{"id":"075a9bf3-1dc2-4410-bccb-c6f51be541b6","arxiv_id":"1908.09476","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In four-band superconductors, unconventional s-wave pairing states survive nonmagnetic disorder much better than single-band theory predicts, with resilience set by the superconducting fitness.","lead":"This paper shows that unconventional s-wave superconducting states in multiband materials are far more resistant to impurities than previously expected. The resilience is controlled by the superconducting fitness, a quantity that can be computed for real candidate superconductors.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unresolved disagreement with Ref. [26] on which Hamiltonian enters the superconducting fitness directly affects the central effective-rate formula (Eq. 24).","rationale":"I read the paper as a serious attempt to generalize disorder pair-breaking to multiband systems with four internal degrees of freedom, and the internal algebra up to Eq. (24) is coherent under the stated assumptions. The reader's CONDITIONAL verdict is appropriate, and I do not propose changing it. However, the single most load-bearing unresolved point is not the isotropic-potential or well-separated-band approximation per se, but the explicit, self-reported disagreement with Ref. [26] over which Hamiltonian defines the superconducting fitness. Since Eq. (24) is the central diagnostic and the authors acknowledge that a concurrent analysis predicts a qualitatively different result (complete insensitivity versus parametric enhancement), the correctness of the central claim is not settled. The reader mentioned this disagreement in the rationale but did not make it the weakest assumption; my concern partially agrees with the reader's assessment. A numerical SCBA calculation for a concrete model, retaining all band terms, would settle the dispute directly and is the natural next check.","tokens_in":11022,"tokens_out":9291,"duration_ms":106607,"concrete_test":"Numerically solve the SCBA gap equation (Eqs. 8–9) for the CuxBi2Se3 model of Eq. (20) without the well-separated-band approximation, using the isotropic impurity potential of Eq. (4) in the linearized regime. Extract the effective scattering rate tau_nu from the suppression of T_c relative to the Abrikosov-Gor'kov form (Eq. 15). Then compare the extracted tau_nu with Eq. (21)/(24) and with the alternative fitness (relative to the impurity Hamiltonian) proposed in Ref. [26]. If the numerical result matches Eq. (24), the disagreement is resolved in favor of the authors; if it matches Ref. [26] or neither, the central claim requires revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—Eq. (24), 1/tau_nu = 1/tau - (1/tau_0)(1 - bar F_C)—identifies disorder robustness with the Fermi-surface-averaged superconducting fitness computed from the normal-state Hamiltonian H_k. The paper's own final section reports that a concurrent independent analysis (Ref. [26]) obtains an effective scattering rate involving the fitness with respect to the impurity Hamiltonian, not the normal-state Hamiltonian, and hence predicts complete insensitivity for states that are fit relative to the impurity potential, rather than the parametric enhancement claimed here. This is not a cosmetic difference: the lambda_l assignments in Tables I and II and the diagnostic bar F_C in Eq. (24) follow from commuting the pairing potential with the specific gamma_l matrices of H_k, whereas the alternative calculation uses a different operator. Both derivations start from the same SCBA framework (Eqs. 8–10), so the source of the discrepancy must be identified before Eq. (24) can be regarded as established. No internal algebraic error was found in the present manuscript, but the unresolved state of the dispute makes the central fitness-based robustness formula and its material-selection diagnostic load-bearing and unverified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies disorder effects on multiband superconductors with four internal degrees of freedom. Using the self-consistent Born approximation, the authors derive effective impurity scattering rates for unconventional s-wave pairing states [Eqs. (16) and (21)] and show that these rates are reduced relative to the normal-state rate. They express the reduction in terms of the Fermi-surface-averaged superconducting fitness [Eq. (24)], apply the formalism to YPtBi and CuxBi2Se3, and conclude that nontrivial s-wave channels are parametrically more robust against nonmagnetic disorder than single-band unconventional states. The paper also shows that this protection extends to other pairing states in the same irreducible representation, with the s-wave state providing an upper bound on the stability.","tokens_in":11185,"tokens_out":18894,"duration_ms":184128,"significance":"If Eq. (24) is established, the paper provides a simple, broadly applicable diagnostic for identifying disorder-robust unconventional pairing channels in multiband systems, connecting the concept of superconducting fitness to the observable T_c suppression. The applications to YPtBi and CuxBi2Se3, including tabulated λ_l assignments and numerical T_c curves, make the proposal concrete. The authors are careful to state their principal assumptions (isotropic impurity potential, well-separated bands) and to flag the disagreement with Ref. [26]. The clean derivation structure and explicit model-specific tables are definite strengths. However, because the disagreement with Ref. [26] directly concerns the content of the central formula, the paper's significance is currently conditional on resolving that discrepancy.","major_comments":[{"comment":"The manuscript explicitly reports that a concurrent independent calculation (Ref. [26]) obtains an effective scattering rate based on the superconducting fitness with respect to the impurity Hamiltonian, yielding complete insensitivity for perfectly fit states, whereas Eq. (24) here uses the fitness with respect to the normal-state Hamiltonian and predicts only parametric enhancement. Because both calculations start from the same SCBA equations (8)-(10), this is a genuine unresolved discrepancy about the central result, not a cosmetic difference. The authors must identify the source of the discrepancy and either prove that the normal-state fitness is the correct quantity or revise Eq. (24) accordingly. This issue is load-bearing: Eq. (24) is the paper's main quantitative diagnostic and underpins the material-specific claims in Sections III and IV. The brief acknowledgment in the final paragraph is not sufficient to establish the central claim.","section":"Section V, final paragraph; Eq. (24)"},{"comment":"The reduction of the anomalous self-energy to Eq. (10) neglects interband contributions on the strength of the 'well separated bands' assumption. This assumption is not quantified, and it is particularly delicate for YPtBi, which is described as a zero-band-gap semimetal with the chemical potential in the lower band. The effective-rate formulas (16) and (21), and the universal form of Eq. (24), depend on this approximation. Please provide an estimate of the omitted interband terms for the two model Hamiltonians, or state explicitly the conditions under which they vanish by symmetry. Without this, the numerical T_c curves in Figs. 3 and 4 should be regarded as illustrative rather than quantitative.","section":"Section II, text following Eq. (10)"},{"comment":"The sentence 'This result follows from the observation that λ_l = +1 (−1) when γ_l ~∆_ν − ~∆_ν γ_{l,*} = 0 (2γ_l ~∆_ν)' is the entire derivation of the central formula. The connection between this commutator condition and the Fermi-surface average of the normalized fitness ~F_C is not shown. Given that Eq. (24) is the headline result, the paper should present the explicit algebra, or an appendix, that converts Eqs. (16) and (21) into Eq. (24), including the exact meaning of γ_{l,*} and the role of H^T_{-k} in the definition of F_C. This is needed both for verifiability and for resolving the relation to Ref. [26].","section":"Section V, derivation of Eq. (24)"}],"minor_comments":[{"comment":"The phrase 'it was it was shown in Ref. [20]' contains a duplicated clause, and 'psuedospin' should be 'pseudospin'.","section":"Introduction, second paragraph"},{"comment":"The captions state that the line τ_ν = τ applies to pairing states in all other nontrivial irreps, but the curves are not individually identified. Please add a legend or explicitly name the pairing state for each curve so the reader can connect the plots to Tables I and II.","section":"Figures 3 and 4 captions"},{"comment":"The quantity α_ν is introduced without fully specifying the notation: please state that ~∆_ν is the s-wave basis state of channel ν and that the Fermi-surface average is taken over the band(s) at the Fermi energy.","section":"Eq. (18)"},{"comment":"Reference [26] is cited as a preprint without journal, volume, or page details; the published version should include the full citation and a more detailed comparison in the text.","section":"Reference list"}],"recommendation":"major_revision","confidential_remarks":"The unresolved disagreement with Ref. [26] is the central obstacle. If the authors cannot resolve it in revision, the paper should not be published in its current form. The editor may wish to have the manuscript reviewed in tandem with Ref. [26] to determine which calculation is correct."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know: this is a careful, self-consistent SCBA treatment of impurity effects on Tc in four-component multiband superconductors, and it produces a clean formula linking the effective scattering rate to the Fermi-surface-averaged superconducting fitness. The catch is that the formula is currently in dispute with a concurrent paper (Ref 26), and the authors openly say so. Until that discrepancy is resolved, Eq. (24) is not established.\n\nWhat the paper does well: the Euclidean Dirac matrix formalism is a clean generalization that keeps track of orbital and spin degrees of freedom explicitly. The derivation of Eqs. (16) and (21) is internally consistent, the assumptions (isotropic impurity potential, well-separated bands, SCBA) are stated, and the applications to YPtBi and CuxBi2Se3 are concrete. The A1u robustness result for CuxBi2Se3 agrees with earlier work, a good sanity check. The YPtBi analysis appears new, as does the explicit connection to the superconducting fitness.\n\nThe soft spots, in order. First, the unresolved disagreement with Ref 26 is real and not cosmetic. Both papers start from the same SCBA equations, but the present paper uses the normal-state Hamiltonian in the fitness, while the concurrent paper uses the impurity Hamiltonian. One predicts parametric enhancement; the other predicts complete insensitivity for perfectly-fit states. The authors flag this but do not identify the source of the discrepancy. That makes the central quantitative claim and the material-selection diagnostic load-bearing and unverified. Second, the isotropic impurity potential assumption is a real limitation, and the authors acknowledge it. If impurities couple anisotropically to orbital or spin degrees of freedom, the effective rates could change substantially. Third, the well-separated-band approximation is fine for the specific materials considered but could break down in more complex systems. These are secondary, not fatal.\n\nWho this is for: people working on disorder in multiband or topological superconductors, especially YPtBi and CuxBi2Se3. A serious referee should engage with this paper, mainly because the dispute with Ref 26 needs adjudication. I would accept it for peer review, with the expectation of either a revision that resolves the discrepancy or a clear statement that the two results are different and why. It deserves reviewer time, but I would not treat Eq. (24) as final until the disagreement is settled.","headline":"A clean SCBA framework for multiband disorder, but the headline fitness formula is disputed by a concurrent paper and needs refereeing.","tokens_in":11766,"tokens_out":5287,"would_cite":true,"duration_ms":48408,"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":"In multiband superconductors with four internal electron degrees of freedom, unconventional s-wave pairing states have a reduced effective disorder-scattering rate, making them more resilient than single-band theory predicts.","keywords":["unconventional superconductivity","disorder","Anderson's theorem","superconducting fitness","multiband superconductivity","topological superconductor","self-consistent Born approximation","s-wave pairing"],"falsifier":"Measure the suppression of the transition temperature in a single-band-Fermi-surface sample of Cu$_x$Bi$_2$Se$_3$ as nonmagnetic disorder is introduced by electron irradiation; if the effective scattering rate extracted from the suppression equals the normal-state rate rather than the reduced value in Eq. (21), the central claim fails. A complementary test is to use impurities with known orbital selectivity and check whether the protection weakens, which would expose the isotropic-potential assumption.","tokens_in":10760,"feed_emoji":"🛡️","tokens_out":16281,"duration_ms":130661,"temperature":0.7,"pith_summary":"This paper argues that in multiband superconductors where electrons carry four internal degrees of freedom, pairing states that are formally s-wave but transform nontrivially under crystal symmetries suffer far less damage from nonmagnetic impurities than standard single-band theory would predict. Using the self-consistent Born approximation, the authors derive a general expression for the critical temperature under disorder and show that the effective scattering rate for such an s-wave channel is reduced from the normal-state value by an amount set by the Fermi-surface average of the superconducting fitness. Applied to the candidate topological superconductors YPtBi and Cu$_x$Bi$_2$Se$_3$, the theory predicts enhanced resilience for the novel s-wave states, and partial protection for any other pairing state in the same irreducible representation according to its similarity to the s-wave state at the Fermi surface. The result matters because it overturns the usual expectation that unconventional pairing is fragile against disorder and provides a computable criterion for identifying disorder-resistant superconducting channels in real materials.","feed_headline":"Unconventional s-wave states resist impurities better than expected","feed_subtitle":"Multiband metals can shield nontrivial s-wave pairing from impurities, guiding topological superconductor search.","key_machinery":"The central machinery is the anomalous impurity self-energy $\\Sigma_2$ evaluated in the self-consistent Born approximation, whose key feature is that the Fermi-surface average of the band-projected pairing potential $P_{k,j}\\tilde\\Delta_k P^T_{-k,j}$ does not vanish because the projection operators $P_{k,\\pm}=(1\\pm\\hat{\\epsilon}_k\\cdot\\vec{\\gamma})/2$ are momentum-dependent matrices. This yields the effective scattering rate formula Eq. (24) in terms of the Fermi-surface-averaged superconducting fitness, and the similarity parameter $\\alpha_\\nu$ of Eq. (18) that controls how much protection other states in the same irrep inherit.","core_discovery":"In a time-reversal- and inversion-symmetric two-band system with four internal degrees of freedom, the anomalous self-energy generated by isotropic nonmagnetic impurities does not vanish for unconventional pairing states, because the band projection operators carry the nontrivial spin-orbital texture of the normal-state Hamiltonian. For an s-wave pairing channel $\\nu$, the effective scattering rate that enters the $T_c$ suppression formula becomes $1/\\tau_\\nu = 1/\\tau - (1/\\tau_0)(1 - \\bar F_C)$, where $\\bar F_C$ is the Fermi-surface average of the normalized superconducting fitness $\\tilde F_C(k)$ that measures the fraction of interband pairing. Thus a state with zero fitness is completely insensitive to nonmagnetic disorder, generalizing Anderson's theorem, while realistic states with small fitness retain a parametrically enhanced robustness. The same anomalous self-energy is nonzero and momentum-independent for every pairing state in the same irreducible representation, so the s-wave channel sets an upper bound on the disorder stability of all states in that irrep, quantified by the overlap $\\alpha_\\nu$ between the state and the s-wave gap at the Fermi surface.","pith_inferences":["A practical next step would be to compute $\\bar F_C$ from first-principles band structures for other candidate multiband superconductors, turning the fitness criterion into a screening tool that ranks materials by expected disorder tolerance before any irradiation experiment.","The predicted difference in $T_c$ suppression between trivial and nontrivial s-wave channels is sharp enough that controlled electron-irradiation experiments on Cu$_x$Bi$_2$Se$_3$ or YPtBi could discriminate between candidate pairing symmetries.","The theory assumes a scalar impurity potential in the orbital basis; real defects often couple selectively to orbitals or spins, so an extension to anisotropic impurity scattering would show whether the protection survives in actual materials or is a property of the simplest model."],"forward_implications":["Nontrivial s-wave channels in any four-degrees-of-freedom multiband superconductor are generically more resilient to nonmagnetic disorder than sign-changing single-band states, with the resilience set by the Fermi-surface-averaged superconducting fitness.","Other pairing states in the same irreducible representation share this protection; for example, the d-wave $E_g$ state in YPtBi is nearly as stable against disorder as the quintet s-wave $E_g$ states because of its high overlap with them at the Fermi surface.","Systems with a nontrivial inversion operator are especially favorable: odd-parity s-wave states commute with three of the five gamma matrices in the generic Hamiltonian, typically giving smaller fitness and hence greater robustness, as seen in the comparison of Cu$_x$Bi$_2$Se$_3$ with YPtBi.","The theory offers a practical diagnostic for materials search: evaluate the Fermi-surface-averaged fitness of the s-wave channel in each irreducible representation; a value $\\bar F_C \\ll 1$ identifies a candidate disorder-resistant unconventional superconductor."],"supporting_citations":[{"why":"Establishes Anderson's theorem for conventional s-wave states, the standard result from which the paper's nontrivial s-wave channels depart.","marker":"[2]"},{"why":"Supplies the single-band disorder theory and the digamma-function $T_c$ suppression formula that the paper generalizes.","marker":"[1]"},{"why":"Introduces the odd-parity s-wave states in Cu$_x$Bi$_2$Se$_3$ and the global transformation showing perfectly fit states are disorder-insensitive.","marker":"[9]"},{"why":"Provides the YPtBi $j=3/2$ pairing states and model parameters used in the first material application.","marker":"[10]"},{"why":"Shows that spin-orbital locking protects the A1u state in Cu$_x$Bi$_2$Se$_3$, the specific case that this paper generalizes.","marker":"[20]"},{"why":"Confirms numerically with a self-consistent T-matrix theory that the mass term controls disorder robustness in Cu$_x$Bi$_2$Se$_3$.","marker":"[21]"},{"why":"Defines the superconducting fitness $F_C$ whose Fermi-surface average appears in the effective scattering rate, Eq. (24).","marker":"[24]"},{"why":"Gives the generic form of the inversion-symmetric four-band normal-state Hamiltonian in Eq. (1).","marker":"[25]"},{"why":"Supplies the $k\\cdot p$ Hamiltonian for Cu$_x$Bi$_2$Se$_3$ used in the second material application.","marker":"[29]"}],"fun_headline_variants":["Unconventional s-wave pairing defies disorder via spin-orbit texture","Multiband metals shield exotic s-wave states from impurities","Superconducting fitness quantifies disorder robustness in multiband systems","Nontrivial s-wave states survive impurities better via interband pairing","Spin-orbit texture boosts disorder tolerance of unconventional s-wave states"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes impurities are identical point scatterers that affect all four internal electron states equally, and that the two bands are well separated in energy; if real disorder couples preferentially to particular orbital or spin states, or the bands are close enough to hybridize, the predicted reduction in the effective scattering rate could be overstated.","fun_headline_variants_meta":{"raw":{"variants":["Unconventional s-wave pairing defies disorder via spin-orbit texture","Multiband metals shield exotic s-wave states from impurities","Superconducting fitness quantifies disorder robustness in multiband systems","Nontrivial s-wave states survive impurities better via interband pairing","Spin-orbit texture boosts disorder tolerance of unconventional s-wave states"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000538,"raw_usage":{"total_tokens":2571,"prompt_tokens":923,"completion_tokens":1648,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":539,"completion_tokens_details":{"reasoning_tokens":1558}},"tokens_in":539,"tokens_out":1648,"duration_ms":12623,"temperature":1.0,"reasoning_tokens":1558,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:10:36.326839+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the suppression of the transition temperature in a single-band-Fermi-surface sample of Cu$_x$Bi$_2$Se$_3$ as nonmagnetic disorder is introduced by electron irradiation; if the effective scattering rate extracted from the suppression equals the normal-state rate rather than the reduced value in Eq. (21), the central claim fails. A complementary test is to use impurities with known orbital selectivity and check whether the protection weakens, which would expose the isotropic-potential assumption.","supporting_citations":[{"cited_title":"Model Hamiltonian for topological insula- tors,","cited_arxiv_id":null,"evidence_quote":"Supplies the $k\\cdot p$ Hamiltonian for Cu$_x$Bi$_2$Se$_3$ used in the second material application."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the single-band disorder theory and the digamma-function $T_c$ suppression formula that the paper generalizes."},{"cited_title":"Odd-Parity Topological Supercon- ductors: Theory and Application to Cu xBi2Se3,","cited_arxiv_id":null,"evidence_quote":"Introduces the odd-parity s-wave states in Cu$_x$Bi$_2$Se$_3$ and the global transformation showing perfectly fit states are disorder-insensitive."},{"cited_title":"Pairing of j = 3/2 Fermions in Half-Heusler Su- perconductors,","cited_arxiv_id":null,"evidence_quote":"Provides the YPtBi $j=3/2$ pairing states and model parameters used in the first material application."},{"cited_title":"Spin-Orbit Locking as a Pro- tection Mechanism of the Odd-Parity Superconducting State against Disorder,","cited_arxiv_id":null,"evidence_quote":"Shows that spin-orbital locking protects the A1u state in Cu$_x$Bi$_2$Se$_3$, the specific case that this paper generalizes."},{"cited_title":"Robust superconductivity with nodes in the superconducting topological insulator Cu xBi2Se3: Zee- man orbital field and nonmagnetic impurities,","cited_arxiv_id":null,"evidence_quote":"Confirms numerically with a self-consistent T-matrix theory that the mass term controls disorder robustness in Cu$_x$Bi$_2$Se$_3$."},{"cited_title":"Tailoring Tc by symmetry principles: The concept of supercon- ducting fitness,","cited_arxiv_id":null,"evidence_quote":"Defines the superconducting fitness $F_C$ whose Fermi-surface average appears in the effective scattering rate, Eq. (24)."},{"cited_title":"Bogoliubov Fermi surfaces: General theory, magnetic order, and topology,","cited_arxiv_id":null,"evidence_quote":"Gives the generic form of the inversion-symmetric four-band normal-state Hamiltonian in Eq. (1)."}],"review_version":1}