{"id":"9bf590e7-32cc-4f52-8cb4-e50017428a01","arxiv_id":"2411.12733","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Proximity-induced pair correlations in 1D superconductor-normal nanowires are computed with Keldysh NEGF, giving algebraic decay in clean wires, a disorder-driven crossover to exponential decay, and a spectral explanation of resonant Cooper pair injection.","lead":"Using a non-equilibrium Green's function method, this paper computes how superconducting pairing leaks into a normal nanowire segment, showing a slow algebraic decay in clean wires and a switch to exponential decay when disorder is added. It then uses the same correlations to explain the enhanced Cooper pair injection seen in a recent resonant-tunneling experiment.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Rigid step-function Δ with no feedback from F(x) contradicts the paper's self-consistency claim; quantitative ξSN values and the disorder crossover may shift under a self-consistent update.","rationale":"The reader identified the rigid, non-self-consistent pair potential as the weakest assumption; this is indeed the most load-bearing concern. The paper's own assertion that self-consistency is imperative for the correct decay of the pair amplitude makes the issue central: if the NEGF method is not self-consistent, its agreement with self-consistent BdG is only qualitative, and the quantitative outputs (ξSN, its scalings, and the disorder crossover) are not guaranteed. The disorder crossover is additionally supported by only a few points and is explicitly conjectural, but that is a secondary evidence gap. A direct self-consistent BdG comparison is the decisive test: it settles whether the rigid-Δ approximation changes the extracted parameters enough to affect the paper's quantitative claims. Since the central framework remains a valid calculation for the fixed-Δ model, the qualitative insights likely survive, so a CONDITIONAL verdict is appropriate and my analysis does not alter it.","tokens_in":16987,"tokens_out":11277,"duration_ms":116672,"concrete_test":"Run a self-consistent BdG calculation on the same 1D lattice (a=2 nm, t0=415 meV, μ=10 meV, Δ=1 meV) where Δ(x)=g(x)F(x) is iterated to convergence, with g(x) chosen so the bulk Δ matches 1 meV. Compare the resulting F(x) in the normal region and the fitted ξSN with Fig. 1(c) and the scaling data in Fig. 2. If ξSN shifts by more than ~20% or the power-law form changes, the rigid-Δ assumption is quantitatively significant and the extracted coherence lengths require revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing assumption is that the pair potential Δ(x)=ΔΘ(-x) in Eq. (1) is a rigid input. The text in Sec. III A states that 'The NEGF method implicitly takes this into account through the self-consistent calculation of the Green's functions of the bath and subsequently of the system.' But that self-consistency refers only to iteratively computing the lead self-energies for a fixed Δ; the computed F(x) is never fed back to update Δ(x). Thus the inverse proximity effect—suppression of the order parameter near the interface and its backaction on Andreev reflection—is omitted. Since the central quantitative results (F0, ξSN, the μ/Δ scalings in Fig. 2, and the disorder crossover in Sec. III D) are obtained with this rigid potential, they could be quantitatively modified if self-consistency is included. Moreover, the disorder crossover itself rests on only four (σD, γ) data points in Fig. 5 and is explicitly conjectured ('we conjecture that this crossover marks an order-disorder phase transition'), so the evidence is thin even within the fixed-Δ model. The benchmark with Ref. [63] is qualitative (same power-law form) and does not validate the specific coefficients.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a Keldysh non-equilibrium Green's function (NEGF) framework to compute the induced pair amplitude F(x) in one-dimensional superconductor-normal hybrids with a rigid step-function pair potential Δ(x)=ΔΘ(-x). The central claims are: (i) F(x) decays algebraically, F(x)=F0 ξSN/(x+ξSN), in clean wires, with an effective coherence length ξSN that grows with μ/Δ and decreases with Δ; (ii) an experimentally motivated disorder profile produces a crossover from power-law to exponential decay of F(x); and (iii) the spectral decomposition of the induced pair amplitude explains the enhanced Cooper-pair injection in the resonant device of Ref. [75]. The paper also discusses the Andreev approximation, the local density of states, and the inverse proximity effect in a qualitative way.","tokens_in":17285,"tokens_out":5101,"duration_ms":51469,"significance":"If the framework is correct, it offers a useful and comparatively general NEGF-based tool for studying proximity-induced correlations in hybrid devices, going beyond the restricted geometries of typical self-consistent BdG calculations and connecting correlations to transport observables. The spectral-resolution analysis of the resonant Cooper-pair injector is a constructive, falsifiable application. The numerical method is based on standard NEGF equations, and the power-law decay of F(x) agrees qualitatively with the self-consistent BdG results of Ref. [63]. However, the paper's claims are partially stronger than its evidence: the pair potential is not computed self-consistently despite the text saying otherwise, the scaling laws in Fig. 2 are phenomenological fits rather than derivations, and the disorder crossover in Sec. III D rests on a small number of points and an explicit conjecture. These issues do not invalidate the method, but they require a clearly revised presentation of what is derived versus what is fitted or conjectured.","major_comments":[{"comment":"The text states that 'The NEGF method implicitly takes this into account through the self-consistent calculation of the Green's functions of the bath and subsequently of the system,' but this is not order-parameter self-consistency. Eq. (1) fixes Δ(x)=ΔΘ(-x) as an input, and the computed pair amplitude F(x) is never fed back to update Δ(x). The inverse proximity effect, i.e., suppression of the superconducting order parameter near the interface and its backaction on Andreev reflection, is therefore not captured. Since F0, ξSN, and the scalings in Fig. 2 are all extracted within this fixed-Δ model, their quantitative values could change under a self-consistent treatment. Please either implement an order-parameter self-consistency loop (or a controlled estimate of its effect) or revise the language in Sec. III A to state clearly that all results are obtained with a rigid, non-self-consistent pair potential.","section":"Sec. III A, Eq. (1)"},{"comment":"The decay form F(x)=F0 ξSN/(x+ξSN) is assumed a priori, and ξSN is obtained by fitting to this form rather than derived from the model. The reported ξSN(μ/Δ)=6.8√(μ/Δ)+28.7 nm and ξSN(Δ)=15.1+26.2/Δ-1.5/Δ² are phenomenological fits, as the authors note, but the surrounding discussion draws conclusions such as 'the growth of the effective coherence length is ∼√(μ/Δ)' as though this were a model prediction. The benchmark with Ref. [63] validates the power-law functional form qualitatively, but it does not validate these specific coefficients or the μ/Δ dependence. Please state explicitly which results are analytic or structural and which are numerical fits, and provide fit-quality measures or uncertainties for the extracted ξSN values.","section":"Eq. (11), Fig. 2"},{"comment":"The claimed disorder-induced crossover from power-law to exponential decay is supported by only four (σD,γ) data points in Fig. 5(b), and the text itself labels the associated order-disorder transition a conjecture. No error bars, confidence intervals for the fitted decay forms, statistical comparison between power-law and exponential fits, or finite-size scaling are given, although the text acknowledges that finite-size scaling is difficult. As presented, the crossover is a suggestive model prediction rather than an established result. Please add quantitative evidence for the crossover, such as fit residuals, multiple system sizes, or a well-defined observable that sharply identifies the transition, or explicitly present this part as a preliminary finding that requires further numerical study.","section":"Sec. III D, Fig. 5"}],"minor_comments":[{"comment":"The caption labels the top row (a-b) and bottom row (c-d), but the text refers to Fig. 2(d) for ξSN versus Δ, which appears to be panel (e) in the caption. Please harmonize the panel labels.","section":"Fig. 2"},{"comment":"The text lists μ/Δ ∈ {1.2,3,5,10,25}, while the Fig. 2 caption lists {1.2,3,5,10,15}. Please correct the discrepancy.","section":"Sec. III B, Fig. 2"},{"comment":"The phrase 'self-consistent calculation of the Green's functions' is misleading and should be replaced with 'iteratively converged Green's functions' or similar, reserving 'self-consistent' for treatments in which Δ is updated using F(x).","section":"Sec. III A"},{"comment":"The vertical dotted line in Fig. 5(c) is described as denoting 'the variance of the disorder profile aσD,' but Eq. (13) defines σ(x) through a standard deviation aσD; please clarify whether the line marks the standard deviation or the variance and use consistent notation.","section":"Sec. III D, Eq. (13)"},{"comment":"The sentence 'instead goes to zero in the opposite µ/∆ → 0 limit [77]' has a dangling bracket; the citation should be placed after the full clause.","section":"Sec. III B"},{"comment":"The manuscript does not include a data or code availability statement. Given that the fitting coefficients in Fig. 2 and the disorder averages in Fig. 5 are central quantitative outputs, a reproducibility statement would strengthen the paper.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of cond-mat.mes-hall and presents a useful numerical framework. The main risk is the overstatement of self-consistency and the thinness of the disorder-crossover evidence; these are fixable in revision. I would recommend inviting a revision in which the authors either implement an order-parameter self-consistent scheme or explicitly frame all results as fixed-Δ calculations, and in which the fitted scalings and the conjectured crossover are presented with appropriate caveats and statistical support."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper gives a Keldysh-NEGF recipe for computing induced pair amplitudes in 1D SN hybrids and applies it in three places: clean-wire decay, spectral decomposition via the correlation Green's function, and a published resonant Cooper-pair injector. The clean-wire power law was already in Rai et al. (Ref. [63]), and the authors say so; the benchmark is qualitative (same functional form), not coefficient-level. The credible new material is Sec. III C and III E: resolving F(x) by energy shows which spectral window carries the proximity correlations and how interference among those windows produces the net decay; the injector analysis ties the enhanced sub-gap conductance to a resonant enhancement of the pair amplitude inside the quantum well. That is a useful way to think about transport in layered devices.\n\nThe soft spots are real but not fatal. The stress-test note is right: the word 'self-consistent' is doing too much work. Delta is a rigid step input; the NEGF bath self-energies are iterated, but F(x) is never fed back into Delta. The paper says the method 'takes this into account' and claims to capture backaction from the normal region on the superconductor. It does not, because Delta(x) is static. That claim should be dialed back. Second, the disorder crossover in Sec. III D rests on four (sigma_D, gamma) points and is explicitly conjectured to be a phase transition. Fine as a motivated guess, but it is not established, and the text should say so more plainly. Third, the xi_SN(mu/Delta) and xi_SN(Delta) scalings are fits with phenomenological coefficients; the authors admit this, which is good, but the results should not be read as derivations. The citation pattern is solid: Ref. [63] is correctly credited and benchmarked, and self-citation is not an issue here.\n\nNo code is provided, so the numerics are unverified, but the method description is detailed enough to attempt a reproduction. The central modeling assumption, rigid Delta, is the conventional first step in this literature, and the NEGF pipeline itself is sound. If inverse proximity effects are significant, the quantitative xi_SN values could shift, though the power-law form is not likely to change.\n\nWho is this for? NEGF toolbuilders and experimental groups interpreting sub-gap conductance in superconductor-semiconductor heterostructures. It is not a breakthrough, but it is a competent methods paper with a credible new decomposition and a clear device application. A serious editor should send this to peer review rather than desk reject it; I would accept with revisions, mainly asking the authors to soften the self-consistency claim, reframe the disorder crossover as suggestive rather than established, and ideally release code or data.","headline":"A workmanlike NEGF framework for induced pair amplitudes; the spectral-resolution and device-application sections are the real contributions, while the self-consistency language and the disorder crossover need scaling back.","tokens_in":17816,"tokens_out":2469,"would_cite":true,"duration_ms":26862,"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":"A Keldysh non-equilibrium Green's function method computes the induced pair amplitude in one-dimensional superconductor-normal hybrids, yielding an algebraic decay law in clean wires, a disorder-driven crossover to exponential decay, and…","keywords":["proximity effect","pair amplitude","Keldysh non-equilibrium Green's functions","superconductor-normal hybrid","Andreev approximation","disorder-induced crossover","Cooper pair injection","one-dimensional nanowire"],"falsifier":"A direct comparison calculation would settle it: solve the same one-dimensional SN junction with a self-consistent BdG scheme in which $\\Delta(x)=g(x)F(x)$ is updated using the computed pair amplitude, at $\\mu=10$ meV and $\\Delta=1$ meV. If the self-consistent pair amplitude does not follow $F(x)=F_0\\xi_{SN}/(x+\\xi_{SN})$ with $\\xi_{SN}\\approx44$ nm, or if the disorder crossover line shifts, the rigid-order-parameter assumption is the reason.","tokens_in":16764,"feed_emoji":"🔬","tokens_out":9076,"duration_ms":79716,"temperature":0.7,"pith_summary":"The paper sets out to show that the Keldysh non-equilibrium Green's function method, through the correlation Green's function $G^n(E)$, is a working microscope for proximity-induced superconductivity in one-dimensional superconductor-normal hybrids. In a clean wire the induced pair amplitude $F(x)=\\langle c^\\dagger_{\\uparrow x}c^\\dagger_{\\downarrow x}\\rangle$ decays algebraically, $F(x)=F_0\\,\\xi_{SN}/(x+\\xi_{SN})$, rather than exponentially, and the extracted length $\\xi_{SN}$ grows with $\\mu/\\Delta$ and shrinks as $1/\\Delta$. The paper also finds that a localized disorder profile at the junction turns this algebraic tail into an exponential one, and that the spectral resolution of $F(x)$ explains why a resonant tunneling device injects Cooper pairs efficiently. If correct, the framework gives device designers a direct way to compute and engineer induced correlations in layered superconducting hybrids.","feed_headline":"Proximity decay in clean nanowires is algebraic, not exponential","feed_subtitle":"A Green's-function method resolves how Cooper-pair correlations leak into normal wires and why disorder flips the decay.","key_machinery":"The load-bearing object is the correlation Green's function $G^n(E)$ defined in Eq. (9), a Nambu-space two-point correlator whose energy integral gives $-i\\langle\\hat\\psi^\\dagger_i\\hat\\psi_j\\rangle$ and hence the pair amplitude $F(x)$ at each site. It is built from the retarded Green's function and the contact in-scattering functions, so it connects the equilibrium many-body correlations to the same contact self-energies used in transport calculations. What it does for the argument is to decompose the induced pairing by energy and position: the interference of the spectrally resolved oscillating correlations is what produces the algebraic tail in clean wires, the exponential tail under disorder, and the resonant enhancement in the Cooper-pair injector.","core_discovery":"On its own terms, the paper's central discovery is that the correlation Green's function $G^n(E)=G^r(E)[\\Sigma^{in}_L(E)+\\Sigma^{in}_R(E)]G^a(E)$ carries the full position- and energy-resolved information about induced pairing. Integrating it over energy recovers the pair amplitude $F(x)$, and resolving it by energy shows that sub-gap states ($|E|<\\Delta$) are the carriers of proximity correlations, while energies far above the chemical potential contribute nothing. The same object yields the clean-wire algebraic decay law $F(x)=F_0\\,\\xi_{SN}/(x+\\xi_{SN})$, with $\\xi_{SN}$ behaving as $\\sim\\sqrt{\\mu/\\Delta}$ in the Andreev regime and $\\sim 1/\\Delta$ as the order parameter grows, and it produces the disorder-driven crossover from power-law to exponential decay. Applied to the resonant Cooper-pair injection device, the calculation shows that resonant tuning makes the in-gap spectral pair amplitude positive and resonantly enhanced inside the quantum well, which is the microscopic counterpart of the enhanced sub-gap conductance seen in transport.","pith_inferences":["The rigid step-function $\\Delta$ means the calculation omits the inverse proximity effect; a natural extension is to update $\\Delta(x)=g(x)F(x)$ self-consistently and check whether the algebraic exponent and $\\xi_{SN}$ values survive.","The conjectured order-disorder crossover, with $f_c/f_d\\sim\\log L$ separating clean and disordered correlation patterns, could be tested by finite-size scaling of $F(x)$ on much longer wires, since the current evidence rests on four data points.","The non-negativity of the in-gap spectral pair amplitude in the resonant case suggests an energy-resolved experimental probe: tunneling spectroscopy into the quantum well might see the same resonance directly, without relying on integrated conductance.","If the method generalizes as claimed, it supplies a ready design loop for Majorana nanowires and other topological hybrids, where optimizing the induced pair amplitude and its decay length in the semiconductor is the central engineering problem."],"forward_implications":["In clean one-dimensional superconductor-normal wires at zero temperature, induced correlations penetrate with an algebraic tail $F(x)=F_0\\xi_{SN}/(x+\\xi_{SN})$; the effective coherence length $\\xi_{SN}$ grows roughly as $\\sqrt{\\mu/\\Delta}$ and falls as $1/\\Delta$.","Disorder localized near the interface, characterized by strength $10^{-\\gamma}$ and spatial spread $\\sigma_D$, drives a crossover from algebraic to exponential decay of $F(x)$, so the disorder profile must be controlled to engineer reproducible proximity devices.","The spectral resolution shows that only energies inside the gap $|E|<\\Delta$ carry the proximity correlations; the integral over energy in Eq. (10) can be truncated in practice.","In the resonant Cooper-pair injection device, the enhanced sub-gap conductance has a microscopic counterpart: a resonantly enhanced, non-negative in-gap pair amplitude inside the quantum well, offering a correlation-based design target for injection efficiency.","Because the method is built on contact self-energies, it extends to transport setups, finite temperatures, and higher-dimensional lattices with other pairing symmetries."],"supporting_citations":[{"why":"Supplies the self-consistent BdG benchmark showing power-law decay of the order parameter in clean nanowires, which the NEGF result reproduces.","marker":"[63]"},{"why":"The experimental report of enhanced Cooper-pair injection via resonant tunneling that the paper's pair-amplitude analysis explains.","marker":"[75]"},{"why":"The earlier BTK-based theoretical model of the same device, which the paper extends by computing induced pair correlations.","marker":"[84]"},{"why":"Defines the Andreev approximation whose regime boundary the paper explores by varying $\\mu/\\Delta$.","marker":"[76]"},{"why":"Provides the analytic transport picture of sub-gap conductance with arbitrary carrier density that the correlation results are said to align with.","marker":"[77]"},{"why":"Anderson's theorem for bulk superconductors, cited to motivate why disorder can still affect proximity-induced pair amplitudes.","marker":"[80]"},{"why":"Prior work on disorder-robust odd-frequency p-wave pairing in normal-metal/superconductor junctions, mentioned as a future direction for the framework.","marker":"[64]"},{"why":"Earlier pair-amplitude dynamics in superconductor-quantum dot hybrids, part of the context for using pair amplitude as a diagnostic.","marker":"[62]"}],"fun_headline_variants":["Disorder flips proximity decay from algebraic to exponential","Clean nanowires show algebraic, not exponential, proximity decay","Green's function method reveals how Cooper pairs leak into normal wires","Resonant tuning enhances sub-gap Cooper-pair injection","Microscopic theory explains enhanced Cooper-pair injection"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the superconducting order parameter is a rigid step function $\\Delta(x)=\\Delta\\Theta(-x)$ fixed before the calculation, so the computed pair amplitude is never fed back to suppress $\\Delta$ near the interface; if the inverse proximity effect is strong, the extracted decay law and $\\xi_{SN}$ would change.","fun_headline_variants_meta":{"raw":{"variants":["Disorder flips proximity decay from algebraic to exponential","Clean nanowires show algebraic, not exponential, proximity decay","Green's function method reveals how Cooper pairs leak into normal wires","Resonant tuning enhances sub-gap Cooper-pair injection","Microscopic theory explains enhanced Cooper-pair injection"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00106,"raw_usage":{"total_tokens":4446,"prompt_tokens":946,"completion_tokens":3500,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":3419}},"tokens_in":562,"tokens_out":3500,"duration_ms":25830,"temperature":1.0,"reasoning_tokens":3419,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:12:26.590872+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct comparison calculation would settle it: solve the same one-dimensional SN junction with a self-consistent BdG scheme in which $\\Delta(x)=g(x)F(x)$ is updated using the computed pair amplitude, at $\\mu=10$ meV and $\\Delta=1$ meV. If the self-consistent pair amplitude does not follow $F(x)=F_0\\xi_{SN}/(x+\\xi_{SN})$ with $\\xi_{SN}\\approx44$ nm, or if the disorder crossover line shifts, the rigid-order-parameter assumption is the reason.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the self-consistent BdG benchmark showing power-law decay of the order parameter in clean nanowires, which the NEGF result reproduces."},{"cited_title":"Bouscher, D","cited_arxiv_id":null,"evidence_quote":"The experimental report of enhanced Cooper-pair injection via resonant tunneling that the paper's pair-amplitude analysis explains."},{"cited_title":"Bouscher, R","cited_arxiv_id":null,"evidence_quote":"The earlier BTK-based theoretical model of the same device, which the paper extends by computing induced pair correlations."},{"cited_title":"Andreev, Thermal conductivity of the intermediate state in superconductors, Sov","cited_arxiv_id":null,"evidence_quote":"Defines the Andreev approximation whose regime boundary the paper explores by varying $\\mu/\\Delta$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Anderson's theorem for bulk superconductors, cited to motivate why disorder can still affect proximity-induced pair amplitudes."},{"cited_title":"L¨ othman, C","cited_arxiv_id":null,"evidence_quote":"Prior work on disorder-robust odd-frequency p-wave pairing in normal-metal/superconductor junctions, mentioned as a future direction for the framework."},{"cited_title":"Heckschen and B","cited_arxiv_id":null,"evidence_quote":"Earlier pair-amplitude dynamics in superconductor-quantum dot hybrids, part of the context for using pair amplitude as a diagnostic."}],"review_version":1}