{"id":"d3f5e1c9-9e87-4ce4-8131-fedb4d164fba","arxiv_id":"1908.07953","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A circuit-theory extension of the Eilenberger formalism computes Josephson current and density of states for disordered multiband superconductors with coexisting spin-density-wave order.","lead":"This paper builds a quasiclassical theory of Josephson junctions between two-band superconductors that host both superconducting and spin-density-wave order, including disorder and arbitrary interface transparency. It derives circuit-theory-like formulas for supercurrent and local density of states that could be used to interpret experiments on iron-based superconductors.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central current trace formula, Eq. (27), drops the spatial phase-gradient term with only a typo-level justification; this unquantified assumption underlies the paper's Josephson results.","rationale":"I agree with the reader that the phase-gradient assumption is the weakest load-bearing step. The central claim that circuit theory follows from the Yip auxiliary-solution method and reproduces the Josephson current from bulk propagators is only as secure as Eq. (27), and Eq. (27) is entered with a sentence that appears to contain a typo and no quantitative justification. This is particularly important because the same trace enters the compact formula Eq. (45) and the claimed Jc ~ Δ²/(eRN M) suppression. At the same time, I would not reject the paper: the auxiliary-solution method is a known technique, the paper checks known limits such as Ambegaokar-Baratoff and Kulik-Omelyanchuk behavior, and the numerical self-consistency in Section IV provides independent evidence that the underlying Eilenberger machinery is being implemented. The correct response is to require the authors to verify or quantify the dropped phase-gradient term before the central trace formula is accepted. This does not change the reader's conditional verdict, so I recommend UNCHANGED.","tokens_in":18475,"tokens_out":19266,"duration_ms":193270,"concrete_test":"Solve the full Eilenberger equations, Eq. (12), directly for a complex order parameter Δ(x)e^{iχ(x)} without applying the unitary transformation Eq. (26), for a representative high-transparency junction (e.g. D=0.9) in the SC+SDW coexistence regime. Compute the interface trace Tr[τ3ρ3σ0Ĝa(0)] and the resulting J(χ), and compare with Eq. (27). If the difference is negligible over the plotted range, the dropped phase-gradient term is harmless; if the difference is visible on the scale of Fig. 6, Eq. (27) and all derived circuit-theory formulas inherit a missing contribution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim rests on the trace formula Eq. (27), which is obtained after applying the unitary transformation Eq. (26) and then dropping the ∂χ/∂x term. The text justifies this by saying that the term 'can be ignored for one does expect the phase to vary substantially across the junction', which presumably should read 'does not expect'. No estimate is provided for the size of the neglected term. In a current-carrying Josephson junction the phase is not strictly uniform in the electrodes: ∂χ/∂x is tied to the supercurrent density and can be sizeable near a high-transparency constriction, over lengths comparable to the coherence length. Because Eq. (27) is used to produce the current-phase relation in Fig. 6, the Ic(V) results in Fig. 7, and the circuit-theory traces in Eqs. (42)-(45), an unquantified omission of this term is the most load-bearing step in the paper. This is not an internal inconsistency, and the reproduction of known limits (e.g. Ambegaokar-Baratoff) gives confidence, but the assumption should be checked rather than asserted.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a quasiclassical Eilenberger theory for Josephson junctions between disordered multiband superconductors with coexisting s± superconducting and spin-density-wave order. The authors use the Yip auxiliary-solution method to satisfy the nonlinear Zaitsev boundary conditions, solve the Eilenberger equations self-consistently for the spatial profiles of the superconducting and magnetic order parameters, and compute the local density of states and the Josephson current as functions of phase and voltage. They then connect their trace formula to Nazarov circuit theory, give averaged expressions in terms of transmission-eigenvalue distributions, and present a compact current-phase relation for the SDW-coexisting case. The central claim is that the Josephson physics of such junctions is determined by bulk quasiclassical propagators plus the transmission distribution, with a universal suppression Jc ~ Δ²/(e RN M) when magnetic order dominates.","tokens_in":18666,"tokens_out":4406,"duration_ms":46176,"significance":"If the central trace formula and its circuit-theory interpretation are correct, this is a useful and fairly general result: it extends Nazarov's circuit theory beyond single-band systems to multiband superconductors with competing orders, and it provides a common framework that reproduces known limits such as Ambegaokar-Baratoff and Kulik-Omelyanchuk. The paper's strengths include a fully self-consistent numerical solution of the Eilenberger equations with nonlinear boundary conditions, explicit phase-diagram and density-of-states computations, and a clear demonstration of how transmission averaging reproduces several previously known current-phase relations. The main weakness is that the central trace formula is asserted rather than derived, and the approximation under which the spatial phase gradient is neglected is not quantified; these issues carry over into all Josephson-current results and the final averaged formulas.","major_comments":[{"comment":"The central current trace formula is introduced as the result of \"somewhat lengthy calculation\" but is not actually derived in the main text or in the appendices. Appendices A–D construct the auxiliary solutions and bulk propagators and state the matching coefficients, but they stop short of showing how these ingredients produce the trace formula Eq. (27). Because Eq. (27) is the starting point for all Josephson-current results in Figs. 6–7 and for Eqs. (42)–(45), the derivation must be supplied in a revised version, either in the main text or in a dedicated appendix with the key intermediate steps.","section":"Section IV.B, Eq. (27)"},{"comment":"The unitary transformation Eq. (26) introduces a term proportional to ∂χ/∂x in the Eilenberger equation, and the paper drops this term with the sentence \"This term, however, can be ignored for one does expect the phase to vary substantially across the junction.\" The sentence appears to contain a typo (presumably \"does not expect\"), but more importantly no estimate is given for the size of the neglected term. In a current-carrying Josephson junction the phase gradient in the electrodes is tied to the supercurrent and is not obviously negligible near a high-transparency interface; since Eq. (27) is obtained after this omission, this unquantified approximation is load-bearing. The authors should justify the short-junction limit explicitly, estimate (ℏvF/2)∂χ/∂x relative to Δ and M, or present a benchmark calculation that retains the phase-gradient term.","section":"Section IV.B, Eqs. (25)–(27)"},{"comment":"The compact current-phase relation Eq. (45) contains an unspecified prefactor cp described only as \"of the order of one.\" Since Eq. (45) is presented as a closed-form result and is used for the scaling claim Jc ~ Δ²/(e RN M), leaving cp undetermined turns a supposedly quantitative formula into a parametric statement. The authors should evaluate cp explicitly from the transmission-average integral, or at least provide the integral representation from which cp follows, so that the result is reproducible and falsifiable.","section":"Section V, Eq. (45)"}],"minor_comments":[{"comment":"The sentence \"This term, however, can be ignored for one does expect the phase to vary substantially across the junction\" appears to contain a typo; it should presumably read \"does not expect.\" As written it asserts the opposite of the intended approximation.","section":"Section IV.B, after Eq. (26)"},{"comment":"The sentence \"the quasiclassical functions which account for the magnetic order do not explicitly enter into the expression for the Josephson current\" is misleading, because M enters explicitly in Eqs. (42)–(44) and also through the bulk propagators in Eq. (27). The intended statement should be rephrased, for example that M enters only through the bulk Green's functions and not through an additional interface-specific magnetic term.","section":"Section IV.B, after Eq. (27)"},{"comment":"The notation SSDWINISSDW is used without a definition or a diagram; a brief explanation of this junction type would help the reader, especially since it is the basis for Eq. (45).","section":"Section V, Eq. (42)–(45)"},{"comment":"The two non-BCS curves in Fig. 7 are not accompanied by the disorder parameters Γ0 and Γπ used in the calculation; specifying these parameters would make the numerical results reproducible.","section":"Fig. 7"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be a solid contribution if the derivation of Eq. (27) is supplied and the phase-gradient approximation is quantified; as written, the central Josephson results rest on an unshown calculation and an unquantified approximation. The unspecified prefactor cp in Eq. (45) should also be removed or computed. These are fixable within the manuscript's scope, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is a working treatment of a Josephson junction between two-band s± superconductors with coexisting SDW order, using Yip's auxiliary solutions to satisfy Zaitsev boundary conditions at arbitrary transparency, and showing that this yields Nazarov-style circuit theory. The spatial profiles of Δ and M, the interface LDOS, the current-phase relation, and Ic(V) in the coexisting phase are all new. The paper also does a good job of checking itself against known limits: the low-transparency current reduces to Ambegaokar-Baratoff, and the transmission-averaged formulas reproduce Kulik-Omelyanchuk-type results in the appropriate limits. That gives me real confidence the formal machinery is not just internally consistent but physically faithful.\n\nThe soft spots are real but mostly addressable. The biggest one is Eq. (27): it is the workhorse trace formula from which the Josephson results follow, and the text says only that it comes from a “somewhat lengthy calculation.” Appendices A–D lay out the auxiliary-solution method, but they do not actually show the derivation of (27). A referee should insist on at least a sketch, because too much of the paper's output hangs on that single formula. The phase-gradient issue is related: the text says the ∂χ/∂x term “can be ignored for one does expect the phase to vary substantially across the junction,” which is surely a typo for “does not expect.” In a short junction with bulk reservoirs this is a standard and usually harmless approximation, but the paper should state it as an assumption and give a rough estimate of the neglected correction, especially for high transparency. That is a minor fix, not a fatal flaw. Also minor: Eq. (45) leaves cp as an order-one prefactor, so the compact current-phase formula is not fully explicit; that undercuts a bit of the practical appeal.\n\nOn the whole the central argument holds up. The method is sound, the checks are meaningful, and the physics—especially the Jc ~ Δ²/(eRNM) suppression when SDW dominates—is sensible. This is a paper for people working on iron-pnictide junctions, proximity effects, or quasiclassical circuit theory. It deserves a serious referee and would likely be a solid PRB-type paper after the derivation of (27) is expanded and the phase-gradient assumption is stated cleanly.\n\nMy recommendation: send it to peer review, with the request that the authors make the derivation of the central trace formula transparent and quantify or explicitly justify the neglect of ∂χ/∂x.","headline":"A useful and mostly sound extension of quasiclassical circuit theory to multiband SC+SDW junctions, but the central current formula deserves a fuller derivation before I would build on it.","tokens_in":19217,"tokens_out":2200,"would_cite":false,"duration_ms":24937,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.45.+c","74.50.+r","74.20.Rp"],"model":"deepseek-v4-flash","headline":"The paper claims that Josephson current in disordered multiband junctions with coexisting spin-density-wave order is fixed by bulk Green's functions and transmission eigenvalues alone, derived from the quasiclassical equations by Yip's…","keywords":["quasiclassical theory of superconductivity","Eilenberger equations","Josephson junctions","multiband superconductors","s±-wave pairing","spin-density-wave order","circuit theory of Andreev reflection","disorder in superconductors"],"falsifier":"Numerically solve the full Eilenberger equations retaining the phase-gradient term $\\partial\\chi/\\partial x$ for a junction whose electrodes carry a spatially varying phase and compare the resulting Josephson current with Eq. (27); a significant change would falsify the claim that the term can be ignored. A complementary experiment is to measure the current-phase relation of a disordered $s^{\\pm}$-wave junction with strong spin-density-wave order and check whether the critical current follows $\\Delta^2/(eR_N M)$ rather than the clean BCS scaling.","tokens_in":18252,"feed_emoji":"⚡","tokens_out":10212,"duration_ms":89144,"temperature":0.7,"pith_summary":"This paper claims that the standard microscopic route to superconducting junctions—the Yip auxiliary-solution method applied to Eilenberger equations with Zaitsev boundary conditions—automatically reproduces the connector rules of Nazarov circuit theory. It then works out a concrete case: a contact between two disordered two-band superconductors with coexisting $s^{\\pm}$-wave pairing and spin-density-wave order. The payoff is a compact current-phase relation in which the supercurrent is fixed by bulk Green's functions and the distribution of transmission eigenvalues alone, and is suppressed as $J_c \\sim \\Delta^2/(eR_N M)$ when the magnetic order dominates. If this is right, the current-phase relationship of such junctions is universal in the same sense as the standard Ambegaokar-Baratoff and Kulik-Omelyanchouk results, and the same machinery extends to other correlated-electron proximity circuits.","feed_headline":"One formula fixes supercurrent in multi-band junctions","feed_subtitle":"Bulk Green's functions plus transmission eigenvalues set the Josephson current, even when magnetic order coexists.","key_machinery":"The load-bearing device is Yip's auxiliary-solution construction: instead of solving the Eilenberger equation with the nonlinear Zaitsev boundary conditions directly, one solves two auxiliary (divergent) propagators on each side of the interface, combines them into the bounded physical solution via Eq. (A3), and finds that the interface matching conditions take the rational circuit-theory form of Eq. (A5). The essential objects are the bulk quasiclassical Green's functions ($g_z$, $f_z$, $s_z$, ...) and the transmission-eigenvalue distribution $\\rho(D)$, with the transparency $D$ entering only through the denominators. This machinery converts the full boundary-value problem into a one-parameter average over transmissions and yields the central trace formula Eq. (27) and its compact current-phase consequences, Eqs. (43) and (45).","core_discovery":"The central claim is that the auxiliary-solution construction converts a nonlinear boundary-value problem into circuit-theory form: the interface values of the quasiclassical propagators are rational functions of the bulk left and right propagators with the transparency $D$ as the only parameter, exactly as in Nazarov's boundary rules. In the coexistence model the Josephson current is therefore a transmission-eigenvalue average of bulk quasiclassical propagators, Eq. (43), which in the symmetric step-function limit reduces to $J(\\chi) = c_p \\frac{\\Delta^2}{eR_N} \\frac{\\sin\\chi}{\\sqrt{M^2+\\Delta^2\\cos^2(\\chi/2)}}$, Eq. (45). The paper further shows numerically that the superconducting order parameter varies only in a narrow neighbourhood of the interface while the spin-density-wave order varies on a longer scale, justifying the step-function approximation used to derive these universal current-phase relations. Disorder and magnetization enter through the anomalous Green's functions and smear the sharp features of the voltage-dependent critical current compared with a clean BCS contact.","pith_inferences":["A natural extension, not carried out in the paper, is to apply the same auxiliary-solution machinery to junctions with other coexisting orders (nematic, charge-density wave) and to three-band models where nematic order can appear.","The phrase in the paper justifying the dropped phase-gradient term ('one does expect the phase to vary substantially') is likely a typo for 'does not expect'; explicitly retaining $\\partial\\chi/\\partial x$ in a numerical solution would test whether the central trace formula survives when the phase varies inside the electrodes.","The predicted $J_c \\sim \\Delta^2/(eR_N M)$ suppression offers a diagnostic: in a disordered iron-pnictide junction, the critical current should dip inside the SDW-dominated part of the phase diagram, tracking the magnetic order parameter."],"forward_implications":["For short junctions between arbitrary superconductors, the Josephson current is fixed by bulk quasiclassical Green's functions and the transmission distribution, with no additional interface fitting parameters.","In the regime where spin-density-wave order dominates, $M \\gg \\Delta$, the critical current scales as $\\Delta^2/(eR_N M)$; this suppression is a concrete, measurable signature inside the magnetic phase.","The self-consistent profiles justify treating the order parameters as step functions, so current-phase relations computed from bulk values describe the full junction.","Disorder-induced interband scattering broadens the singular features in $I_c(V)$, so junction spectroscopy of such materials cannot be interpreted with clean BCS formulas.","The known universal Josephson results—Ambegaokar-Baratoff, Kulik-Omelyanchouk, Kupriyanov-Lukichev, and chaotic-cavity limits—all reappear as special cases of the same transmission averaging."],"supporting_citations":[{"why":"Supplies the auxiliary-solution method that converts the nonlinear boundary-value problem into circuit-theory form.","marker":"[31]"},{"why":"Gives the nonlinear Zaitsev boundary conditions at the interface that the auxiliary solutions must satisfy.","marker":"[32]"},{"why":"Original circuit theory of Andreev reflection whose connector rules are reproduced here.","marker":"[9]"},{"why":"Extended circuit-theory formulation used to identify the interface matching relations.","marker":"[10]"},{"why":"Gives the transmission-eigenvalue distribution for a diffusive connector (p=1) used in the averaging.","marker":"[12]"},{"why":"Gives the p=3/2 transmission distribution for randomly scattered interfaces used in the SINIS current-phase relation.","marker":"[13]"},{"why":"Establishes the connection between averaging over transmission eigenvalues and Usadel-equation SINIS results.","marker":"[54]"},{"why":"Disordered SNS Josephson current recovered from the p=1 average, one of the benchmarks the paper reproduces.","marker":"[52]"},{"why":"Earlier formulation of the two-band superconducting/spin-density-wave coexistence model and self-consistency equations used in the calculations.","marker":"[36]"}],"fun_headline_variants":["Circuit theory simplifies multiband supercurrent to one formula","Universal supercurrent from circuit rules in multiband junctions","Transmission eigenvalues set Josephson current in multiband contacts","One equation captures supercurrent in disordered multiband superconductors","Multiband Josephson current follows Nazarov-style circuit rules"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole calculation rests on assuming that the superconducting phase is essentially uniform inside each electrode so that the phase-gradient term generated by the unitary transformation can be dropped; if the phase actually varies across the junction, the central trace formula and every Josephson result built on it miss additional contributions.","fun_headline_variants_meta":{"raw":{"variants":["Circuit theory simplifies multiband supercurrent to one formula","Universal supercurrent from circuit rules in multiband junctions","Transmission eigenvalues set Josephson current in multiband contacts","One equation captures supercurrent in disordered multiband superconductors","Multiband Josephson current follows Nazarov-style circuit rules"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000634,"raw_usage":{"total_tokens":2932,"prompt_tokens":962,"completion_tokens":1970,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":578,"completion_tokens_details":{"reasoning_tokens":1885}},"tokens_in":578,"tokens_out":1970,"duration_ms":110028,"temperature":1.0,"reasoning_tokens":1885,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:54:01.465305+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically solve the full Eilenberger equations retaining the phase-gradient term $\\partial\\chi/\\partial x$ for a junction whose electrodes carry a spatially varying phase and compare the resulting Josephson current with Eq. (27); a significant change would falsify the claim that the term can be ignored. A complementary experiment is to measure the current-phase relation of a disordered $s^{\\pm}$-wave junction with strong spin-density-wave order and check whether the critical current follows $\\Delta^2/(eR_N M)$ rather than the clean BCS scaling.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the auxiliary-solution method that converts the nonlinear boundary-value problem into circuit-theory form."},{"cited_title":"Yip, Jour","cited_arxiv_id":null,"evidence_quote":"Gives the nonlinear Zaitsev boundary conditions at the interface that the auxiliary solutions must satisfy."},{"cited_title":"Argaman, Europhysics Letters (EPL) 38, 231 (1997)","cited_arxiv_id":null,"evidence_quote":"Gives the transmission-eigenvalue distribution for a diffusive connector (p=1) used in the averaging."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the connection between averaging over transmission eigenvalues and Usadel-equation SINIS results."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Disordered SNS Josephson current recovered from the p=1 average, one of the benchmarks the paper reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier formulation of the two-band superconducting/spin-density-wave coexistence model and self-consistency equations used in the calculations."}],"review_version":1}