Abstract
The suggestion that non-reciprocal critical current (NRC) may be an intrinsic property of non-centrosymmetric superconductors has generated renewed theoretical and experimental interest motivated by an analogy with the non-reciprocal resistivity due to the magnetochiral effect in uniform materials with broken spatial and time-reversal symmetry. Theoretically it has been understood that terms linear in the Cooper pair momentum do not contribute to NRC, although the role of higher-order terms remains unclear. In this work we show that critical current non-reciprocity is a generic property of multilayered superconductor structures in the presence of magnetic field-generated diamagnetic currents. In the regime of an intermediate coupling between the layers, the Josephson vortices are predicted to form at high fields and currents. Experimentally, we report the observation of NRC in nanowires fabricated from InAs/Al heterostructures. The effect is independent of the crystallographic orientation of the wire, ruling out an intrinsic origin of NRC. Non-monotonic NRC evolution with magnetic field is consistent with the generation of diamagnetic currents and formation of the Josephson vortices. This extrinsic NRC mechanism can be used to design novel devices for superconducting circuits.
Similar content being viewed by others
Introduction
Diodes are the most basic elements of semiconductor electronics and development of superconducting diodes can extend the functionality of the superconducting circuitry. A non-reciprocal critical current (NRC) in a multiply-connected superconductors is a well known effect and can be readily observed in, e.g., an asymmetric superconducting rings1. An implicit suggestion that NRC may be an intrinsic property of non-centrosymmetric superconductors2 generated a renewed theoretical and experimental interest motivated by an analogy with the non-reciprocal resistivity due to the magnetochiral effect, which can appear in uniform materials with broken spatial and time-reversal symmetry3. However, a direct analogy between corrections to resistivity and superconducting current is misleading because the anisotropy of scattering is caused by the spin-orbit effects, while the proposed origin of nonreciprocity in singlet-pairing superconductors is a spin-independent Lifshitz invariant4,5. It has been demonstrated in the literature6 that in uniform singlet superconductors in constant magnetic field the Lifshitz invariants can be eliminated by a gauge (Galilean) transformation from both the Ginzburg-Landau (GL) equation and the expression for the supercurrent, so that linear in Cooper pair momentum terms do not lead to nonreciprocity. Phenomenological treatment shows that cubic in the Cooper pair momentum terms can lead to NRC corrections7,8. It has been also suggested that the Rashba terms in the electron spectrum contribute to nonreciprocity within the formalism of quasiclassical Eilenberger equations9. A symmetry analysis and microscopic calculations of the Cooperon propagator in the presence of the Zeemann effect and a linear or cubic Dresselhaus spin-orbit interactions in the electron spectrum show that the NRC magnitude and sign depend on the crystallographic orientation of the supercurrent flow (Y.L.G., J.I.V., A.S. & L.P.R., manuscript in preparation). Such anisotropy should characterize non-reciprocity in superconductor/semiconductor heterostructures and cubic uniform singlet noncentrosymmetric superconductors. Another suggested mechanism of NRC is the formation of non-uniform currents in superconducting multilayers10,11.
In this paper we show that NRC naturally arises in the presence of the magnetic-field–generated diamagnetic currents when neighboring layers in multilayer superconductors are strongly coupled. The total current through a multilayer structure is divided between the layers as an inverse ratio of their kinetic inductances. Initially only one of the layers reaches the maximum (critical) current as the total current increases. A further current increase in the superconducting state requires generating a phase difference between the layers, which adds Josephson energy penalty to the total energy, and in the case of strong interlayer coupling (large Josephson energy) the whole system transitions into a normal state. Field-generated diamagnetic currents will either increase or decrease an external current where transition to the normal state occurs, thus leading to the NRC. In the regime of intermediate interlayer coupling strengths, an interlayer phase difference can change by 2π, leading to the formation of Josephson vortices. Experimentally, we report observation of NRC in nanowires fabricated from InAs/Al heterostructures. The observed non-monotonic evolution of NRC with magnetic field is consistent with the formation of Josephson vortices. Our findings show that the extrinsic contribution to NRC is generic to multilayer superconductors, and may provide a compelling explanation to the NRC observed in Refs. 2 and12, in the latter work magnetic flux produced by a Co layer generates opposite diamagnetic currents in the adjacent Nb and V layers.
The term “superconducting diode effect” has been used to describe NRC in different systems, including thin superconducting films13,14,15,16,17,18 and Josephson junctions19,20,21,22,23,24,25,26,27. In the former experiments the presence of out-of-plane magnetic field and formation of vortices is essential for the observation of NRC, in this case the critical current is determined by the strength and symmetry of the flux pinning potential. In the latter case the critical current in Josephson junctions is determined by the overlap of Andreev states. In this paper we restrict our discussion to the origin of NRC in long nanowires, where critical current is determined by the de-pairing velocity of Cooper pairs (the Bardeen limit28).
Results
We have studied switching currents Isw defining a transition from superconducting to normal state in nanowires fabricated from Al/InGaAs/InAs/InGaAs heterostructures29, where patterned Al top layer forms a nanowire and induces superconductivity in a high mobility InAs quantum well via the proximity effect. An AFM micrograph of a typical device is shown in the inset in Fig. 1. A typical current-voltage characteristic exhibits a sharp switching transition limited by the current resolution ( < 5 nA for the fastest sweep rates used in our experiments). A histogram of switching currents \({I}_{sw}^{\pm }\) for positive (+) and negative (-) current sweeps is shown in Fig. 1a for 10,000 sweeps. Field dependence of average values \(\langle {I}_{sw}^{+}\rangle\) and \(\langle {I}_{sw}^{-}\rangle\) is plotted in Fig. 1b for the in-plane field B⊥ perpendicular to the wire. The \(\langle {I}_{sw}^{+}\rangle\) and \(\langle {I}_{sw}^{-}\rangle\) can be separated into a symmetric \(\langle {I}_{sw}\rangle=(\langle {I}_{sw}^{+}\rangle+\langle {I}_{sw}^{-}\rangle )/2\) and asymmetric \(\Delta I=\langle {I}_{sw}^{+}\rangle -\langle {I}_{sw}^{-}\rangle\) parts, the latter being the non-reciprocal component of the supercurrent. Both 〈Isw〉 and ΔI are nonmonotonic functions of magnetic field. As shown in the Supplement, a minima of 〈Isw〉 at low fields vanishes above 350 mK (0.3 TC), while there is no change in ΔI at least up to 750 mK (>0.6TC). This difference in energy scales for the appearance of NRC and non-monotonic evolution of 〈Isw〉 indicates that these are unrelated phenomena, and below we focus on the origin of NRC. Some devices were fabricated with a top gate, which allows electrostatic control of the electron density in the InAs layer not covered by Al; we found that depletion of the 2D electron gas in the exposed InAs results in a slight increase of 〈Isw〉 but does not affect ΔI. Similar field effect has been observed previously in superconductor nanodevices30 and was attributed to the presence of quasiparticles31, a conclusion consistent with the observed gate dependence of the 〈Isw〉.
a Histograms of switching currents for 10,000 positive \({I}_{sw}^{+}\) and negative \({I}_{sw}^{-}\) current sweeps performed at T = 30 mK and B⊥ = 100 mT. Inset shows a typical current-voltage characteristic. b Average switching current for positive \(\langle {I}_{sw}^{+}\rangle\) and negative \(\langle {I}_{sw}^{-}\rangle\) sweeps, non-reciprocal difference \(\Delta I=\langle {I}_{sw}^{+}\rangle -\langle {I}_{sw}^{-}\rangle\) and an average of all sweeps 〈Isw〉 is plotted as a function of in-plane magnetic field B⊥. In (c) enlarged ΔI data is colored to signify non-monotonic field dependence and multiple sign changes. d Dependence of ΔI on in-plane field orientation is measured at a constant B = 100 mT. Blue line is a fit with a sine function. Insert shows an AFM image of a 3μm-long wire connected to wide contacts, yellow areas are Al, in darker areas Al is removed and InAs is exposed.
Unlike the linear in Cooper pair momentum terms, higher order terms cannot be removed by gauge transformation and it was shown that the presence of terms ~ α3Q3Δ2 cubic in the Cooper pair momentum in an expansion of the Ginsburg-Landau coefficients can generally lead to non-zero ΔI which is a non-monotonic function of B and can even change sign7,9 (here Q = − iℏ ∇ − 2eA is a generalized Cooper pair momentum, A is electromagnetic vector-potential). However, for proximitized InAs layer, a generation of the terms higher order in the Cooper pair momentum in the presence of the Rashba spin-orbit and Zeeman interactions coexists with a similar generation of such terms due to the Dresselhaus spin-orbit interactions. The importance of the Dresselhaus-like terms in the electron spectrum is not limited to proximity structures, and they can play significant role in any noncentrosymmetric material. Investigation of realistic cubic terms in the Cooper pair momentum showed (Y.L.G., J.I.V., A.S. & L.P.R., manuscript in preparation) that nonreciprocity becomes highly anisotropic as a result of Dresselhaus-induced contribution. For comparison with experiments, it is instructive to express the odd in Cooper pair momentum part of the kinetic energy in coordinates rotated by π/2 with respect to the principal crystallographic axes of InAs, where \(\hat{x}||[1\overline{1}0]\) and \(\hat{y}||[110]\), see insert in Fig. 2a. In these coordinates, the cubic in the Cooper pair momentum kinetic term originating from the cubic Dresselhaus electron spin-orbit interaction reads
where coefficient κ contains the Dresselhaus constant βD and other material parameters. The resulting NRC correction to the supercurrent is
This correction is independent of the sign of I and is added or subtracted to the B = 0 current value depending on the direction of the current flow. Here Bx and By enter symmetrically for wires oriented along x and y. However, in the configuration with the current \(I\parallel \hat{x}\) and magnetic field By and the configuration with \(I\parallel \hat{y}\) and Bx, this expression has opposite signs for the same mutual orientation of I and B, see inset in Fig. 2a. Thus, the Dresselhaus-induced contribution results in NRC with opposite sign for wires oriented along \([1\overline{1}0]\) and [110] crystallographic axis. The cubic (and generally all odd) in Cooper pair momentum terms originating from the Rashba electronic interactions, when added with the Dresselhaus-induced terms, will produce anisotropy in the absolute value of NRC, and, in particular, different values of non-reciprocal asymmetrical component of the current for those two directions. Theoretical investigation of electronic spectra of these systems32 suggests that in narrow InAs quantum wells cubic Dresselhaus terms are larger than the Rashba terms. The lower limit for the value of the Dresselhaus contribution can be extracted from the total spin-orbit anisotropy (which is defined by the ratio between a linear Rashba, and a linear and cubic Dresselhaus terms in electronic spectrum), which was measured to be 70% in spin-galvanic and circular photogalvanic experiments33 and > 15% in transport experiments34,35. Such anisotropies must result in the corresponding crystallographic anisotropy of the NRC, which is not observed in our experiments, Fig. 2a. Therefore, we conclude that the NRC we observed is not intrinsic. The observed NRC does not depend on the wire length, Fig. 2b, which rules out trivial effects related to the formation of spurious loops due to the presence of wire/contact boundaries.
a NRC is plotted for two 2 μm-long wires oriented along [110] and \([1\overline{1}0]\) crystallographic axes. Insets define mutual orientataion of wires and fields. b NRC for 2, 3, and 5 μm wires. The top and bottom curves are shifted vertically by 0.2 μA. Brackets with arrows indicate a maximum ΔB needed to insert a flux ϕ0 = h/2e in the area defined by the corresponding wire lengths, as indicated by a dashed loop in the inset. An effective length for the period marked by a magenta bracket is l = 0.5μm for the same loop.
While recent interest in NRC has been motivated by a possibility of the intrinsic origin of the effect, NRC naturally arises in multiply-connected superconductors. In superconducting loops, the critical current is modulated by an external flux ϕ = BSloop piercing the loop. In a loop with asymmetric arms, the current maximum is shifted from B = 0, and the sign of the shift depends on the direction of the current as shown in Fig. 3a. A nonreciprocal component of the switching current ΔI is linear in B in the vicinity of B = 0, reaches extrema at ϕ ≈ ϕ0/4, changes sign and oscillates with a period Δϕ = ϕ0. Thus, an asymmetric loop is the simplest “superconducting diode”. There is a clear similarity between ΔI measured in an asymmetric superconducting loop and in an Al/InAs nanowire as emphasized in Fig. 3b, suggesting that nonmonotonic NRC in our nanowires may be due to emerging current loops.
a An average switching current for positive \(\langle {I}_{sw}^{+}\rangle\) and negative \(\langle {I}_{sw}^{-}\rangle\) sweeps, non-reciprocal difference \(\Delta I=\langle {I}_{sw}^{+}\rangle+\langle {I}_{sw}^{-}\rangle\) and an average of all sweeps 〈Isw〉 plotted as a function of out-of-plane magnetic field Bo for a loop shown in the insert in (b). Note that 〈Isw〉 is maximal while ΔI = 0 when the flux ϕ = nϕ0. In (b) ΔI for the nanowire and the loop are plotted together as a function of a reduced flux ϕ/ϕ0, where we used Swire = 0.0052μm2 for the effective area in the wire and Sloop = 2.59μm2 in the loop.
External magnetic field generates circular diamagnetic currents in a superconductor, as shown schematically in Fig. 4a,b, and these currents affect the critical current. In homogeneous superconductors the presence of diamagnetic currents will not result in the critical current non-reciprocity, but in a heterogeneous superconductor, in general, their presence will lead to NRC. Qualitatively, the origin of NRC can be understood from a simplified model of a superconductor heterostructure represented as two coupled zero-thickness superconducting wires separated by a distance d, Fig. 4c. The total energy of the two-wire system can be written as a sum of kinetic and Josephson energies,
where \({{{{{{{{\mathcal{E}}}}}}}}}_{k}={L}_{1}{I}_{1}^{2}+{L}_{2}{I}_{2}^{2},\,{{{{{{{{\mathcal{E}}}}}}}}}_{J}\) is the Josephson coupling, Δϕ = ϕ1(x) − ϕ2(x) is the phase difference between superconducting condensates, and Li are the kinetic inductances per unit length in wires labeled by an index i = 1, 2. Supercurrents in each wire \({I}_{i}={(2e{L}_{i})}^{-1}(\hslash {\partial }_{x}{\phi }_{i}-2e{A}_{x})\) should satisfy charge conservation constraint I1(x) + I2(x) = Iext, where Iext is the applied external current. Detailed solution for this model can be found in the Supplementary Materials, and we outline now the main results. For small external currents (I1 < Ic1 and I2 < Ic2, where Ici are the critical currents in the wires) it is energetically favorable to keep the phase difference Δϕ constant (Δϕ = 0 for \({{{{{{{{\mathcal{E}}}}}}}}}_{J}\, > \,0\)). Then, the currents can be expressed as I1 = I01 + Idia and I2 = I02 − Idia, where I01, I02 ∝ Iext with I01/I02 = L2/L1 = η−1 and Idia = Byd/(L1 + L2). Dependence of I1 and I2 on Iext for By > 0, By = 0 and By < 0 is plotted schematically in Fig. 4d. As Iext increases and one of the currents (I1 in our example) reaches the critical value Ic1, further external current increase requires an increase of ∣Δϕ∣ because the excess current has to flow through the remaining superconducting wire with the current I2. In the case of weak interwire coupling, \({{{{{{{{\mathcal{E}}}}}}}}}_{J}\ll {{{{{{{{\mathcal{E}}}}}}}}}_{k}\), deviation of Δϕ from zero does not lead to a significant energy penalty and the critical current of the whole system \({I}_{c}^{{\prime} }={I}_{c1}+{I}_{c2}\) does not depend on the magnetic field direction. In the opposite regime of strong coupling, \({{{{{{{{\mathcal{E}}}}}}}}}_{J}\gg {{{{{{{{\mathcal{E}}}}}}}}}_{k}\), the energy cost associated with the formation of Josephson currents (the last term in Eq. (3)) is prohibitively high and the whole system transitions to a normal state at Iext ≈ (1 + η)(Ic1 − Idia), resulting in ΔI = − 2(η + 1)Idia(By) (this equation is correct for β > η + (η + 1)Idia/Ic1, where β = Ic2/Ic1, NRC for other scenarios is listed in the Supplementary Materials). Thus a superconducting diode effect is a generic property of coupled multilayer superconductors.
Diamagnetic currents in (a) a uniform superconductor and (b) a heterostructure. c A two-layer heterostructure is modeled as two zero thickness wires separated by a distance d with coupling described by the Josephson energy \({{{{{{{{\mathcal{E}}}}}}}}}_{J}\). d Schematic of current distribution between the wires Ii = I0i − (−1)iIdia, i = 1, 2, and the phase difference Δϕ as a function of an external current Iext = I1 + I2 for B = 0 (solid lines), B > 0 (dotted lines) and B < 0 (dashed lines). For weakly coupled wires \({{{{{{{{\mathcal{E}}}}}}}}}_{J}\ll {{{{{{{{\mathcal{E}}}}}}}}}_{k}\), the critical current is field-independent \({I}_{c}^{{\prime} }={I}_{c1}+{I}_{c2}\), see the text; the critical current is reduced and acquires a linear-in-B correction in a strong coupling regime \({{{{{{{{\mathcal{E}}}}}}}}}_{J}\gg {{{{{{{{\mathcal{E}}}}}}}}}_{k}\) due to the phase locking Δϕ = 0. e In the intermediate coupling regime \({{{{{{{{\mathcal{E}}}}}}}}}_{J} \sim {{{{{{{{\mathcal{E}}}}}}}}}_{k}\) Josephson vortices may form generating a 2π phase twist, in this case NRC becomes a non-monotonic function of B. f Calculated NRC ΔI is plotted as a function of flux Φ = SvBy for several Bx, Eq. (S12).
As By and diamagnetic currents increase, the phase locking condition Δϕ = 0 along the length of the wires leads to a significant increase of \({{{{{{{{\mathcal{E}}}}}}}}}_{k}\). At a critical field Bc = (3/π2)Φ0/(lJd) it becomes energetically favorable to reduce the overall energy by twisting the phase difference by 2π locally forming a Josephson vortex (\({l}_{J}\,\approx \,{\Phi }_{0}/(2\pi \sqrt{2{{{{{{{{\mathcal{E}}}}}}}}}_{J}{L}_{2}})\) and Φ0 = h/2e is the flux quantum). Evolution of the phase difference \(\Delta \phi (x)=4\arctan [\exp (x/{l}_{J})]\), vortex-induced currents in the wires Iv(x), and interwire Josephson current density jJ(x) across a vortex are shown schematically in Fig. 4e. The maximum of Iv(x) at the center of the vortex determines the Josephson vortex contribution to NRC. In the absence of quantum fluctuations formation of a vortex is accompanied by an abrupt re-distribution of currents between the wires, which results in a sawtooth NRC dependence on the magnetic field. Generation of multiple Josephson vortices does not modify NRC compared to a single vortex case unless the vortices significantly overlap so that the maximum of Iv(x) exceeds its single-vortex value.
In Fig. 4f we plot ΔI(By) for several Bx using Eq. (S12) in the Supplemental Material. A gradual change of ΔI near Φ0/2 is due to quantum fluctuations of the winding number due to strong coupling of the vortex to current-carrying wires. This smearing is similar to the gradual change of a critical current in a ring connected to superconducting leads (Fig. 3), as compared to an abrupt reversal of persistent currents at Φ0/2 in isolated rings36. The period of oscillations of ΔI corresponds to the flux threading an effective vortex area Sv = (π2/3)lJd = lvd. The period ΔB⊥ = 400 mT translates into the length lv ≈ 500nm, where Δϕ substantially deviates from zero. We estimate \({l}_{v}\, < \,{\xi }_{InAs}=\sqrt{{\xi }_{InAs}^{0}{l}_{InAs}^{m}}\approx 750\) nm and expect proximity-induced superconductivity in InAs to be preserved in the presence of a vortex. Here we use \({\xi }_{InAs}^{0}=\hslash {v}_{F}/\pi {\Delta }^{*}\,\approx \,1.8\,\mu\)m, induced gap in InAs Δ* ≈ Δ = 1.796kBTc = 230μeV (induced gap is close to the gap of Al in these heterostructures37), and the mean free path in uncapped InAs 2D gas \({l}_{InAs}^{m}\,\approx \,300\) nm.
Finally, we use the two-wire model to estimate the temperature and in-plane field \({B}_{\parallel }\parallel \hat{x}\) dependences of NRC assuming that both parameters affect the Cooper pair density n2 in InAs. In the vicinity of B⊥ = 0 the amplitude of \(\Delta I\propto {L}_{2}^{-1}\propto {n}_{2}\) and is expected to decrease with an increase of T or B∥. The critical field \({B}_{c}\propto \sqrt{{{{{{{{{\mathcal{E}}}}}}}}}_{J}{L}_{2}}\) depends on \({{{{{{{{\mathcal{E}}}}}}}}}_{J}\propto {n}_{2}\), and the period of oscillations is expected to be T- and B∥-independent, Fig. 4d. Josephson coupling \({{{{{{{{\mathcal{E}}}}}}}}}_{J}\) is exponentially sensitive to the thickness of the InGaAs spacer and we expect slight variations of the period ΔB⊥ between the samples. These qualitative estimates are consistent with experimental observations, see Figs. S2 and S3 in the Supplemental Material.
Methods
Materials
The wafer was grown using Molecular Beam Epitaxy (MBE) on an InP substrate. The heterostructure consists of 1 μm graded InxAl1−xAs insulating buffer followed by a In0.75Ga0.25As(4nm)/InAs(7nm)/In0.75Ga0.25As(10nm) multilayer structure capped in-situ with 7nm of Al. The two-dimensional electron gas has a peak mobility of 28000 cm2/Vs at a density 8 × 1011 cm−2.
Sample fabrication
The nanowires were fabricated using standard electron beam lithography. The mesas were defined by first removing the top Al layer with Al etchant Transene D and then a deep wet etching using H3PO4:H2O2:H2O:C6H8O7 (1ml:8ml:85ml:2g). Nanowires are defined in the second step of lithography by patterning the Al layer. Some devices have a top electrostatic gate, in these devices a 20 nm HfO2 is grown by atomic layer deposition followed by a deposition of a Ti/Au (10/100 nm) gate.
Measurements
Current-voltage sweeps were performed using a homemade high-speed high resolution DAC/ADC (digital-to-analog and analog-to-digital converter) setup. The sweeps were automatically interrupted at the superconductor-normal transition (Isw) in order to minimize device heating. Current sweep rate and delay between sweeps have been optimized to obtain <5 nA current resolution and to keep device temperature <50 mK at the base temperature of the fridge. The data has been corrected for an instrumental cooldown-dependent constant current offset (generated in the circuit by uncompensated voltages in the system and limited by a 100 kΩ current source resistor) to insure that ΔI = 0 at B = 0.
Data availability
All the data that support the findings of this study are available from the corresponding author on reasonable request. Source data are provided with this paper.
Code availability
The code used for the numerical simulation has been provided with this paper as Supplementary Data 1.
References
Burlakov, A. A., Gurtovoi, V. L., Il’in, A. I., Nikulov, A. V. & Tulin, V. A. Superconducting quantum interference device without Josephson junctions. JETP Lett. 99, 169 (2014).
Ando, F. et al. Observation of superconducting diode effect. Nature 584, 373 (2020).
Rikken, G. L. J. A. & Wyder, P. Magnetoelectric Anisotropy in Diffusive Transport. Phys. Rev. Lett. 94, 016601 (2005).
Levitov, L. S., Nazarov, Y. V. & Eliashberg, G. M. Magnetostatics of Superconductors without an inversion center. JETP Lett. 41, 445 (1985).
Edelstein, V. M. The Ginzburg - Landau equation for superconductors of polar symmetry. J. Phys. Condens. Matter 8, 339 (1996).
Agterberg, D. F. Magnetoelectric effects, helical phases, and FFLO phases. In Non-centrosymmetric Superconductors (eds. Bauer, E. & Sigrist, M.) 155–170 (Springer, 2012).
Daido, A., Ikeda, Y. & Yanase, Y. Intrinsic Superconducting Diode Effect. Phys. Rev. Lett. 128, 037001 (2022).
He, J. J., Tanaka, Y. & Nagaosa, N. A phenomenological theory of superconductor diodes. N. J. Phys. 24, 053014 (2022).
Ilić, S. & Bergeret, F. S. Theory of the supercurrent diode effect in Rashba superconductors with arbitrary disorder. Phys. Rev. Lett. 128, 177001 (2022).
Vodolazov, D. Y. et al. Peculiar superconducting properties of a thin film superconductor-normal metal bilayer with large ratio of resistivities. Supercond. Sci. Technol. 31, 115004 (2018).
Zinkl, B., Hamamoto, K. & Sigrist, M. Symmetry conditions for the superconducting diode effect in chiral superconductors. Phys. Rev. Res. 4, 033167 (2022).
Narita, H. et al. Field-free superconducting diode effect in noncentrosymmetric superconductor/ferromagnet multilayers. Nat. Nanotechnol. 17, 823 (2022).
Bauriedl, L. et al. Supercurrent diode effect and magnetochiral anisotropy in few-layer NbSe2. Nat. Commun. 13, 4266 (2022).
Lin, J.-x. et al. Zero-field superconducting diode effect in small-twist-angle trilayer graphene. Nat. Phys. 18, 1221–1227 (2022).
Lyu, Y.-Y. et al. Superconducting diode effect via conformal-mapped nanoholes. Nat. Commun. 12, 2703 (2021).
Suri, D. et al. Non-reciprocity of vortex-limited critical current in conventional superconducting micro-bridges. Appl. Phys. Lett. 121, 102601 (2022).
Ustavschikov, S. S. et al. Diode Effect in a Superconducting Hybrid Cu/MoN Strip with a Lateral Cut. J. Exp. Theor. Phys. 135, 226 (2022).
Hou, Y. et al. Ubiquitous Superconducting Diode Effect in Superconductor Thin Films. arXiv:2205.09276 (2022).
Diez-Merida, J. et al. Magnetic Josephson Junctions and Superconducting Diodes in Magic Angle Twisted Bilayer Graphene, arXiv:2110.01067 (2021).
Baumgartner, C. et al. Supercurrent rectification and magnetochiral effects in symmetric Josephson junctions. Nat. Nanotechnol. 17, 39 (2022).
Pal, B. et al. Josephson diode effect from cooper pair momentum in a topological semimetal. Nat. Phys. 18, 1228–1233 (2022).
Jeon, K. R. et al. Zero-field polarity-reversible Josephson supercurrent diodes enabled by a proximity-magnetized Pt barrier. Nat. Mater. 21, 1008–1013 (2022).
Gupta, M. et al. Superconducting Diode Effect in a Three-terminal Josephson Device. arXiv:2206.08471 (2022).
Golod, T. & Krasnov, V. M. Demonstration of a superconducting diode-with-memory, operational at zero magnetic field with switchable nonreciprocity. Nat. Commun. 13, 1 (2022).
Wu, H. et al. The field-free Josephson diode in a van der Waals heterostructure. Nature 604, 653 (2022).
Chiles, J. et al. Non-Reciprocal Supercurrents in a Field-Free Graphene Josephson Triode, arXiv:2210.02644 (2022).
Turini, B. et al. Josephson Diode Effect in High Mobility InSb Nanoflags. Nano Leters 22, 8502–8508 (2022).
Bardeen, J. Critical Fields and Currents in Superconductors. Rev. Mod. Phys. 34, 667 (1962).
Shabani, J. et al. Two-dimensional epitaxial superconductor-semiconductor heterostructures: A platform for topological superconducting networks. Phys. Rev. B 93, 155402 (2016).
Simoni, G. D., Paolucci, F., Solinas, P., Strambini, E. & Giazotto, F. Metallic supercurrent field-effect transistor. Nat. Nanotechnol. 13, 802 (2018).
Golokolenov, I., Guthrie, A., Kafanov, S., Pashkin, Y. & Tsepelin, V. On the origin of the controversial electrostatic field effect in superconductors. Nat. Commun. 12, 2747 (2021).
Winkler, G. W. et al. Unified numerical approach to topological semiconductor-superconductor heterostructures. Phys. Rev. B 99, 245408 (2019).
Giglberger, S. et al. Rashba and dresselhaus spin splittings in semiconductor quantum wells measured by spin photocurrents. Phys. Rev. B 75, 035327 (2007).
Farzaneh, S. M. et al. Magneto-anisotropic weak antilocalization in near-surface quantum wells. arXiv:2208.06050 (2022).
Baumgartner, C. et al. Effect of rashba and dresselhaus spin–orbit coupling on supercurrent rectification and magnetochiral anisotropy of ballistic josephson junctions. J. Phys.: Condens. Matter 34, 154005 (2022).
Geim, A. K. et al. Phase transitions in individual sub-micrometre superconductors. Nature 390, 259 (1997).
Nichele, F. et al. Scaling of Majorana Zero-Bias Conductance Peaks. Phys. Rev. Lett. 119, 1 (2017).
Acknowledgements
The Al/InAs heterostructures were provided by Michael Manfra Group at Purdue University. The authors thank Lev Ioffe and Igor Aleiner for stimulating discussions. Experimental part was initially supported by the U.S. Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering under Award DE-SC0008630, the work was completed with the support by NSF award DMR-DMR-2005092 (A.S. and L.P.R.). Theoretical work is supported by the U.S. Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering under Award DE-SC0010544 (Y.L-G) and the Office of the Under Secretary of Defense for Research and Engineering under award number FA9550-22-1-0354 (J.I.V.).
Author information
Authors and Affiliations
Contributions
L.P.R. conceived, A.S. performed experiments, J.I.V. and Y.L.G. developed the theory. All authors participated in writing the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Elia Strambini and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Sundaresh, A., Väyrynen, J.I., Lyanda-Geller, Y. et al. Diamagnetic mechanism of critical current non-reciprocity in multilayered superconductors. Nat Commun 14, 1628 (2023). https://doi.org/10.1038/s41467-023-36786-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-36786-5
This article is cited by
-
Link between supercurrent diode and anomalous Josephson effect revealed by gate-controlled interferometry
Nature Communications (2024)
-
Unconventional superconductivity in chiral molecule–TaS2 hybrid superlattices
Nature (2024)
-
DNA as a perfect quantum computer based on the quantum physics principles
Scientific Reports (2024)
-
Superconducting diode effect sign change in epitaxial Al-InAs Josephson junctions
Communications Physics (2024)
-
Sign reversal of the Josephson inductance magnetochiral anisotropy and 0–π-like transitions in supercurrent diodes
Nature Nanotechnology (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.