Abstract
Supercurrent diodes are nonreciprocal electronic elements whose switching current depends on their flow direction. Recently, a variety of composite systems combining different materials and engineered asymmetric superconducting devices have been proposed. Yet, ease of fabrication and tunable sign of supercurrent rectification joined to large efficiency have not been assessed in a single platform so far. We demonstrate that all-metallic superconducting Dayem nanobridges naturally exhibit nonreciprocal supercurrents under an external magnetic field, with a rectification efficiency up toâ~â27%. Our niobium nanostructures are tailored so that the diode polarity can be tuned by varying the amplitude of an out-of-plane magnetic field or the temperature in a regime without magnetic screening. We show that sign reversal of the diode effect may arise from the high-harmonic content of the current phase relation in combination with vortex phase windings present in the bridge or an anomalous phase shift compatible with anisotropic spin-orbit interactions.
Similar content being viewed by others
Introduction
Non-reciprocal charge transport is an essential element in modern electronics as a building block for multiple components such as rectifiers, photodetectors, and logic circuits. For instance, pn-junctions and Schottky-barrier devices are archetypal semiconductor-based examples of systems, known as diodes, with direction-selective charge propagation. Their operation stems from the spatial asymmetry of the heterojunction that provides inversion symmetry breaking. Likewise, dissipationless rectification refers to the asymmetric switching of the critical current (Isw) required to turn a superconductor into the normal state depending on the current bias polarity. Breaking both inversion and time-reversal symmetry, which are preserved in conventional superconductors, is the foundational aspect to enable the diode effect, as recently observed in superconducting materials1,2,3,4 and heterostructures5,6,7,8,9,10,11,12. Recent experimental findings have boosted a number of theoretical investigations in superconductors13,14,15,16 and Josephson junctions (JJs)17,18,19. In particular, several mechanisms have been proposed to account for the supercurrent diode effect (SDE). On the one hand, those based on intrinsic depairing currents focus on finite momentum pairing that arise from the combination of spin-orbit coupling (SOC) and Zeeman field13,14,15,16,20, or from Meissner currents18. On the other hand, other works underline the role of Abrikosov vortices, magnetic fluxes, screening currents or self-field effects as key elements for setting out nonreciprocal charge transport in superconductors21,22,23,24,25,26, such as in systems with trapped Abrikosov vortices27,28,29 or in micron-sized Nb-based strips30,31,32,33,34.
Till now, most of the research efforts have aimed to the realization of a SDE that maximizes the rectification efficiency, while the change of its polarity has been reported in a few cases only6,10,11,35,36,37. The SDE sign reversal has been interpreted as a consequence of finite momentum pairing6,11,35,36 requiring in-plane magnetic fields or diamagnetic currents and Josephson vortices10,29, as well as ascribed to vortex ratchet and asymmetric pinning effects38,39,40,41. All these outcomes suggest the need for an effective mastering of the polarity change of the SDE and its implementation in a simple and monolithic platform suitable for nanoscale miniaturization, not accomplished yet.
Here, we experimentally demonstrate a sign reversal tunable SDE in elemental superconducting weak links made on niobium (Nb). Nano-sized constrictions of Nb realize Dayem bridges whose switching currents for positive and negative sweep direction, \({I}_{sw}^{+}\) and \({I}_{sw}^{-}\), respectively, differ in the absolute value. This difference can be tuned both in amplitude and sign by an out-of-plane magnetic field (Bz), without inverting the polarity of Bz. Thermal effects can lead to two different energy scales for the maximal amplitude and the sign reversal of the diode efficiency which is inconsistent with physical scenarios related to Meissner currents or self-field effects. We show that sign reversal of the non-reciprocal response may arise from the phase shift due to the vortex phase winding or from spin-orbit effects due to the material granularity, in either case jointly with a few-harmonic content of the current-phase relation (CPR) of the weak link.
Results
Metallic diode architectures
We analyze two different geometries of Nb Dayem bridges, i.e., weak links made of a constant-thickness and all-metallic constriction between two superconducting banks42. The schematics of the electronic circuitry and false-color scanning electron micrographs of the devices are shown in Fig. 1a. In the first type of samples, 25-nm-thick micrometer-wide banks are connected via a link whose length l isâ~â80 nm and width wâ~â180 nm. The second type consists of 55-nm-thick banks connected via a quasi-one-dimensional wire with lâââ1 μm, and wâââ80 nm. Hereafter, we shall refer to the first and second types of bridges as short and long, respectively. Both device families are patterned through a single electron-beam lithography step followed by sputter deposition of the Nb thin film and lift-off. A 4-nm-thick Ti layer is pre-sputtered for adhesion purposes.
The differential resistance Râ=âdV/dI versus temperature T of two representative bridges is shown in Fig. 1b. The first abrupt reduction of R marks the critical temperature of the Nb films TTFâââ8.1(7.9) K for the 55(25)-nm-thick sample, suggesting minor impact of the film thickness on the superconducting state. The resistance drops to zero at the critical temperature of the weak link (Tc), which strongly depends, along with its normal-state resistance RN, on the geometry42. While the short bridge exhibits RNâ~â40âΩ, in the long one has RNâ~â270âΩ.
Below Tc, dissipationless transport occurs in the bridges owing to Cooper pairs supercurrent. The temperature dependence of the switching current Isw of both devices is displayed in Fig. 1c. From the fit to the Bardeen equation \({I}_{sw}(T)={I}_{sw}^{0}{[1-{(\frac{T}{{T}_{c}})}^{2}]}^{\frac{3}{2}}\), we extract a zero-temperature switching current \({I}_{sw}^{0}\simeq 720\,\)μA and a critical temperature \({T}_{c}^{S}\simeq 4.3\) K for the shor bridge. Similarly, for the long weak link we obtain \({I}_{sw}^{0}\simeq 42\,\mu\)A and \({T}_{c}^{L}\simeq 2.1\) K. From these values, we determine a zero-temperature BCS energy gap \({\Delta }_{0}=1.764{k}_{B}{T}_{c}^{S(L)}\simeq 650(320)\,\mu\) eV for the short(long) bridge, where kB is the Boltzmann constant. For the long bridge, we deduce a superconducting coherence length \({\xi }_{0}=\sqrt{\hslash l/({R}_{N}wt{e}^{2}{N}_{F}{\Delta }_{0})}\simeq 11\) nm, where t is the film thickness, NFâââ5.33âÃâ1047Jâ1mâ3 is the density of states at the Fermi level of Nb43, and e is the electron charge. Similarly, we can evaluate the London penetration depth \({\lambda }_{L}=\sqrt{\hslash {R}_{N}wt/(\pi l{\mu }_{0}{\Delta }_{0})}\simeq\) 790 nm, where μ0 is the vacuum magnetic permeability. Since w,âtââªâλL, the bridges can be uniformly penetrated by an external magnetic field.
The current vs voltage (IV) characteristics of the short and long bridges are shown in Fig. 1d, e, respectively, for selected values of bath temperature. The devices show an abrupt transition to the normal state at the switching current Isw and display the typical hysteresis of metallic junctions which originates from Joule heating induced in the bridge when the bias current is swept back from the resistive to the dissipationless state44.
Short Dayem bridge diode performance
Let us now discuss how the short Dayem bridge in Fig. 1a can be used as a supercurrent diode. Non-reciprocal dissipationless transport is revealed by comparing the switching currents while sweeping the biasing current from zero to positive values (\({I}_{sw}^{+}\)) or from zero to negative values (\({I}_{sw}^{-}\)) in the presence of an out-of-plane magnetic field Bz. The switching currents at Tâ=â0.3 K are reported in Fig. 2a. The magnetic field increasingly reduces the superconducting gap and thereby both the switching currents. A linear decrease in Bz of \({I}_{sw}^{+}\) and \(| {I}_{sw}^{-}|\) is observed up to ~0.07 T. At larger fields, the dependence of the switching currents on Bz is sublinear. Notably, both Isw(Bz) are not antisymmetric with respect to the magnetic field (Isw(Bz)ââ âIsw(âââBz)) while the symmetry relation \({I}_{sw}^{+}({B}_{z})\simeq -{I}_{sw}^{-}(-{B}_{z})\) is respected, within the small experimental fluctuations, as theoretically expected. This symmetry relation is further confirmed in the switching currents difference \(\Delta {I}_{sw}\equiv {I}_{sw}^{+}-| {I}_{sw}^{-}|\) displaying an odd-in-Bz superconducting diode effect (ÎIsw(Bz)âââââÎIsw(âââBz)) as shown in Figure 2b. ÎIsw is characterized by a maximum at Bmaxâââ0.05 T and a sign inversion at BRâââ0.1 T where \({I}_{sw}^{+}\) and \(| {I}_{sw}^{-}|\) have a crossing (see Fig. 2a). From now on, Bmax indicates the position in field of the rectification peak. Two IV curves, recorded for magnetic fields lower and larger than BR, are plotted in Fig. 2c to emphasize the sign change in the rectification.
Nonreciprocal transport can be conveniently quantified by the rectification efficiency defined as \(\eta =\frac{{I}_{sw}^{+}-| {I}_{sw}^{-}| }{{I}_{sw}^{+}+| {I}_{sw}^{-}| }\). Figure 2d shows the evolution of η versus Bz and T. η(Bz) is substantially unaffected by thermal effects up to Tâââ1.75 K \(=0.41\,{T}_{c}^{S}\) where a maximum rectification ηmaxâ~â27% is obtained. The evolution of ηmax in temperature is displayed in Fig. 2e (i). In addition, we parametrize the diode sensitivity to the magnetic field in the vicinity of the abrupt sign change as Îâ=âηmax/(â£BmaxâââBRâ£). A maximum value Îâ~â650 Tâ1 is achieved around 2.25 K (see Fig. 2e (i)) At higher temperatures, the quantities ηmax Î, and the characteristic magnetic fields Bmax and BR related to rectification (see Fig. 2e (ii)), all decrease in a similar fashion. The full profile of the rectification efficiency versus Bz is better visualized in Fig. 2f where η(Bz) is plotted for a few selected values of temperature.
Long Dayem bridge diode performance
Next, we characterize the long nanobridge shown in Fig. 1a. Figure 3a reports the decay of \({I}_{sw}^{+}\) and \(| {I}_{sw}^{-}|\) as a function of Bz. At first, we notice that the switching currents are damped down toâ~â60% of their zero-field value at Bzâââ0.3 T, whereas in the previous sample, the same damping is achieved for lower fields (Bzâââ0.07 T see Fig. 2a). Figure 3b displays ÎIsw versus Bz. For low magnetic fields, ÎIsw(Bz) exhibits a linear relation. While increasing â£Bz⣠further, ÎIsw bends and then inverts its trend: an abrupt jump realizes a sign reversal at â£BRâ£â~â0.34 T. Then, a relative peak at â£Bmaxâ£âââ0.38 T, marked by a red dashed line, represents the field at which maximum rectification efficiency is achieved as before. Next to the sign change, ÎIsw leaves the clean trend and looks noisy. Such small jumps are reproducible, thus ruling out a stochastic nature of the underlying processes. Finally, ÎIsw oscillates at higher magnetic fields, as shown in the inset of Fig. 3b. Two IV curves, for fields lower and larger than BR, are plotted in Fig. 3c to highlight that the rectification sign changes from negative to positive as the field Bzâ>â0 increases contrary to the short bridge. This change in symmetry is attributed to the vortex nucleation as discussed later.
The magnetic field and temperature dependence of the rectification efficiency η is presented as a color plot in Fig. 3d. The sign change and the maximum rectification are affected by temperature in a different way as compared to the short constriction. The linear increase of the rectification at low fields smears out with temperature, reducing BR until it vanishes at Tâââ1.1 K \(=0.5\,{T}_{c}^{L}\). Rectification parameters displayed in Fig. 3e shows that Bmax is more robust in temperature than BR: it is still observable at Tâââ1.8 K \(=0.8\,{T}_{c}^{L}\). The sudden change of sign is quantified by a maximum Îâ~â360 Tâ1 at 0.15 K. As before, the profile of rectification efficiency as a function of Bz is shown in Fig. 3f for a few selected values of temperature. The difference in the temperature trend between BR and Bmax (see Fig. 3e (ii)) suggests two different energy scales responsible for the sign reversal and the maximum rectification, as confirmed by measurements obtained in other devices with same nominal dimensions. (see Supplementary Figure 1). Rectification on the second sample exhibits similar η(Bz) lineshape with almost identical Bmax and ηmax values and temperature dependence. In this sample, low-field features fade more rapidly with temperature, which appears to be sample-dependent.
Modeling the sign reversal of the diode effect
Previous SDE reported in conventional superconductor are normally attributed to Abrikosov vortices and self-field effects29,34 or Meissner currents30,31,32,33. We have neglected those mechanisms since they are weak in amplitude. Such a conclusion can be directly deduced from the stark contrast of the profile of the switching current vs B-field data points. In both scenarios Isw(Bz) reaches a maximum at nonzero magnetic fields. Moreover, amplitude asymmetry does not change sign as a function of the magnetic field without reversing its orientation before the critical current vanishes. Instead, the switching current-field profile in our devices decays monotonically at low fields with no signature of the characteristic skewness. For this reason, we propose two physical scenarios compatible with our devices that may explain our experimental findings. Both of them rely on non-sinusoidal CPRs, typical of superconducting nanobridges42, combined with an inversion-symmetry breaker. In model I, this is represented by a supercurrent vortex, while in model II by spin-orbit couplings. An out-of-plane magnetic field parametrized by B* is considered, where B* is given by the zero rectification field η(B*)â=â0 and by the vanishing of high-harmonic amplitude.
In model I, the Dayem nanobridge is schematized as a one-dimensional chain of weak links of width w formed by the Nb grains. A supercurrent vortex can nucleate in one of these weak links45,46, as sketched in Fig. 4a, which induces a phase winding in the superconducting order parameter. These vortices have a typical size of the order of ξ, so only a few of them can be accommodated within the bridge. Notice that such vorticity is not a screening of the B-field, since the small dimensions of the bridge (wââªâλL) allow full penetration of Bz. In this framework, the CPR is affected by two-phase shifts: the conventional vector potential associated with Bz and the phase winding of the vortex. It is indeed the interplay of these two contributions that is responsible for a sign change in the rectification parameter. Though on a different length scale, this physical scenario is similar to that of Josephson phase vortices29,47. The rectification parameter η is then evaluated by determining the maximum and minimum values of the Josephson current with respect to the phase bias, see Methods for details.
Figure 4b, reports the evolution of η in B for different amplitudes of the second harmonic, I2, of the CPR. The magnitude of η scales with I2 showing multiple nodes, whose position in B is independent of I2. The sign change also depends on the position of the vortex, as displayed in Fig. 4c, where η(B) is evaluated for a vortex nucleated at different distances from the lateral edge of the bridge (xν). This behavior suggests a phase-shift competition dominated at low fields by the vortex phase slip and at large fields by the vector potential.
Another scenario able to describe the sign reversal in the diode rectification can be envisioned by combining the colored CPR of the nanobridge with an anomalous phase shift14,48,49,50,51 induced by spin-orbit interactions and magnetic fields. Nb is a heavy metal and thus possesses an atomistic spin-orbit interaction of the order ofâ~â100 meV52. We extend such concept to our weak links where we call for a spin-orbit coupling of the order of a few meV, i.e., 2 orders of magnitude weaker than the atomic limit53. In particular, we can expect that mirror symmetry can be locally or globally broken in polycrystalline films54, thereby leading to spin-orbit interaction of both Rashba and Dresselhaus types (see Methods). Figure 4d shows a sketch of the bridge modeled as an effective \(SS{\prime} S\) structure, where the S and \(S{\prime}\) components have different amplitudes of the superconducting gap and different spin-orbit couplings breaking horizontal and vertical mirror symmetries. The anisotropic spin-orbit interaction generates an anomalous phase shift in the CPR that varies with the magnetic field, as explicitly shown in Supplementary Figure 2. The anomalous phase is then introduced in the CPR via a phenomenological parameter ÎB providing a first-order cosine component in the Fourier expansion, i.e., \(I={\sum }_{n}{I}_{n}\sin (\varphi )+{\Gamma }_{B}\cos (\varphi )\) (see Methods for details). The anomalous phase Ï0 is then related to the amplitude of ÎB. We assume a linear damping of the harmonics Inâ=âIn,0(1âââB/B*), which defines the scale parameter B*.
From model II, the diode sign reversal takes place only in the presence of a sizable third harmonic component. Figure 4e reports η(B) for some values of I3,0. By increasing I3,0, the sign inversion gets more pronounced, whereas the maxima and minima of the rectification (\({B}_{max}^{s}\simeq 0.74{B}^{* }\), \({B}_{{{{{{{\rm{min}}}}}}}}^{s}\simeq 0.89{B}^{* }\)) are barely affected by the weight of the harmonic. Moreover, by including more harmonics in the CPR, the lineshape of η(B) is modified. For instance, Figure 4f shows that a fourth-order harmonic affects the magnetic field dependence by substantially removing the sign change.
Discussion
The comparison between our experimental findings and the proposed models reveals some important features supporting the proposed mechanisms. In particular, for both bridges η has an almost monotonic damping in temperature that can be explained in both models by the reduction of high-order harmonics. This is expected in long metallic weak links where the CPR evolves from highly distorted to sinusoidal-like shapes at large temperatures42,55. Moreover, as shown in Figs. 2e and 3e, Bmax is temperature resilient until Tâ¥0.5âTc. This feature is fairly captured by both models, as the maximum rectification looks almost independent of the harmonic content (see Fig. 4b for model I and Fig. 4e for model II). However, long bridges exhibit features that are mostly accounted by model I, while short ones are more compatible with model II. For example, in long (short) bridges the sign reversal is present below (above) Bmax as shown in Fig. 3b and Fig. 2b to be compared with Fig. 4b and e, respectively. Multiple sign reversal nodes appear in high fields only for long bridges, as shown in Fig. 3b and well described by the interferometric mechanism of model I. While the rectification lineshape given by model II presents only one inversion node. Moreover, the quick damping of the rectification inversion observed at low fields (â<â0.3 T) in Fig. 3f is captured by the vortex dynamics described in Fig. 4c. The relative size of the vortex ξ/w is temperature-dependent and influences the vortex position in the bridge. Thus, it is plausible to expect that variations of the vortex size mostly affect the rectification lineshape at low fields while remaining substantially unchanged for larger fields, as shown in Fig. 4c.
Finally, it is interesting to note that by extending the proposed models to in-plane magnetic fields, a sizable supercurrent rectification is anticipated but without sign reversal. In particular, for model I, no phase shift is expected from the spatial dependence of the vector potential, since the orbital coupling between an in-plane field and the electron momentum becomes negligible. Thus, the source of phase interference with the vortex winding is eliminated. For model II, the anomalous phase and the harmonic content would be differently affected by an in-plane Zeeman field compared to the out-of-plane orientation not resulting in sign reversal.
Conclusions
In summary, we have demonstrated the implementation of supercurrent diodes in Nb Dayem nanobridges. By breaking the time-reversal symmetry with an out-of-plane magnetic field, we demonstrate that both the amplitude and sign of the rectification can be tuned without inverting the polarity of the applied field. The evolution of the critical current on B field decays monotonically without showing any skewness, which rules out self-field effects or Meissner currents. We have developed two theoretical models to account for the sources of time- and inversion-symmetry breaking, one based on a vortex phase winding, and one that takes into account the spin-orbit interactions present in polycrystalline heavy materials. Yet, a quantitative description of the supercurrent diode effect in metallic nanoconstrictions should account for both scenarios, which complement each other and can coexist. Furthermore, the fabrication process is simple when compared to that of other platforms, a compelling step towards scalability. Analogous nanobridges can be realized from several elemental superconductors currently at the base of other architectures, such as nanocryotrons56, rapid single-flux quanta (RSFQ)57 and memories58, which would ease a potential integration.
Finally, the sharp sign reversal of the diode rectification allows us to envisage applications of Dayem nanobridges as B-field threshold detectors. When biased in the vicinity of the rectification node, small variations of an environmental magnetic field would result in modifications of the sign of the rectification parameter.
Methods
Sample fabrication
Nb strips and constrictions are patterned by e-beam lithography on AR-P 679.04 (PMMA) resist. PMMA residuals are removed by O2-plasma etching after developing. Nb thin films were deposited by sputtering with a base pressure of 2âÃâ10â8 Torr in a 4 mTorr Ar (6N purity) atmosphere and liftoff by acetone or AR-P 600.71 remover. A thin Ti layer was previously sputtered to improve Nb adhesion and base pressure in the deposition chamber, resulting in a nominal thickness Ti(4nm)/Nb(25nm) and Ti(4nm)/Nb(55nm) for the so-called short and long nanobridges, respectively.
Transport measurements
Transport measurements were carried out in filtered (two-stage RC and Ï filters) cryogen-free 3He-4He dilution refrigerators by a standard 4-point probe technique. DC current-voltage characteristics were measured by sweeping a low-noise current bias positively and negatively, and by measuring the voltage drop across the weak links with a room-temperature, low-noise preamplifier for each current value every ~20 ms. Current step size ÎI was adapted depending on the switching current keeping the values lower than 0.002Isw(Bâ=â0) ~1μ A and 100 nA for short and long devices, respectively. The switching currents are extracted from the maximum of the derivative dV/dI and averaged from 5-10 reiterations of the IV curves. Error bars account for the standard deviation and propagation of errors in Isw, ÎIsw, and η. For the characteristic fields, the error is given by the magnetic field step ÎBzâ<â10 mT. Joule heating is minimized by automatically switching the current off once the device turns into the normal state. A delay between sweeps was optimized to keep the stability of the fridge temperature lower than 50 mK. Furthermore, no changes in the switching currents (up to the accuracy given by the standard deviation) were observed in different cooling cycles, by changing the order of the sweeps or by adding an extra delay in the acquisition protocol concluding that hysteretic behavior or local heating is negligible.
Theoretical models
The current phase relation for the model I is given by \(I={\sum }_{n = 1,2}\int\nolimits_{0}^{w}{I}_{n}\sin \left(\phi (x)+n\varphi \right)dx\). The critical current in both directions is evaluated by determining the maximum and minimum values with respect to the phase bias Ï. Here, the spatial-dependent phase difference Ï(x) is given by the magnetic field and the phase winding contribution due to the vortex47. We have assumed that the supercurrent has a subdominant second harmonic contribution, as expected in long, diffusive weak links42,55. The spatial dependence of the phase along the transverse direction x is expressed as \(\phi (x)=\frac{2\pi {d}_{b}{B}_{z}x}{{{{\Phi }}}_{0}}+{\phi }_{v}(x)\) where db is a characteristic length of the weak link, related to the width of the junction, and \({\phi }_{v}(x)={\gamma }_{v}\arctan [\frac{(x-{x}_{v})}{{y}_{v}}]\) is the spatially-inhomogeneous phase offset due to the vortex structure. γv indicates the sign and amplitude of the winding and (xv,âyv) the position of the vortex core with respect to the boundaries (see Fig. 4a). We also assume that the second harmonic amplitude of the supercurrent is vanishing at Bâ=âB*.
Concerning model II, the current phase relation is given by \(I={\sum }_{n}{I}_{n}\sin (\varphi )+{\Gamma }_{B}\cos (\varphi )\). The critical currents are evaluated by determining the maximum and minimum values of the Josephson current. For the examined model, we assume that the amplitudes of the n-harmonic In get suppressed by B with a linear rate, i.e. \({I}_{n}={I}_{n,0}(1-\frac{B}{{B}^{* }})\) for nâ¥1. The linear trend is compatible with the observed behavior of the overall supercurrent amplitude as shown in Figs. 2a and 3a for the range of applied field where the rectification is nonvanishing. Furthermore, the suppression of the amplitude of the harmonics with the magnetic field can be related to the reduction of the effective transmission across the grains due to depairing and magnetic interference. The decrease in the transmission implies a reduction of the non-harmonic amplitudes42,55. In this model, we have performed a real space simulation for the examined geometry. By solving the Bogoliubov-de Gennes equations on a finite-size slab in the presence of an out-of-plane magnetic field, we demonstrate that an anomalous phase shift can be obtained. The simulation is performed for a system size NxâÃâNy with Nxâ=â150 and Nyâ=â100. The employed tight-binding model includes a nearest neighbour hopping amplitude t, and the conventional spin-singlet local pairing amplitude. We apply a phase bias across the weak link and determine the free energy as shown in the Supplementary Figure 2a. The resulting anomalous phase increases with the magnetic field and depends on the strength of the Rashba and Dresselhaus interactions as shown in the Supplementary Figure 2b. The linear Rashba term on a lattice for a two-dimensional geometry is expressed as \({H}_{R}={\alpha }_{R}[\sin ({k}_{x}){\sigma }_{y}-\sin ({k}_{y}){\sigma }_{x}]\) while the Dresselhaus term is given by \({H}_{D}={\alpha }_{D}[\sin ({k}_{x}){\sigma }_{x}-\sin ({k}_{y}){\sigma }_{y}]\) with Ïi(iâ=âx,ây) being the Pauli matrices associated with the spin angular momentum. We notice that HR breaks the horizontal mirror symmetry while HD breaks both the vertical and horizontal mirror symmetries. This reduction of mirror symmetry is expected to be locally or globally broken in granular films54. Notably, the presence of the Dresselhaus term is crucial to induce an anomalous phase shift in the supercurrent in the presence of an out-of-plane magnetic field.
Concerning the size of the Rashba and Dresselhaus couplings, our analysis shows a sizable, nonvanishing anomalous phase in the presence of an out-of-plane magnetic field for values from 0.01 to 0.3, in units of the effective electron hopping amplitude t. We notice that t is an effective electronic amplitude that describes the low energy properties of the superconductor which can be taken in the order of 20 meV. For such value, Rashba/Dresselhaus couplings able to generate the anomalous phase will range from 0.2 to 6 meV. Although there is no experimental estimations of Rashba and Dresselhaus coupling for Nb thin films, these values are plausible for heavy metals like Nb, since Rashba spin splitting scales with the size of the atomic spin-orbit coupling59,60. Nb has a multiorbital electronic structure61, in this case, the breaking of the inversion or mirror symmetry primarily leads to an orbital Rashba coupling62,63 among orbitals with different mirror parity. Then, the spin Rashba interaction arises through the atomic spin-orbit coupling. The large amplitude of the atomic spin-orbit coupling in Nb can thus yield a sizable spin Rashba interaction 59,60,64,65.
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Code availability
The code that support the findings of this study are available from the corresponding author upon reasonable request.
References
Wakatsuki, R. et al. Nonreciprocal charge transport in noncentrosymmetric superconductors. Sci. Adv. 3, e1602390 (2017).
DÃez-Mérida, J. et al. Symmetry-broken Josephson junctions and superconducting diodes in magic-angle twisted bilayer graphene. Nat. Commun. 14, 2396 (2023).
Lin, J.-X. et al. Zero-field superconducting diode effect in small-twist-angle trilayer graphene. Nat. Phys. 18, 1221â1227 (2022).
Paolucci, F., De Simoni, G. & Giazotto, F. A gate-and flux-controlled supercurrent diode effect. Appl. Phys. Lett. 122, 042601 (2023).
Ando, F. et al. Observation of superconducting diode effect. Nature 584, 373â376 (2020).
Pal, B. et al. Josephson diode effect from cooper pair momentum in a topological semimetal. Nat. Phys. 18, 1228â1233 (2022).
Bauriedl, L. et al. Supercurrent diode effect and magnetochiral anisotropy in few-layer NbSe2. Nat. Commun. 13, 4266 (2022).
Baumgartner, C. et al. Supercurrent rectification and magnetochiral effects in symmetric Josephson junctions. Nat. Nanotechnol. 17, 39â44 (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).
Sundaresh, A., Väyrynen, J. I., Lyanda-Geller, Y. & Rokhinson, L. P. Diamagnetic mechanism of critical current non-reciprocity in multilayered superconductors. Nat. Commun. 14, 1628 (2023).
Lotfizadeh, N. et al. Superconducting diode effect sign change in epitaxial Al-InAs Josepshon junctions. arXiv:2303.01902 (2023).
Turini, B. et al. Josephson diode effect in high-mobility InSb nanoflags. Nano Lett. 22, 8502â8508 (2022).
Daido, A., Ikeda, Y. & Yanase, Y. Intrinsic superconducting diode effect. Phys. Rev. Lett. 128, 037001 (2022).
IliÄ, S. & Bergeret, F. S. Theory of the supercurrent diode effect in Rashba superconductors with arbitrary disorder. Phys. Rev. Lett. 128, 177001 (2022).
Yuan, N. F. Q. & Fu, L. Supercurrent diode effect and finite-momentum superconductors. Proc. Natl Acad. Sci. 119, e2119548119 (2022).
He, J. J., Tanaka, Y. & Nagaosa, N. A phenomenological theory of superconductor diodes. New J. Phys. 24, 053014 (2022).
Misaki, K. & Nagaosa, N. Theory of the nonreciprocal Josephson effect. Phys. Rev. B 103, 245302 (2021).
Davydova, M., Prembabu, S. & Fu, L. Universal Josephson diode effect. Sci. Adv. 8, eabo0309 (2022).
Zhang, Y., Gu, Y., Li, P., Hu, J. & Jiang, K. General theory of Josephson diodes. Phys. Rev. X 12, 041013 (2022).
Scammell, H. D., Li, J. I. A. & Scheurer, M. S. Theory of zero-field superconducting diode effect in twisted trilayer graphene. 2D Materials 9, 025027 (2022).
Wambaugh, J. F., Reichhardt, C., Olson, C. J., Marchesoni, F. & Nori, F. Superconducting fluxon pumps and lenses. Phys. Rev. Lett. 83, 5106â5109 (1999).
Vodolazov, D. Y. & Peeters, F. M. Superconducting rectifier based on the asymmetric surface barrier effect. Phys. Rev. B 72, 172508 (2005).
Villegas, J. E., Gonzalez, E. M., Gonzalez, M. P., Anguita, J. V. & Vicent, J. L. Experimental ratchet effect in superconducting films with periodic arrays of asymmetric potentials. Phys. Rev. B 71, 024519 (2005).
Van de Vondel, J., de Souza Silva, C. C., Zhu, B. Y., Morelle, M. & Moshchalkov, V. V. Vortex-rectification effects in films with periodic asymmetric pinning. Phys. Rev. Lett. 94, 057003 (2005).
Carapella, G., Granata, V., Russo, F. & Costabile, G. Bistable Abrikosov vortex diode made of a Py-Nb ferromagnet-superconductor bilayer structure. Appl. Phys. Lett. 94, 242504 (2009).
Cerbu, D. et al. Vortex ratchet induced by controlled edge roughness. New J. Phys. 15, 063022 (2013).
Golod, T. et al. Reconfigurable Josephson Phase Shifter. Nano Lett. 21, 5240â5246 (2021).
Lyu, Y.-Y. et al. Superconducting diode effect via conformal-mapped nanoholes. Nat. Commun. 12, 2703 (2021).
Golod, T. & Krasnov, V. M. Demonstration of a superconducting diode-with-memory, operational at zero magnetic field with switchable nonreciprocity. Nat. Commun. 13, 3658 (2022).
Suri, D. et al. Non-reciprocity of vortex-limited critical current in conventional superconducting micro-bridges. Appl. Phys. Lett. 121, 102601 (2022).
Chahid, S., Teknowijoyo, S., Mowgood, I. & Gulian, A. High-frequency diode effect in superconducting \({{{{{{{{{{\rm{Nb}}}}}}}}}}}_{3}{{{{{{{{{\rm{Sn}}}}}}}}}}\) microbridges. Phys. Rev. B 107, 054506 (2023).
Hou, Y. et al. Ubiquitous superconducting diode effect in superconductor thin films. Phys. Rev. Lett. 131, 027001 (2023).
Satchell, N., Shepley, P. M., Rosamond, M. C. & Burnell, G. Supercurrent diode effect in thin film Nb tracks. J. Appl. Phys. 133, 203901 (2023).
Vávra, O., Pfaff, W., Monaco, R., Aprili, M. & Strunk, C. Current-controllable planar S-(S/F)-S Josephson junction. Appl. Phys. Lett. 102, 072602 (2013).
Costa, A. et al. Sign reversal of the Josephson inductance magnetochiral anisotropy and 0-Ï-like transitions in superconducting diodes. Nat. Nanotechnol.1â7 (2023).
Kawarazaki, R. et al. Magnetic-field-induced polarity oscillation of superconducting diode effect. Appl. Phys. Express. 15, 113001 (2022).
Gupta, M. et al. Gate-tunable superconducting diode effect in a three-terminal Josephson device. Nat. Commun. 14, 3078 (2023).
Gillijns, W., Silhanek, A. V., Moshchalkov, V. V., Reichhardt, C. J. O. & Reichhardt, C. Origin of reversed vortex ratchet motion. Phys. Rev. Lett. 99, 247002 (2007).
He, A., Xue, C. & Zhou, Y.-H. Switchable reversal of vortex ratchet with dynamic pinning landscape. Appl. Phys. Lett. 115 (2019).
Ideue, T., Koshikawa, S., Namiki, H., Sasagawa, T. & Iwasa, Y. Giant nonreciprocal magnetotransport in bulk trigonal superconductor PbTaSe2. Phys. Rev. Res. 2, 042046 (2020).
Jiang, J. et al. Reversible ratchet effects in a narrow superconducting ring. Phys. Rev. B 103, 014502 (2021).
Likharev, K. K. Superconducting weak links. Rev. Mod. Phys. 51, 101â159 (1979).
Jani, A. R., Brener, N. E. & Callaway, J. Band structure and related properties of bcc niobium. Phys. Rev. B 38, 9425â9433 (1988).
Skocpol, W. J., Beasley, M. R. & Tinkham, M. Self-heating hotspots in superconducting thin-film microbridges. J. Appl. Phys. 45, 4054â4066 (2003).
Aranson, I., Gitterman, M. & Shapiro, B. Y. Motion of vortices in thin superconducting strips. J. Low Temp. Phys. 97, 215â228 (1994).
Aranson, I., Gitterman, M., Shapiro, B. Y. & Vinokur, V. Nucleation, growth and kinetics of the vortex phase in thin superconducting strips. Phys. Scrip. 1996, 125 (1996).
Krasnov, V. M. Josephson junctions in a local inhomogeneous magnetic field. Phys. Rev. B 101, 144507 (2020).
Assouline, A. et al. Spin-orbit induced phase-shift in Bi2Se3 Josephson junctions. Nat. Commun. 10, 126 (2019).
Strambini, E. et al. A Josephson phase battery. Nat. Nanotechnol. 15, 656â660 (2020).
Mayer, W. et al. Gate controlled anomalous phase shift in Al/InAs Josephson junctions. Nat. Commun. 11, 212 (2020).
Margineda, D., Claydon, J. S., Qejvanaj, F. & Checkley, C. Observation of anomalous Josephson effect in nonequilibrium Andreev interferometers. Phys. Rev. B 107, L100502 (2023).
Machin, D. J. & Sullivan, J. F. The preparation and properties of complexes derived from niobium(IV) halides. J. less-common met. 19, 405â411 (1969).
Gül, Ã. et al. Andreev Reflection in the Fractional Quantum Hall State. Phys. Rev. X 12, 021057 (2022).
Park, J.-S. Stabilization and self-passivation of symmetrical grain boundaries by mirror symmetry breaking. Phys. Rev. Mater. 3, 014602 (2019).
Golubov, A. A., Kupriyanov, M. Y. & Ilâichev, E. The current-phase relation in Josephson junctions. Rev. Mod. Phys. 76, 411â469 (2004).
McCaughan, A. N. & Berggren, K. K. A superconducting-nanowire three-terminal electrothermal device. Nano Lett. 14, 5748â5753 (2014).
Likharev, K. K. & Semenov, V. K. Rsfq logic/memory family: A new Josephson-junction technology for sub-terahertz-clock-frequency digital systems. IEEE Trans. Appl. Supercond. 1, 3â28 (1991).
Ligato, N., Strambini, E., Paolucci, F. & Giazotto, F. Preliminary demonstration of a persistent Josephson phase-slip memory cell with topological protection. Nat. Commun. 12, 5200 (2021).
Sunko, V. et al. Maximal Rashba-like spin splitting via kinetic-energy-coupled inversion-symmetry breaking. Nature 549, 492â496 (2017).
Bihlmayer, G., Noël, P., Vyalikh, D. V., Chulkov, E. V. & Manchon, A. Rashba-like physics in condensed matter. Nat. Rev. Phys. 4, 642â659 (2022).
RüÃmann, P. & Blügel, S. Density functional Bogoliubov-de Gennes analysis of superconducting Nb and Nb(110) surfaces. Phys. Rev. B. 105, 125143 (2022).
Park, S. R., Kim, C. H., Yu, J., Han, J. H. & Kim, C. Orbital-Angular-Momentum Based Origin of Rashba-Type Surface Band Splitting. Phys. Rev. Lett. 107, 156803 (2011).
Mercaldo, M. T., Solinas, P., Giazotto, F. & Cuoco, M. Electrically Tunable Superconductivity Through Surface Orbital Polarization. Phys. Rev. Appl. 14, 034041 (2020).
Ishizaka, K. et al. Giant Rashba-type spin splitting in bulk BiTeI. Nat. Mater. 10, 521â526 (2011).
Yaji, K. et al. Large Rashba spin splitting of a metallic surface-state band on a semiconductor surface. Nat. Commun. 1, 17 (2010).
Acknowledgements
This work was funded by the EUâs Horizon 2020 Research and Innovation Framework Program under Grant Agreement No. 964398 (SUPERGATE), No. 101057977 (SPECTRUM), and by the PNRR MUR project PE0000023-NQSTI.
Author information
Authors and Affiliations
Contributions
D.M. fabricated the samples, conducted the experiments, and analyzed data with inputs from A.C., E.S. and F.G. The theoretical models describing the experiment was developed by Y.F., M.T.M., and M.C. The manuscript was written by D.M., A.C., E.S., M.C., and F.G. with inputs from all the authors. E.S. and F.G. conceived the experiment. F.G. supervised and coordinated the project. All authors discussed the results and their implications equally at all stages.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Communications Physics thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Margineda, D., Crippa, A., Strambini, E. et al. Sign reversal diode effect in superconducting Dayem nanobridges. Commun Phys 6, 343 (2023). https://doi.org/10.1038/s42005-023-01458-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s42005-023-01458-9
This article is cited by
-
Superconducting diode effect sign change in epitaxial Al-InAs Josephson junctions
Communications Physics (2024)
-
Large superconducting diode effect in ion-beam patterned Sn-based superconductor nanowire/topological Dirac semimetal planar heterostructures
Nature Communications (2024)