{"id":"7acfab35-fe8b-46a5-be12-7ca59d25d3fd","arxiv_id":"1908.02514","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Spin-nematic (quadrupolar) order turns the fermion-quadrupolar interaction into a four-fermion pairing interaction that yields p-wave superconductivity; the sign of the coupling fixes equal-spin versus opposite-spin pairing, and the critical coupling vanishes at zero chemical potential.","lead":"This physics paper proposes that a special spin ordering called spin-nematic order, in which spins align like tiny rods without creating any magnet, can by itself make electrons pair up and conduct without resistance. The model predicts a pairing pattern set by the sign of the interaction and offers a candidate explanation for superconductivity in the iron compound FeSe.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (25)-(27) imply κcr=0 for all μ inside the band, so the finite quantum-critical coupling in Fig. 2 is likely a numerical artifact; the central transition claim fails.","rationale":"The paper's goal is to show static quadrupolar order alone produces p-wave pairing with a μ-dependent critical coupling. The algebraic step from the classical nematic background to the four-fermion interaction (18) is plausible and the mean-field decoupling is standard; I see no error in the sign selection of the parallel vs antiparallel channel. However, the central quantitative result—a quantum critical coupling that grows with μ—appears to be an artifact. The gap equation (26) is the standard BCS equation; its linearization (27) contains ∫ e_k^2/|ξ_k|. For the cubic tight-binding dispersion defined in Eq. (19), any μ inside the band has a Fermi surface. The T1u form factor e_k = sinkx+sinky+sinkz is not identically zero on that surface, so the integral diverges logarithmically. Thus an arbitrarily small |κ| yields pairing; there is no finite critical coupling. The figure's finite κcr must arise either from a wrong dispersion ('-2t ε_k - μ' double-counts t) or from a coarse-mesh regularization of a divergent integral. This directly contradicts the abstract's 'quantum critical transition' and the conclusion's 'critical value depends on the chemical potential.' The reader's conditional verdict focused on the rigid-nematic background and missing numerics; while related to numerical divergence, my concern identifies a false positive in the central phase diagram. If the divergence check confirms κcr = 0 for in-band μ, the paper's main claim cannot stand as written, though the pairing mechanism itself (static nematic order → p-wave superconductivity) may still be salvageable. For this reason I recommend REJECT rather than CONDITIONAL.","tokens_in":7122,"tokens_out":18103,"duration_ms":172416,"concrete_test":"Evaluate I(μ) = ∫_{BZ} d^3k/(2π)^3 (sinkx+sinky+sinkz)^2 / | -2t(coskx+cosky+coskz) - μ | for μ/t = 4 (and μ/t = 0) on a simple cubic lattice with momentum grids N^3 for N = 32, 64, 128, 256. If I(μ) grows with N (logarithmically) so that 1/I → 0, then κcr = 0 for in-band μ and Fig. 2's finite critical coupling is a grid artifact. If instead I(μ) saturates at a finite value, re-check the derivation: re-derive Eq. (25) from the free energy Eq. (22) using ξ_k = ε_k - μ, and confirm whether the denominator should be |ξ_k|; the text's '-2t ε_k - μ' is not consistent with the ε_k defined in Eq. (19).","verdict_should_be":"REJECT","load_bearing_attack":"The most load-bearing problem is not the rigid-nematic assumption but the gap equation itself. Eq. (19) defines ε_k = -2t(coskx+cosky+coskz), so the single-particle energy is ξ_k = ε_k - μ. Eqs. (25) and (27) instead use (-2t ε_k - μ), which either double-counts the hopping amplitude or is dimensionally inconsistent. With the correct ξ_k, the linearized gap equation is 1 = |κ| ∫_{BZ} d^3k/(2π)^3 (sinkx+sinky+sinkz)^2 / |ξ_k|. For any μ inside the tight-binding band, the Fermi surface has nonzero measure and the T1u form factor does not vanish identically on it, so the integral diverges logarithmically; Cooper's theorem then gives κcr = 0 for every in-band μ, not only μ = 0. The finite κcr that grows with μ shown in Fig. 2 is therefore either computed with the wrong dispersion or is a numerical artifact of an unregularized divergent integral on a finite grid. This removes the advertised quantum critical transition from normal spin-nematic to coexisting superconducting-nematic phase and the μ-dependent critical coupling in the abstract and Conclusion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a spin-fermion model on a cubic lattice, combining a spin-1 Heisenberg/biquadratic Hamiltonian, a tight-binding fermion band, and a fermion-quadrupolar interaction. In the ferroquadrupolar/nematic phase the quadrupolar operator is replaced by its classical expectation value, reducing the interaction to a four-fermion pairing interaction. A Hartree-Fock treatment with a T1u gap ansatz yields spin-triplet p-wave superconductivity, with the sign of the coupling selecting equal-spin versus opposite-spin pairing. The paper claims that at zero chemical potential an arbitrarily small coupling produces superconductivity, while for nonzero chemical potential there is a quantum critical coupling separating the normal nematic phase from a phase where superconductivity coexists with nematicity, and it connects the low-pressure regime of FeSe to this scenario.","tokens_in":7424,"tokens_out":9259,"duration_ms":100495,"significance":"The proposed mechanism—static quadrupolar order generating odd-parity pairing—is conceptually interesting and the paper makes a specific, falsifiable prediction for the chemical-potential dependence of the critical coupling. The model is simple and the T1u channel selection is clearly identified, and the paper does not fit any parameter to data. However, the central quantitative prediction is not reliable because the gap equations use a single-particle dispersion inconsistent with the Hamiltonian, and the advertised finite critical coupling for in-band chemical potentials appears to be an artifact of that inconsistency.","major_comments":[{"comment":"The single-particle dispersion in the gap equations is inconsistent with the Hamiltonian. Equation (19) defines ε_k = -2t(coskx+cosky+coskz), so the quasiparticle energy should contain ξ_k = ε_k - μ. Equation (23) omits μ entirely, and Eqs. (25)-(27) use (-2t ε_k - μ), which is neither ξ_k nor any dimensionally consistent expression. All subsequent results, including Figs. 1 and 2, need to be recomputed with the correct dispersion.","section":"Superconductivity, Eqs. (19)-(27)"},{"comment":"With the correct dispersion, the linearized gap equation is 1 = |κ| ∫ d^3k/(2π)^3 (sin kx + sin ky + sin kz)^2 / |ε_k - μ|. For any μ inside the band (-6t < μ < 6t), the Fermi surface has nonzero measure, and the integrand diverges logarithmically near it, giving κcr = 0 rather than a finite value. The finite κcr(μ) shown in Fig. 2, and the quantum critical transition claimed in the abstract and Conclusion, are therefore not supported by the corrected equation. This is a load-bearing error in the central claim.","section":"Eq. (27) and Fig. 2"},{"comment":"The reduction of the fermion-quadrupolar interaction to the four-fermion interaction (18) is not shown; since this is the basis of the pairing channel, the algebra should be provided. More importantly, the replacement of S^ba by its classical expectation value in Eq. (17) is made once and never updated: the fermions have no feedback on the nematic order, and dynamical quadrupolar fluctuations are discarded. Given that the paper claims coexistence of nematicity and superconductivity and motivates the model with FeSe, where nematic fluctuations are strong, the rigidity of the nematic background is a significant limitation that should be justified or relaxed.","section":"Fermion-Quadrupolar Interaction, Eqs. (16)-(18)"},{"comment":"No numerical grid, momentum cutoff, or regularization is specified for the integrals in Eqs. (25)-(27). Because the κcr integral is divergent at the Fermi surface for in-band μ, any finite value obtained from it is regularization-dependent, and the results in Figs. 1 and 2 cannot be reproduced as presented.","section":"Figs. 1 and 2"}],"minor_comments":[{"comment":"The statement that the system (25) has no solution with both Δp and Δap nonzero is asserted without proof; either provide a derivation or a reference.","section":"After Eq. (24)"},{"comment":"Reference [5] is misformatted (\"Aarts Steglich\" should be Steglich et al.), and several other references contain typographical or formatting errors.","section":"References"},{"comment":"The caption is ambiguous: \"upper graph/lower one\" combined with two colors makes it difficult to identify which curve corresponds to which sign of the coupling and which value of μ; please label the curves directly.","section":"Fig. 1 caption"},{"comment":"The identity <c_{k↑}c_{-k↑}> = <c_{k↓}c_{-k↓}> is used without justification; in a nematic background that breaks spin-rotation symmetry, this equality should be verified explicitly.","section":"Eq. (20)"},{"comment":"The statement that the low-pressure region of FeSe is \"well described by the theory\" is qualitative only; no direct comparison with the measured Tc or pressure dependence is provided.","section":"Conclusion"}],"recommendation":"reject","confidential_remarks":"The manuscript contains a promising idea, but the central quantitative claim is invalidated by the dispersion error in the gap equations; the corrected equations remove the finite κcr for in-band μ, so the advertised quantum critical transition is an artifact. I recommend rejection rather than major revision because the central scenario, as stated, cannot survive the corrected gap equation. The report gives the authors the concrete technical path if they wish to reformulate the model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nTwo things to know. The paper's core construction is new and worth a look: it builds a fermion-quadrupolar interaction out of the triplet-pairing d-vector, then shows that in a nematic background this reduces to a four-fermion pairing term with a sign-selected spin structure. The sign rule—κ>0 equal-spin, κ<0 opposite-spin—is concrete and testable. That part is good. But the paper's central phase diagram claim, the finite μ-dependent critical coupling, is wrong. It comes from a dispersion error.\n\nLook at Eqs. (25)-(27). The paper defines ε_k = -2t(coskx+cosky+coskz), so the single-particle energy is ξ_k = ε_k - μ. The gap equations instead use (-2t ε_k - μ). That double-counts t and is dimensionally off. With the correct ξ_k, the linearized gap equation is 1 = |κ| ∫ e_k^2 / |ξ_k|. For any μ inside the band, the Fermi surface has nonzero measure and e_k = sinkx+sinky+sinkz does not vanish identically on it, so the integral diverges logarithmically. Cooper's theorem then gives κcr=0 for every in-band μ, not just μ=0. The finite κcr that grows with μ in Fig. 2 is almost certainly a numerical artifact of an unregularized divergent integral or the wrong dispersion. That removes the advertised quantum critical transition from the normal nematic to the coexisting superconducting-nematic phase.\n\nThe reader's other concerns are real but secondary. The no-coexistence claim after Eq. (25) is asserted, not proved. The numerics lack mesh and regulator details. And treating the nematic order as a rigid background ignores fluctuations, which is a strong assumption for FeSe. But those could be fixed or tempered. The dispersion error is load-bearing and needs a full re-derivation.\n\nGive the paper credit for not fitting parameters: κcr(μ) is computed, not tuned. And the reduction from (16) to (18) is genuinely clever. But as written, the main quantitative result does not survive.\n\nWho is this for? Someone interested in nematicity-induced pairing mechanisms could mine the construction and the sign rule. The paper deserves a serious referee, not because the conclusions are right, but because the idea is novel and the error is identifiable and fixable. I'd send it out, with a directive to redo the gap equations with the correct dispersion and either prove or drop the coexistence claim.\n\nWould I cite it? Not as it stands. Would I bring it to reading group? Maybe, as a cautionary tale about dispersion mistakes.\n\nBest,\n[Name]","headline":"Novel mechanism, but the central κcr(μ) result is an artifact of a wrong single-particle dispersion; the idea deserves peer review, the phase diagram does not.","tokens_in":7974,"tokens_out":6224,"would_cite":false,"duration_ms":65605,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":null,"created_at":"2026-08-14T14:43:31.196637+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":null,"supporting_citations":[],"review_version":1}