Experimental evidence of a \varphi Josephson junction

Experimental evidence of a Josephson junction

H. Sickinger Physikalisches Institut and Center for Collective Quantum Phenomena in LISA, Universität Tübingen, Auf der Morgenstelle 14, D-72076 Tübingen, Germany    A. Lipman The Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel    M. Weides Peter Grünberg Institute and JARA-Fundamentals of Future Information Technology, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany Current address: Physikalisches Institut, Karlsruher Institut für Technologie, 76131 Karlsruhe, Germany    R. G. Mints The Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel    H. Kohlstedt Nanoelektronik, Technische Fakultät, Christian-Albrechts-Universität zu Kiel, D-24143 Kiel, Germany    D. Koelle    R. Kleiner    E. Goldobin Physikalisches Institut and Center for Collective Quantum Phenomena in LISA, Universität Tübingen, Auf der Morgenstelle 14, D-72076 Tübingen, Germany
July 14, 2019

We demonstrate experimentally the existence of Josephson junctions having a doubly degenerate ground state with an average Josephson phase . The value of can be chosen by design in the interval . The junctions used in our experiments are fabricated as 0- Josephson junctions of moderate normalized length with asymmetric 0 and regions. We show that (a) these Josephson junctions have two critical currents, corresponding to the escape of the phase from and states; (b) the phase can be set to a particular state by tuning an external magnetic field or (c) by using a proper bias current sweep sequence. The experimental observations are in agreement with previous theoretical predictions.

Josephson junction
74.50.+r, 85.25.Cp

Josephson junctions (JJs) with a phase shift of in the ground stateBulaevskiĭ et al. (1977) attracted a lot of interest in recent yearsBaselmans et al. (1999); Ryazanov et al. (2001a); Kontos et al. (2002); Weides et al. (2006a); van Dam et al. (2006); Gumann et al. (2007). In particular, these JJs can be used as on-chip phase batteries for biasing various classicalOrtlepp et al. (2006) and quantumFeofanov et al. (2010) circuits, allowing for removing external bias lines and reducing decoherence. Currently, it is possible to fabricate simultaneously both 0 and JJs using various technologies such as superconductor-ferromagnet heterostructuresRyazanov et al. (2001b); Weides et al. (2006b); Frolov et al. (2008); Weides et al. (2007); Pfeiffer et al. (2008) or JJs based on d-wave superconductorsVan Harlingen (1995); Tsuei and Kirtley (2000); Smilde et al. (2002); Gürlich et al. (2009).

It would be remarkable to have a phase battery providing an arbitrary phase shift , rather than just 0 or . The simplest idea is to combine 0 and JJs to obtain a JJ. However, this is not as straightforward as it may seem. The balance between 0 and parts is complicated as shown in the pioneering workBulaevskii et al. (1978), where conditions for having a non-trivial -state were derived. Long artificialBuzdin and Koshelev (2003); Pugach et al. (2010) and naturalMints (1998); Il’ichev et al. (1999); Mints et al. (2002) arrays of -0--0-- JJs with short segments were analyzed in detail and suggested as systems, where a JJ could be realized. More recently, only one period of such an array, i.e., one 0- JJ, was analyzed in an external magnetic fieldGoldobin et al. (2011). In these works, the Josephson phase is considered as a sum of a constant (or slowly varying) phase and a deviation from the average phase. Then, for the average phase one obtains an effective current-phase relation (CPR) for the supercurrentGoldobin et al. (2011)


where the averaged value of the critical current . The CPR (1) exactly corresponds to a JJ at zero normalized magnetic field , if , cf. Fig.1(a). Here is the total length of JJ, while and are the lengths of 0 and parts, accordingly, is the width of JJ.

Figure 1: (Color online) (a) effective CPR and (b) effective Josephson energy calculated numerically (thick lines) in comparison to those given by approximate analytical formulas (1) and (3) (thin lines). For the curves look mirror reflected with respect to the -axis, i.e., . (c) numerically calculated effective Josephson energy at tilted by an applied bias current (energy supplied by the current source ). Several important situations are shown: ground state (same curve as in (b)), (retrapping), (escape from well), and (escape from well).

It is worth noting that the term “ JJ”, introduced in Ref. Buzdin and Koshelev, 2003, refers to a JJ with a degenerate ground state phase . In the particular case of Eq. (1) at one has . The coefficients and are defined asLipman et al. (2012)


where , are the lengths normalized to the Josephson length and , are the critical current densities of 0 and parts, respectively. Here is calculated using the average value of the critical current .

The physics of JJs with a CPR given by Eq. (1) is quite unusualGoldobin et al. (2007). In particular, one should observe two critical currentsGoldobin et al. (2007, 2011) at , corresponding to the escape of the phase from the left () or the right () well of the double-well Josephson energy potential


where . The critical currents are different because the maximum slope (maximum supercurrent in Fig. 1(a)) on the rhs (positive bias) of the well is smaller than the maximum slope (maximum supercurrent in Fig. 1(a)) on the rhs of the well, see Fig. 1(b).

In this letter we present experimental evidences of a JJ made of one 0 and one segment, see Fig. 2(a).

Figure 2: (Color online) (a) Sketch (cross section) of the investigated SIFS 0- JJ with a step in the F-layer thickness. (b) Optical image (top view, colored manually) of an investigated sample having overlap geometry. The junction area (0-segment and -segment) is . The “etched F-layer” area shows where the thickness of the F-layer was reduced from to to produce a 0- JJ.
Figure 3: (Color online) Domain of existence of state. The shows the position of the investigated JJ at .

The samples were fabricated as NbAl-AlONiCuNb heterostructuresWeides et al. (2006b, 2010). These superconductor-insulator-ferromagnet-superconductor (SIFS) JJs have an overlap geometry as shown in Fig. 2(b). Each junction consists of two parts, a conventional 0-segment and a -segment. It is well knownKontos et al. (2002); Oboznov et al. (2006); Weides et al. (2006a); Buzdin (2005) that the critical current in SFS or SIFS JJs strongly depends on the thickness of the F-layer and can become negative within some range of ( junction). Therefore, to produce the 0 and the segments, the F-layer has different thicknesses and , as shown in Fig. 2(a). To achieve this, the F-layer of thickness , corresponding to ( JJ) was fabricated first. Then the area indicated in Fig. 2(b) was etched down to , corresponding to (0 JJ). Usually one obtains asymetric valuesKemmler et al. (2010), i.e. .

We have studied 3 samples. One of them has very little asymmetry of about 5%, and , which corresponds to . In this case even 5% of asymmetry brings the sample out of the domain, see Fig. 3 for a qualitative picture. The other two JJs have and are deep inside the -domain in parameter space. Both samples show similar results. Here we present the results obtained on one of them. Its parameters are summarized in Tab. 1 and its position within the domain is indicated in Fig. 3. The values of and cannot be measured directly. In our case we have measured the dependence (see below) and then simulated it numerically using and as fitting parameters. The best fitting was obtained for the values specified in Tab. 1.

parameter -part -part whole JJ
physical length
, , ()
normalized lengths 0.68 0.68 1.36
Table 1: Junction parameters at . The values of , and the normalized length are obtained from the fits.

According to theoretical predictionsGoldobin et al. (2007, 2011) for a JJ at zero magnetic field one expects two critical currents (for each bias direction), corresponding to the escape of the phase from and wells of , respectively for bias current and vice versa for . The current is always observed. To observe one has to have low damping so that retrapping in the well is avoided (for positive ). The characteristic (IVC) of the investigated 0- JJ at measured at show only one critical current. Therefore the experiments were performed in a He cryostat at down to where the damping in SIFS JJs reduces drasticallyPfeiffer et al. (2008). In the temperature range from down to both and are clearly visible in the IVCs as shown in Fig. 4.

EarlierGoldobin et al. (2007) we proposed a technique that allows to choose which critical current one traces in the IVC, i.e. from which well the phase escapes. The control is done by choosing a proper bias sweep sequence. For example, if the junction is returning from the positive voltage state, the potential is tilted so that the phase slides to the right. When the tilt becomes small enough the phase will be trapped, presumably in the right well, cf. Fig. 1(c). However, this natural assumption is not always true (see below). Then, if the phase is trapped at , we can sweep the bias (a) in the positive direction and will observe escape from (to the right) at or (b) in the negative direction to observe escape from (to the left) at .

In experiment, at , when the damping is very low, the currents and are traced in random order. Recording one IVC after the other, we were able to obtain IVCs with all 4 possible combinations: (a) , (b) , (c) , (d) . Choosing a specific sweep sequence as described above does not make the outcome ( or ) predictable. We believe that in this temperature range the damping is so low that, upon returning from the positive voltage state, the phase does not simply stop in the well, but can also reflect from the barrier and find itself in a well, cf. Fig. 1(c). The absence of determinism suggests that most probably we are dealing with a system exhibiting chaotic dynamics. This issue will be investigated elsewhere. Nevertheless, at we managed to achieve deterministic behavior as described above, see Fig. 4.

Figure 4: (Color online) Current-voltage characteristics showing lower and higher critical currents measured at . At this temperature the behavior is deterministic: if one sweeps as one always observes the critical currents ; if one sweeps as , one always observes ; finally, if one sweeps as , one always observes .

Another fingerprint of a JJ is its dependence, which (a) should have the main cusp-like minima shifted off from point-symmetrically with respect to the origin, see Fig. 4 of Ref. Goldobin et al., 2011; (b) should show up to four branches in total (for both sweep polarities) at low magnetic fieldGoldobin et al. (2011). In essence, the latter feature results from the escape of the phase from two different energy minima in two different directions and is an extension of the two-critical-currents story to the case of non-zero magnetic field. Instead, the feature (a) alone cannot serve as a proof of a JJ as it is a common feature of every asymmetric 0- JJ even if its ground state is 0 or Goldobin et al. (2011).

The experimentally obtained dependence at is shown in Fig. 5. First, in the whole temperature range we observe the main minima shifted point-symmetrically to a finite field, which is . Second, at one observes the crossing of two branches and two critical currents for each bias current polarity. The left hand branches correspond to the escape of the phase out of the -well, while the right branch corresponds to the escape from the -well. Note, that at the crossing of branches is not visible. One just observes a cusp-like minimum. In this case, both cannot be seen together and the existence of two states, as should be present for a -JJ, cannot be proven.

To trace the intersection of the branches better we were applying a special value of the magnetic field during reset from the voltage state back to state, to “prepare” a specific state ( or ) of the system. Then the field was set back to the “current” value and the bias was ramped up to trace . By doing this for different values of , one is able to trace the intersection of the branches much better than without this technique, see Fig. 5. Still, in experiment we were not able to trace the branches for positive and negative up to the point where they meet, as shown in Fig. 4(a) and (b) of Ref. Goldobin et al., 2011, probably because of retrapping.

Figure 5: (Color online) Experimentally measured dependence (black symbols) and numerically calculated curve (red/gray smaller symbols) that provides the best fit.

To extract some parameters of our JJ from the curve, we performed numerical simulations, by solving the full sine-Gordon equation for a 0- JJ, using the normalized length and the critical currents and as fitting parameters. The simulations were performed using StkJJGoldobin (2011). The objective was to obtain the best fit close to the origin, especially the point where the branches of cross. The best fit is shown in Fig. 5 and was obtained for and the values and from Tab. 1. Measured data in Fig. 5 are shifted along axis by to compensate for average remanent magnetization of the F-layer. To obtain an almost perfect fit we have also assumed a difference in the constant remanent magnetizations in the 0 and in the parts of , similar to earlier studiesKemmler et al. (2010); Weides et al. (2010). One can also see that one of the experimental branches after the main maximum runs parallel to the simulated curve (the experimental is not point symmetric). This shift stays the same even if we cycle the sample through of Nb, but changes if we cycle through the Curie temperature of the F-layer. We conclude that this is related to the non-uniform remanent magnetization of the F-layer.

The parameters obtained from the fit allow to see the location of our JJ on the plane, see Fig. 3. One sees that the JJ is not really short and lays quite deep in the domain. In this region the analytic resultsGoldobin et al. (2011); Lipman et al. (2012) and, in particular, Eqs. (1) and (3) are valid only qualitatively, cf., Fig. 3 in Ref. Goldobin et al., 2011. Therefore we calculated the CPR of our JJ numerically. We started with a 0- JJ with parameters , , and obtained from the fit, assumed some applied bias current and found all static solutions numerically. Then for each of these solutions we calculated and plotted all those on a plot. By repeating this procedure for different , we obtain an effective CPR shown in Fig. 1(a) together with the curve produced by the analytical formula (1). One can see that the exact effective CPR calculated numerically and the approximate CPR calculated analytically are qualitatively similar. From the numerical CPR the value of for our particular JJ is extracted.

In conclusion, we have demonstrated experimentally the realization of a Josephson junction based on a - SIFS junction. We have observed experimentally a shift of the minimum according to the effective CPR (1), as predictedGoldobin et al. (2011), as well as two critical currents close to zero magnetic field for each bias current polarity. We showed that one can choose between the states of the system using an externally applied magnetic field, which removes their degeneracy. We also showed that one can bring the system into one of the two states by properly tilting the potential using the bias current. Depending on the damping (temperature) one can achieve deterministic behavior (when damping is small enough to see the lower critical current, but large enough to trap the phase in a particular well) as well as random behavior at very low damping. The obtained JJ has at ).

We acknowledge financial support by the German Israeli Foundation (Grant No. G-967-126.14/2007) and DFG (via SFB/TRR-21, project A5 as well via project Go-1106/3). H.S. gratefully acknowledges support by the Evangelisches Studienwerk e.V. Villigst


  • Bulaevskiĭ et al. (1977) L. N. Bulaevskiĭ, V. V. Kuziĭ,  and A. A. Sobyanin, JETP Lett. 25, 290 (1977), [Pis’ma Zh. Eksp. Teor. Fiz. 25, 314 (1977)].
  • Baselmans et al. (1999) J. J. A. Baselmans, A. F. Morpurgo, B. J. V. Wees,  and T. M. Klapwijk, Nature (London) 397, 43 (1999).
  • Ryazanov et al. (2001a) V. V. Ryazanov, V. A. Oboznov, A. Y. Rusanov, A. V. Veretennikov, A. A. Golubov,  and J. Aarts, Phys. Rev. Lett. 86, 2427 (2001a).
  • Kontos et al. (2002) T. Kontos, M. Aprili, J. Lesueur, F. Genêt, B. Stephanidis,  and R. Boursier, Phys. Rev. Lett. 89, 137007 (2002).
  • Weides et al. (2006a) M. Weides, M. Kemmler, E. Goldobin, D. Koelle, R. Kleiner, H. Kohlstedt,  and A. Buzdin, Appl. Phys. Lett. 89, 122511 (2006a)cond-mat/0604097 .
  • van Dam et al. (2006) J. A. van Dam, Y. V. Nazarov, E. P. A. M. Bakkers, S. De Franceschi,  and L. P. Kouwenhoven, Nature (London) 442, 667 (2006).
  • Gumann et al. (2007) A. Gumann, C. Iniotakis,  and N. Schopohl, Appl. Phys. Lett. 91, 192502 (2007).
  • Ortlepp et al. (2006) T. Ortlepp, Ariando, O. Mielke, C. J. M. Verwijs, K. F. K. Foo, H. Rogalla, F. H. Uhlmann,  and H. Hilgenkamp, Science 312, 1495 (2006).
  • Feofanov et al. (2010) A. K. Feofanov, V. A. Oboznov, V. V. Bol’ginov, J. Lisenfeld, S. Poletto, V. V. Ryazanov, A. N. Rossolenko, M. Khabipov, D. Balashov, A. B. Zorin, P. N. Dmitriev, V. P. Koshelets,  and A. V. Ustinov, Nat. Phys. 6, 593 (2010).
  • Ryazanov et al. (2001b) V. V. Ryazanov, V. A. Oboznov, A. V. Veretennikov,  and A. Y. Rusanov, Phys. Rev. B 65, 020501 (2001b).
  • Weides et al. (2006b) M. Weides, M. Kemmler, H. Kohlstedt, R. Waser, D. Koelle, R. Kleiner,  and E. Goldobin, Phys. Rev. Lett. 97, 247001 (2006b)cond-mat/0605656 .
  • Frolov et al. (2008) S. M. Frolov, M. J. A. Stoutimore, T. A. Crane, D. J. Van Harlingen, V. A. Oboznov, V. V. Ryazanov, A. Ruosi, C. Granata,  and M. Russo, Nat. Phys. 4, 32 (2008).
  • Weides et al. (2007) M. Weides, C. Schindler,  and H. Kohlstedt, J. Appl. Phys. 101, 063902 (2007).
  • Pfeiffer et al. (2008) J. Pfeiffer, M. Kemmler, D. Koelle, R. Kleiner, E. Goldobin, M. Weides, A. K. Feofanov, J. Lisenfeld,  and A. V. Ustinov, Phys. Rev. B 77, 214506 (2008)0801.3229 .
  • Van Harlingen (1995) D. J. Van Harlingen, Rev. Mod. Phys. 67, 515 (1995).
  • Tsuei and Kirtley (2000) C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 (2000).
  • Smilde et al. (2002) H.-J. H. Smilde, Ariando, D. H. A. Blank, G. J. Gerritsma, H. Hilgenkamp,  and H. Rogalla, Phys. Rev. Lett. 88, 057004 (2002).
  • Gürlich et al. (2009) C. Gürlich, E. Goldobin, R. Straub, D. Doenitz, Ariando, H.-J. H. Smilde, H. Hilgenkamp, R. Kleiner,  and D. Koelle, Phys. Rev. Lett. 103, 067011 (2009).
  • Bulaevskii et al. (1978) L. N. Bulaevskii, V. V. Kuzii,  and A. A. Sobyanin, Solid State Commun. 25, 1053 (1978).
  • Buzdin and Koshelev (2003) A. Buzdin and A. E. Koshelev, Phys. Rev. B 67, 220504(R) (2003)cond-mat/0305142 .
  • Pugach et al. (2010) N. G. Pugach, E. Goldobin, R. Kleiner,  and D. Koelle, Phys. Rev. B 81, 104513 (2010).
  • Mints (1998) R. G. Mints, Phys. Rev. B 57, R3221 (1998).
  • Il’ichev et al. (1999) E. Il’ichev, V. Zakosarenko, R. P. J. IJsselsteijn, H. E. Hoenig, H.-G. Meyer, M. V. Fistul,  and P. Müller, Phys. Rev. B 59, 11502 (1999).
  • Mints et al. (2002) R. G. Mints, I. Papiashvili, J. R. Kirtley, H. Hilgenkamp, G. Hammerl,  and J. Mannhart, Phys. Rev. Lett. 89, 067004 (2002).
  • Goldobin et al. (2011) E. Goldobin, D. Koelle, R. Kleiner,  and R. G. Mints, Phys. Rev. Lett. 107, 227001 (2011)1110.2326 .
  • Lipman et al. (2012) A. Lipman, E. Goldobin, D. Koelle, R. Kleiner,  and R. Mints, “Josephson junction with a magnetic-field tunable current-phase relation,”  (2012), unpublished.
  • Goldobin et al. (2007) E. Goldobin, D. Koelle, R. Kleiner,  and A. Buzdin, Phys. Rev. B 76, 224523 (2007)0708.2624 .
  • Weides et al. (2010) M. Weides, U. Peralagu, H. Kohlstedt, J. Pfeiffer, M. Kemmler, C. Gürlich, E. Goldobin, D. Koelle,  and R. Kleiner, Supercond. Sci. Technol. 23, 095007 (2010).
  • Oboznov et al. (2006) V. A. Oboznov, V. V. Bol’ginov, A. K. Feofanov, V. V. Ryazanov,  and A. I. Buzdin, Phys. Rev. Lett. 96, 197003 (2006).
  • Buzdin (2005) A. I. Buzdin, Rev. Mod. Phys. 77, 935 (2005).
  • Kemmler et al. (2010) M. Kemmler, M. Weides, M. Weiler, M. Opel, S. T. B. Goennenwein, A. S. Vasenko, A. A. Golubov, H. Kohlstedt, D. Koelle, R. Kleiner,  and E. Goldobin, Phys. Rev. B 81, 054522 (2010).
  • Goldobin (2011) E. Goldobin, ‘‘StkJJ – User’s Reference,” http://www.geocities.com/SiliconValley/Heights/7318/StkJJ.htm (2011).
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