{"id":"ffc71c01-df5b-4a63-8088-a9437c363a35","arxiv_id":"2412.18864","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The Hatsugai-Kohmoto model with pairing interactions shows strong pair-density-wave fluctuations at momentum (π,π) for low doping and intermediate interaction, competing with ordinary superconductivity.","lead":"This paper calculates pairing instabilities in an exactly solvable toy model of a doped Mott insulator, the Hatsugai-Kohmoto model with added pairing interactions, and finds finite-momentum pair-density-wave fluctuations at wavevector (π,π) in a low-doping, intermediate-interaction regime. The result offers a simple analytic setting for a debated kind of superconducting order in high-temperature cuprates, where fluctuating pair-density waves may be a key ingredient.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central PDW/SC phase boundary rests on an untested constant-density-of-states approximation; the true square-lattice DOS may change the competition.","rationale":"The reader's weakest-assumption analysis identifies the same load-bearing step: the constant-DOS approximation. I agree that this is the single most concerning point because the central claim is not just the existence of a finite-Q pairing susceptibility but its dominance over the q=0 channel in a specific part of the phase diagram. That dominance is decided by comparing two integrals that are both evaluated with the flat DOS, and the paper explicitly acknowledges the approximation without testing it. A direct numerical or analytic-DOS recomputation would settle the issue. Secondary issues noted in the paper, such as the likely typo in Eq. (14) and the partially derived appendix formulas, are presentation problems rather than independent threats to the central claim, and the arbitrary choice v=v' is transparent and can be varied once the underlying χ0 are computed correctly. The verdict of CONDITIONAL remains appropriate: the claim is plausible and internally coherent, but its quantitative phase boundary is not yet robustly established.","tokens_in":17109,"tokens_out":16956,"duration_ms":165013,"concrete_test":"Recompute T_c^PDW and T_c^SC for v=v'=0.75 at representative (u,x) values using the exact square-lattice DOS, either by direct k-grid evaluation of Eqs. (26)/(30) and (27) with the finite-T occupation factors and Eq. (6), or by substituting the analytic DOS ρ(ω) = K(1-(ω/4t)^2)/(2π^2 t) into Eqs. (19), (27), and (30), then re-extract the PDW/SC boundary shown in Fig. 4. If the low-x intermediate-u window with T_c^PDW > T_c^SC persists and x_c remains within roughly 20% of 0.25, the flat-DOS concern is resolved; if the window disappears or shifts substantially, the central claim must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that χ0(0,Q) diverges first in a low-doping, intermediate-U window (Sec. IV, Fig. 4) is computed using the replacement ρ(ω)=1/W for the square-lattice density of states, stated explicitly in Sec. III as 'for simplicity'. This flat DOS enters every quantitative step that produces the phase diagram: the pair-binding equation (19), the SC susceptibility (27), the PDW susceptibility (30)-(32), and the chemical-potential versus doping relations (33)-(34). The real square-lattice DOS has a logarithmic van Hove singularity at band center and sharp band edges, which change the relative weight of the Ω0, Ω1, and Ω2 momentum regions and also shift μ(T,x,u) at fixed doping. Because the SC and PDW integrals probe different energy ranges, the ratio χ0(0,Q)/χ0(0,0) can change when the true DOS replaces the flat one. The paper offers no numerical k-sum test or error estimate showing that the PDW-dominated window and the quoted x_c≈0.25 survive. If the window shrinks or disappears under the actual DOS, the main conclusion is an artifact of the approximation rather than a property of the HK-SC model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the Hatsugai-Kohmoto (HK) model on a square lattice with an added finite-momentum pairing interaction (the HK-SC model). It derives Cooper-pair binding energies and static pair susceptibilities in the zero-center-of-mass (SC) and Q=(π,π) (PDW) channels. The central claim is that, for intermediate U and low doping up to a critical doping of order x_c≈0.25, the PDW susceptibility diverges at a finite temperature and dominates the SC channel, while at larger doping SC wins. The authors compute closed-form expressions for the pair-binding energies, critical temperatures, and chemical potentials using a constant density of states, and interpret the resulting phase diagram as evidence that PDW fluctuations emerge generically on doping a Mott insulator, with possible relevance to underdoped cuprates.","tokens_in":17288,"tokens_out":5210,"duration_ms":51739,"significance":"If the central claim survives scrutiny, this would be a valuable addition to the field: it provides an analytically tractable two-dimensional model in which PDW fluctuations compete with uniform superconductivity, and it connects the exactly solvable HK model to a currently active experimental and theoretical topic. The paper is constructive rather than phenomenological: it contains no parameter fitting, the susceptibility calculation is transparently built from previously established HK Green's functions, and the analytic expressions are sufficiently detailed to be checked and reused. The significance is somewhat conditional, however, because the quantitative phase diagram is obtained from an approximate density of states, and because the treatment of the pairing-perturbed model is not literally exact.","major_comments":[{"comment":"The central quantitative claim that PDW fluctuations dominate in a low-doping, intermediate-U window rests on the replacement of the square-lattice density of states by a constant, ρ(ω)=1/W, introduced in Sec. III immediately after Eq. (15) and used in Eq. (19), Eq. (27), Eqs. (30)-(32), and the chemical-potential relations (33)-(34). The real square-lattice DOS has van Hove singularities and sharp band edges, which change the relative weight of the Ω0, Ω1, and Ω2 energy regions and also shift μ(T,x,u) at fixed doping. Since the q=0 and q=Q susceptibilities probe different energy ranges, the ratio χ0(0,Q)/χ0(0,0) could change substantially when the true DOS replaces the flat one. The paper provides no numerical k-sum test, no comparison with the actual square-lattice DOS, and no error estimate showing that the PDW-dominated window and the quoted critical doping x_c≈0.25 survive. Because this approximation enters every equation that produces the phase diagram in Fig. 4, this is a load-bearing point that needs to be addressed before the central claim can be accepted.","section":"Sec. III, Eq. (19); Sec. IV, Eqs. (27), (30)-(32)"},{"comment":"Equation (14) appears to contain an error in the Ω1 terms: the factors (1−n^q_0(k)) should presumably be (1−n^q_1(k)), since those terms involve β_k with k∈Ω1 and since the variational equation (15) for the same channel uses (1−n^q_1(k)) and n^q_1(k)+1/2. As written, Eq. (14) is inconsistent with Eq. (15) and with the derivation of the pair-binding energy. The authors should correct this expression and re-verify that the self-consistent equation (15) follows from the corrected variational calculation.","section":"Eq. (14)"},{"comment":"The claim of exactness is overstated for the HK-SC model. The bare HK model is exactly solvable because it is diagonal in momentum, but once the pairing term V_q is added the Hamiltonian is no longer k-diagonal. Equation (22) is a Dyson/RPA-type resummation for the pairing susceptibility, and the subsequent calculation explicitly assumes T≪U,W and the constant-DOS approximation. The text at the start of Sec. IV says 'We will perform here an exact calculation of Eq. (21)', and the title and abstract describe the doped Mott insulator as 'exactly solvable'. I recommend either proving that Eq. (22) is exact for the HK-SC model or clearly labeling the susceptibility computation as an approximate controlled treatment. This distinction matters because the critical temperatures T_c^SC and T_c^PDW are inferred from the pole condition of Eq. (23) within this approximation.","section":"Sec. IV, Eq. (22), and title/abstract"}],"minor_comments":[{"comment":"The sentence 'In Appendix V, we show some expressions of μ(T,x,U)' appears to be a typo; there is no appendix labeled V in the manuscript, only the unnumbered 'APPENDIX: RESULTS FOR T_c^PDW AND μ(T)'.","section":"Sec. II"},{"comment":"The expression for χ0(0,Q) in Eq. (32) contains ambiguous notation, including 'tanh^{-1}(βU/2)', unbalanced parentheses, and line breaks that make the formula difficult or impossible to verify. Please provide a cleaner closed form, or include a supplementary notebook or code snippet that generates the plotted curves.","section":"Appendix, Eq. (32)"},{"comment":"The notation μ'_c appears in the Appendix without being defined in the main text, and the sentence 'We note that μ′c in Eqs. (31) and (34)' seems to refer to a quantity in Eq. (34) that is not denoted μ'_c. Please make the notation for the chemical potential at the critical temperature consistent throughout Eqs. (31)-(34).","section":"Appendix, Eqs. (31)-(34)"},{"comment":"Figure 4 is labeled 'schematic phase diagram' with a dashed guide line. Since the precise location of the PDW/SC boundary is a central quantitative result, it would be helpful to either show computed boundary curves for the parameter values used in the paper, or explicitly state that the drawn boundary is only illustrative and not extracted from the calculations.","section":"Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The main risk is the flat-DOS approximation. If the authors can supply a numerical check with the true square-lattice DOS and the PDW-dominated window survives, the paper would be a solid contribution. If such a check shows the window shrinks or disappears, the central claim would need to be substantially weakened. I would also ask the editor to encourage the authors to soften the 'exact' language for the HK-SC susceptibility calculation, as this wording is likely to attract criticism."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThe short version: this paper does something new—it extends the Hatsugai-Kohmoto exact Green's function to the finite-momentum pairing channel at Q=(π,π) and computes the PDW susceptibility and the PDW/SC competition. That's a real extension of the known q=0 HK superconductivity results, and the analytic treatment is careful within its stated assumptions. If you work on HK models or on PDW mechanisms in exactly solvable settings, this is worth reading.\n\nWhat's good: the derivation of the pair-binding energy and susceptibilities is mostly clear, the use of the exact HK Green's function is appropriate, and the authors are upfront that the constant-DOS replacement is a simplification. The qualitative picture—PDW fluctuations win at low doping and intermediate U, SC wins further out—is internally consistent, and the appendix gives explicit formulas, even if they are ugly.\n\nThe soft spots are real but not fatal. The biggest one is the constant density of states. The entire phase diagram, including the quoted x_c≈0.25, comes from replacing the square-lattice DOS with a flat band. The real DOS has van Hove singularities and sharp edges, and the SC and PDW integrals probe different energy ranges, so the ratio can shift. The paper offers no check (numerical k-sum or otherwise) that the PDW window survives with the true DOS. That is a genuine gap, and it is the main reason I would not take the quantitative boundaries at face value.\n\nSecond, the word 'exact' is doing more work than the analysis supports. The susceptibility is computed with the exact Green's function of the bare HK model, but the pairing interaction is then handled through a ladder/Dyson equation with a finite V_q. That is standard BCS-type treatment, not an exact solution of the HK-SC Hamiltonian. The title is stronger than the method.\n\nThird, there is a likely typo in Eq. (14): the Omega1 terms use n^q_0(k) where n^q_1(k) is presumably intended, and the appendix expressions are not fully derived. These are minor and fixable.\n\nOverall, the central qualitative claim is plausible and deserves a serious referee. The flat-DOS dependence should be tested before publication; if it survives, this is a solid subfield contribution. I would send it to review, with a request to address the DOS robustness and to soften the 'exact' language.","headline":"A solid extension of the HK model to finite-momentum pairing, with a plausible PDW-vs-SC phase diagram that hinges on an untested constant-DOS approximation.","tokens_in":17876,"tokens_out":2794,"would_cite":true,"duration_ms":23192,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Adding a pairing interaction to the exactly solvable Hatsugai-Kohmoto model produces a finite-temperature pair-density-wave instability at center-of-mass momentum $Q=(\\pi,\\pi)$ that beats uniform superconductivity at low doping and…","keywords":["Hatsugai-Kohmoto model","pair-density wave","doped Mott insulator","exactly solvable model","pair susceptibility","superconductivity","non-Fermi liquid","cuprate superconductors"],"falsifier":"Repeat the same exact susceptibility calculation using the actual square-lattice density of states (or a momentum-resolved numerical method) and check whether $\\chi_0(0,Q)$ still diverges at finite $T$ and still exceeds $\\chi_0(0,0)$ for $x\\lesssim 0.25$ and intermediate $U$; if it does not, the PDW-dominated phase diagram is an artifact of the constant-DOS approximation.","tokens_in":16885,"feed_emoji":"⚛️","tokens_out":9834,"duration_ms":78824,"temperature":0.7,"pith_summary":"The paper studies the Hatsugai-Kohmoto model on a square lattice—an exactly solvable model that is a Mott insulator at half-filling and a non-Fermi liquid when doped—augmented by a BCS-like pairing interaction. It claims that, at low doping up to a critical value around $x_c\\approx 0.25$ and for intermediate interaction strength $U$, the static pair susceptibility at center-of-mass momentum $Q=(\\pi,\\pi)$ diverges at a finite critical temperature, signaling strong pair-density-wave (PDW) fluctuations. In that doping window the PDW channel beats the ordinary $q=0$ superconducting channel; further doping favors uniform superconductivity instead. The interest is that the HK model's Mott transition has been argued to lie in the same universality class as the Hubbard model's, so this exact result offers a potential picture of PDW fluctuations on doping a Mott insulator, a scenario proposed for underdoped cuprates.","feed_headline":"Pair-density waves beat uniform superconductivity at low doping","feed_subtitle":"Exactly solvable Mott model puts (π,π) pairing ahead of uniform superconductivity at low doping.","key_machinery":"The central object is the static pair susceptibility $\\chi_0(0,q)$ of the HK model with a pairing interaction $V_q$, whose divergence condition $\\chi_0(0,q)=1/V_q$ locates the superconducting ($q=0$) and PDW ($q=Q=(\\pi,\\pi)$) instabilities. The machinery is the exact two-band structure of the HK Green's function, the constant density of states $\\rho(\\omega)=1/W$ (taken for simplicity), and the identity $\\xi_{k+Q}=-\\xi_k-2\\mu$, which makes the PDW pair-energy denominator independent of $k$ and renders the susceptibility integrals exact.","core_discovery":"Working with the exact Green's function of the HK model and a Dyson equation for the pair susceptibility, the authors show that the $Q=(\\pi,\\pi)$ pairing channel has a Cooper-like bound state with negative pair-binding energy $\\varepsilon'$ for low doping and intermediate $U$, and that the corresponding static susceptibility $\\chi_0(0,Q)$ diverges at a finite $T_c^{\\mathrm{PDW}}$ once the pairing strength $V'$ exceeds a small threshold. Because $\\xi_{k+Q}=-\\xi_k-2\\mu$, the two-particle energy denominators in the PDW susceptibility become $k$-independent, so the momentum integrals can be done exactly. The result is a phase diagram in which PDW fluctuations dominate over $q=0$ superconductivity in an underdoped, intermediate-$U$ region, while uniform SC takes over at higher doping; for $V'\\to 0$ the PDW susceptibility is enhanced but does not diverge at finite $T$, so the model realizes a fluctuating PDW regime rather than long-range PDW order in that limit.","pith_inferences":["A direct test would repeat the exact susceptibility integrals with the actual square-lattice density of states (including van Hove singularities) instead of $\\rho(\\omega)=1/W$; if the PDW-dominated window survives, the phase diagram becomes a quantitative benchmark for approximate methods in the Hubbard model.","Because the $Q=(\\pi,\\pi)$ PDW here is tied to $\\eta$-pairing, one could probe the fluctuating PDW regime through dynamical charge correlations or spectral functions of the HK model, looking for signatures that could be compared with finite-size numerics.","The exact $T_c^{\\mathrm{PDW}}$ formulas could be used as a controlled starting point to study how $k$-dependent pairing interactions, disorder, or the orbital-HK extension shift the PDW-versus-SC competition.","If the fluctuating-PDW regime is the relevant 'mother state' for underdoped cuprates, a testable next step is to check whether the model reproduces pseudogap-like spectral features, such as Fermi arcs, in the same doping and interaction window."],"forward_implications":["At finite pairing strength $V'$ and low doping up to $x_c\\approx 0.25$ with intermediate $U$, the static pair susceptibility at $Q=(\\pi,\\pi)$ diverges at a finite $T_c^{\\mathrm{PDW}}$, establishing a genuine finite-temperature PDW instability in the HK-SC model.","In this underdoped window the PDW channel dominates the $q=0$ SC channel; at higher doping the uniform superconducting instability becomes the leading one.","For $V'\\to 0$ the $Q=(\\pi,\\pi)$ susceptibility is strongly enhanced as $T$ drops but does not diverge at finite temperature, so the model describes a fluctuating PDW regime rather than true long-range PDW order.","Increasing $U$ up to the bandwidth strengthens PDW pair binding while weakening uniform SC, which is why intermediate correlation strength favors PDW; $U>W$ suppresses both phases.","If the HK model's Mott transition indeed shares a universality class with the Hubbard model, the exact PDW-fluctuation phase diagram becomes a candidate qualitative explanation for fluctuating-PDW behavior in doped Mott insulators such as underdoped cuprates."],"supporting_citations":[{"why":"Defines the Hatsugai-Kohmoto model, the exactly solvable momentum-space interacting system whose pairing instabilities this paper computes.","marker":"[37]"},{"why":"Supplies the Dyson-equation treatment of the pair susceptibility and the exact q=0 superconducting instability that the present work extends to the PDW channel.","marker":"[42]"},{"why":"Argues the HK metal-insulator transition lies in the same universality class as the Hubbard Mott transition, motivating the cuprate relevance of the PDW result.","marker":"[59]"},{"why":"Cooper's pair-binding argument underpins the variational calculation of the PDW pair-binding energy.","marker":"[62]"},{"why":"Provides the BCS exponential form for the q=0 superconducting critical temperature that the paper recovers.","marker":"[63]"},{"why":"Connects the Q=(π,π) PDW to η-pairing off-diagonal long-range order in the Hubbard model.","marker":"[65]"},{"why":"One of the recent results showing that the Q=(π,π) PDW susceptibility is a close competitor to SC in correlated models, guiding the choice of pairing momentum.","marker":"[34]"},{"why":"Supports the close-competitor status of the Q=(π,π) PDW and the fluctuating-PDW description of the pseudogap.","marker":"[64]"},{"why":"Defines the physics of pair-density waves and their proposed relevance to cuprate superconductors.","marker":"[1]"},{"why":"Argues PDW fluctuations are relevant for the cuprate pseudogap, the interpretation the paper borrows for its fluctuating PDW regime.","marker":"[8]"}],"fun_headline_variants":["Exact model shows pair-density waves win at low doping","PDW fluctuations outshine uniform superconductivity in Mott model","Solvable doped Mott insulator reveals strong PDW pairing","π-π pairing dominates in exactly solvable Mott insulator"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes a flat, featureless density of states in place of the real square-lattice one, and the values of the pair-binding energies, critical temperatures, and the critical doping $x_c\\approx 0.25$ depend on that choice; with the true band structure the PDW-versus-SC competition could change.","fun_headline_variants_meta":{"raw":{"variants":["Exact model shows pair-density waves win at low doping","PDW fluctuations outshine uniform superconductivity in Mott model","Solvable doped Mott insulator reveals strong PDW pairing","π-π pairing dominates in exactly solvable Mott insulator"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000553,"raw_usage":{"total_tokens":2654,"prompt_tokens":979,"completion_tokens":1675,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":1609}},"tokens_in":595,"tokens_out":1675,"duration_ms":11480,"temperature":1.0,"reasoning_tokens":1609,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:23:08.873028+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the same exact susceptibility calculation using the actual square-lattice density of states (or a momentum-resolved numerical method) and check whether $\\chi_0(0,Q)$ still diverges at finite $T$ and still exceeds $\\chi_0(0,0)$ for $x\\lesssim 0.25$ and intermediate $U$; if it does not, the PDW-dominated phase diagram is an artifact of the constant-DOS approximation.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Argues the HK metal-insulator transition lies in the same universality class as the Hubbard Mott transition, motivating the cuprate relevance of the PDW result."},{"cited_title":"Liu, X.-X","cited_arxiv_id":null,"evidence_quote":"One of the recent results showing that the Q=(π,π) PDW susceptibility is a close competitor to SC in correlated models, guiding the choice of pairing momentum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the physics of pair-density waves and their proposed relevance to cuprate superconductors."}],"review_version":1}