Zurek-Kibble Mechanism for the Spontaneous Vortex Formation in NbAl/Alox/Nb Josephson Tunnel Junctions: New Theory and Experiment

  • R. Monaco, J. Mygind, M. Aaroe, R. J. Rivers, and V. P. Koshelets, Istituto di Cibernetica del C.N.R., 80078, Pozzuoli, Italy and Unitá INFM-Dipartimento di Fisica, Universitá di Salerno, 84081 Baronissi, Italy, Department of Physics, B309, Technical University of Denmark, DK-2800 Lyngby, Denmark, Blackett Laboratory, Imperial College London, London SW7 2AZ, United Kingdom, and Institute of Radio Engineering and Electronics, Russian Academy of Science, Mokhovaya 11, Building 7, 125009, Moscow, Russia
Phys. Rev. Lett. 96, 180604
Log-log plot of the measured frequency f1 of trapping single fluxons versus the quenching time τQ. Each point corresponds to many thermal cycles. The vertical error bars gives the statistical error while the relative error bars in τQ amounting to ±10% are as large as the dots’ width. The solid line is the fit to an allometric relationship f1=aτQb which yields a=0.0092±10% (taking τQ in seconds) and b=0.51±6%.Log-log plot of the measured frequency f1 of trapping single fluxons versus the quenching time τQ. Each point corresponds to many thermal cycles. The vertical error bars gives the statistical error while the relative error bars in τQ amounting to ±10... Show more

Because phase transitions take place in a finite time, causality guarantees that correlation lengths remain finite, even when the transitions are continuous. The Zurek-Kibble (ZK) scenario [1–3] proposes that such transitions take effect as fast as possible, i.e., the domain structure after the quenching of the system initially reflects the causal horizons. As a result, once the transition is implemented, correlation lengths scale as (positive) powers of the quench time τQ (inverse quench rate). This proposal is amenable to direct testing for transitions whose domain boundaries carry topological charge. Such is the case for Josephson tunnel junctions (JTJs), where the topological charge is a fluxon, i.e., a supercurrent vortex carrying magnetic flux Φ0=h/(2e) in the plane of the oxide layer between the two superconductors that make up the JTJ. In the case considered in this Letter of annular JTJs, obtained by the superposition of two superconducting rings, it corresponds to a magnetic flux that threads just one of the two rings an odd number of times.

In 2000 an idealized model was proposed [4] to test the ZK scenario for JTJs, assuming that causal horizons are constrained by the velocity of electromagnetic waves in the JTJ, the Swihart velocity [5,6]. As a result, the probability f1 for spontaneously producing one fluxon in the thermal quench of a symmetric annular JTJ of circumference C was predicted to scale as

In Eq. (1), ξ¯ is the Zurek-Kibble causal length, the correlation length of the relative phase angle at the time of defect formation. It is defined through Eq. (1) in terms of the cold correlation length ξ0, the relaxation time τ0 of the long wavelength modes, and the quench time τQ, in turn defined by: TC/τQ=-(dT/dt)T=TC. Equation (1) holds for C<ξ¯. Under the assumptions of (a) weak coupling of the superconductors and (b) exact critical slowing of the Swihart velocity at the critical temperature T=Tc, we predicted σ=0.25 [4].

In 2001 our first proof-of-principle experiment [7] was performed, to test Eq. (1). The experiment consisted of taking an annular JTJ isolated from its surroundings and making it undergo a forced phase transition by heating it above its superconducting critical temperature and letting it cool passively back towards the liquid-He temperature. Once the thermal cycle is over, the junction I-V curve is inspected and any trapped fluxon can be detected by the appearance of current peaks at discrete voltages in the I-V characteristic of a JTJ. By detecting the voltage position of these peaks the number of trapped fluxons (or antifluxons) can be determined.

The experiment of Ref. [7] was remarkably successful, as the only experiment of the several [8–17] previously performed on condensed matter systems sensitive enough to show quantitative ZK scaling behavior, although qualitative support for scaling had been demonstrated with He3 and high- Tc superconductors. However, this preliminary experiment suffered from severe restrictions, in particular, a limited range of cooling rates, statistically poor data, and above all insufficient shielding of the earth’s magnetic field with the possibility of systematic errors. The outcome was σ commensurate with 0.25.

In this Letter we shall present results from a new experiment designed to circumvent these problems, that we shall discuss later. Our experiment shows extremely reliable scaling behavior of the form (1), but with σ=0.5 to high accuracy. This is obviously at total variance with our earlier prediction. However, as we shall show, the value σ=0.25 depended critically on both assumptions (a) and (b) cited earlier being exactly satisfied. For realistic JTJs these assumptions are only approximately satisfied in the immediate vicinity of T=Tc. Since defects form very close to Tc, we shall see that this is sufficient for the value of σ to jump from 0.25 to σ=0.5. Thus, rather than its negation, our experiment arguably provides an even more robust demonstration of the ZK scenario.

To reiterate, the theory in Ref. [4] was developed for JTJs whose electrodes are weak coupling superconductors; for such JTJs the temperature dependence of the critical current density Jc(T) is given by the Ambegaokar-Baratoff equation [18]:

where Δ(T) is the superconducting gap energy and ρN is JTJ normal resistance per unit area squared. Equation (2) provides a linear decrease of Jc near Tc:

where T(t=0)=Tc and the dimensionless quantity α is approximately equal to 2Δ(0)/kBTC=3.5.

However, our JTJs are based on Nb, a strong-coupling superconductor, for which Eq. (2) is not necessarily valid. In practice, high quality and reproducible barriers are achieved by depositing a thin Al overlay onto the Nb base electrode which will be only partially oxidized, leaving a Nb-Al bilayer underneath having a non-BCS temperature dependence of the energy gap and of the density of states. The proximity effect in Nb-Al/Alox/Nb JTJs has been extensively studied and it is known to influence the electrical properties of the junctions, such as the current-voltage characteristic and the temperature dependence of the critical current density. Specifically, the proximity effect in superconductor–normal-metal–insulator–superconductor junctions can lead to dominance by an otherwise subdominant temperature dependence of the critical current density [19] in the vicinity of Tc of the form

where α is a constant depending intricately on the degree of proximity. The last equation models the tail shaped dependence of Jc vs T near Tc; it has been theoretically derived and experimentally confirmed by Golubov et al. [20] in 1995. Rephrasing the arguments of Ref. [4] with Eq. (4) replacing Eq. (3), the Josephson penetration depth λJ(T), which plays the role of system equilibrium coherence length ξ(T), near Tc diverges linearly with time


ds being electrode thickness.

ξ˙(t)=dξ(t)/dt<0 measures the rate at which the defects contract, i.e., the speed of interfaces between ordered and disordered ground states. Since ξ˙(t) decreases with time t>0, the earliest possible time t at which defects could possibly appear is determined by causality,

where t¯ is the causal time and c¯ is the Swihart velocity.

As we said, in Ref. [4] we had also assumed that the Swihart velocity vanished at Tc. For realistic JTJs, this is not so; it just becomes very small. Swihart [5] has demonstrated that for a thin-film superconducting strip transmission line the solution for the velocity varies continuously as one passes through the critical temperature into the normal state. As a result, we assume c¯(t)=c¯nn near the transition temperature where c¯nn is the speed of light in a microstrip line made of normal metals. In the case of a microstrip line made by two electrodes having the same thickness ds and the same skin depth δ, with dsδ, separated by a dielectric layer of thickness dox and dielectric constant ϵ, c¯nn(2/δ)doxds/ϵμ [5]. Its value depends on the temperature very weakly, but depends on the frequency f through δ=ρ/πμf, ρ being the normal-metal residual resistivity.

The match c¯(Tc)=c¯nn is certainly realistic and we still have approximate critical slowing down insofar as c¯nn is much smaller than the zero temperature Swihart velocity c¯0dox/2λL0ϵμ, i.e., when the zero temperature London penetration depth λL0δ2/ds. For 300 nm thick Nb electrodes ( ρ=3.8μΩcm and λL0=90nm), δ1mm at say f=10kHz, so the last inequality is fully satisfied. At the same frequency, for a value of the specific barrier capacitance cs=ϵ/dox=0.027F/m2 typical of low current density Nb-Al/Alox/Nb JTJs, we get c¯nn=6×103m/s and c¯0=1.6×107m/s [21].

The solution of the causality equation Eq. (7) with a nonvanishing Swihart velocity yields: t¯=ξ0τQ/c¯nn, again discretely different from its idealized counterpart [4]. Inserting the value of t¯ into Eq. (5) we obtain the new Zurek length ξ¯

where τ0=ξ0/c¯nn ( τ0=O(1ns). We reach the important conclusion for realistic JTJs that the probability f1 for spontaneously producing one fluxon in the quench is still predicted to scale with the quench time τQ according to Eq. (1), but the critical exponent is now σ=0.5, rather than σ=0.25. Detailed calculations will be reported elsewhere [22].

However, it is worth observing that by varying τQ in the experimentally achieved four-decade range 1ms<τQ<10s, we get 10μs<t¯<1ms that is always much larger than τ0; it means that by the time the Josephson phase freezes the Josephson effect is well established. Further, in the same τQ range the normalized freezing temperature T¯/TC at which the defect is formed is 0.99<T¯/TC<0.9999. It would be really hard, if not impossible, to measure the temperature dependence of Jc and c¯ so close to TC. We have then to resort to theoretical predictions.

The new experiment has a faster and more reliable single heating system, obtained by integrating a meander line 50μm wide, 200 nm thick, and 8.3 mm long Mo resistive film in either ends of the 4.2mm×3mm×0.35mm Si chip containing the Nb/AlOx/Nb trilayer JTJs. These resistive elements have a nominal resistance of 50Ω at LHe temperatures and, due to their good adhesion with the substrate, are very effective in dissipating heat in the chip. In fact, voltage pulses a few μs long and a few volts high applied across the integrated heater provided quench times as low as 1 ms, that is more than 2 orders of magnitude smaller than for the previous situation [7]. Further, automatization of thermal cycles was implemented that allowed for much more robust statistics to be achieved. At the end of each thermal quench the junction I-V curve is automatically stored and an algorithm has been developed for the detection of the trapped fluxons. Finally, all the measurements have been carried out in a magnetic and electromagnetically shielded environment. The low frequency magnetic shielding was achieved by using μ-metal, cryoperm, and lead cans. During the quench the JTJ was also electrically isolated.

The experimental results reported here are restricted to one of two annular JTJs with similar geometry to the sample used in the earlier experiment ( C=500μm and ds=300nm), but with about 50 times lower critical current density ( Jc(0)60A/cm2): this leads to a value of ξ0=17μm (we have assumed α=α=3.5 in Eq. (6)). An extensive description of the chip layout, the experimental setup and the quenching data will be given elsewhere [22] while the details of the fabrication process can be found in Ref. [23]. The quench time τQ was continuously varied over more that 4 orders of magnitude (from 20 s down to 1 ms) by varying the helium exchange gas pressure inside the vacuum can and/or the width and the amplitude of the voltage pulse across the integrated resistive element.

Figure 1 shows on a log-log plot the measured frequency f1=n1/N of single fluxon trapping, obtained by quenching the sample N times for each value of a given quenching time τQ, n1 being the number of times that the inspection of the low temperature JTJ current-voltage characteristics at the very end of each thermal cycle showed that one defect was spontaneously produced. N ranged between 100 and 2600 and n1 was never smaller then 10, except for the rightmost point for which n1=3 (and N=1800). The sample has undergone a total of more than 100 000 thermal cycles without any measurable change of its electrical parameters. The vertical error bars gives the statistical error f1/n1. The measurement of τQ by fitting solutions to the heat equation follows that of [7] in its high accuracy. The relative error bars in τQ amounting to ±10% are as large as the dot’s width. To test Eq. (1), we have fitted the data with an allometric function f1=aτQ-b, with a and b as free fitting parameters. A linear regression of logf1 vs logτQ, represented by the continuous line in Fig. 1, yields b=0.51±6%, in excellent agreement with the predicted value 0.5. The same fit gives a=0.0092±10% (taking τQ in seconds) that is 6–7 times larger than the predicted value C/ξ0c¯nn (with C=500μm, ξ0=17μm and c¯nn=6×103m/s). As a bound we only expect agreement in the overall normalization a to somewhat better than an order of magnitude. Empirically, the different condensed matter experiments have shown that the ratio aobserved/apredicted varies widely from system to system; O(1) for superfluid He3 [8,9], very small for high- Tc superconductors [13]. We point out that in our case the value of the prefactor a is dependent, although weakly, on the choice of α and f, being C/ξ0c¯nn(α/f)1/4. The choice f=10kHz was determined assuming that 1/a212kHz is the characteristic frequency of our system at the time of the thermal quench. As we noted earlier, the new scaling exponent σ=0.5, obtained by applying causality arguments to realistic JTJs, is twice as large than that observed earlier [7], also in samples made with the same Nb-Al/Alox/Nb technology. The reason for this discrepancy resides in the fact that at that time we were unaware of the high sensitivity of f1 to external magnetic fields which, although small, were most likely present. In the present experiment much care has been devoted both to avoid magnetic or current carrying materials in the cryoprobe during the quench and to screen the chip environment from 50 Hz noise as well as from dc magnetic field.

The data of Fig. 1 resolve another issue. It has been disputed that the Swihart velocity, which gives the behavior above, is the relevant velocity for field ordering along the JTJ oxide. If we had taken the relevant velocity to be that of phase ordering in the individual superconductors, as invoked by Zurek [1,2] when considering the spontaneous flux produced on quenching annuli of simple superconductors, at the same level of approximation we would have predicted b=0.25 and the prefactor orders of magnitude larger.

We have not observed the change in behavior predicting f1 to change from the linear behavior with C/ξ¯ of Eq. (9) to a random walk behavior in the phase,

once f1 is sufficiently large. A first guess would suggest that the transition from Eq. (1) to Eq. (9) would occur when f11. However, future experiments to be carried out on ring-shaped JTJs having the circumference larger than that of the sample reported in this paper and with τQ in the same range of this experiment should clearly reveal the transition from Eq. (1) to Eq. (9).

In summary, we see this experiment as providing strong corroboration of Zurek-Kibble scaling over a wide range of quenching time τQ in accord with our predictions for Nb-Al/Alox/Nb JTJs. As such, it replaces the experiment reported in Ref. [7] by being more realistic theoretically and more sophisticated experimentally. We stress that this experiment is the only one to date to have confirmed the Zurek-Kibble scaling exponent for a condensed matter system. Further experiments can be devised to investigate the transition to the random walk regime and the effect of the thermal gradients. Such experiments are in the process of being performed.


  1. W. H. Zurek, Nature (London) 317, 505 (1985); ; Acta Phys. Pol. B24, 1301 (1993).
  2. W. H. Zurek, Phys. Rep. 276, 177 (1996).
  3. T. W. B. Kibble, Phys. Rep. 67, 183 (1980).
  4. E. Kavoussanaki, R. Monaco, and R. J. Rivers, Phys. Rev. Lett. 85, 3452 (2000); ; R. Monaco, R. J. Rivers, and E. Kavoussanaki, J. Low Temp. Phys. 124, 85 (2001).
  5. J. C. Swihart, J. Appl. Phys. 32, 461 (1961).
  6. A. Barone and G. Paternó, Physics and Applications of the Josephson Effect (John Wiley & Sons, New York, 1982).
  7. R. Monaco, J. Mygind, and R. J. Rivers, Phys. Rev. Lett. 89, 080603 (2002); ; Phys. Rev. B 67, 104506 (2003).
  8. C. Bauerle et al., Nature (London) 382, 332 (1996).
  9. V. M. H. Ruutu et al., Nature (London) 382, 334 (1996).
  10. P. C. Hendry et al., Nature (London) 368, 315 (1994).
  11. M. E. Dodd et al., Phys. Rev. Lett. 81, 3703 (1998); ; J. Low Temp. Phys. 115, 89 (1999).
  12. R. Carmi and E. Polturak, Phys. Rev. B 60, 7595 (1999).
  13. A. Maniv, E. Polturak, and G. Koren, Phys. Rev. Lett. 91, 197001 (2003).
  14. R. Carmi, E. Polturak, and G. Koren, Phys. Rev. Lett. 84, 4966 (2000).
  15. J. R. Kirtley, C. C. Tsuei, and F. Tafuri, Phys. Rev. Lett. 90, 257001 (2003).
  16. S. Casado, W. González-Viñas, H. Mancini, and S. Boccaletti, Phys. Rev. E 63, 057301 (2001).
  17. S. Ducci, P. L. Ramazza, W. González-Viñas, and F. T. Arecchi, Phys. Rev. Lett. 83, 5210 (1999).
  18. V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 486 (1963); ; 11, 104(E) (1963).
  19. N. L. Rowell and H. J. T. Smith, Can. J. Phys. 54, 223 (1976).
  20. A. A. Golubov et al., Phys. Rev. B 51, 1073 (1995).
  21. This value comes from the Swihart formula c¯0=dox/2λL0ϵμ, with a +20% correction due to the presence of the idle region.

  22. R. Monaco, J. Mygind, M. Aaroe, V. P. Koshelets, and R. J. Rivers (to be published).
  23. P. N. Dmitriev, I. L. Lapitskaya, L. V. Filippenko, A. B. Ermakov, S. V. Shitov, G. V. Prokopenko, S. A. Kovtonyuk, and V. P. Koshelets, IEEE Trans. Appl. Supercond. 13, 107 (2003).

About the Authors

Image of R. Monaco
Image of J. Mygind
Image of M. Aaroe
Image of R. J. Rivers
Image of V. P. Koshelets

Related Articles

Focus: Better Microwaves from a Metamaterial

Focus: Better Microwaves from a Metamaterial

A system that produces radiation by sending electrons through a metamaterial waveguide could potentially generate high-power microwaves with high efficiency. Read More »

Arts &amp; Culture: Science and Filmmaking in an Equal Relationship

Arts & Culture: Science and Filmmaking in an Equal Relationship

Imagine Science—a yearly film festival that takes place in New York, Paris, and Abu Dhabi—is forging new ways to bring together scientists and filmmakers. Read More »

Synopsis: Explaining Aftershock Clustering

Synopsis: Explaining Aftershock Clustering

A study of bursting phenomena like earthquakes suggests that events appear to cluster in time because of the way that small events like aftershocks are identified. Read More »

More Articles