{"id":"1d4a27e5-d52c-4bec-b396-59d0f05dd635","arxiv_id":"2502.04739","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"In a tight-binding extended Hubbard model, mixed singlet-triplet superconducting phases emerge below Tc, producing distinctive multi-peak density-of-states signatures that change with temperature.","lead":"This paper calculates how a model superconductor can switch its pairing symmetry as the temperature drops, and shows that the energy gap measured at low temperatures can look very different from the gap just below the transition temperature. The result warns experimentalists that measuring a superconducting gap at the lowest temperature may misidentify the mechanism that causes pairing.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unitarity constraint (Eq. 12) excludes non-unitary mixed phases; if s+d+p were allowed, the phase diagram and the claimed DOS fingerprint could change.","rationale":"The reader's weakest assumption is the unitarity condition Eq. (12), and I agree this is the most load-bearing concern. The paper is transparent about this restriction, but the central claims of prevalence and a distinctive spectroscopic fingerprint are tied to it. The proposed numerical test directly addresses whether allowing non-unitary mixed phases changes the phase diagram and DOS. The reader's CONDITIONAL verdict already reflects this limitation, so no verdict adjustment is needed. I do not see a more serious internal error: the self-consistent equations, free-energy comparison, and DOS calculations are internally consistent within the unitary mean-field framework, and the paper explicitly flags the main caveat.","tokens_in":17581,"tokens_out":17764,"duration_ms":187483,"concrete_test":"Recompute the T=0 phase diagram and the representative cooling curves in Figs. 8–10 without imposing Eq. (12): solve the full 2x2 BCS gap equations for complex Δ_s and Δ_t with arbitrary relative phase, evaluate the free energy for all stationary points, and carry the global minimum into the DOS calculation. Specifically, at the □ point of Fig. 7(a) and the 9 point of Fig. 7(c), test whether a non-unitary s+d+p solution has lower free energy than s+d+ip and whether the DOS peak count or peak positions change. If s+d+ip remains the global minimum and its multi-peak DOS persists, the concern is retired; if not, the central claim is conditional on the unitarity assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that mixed-symmetry s+d+ip phases are prevalent and show a distinctive multi-peak DOS rests on the unitarity assumption in Eq. (12), which forces Re[(Δ_s)^* Δ_t] = 0 (Eq. 13) and thus a π/2 phase difference between singlet and triplet components. The paper explicitly excludes non-unitary phases such as s+d+p, citing that most superconductors have unitary gaps, but this is a modeling choice rather than a derivation. Because the mean-field free energy in Eq. (14) and the quasiparticle energy E_k are written in the unitary form |Δ_k|^2 = |Δ_s|^2 + |Δ_t|^2, the calculation never competes non-unitary candidates against unitary ones; the prevalence of s+d+ip over s+d+p is therefore imposed, not derived. A non-unitary mixed phase would have a two-branch quasiparticle spectrum and generically a different DOS, so the claimed multi-peak fingerprint could be an artifact of the unitary restriction. This does not make the paper internally inconsistent, but it does limit the generality of the abstract's 'prevalence' and 'complete picture' claims.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the extended Hubbard model with on-site U and nearest-neighbor V in one and two dimensions, using BCS mean-field theory to solve self-consistently for superconducting order parameters of s-, p-, and d-wave symmetry. The authors map Tc and T=0 phase diagrams for several fillings, showing that mixed-symmetry phases (s+ip in 1D; s+id and s+d+ip in 2D) become stable below Tc even when a pure-symmetry phase wins at Tc. They compute the density of states as a function of temperature and identify the origins of the peaks (Van Hove and BCS coherence peaks). The central result is that the 2D s+d+ip phase shows a distinctive multi-peak DOS, which the authors propose as a spectroscopic fingerprint for mixed-symmetry superconductivity. The paper emphasizes that gap symmetry can change with temperature, so low-temperature measurements alone may misidentify the pairing mechanism.","tokens_in":17805,"tokens_out":10770,"duration_ms":102869,"significance":"The paper is a careful and clearly written mean-field study that extends earlier work on the extended Hubbard model. Its strengths include full self-consistent solution of the coupled gap equations, phase diagrams over a wide parameter range, and analytic expressions for the DOS peak positions (Eqs. 20, 22-27) that allow the peak structure to be traced to specific features of the band structure and gap. The identification of a multi-peak DOS unique to the s+d+ip phase, if robust, provides a concrete experimental signature for mixed-symmetry superconductivity in quasi-2D materials. However, the analysis is deliberately restricted to unitary gaps, which excludes non-unitary mixed phases such as s+d+p; the paper does not assess whether such phases would be competitive. The 'complete picture' claim in the abstract is therefore stronger than what is demonstrated. Overall, this is a useful contribution that would benefit from a more cautious presentation of its generality.","major_comments":[{"comment":"The unitarity condition (Eq. 12) forces any coexisting singlet and triplet components to have a π/2 phase difference (Eq. 13), so non-unitary mixed phases such as s+d+p are excluded from the calculation. Because the quasiparticle energy Ek and the free energy (Eq. 14) are written in the unitary form, the self-consistent search never compares non-unitary candidates against unitary ones. The phase diagrams in Fig. 10 and the claim that s+d+ip is 'the most stable mixed-symmetry phase' (Discussion, Section V) are therefore conditional on this assumption. While the paper acknowledges the restriction, it does not justify that non-unitary phases would not appear in some parameter regions, nor does it discuss how the DOS fingerprint would change for a non-unitary two-branch spectrum. This is load-bearing for the abstract's 'prevalence' and 'complete picture' claims. I recommend either (a) softening the generality claims and explicitly labelling the results as valid for unitary gaps, or (b) adding a representative calculation for a non-unitary phase (e.g., s+d+p) to show that it does not alter the phase diagram or the DOS fingerprint.","section":"Section II, Eqs. (12)-(13); Section IV.B/D"},{"comment":"The phase diagrams are computed on a finite U-V grid, but the paper does not state the grid spacing, the convergence tolerance for the self-consistent equations, or the criterion used to identify the free-energy minimum when multiple stationary points coexist. Transition temperatures such as Tc1 = 0.803t/kB and Tc2 = 0.525t/kB (Fig. 8) are quoted to three significant figures without uncertainty. Because a central claim is that the s+d+ip phase occupies the largest region of parameter space at ne = 0.75, the lack of resolution information makes it difficult to assess the robustness of the phase boundaries. Please provide the numerical parameters (step sizes, tolerances, number of k-points) and estimate the resulting uncertainty in the phase boundaries.","section":"Section IV.C, Figs. 7 and 10"}],"minor_comments":[{"comment":"The U,V values for the ne = 0.50 DOS calculation are given as U/t = -0.50, V/t = -3.00 in the text, but the Fig. 12 caption states U/t = -2.50, V/t = -3.60. Please correct this inconsistency.","section":"Section IV.D and Fig. 12 caption"},{"comment":"The Gaussian broadening ν/t = 0.03 is arbitrary and, as the authors note, causes some peaks to merge (e.g., Fig. 13(c)). A brief justification of the chosen width or a discussion of the sensitivity of the multi-peak fingerprint to ν would strengthen the analysis.","section":"Section III.D and Appendix A"},{"comment":"The statement 'A complete picture of the superconducting symmetry can only be attained if measurements are made over the entire temperature range' is stronger than what the model shows; the model shows that symmetry changes can occur, but not that they are ubiquitous. Consider rewording to 'may need to be'.","section":"Abstract and Section I"},{"comment":"The claim that the extended s-wave and p-wave order parameters appear at the same temperature in the d → s+d+ip transition is surprising and could benefit from a brief explanation or a reference to a Ginzburg-Landau analysis.","section":"Section IV.B"}],"recommendation":"major_revision","confidential_remarks":"The paper is a competent mean-field study, but the unitarity assumption is a significant limitation that should be clearly acknowledged. The authors are transparent about it, but the abstract overstates the generality. If they are willing to soften the claims and fix the numerical and consistency issues, a revised version would be suitable for publication. I do not see any fatal flaws."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know before you read it. First, this is a systematic BCS mean-field study of the extended Hubbard model with coexisting singlet and triplet order parameters, and it shows that mixed-symmetry phases appear generically below Tc. Second, the genuinely new content is the T=0 phase diagram in 2D with five order parameters and the multi-peak DOS structure in the s+d+ip phase, which they argue is a spectroscopic fingerprint.\n\nWhat it does well: the calculation is transparent, with the self-consistent equations written out and the phase diagrams carefully presented. The DOS analysis is honest about the origin of each peak (Van Hove versus BCS coherence) and the appendices give the full set of plots. They do not overclaim that this identifies a real material; they position it as a caution for experimental interpretation. Prior work, especially Hutchinson and Marsiglio (Ref. [34]), is properly acknowledged.\n\nSoft spots, in proportion. The unitarity assumption (Eq. 12) genuinely restricts the set of allowed phases: non-unitary mixed phases like s+d+p are excluded by construction, so the abstract's 'prevalence' claim is really 'prevalence among unitary gaps'. The paper says this explicitly, but the wording of the abstract glosses over it. That is a scope limitation, not an internal error. The Gaussian broadening is standard, and the phase boundaries carry no error bars, but for this kind of numerical mean-field study that is normal. The DOS 'fingerprint' is model-dependent; no argument is given that it survives beyond mean-field or in a more realistic band structure. The final advice about measuring over the whole temperature range is sensible, but the word 'only' in the abstract is too strong, since phase-sensitive probes can sometimes determine symmetry at a single temperature.\n\nWho this is for: anyone working on unconventional pairing or interpreting tunneling spectra in quasi-2D superconductors. It is a useful reminder that the gap symmetry measured at low T need not be the symmetry selected at Tc. It deserves a serious referee.\n\nRecommendation: send it to peer review. Ask the authors to make the unitarity restriction more prominent in the abstract or introduction, and to soften the 'only' in the final sentence, but do not reject over either.","headline":"A solid mean-field mapping of mixed-singlet-triplet phases below Tc in the extended Hubbard model, with a real but model-dependent DOS fingerprint; the unitarity restriction is a clear scope limit, not a hidden flaw.","tokens_in":18328,"tokens_out":1989,"would_cite":true,"duration_ms":22128,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.20.Fg","74.25.Jb"],"model":"deepseek-v4-flash","headline":"In a generic tight-binding superconductor, pure singlet or triplet phases give way to mixed-symmetry phases on cooling, and the mixed s+d+ip phase leaves a multi-peak tunneling fingerprint.","keywords":["mixed-symmetry superconductivity","singlet-triplet pairing","extended Hubbard model","density of states","Van Hove singularity","BCS mean-field theory","phase diagram","time-reversal symmetry breaking"],"falsifier":"Take a candidate $s{+}d{+}ip$ superconductor and measure its tunneling density of states from just below $T_c$ down to $T=0$: if the predicted split peaks (Eqs. 23-26) do not appear, or if a non-unitary phase like $s{+}d{+}p$ is stabilized instead, the central claim fails. A simpler numerical check is to relax the unitarity condition in the self-consistent equations and see whether the $s{+}d{+}p$ phase occupies any substantial region of the phase diagram.","tokens_in":17404,"feed_emoji":"🔬","tokens_out":6946,"duration_ms":65559,"temperature":0.7,"pith_summary":"This paper argues that in a translationally invariant extended Hubbard model, superconducting phases with mixed spin-singlet and spin-triplet symmetry emerge naturally below the critical temperature, not just pure s-, p-, or d-wave phases. Using BCS mean-field theory, it maps phase diagrams in one and two dimensions and follows specific coupling points as temperature drops, finding multiple superconducting transitions in which a pure phase stable just below Tc gives way to a mixed-symmetry phase at low temperature. The central new result is that in 2D the s+d+ip phase shows a multi-peak density-of-states structure created by splitting the Van Hove and BCS coherence peaks, a feature absent in the pure phases. A sympathetic reader would care because the symmetry measured at low temperature need not be the symmetry that emerged at Tc, so a complete picture of the superconducting order requires measurements over the entire temperature range.","feed_headline":"Mixed-symmetry superconductivity leaves multi-peak tunneling spectra","feed_subtitle":"In a 2D Hubbard model, pure d-wave order transforms into s+d+ip below Tc and splits the DOS peaks.","key_machinery":"The central object is the mean-field gap matrix of the extended Hubbard model, whose singlet and triplet parts are expanded in symmetry-allowed basis functions; in 2D, $\\Delta_k^{(s)}=\\Delta_0+\\Delta_{s^*}s_k+\\Delta_{d_{x^2-y^2}}d_k$ and $\\Delta_k^{(t)}=\\Delta_{p_x}\\sin k_x+\\Delta_{p_y}\\sin k_y$. The argument is carried by the unitarity condition $\\Delta_k^\\dagger\\Delta_k=|\\Delta_k|^2\\mathbf{1}$, which turns the quasiparticle energy into a quadrature sum of singlet and triplet magnitudes and forces a $\\pi/2$ phase difference between them when they coexist. This is what makes a phase such as $s{+}d{+}ip$ possible while excluding $s{+}d{+}p$. The density-of-states fingerprint then follows from saddle points of the quasiparticle dispersion, with peak positions given analytically in terms of the competing order parameters.","core_discovery":"On the paper's own terms, the discovery is that coexistence of singlet and triplet order parameters in the s+d+ip phase of the two-dimensional extended Hubbard model produces a characteristic multi-peak density of states. Because the unitarity condition forces the singlet and triplet components to differ in phase by pi/2, the quasiparticle energy $E_k=(\\xi_k^2+|\\Delta_k^{(s)}|^2+|\\Delta_k^{(t)}|^2)^{1/2}$ contains both components in quadrature, and the saddle points of $E_k$ along the $k_x=0$ and $k_y=0$ lines no longer coincide. The result is that the single logarithmic Van Hove or BCS coherence peaks of the pure d-wave phase split into two, three, or four peaks, with locations set by the s-wave and d-wave amplitudes through Eqs. (23)-(26). At half filling the mixed phase is fully gapped; below half filling the peak structure is two, three, or occasionally four peaks per frequency band. The paper also establishes that several ground-state phases in both 1D and 2D are mixed s+ip or s+d+ip rather than pure, and that these mixed phases spontaneously break time-reversal symmetry.","pith_inferences":["If non-unitary gaps were admitted, phases like $s{+}d{+}p$ (without the imaginary $i$) might appear, and the multi-peak fingerprint would need revision; this is a testable extension of the model.","Because mixed phases carry a $\\pi/2$ phase difference, they break time-reversal symmetry, suggesting one could look for accompanying spontaneous signatures such as polar Kerr rotation or chiral edge currents.","The paper's peak-count analysis could be applied to quasi-2D organic superconductors where three-peak tunneling spectra have been reported, checking whether the peak positions follow Eqs. (23)-(26).","The formalism naturally extends to longer-range interactions, which would add order parameters in other irreducible representations of the square-lattice point group and potentially produce additional peak structures."],"forward_implications":["A pure-symmetry phase identified just below $T_c$ can be a different superconducting phase at low temperature, so symmetry identification requires measurements over the full superconducting temperature range.","The multi-peak density-of-states structure, unique to the $s{+}d{+}ip$ phase among the phases studied, gives a spectroscopic fingerprint for mixed singlet-triplet pairing.","The $s{+}id$ phase is fully gapped and shows no peak splitting, so a full gap with a single coherence peak does not rule out mixed symmetry.","At $n_e=0.75$ the $s{+}d{+}ip$ phase occupies the largest region of parameter space, making mixed phases a robust outcome rather than a fine-tuned one.","In optical-lattice quantum simulators of the Hubbard model, the predicted temperature-dependent density of states could be measured directly."],"supporting_citations":[{"why":"Supplies the mean-field formalism and prior mixed-symmetry solutions of the extended Hubbard model that this paper extends to full phase diagrams and densities of states.","marker":"[34]"},{"why":"Provides earlier exotic superconducting states in the extended attractive Hubbard model that motivate the mixed-symmetry search.","marker":"[33]"},{"why":"Gives the Ginzburg-Landau understanding of the pi/2 phase difference between coexisting order parameters in tight-binding models.","marker":"[32]"},{"why":"An earlier concrete calculation showing mixed-symmetry order parameters can appear below Tc.","marker":"[31]"},{"why":"Provides the analytic density-of-states saddle-point formulas for d-wave superconductors used to locate the Van Hove and coherence peaks.","marker":"[45]"},{"why":"Reports experimental evidence for mixed-symmetry superconductivity with multi-peak tunneling spectra in an organic metal, the comparison system for the predicted fingerprint.","marker":"[46]"},{"why":"Supplies the context on non-unitary superconductivity used to justify restricting the study to unitary gap parameters.","marker":"[42]"},{"why":"Establishes the phenomenological framework in which different order-parameter symmetries can mix below Tc.","marker":"[17]"},{"why":"Lays out the symmetry classification of superconducting order parameters on which the irreps analysis relies.","marker":"[16]"}],"fun_headline_variants":["Mixed singlet-triplet gaps split DOS peaks","Superconducting symmetry shifts with temperature","s+d+ip order leaves multi-peak fingerprints","Temperature reveals hidden superconducting symmetry","Mixed parity gaps reshape density of states"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the superconducting gap is unitary, so any coexisting singlet and triplet components must differ in phase by exactly $\\pi/2$; lift that constraint and non-unitary phases such as $s{+}d{+}p$ could change both the phase diagram and the density-of-states peaks.","fun_headline_variants_meta":{"raw":{"variants":["Mixed singlet-triplet gaps split DOS peaks","Superconducting symmetry shifts with temperature","s+d+ip order leaves multi-peak fingerprints","Temperature reveals hidden superconducting symmetry","Mixed parity gaps reshape density of states"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000502,"raw_usage":{"total_tokens":2453,"prompt_tokens":945,"completion_tokens":1508,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":561,"completion_tokens_details":{"reasoning_tokens":1444}},"tokens_in":561,"tokens_out":1508,"duration_ms":11314,"temperature":1.0,"reasoning_tokens":1444,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T21:39:30.847694+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a candidate $s{+}d{+}ip$ superconductor and measure its tunneling density of states from just below $T_c$ down to $T=0$: if the predicted split peaks (Eqs. 23-26) do not appear, or if a non-unitary phase like $s{+}d{+}p$ is stabilized instead, the central claim fails. A simpler numerical check is to relax the unitarity condition in the self-consistent equations and see whether the $s{+}d{+}p$ phase occupies any substantial region of the phase diagram.","supporting_citations":[{"cited_title":"Nayak and S","cited_arxiv_id":null,"evidence_quote":"Supplies the mean-field formalism and prior mixed-symmetry solutions of the extended Hubbard model that this paper extends to full phase diagrams and densities of states."},{"cited_title":"Kuboki, Effect of Band Structure on the Symmetry of Superconducting States, Journal of the Physical Society of Japan 70, 2698–2702 (2001)","cited_arxiv_id":null,"evidence_quote":"Provides earlier exotic superconducting states in the extended attractive Hubbard model that motivate the mixed-symmetry search."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Ginzburg-Landau understanding of the pi/2 phase difference between coexisting order parameters in tight-binding models."},{"cited_title":"Kheirkhah, Z","cited_arxiv_id":null,"evidence_quote":"An earlier concrete calculation showing mixed-symmetry order parameters can appear below Tc."},{"cited_title":"Zhou and H","cited_arxiv_id":null,"evidence_quote":"Reports experimental evidence for mixed-symmetry superconductivity with multi-peak tunneling spectra in an organic metal, the comparison system for the predicted fingerprint."},{"cited_title":"Hashimoto, I","cited_arxiv_id":null,"evidence_quote":"Supplies the context on non-unitary superconductivity used to justify restricting the study to unitary gap parameters."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the phenomenological framework in which different order-parameter symmetries can mix below Tc."},{"cited_title":"On the Superfluidity of Liquid He3","cited_arxiv_id":null,"evidence_quote":"Lays out the symmetry classification of superconducting order parameters on which the irreps analysis relies."}],"review_version":1}