{"id":"65526f97-aceb-4049-a6f9-01a47f1c906a","arxiv_id":"2506.15770","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Simulations of the doped Hubbard model show a momentum-dependent Mott-gap dichotomy and mismatched pseudogap temperatures from different probes, supporting a crossover picture of the cuprate pseudogap.","lead":"This paper uses exact quantum Monte Carlo simulations of the Hubbard model to map how electrons behave in a model relevant to high-temperature superconductors, finding clear differences between hole and electron doping. It reports that the mysterious pseudogap, a state where parts of the Fermi surface are gapped, looks like a smooth crossover caused by strong interactions rather than a sharp phase transition.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The no-pocket and T*_A evidence rests on un-error-barred MaxEnt self-energies and a three-point spline peak; either artifact would weaken the crossover conclusion.","rationale":"I read the paper as a serious DQMC study whose asymmetry and dichotomy results are well supported: DQMC is unbiased, the -beta G(k,tau=beta/2) proxy independently confirms the Fermi-surface trends, and data and analysis code are released. The concern is not with the numerical core but with the interpretive step from those data to 'smooth crossover, no pockets.' The reader conditioned on MaxEnt reliability and 8x8 finite-size effects; I agree and locate the most fragile point more precisely. The self-energy zero-crossing that rules out pockets is a derived object: MaxEnt gives A(k,omega), Kramers-Kronig gives Re G, Dyson inversion gives Sigma, and splines interpolate it. Each step can smooth or shift a narrow feature, and the paper states that error bars are not shown for spectral functions and self-energies. A pocket back side would be a weak, narrow feature near the antinode, exactly where the method is least reliable. Separately, T*_A relies on a cubic-spline peak through a very small number of temperature points; without error bars, the 'peak' at T/t = 1/3 is not established. Both issues are correctable with additional computation: an injected-pocket recovery test and denser temperature sampling with bootstrapped errors. If those tests pass, the crossover claim is substantially strengthened; if they fail, the claim should be weakened to 'consistent with a crossover within the accessible temperature and momentum resolution.' The reader's CONDITIONAL verdict already captures this need, so I recommend no change.","tokens_in":34538,"tokens_out":6183,"duration_ms":76732,"concrete_test":"Using the released DQMC data and code, rerun the n=0.95, U=6, t'=-0.25, T=t/4 analysis with a synthetic antinodal pocket injected into G(k,tau) (a second pole in Sigma with width comparable to or smaller than the twist momentum spacing), then apply the identical MaxEnt, Kramers-Kronig, and spline pipeline. If the pipeline does not recover the second zero-crossing in Re Sigma + epsilon - mu, the no-pocket conclusion is not falsifiable by this method. Independently, add two additional temperature points near T/t = 0.29 and 0.22 for antinodal A(k_F,0) and recompute T*_A with bootstrap error bars; if the peak at T/t = 1/3 shifts or disappears, the probe-dependent-T* argument is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is a smooth crossover rather than a symmetry-breaking transition, and the two load-bearing pieces of evidence are (i) the absence of Fermi pockets from the sign structure of Re Sigma_k(omega=0) + (epsilon_k - mu), and (ii) the probe dependence of T*, anchored by T*_A defined as the antinodal spectral-weight peak. Step (i) is the more fragile. The self-energy is obtained by MaxEnt continuation of G(k,tau), Kramers-Kronig inversion, Dyson inversion, and then 2D spline interpolation over momenta from 15 twisted 8x8 boundary conditions. A small pocket back side would appear as a second zero-crossing of Re Sigma + epsilon - mu in a narrow momentum window near the antinode, exactly where A(k,0) is smallest and where both MaxEnt bias and spline smoothing are most likely to erase it. Thus 'no signature of pocket formation' is a resolution/continuation upper limit, not a demonstrated absence. Step (ii) is also under-resolved: for n=0.95, T*_A = t/3 is identified from a cubic-spline peak through roughly three temperature points (T/t <= 0.5, 0.33, 0.25) with no error bars, so the antinodal nonmonotonicity is not independently confirmed. If either artifact is real, the evidence that different probes yield different T* and that the pseudogap is a smooth crossover loses its footing.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses determinant quantum Monte Carlo (DQMC) with twisted boundary conditions and maximum-entropy analytic continuation to compute single-particle spectral functions, self-energies, NMR-like spin susceptibilities, transport, and thermodynamics of the two-dimensional Hubbard model with U/t = 6, t'/t = -0.25, and t''/t = 0 or 0.15, on 8x8 lattices for both hole and electron doping. It reports a systematic electron-hole asymmetry: the electron-doped side exhibits stronger antiferromagnetic correlations, more coherent quasiparticles, and hot spots on the Fermi surface, while the hole-doped side develops a nodal-antinodal dichotomy and a pseudogap without signatures of Fermi pockets. The authors attribute this dichotomy to the momentum dependence of the Mott gap, as seen in both the spectral function and a Re Sigma ~ 1/omega self-energy, and they compare pseudogap temperatures extracted from the antinodal spectral weight, Knight shift, and (T1T)^-1, finding that these temperatures do not coincide. They conclude that the pseudogap is a smooth crossover driven by strong correlations, not a symmetry-breaking phase transition, even without disorder.","tokens_in":34847,"tokens_out":5034,"duration_ms":60343,"significance":"If the central claims hold, this is a valuable contribution to the pseudogap debate. The paper's strengths include systematic DQMC simulations with twisted-boundary-condition momentum resolution, a cross-check of the MaxEnt Fermi surfaces against the continuation-free proxy -beta G(k, tau = beta/2), direct comparison of spectroscopy-like and NMR-like probes within the same model, and public release of data and analysis routines. The electron-hole asymmetry and the explicit statement that the pseudogap temperature is probe-dependent are interesting and testable, and the comparison across independent probes is not circular. However, the load-bearing negative claim of no Fermi pockets and the quantitative crossover temperatures rest on MaxEnt spectral functions and self-energies without error bars, and on finite-size assumptions that are acknowledged but not checked for the spectral quantities. The significance is therefore conditional on resolving these points.","major_comments":[{"comment":"The central negative claim that no Fermi pockets appear is a resolution-limited statement. The quantity Re Sigma_k(omega = 0) + (epsilon_k - mu) is obtained through MaxEnt continuation, Kramers-Kronig inversion, Dyson inversion, and two-dimensional spline interpolation over momenta from 15 twisted 8x8 boundary conditions; Appendix B states that error bars are not shown for single-particle spectral functions and self-energies. A small pocket back side would appear as a second zero crossing in a narrow momentum window near the antinode, precisely where A(k,0) is smallest and where MaxEnt bias and spline smoothing are most likely to erase it. Please provide a quantitative test of this resolution limit, for example by applying the same continuation and interpolation pipeline to synthetic G(k,tau) generated from spectra with and without a small pocket, or by showing convergence with respect to the number of twists and lattice size.","section":"Section III, Figs. 8-10; Appendix B"},{"comment":"The value T*_A = t/3 for n = 0.95 is identified from a cubic-spline peak of the antinodal spectral weight versus temperature using essentially three temperature points (T/t = 1/2, 1/3, 1/4) and no error bars. Because the probe-dependence conclusion in Fig. 14 hinges on this peak, please provide bootstrap confidence intervals or additional temperature points to confirm that the antinodal nonmonotonicity is real and that the peak position is robust.","section":"Section IV, Figs. 7(d) and 14"},{"comment":"The paper acknowledges that finite-size effects are more pronounced for electron doping and states that they 'are likely to extend to other properties not explicitly analyzed here, such as the single-particle spectral function, self-energy, and (T1T)^-1,' yet no finite-size check is presented for the no-pocket result or for T*_A. Given that the no-pocket and crossover conclusions are the central claims, please add a 10x10 or 12x12 comparison for the hole-doped n = 0.95 case, or an equivalent twist-convergence study, to demonstrate that the real-frequency self-energy structure is converged.","section":"Appendix C4 and Figs. C8, 9-10"},{"comment":"The text states that for t''/t = 0 the Knight-shift peak 'occurs at much lower temperatures, even outside the range of our study,' but Fig. 14 reports T*_Ks values for t''/t = 0. Please clarify how T*_Ks was extracted and reconcile these statements, since the summary comparison of T* values depends on the consistency of these definitions.","section":"Section IV and Fig. 14"}],"minor_comments":[{"comment":"The abstract says 'transition towards the pseudogap' while the conclusion says the pseudogap is a smooth crossover; please use consistent terminology to avoid implying a phase transition.","section":"Abstract and Section V"},{"comment":"The caption says dashed lines highlight the change in sign of Re Sigma_k(omega = 0) + (epsilon_k - mu); please state explicitly whether these dashed lines are zero contours or merely guides to the eye.","section":"Fig. 9 caption"},{"comment":"The Knight shift and (T1T)^-1 are plotted in arbitrary units; please state in the captions that constant prefactors are omitted and whether the quantities are normalized per site.","section":"Appendix A3 and Figs. 12-13"},{"comment":"The notation for twisted-boundary-condition Green's functions is compact and could be confusing; a one-sentence reminder that the twist angle enters through the redefined operators in Eq. (A18) would improve readability.","section":"Appendix A2 and Eq. (A21)"},{"comment":"The footnote-like reference [132] embedded in a sentence about interpolation would be clearer as a regular parenthetical or a numbered footnote in the standard style.","section":"Appendix A1, footnote [132]"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this if you want the current best DQMC evidence on hole-electron asymmetry and probe-dependent T* in the Hubbard model. The paper is honest about its own limitations, and the central crossover claim is plausible but not airtight.\n\nWhat earns credit: the twisted-boundary spectral maps are a real step up in momentum resolution; the comparison of T* from spectra, Knight shift, and T1 across a doping ladder is new and well-organized; and the asymmetry in coherence between hole and electron doping is backed by transport, thermodynamics, and spin structure, not just spectra. Data and code are public. The use of -beta G(k,tau=beta/2) as a MaxEnt-free cross-check is good practice.\n\nThe soft spots land where the stress-test says. The no-pocket conclusion rests on MaxEnt self-energies with no error bars (stated in Appendix B), followed by Kramers-Kronig and Dyson inversion and 2D spline interpolation over 15 twisted 8x8 boundary conditions. A small pocket back side near the antinode, where spectral weight is weakest, could be erased by continuation bias or spline smoothing. The authors say 'no signature of pocket formation', which is carefully phrased, but the later statement that the results contradict pocket-based arc explanations is stronger than the evidence supports. Second, T*_A for the key doping n=0.95 is a cubic-spline peak through three temperature points, with no error bars; the nonmonotonicity is visible in the raw curves, but the quoted value is not precisely determined. Third, finite-size effects are explicitly said to likely extend to spectral quantities, yet the main maps are 8x8. None of this kills the paper; it means the headline should be read as 'consistent with a crossover' rather than 'proven'.\n\nWorth a serious referee. The experimental backbone is solid, the parameter scans are systematic, and the comparison of T* definitions is a useful contribution to a long-running debate. A referee should push for a sensitivity analysis on the no-pocket claim and for error estimates on the spectral functions, but this is a paper that belongs in the literature. I would cite it for the T* comparison and the asymmetry results.","headline":"Solid DQMC study with honest caveats; crossover claim is plausible but the no-pocket and T*_A evidence are resolution-limited.","tokens_in":35406,"tokens_out":3184,"would_cite":true,"duration_ms":39048,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.27.+a","74.72.Kf"],"model":"deepseek-v4-flash","headline":"This paper argues, from unbiased numerically exact simulations of the doped Hubbard model, that the cuprate pseudogap is a smooth crossover driven by strong correlations — not a symmetry-breaking phase transition.","keywords":["Hubbard model","pseudogap","determinant quantum Monte Carlo","maximum entropy analytic continuation","Fermi arcs","nodal-antinodal dichotomy","Mott gap","cuprate superconductors"],"falsifier":"Run the same model at temperatures below $T/t = 1/4$ on larger clusters with a real-frequency method that avoids maximum-entropy reconstruction: if the antinodal spectral weight develops a kink that sharpens as the lattice grows, if a Fermi pocket with a resolvable back side appears in the self-energy, or if the three $T^*$ values converge to a single lattice-size-independent number, the crossover conclusion would fail, whereas a smooth, size-independent evolution of the antinodal weight is what the paper's picture predicts.","tokens_in":34323,"feed_emoji":"⚛","tokens_out":22470,"duration_ms":202688,"temperature":0.7,"pith_summary":"This paper sets out to settle a foundational debate in cuprate physics: whether the pseudogap — the suppression of low-energy electronic states seen below a temperature $T^*$ in hole-doped cuprates — is a distinct phase of matter or a smooth, disorder-free crossover. Using determinant quantum Monte Carlo, a numerically exact simulation of the Hubbard model, the authors compute momentum-resolved single-particle spectra, the simulated analogue of angle-resolved photoemission spectroscopy. They find that hole doping produces Fermi arcs and a nodal–antinodal dichotomy that mirrors the momentum dependence of the Mott gap, showing up as a pole-like $\\mathrm{Re}\\,\\Sigma \\sim 1/\\omega$ self-energy, whereas electron doping yields coherent quasiparticles and antiferromagnetic hot spots with no pseudogap. Because spectroscopy and simulated nuclear magnetic resonance give different pseudogap temperatures, and the simulations contain no disorder, the paper concludes that the pseudogap onset is a smooth crossover driven by strong correlations rather than a phase transition. If correct, this tells experimenters that probe-dependent $T^*$ values are expected, and it redirects theory from broken-symmetry order toward strong-correlation physics as the explanation.","feed_headline":"Simulations find the pseudogap is a smooth crossover, not a transition","feed_subtitle":"Probe-dependent pseudogap temperatures in disorder-free simulations point to strong-correlation physics, not broken symmetry.","key_machinery":"The load-bearing object is the single-particle self-energy $\\Sigma_k(\\omega)$, extracted from the simulated spectral function $A(\\mathbf{k},\\omega)$ through Dyson's equation $G^R(\\mathbf{k},\\omega) = 1/(\\omega - (\\epsilon_{\\mathbf{k}}-\\mu) + i\\delta - \\Sigma_{\\mathbf{k}}(\\omega))$. An explicit two-pole calculation in the paper's appendix shows that a pole-like real part $\\mathrm{Re}\\,\\Sigma \\sim (\\omega-\\omega_0)^{-1}$ opens a gap of size roughly $2a^{-1/2}$ at the Fermi level, while a linear negative slope $\\mathrm{Re}\\,\\Sigma \\sim -b\\omega$ merely renormalizes the band and leaves it gapless; the pseudogap appears where the pole-like remnant of the Mott gap sits closest to the Fermi level — near the antinode for hole doping, near the node for electron doping. Methodologically, the argument rides on three tools: determinant quantum Monte Carlo (a numerically exact simulation of the interacting Hubbard model), maximum-entropy analytic continuation (which converts imaginary-time Green's functions and spin susceptibilities into real-frequency spectra), and twisted boundary conditions on an $8\\times 8$ cluster (which raise the momentum resolution to that of a $64 \\times 64$ lattice).","core_discovery":"On the paper's own terms, the central discovery is that the hole-doped pseudogap of the Hubbard model is a smooth crossover driven by strong correlations, and that its momentum-space signatures — Fermi arcs and a nodal–antinodal dichotomy — follow directly from the momentum dependence of the Mott gap rather than from antiferromagnetic order, Fermi pockets, or a Lifshitz transition. The decisive evidence sits in the self-energy: at low hole doping the real part $\\mathrm{Re}\\,\\Sigma_k(\\omega)$ develops a pole-like $1/(\\omega-\\omega_0)$ remnant of the Mott gap whose position relative to the Fermi level shifts with momentum, suppressing coherent weight at the antinode while the node stays coherent; electron doping reverses the dichotomy, and the same simulations show no enclosed pockets anywhere in the Brillouin zone. The temperature dependence of the antinodal spectral weight defines a spectroscopic pseudogap onset $T^*_A \\sim t/3$, while simulated NMR probes — the Knight shift and $(T_1T)^{-1}$ evaluated with oxygen-site form factors — give their own onset temperatures $T^*_{K_s}$ and $T^*_{T_1}$ that do not coincide with $T^*_A$. Because no disorder enters the calculation, the paper concludes that the pseudogap is a true smooth crossover even in the clean limit, describing a high-temperature pseudogap that is visible only for dopings near half-filling.","pith_inferences":["Testable extension: map the three onset temperatures $T^*_A$, $T^*_{K_s}$, and $T^*_{T_1}$ across a grid of $U$, $t'$, and $t''$ values — the crossover picture predicts the probe-dependence survives over a wide parameter range, since the mechanism is the momentum-dependent Mott gap rather than a special hopping geometry.","The pole-like self-energy implies a zero of the Green's function near the same frequency, so re-analyzing the same simulation data through the zeros of $G(\\mathbf{k},\\omega)$ could provide a sharper, MaxEnt-robust diagnostic for the crossover than the peak-based definition of $T^*_A$.","The simulations stop at $T/t \\ge 1/4$; the crossover conclusion implicitly invites a direct check that the antinodal spectral weight stays smooth as the lattice grows and the temperature drops, with no emerging kink or pocket back side.","If the crossover view is right, the two temperature scales seen in NMR experiments on real cuprates are likely two smooth onsets of the same Mott-controlled physics rather than signatures of two distinct orders — a reinterpretation that simultaneous ARPES and Knight-shift measurements on the same sample could test."],"forward_implications":["Different probes will legitimately define different pseudogap temperatures $T^*$: the spectroscopic onset, the NMR Knight-shift onset, and the spin-lattice-relaxation onset are distinct scales, so apparent disagreements among ARPES, NMR, and transport experiments are what a correlation-driven crossover predicts rather than evidence of a hidden transition.","Fermi arcs are not the visible halves of damped Fermi pockets: the self-energy analysis finds no enclosed pockets, ruling out the arc-from-pocket-back-side mechanism at the studied dopings and temperatures.","Weak-coupling fermiology — band structure plus antiferromagnetic nesting — correctly describes the electron-doped side (coherent quasiparticles, hot spots) but fails for hole doping, where the antinodal self-energy is governed by the remnant Mott gap; hole-doped pseudogap physics is a strong-coupling effect.","Hole doping reproduces the linear-in-$T$ resistivity associated with strange metallicity, while electron doping curves toward Fermi-liquid behavior, connecting the coherence asymmetry between the two dopings to transport anomalies.","The pseudogap seen here is the high-temperature pseudogap, appearing only near half-filling; the paper leaves open a distinct lower-temperature pseudogap away from half-filling with a much smaller energy scale."],"supporting_citations":[{"why":"Supplies the determinant quantum Monte Carlo algorithm, the numerically exact method behind all of the paper's many-body simulations.","marker":"[86]"},{"why":"Establishes the DQMC approach for the two-dimensional Hubbard model that this study applies to hole and electron doping.","marker":"[87]"},{"why":"Provides the maximum-entropy analytic continuation that turns imaginary-time Green's functions into the real-frequency spectral functions and self-energies.","marker":"[88]"},{"why":"Supplies the analytic continuation methodology used for the spectral, transport, and NMR quantities extracted from imaginary-time data.","marker":"[89]"},{"why":"Introduces the twisted-boundary-conditions scheme on small Hubbard clusters and the strong-coupling pseudogap mechanism the paper builds on for momentum-resolved spectra.","marker":"[84]"},{"why":"Gives the prior DQMC simulation of the NMR response in the pseudogap regime whose form-factor treatment the paper adopts and whose conclusions its T* values contrast with.","marker":"[81]"},{"why":"Reports the linear-in-T resistivity of the doped Hubbard model that the present hole-doped transport results reproduce as a signature of strange metallicity.","marker":"[108]"},{"why":"States the Fermi-arc-from-damped-Fermi-pocket scenario that the paper's self-energy analysis, finding no pockets, directly contradicts.","marker":"[44]"},{"why":"Proposes the pseudogap-Lifshitz-transition link that the paper finds no evidence for at the relevant dopings.","marker":"[45]"}],"fun_headline_variants":["Pseudogap is a smooth crossover, not a transition","Hubbard model shows pseudogap: smooth crossover, no pockets","Fermi arcs from Mott gap momentum dependence, not AF order","Electron vs hole doping: asymmetric coherent weights and T*"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the computer reconstruction of real-frequency spectra from the imaginary-time simulation data is quantitatively reliable — including the pole-shaped feature in the self-energy that carries the argument — and that an $8\\times 8$ lattice with twisted boundary conditions stands in for the infinite system, a premise the paper itself flags by noting that no error bars are shown for spectral functions and that finite-size effects are more pronounced for electron doping and likely extend to the spectra.","fun_headline_variants_meta":{"raw":{"variants":["Pseudogap is a smooth crossover, not a transition","Hubbard model shows pseudogap: smooth crossover, no pockets","Fermi arcs from Mott gap momentum dependence, not AF order","Electron vs hole doping: asymmetric coherent weights and T*"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000327,"raw_usage":{"total_tokens":1844,"prompt_tokens":976,"completion_tokens":868,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":592,"completion_tokens_details":{"reasoning_tokens":795}},"tokens_in":592,"tokens_out":868,"duration_ms":8931,"temperature":1.0,"reasoning_tokens":795,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T23:51:56.474850+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same model at temperatures below $T/t = 1/4$ on larger clusters with a real-frequency method that avoids maximum-entropy reconstruction: if the antinodal spectral weight develops a kink that sharpens as the lattice grows, if a Fermi pocket with a resolvable back side appears in the self-energy, or if the three $T^*$ values converge to a single lattice-size-independent number, the crossover conclusion would fail, whereas a smooth, size-independent evolution of the antinodal weight is what the paper's picture predicts.","supporting_citations":[{"cited_title":"Blankenbecler, D","cited_arxiv_id":null,"evidence_quote":"Supplies the determinant quantum Monte Carlo algorithm, the numerically exact method behind all of the paper's many-body simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the DQMC approach for the two-dimensional Hubbard model that this study applies to hole and electron doping."},{"cited_title":"Jarrell and J","cited_arxiv_id":null,"evidence_quote":"Provides the maximum-entropy analytic continuation that turns imaginary-time Green's functions into the real-frequency spectral functions and self-energies."},{"cited_title":"Gunnarsson, M","cited_arxiv_id":null,"evidence_quote":"Supplies the analytic continuation methodology used for the spectral, transport, and NMR quantities extracted from imaginary-time data."},{"cited_title":"Strong-coupling mechanism of the pseudogap in small Hubbard clusters","cited_arxiv_id":"2010.12601","evidence_quote":"Introduces the twisted-boundary-conditions scheme on small Hubbard clusters and the strong-coupling pseudogap mechanism the paper builds on for momentum-resolved spectra."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the prior DQMC simulation of the NMR response in the pseudogap regime whose form-factor treatment the paper adopts and whose conclusions its T* values contrast with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the linear-in-T resistivity of the doped Hubbard model that the present hole-doped transport results reproduce as a signature of strange metallicity."}],"review_version":1}