{"id":"bb1c7a0e-b561-4f9d-9ec7-e016cd1eed0c","arxiv_id":"1908.08766","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A generalized Anderson theorem protects nodal superconductivity in multiorbital materials when impurity scattering does not mix the internal degrees of freedom.","lead":"This paper proposes a generalized Anderson theorem: nodal superconductors with extra internal degrees of freedom are protected from nonmagnetic impurities as long as the pairing is isotropic and the impurities do not mix those degrees of freedom. It also reports thermal conductivity data on the doped topological insulator CPSBS showing that a nodal superconductor survives despite a scattering rate far larger than its gap.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Protection requires V proportional to τ0, not merely orbital-diagonal scattering: an inversion-even τ3 term allowed by symmetry is omitted from Eq. (8) and gives nonzero fitness for the E_u gap.","rationale":"The generalized Anderson theorem is a clean algebraic result, and the experimental evidence for nodal quasiparticles in CPSBS is persuasive: sample II has essentially full superconducting volume and still shows residual heat conduction incompatible with a fully gapped state. The reader's conditional verdict is appropriate. However, the reader's weakest assumption — that impurities must break inversion symmetry for the protection to fail — is not the sharpest point. Even granting full inversion symmetry of the impurity potential, the allowed τ3⊗σ0 term is diagonal in orbital space and does not mix internal degrees of freedom, yet it does not commute with inter-orbital odd-parity order parameters such as the E_u state of Eq. (7). Consequently, the statement that all nonmagnetic scattering not mixing internal degrees of freedom is harmless is internally too broad; the actual condition is that V commute with Δ, which for inter-orbital pairing requires V ∝ τ0. This does not invalidate the theorem for the idealized case V = V0 τ0⊗σ0, nor does it overturn the nodal experimental data, but it makes the application to CPSBS dependent on the unstated assumption V11 = V22. That assumption is testable by symmetry analysis or DFT-based impurity modeling. The correct response is therefore to keep the conditional verdict, with the condition sharpened to include verification that the τ3 component of the impurity potential is negligible.","tokens_in":21635,"tokens_out":15138,"duration_ms":180375,"concrete_test":"Recompute the effective scattering rate in Eq. (5) with V = V0 τ0⊗σ0 + V3 τ3⊗σ0 and the E_u order parameter of Eq. (7), evaluating FC from Eq. (6); the result is nonzero and proportional to V3^2, demonstrating that Eq. (8)'s reduction to V0 alone is not symmetry-enforced. To settle whether this matters for CPSBS, compute V3/V0 for a realistic Cu interstitial using a DFT Wannier projection of the impurity potential onto the P1z+ and P2z− orbitals; if V3 is comparable to V0, the claimed protection of the nodal E_u state is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that superconductors with momentum-dependent band-basis gaps are generically protected from nonmagnetic scattering that does not mix internal degrees of freedom. The argument reduces the impurity potential to V = V0 τ0⊗σ0 + Vs S·(τ0⊗σ) using parity: the overlap between opposite-parity orbitals vanishes for a symmetric potential, so V has no τ1 or τ2 terms. But this conclusion is too strong. Inversion symmetry forbids only off-diagonal orbital terms; it allows τ3⊗σ0, which is diagonal in orbital space and does not mix the internal degrees of freedom. Table S1 indeed lists (3,0) as an allowed A1g term in the normal-state Hamiltonian, and nothing forces the diagonal impurity matrix elements on P1z+ and P2z− to be equal. Writing V = V0 τ0⊗σ0 + V3 τ3⊗σ0, the fitness function for the E_u order parameter Δ = Δ0[iτ2⊗σ1(iσ2)] of Eq. (7) is nonzero because [τ3, τ2] ≠ 0. Substituting this V into Eq. (5) gives ΓEff ∝ V3^2 > 0, so the predicted zero pair-breaking rate is lost. The same applies to the other inter-orbital odd-parity channels (A1u, A2u). Thus the theorem as stated is not valid for all nonmagnetic scattering that fails to mix internal degrees of freedom; it holds only when V is strictly proportional to the identity in the internal-orbital space. The assumption V11 = V22 is not derived from inversion symmetry and is not stated. Since the CPSBS explanation requires exactly this unargued equality, the application of the theorem to real impurities is less secure than the text suggests, even before considering inversion-broken impurity sites.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a generalized Anderson theorem for superconductors with multiple internal degrees of freedom. The central idea is that when the pairing interaction is isotropic in the local orbital basis, a nonmagnetic impurity potential that is proportional to the identity in the internal space commutes with every allowed gap matrix, so the superconducting fitness function vanishes and the effective pair-breaking rate ΓEff is zero, even if the gap is nodal in the band basis. The authors apply this framework to the Bi2Se3-based superconductor CPSBS, combining specific-heat and ultra-low-temperature thermal-conductivity data to argue that nodal superconductivity persists although the normal-state scattering rate is more than an order of magnitude larger than the gap. The paper also derives ΓEff from the Abrikosov-Gorkov equations and expresses it in terms of the fitness function. The algebraic core—V ∝ τ0⊗σ0 commutes with all gap matrices—is sound, but the broader claim that any nonmagnetic scattering that does not mix internal degrees of freedom is harmless is not justified by the symmetry argument presented.","tokens_in":22010,"tokens_out":12402,"duration_ms":131402,"significance":"If the theorem as stated were correct, it would provide a general mechanism protecting unconventional multiorbital superconductors against disorder and would resolve the long-standing puzzle of robustness in Bi2Se3-based materials. The experimental data on CPSBS are of high quality and the use of the superconducting fitness function is elegant. The machine-checkable algebraic statement that commutativity implies zero pair-breaking is correct and is a useful contribution. However, the central claim is overstated: inversion symmetry does not imply V ∝ identity, and the allowed diagonal orbital term τ3⊗σ0 breaks the protection for the odd-parity inter-orbital order parameters. This severely limits the generality of the proposed theorem and undermines the claim that the CPSBS robustness is explained by the generalized Anderson mechanism as presented.","major_comments":[{"comment":"The reduction of the impurity potential to V = V0 τ0⊗σ0 + Vs S·(τ0⊗σ) is not derived from the stated assumption of a symmetric impurity potential. Inversion symmetry forbids only the off-diagonal orbital terms τ1⊗σ0 and τ2⊗σ0; it allows a diagonal orbital term V3 τ3⊗σ0. Since the two effective orbitals P1z+ and P2z− are inequivalent—indeed Table S1 lists the (3,0) A1g term in the normal-state Hamiltonian—a local impurity potential will generically produce different on-site matrix elements on the two orbitals, so V3 ≠ 0. For the Eu order parameter of Eq. (7), the commutator [τ3, τ2] is nonzero, giving FC ≠ 0 and ΓEff ∝ V3^2; the same applies to the A1u and A2u odd-parity channels. Therefore the theorem as stated does not hold for all nonmagnetic scattering that fails to mix internal degrees of freedom; it holds only when V is strictly proportional to the identity in orbital space. The assumption V11 = V22 is not stated anywhere and is not a consequence of inversion symmetry. Since the CPSBS explanation relies on exactly this equality, this issue is load-bearing and must be addressed, either by proving V3 = 0 from a microscopic model or by reformulating the theorem with V ∝ identity as an explicit assumption and discussing the physical consequences for the experimental interpretation.","section":"Robust superconductivity in the Bi2Se3-based materials, Eq. (8)"},{"comment":"The statement that 'the conclusion of zero scattering rate is valid for any SC order parameter possible for the Bi2Se3-based materials, because the identity matrix τ0⊗σ0 commutes with any Δ' is only valid if the impurity potential has already been restricted to be proportional to τ0. The paper's broader claims—for example, that superconductors with momentum-dependent band-basis gaps are 'generically protected from nonmagnetic scattering that do not mix the internal DOF'—are false under the symmetries of the problem, because τ3⊗σ0 is diagonal in orbital space, does not mix internal degrees of freedom, and yet gives a nonzero fitness function for the odd-parity inter-orbital order parameters. The authors should either provide a microscopic justification for the vanishing of the τ3 component, or explicitly restrict the theorem to impurity potentials proportional to the identity in the internal space and clearly state that this is a necessary condition, not a consequence of symmetry.","section":"Abstract and main text"},{"comment":"The derivation of the effective scattering rate and the conclusion that ΓEff = 0 when FC = 0 are obtained within the self-consistent Born approximation. The experimental regime is ℏΓ/Δ0 ≳ 60, which lies far outside the weak-scattering limit. The paper states that the protection holds 'even for arbitrarily large τ1 and τ2 in the normal state,' but this statement is only established within the Born approximation. To support the strong-disorder conclusion, the authors should either extend the argument beyond the Born approximation (e.g., via a T-matrix treatment) or qualify the claim to the regime of validity of the derivation.","section":"Supplementary Materials, Eqs. (27)–(32)"}],"minor_comments":[{"comment":"There is a typo: 'CPSPB' should be 'CPSBS' in the sentence discussing the Eu channel.","section":"Main text, after Eq. (8)"},{"comment":"The matrix elements in the unitary transformation are typeset ambiguously; for instance 'h11hv−ih12h30/hp' should read '(h11 hv − i h12 h30)/hp'. Please correct the parentheses throughout the matrix.","section":"Supplementary Materials, Eq. (S9)"},{"comment":"The notation 'the horizontal bar indicates impurity averaging' is unclear; please specify explicitly that the overline denotes the configurational average over impurity positions.","section":"Main text, Eq. (5)"},{"comment":"The row 'Dimension' contains the entry '2D?' with a question mark for one material; please resolve this ambiguity.","section":"Supplementary Materials, Table S4"},{"comment":"The abstract states that the scattering-rate energy scale is 'orders of magnitude larger' than the gap, while the main text says 'more than an order of magnitude'; please harmonize the wording.","section":"Abstract and main text"},{"comment":"In the phrase 'nonmagnetic scattering that do not mix the internal DOF,' the verb should agree with the singular subject 'scattering'—use 'does not mix.'","section":"Main text, Abstract"}],"recommendation":"major_revision","confidential_remarks":"The manuscript has a strong experimental component and a clean algebraic observation, but the theoretical claim as stated is not supported. The τ3⊗σ0 issue is a genuine correctness problem for the generalized theorem and for the interpretation of the CPSBS data. I recommend major revision rather than rejection because the underlying algebra (commutativity implies zero pair-breaking) is sound and the authors may be able to reformulate the theorem with the necessary assumption of V ∝ identity, or provide a microscopic argument for why V3 is negligible. The present version, however, should not be accepted without addressing this point."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThis paper is worth reading, but the headline theorem is overstated. What is genuinely new: it takes the superconducting fitness function and shows that the pair-breaking rate in an Abrikosov-Gor'kov treatment can be written as a Fermi-surface average of |FC|^2. For isotropic pairing in the orbital basis, any impurity potential proportional to identity in the internal space gives zero pair breaking, so nodal states can survive strong disorder. The experimental part on CPSBS is also solid: the thermal conductivity in a sample that is essentially 100% superconducting shows residual mobile quasiparticles at 50 mK, which is hard to explain without nodes, and the Drude estimate puts hbar*Gamma well above the gap. The algebraic connection between FC and the effective scattering rate looks right.\n\nThe soft spot is the impurity-potential assumption. The text claims that inversion symmetry plus opposite-parity orbitals forces V = V0 tau0 sigma0, but inversion symmetry only forbids inter-orbital terms tau1 and tau2. It allows tau3 sigma0 (diagonal but unequal), and nothing in the paper rules out different diagonal matrix elements on the two orbitals. The stress-test is correct: for the E_u order parameter, [tau3, tau2] is nonzero, so a V3 tau3 term gives nonzero fitness and pair breaking. The statement that tau0 sigma0 commutes with any Delta is true, but only if V is actually proportional to tau0 sigma0, and the paper has not shown that. For a delta-function impurity, the diagonal matrix elements on the two parity orbitals are generally not equal, so the protection is not generic. The theorem should be stated as: protected when V is proportional to tau0 sigma0, not for all nonmagnetic scattering that fails to mix orbitals.\n\nThis is load-bearing for the claimed explanation of CPSBS. It is addressable: one could compute V1 and V2 for realistic impurity models, or argue that the intercalant disorder is screened and long-wavelength so the potential is nearly constant across orbitals. Until then, the application to real impurities is less secure than the text suggests. Minor issue: the scattering-rate estimate lacks uncertainty propagation, and the evidence for nodes, while very suggestive, is not a formal proof.\n\nWho is this for? Anyone working on doped topological insulators or multiorbital superconductors. I would send it to a serious referee rather than desk reject, but the referee should push on the impurity assumption and the generality claim. I would bring it to reading group and would likely cite the fitness formalism after the assumption is cleaned up.","headline":"The fitness-function machinery is genuinely useful and the CPSBS thermal conductivity is impressive, but the generalized Anderson theorem as stated rests on an unjustified V proportional to identity assumption and is not generic for nonmixing impurities.","tokens_in":22503,"tokens_out":2694,"would_cite":true,"duration_ms":31689,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.20.Fg","74.25.F-","74.70.-b"],"model":"deepseek-v4-flash","headline":"This paper argues that nodal superconductivity in multiorbital materials is protected from nonmagnetic disorder when pairing is isotropic in the orbital basis and impurities do not mix internal degrees of freedom, with…","keywords":["superconducting fitness","Anderson's theorem","nodal superconductivity","disorder robustness","topological superconductors","multiorbital superconductivity","Bi2Se3-based superconductors","CPSBS"],"falsifier":"Introduce into CPSBS a controlled concentration of defects that are known to sit off inversion centers, and measure whether $T_c$ is suppressed along the generalized Abrikosov-Gor'kov curve as the defect density rises; observation of such a suppression would show that inversion-breaking scattering activates the fitness function. A null result—unchanged $T_c$ despite clear inversion-breaking scattering—would contradict the paper's central claim.","tokens_in":21471,"feed_emoji":"🛡️","tokens_out":9894,"duration_ms":89134,"temperature":0.7,"pith_summary":"Conventional wisdom holds that nonmagnetic impurities destroy nodal superconductors, yet Bi2Se3-based superconductors remain robust under strong disorder. This paper explains the puzzle by generalizing Anderson's theorem to superconductors with extra internal degrees of freedom such as orbitals: when the pairing interaction is isotropic in the microscopic orbital basis, the momentum-dependent gap that appears in the band basis is protected from any nonmagnetic scattering that does not mix the internal degrees of freedom. The pair-breaking rate is controlled by a superconducting fitness function, and for the Bi2Se3 family the relevant scattering potential is proportional to the identity matrix in orbital and spin space, which commutes with every allowed gap matrix and gives zero scattering rate. The paper reports that in CPSBS, thermal conductivity down to 50 mK reveals residual mobile quasiparticles—nodes—while the estimated scattering rate is more than an order of magnitude larger than the superconducting gap, exactly the regime in which ordinary nodal superconductivity would be suppressed.","feed_headline":"When disorder cannot mix orbitals, nodal superconductivity survives","feed_subtitle":"CPSBS heat transport shows intact nodal order even when scattering is 60 times the gap.","key_machinery":"The central object is the superconducting fitness function $F_C(\\mathbf{k}-\\mathbf{k}') = V(\\mathbf{k}-\\mathbf{k}')\\Delta - \\Delta V^*(\\mathbf{k}-\\mathbf{k}')$, a modified commutator between the impurity scattering potential and the gap matrix; the effective pair-breaking rate is its normalized Fermi-surface trace average. The argument uses the fact that in Bi2Se3-based superconductors the nonmagnetic scattering potential has the form $V = V_0\\tau_0\\otimes\\sigma_0$, because the two effective orbitals have opposite parity and the overlap integral of a symmetric impurity potential between them vanishes. Since the identity in orbital and spin space commutes with every allowed gap matrix of the form $\\tau_a\\otimes\\sigma_b$, the fitness function is identically zero, so the impurity-induced pair-breaking rate vanishes even though the gap is nodal in the band basis. A unitary transformation from the orbital basis to the band basis is the step that converts the constant orbital-basis gap into a momentum-dependent nodal structure, and the node positions are set purely by normal-state parameters such as $h_{20}(\\mathbf{k})\\sim B_0 k_z$ and $h_{12}(\\mathbf{k})\\sim A_0 k_x$.","core_discovery":"The central claim is a generalized Anderson's theorem: in a superconductor with multiple internal degrees of freedom, an isotropic pairing interaction written in the local orbital basis produces a momentum-dependent gap once projected onto the band basis, and that gap is immune to nonmagnetic impurities as long as the impurity potential does not mix the internal degrees of freedom. The argument is carried by the superconducting fitness function $F_C = V\\Delta - \\Delta V^*$, whose Fermi-surface average determines the effective pair-breaking rate $\\hbar\\Gamma_{\\rm eff} = \\frac{1}{4}\\langle\\mathrm{Tr}[\\tilde F_C^\\dagger\\tilde F_C]\\rangle_k$. For the Bi2Se3-based materials, the two effective orbitals have opposite parity, so a symmetric impurity potential reduces to $V = V_0\\tau_0\\otimes\\sigma_0$, the identity in both orbital and spin space; because this commutes with every allowed gap matrix of the form $\\tau_a\\otimes\\sigma_b$, the fitness function vanishes and the scattering rate is zero. The paper presents CPSBS as the extreme demonstration: thermal conductivity data show unambiguous nodal quasiparticle heat transport, while the Drude estimate gives $\\hbar\\Gamma \\gtrsim 30$ meV compared with $\\Delta_0 \\simeq 0.5$ meV, so strong scattering coexists with a nodal superconducting state instead of destroying it.","pith_inferences":["A direct extension would be to test the boundary of the protection by deliberately introducing inversion-breaking defects, which should add orbital-mixing terms to the impurity potential and restore pair-breaking; the predicted $T_c$ suppression as a function of such defect density is quantitatively given by the fitness formula.","The same fitness logic suggests a screening criterion for other multiorbital or multi-valley superconductors: if the dominant impurity potential is block-diagonal in the internal degrees of freedom and the pairing is isotropic in the local basis, nodal gaps should be disorder-robust.","Because bulk nodal quasiparticles dominate the low-temperature heat transport in nodal topological superconductors, the surface Majorana contribution remains tiny in bulk samples; isolating a quantized Majorana thermal signal would require ultra-clean thin geometries or surface-sensitive measurements."],"forward_implications":["Any superconductor whose nonmagnetic impurity potential is proportional to the identity in all internal spaces has zero effective pair-breaking rate, so its $T_c$ is not suppressed by disorder even if the band-basis gap has nodes.","The Abrikosov-Gor'kov-type formula $\\log(T_c/T_c^0)=\\Psi(1/2)-\\Psi(1/2+\\hbar\\Gamma_{\\rm eff}/2\\pi k_B T_c)$ remains the correct description of disorder-suppressed $T_c$ in multiorbital systems provided $\\hbar\\Gamma_{\\rm eff}$ is computed from the fitness function.","The observed robustness of superconductivity in Bi2Se3-based materials against disorder is not evidence against unconventional pairing; it is the expected behaviour for isotropic pairing in the orbital basis with parity-preserving scattering.","In CPSBS, the absence of universal thermal conductivity—with $\\kappa_0/T$ about three orders of magnitude below the clean-limit estimate—shows that strong scattering persists in the superconducting state, and the survival of the nodal signature is a direct consequence of the generalized theorem."],"supporting_citations":[{"why":"Anderson's original theorem for dirty superconductors, which this work generalizes to multiple internal degrees of freedom.","marker":"[20]"},{"why":"Supplies the two-orbital model with isotropic pairing for Bi2Se3-based superconductors and the order-parameter symmetry analysis used throughout.","marker":"[22]"},{"why":"Introduced the superconducting fitness function as a measure of incompatibility between the normal-state Hamiltonian and the gap matrix in multiorbital systems.","marker":"[25]"},{"why":"Formalized the concept of superconducting fitness that the present paper extends to impurity scattering.","marker":"[26]"},{"why":"Previous study of CPSBS establishing the nematic superconducting state with nodes, providing the sample-specific context and specific-heat data used here.","marker":"[19]"},{"why":"Reports CPSBS as a superconductor derived from a topological insulator heterostructure and supplies the normal-state resistivity used in the Drude estimate.","marker":"[18]"},{"why":"Line-nodal gap theory used to fit the electronic specific heat and extract the superconducting volume fraction of the samples.","marker":"[50]"},{"why":"Documented robust odd-parity superconductivity against disorder in NbxBi2Se3, a motivating experimental observation for the generalized theorem.","marker":"[16]"}],"fun_headline_variants":["Generalized Anderson theorem shields nodal superconductors from disorder","Nodal superconductivity survives disorder via generalized Anderson theorem","Impurity immunity in Bi2Se3 superconductors explained by new theorem","Beyond Anderson: nodal order persists under strong scattering","CPSBS shows nodal order immune to scattering 60 times gap"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the impurity potential is inversion-symmetric, so scattering never mixes the two opposite-parity orbitals; if real defects sit off inversion centers or have an odd-parity component, the predicted zero scattering rate fails.","fun_headline_variants_meta":{"raw":{"variants":["Generalized Anderson theorem shields nodal superconductors from disorder","Nodal superconductivity survives disorder via generalized Anderson theorem","Impurity immunity in Bi2Se3 superconductors explained by new theorem","Beyond Anderson: nodal order persists under strong scattering","CPSBS shows nodal order immune to scattering 60 times gap"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000206,"raw_usage":{"total_tokens":1420,"prompt_tokens":989,"completion_tokens":431,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":605,"completion_tokens_details":{"reasoning_tokens":349}},"tokens_in":605,"tokens_out":431,"duration_ms":4926,"temperature":1.0,"reasoning_tokens":349,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:30:28.883742+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Introduce into CPSBS a controlled concentration of defects that are known to sit off inversion centers, and measure whether $T_c$ is suppressed along the generalized Abrikosov-Gor'kov curve as the defect density rises; observation of such a suppression would show that inversion-breaking scattering activates the fitness function. A null result—unchanged $T_c$ despite clear inversion-breaking scattering—would contradict the paper's central claim.","supporting_citations":[{"cited_title":"Impurity scattering in superconductors,","cited_arxiv_id":null,"evidence_quote":"Supplies the two-orbital model with isotropic pairing for Bi2Se3-based superconductors and the order-parameter symmetry analysis used throughout."},{"cited_title":"Identifying detrimental effects for multiorbital supercon- ductivity: Application to sr 2ruo4,","cited_arxiv_id":null,"evidence_quote":"Formalized the concept of superconducting fitness that the present paper extends to impurity scattering."},{"cited_title":"Topological superconductors: a review,","cited_arxiv_id":null,"evidence_quote":"Reports CPSBS as a superconductor derived from a topological insulator heterostructure and supplies the normal-state resistivity used in the Drude estimate."},{"cited_title":"Anomalous suppression of the superfluid density in the CuxBi2Se3 superconductor upon progressive cu intercalation,","cited_arxiv_id":null,"evidence_quote":"Documented robust odd-parity superconductivity against disorder in NbxBi2Se3, a motivating experimental observation for the generalized theorem."}],"review_version":1}