Superconducting orbital diode effect in SN bilayers
Pith reviewed 2026-05-10 16:00 UTC · model grok-4.3
The pith
In thin SN bilayers, a nonideal interface with finite resistance can strengthen the orbital superconducting diode effect beyond the ideal case.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the diffusive regime with weak intralayer inhomogeneities, the superconducting diode effect in an SN bilayer depends nonmonotonically on the resistance of the SN interface; when the bilayer thickness is much smaller than the coherence length, this dependence allows a nonideal interface to produce a stronger diode effect than an ideal one.
What carries the argument
The orbital mechanism arising from the inhomogeneous distribution of Meissner currents due to spatially varying superfluid density, modulated by finite resistance at the SN interface.
If this is right
- The diode effect strength varies nonmonotonically with interface resistance for thin bilayers.
- Finite interface resistance can produce higher diode efficiency than a perfect interface.
- The asymmetry appears in both critical current and kinetic inductance.
- The effect holds only within the limit of weak intralayer inhomogeneities.
Where Pith is reading between the lines
- Practical SN devices with imperfect interfaces may show stronger diode behavior than ideal-interface calculations predict.
- Varying interface transparency in experiments on thin bilayers could map the predicted nonmonotonic curve.
- Interface engineering might be applied to other hybrid superconducting structures to optimize current-direction asymmetry.
Load-bearing premise
The analysis assumes only weak intralayer inhomogeneities in the diffusive regime.
What would settle it
Measure the ratio of critical currents for opposite directions in fabricated SN bilayers with tunable interface resistance and thickness much less than the coherence length to check for nonmonotonic enhancement at finite resistance.
Figures
read the original abstract
We study the superconducting diode effect (SDE) in a diffusive superconductor - normal metal (SN) bilayer subjected to an in-plane magnetic field. The supercurrent flows along the layers, perpendicular to the field. The SDE, manifested as an asymmetry in the critical (depairing) currents and kinetic inductance for opposite current directions, arises from an orbital mechanism due to the inhomogeneous distribution of the Meissner currents caused by a spatially varying superfluid density. Recently, Levichev et al. [Phys. Rev. B 108, 094517 (2023)] demonstrated the realization of this effect in such a structure, supporting numerical calculations for an ideal interface with an experiment. In this work, we investigate the influence of a nonideal interface with finite resistance on the SDE. Employing an analytical approach, we focus on limiting cases corresponding to weak intralayer inhomogeneities. We find that the strength of the SDE depends nonmonotonically on the interface resistance when the bilayer thickness is small compared to the coherence length. Remarkably, a nonideal interface can enhance the SDE compared to the ideal case.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analytically investigates the superconducting orbital diode effect (SDE) in diffusive SN bilayers under an in-plane magnetic field, focusing on the role of a nonideal interface with finite resistance. Building on prior numerical work for ideal interfaces, it restricts to limiting cases of weak intralayer inhomogeneities and derives that the SDE strength depends nonmonotonically on interface resistance when the bilayer thickness d is much smaller than the coherence length ξ; remarkably, nonideal interfaces can enhance the SDE relative to the ideal case.
Significance. If the analytical results hold within the stated approximations, the work provides controlled limiting-case expressions that clarify how interface resistance modulates the orbital SDE arising from inhomogeneous Meissner currents and spatially varying superfluid density. This extends the numerical findings of Levichev et al. and offers potential guidance for tuning SDE in thin hybrid structures, though the enhancement effect is tied to the weak-inhomogeneity diffusive regime.
major comments (1)
- The central nonmonotonic dependence and enhancement claim are derived only under the weak-inhomogeneity approximation (abstract and §3); the manuscript should explicitly state the range of validity (e.g., quantitative bound on the inhomogeneity parameter) and note that the result does not necessarily extrapolate beyond d ≪ ξ or to strong inhomogeneities, as this is load-bearing for the 'remarkably' enhancement statement.
minor comments (2)
- Notation for the interface resistance parameter (likely denoted R or similar in the boundary conditions) should be introduced earlier and used consistently with the Usadel equations referenced from prior work.
- Figure captions (if present) would benefit from explicit labels indicating which curves correspond to ideal vs. nonideal interfaces and the specific limiting-case parameters used.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive suggestion. We agree that the validity range of the weak-inhomogeneity approximation should be stated more explicitly and will revise the text accordingly.
read point-by-point responses
-
Referee: The central nonmonotonic dependence and enhancement claim are derived only under the weak-inhomogeneity approximation (abstract and §3); the manuscript should explicitly state the range of validity (e.g., quantitative bound on the inhomogeneity parameter) and note that the result does not necessarily extrapolate beyond d ≪ ξ or to strong inhomogeneities, as this is load-bearing for the 'remarkably' enhancement statement.
Authors: We agree with the referee that the nonmonotonic dependence on interface resistance and the possibility of enhancement are obtained within the weak-inhomogeneity limit. In the revised manuscript we will add an explicit statement in the abstract and at the beginning of §3 clarifying the range of validity: the analytical expressions hold when the inhomogeneity parameter (defined via the relative variation of the superfluid density across the bilayer) remains small enough for the perturbative expansion to be controlled, and we will note that the results are derived under the assumption d ≪ ξ. We will also caution that the enhancement effect is not guaranteed to persist outside this regime or for strong inhomogeneities, thereby qualifying the 'remarkably' phrasing. revision: yes
Circularity Check
No significant circularity; derivation self-contained in standard Usadel framework
full rationale
The paper derives the nonmonotonic SDE dependence analytically from the standard diffusive Usadel equations with interface boundary conditions, restricted to explicit limiting cases of weak intralayer inhomogeneities when d ≪ ξ. No parameters are fitted to data and then relabeled as predictions, no self-citation chain bears the central claim, and no ansatz or uniqueness theorem is smuggled in. The result follows directly from solving the governing equations under the stated approximations without reduction to inputs by construction.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Diffusive regime applies to the SN bilayer
- ad hoc to paper Weak intralayer inhomogeneities allow analytical treatment
Reference graph
Works this paper leans on
-
[1]
Despite the fact that the superfluid densityn(x) in this case is a steplike function, theτcorrections make the ratio between the densitiesn S0/nN0 sensitive to the current and magnetic field, which implies variation of the profile [see Eqs. (26) and (59)]. D. Conditions of applicability Now we discuss the applicability conditions of the above solutions. T...
-
[2]
Limit of thin layers To begin with, we consider the case of thin layers, dS(N) ≪ξ S(N) . We also assume that the respective parameters of the layers are of the same order:d S ∼ dN ∼d,σ S ∼σ N, andτ S ∼τ N ∼τ. In this case, the primary corrections responsible for the SDE are δnS,d ∼δn S,τ(dS/ξS)2 andδn N,τ ∼δn S,τ(τ∆) −1. The corresponding asymmetric corre...
-
[3]
Limit of moderately thick layers Now, we consider the case, in which at least one of the layers is moderately thick,d S(N) ∼ξ S (recall thatξ S ≪ ξN in the strongly resistive limit, which we consider now). In this regime, thedcorrection in the S layer dominates over the corrections in the N layer,δn S,d ≫δn N,τ , δnN,d. Consequently, the contribution to t...
-
[4]
=L k(x= 1, B= 0) [note thatI= 0 corresponds to x= 1 in Eq. (97)]. AtB= 0, the dependence is symmetric, Lk(I) =L k(−I), and thus the SDE is absent. In contrast, atB̸= 0 (the green curve corresponds toB= 0.6B c), this reciprocity is broken, resulting inL k(I)̸=L k(−I). Equation (97) together with Eq. (85) parametrically de- fines theL k(I, B) dependence (wi...
-
[5]
Derivation The kinetic inductance is defined through the relation between the kinetic energy per unit length and the total current. In the simplest linear case and in the absence of a magnetic field, this relation takes the form: Ek =L kI2/2.(A2) Linearity here implies that the density does not depend on the current or momentum, i.e.,n x = const q, leadin...
-
[6]
Relation to experiment We can rewrite Eq. (A3) for the kinetic energy as a work for acceleration of superconducting electrons with powerIEduring timet: Ek = Z t 0 IEdt,E=L k dI dt .(A9) Here,Eis the electric field (voltage per unit length) that accelerates electrons. The relation betweenEandIcor- responds to the kinetic contribution to the impedance Zk =−...
-
[7]
Effective self-consistency equation Consideration in Sec. III, which led to Eq. (40), should be supplemented by discussion of the self-consistency equation. Its effective form can be obtained from Eq. (9) with homogeneousθ S(x) =θand ∆(x) = ∆. In terms of the effective order parameterE g [see Eq. (41)], the effective self-consistency equation can be writt...
-
[8]
Critical pair-breaking parameter Equation (40) is valid for the pair-breaking parameter less than the critical value Γ c(T). At the critical value, a phase transition occurs and superconductivity disap- pears:θ(Γ→Γ c)→0 andE g(Γ→Γ c)→0. In this limit, Eq. (B1) produces the following equation for the critical pair-breaking parameter: ln Tc0 T =ψ 1 2 + Γc(T...
-
[9]
(B1) with an integral and calculate it explicitly
Zero-temperature limit At zero temperature, we can replace the sum in Eq. (B1) with an integral and calculate it explicitly. We ob- tain the dependence of the effective order parameter on the pair-breaking parameter,E g(Γ), determined by the following equation: ln(Eg/Eg0) +f(Γ/E g) = 0,(B5) with theffunction defined as f(z) = πz/4, z <1, arccoshz+ z...
-
[10]
(40) and solve the self-consistency Eq
Vicinity of the phase transition In the vicinity of the phase transition (θ→0,E g →0), we can find an explicit functionθ(ω n >0) from Eq. (40) and solve the self-consistency Eq. (B1): θ(ωn)≈ Eg Γ +ω n − E3 g ωn 2(Γ +ω n)4 ,(B9) E2 g(T,Γ) = eE2 g(T) 1−Γ/Γ c(T) .(B10) Here, the temperature-dependent prefactor is given by eE2 g(T) =TΓ c(T)Ξ Γc(T) 2πT ,(B11) ...
-
[11]
I. O. Kulik and I. K. Yanson,Josephson Effect In Super- conducting Tunneling Structures(John Wiley & Sons, New York, 1972)
work page 1972
-
[12]
A. Barone and G. Paterno,Physics and Applications of the Josephson Effect(Wiley, New York, 1982)
work page 1982
-
[13]
P. J. W. Moll and V. B. Geshkenbein, Evolution of su- perconducting diodes, Nat. Phys.19, 1379 (2023)
work page 2023
-
[14]
L. S. Levitov, Yu. V. Nazarov, and G. M. Eliashberg, Magnetostatics of superconductors without an inversion center, Pis’ma Zh. Eksp. Teor. Fiz.41, 365 (1985), [JETP Lett.41, 445 (1985)]
work page 1985
-
[15]
V. M. Edelstein, The Ginzburg–Landau equation for su- perconductors of polar symmetry, J. Phys.: Condens. Matter8, 339 (1996)
work page 1996
-
[16]
V. M. Krasnov, V. A. Oboznov, and N. F. Pedersen, Fluxon dynamics in long Josephson junctions in the pres- ence of a temperature gradient or spatial nonuniformity, Phys. Rev. B55, 14486 (1997)
work page 1997
-
[17]
J. B. Majer, J. Peguiron, M. Grifoni, M. Tusveld, and J. E. Mooij, Quantum ratchet effect for vortices, Phys. Rev. Lett.90, 056802 (2003)
work page 2003
-
[18]
J. E. Villegas, S. Savel’ev, F. Nori, E. M. Gonzalez, J. V. Anguita, R. Garc´ ıa, and J. L. Vicent, A superconducting reversible rectifier that controls the motion of magnetic flux quanta, Science302, 1188 (2003)
work page 2003
-
[19]
D. Y. Vodolazov, B. A. Gribkov, S. A. Gusev, A. Yu. Klimov, Yu. N. Nozdrin, V. V. Rogov, and S. N. Vdovichev, Considerable enhancement of the critical cur- rent in a superconducting film by a magnetized magnetic strip, Phys. Rev. B72, 064509 (2005)
work page 2005
-
[20]
C. C. de Souza Silva, J. Van de Vondel, M. Morelle, and V. V. Moshchalkov, Controlled multiple reversals of a ratchet effect, Nature440, 651 (2006)
work page 2006
-
[21]
M. Morelle and V. V. Moshchalkov, Enhanced criti- cal currents through field compensation with magnetic strips, Appl. Phys. Lett.88, 172507 (2006)
work page 2006
-
[22]
A. Yu. Aladyshkin, D. Yu. Vodolazov, J. Fritzsche, R. B. G. Kramer, and V. V. Moshchalkov, Reverse- domain superconductivity in superconductor- ferromagnet hybrids: Effect of a vortex-free channel on the symmetry ofI-Vcharacteristics, Appl. Phys. Lett. 97, 052501 (2010)
work page 2010
-
[23]
M. A. Silaev, A. Yu. Aladyshkin, M. V. Silaeva, and A. S. Aladyshkina, The diode effect induced by domain-wall superconductivity, J. Phys.: Condens. Matter26, 095702 (2014)
work page 2014
-
[24]
T. Yokoyama, M. Eto, and Yu. V. Nazarov, Anomalous Josephson effect induced by spin-orbit interaction and Zeeman effect in semiconductor nanowires, Phys. Rev. B 89, 195407 (2014)
work page 2014
-
[25]
R. Wakatsuki, Y. Saito, S. Hoshino, Y. M. Itahashi, T. Ideue, M. Ezawa, Y. Iwasa, and N. Nagaosa, Nonre- ciprocal charge transport in noncentrosymmetric super- conductors, Sci. Adv.3, e1602390 (2017)
work page 2017
-
[26]
C.-Z. Chen, J. J. He, M. N. Ali, G.-H. Lee, K. C. Fong, and K. T. Law, Asymmetric Josephson effect in inversion symmetry breaking topological materials, Phys. Rev. B 98, 075430 (2018)
work page 2018
- [27]
-
[28]
N. Nagaosa and Y. Yanase, Nonreciprocal transport and optical phenomena in quantum materials, Annu. Rev. 20 Condens. Matter Phys.15, 63 (2024)
work page 2024
-
[29]
D. Shaffer and A. Levchenko, Theories of superconduct- ing diode effects, arXiv:2510.25864
-
[30]
M. Yu. Levichev, I. Yu. Pashenkin, N. S. Gusev, and D. Yu. Vodolazov, Finite momentum superconductivity in superconducting hybrids: Orbital mechanism, Phys. Rev. B108, 094517 (2023)
work page 2023
-
[31]
C. Baumgartner, L. Fuchs, A. Costa, S. Reinhardt, S. Gronin, G. C. Gardner, T. Lindemann, M. J. Manfra, P. E. Faria Junior, D. Kochan, J. Fabian, N. Paradiso, and C. Strunk, Supercurrent rectification and magne- tochiral effects in symmetric Josephson junctions, Nat. Nanotechnol.17, 39 (2022)
work page 2022
-
[32]
K. D. Usadel, Generalized diffusion equation for super- conducting alloys, Phys. Rev. Lett.25, 507 (1970)
work page 1970
- [33]
-
[34]
L. N. Cooper, Superconductivity in the neighborhood of metallic contacts, Phys. Rev. Lett.6, 689 (1961)
work page 1961
-
[35]
P. G. de Gennes, Boundary effects in superconductors, Rev. Mod. Phys.36, 225 (1964)
work page 1964
-
[36]
W. L. McMillan, Tunneling model of the superconducting proximity effect, Phys. Rev.175, 537 (1968)
work page 1968
-
[37]
A. A. Golubov and M. Yu. Kupriyanov, Josephson effect in SNINS and SNIS tunnel structures with finite trans- parency of the SN boundaries, Zh. Eksp. Teor. Fiz.96, 1420 (1989), [Sov. Phys. JETP69, 805 (1989)]
work page 1989
-
[38]
A. A. Golubov, Proximity effect in dirty N/S multilayers, inSuperconducting Superlattices and Multilayers, Vol. 2157, edited by I. Bozovic (SPIE, 1994) pp. 353–362
work page 1994
-
[39]
Ya. V. Fominov and M. V. Feigel’man, Superconductive properties of thin dirty superconductor–normal-metal bi- layers, Phys. Rev. B63, 094518 (2001)
work page 2001
-
[40]
A. I. Larkin and Yu. N. Ovchinnikov, Vortex motion in superconductors, inNonequilibrium Superconductivity, edited by D. N. Langenberg and A. I. Larkin (Elsevier, New York, 1986) p. 530
work page 1986
-
[41]
A. D. Zaikin and G. F. Zharkov, Theory of wide dirty SNSjunctions, Fiz. Nizk. Temp.7, 375 (1981), [Sov. J. Low Temp. Phys.7, 184 (1981)]
work page 1981
-
[42]
M. Yu. Kuprianov and V. F. Lukichev, Influence of boundary transparency on the critical current of “dirty” SS′S structures, Zh. Eksp. Teor. Fiz.94, 139 (1988), [Sov. Phys. JETP67, 1163 (1988)]
work page 1988
-
[43]
To avoid confusion, we note that in Ref. [20], the quan- tityq [20] 0 , chosen as an independent variable characteriz- ing the distribution of the condensate momenta inside the SN bilayer, is the thickness-averaged momentum: q[20] 0 = R qs(x)dx/d. AtB̸= 0 andI= 0, the current distribution is such thatq [20] 0 ̸= 0, therefore this state was referred to as ...
-
[44]
Theτ S(N) parameters defined in Eq. (25), can be es- timated asτ S(N) ∼ T −1dS(N) /vS(N) , whereTis the transparency of the interface andv S(N) are the Fermi velocities in the layers. Note thatτ S(N) do not depend on the mean free pathsl S(N) in the layers and contain the ballistic timesd S(N) /vS(N) of passage across the lay- ers. If the transverse resis...
-
[45]
A. A. Abrikosov and L. P. Gor’kov, Contribution to the theory of superconducting alloys with paramagnetic im- purities, Zh. Eksp. Teor. Fiz.39, 1781 (1960), [Sov. Phys. JETP12, 1243 (1961)]
work page 1960
-
[46]
Using the relationσ= 2e 2νDbetween the conductivity and the diffusion constant, we can rewrite Eq. (43) in terms of the total (volume-integrated) density of states νS(N) VS(N) or in terms of the level spacingsδ S(N) : ⟨F⟩= FSνSVS +F N νN VN νSVS +ν N VN = FSδ−1 S +F N δ−1 N δ−1 S +δ −1 N
-
[47]
J. H. Claassen, J. E. Evetts, R. E. Somekh, and Z. H. Barber, Observation of the superconducting proximity effect from kinetic-inductance measurements, Phys. Rev. B44, 9605 (1991)
work page 1991
-
[48]
A. A. Neilo, S. V. Bakurskiy, N. V. Klenov, I. I. Soloviev, and M. Yu. Kupriyanov, Magnetic control of the kinetic inductance in elements of superconducting electronics, JETP Lett.121, 58 (2025)
work page 2025
-
[49]
Maki, Gapless superconductivity, inSuperconductiv- ity, Vol
K. Maki, Gapless superconductivity, inSuperconductiv- ity, Vol. 2, edited by R. D. Parks (Marcel Dekker, New York, 1969) Chap. 18
work page 1969
-
[50]
Ya. V. Fominov, N. M. Chtchelkatchev, and A. A. Gol- ubov, Nonmonotonic critical temperature in supercon- ductor/ferromagnet bilayers, Phys. Rev. B66, 014507 (2002)
work page 2002
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.