Josephson coupling in the dissipative state of a thermally hysteretic \mu-SQUID

Josephson coupling in the dissipative state of a thermally hysteretic -Squid

Sourav Biswas Department of Physics, Indian Institute of Technology Kanpur, Kanpur 208016, India    Clemens B. Winkelmann Univ. Grenoble Alpes, CNRS, Grenoble INP, Institut Néel, 38000 Grenoble, France    Hervé Courtois Univ. Grenoble Alpes, CNRS, Grenoble INP, Institut Néel, 38000 Grenoble, France    Anjan K. Gupta Department of Physics, Indian Institute of Technology Kanpur, Kanpur 208016, India
July 26, 2019

Micron-sized superconducting interference devices (-SQUIDs) based on constrictions optimized for minimizing thermal runaway are shown to exhibit voltage oscillations with applied magnetic flux despite their hysteretic behavior. We explain this remarkable feature by a significant supercurrent contribution surviving deep into the resistive state, due to efficient heat evacuation. A resistively shunted junction model, complemented by a thermal balance determining the amplitude of the critical current, describes well all experimental observations, including the flux modulation of the (dynamic) retrapping current and voltage by introducing a single dimensionless parameter. Thus hysteretic -SQUIDs can be operated in the voltage read-out mode with a faster response. The quantitative modeling of this regime incorporating both heating and phase dynamics paves the way for further optimization of -SQUIDs for nano-magnetism.

I Introduction

A superconducting quantum interference device (SQUID), in which two Josephson junctions form a closed loop, exhibits a modulation of the critical current , as a function of the magnetic flux through the loop with a period . It is the most sensitive magnetic field transducer to date.granatasquid (); recentsquidreview () Miniaturized SQUIDs have been used for probing magnetic properties of nanoparticles mag (); mag1 (); cobultnano () and surfaces with sub-m resolution.scanningsquidmicro1 (); scanningsquidmicro2 () The coupling of a nano-particle’s magnetic flux to a - or nano-SQUID is far better granatasquid () than to a conventional SQUID, leading to a magnetic moment resolution down to below 1 .mag1 (); mag2 () Thus, optimizing such -SQUIDs in terms of sensitivity, ease of fabrication and operation, or operating temperature and magnetic field range is the focus of a large panel of recent works.highqualitysquid (); siliconsquid (); nbnanosquid (); diamondsquid (); lownoisesquid (); NbNsquid (); giazotto () However, hysteresis in the current-voltage characteristics (IVCs) of -SQUIDs severely limits their flux resolution and speed. Hysteresis in conventional SQUIDs based on superconductor-insulator-superconductor Josephson junctions is well understood in the frame of the resistively and capacitively shunted junction (RCSJ) model.squidbook (); tinkhambook () In contrast, hysteresis in weak-link (WL) based -SQUIDs arises from the Joule heating leading to a self-sustained hot-spot.scocpol (); tinkhamnanowire (); herveprl (); hazraprb () Various strategies to avoid this phenomenon have been exploredbaranovprb (); nikhilsust (); lahmapl () but with incomplete success. Recently, high sensitivity non-hysteretic nano-SQUIDs, based on high critical temperature cuprate grain boundary junctions shunted with gold films, have been demonstrated for nano-magnetism.hightc1 (); hightc2 () Although such devices have some advantage as heating effects seem to be negligible there, their fabrication is quite challenging.

In a WL biased with a current close to the critical current , the transition to the dissipative state is triggered by a phase slipgregory (); nayanaprl (), which changes the quantum phase difference between the leads by . The ensuing voltage peak, and thus heating, generally suffices to create an avalanche of phase-slips, driving the local WL temperature above the bath temperature . When the bias current is ramped down, superconductivity is recovered only at the so-called retrapping current () leading to hysteresis. The hot-spot model by Skocpol, Beasley and Tinkham (SBT)scocpol () considers that in the finite-voltage state the temperature is above the critical temperature . The Josephson coupling is then lost so that no SQUID-type behavior is observed in this state for most of the temperature range.nikhilprl (); hazrair (); hazraprb () Still, it has been observed in some devices based on WLswangjap (); tinkhamnanowire (); wangapl () and SNS junctions.anger (); ronzani (); krasnov () In most cases, the theoretical modeling neglected thermal effects and relied on a conventional RCSJ model, but with an effective capacitancesong () well above the actual geometric one. Eventually, the SBT model was also extended to the case of a WL temperature remaining below , still ignoring the WL phase dynamics.wangjap (); wangapl ()

Recently, some of us proposed a dynamic thermal modelanjanjap () of WLs, incorporating both the overheating of the WL at a temperature and a resistively-shunted junction (RSJ) type phase dynamics. If remains below then the Josephson coupling across the WL is not fully destroyed. The Josephson current, together with the normal current, persists over a portion of the finite voltage branch of the IVCs,twofluidmodel () thus leading to SQUID-type voltage oscillations. This dynamic thermal model (DTM) is in between the SBT model, where the Josephson coupling does not exist at non-zero voltages, and the R(C)SJ model that ignores the thermal effects. A similar approach was used to describe the radio-frequency response of a SNS junction.PRB-DeCecco () Alternatively, hysteresis in the phase slip-controlled regime can be described using a more elaborate non-equilibrium approach using time-dependent Ginzburg-Landau equations.peeter ()

In this paper, we report temperature- and magnetic field-dependent transport in -SQUIDs with a geometry optimized for a moderate Joule heating. Despite being clearly hysteretic, the devices exhibit voltage oscillations with the magnetic flux. The related non-zero flux-sensitive supercurrent contribution surviving well above the critical current, which cannot be understood using the static SBT model, is analyzed with the dynamic thermal model, which quantitatively captures every observation. We eventually discuss the flux-sensitivity of the studied hysteretic -SQUIDs in the dynamic regime.

Ii Dynamic Thermal Model

In the dissipative state, phase slipsnayanaprl (); bezryadin-phase-slip (); bezryadin () occur at a rate , where is the zero-magnetic field critical current, taken here at the bath temperature. Each phase-slip deposits a Joule heat leading to a temperature rise in the WL region. The Joule heat is generated over a length scale determined by the inelastic quasi-particle diffusion length,SBT-phase-slip (); dolan-PS () much longer than the WL dimensions studied here. Thus we assume a uniformnayanaprl () temperature over the entire WL, which determines its critical current . The characteristic time for the thermal balance is where is the WL heat capacity and is the thermal conductance to the bath. Under the quasi-static approximation, the instantaneous WL temperature is dictated by a thermal balance between the Joule heat and the conduction to the bath .

If the temperature remains below at any instant of time , the bias current is dynamically shared between a super-current and a complementary normal current. This gives rise to a time-dependent voltage related to the phase difference through the Josephson relation . The heat balance equation governing the dynamics of temperature is written as . These two equations can be re-arranged in terms of dimensionless variables as:anjanjap ()


Here, the currents and are, respectively, the currents and in units of the zero-magnetic field critical current . The time unit is and . The dimensionless temperature is defined as:


We also define the dimensionless parameter


as the ratio of Joule heat generation (at ) and heat evacuation (at ).

Figure 1: Device state diagram with the grey shaded region as the dynamic regime where the WL has a finite voltage but carries a non-zero supercurrent. The region on the left of this grey-region is the zero-voltage state where all the bias current is carried as supercurrent while the right-region has no supercurrent. The red (green) line depicts the dependence of (). Inset: Equivalent circuit of the DTM. (b) Variation of the dimensionless time-averaged WL temperature with the bias current in the dynamic regime for different values. (c) Ratio of the time-averaged supercurrent and the critical current as a function of at the bias current values of and .

The static retrapping current notation () is defined by the WL being right at the critical temperature (). From Eq. 2, the thermal balance gives:


At a larger current we have so that there is no Josephson coupling and hence . A non-zero supercurrent can be carried by the WL only for .

The dynamic retrapping current is the current below which the dynamic state ceases to exist and the zero-voltage state becomes stable against any phase-slip. Here and in the following, we consider WLs featuring a linear temperature dependence of the critical current and a sinusoidal current-phase relation . One can then obtain by solving Eqs. (1,2):anjanjap ()


From numerical simulations, we find that a non-sinusoidal , within a regime of single-valued current, negligibly affects this relation. Moreover, the elevated WL temperature in the dynamic state gives rise to an increase in coherence length and thus the current-phase relation is close to sinusoidal.

The extent of the dynamic regime, defined by the current bias window , depends on the value of the dimensionless parameter . Figure 1(a) depicts the device state diagram found using the dependence of and . For large values of , i.e. poor heat evacuation and/or high , and are both below and very close to each other. The dynamic regime then occurs in a bias current window of vanishing width, making its observation in IVCs practically impossible.nikhilprl (); nikhilsust () In this limit, the physics is well described by the SBT and other static thermal models.scocpol (); tinkhamnanowire (); nikhilprl () In contrast, for of about unity or smaller, the dynamic regime spans over a significant current range. The static retrapping current can then significantly exceed and the dynamic retrapping current is close to . For extremely good heat evacuation , one has , and RSJ model is recovered.

In the dynamic regime between and , the WL temperature oscillates with time about an average value. However, for a large value of , the magnitude of these oscillations is negligible compared to the average WL temperature.anjanjap () This is always the case as , which can range from tens of ns to s, is greater than , which is of ps order. The WL can thus be considered at a constant (time-averaged) temperature .anjanjap (); PRB-DeCecco (); nayanaprl () The corresponding time-averaged reduced temperature can be calculated as a function of the current bias from Eqs. (1, 2):anjanjap ()


with as the single parameter. Figure 1(b) shows how (or equivalently ) decreases with the current bias for various values of the parameter , starting from 1 (or ) at . At every bias, one can thus calculate the critical current , and the related time-averaged voltage as:


and the time-averaged supercurrent as:


Figure 1(c) shows as a function of the parameter at bias current values equal to and . For , the time-averaged supercurrent is zero when the bias current reaches .scocpol (); nikhilprl () In practice, as soon as exceeds about 2, the WL switches almost immediately to a fully normal state with and (almost) zero supercurrent. For small values, the supercurrent is comparable to the full critical current .

Iii Experimental Details

The fabrication of the -SQUIDs starts with the deposition of a Nb thin film with thickness 40 nm on a Si substrate. A resist layer was afterwards patterned using laser lithography for the outer leads and contact pads, and electron-beam lithography for the smaller structures. A 25-nm thick Al layer was then deposited followed by lift-off. The latter acting as a protective mask, the Nb devices were obtained by a SF reactive-ion etch. The Al mask is eventually chemically etched. Figure 2(b) shows the loop of a SQUID with the two constrictions, of nominal width and length 40 and 160 nm respectively, in parallel. The critical current was tuned down (to the 100 A range) by trimminglahmapl () down the Nb thickness to 20 ( 2) nm in subsequent reactive-ion etching steps without Al mask. The etching process can also lead to a reduction in due to the appearance of a non-superconducting layer on top, bottom and also on the sides.gubin (); siegel (); physicac () However the thickness of such a layercooper (); proximity () is estimated to be 2 nm only, which leaves a large and effective superconducting channel at the core of the WLs.

Electrical transport measurements were carried out in a closed-cycle refrigerator with base temperature 1.3 K. The electrical signals are thoroughly filtered, both at room temperature ( filters) and at base temperature (copper-powder filters). Home-made ground-isolated current sources and voltage amplifiers were used. From the temperature-dependent four-probe transport measurements of a first device, we find the onset of superconductivity at 8.6 K and a sheet resistance = 5.8 in the normal state. In the following, we will present experimental data mainly from one sample. Another sample featured a similar behavior, see Appendix A.

Iv Results

Zero-field current-voltage characteristics (IVCs) of a -SQUID at various bath temperatures are shown in Fig. 2(a). The critical current varies strongly with temperature and exceeds 100 A below 2 K, see Fig. 2(b). The critical current density () at 1.3 K, found as 21.1 MA/cm, is close to the Ginzburg Landau depairing current density = 36 MA/cm (at zero temperature), estimated using parameters for bulk and clean Nb.tinkhambook () This value is similar to our earlier devices.nikhilprl () Together with the linear dependence (up to ) of the critical current with the bath temperature , this confirms the intrinsic superconducting nature of the WLs as opposed to that of SNS WLs.SNS_Ic (); PRB-DeCecco (); nikhilprl () Below the crossover temperature K, hysteresis is observed with a retrapping to the zero-resistance state at a well-defined dynamic re-trapping current anjanjap () , see Fig. 2(a). A fit of its temperature dependence to Eq. 6 provides the value for the critical temperature of the WL = 6.0 K. At bias currents significantly above , further thermal instabilities occur in larger portions of the device as evidenced from additional re-trapping currents, see Appendix A. Our study here is focussed on the hysteretic regime and at bias currents below or in the vicinity of the critical current .

Figure 2: (a) Current-voltage characteristics at different temperatures in the hysteretic regime. Red, black, yellow and blue colors correspond to = 1.67, 2.27, 2.77 and 2.98 K, respectively. (b) Temperature dependence of the critical current and the (dynamic) retrapping current . The crossover temperature is close to 3 K. The solid lines are fits below . The critical current is fit to 211.2, which gives the zero-temperature critical current 211.2 A. The dynamic retrapping is fit to Eq. 6, which gives the WL critical temperature = 6.0 K. Inset: Scanning electron micrograph of the SQUID loop.

In the dissipative branch, the IVC slope varies with the bias current from 5.2 (just above ) to 7 (above ) at 1.67 K, see Fig. 2(a). In the low-bias regime of interest here, the differential resistance is always close to independent of bias current and temperature, see Appendix A. This value is significantly below the value of 11.5 for the resistance of two WLs in parallel estimated from the sheet resistance. This indicates that the WLs are not fully resistive in this dissipative state.

When applying a perpendicular magnetic field , and thus a flux , the critical current displays a -periodic modulation, taking an effective SQUID loop area S = 1.6 m. The flux modulation of the critical current, starting from its maximal zero-field value , is not complete and has a rather triangular shape, see Figs. 3(a,b). This behavior cannot be explained solely by asymmetric critical currents between the two arms. Moreover, self-flux related effects,squidbook () related to the loop inductance , are negligible here, as we estimate 0.1. For a WL with a length , the supercurrent-phase relation is non-sinusoidal.wolf (); hazrananosquid (); iphi1 (); likharev () Numerical calculations using a non-sinusoidal relation indeed yield an incomplete cancellation of at (Appendix B), similar to the experimental behavior.

Figure 3: (a,b) Oscillations of the critical and the dynamic retrapping currents with the magnetic field in the hysteretic regime at 1.3 and 2.42 K, respectively. (c) Voltage oscillation with magnetic field for = 39 to 53 A with 2 A intervals at 2.42 K. (e) The same for = 25 to 50 A at 3.26 K. (d,f) IVCs at different flux values ( and ) at 2.42 and 3.26 K, respectively.

Strikingly, the retrapping current also shows oscillations with the magnetic flux, see Figs. 3(a,b), in contradiction with the SBT picture of a fully normal state of the device in the dissipative state. A similar feature was observed in SQUIDs based on normal metal WLs anger (); ronzani (); krasnov () but not satisfactorily explained. It constitutes a first indication that superconductive coupling is not fully suppressed by the electron heating in the dissipative branch of our hysteretic devices.

Moreover, the SQUID voltage also shows an oscillatory dependence on magnetic field for a wide bias current window, see Figs. 3(c-f). This is completely opposed to the SBT model behavior, in which the dissipative state displays no signature of Josephson coupling across the SQUID WLs. The initial jump in voltage, seen for (see Fig. 3(c)), occurs due to the first arrival to the resistive branch. At a fixed temperature, the IVCs at different flux-values are found to merge on the linear branch, see Fig. 3(d,f), beyond a particular bias current. We identify this current as the static retrapping current ,anjanjap () as discussed in the model section. The oscillations consistently disappear at a bias current beyond . At lower temperatures ( = 1.3 and 2 K) in the hysteretic regime, oscillations are observed over a narrow bias current span just above , see Fig. 4 (a,b).

The flux-to-voltage transduction function , defined as the maximum of , is found to be 27 V/ in the dynamic regime at = 2.42 K. With a voltage noise of 1 nV in our circuit, this gives a flux noise density 37 . In the non-hysteretic regime, thanks to higher flux-to-voltage transduction , the sensitivity reaches 3 at = 3.26 K. The latter value is similar to the ones reported in non-hysteretic -SQUIDs with room temperature amplifiers.granatasquid (); troeman-2007 ()

Figure 4: (a) oscillation at the lowest temperature of 1.3 K at a bias current 48.3 A just above the dynamic retrapping current . (b) Same at 2 K for three different bias currents (45, 46 and 47 A) very close to . At these lower temperatures, voltage values are extracted from IVCs at different magnetic field values.

V Discussion

For modeling the -SQUID, we assume both WLs to be identical, with a temperature-dependent critical current , a temperature-independent normal-state resistance , and a heat-loss coefficient . The two WLs’ phases and maintain a constant difference , forced by the magnetic flux through the SQUID loop. Consequently, the two WLs’ average temperatures are identical in the dynamic steady state. Eventually, the SQUID behaves as a single WL with normal resistance , heat loss coefficient and critical current . The flux-modulation of the critical current alters the expression for the dynamic retrapping current into


At zero flux, this expression matches with Eq. 6. In the limit of a small , one recovers the usual modulation. In contrast, is independent of the flux.

For every bath temperature, we use Eq. 6 of the DTM with the measured zero-flux dynamic retrapping current and critical current to extract the value of . Figure 5(a) shows how the parameter varies with the bath temperature over the hysteretic regime, from close to zero at to about 8 at 1.3 K. For small values, and are (almost) indistinguishable in IVCs. Thus the error bars in increases when is increased towards , see Fig. 5(a) and the method cannot be used beyond. As discussed below, the variation of the residual supercurrent with the bias current is then a more appropriate method to extract the value of . The temperature coefficient of the measured critical currrent below being known (see Fig. 2(b)), we use Eq. 4 with as the single free parameter to fit the curve. We obtain 2.6 nW/K. Alternatively, one can also use the value of obtained from the fit of the dynamic retrapping current as a function of the bath temperature , which gives k 2.8 nW/K. Using a typical value of heat-transfer coefficient as 5 W/cm.K,hazraprb (); nikhilprl () and these two close values of , the effective heat loss area is estimated to be m which is larger than but still comparable to the m area of the WL.

Using the variation of the parameter as a function of the bath temperature , Eq. 7 provides us with the behavior of the WL temperature at a bias current equal to the corresponding dynamic retrapping current , see Fig. 5(a). At the crossover temperature where is small, the WL is at thermal equilibrium with bath, i.e. as in the isothermal RSJ model. Towards low temperature, the WL temperature increases towards .

Figure 5: (a) Black symbols: change in with bath temperature for the measured first device as found using and in the hysteretic regime (below 3 K). The blue solid line is a fit as per Eq. 4 in DTM with = 2.6 nW/K. Red symbols and line represent the WL temperature as a function of the bath temperature and at a bias current equal to the corresponding . (b) Variation of the modulation amplitudes and with bath temperature. (c) Symbols: measured supercurrent modulation amplitude in units of the critical current modulation amplitude in the dissipative state as a function of bias current ranging from to at different bath temperatures. Solid lines: best fit to the DTM with fit parameters listed in Table I. The value of , setting the x-axis scale, depends on the bath temperature. Arrows and dotted lines indicate the positions of the static and dynamic retrapping current respectively.

A small value, required for observation of the dynamic regime, necessitates a small critical current and/or a large thermal conductance to the bath . Compared to earlier similar devices nikhilprl () for which we estimate to be about 20 at 4.2 K, we enhanced by widening the leads right outside the SQUID loop, while still keeping a short and narrow neck between the SQUID loop and the wide leads to avoid vortices between the two WL’s current path.bezryadin () As described earlier, we reduced approximately by an order of magnitude by trimming the Nb thickness down, which overall dominates the decrease in value. The value = 2 obtained at 2.1 K, see Fig. 5(a), approximately defines the low temperature limit for practical operation of the SQUID in the voltage-modulation mode, significantly below the hysteresis temperature . At lower temperatures, the bias-current range of dynamic regime is narrow and the voltage oscillations of small amplitude.

In the experiment, and as discussed above, the critical current is not fully modulated by the flux, which implies the same for the supercurrent. Thus one cannot compare directly the supercurrent calculated from the model to the one deduced from the measured voltage oscillations. We consider the amplitude of the supercurrent modulation in units of the critical current modulation by the flux, i.e. . From the experimental data, we calculate as being equal to where and are the modulation amplitude of, respectively, the voltage and the critical current. As for the theory, we calculate:


with the temperature being found using Eq. 7 with as the single adjustable parameter.

Figure 5(c) shows the experimental (symbols) and theoretical (lines) values of the ratio as a function of the normalized bias current capturing most of the dynamic regime. A very good quantitative agreement is obtained. The fit values of listed in Table 1 agree well with those deduced, and plotted in Fig. 5(a), from the analysis of the dynamic retrapping current. The comparison made here is fully justified only in the case of a sinusoidal . Extending this to the case of a non-sinusoidal current-phase relation is intuitive but still not fully theoretically established. Still, the successful comparison of experimental data and theoretical calculation demonstrates that the DTM describes accurately the transition between the isothermal Josephson junction behavior and the electronically-overheated and hysteretic -SQUID behavior.

T (K) I (A) I (A) from value from Fig. 5 fit
2.42 58 44 1.13 1.1
2.62 46 40 0.66 0.74
2.88 40 38 0.35 0.5
3.26 31.6 31.6 - 0.28
Table 1: Sample parameters including the calculated and fitted values of for different temperatures.

The insights gained from the above study, in particular the key role of the parameter , provide a guideline for designing devices with improved performance. While for a wide dynamic regime is obtained, featuring both hysteretic behavior and SQUID voltage-oscillations, one needs to reach so that hysteresis disappears and the voltage modulations reach a significant fraction of . This is illustrated for a device in Fig. 6, with narrower WLs as compared to the previous one, resulting in a smaller critical current and thus a small , even at 1.3 K.

Figure 6: oscillations of another device with a small critical current = 51 A and thus a small = 0.36 at 1.3 K. Here the bias current ranges from 50 to 70 A.

In this device, the flux-to-voltage transduction function is 1 mV/ at 1.3 K. With an estimated voltage noise of 1 nV in our circuit, we find a flux noise density 1 , i.e. significantly below the values of 3 previously reported in non-hysteretic -SQUIDs using room temperature amplifiers.granatasquid (); troeman-2007 () The use of a low-temperature current amplifier,mag2 () while voltage biasing the -SQUID, is expected to further improve the sensitivity.

Vi Conclusion

In conclusion, we discussed the crossover from the fully overheated WL, i.e. the SBT regime, where the supercurrent is either on or off, to the (isothermal) RSJ case, where the supercurrent contribution decays progressively when the bias current exceeds the critical current. This physics is relevant not only for WLs but also for Josephson junctions based on nanowires, two-dimensional materials and topological insulators, where a large supercurrent density can appear, implying a large power density at the resistive switch, together with poor heat evacuation, thus creating hysteresis. A single parameter reflects the balance between the heat evacuation from the WL and the injected heat, it can be tuned by trimming the critical current and/or varying the thermal coupling to the bath. This balance determines for the amplitude of voltage modulation in the phase dynamic regime. The existence of voltage oscillations makes hysteretic -SQUIDs useful as flux-to-voltage transducers for probing magnetism at the nanoscale with a wide bandwidth.


We are indebted to Thierry Crozes for help in the device fabrication at the Nanofab platform at Néel Institute. Help from Avijit Duley in calculating the current-phase relation is thankfully acknowledged. SB acknowledges financial support from CSIR, Government of India. AKG acknowledges a research grant from the SERB-DST of the Government of India. AKG also thanks Université Grenoble Alpes for a visiting position. CW and HC acknowledge financial support from the LabEx LANEF (ANR-10-LABX-51-01) project ÒUHV-NEQÓ. We acknowledge financial support from CEFIPRA through project 5804-2.

Institute of Engineering Univ. Grenoble Alpes

Appendix A: Additionnal Experimental Data

Figure 7(a,b) shows the large scale SEM image of the -SQUID together with its resistance versus bath temperature () curve and non-hysteretic IVCs. The WL of 6 K cannot be found from R-T plot as the bias current of 10 A is too high.

Figure 7: (a) Larger area SEM image of the -SQUID showing the current and voltage leads. (b) plot at a fixed = 10 A showing the superconducting transition at an onset temperature of 8.6 K. Inset shows some IVCs in the non-hysteretic regime. (c) Hysteretic IVCs over larger range of bias current showing thermal instability at .

Figure 7(c) shows hysteretic I-V characteristics over a large bias range. We see multiple re-trapping currents with the larger magnitude ones representing thermal instabilities in wider portions (see arrows labeled as A and B in Fig. 7(a) of the device as evidenced from the resistance values above the respective re-trapping currents. The physics of thermal instability at higher bias currents (above ), see Fig. 7(c), beyond the dynamic regime has been already reported by some of us.nikhilprl (); nikhilsust ()

Figure 8: The variation of the differential resistance with in the range . The dotted line shows the saturation to K for .

Similar results on transport measurements observed in another (third) -SQUID are shown in Fig. 9(a,b). The crossover temperature of this device is found to be 2.75 K. In the hysteretic regime (below 2.75 K), the retrapping current oscillates along with with (see Fig. 9(a) for = 2.16 K). Estimated resistance from the linear ohmic branch of IVCs beyond is 4.1 . Period of oscillation is same, i.e. 1.25 mT, as that of the first device. The oscillations at 2.16 K are shown in Fig. 9(b) for bias currents ranging from = 25 to 38 A.

Figure 9: (a) oscillation with B at = 2.16 K for a third -SQUID with = 2.75 K. (b) V-B oscillations at this temperature for = 25 to 38 A.

Appendix B: Non-Sinusoidal Relation and Modulation

As opposed to the short () WLs, where the supercurrent-phase relation is sinusoidal (),likharev () the longer WLs () exhibit a non-sinusoidal . Using Ginzburg-Landau theory (which is valid close to ), the relation for a single WL of different lengths is calculated and shown in Fig. 10.wolf () For a SQUID with identical WLs, the magnetic flux gives rise to a phase difference between the two WLs. Adding the two WL’s supercurrents with this phase difference, we obtain the total supercurrent . Figure 10(b) shows the total as a function of at . For a short WL-based SQUID, with perfectly sinusoidal , this total at is identically zero, which is clearly not the case for long WL based SQUIDs.

Figure 10: (a) relation for different values. For , is no longer single valued. (b) as a function of at . Here is SQUID critical current at zero flux. (c) variation with flux. The green, red and blue curves, respectively, represent 1, 2 and 3. The black line in (b) shows for perfectly sinusoidal relation at .

Figure 10(c) shows the oscillations for 1, 2 and 3 with as the maximum (with respect to ) value of at a given . The modulation amplitude decreases with increasing . This demonstrates how the modulation of a SQUID with flux is limited by the non-sinusoidal relation of the WLs.

Appendix C: Rcsj Model

Here we attempt the fitting of the super-current relative modulation with the RCSJ model. This is anyway not very plausible due to the lack of sharp cut-off in which is quite apparent in the experiments.

According to the RCSJ model,squidbook () the current , as shared between resistances, capacitances and the Josephson junctions of the SQUID, can be written as:


Here, and are the resistance and capacitance of each of the two junctions. We have again assumed the screening parameter . using the same dimensionless quantities ate Section I, we get


Here is the effective Stewart-McCumber parameter for the SQUID.

Figure 11: Variation of residual super-current () in the resistive branch at = 2.42, 2.62, 2.88 and 3.26 K as calculated from experimental data (shown by symbols). Dotted lines represent the best fit to the RCSJ model. in x-axis is the critical current value at zero flux for respective .
T (K) I (A) I (A) calculated fitted
2.42 58 44 2.13 5.4
2.62 46 40 1.5 3.2
2.88 40 38 0.8 2.33
3.26 31.6 31.6 - 1.0
Table 2: Comparison of measured , with calculated (from data) and fitted for different temperatures.

We get the IVCs from the numerical steady-state solutions of Eq. 13 and . Using these solutions, we find that depends on . This fact is used to extract the values for experimentally measured . Using the expression , we have tried to fit the measured to the RCSJ model in Fig. 11. The fitted values are listed in Table 2 together with the values extracted from the values of the dynamic retrapping current . We see from the fit that RCSJ does not fit well as compared to the DTM.


  • (1) C. Granata and A. Vettoliere, Nano superconducting quantum interference device: A powerful tool for nanoscale investigations, Phys. Rep. 614, 1-69 (2016).
  • (2) M. J. Martínez-Pérez and D. Koelle, NanoSQUIDs: Basics & recent advances, Phys. Sci. Rev. 2, 20178001 (2017).
  • (3) W. Wernsdorfer, Single nanoparticle measurement techniques, J. Appl. Phys. 87, 5094 (2000).
  • (4) W. Wernsdorfer, Classical and quantum magnetization reversal studied in nanometer-sized particles and clusters, Adv. Chem. Phys. 118, 99 (2001).
  • (5) M. J. Martínez-Perez, B. Muller, D. Schwebius, D. Korinski, R. Kleiner, J. Sese, D. Koelle, NanoSQUID magnetometry of individual cobalt nanoparticles grown by focused electron beam induced deposition, Supercond. Sci. and Technol. 30, 024003 (2017).
  • (6) C. Veauvy, K. Hasselbach and D. Mailly, Scanning -superconduction quantum interference device force microscope, Rev. Sci. Instrum, 73, 3825 (2002).
  • (7) A. Finkler, Y. Segev , Y. Myasoedov, M. L. Rappaport, L. NeÕeman, D. Vasyukov, E. Zeldov, M. E. Huber, J. Martin and A. Yacoby, Self-Aligned Nanoscale SQUID on a Tip, Nano. Lett 10, 1046 (2010).
  • (8) D. Vasyukov, Y. Anahory, L. Embon, D. Halbertal, J. Cuppens, L. Neeman, A. Finkler, Y. Segev, Y. Myasoedov, M. L. Rappaport, M. E. Huber and E. Zeldov, A scanning superconducting quantum interference device with single electron spin sensitivity, Nature Nanotech. 8, 639 (2013).
  • (9) L. Chen, H. Wang, X. Liu, L. Wu and Z. Wang, A High-Performance Nb Nano-Superconducting Quantum Interference Device with a Three-Dimensional Structure, Nano. Lett 16, 7726 (2016).
  • (10) J. E. Duvauchelle, A. Francheteau, C. Marcenat, F. Chiodi, D. Débarre, K. Hasselbach, J. R. Kirtley, and F. Lefloch, Silicon superconducting quantum interference device, Appl. Phys. Lett. 107, 072601 (2015).
  • (11) R. Wolbing, J. Nagel, T. Schwarz, O. Kieler, T. Weimann, J. Kohlmann, A. Zorin, M. Kemmler, R. Kleiner, D. Koelle, Nb nano superconducting quantum interference devices with high spin sensitivity for operation in magnetic fields up to 0.5 T, Appl. Phys. Lett. 102, 192601 (2013).
  • (12) S. Mandal, T. Bautze, O. A. Williams, C. Naud, E. Bustarret, F. Omnès, P. Rodière, T. Meunier, C. Bauerle, L. Saminadayar, The Diamond SQUID, ACS Nano. 5, 7144 (2011).
  • (13) T. Schwarz, J. Nagel, R. Wölbing, M. Kemmler, R. Kleiner, and D. Koelle, Low-Noise Nano Superconducting Quantum Interference Device Operating in Tesla Magnetic Fields, ACS Nano. 7, 844 (2013).
  • (14) R. Russo, E. Esposito, A. Crescitelli, E. Di Gennaro, C. Granata, A. Vettoliere, R. Cristiano and M. Lisitskiy, NanoSQUIDs based on niobium nitride films, Supercond. Sci. and Technol. 30, 024009 (2017).
  • (15) A. Ronzani, C. Altimiras and F. Giazotto, Highly Sensitive Superconducting Quantum-Interference Proximity Transistor, Phys. Rev. Appl. 2, 024005 (2014).
  • (16) J. Clarke and A. I. Braginski (Editors), The SQUID Handbook, 2004.
  • (17) M. Tinkham, Introduction to Superconductivity 2nd ed. (Mc. Graw-Hill, New York, 1996).
  • (18) W. J. Skocpol, M. R. Beasley and M. Tinkham, Self-heating hotspots in superconducting thin-film microbridges, J. Appl. Phys. 45, 4054 (1974).
  • (19) M. Tinkham, J. U. Free, C. N. Lau and N. Markovic, Hysteretic I-V curves of superconducting nanowires, Phys. Rev. B 68, 134515 (2003).
  • (20) H. Courtois, M. Meschke, J. T. Peltonen and J. P. Pekola, Origin of Hysteresis in a Proximity Josephson Junction, Phys. Rev. Lett. 101, 067002 (2008).
  • (21) D. Hazra, L. M. A. Pascal, H. Courtois and A. K. Gupta, Hysteresis in superconducting short weak links and -SQUIDs, Phys. Rev. B 82, 184530 (2010).
  • (22) V. V. Baranov, A. G. Balanov and V. V. Kabanov, Dynamics of resistive state in thin superconducting channels, Phys. Rev. B 87, 174516 (2013).
  • (23) N. Kumar, C. B. Winkelmann, S. Biswas, H. Courtois and A. K. Gupta, Controlling hysteresis in superconducting constrictions with resistive shunt, Supercond. Sci. and Technol. 28, 072003 (2015).
  • (24) S. K. H. Lahm and D. L. Tilbrook, Development of a niobium nanosuperconducting quantum interference device for the detection of small spin populations, Appl. Phys. Lett. 82, 1078 (2003).
  • (25) T. Schwarz, R. Wölbing, C. F. Reiche, B. Müller, M. J. Martínez-Pérez, T. Mühl, B. Büchner, R. Kleiner and D. Koelle, Low-Noise YBaCuO Nano-SQUIDs for Performing Magnetization-Reversal Measurements on Magnetic Nanoparticles, Phys. Rev. Appl. 3, 044011 (2015).
  • (26) R. Arpaia, M. Arzeo, S. Nawaz, S. Charpentier, F. Lombardi and T. Bauch, Ultra low noise YBaCuO nano superconducting quantum interference devices implementing nanowires, Appl. Phys. Lett. 104, 072603 (2014).
  • (27) G. Kimmel, A. Glatz, and I. S. Aranson, Phase slips in superconducting weak links, Phys. Rev. B 95, 014518 (2017).
  • (28) N. Shah, D. Pekker and P. M. Goldbart, Inherent Stochasticity of Superconductor-Resistor Switching Behavior in Nanowires, Phys. Rev. Lett. 101, 207001 (2008).
  • (29) N. Kumar, T. Fournier, H. Courtois, C. B. Winkelmann and A. K. Gupta, Reversibility of superconducting Nb weak links driven by the proximity effect in a quantum interference device, Phys. Rev. Lett. 114, 157003 (2015).
  • (30) D. Hazra, J. R. Kirtley and K. Hasselbach, Retrapping Current in Bridge-Type Nano-SQUIDs, Phys. Rev. Appl. 4, 024021 (2015).
  • (31) L.-K. Wang, D.-J. Hyun, and B. S. Deaver, Heating and flux flow in niobium variable-thickness bridges, J. Appl. Phys. 49, 5602 (1978).
  • (32) L.-K. Wang, A. Callegari, B. S. Deaver, D. W. Barr, and R. J. Mattauch Microwave mixing with niobium variable thickness bridges, Appl. Phys. Lett. 31, 306 (1977).
  • (33) V. M. Krasnov, T. Golod, T. Bauch, and P. Delsing, Anticorrelation between temperature and fluctuations of the switching current in moderately damped Josephson junctions, Phys. Rev. B 76, 224517 (2007).
  • (34) L. Angers, F. Chiodi, G. Montambaux, M. Ferrier, S. Guéron, H. Bouchiat, and J. C. Cuevas, Proximity dc squids in the long-junction limit, Phys. Rev. B 77, 165408 (2008).
  • (35) A. Ronzani, M. Baillergeau, C. Altimiras, and F. Giazotto, Micro-superconducting quantum interference devices based on V/Cu/V Josephson nanojunctions, Appl. Phys. Lett. 103, 052603 (2013).
  • (36) Y. Song, Origin of capacitance in superconducting microbridges, J. Appl. Phys. 47, 2651 (1976).
  • (37) A. K. Gupta, N. Kumar and S. Biswas, Temperature and phase dynamics in superconducting weak links, J. Appl. Phys. 116, 173901 (2014).
  • (38) J. Bardeen, Two-Fluid Model of Superconductivity, Phys. Rev. Lett. 1, 399 (1958).
  • (39) A. De Cecco, K. Le Calvez, B. Sacépé, C. B. Winkelmann, and H. Courtois, Interplay between electron overheating and ac Josephson effect, Phys. Rev. B 93, 180505(R) (2016).
  • (40) G. Berdiyorov, K. Harrabi, F. Oktasendra, K. Gasmi, A. I. Mansour, J. P. Maneval, and F. M. Peeters, Dynamics of current-driven phase-slip centers in superconducting strips, Phys. Rev. B 90, 054506 (2014).
  • (41) S. L. Chu, A. T. Bollinger, and A. Bezryadin, Phase slips in superconducting films with constrictions, Phys. Rev. B 70, 214506 (2004).
  • (42) M. W. Brenner, D. Roy, N. Shah, and A. Bezryadin, Dynamics of superconducting nanowires shunted with an external resistor, Phys. Rev. B 85, 224507 (2012).
  • (43) W. J. Skocpol, M. R. Beasley and M. Tinkham, Phae-slip centers and nonequilibrium processes in tin microbridges, J. Low Temp. Phys. 16, 145 (1974).
  • (44) G. J. Dolan and L. D. Jackel, Voltage measurements with nonequilibrium region near phase-slip centers, Phys. Rev. Lett. 39, 1628 (1977).
  • (45) Note that in the original paper introducing the DTM,anjanjap () we used a different notation, , for this current.
  • (46) A. I. Gubin, K. S. Il’in, S. A. Vitusevich, M. Siegel, and N. Klein, Dependence of magnetic penetration depth on the thickness of superconducting Nb thin films, Phys. Rev. B 72, 064503 (2005).
  • (47) I. Charaev, T. Silbernagel, B. Bachowsky, A. Kuzmin, S. Doerner, K. Il’in, A. Semenov, D. Roditchev, D. Yu. Vodolazov and M. Siegel, Proximity effect model of ultranarrow NbN strips, Phys. Rev. B 96, 184517 (2017).
  • (48) K. Il’in, D. Rall, M. Siegel, A. Engel, A. Schilling, A. Semenov, H.-W. Huebers, Influence of thickness, width and temperature on critical current density of Nb thin film structures, Physica C 470, 953 (2010).
  • (49) L. N. Cooper, Superconductivity in the Neighborhood of Metallic Contacts, Phys. Rev. Lett. 6, 689 (1961).
  • (50) Ya. V. Fominov and M. V. Feigel’man, Superconductive properties of thin dirty superconductor-normal-metal bilayers, Phys. Rev. B 63, 094518 (2001).
  • (51) P. Dubos, H. Courtois, B. Pannetier, F. K. Wilhelm, A. D. Zaikin, and G. Schön, Josephson critical current in a long mesoscopic S-N-S junction, Phys. Rev. B 63, 064502 (2001).
  • (52) E. de Wolf and R. de Bruyn Ouboter, The Josephson supercurrent-phase relation, Physica B 176, 133 (1992).
  • (53) D. Hazra, J. R. Kirtley and K. Hasselbach, Nano-superconducting quantum interference devices with suspended junctions, Appl. Phys. Lett. 104, 152603 (2014).
  • (54) K. Hasselbach, D. Mailly and J. R. Kirtley, Micro-superconducting quantum interference device characteristics, J. Appl. Phys. 91, 4432 (2002).
  • (55) K. K. Likharev, Superconducting weak links, Rev. Mod. Phys. 51, 101 (1979).
  • (56) A. G. P. Troeman, H. Derking, B. Boerger, J. Pleikies, D. Veldhuis and H. Hilgenkamp, NanoSQUIDs Based on Niobium Constrictions, Nano Lett. 7, 2152 (2007).
  • (57) D. Pekker, A. Bezryadin, D. S. Hopkins and P. M. Goldbart, Operation of a superconducting nanowire quantum interference device with mesoscopic leads, Phys. Rev. B 72, 104517 (2005).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description