{"id":"63c2fd36-a759-4623-a4fc-558e18d2f87c","arxiv_id":"2505.05761","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A solvable correlated model produces a d-wave superconducting dome peaked at the pseudogap quantum critical point, with d-wave pairing robust against phase fluctuations.","lead":"The authors calculate superconductivity in a solvable model of cuprates and find a dome-shaped superconducting region with optimal doping at the point where the pseudogap closes. The model suggests that a partially flat band protects d-wave pairing from phase fluctuations, which may help explain why d-wave superconductivity dominates in cuprates.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The finite-temperature effective dispersion used for all superconducting calculations may not preserve the flat-band condition that underpins the d-wave robustness claim.","rationale":"The paper's central claim is that pseudogap correlations and a partially flat band produce a d-wave superconducting dome with optimal doping near the pseudogap quantum critical point, and that d_{x^2-y^2}-wave pairing is robust against phase fluctuations. All superconducting calculations feed the effective dispersion ξ_k of Eq. (4) into the gap equation, phase stiffness, and density-density correlation, so the accuracy of that replacement is genuinely load-bearing. The reader identified Eq. (4) as the weakest assumption; I agree, and the mechanism analysis sharpens why: the phase-fluctuation robustness argument depends on the flatness of the band near (π,0), but finite-temperature occupation factors can add a U∇_k n_{k+Q} contribution to the centroid velocity that is not present in either exact spectral branch. The paper's own checks are supportive but incomplete: the DOS comparison is shown at two low temperatures, and the exact-vs-effective mean-field T_c comparison removes phase fluctuations, which are exactly what the flat band is supposed to protect against. I also note the authors' acknowledged limitations, including the long-tail overdoped superconductivity and the absence of the antiferromagnetic phase, and the metadata title/author mismatch, but those are secondary to the technical concern. The concern is real and testable, so a conditional verdict remains appropriate; if the proposed two-branch calculation shows the dome and d-wave dominance persist, the paper's central claim would be materially strengthened. Therefore I leave the reader's CONDITIONAL verdict unchanged.","tokens_in":28473,"tokens_out":10943,"duration_ms":123788,"concrete_test":"Recompute Fig. 2(a) at p ≈ 0.374 and p = 0.32 using the exact two-branch spectral function instead of the centroid ξ_k: keep both poles at ε_k and ε_k+U with weights 1−n_{k+Q} and n_{k+Q}, and apply the Doppler shift with v_k = ∂ε_k/∂k to each branch. If the dome, the d-wave/d_xy ratio, or the underdoped discontinuity shifts by more than about 10%, Eq. (4) is the cause. A cheaper diagnostic: at T = 100 and 150 K, compare |∇_k ξ_k| near (π,0) with |∂ε_k/∂k|; a centroid slope much larger than the branch slope would indicate the flat-band protection is partially an artifact of the effective dispersion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (4) replaces the exact two-peak spectral function A(k,ω) = (1−n_{k+Q})δ(ω−ε_k) + n_{k+Q}δ(ω−ε_k−U) by a single centroid band ξ_k = (1−n_{k+Q})ε_k + n_{k+Q}(ε_k+U). At T=0 this is the exact occupied or empty branch, but at finite temperature both branches carry spectral weight. The central d-wave robustness mechanism (main text near Eqs. 5 and 6) relies on the nearly flat band at (π,0): v_k ≈ 0 and a large d-wave gap there make the unpairing condition |p_a·v_k| < sqrt(ξ_k²+|Δ_k|²) almost always satisfied. However, the temperature-dependent n_{k+Q} adds a term U∇_k n_{k+Q} to ∇_k ξ_k, so the centroid band can acquire a non-negligible slope precisely in the momentum region and temperature window that determine T_c. The validation in the End Matter and Fig. 3 checks only the DOS at T = 0.003 K and 69.6 K and the mean-field T_c without phase fluctuations; neither check probes the group velocities or the phase-fluctuation unpairing regions that are the physical basis for d-wave dominance and for the reported underdoped discontinuity on the dome. The most load-bearing assumption is therefore not tested in the regime where the headline claim lives.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies an exactly solvable Hubbard-like model with a pseudogap phase and a partially flat band, and uses it to compute self-consistent d-wave superconductivity with both thermal and zero-point superconducting phase fluctuations in the presence of long-range Coulomb interactions. The central results are a dome-shaped d_{x^2-y^2}-wave superconducting region with optimal doping near the quantum critical point between the pseudogap and metallic phases, a strong suppression of d_{xy}- and s-wave pairing by phase fluctuations, and a discontinuity on the underdoped side of the dome. The authors derive the effective normal-state dispersion ξ_k = (1−n_{k+Q})ε_k + n_{k+Q}(ε_k+U) in Eq. (4), use it for all superconducting calculations, and justify it in the End Matter by comparing the density of states and the zero-phase-fluctuation mean-field T_c against the exact solution.","tokens_in":28727,"tokens_out":4697,"duration_ms":57479,"significance":"If the central claim holds, the paper offers a controlled, analytically solvable route to a cuprate-like phase diagram, including a superconducting dome and robust d-wave symmetry, while explicitly including long-range phase fluctuations and Coulomb interactions, which are usually difficult to treat. The manuscript is commendable for using an exactly solvable correlated model, for providing a detailed path-integral derivation of the phase-fluctuating gap equation in the supplemental material, and for testing the effective dispersion against the exact density of states and the exact mean-field T_c. The main caveat is that the headline results rely on an approximate finite-temperature dispersion whose validation does not directly cover the phase-fluctuation regime where the reported d-wave robustness and the underdoped discontinuity live.","major_comments":[{"comment":"The finite-temperature effective dispersion is load-bearing for the phase-fluctuation calculation, but the validation provided does not test the quantity that controls the d-wave robustness claim. At finite temperature the exact spectral function A(k,ω) = (1−n_{k+Q})δ(ω−ε_k) + n_{k+Q}δ(ω−ε_k−U) has weight in both branches, and the centroid dispersion ξ_k = ε_k + n_{k+Q}U has gradient ∇_k ξ_k = ∇_k ε_k + U∇_k n_{k+Q}; the second term is not controlled by the DOS comparisons in Fig. 3(a) or by the mean-field T_c comparisons in Fig. 3(b), both of which are insensitive to the group velocities that enter the Doppler-shift term p_a·v_k. Since the unpairing criterion |p_a·v_k| < sqrt(ξ_k^2+|Δ_k|^2) near (π,0) is the stated physical reason for d-wave robustness, the authors should either repeat the phase-fluctuation calculation using both spectral branches explicitly or provide a direct check of the group velocities and phase-fluctuation amplitude ⟨p_a^2⟩ in the doping and temperature window where the dome and the discontinuity are obtained.","section":"Eq. (4) and End Matter, 'Justification of the effective electronic energy dispersion'"},{"comment":"The advertised claim that the superconducting dome is 'naturally generated' by the interplay between superconducting order and pseudogap correlations is overstated, because the dome already exists at the BCS level and is explicitly attributed by the authors to the doping-dependent density of states peaked at p = 0.353, a peak inherited from the model's flat band at the quantum critical point. Phase fluctuations modify the dome (including the underdoped discontinuity) but do not create the dome shape. To support the mechanism claimed in the abstract, the authors should either reframe the central claim as 'the pseudogap-induced partially flat band produces a DOS peak that controls the dome location' or provide a control calculation, such as a band structure without the flat band, showing that the QCP location rather than the input DOS determines the optimal doping.","section":"Abstract and Fig. 2(b)"},{"comment":"The first-order-like transition from a strong-phase-fluctuating state to a weak-phase-fluctuating state on the underdoped side of the dome is a central new feature, but its numerical robustness is not documented. The discontinuity appears in a narrow doping window between p = 0.31 and p = 0.32, and its location and sharpness depend on the phase-fluctuation integral, including the momentum cutoff q_c = 1/ξ0 introduced in the End Matter. The authors should report a convergence check in the k-grid, the q-integration cutoff, and the iterative solution procedure, and should state whether the discontinuity is a true solution-branch change or a numerical jump in the self-consistent iteration.","section":"Fig. 2(a) and text near the underdoped discontinuity"}],"minor_comments":[{"comment":"The full-text title, 'An Exactly Solvable Model of Phase-Fluctuating Superconductivity in Cuprates: The Role of Partially Flat Bands', differs from the title under which the paper is submitted; the two should be harmonized.","section":"Title and abstract"},{"comment":"The upper cutoff q_c = 1/ξ0 for the bosonic phase-fluctuation integral is introduced without any sensitivity analysis; a brief statement of how the dome, the T_c values, and the discontinuity depend on this cutoff would make the calculation more reproducible.","section":"Eq. (7) and End Matter, phase-fluctuation cutoff"},{"comment":"The pseudogap onset temperature T* is defined through a Lorentzian broadening Γ = 0.01t and a 10% spectral-weight threshold; the sensitivity of the reported T* line, and hence of its intersection with the superconducting dome, to these numerical choices should be stated.","section":"Supplemental Material, Section I"},{"comment":"The statement that U = 1.15t in the exactly solvable model corresponds to U = 8t in the conventional Hubbard model is not derived; a brief explanation of the mapping would help readers assess the physical regime.","section":"Discussion and Ref. [71]"},{"comment":"The notation for the pairing symmetries is not fully consistent (d_{x^2-y^2}-wave appears as 'dx2−y2-wave' and 'd_{x^2-y^2}-wave' in different places); please unify the notation.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and addresses a timely problem. The main risk is that the central d-wave robustness result rests on an effective finite-temperature dispersion whose group velocities are not validated in the phase-fluctuation regime; however, this is addressable with additional calculations or a more careful bounding analysis. I would not recommend rejection. The title mismatch between the abstract and the full text, and the fact that the manuscript is dated August 2025 under a May 2025 arXiv identifier, should be checked by the editors."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is applying the Yang-Chen phase-fluctuating superconductivity formalism to the Worm et al. exactly solvable pseudogap model, and working out d-wave pairing in that setting. The main results are a superconducting dome peaking near the pseudogap quantum critical point, and a claim that d_{x^2-y^2}-wave pairing is robust against phase fluctuations while d_xy and s-wave are suppressed. That calculation is not in the prior literature, which had the normal state (Worm et al.) or s-wave pairing (Phillips et al.). So the paper does produce a new result.\\n\\nWhat it does well: the self-consistent scheme is laid out carefully, the gap equation and phase stiffness are derived in the supplement, and the authors provide two honest checks of their key approximation—comparing the effective dispersion's density of states against the exact model at two temperatures, and comparing mean-field T_c against an exact diagonalization of the 16-state problem. They also state plainly, in the main text, that the dome shape comes from the doping-dependent normal-state DOS peaked at the QCP. That undercuts the abstract's \"interplay\" framing, but it is not hidden. The treatment of zero-point fluctuations via a renormalized pairing strength is a sensible way to avoid double counting.\\n\\nThe soft spots are real but proportionate. The effective dispersion in Eq. (4) is load-bearing, and the stress-test concern lands: at finite temperature, the centroid band includes U\\nabla_k n_{k+Q}, which can destroy the near-flatness at (\\pi,0) that the d-wave robustness argument depends on. The validation checks show DOS agreement and mean-field T_c agreement, but neither probes the group velocities or the unpairing regions that determine whether phase fluctuations actually favor d-wave. That is a gap between the evidence and the central claim. The circularity concern is partly justified but mild, since the authors say the dome follows from the DOS peak; the abstract overstates by calling it \"natural generation by interplay.\" There is also a metadata problem worth flagging: the arXiv listing gives different authors and a different title from the full text. That has to be fixed regardless of scientific content. No code or data is shipped, and experimental comparisons are qualitative, which is fine for this kind of theory paper but should be said clearly.\\n\\nWho is this for? People who work on solvable correlated models for the cuprate pseudogap and on phase-fluctuation mechanisms. It deserves a serious referee: the framework is controlled, the derivations are explicit, and the d-wave robustness claim is interesting enough to warrant scrutiny. I would send it to review, and ask the authors to test the effective dispersion by computing group velocities in the regions that matter, or when feasible to run the gap equation with the full two-branch spectral function. Fix the metadata, soften the abstract, and this is a reasonable contribution to a specialized journal.","headline":"A credible and clearly written phase-fluctuating d-wave calculation on a solvable pseudogap model; the dome is largely inherited from the input flat band, and the effective dispersion assumption needs stronger checking before the d-wave robustness claim is taken as established.","tokens_in":29319,"tokens_out":1762,"would_cite":true,"duration_ms":21242,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":null,"created_at":"2026-08-15T22:57:22.153916+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":null,"supporting_citations":[],"review_version":1}