{"id":"2c682102-0f16-4ea9-8b49-cbc39be392aa","arxiv_id":"2501.01155","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Impurity-induced in-gap states in superconductors are tied to sign-changing order parameters, and the effect can be used to test pairing symmetry in La3Ni2O7.","lead":"This paper uses T-matrix theory to show that impurity-induced bound states appear in unconventional superconductors exactly when the superconducting gap changes sign along the Fermi surface. It applies the idea to the nickelate superconductor La3Ni2O7, suggesting that STM scans of impurities could reveal whether its pairing is interlayer sign-changing or sign-preserving.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (8) is not Det[I − U G0] as defined: for U = V τ_z the exact zero-energy denominator is (1 + VΣP)² + (VΣQ)², so the dropped 1 + 2VΣP terms change the pole condition and the analytic link between sign reversal and in-gap states is not established.","rationale":"I read the paper in good faith: its central claim rests on two pillars, a numerical T-matrix LDOS calculation (Eqs. 5–7) and an analytic derivation (Eqs. 8 and 12). The numerics are standard and the qualitative physics — mid-gap impurity states for sign-changing d-wave gaps, none for sign-preserving |Δ|, and a clean threshold at Δ_s = Δ_c^s — is consistent with prior literature, so the central claim is not overturned. However, the analytic pillar fails an independent check. For the stated U = V τ_z, G0(k,0) = −(ε_k τ_z + Δ_k τ_x)/E_k², so G0(r,r,0) = −(ΣP)τ_z − (ΣQ)τ_x and Det[I − U G0(0)] = (1 + VΣP)² + (VΣQ)², not V²[(ΣP)² + (ΣQ)²]. The omitted 1 + 2VΣP terms are O(1) precisely where the resonance is claimed (V = 20 requires ΣP ≈ −0.05), so this is not a controlled approximation. The paper's mechanism — 'ΣP ≈ 0 and ΣQ ≈ 0 make A(0) vanish' — is thus incorrect; the exact condition couples the sign-change requirement (ΣQ ≈ 0) to a band-structure-dependent condition (ΣP ≈ −1/V), with the actual pole at finite ω where (V g_0(ω))² balances the residual. For nodal d-wave, the large pole strength C = Σ1/E² keeps that pole near zero over a wide doping range, so the d-wave conclusion is qualitatively robust; but for the fully gapped interlayer pairing the pole energy depends on ΣP' through |1 + VΣP'|/(VΣ_k 1/E_k²) in a way the paper never derives, weakening the Sec. IV claim that in-gap peaks 'generally exist' and leaving the La₃Ni₂O₇ two-peak prediction unquantified against doping and pocket asymmetry. This is an internal inconsistency between the stated determinant and the printed formula, not a disagreement with consensus; the paper's own acknowledged limitations (non-interacting model, point impurity, no disorder) are separate and do not bear on this point. The concrete test — evaluating both determinant expressions numerically and comparing with Fig. 1(b), 2(d), and 4(b) — settles it. The reader already identified this weakness ('identity and cross terms neglected'), and my independent derivation confirms and sharpens it, so the CONDITIONAL verdict stands, now with the explicit condition that Eqs. (8) and (12) be replaced by the exact determinant expansion and used to quantify the regime over which the sign-change signatures persist.","tokens_in":16043,"tokens_out":34434,"duration_ms":300876,"concrete_test":"Recompute the T-matrix denominator at ω = 0 for the single-band d-wave model with U = diag(V,−V): form G0(r,r,0) = −Σ_k (ε_k τ_z + Δ_k τ_x)/E_k² from Eq. (4), evaluate Det[I − U G0] = (1 + VΣP)² + (VΣQ)² at μ = 0.2 and μ = 2, and compare with Eq. (8) and with the Re A(ω) curves in Figs. 1(b) and 2(d). Then repeat for the interlayer model using Eqs. (12)–(15), locating the poles of the exact determinant; if the pole energies and the ΣP ≈ −1/V dependence differ from the paper's stated criterion, Eqs. (8) and (12) must be corrected and the claimed 'general existence' of the resonant peaks re-quantified over doping and V.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III defines A(ω) = Det[I − U G0(r,r,ω)] (Eq. 5) and then states, as the derivation of the central claim, that Re A(0) = V²[(ΣP)² + (ΣQ)²] (Eq. 8). Direct evaluation for the stated single-band Nambu problem gives a different result. With U = diag(V,−V) = V τ_z and G0(k,0) = −(ε_k τ_z + Δ_k τ_x)/(ε_k² + Δ_k²), the local Green's function is G0(r,r,0) = −(ΣP)τ_z − (ΣQ)τ_x, so Det[I − U G0(0)] = (1 + VΣP)² + (VΣQ)² = 1 + 2VΣP + V²[(ΣP)² + (ΣQ)²]. The identity and cross term are not negligible: the vanishing condition is ΣQ ≈ 0 together with ΣP ≈ −1/V (≈ −0.05 at V = 20), not ΣP ≈ 0 as the text claims. Away from that special ΣP value, the exact Re A(0) is O(1)–O(10²) even for a sign-changing gap with ΣQ = 0; the resonance, if present, then sits at the finite energy set by (V g_0(ω))² = (1 + VΣP)² + (VΣQ)², with g_0(ω) = Σ_k ω/(ω² − E_k²), a condition the paper never derives. The same defective structure appears in Eq. (12) for the interlayer case, so the Sec. IV statement that in-gap resonant peaks 'generally exist' for interlayer pairing is analytically unsupported, and the doping/model dependence of the La₃Ni₂O₇ two-peak prediction is not quantified. Because Eqs. (8) and (12) are the paper's analytical demonstration of the central claim, the derivation as printed is incorrect; the evidence for the qualitative link is then only the numerics (Figs. 1, 4, 5) and the Δ_s = Δ_c^s crossover (Fig. 2).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies single-impurity scattering in unconventional superconductors using the T-matrix formalism. It considers intralayer d-wave pairing on a square lattice and interlayer s-wave pairing in a bilayer model, and it claims to establish an analytical link between the sign change of the superconducting order parameter along the Fermi surface and the appearance of in-gap impurity resonances. The authors extend the analysis to a two-orbital bilayer model of La3Ni2O7 and predict two resonant peaks symmetric about the Fermi energy for interlayer pairing. The numerical part follows standard T-matrix methodology, but the analytical derivation of the central claim, especially Eqs. (8) and (12), contains a determinant identity error that undermines the stated proof.","tokens_in":16569,"tokens_out":8645,"duration_ms":86968,"significance":"If the central claim were established, the paper would provide a practical spectroscopic tool for distinguishing pairing symmetries in bilayer nickelates, and the two-peak prediction for interlayer pairing would be a concrete, falsifiable experimental signature. The numerical implementation is standard, the control case with an absolute-value gap is a useful test, and the realistic two-orbital model with self-consistent pairing is a valuable extension. However, the analytical centerpiece of the paper, Eq. (8) and its interlayer counterpart Eq. (12), is incorrect as written: the unit matrix in Det[I - U G0] is omitted. Because the claimed link between sign reversal and in-gap states is presented as a derived result rather than merely as a numerical observation, the paper's significance currently rests on an invalid derivation. The numerics may still support the qualitative conclusion, but the analytical argument must be corrected before the claim can be accepted.","major_comments":[{"comment":"The stated expression for A(0) is not the determinant defined in Eq. (5). For the single-band Nambu problem with U = V τ_z and G0(r,r,0) = -ΣP τ_z - ΣQ τ_x, where P_k = ε_k/(ε_k²+Δ_k²) and Q_k = Δ_k/(ε_k²+Δ_k²), direct evaluation gives Det[I - U G0(0)] = (1+VΣP)² + (VΣQ)². The paper's Eq. (8) omits the (1+VΣP)² term. Consequently, the statement in the text that small ΣP and small ΣQ imply Re A(0) ≈ 0 is incorrect; if ΣP = ΣQ = 0, the exact determinant equals 1. The pole condition is actually 1+VΣP = 0 together with ΣQ = 0, not V²[(ΣP)²+(ΣQ)²] = 0. This invalidates the analytical derivation of the claimed intrinsic link in Section III, including the interpretation of Fig. 1(b) and the conclusion in the abstract that in-gap states are directly linked to the sign reversal of the order parameter.","section":"Section III, Eq. (8)"},{"comment":"The same unit-matrix omission appears in the interlayer expression. For the 4x4 Nambu problem of the bilayer model, the exact determinant of I - U G0 contains terms of the form (1+...)² that are not captured by V²[(ΣP')²+(ΣQ')²]. Therefore the sentence immediately below Eq. (15) claiming that 'our analytical calculation indicates that in-gap resonant peaks generally exist' is not supported by the derivation. At most, the symmetry arguments show that ΣQ' can be small; they do not show that the full determinant vanishes. The prediction of two symmetric resonant peaks for interlayer pairing in Figs. 4 and 5 is thus supported only by the numerics, and the analytical rationale, as well as the dependence on V, t⊥, and chemical potential, must be re-derived with the correct determinant.","section":"Section IV, Eq. (12)"},{"comment":"The crossover at Δ_s = Δ_s^c is interpreted through the concavity of Re A(ω) and the disappearance of sign change along the Fermi surface. Since this interpretation relies on the incorrect Eq. (8), the analytical explanation of why in-gap features disappear at the critical value needs to be redone with the full determinant, including the frequency-dependent terms from G0(ω). The numerical observation of a crossover is not in question, but the claim that it establishes an 'intrinsic relationship' between sign reversal and in-gap states is not established by the presented derivation.","section":"Section III, Figs. 2 and 3"},{"comment":"The same incorrect form of A(0) is used for the p+ip case in Eq. (S1) of the Supplemental Material. Since the Supplemental explicitly repeats Eq. (8), the correction is not a local issue but affects the global analytical framework of the paper. Any revised version must propagate the corrected determinant through all analytical statements, including those in the Supplemental.","section":"Section III and Supplemental S-2"}],"minor_comments":[{"comment":"The section heading reads 'La2Ni3O7' while the rest of the paper, including the title and abstract, refers to La3Ni2O7. Please correct the compound formula.","section":"Section IV heading"},{"comment":"The notation A(ω) is used both for the full determinant and for its real part; Eq. (8) is written as if A(0) equals Re A(0), but the text later discusses Re A(ω) and Im A(ω) separately. Please define Re A(ω) and Im A(ω) explicitly and write the determinant in a way that makes the real and imaginary parts unambiguous.","section":"Section III, around Eq. (8)"},{"comment":"The critical s-wave component Δ_s^c = Δ_d(1-|μ/4|) is stated without derivation in the main text and attributed to Ref. [55]. Since this quantity plays a central role in the numerical crossover analysis, a short derivation or at least a clear statement of its origin in the main text would improve readability.","section":"Section III, Δ_s^c formula"},{"comment":"The definition of the impurity potential U in Eq. (5) uses U_{1+N_l,1+N_l} = -Vimp, which is clear for the single-band bilayer model. However, the realistic two-orbital model in Supplemental S-3 sets U11 = U22 and U55 = U66, implying the impurity acts on both orbitals of layer 1. Please state explicitly that this is the intended generalization and clarify whether the main-text notation N_l refers to layers or to orbitals in the two-orbital case.","section":"Section II, Eq. (5) and Supplemental S-3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses an active and important question, and the numerical T-matrix calculations are competently executed. The central problem is that Eqs. (8) and (12) are analytically wrong, and the paper's main conceptual claim is presented as a proof based on these equations. The error is correctable in principle: the authors should replace the incorrect determinant with the exact expression, re-examine the pole condition, and then either rephrase the claims as conditional on V and model parameters or demonstrate that the corrected condition holds for the parameter ranges considered. I do not see grounds for outright rejection, but the revision must be substantial because the analytical framework, the abstract, and the interpretation of the figures all depend on the defective equations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThe paper has a genuinely useful target—impurity spectroscopy for La3Ni2O7—and the numerical results for interlayer pairing are clean and interesting. But the analytic derivation that is supposed to establish the link between sign-changing gaps and in-gap states has a load-bearing error: Eq. (8) is not Det[I − U G0] as defined.\n\nWhat is actually new: the application to interlayer-paired nickelates, the prediction of two resonant peaks symmetric about the Fermi energy, and the control case with |Δ| (which cleanly shows no in-gap states for a sign-preserving gap). The T-matrix framework itself is standard, and the qualitative link was already known from cuprates and iron-based superconductors. The paper is honest about its limitations: non-interacting model, point impurity, no finite disorder.\n\nThe soft spot is the central analytic claim. For the stated Nambu U = diag(V, −V), direct algebra gives Det[I − U G0(0)] = (1+VΣP)^2 + (VΣQ)^2, not V^2[(ΣP)^2+(ΣQ)^2]. The dropped 1+2VΣP terms are not negligible: for V=20, the pole condition shifts from ΣP≈0 to ΣP≈−0.05, and away from that point the determinant is O(1)–O(100) even when ΣQ=0. So the paper does not prove, analytically, that a sign-changing gap guarantees a resonance. The same defective structure appears in Eq. (12) for interlayer pairing. The numerics in Figs. 1, 4, and 5 may still be correct, and the physical conclusion may survive, but the derivation as printed is incorrect.\n\nWho is this for? Groups doing STM on nickelates and theorists modeling impurity states. It deserves a serious referee, but the manuscript needs a major revision: fix the determinant expressions, show the exact expansion, and quantify the parameter ranges where the cancellation holds. If those are corrected, the paper could be a useful contribution.\n\nBest,\n[Name]","headline":"The numerics for interlayer pairing are promising, but the paper's central analytic derivation of the resonance condition has a miscomputed determinant and needs a major fix.","tokens_in":17124,"tokens_out":3051,"would_cite":false,"duration_ms":26715,"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":"Impurity-induced in-gap states occur exactly when a superconductor's order parameter changes sign along its Fermi surface, and the paper uses this to propose a spectroscopic test for pairing symmetry in La$_3$Ni$_2$O$_7$.","keywords":["unconventional superconductivity","impurity bound states","T-matrix theory","pairing symmetry","nickelate superconductors","interlayer pairing","scanning tunneling microscopy","sign-changing order parameter"],"falsifier":"Measure the local density of states next to a deliberately placed impurity in a fully gapped La$_3$Ni$_2$O$_7$ sample (film or crystal) by scanning tunneling spectroscopy. If the interlayer-pairing scenario is right, a strong scalar impurity should produce two sharp resonant peaks symmetric about the Fermi energy inside the gap; if instead only a single mid-gap resonance appears, or no in-gap resonance is observed, the specific cancellation mechanism proposed here would be ruled out.","tokens_in":15748,"feed_emoji":"🔬","tokens_out":9671,"duration_ms":80505,"temperature":0.7,"pith_summary":"The paper sets out to prove a direct link between the sign change of a superconducting order parameter along the Fermi surface and the appearance of bound states around a single nonmagnetic impurity. Working within the $T$-matrix description of impurity scattering, it shows that a sign-changing gap makes the impurity $T$-matrix denominator nearly vanish at zero energy, producing sharp in-gap resonances, while a sign-preserving gap with the same magnitude does not. It then applies the same logic to interlayer $s$-wave pairing in a bilayer: written in bonding and antibonding bands, that pairing is $s_\\pm$, so it too produces impurity resonances. For the bilayer nickelate La$_3$Ni$_2$O$_7$ with interlayer pairing, the predicted signature is two sharp peaks symmetrically placed about the Fermi energy. This gives a concrete spectroscopy-based criterion for distinguishing interlayer pairing from intralayer $d$-wave or $s_\\pm$ pairing in the nickelates.","feed_headline":"Two in-gap peaks would expose sign-changing pairing in nickelates","feed_subtitle":"Impurity spectra can distinguish interlayer pairing from d-wave pairing in La3Ni2O7.","key_machinery":"The load-bearing object is the impurity $T$-matrix denominator $A(\\omega)=\\det[\\hat I-\\hat U\\hat G_0(\\mathbf{r},\\mathbf{r},\\omega)]$: resonant in-gap states occur where its real and imaginary parts vanish at the same low energy. At zero energy it reduces to $A(0)=V_{\\rm imp}^2[(\\sum_{\\bf k}P_{\\bf k})^2+(\\sum_{\\bf k}Q_{\\bf k})^2]$ (and to the interlayer counterpart in Eqs. (12)-(15)), so the existence of a bound state is controlled by two weighted Brillouin-zone sums over the Fermi surface. The argument is that $P_{\\bf k}$ changes sign across the Fermi surface and largely cancels, while $Q_{\\bf k}$ cancels whenever the gap sign reverses under a symmetry of the band structure. That cancellation is the mechanism: sign-changing order parameters suppress the denominator and create the resonances; sign-preserving order parameters leave the denominator finite.","core_discovery":"On the paper's own terms, the central discovery is that the impurity-induced bound state is not a generic property of a particular gap shape but is tied specifically to the order parameter changing sign as one moves along the normal-state Fermi surface. For intralayer pairing, the zero-energy real part of the $T$-matrix denominator is $A(0)=V_{\\rm imp}^2[(\\sum_{\\bf k}P_{\\bf k})^2+(\\sum_{\\bf k}Q_{\\bf k})^2]$, with $P_{\\bf k}=\\varepsilon_{\\bf k}/(\\varepsilon_{\\bf k}^2+\\Delta_{\\bf k}^2)$ and $Q_{\\bf k}=\\Delta_{\\bf k}/(\\varepsilon_{\\bf k}^2+\\Delta_{\\bf k}^2)$. $P_{\\bf k}$ cancels between hole and electron regions and $Q_{\\bf k}$ cancels under the symmetry that reverses the gap sign, so a sign-changing $d$-wave gap gives $A(0)\\simeq0$ and a mid-gap resonance; replacing $\\Delta_{\\bf k}$ by its absolute value removes the sign change and removes the in-gap feature. Adding an $s$-wave component that shifts the gap nodes off the Fermi surface continuously weakens and then eliminates the resonances at a critical value $\\Delta^c_s=\\Delta_d(1-|\\mu|/4)$. For interlayer $s$-wave pairing, transforming to bonding and antibonding bands makes the gap opposite in sign on the two Fermi pockets, and the analogous cancellation of $Q'_{\\bf k}$ yields two resonant peaks symmetric about zero energy; the same two-peak signature survives in a two-orbital model of La$_3$Ni$_2$O$_7$ with self-consistently determined interlayer pairing.","pith_inferences":["The zero-denominator condition suggests the probe is not limited to the two cases computed here: any sign-changing gap whose Fermi-surface sums cancel should show impurity resonances, and the $p+ip$ example in the paper's supplemental material is a sign that this is a symmetry-independent rule.","A testable extension: the strength and peak position of the two symmetric resonances should depend on the relative weights of the bonding and antibonding Fermi pockets, so doping or pressure that unbalances the pockets is a way to turn the predicted signature off and verify the mechanism.","With ambient-pressure thin films now available, one can look for the two-peak signature directly; a null result in a fully gapped film would be evidence against dominant interlayer pairing of the type modeled here."],"forward_implications":["For intralayer $d_{x^2-y^2}$ pairing near half-filling, a nonmagnetic impurity produces a sharp mid-gap resonant peak, and the same gap with its sign made positive everywhere produces no in-gap structure.","As an $s$-wave component is added to a $d$-wave gap, the mid-gap peak splits and weakens, and it disappears entirely once the $s$-wave component reaches the value at which the gap no longer changes sign along the Fermi surface.","In a bilayer with dominant interlayer $s$-wave pairing, the impurity response consists of two sharp resonant peaks positioned symmetrically about the Fermi energy, because the band-basis gap is $s_\\pm$.","In a two-orbital model of La$_3$Ni$_2$O$_7$ with self-consistent interlayer pairing, the two symmetric peaks survive, so STM on a single impurity can distinguish interlayer pairing from intralayer $d$-wave (one mid-gap peak) or intralayer $s_\\pm$ (weaker in-gap states)."],"supporting_citations":[{"why":"supplies the T-matrix formalism for impurity scattering used throughout the paper.","marker":"[2]"},{"why":"originates the midgap-state argument for sign-changing d-wave pairing that the paper re-derives in T-matrix language.","marker":"[38]"},{"why":"attributes impurity bound states to Andreev scattering, the qualitative picture the paper makes quantitative.","marker":"[5]"},{"why":"provides the two-orbital bilayer tight-binding parameters used for the La3Ni2O7 model.","marker":"[40]"},{"why":"establishes that interlayer interaction gives s± pairing, the pairing state the impurity probe targets.","marker":"[54]"},{"why":"contains the computed impurity spectra for intralayer d-wave and s± pairing that the paper contrasts with the interlayer result.","marker":"[23]"},{"why":"gives the critical s-wave component formula and the 8x8 two-orbital model details needed for the nickelate calculation.","marker":"[55]"}],"fun_headline_variants":["Sign-changing gaps leave distinct impurity-state fingerprints","Two impurity peaks reveal interlayer pairing in La3Ni2O7","Impurity resonances distinguish d-wave from interlayer pairing","Nickelate impurity spectra betray sign-changing order parameter"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument would collapse if the cancellation that makes $A(0)$ small is spoiled by unequal weights of the two Fermi pockets, by a Fermi surface on which the gap nodes do not cross, or by a non-negligible imaginary part of $A(\\omega)$ inside the gap.","fun_headline_variants_meta":{"raw":{"variants":["Sign-changing gaps leave distinct impurity-state fingerprints","Two impurity peaks reveal interlayer pairing in La3Ni2O7","Impurity resonances distinguish d-wave from interlayer pairing","Nickelate impurity spectra betray sign-changing order parameter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000288,"raw_usage":{"total_tokens":1819,"prompt_tokens":1201,"completion_tokens":618,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":817,"completion_tokens_details":{"reasoning_tokens":552}},"tokens_in":817,"tokens_out":618,"duration_ms":6919,"temperature":1.0,"reasoning_tokens":552,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:35:18.748190+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the local density of states next to a deliberately placed impurity in a fully gapped La$_3$Ni$_2$O$_7$ sample (film or crystal) by scanning tunneling spectroscopy. If the interlayer-pairing scenario is right, a strong scalar impurity should produce two sharp resonant peaks symmetric about the Fermi energy inside the gap; if instead only a single mid-gap resonance appears, or no in-gap resonance is observed, the specific cancellation mechanism proposed here would be ruled out.","supporting_citations":[{"cited_title":"Hu, Midgap surface states as a novel signature for d 2 xa-x 2 b -wave superconductivity, Phys","cited_arxiv_id":null,"evidence_quote":"originates the midgap-state argument for sign-changing d-wave pairing that the paper re-derives in T-matrix language."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"attributes impurity bound states to Andreev scattering, the qualitative picture the paper makes quantitative."},{"cited_title":"Effective perpendicular electric field as a probe for interlayer pairing in ambient-pressure superconducting La$_{2.85}$Pr$_{0.15}$Ni$_{2}$O$_{7}$ thin films","cited_arxiv_id":"2503.23861","evidence_quote":"establishes that interlayer interaction gives s± pairing, the pairing state the impurity probe targets."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the critical s-wave component formula and the 8x8 two-orbital model details needed for the nickelate calculation."}],"review_version":1}