InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Engineering » Electrical and Electronic Engineering » "Superconductors - Materials, Properties and Applications", book edited by Alexander Gabovich, ISBN 978-953-51-0794-1, Published: October 17, 2012 under CC BY 3.0 license. © The Author(s).

Effective Models of Superconducting Quantum Interference Devices

By R. De Luca
DOI: 10.5772/48483

Article top

Overview

Figure 1. Schematic representation of a two-junction quantum interferometer

Figure 2. Average phase difference for cos(πψex)=0.3 and: a)iB=1.6; b)iB=0.4; c)iB=0.6. Dotted blue lines represent ϕ0(τ) for β = 0.02, full red lines represent ϕ(τ) as calculated to first order for β = 0.02, and the dashed black lines represents the numerical solution of the complete system. In a) the first order approximation of the solution is not shown, for clarity reasons.

Figure 3. First order correction to ϕ0(τ) as calculated by the procedure described in the text by setting β=0.02, cos(πψex)=0.3and: a)iB=1.6; b)iB=0.4; c)iB=0.6.

Figure 4. Circulating current iS as a function of the normalized time for null values of the initially trapped flux and for iB=2.2 and: a) ψex=0.1(orange), ψex=0.3(cyan) ; b)ψex=0.5; c) ψex=0.7(cyan), ψex=0.9(orange).

Figure 5. Circulating current iS in the presence of an oscillating applied magnetic field shown as a function of the normalized time for A=0, B=0.1, iB=2.5, β=0.01 and for normalized frequencies equal to: a)ω˜=0.03; b)ω˜=0.06; c)ω˜=0.09. No initially trapped flux in the superconducting loop is present.

Figure 6. a) Time average value 〈ic〉 of the critical current ic as a function of the amplitude B of the oscillating magnetic flux for A=0.1 (blue) and A=0.5 (orange). (b) Time average value 〈ic〉 of the critical current ic as a function of A for B=0.1 (blue) and B=0.5 (orange). Both figures are obtained by taking null values of the initially trapped flux.

Figure 7. Schematic representation of a one-dimensional array of equal Josephson junctions.

Figure 8. Critical current iC=|BN+1| as a function of the applied magnetic flux for an array of 10 (a) and 15 (b) Josephson junctions.

Figure 9. a) I-V characteristics of an array of 10 Josephson junctions for the following values of the applied magnetic fluxΨex: 0.05 (blue); 0.075 (cyan); 0.10 (orange); 0.15 (red). (b) I-V characteristics of an array of 15 Josephson junctions for the following values of the applied magnetic fluxΨex: 0.025 (blue); 0.05 (cyan); 0.075 (orange); 0.10 (red).

Figure 10. Average voltage vs. applied magnetic flux for an array of 10 (a) and 15 (b) Josephson junctions for three different values of the normalized bias currentiB. In (a) iB=4.0, 10.0, 16.0(from bottom to top). In (b) iB=4.0, 10.0, 15.0(from bottom to top).

Figure 11. Schematic representation of a Josephson junction array with alternate parameters. In the branches of the parallel connections overdamped 0-junctions (green circles) and π-junctions (red diamonds) are alternately present.

Figure 12. a) Envelope function N E(Ψex) as a function of2 F(Ψex). (b) Critical current as a function of the normalized applied flux for a homogeneous parallel array with N=12, ε=1, σ=1.

Figure 13. Critical current as a function of Ψex for a non-homogeneous parallel array with N=12, ε=1.3, σ=1.5 in the following ranges: (a) [-1.25, 1.25]; (b) [-2.5, 2.5]. In (a) the envelope function N E(Ψex) is shown as a dotted line.

Figure 14. Critical current as a function of Ψex for a non-homogeneous parallel array with N=12, ε=1.3, σ=2in the following ranges: (a) [-1.75, 1.75]; (b) [-5.5, 5.5]. In (a) the envelope function N E(Ψex) is shown as a dotted line.

Figure 15. Critical current ic as a function of the oscillating amplitude B. (a) The normalized frequency is fixed at ω˜=1.0 and various values of the d. c. component A are represented: A=0.1 (orange); A=0.2 (brown); A=0.3 (cyan). (b) The d. c. component A is fixed at A=0.1 and the normalized frequency attains the following values: ω˜=0.5(orange); ω˜=1.0(brown); ω˜=1.5(cyan).

Figure 16. Critical current ic as a function of the normalized frequency ω˜ for A=0.1 and for various values of the oscillating amplitude B: B=0.4 (orange); B=0.8 (brown); B=1.2 (cyan).

Figure 17. Voltage v versus the d. c. component A of the applied flux for iB=2.5 and B=0.0 (orange); B=0.15 (brown); B=0.30 (cyan). In (a) the normalized frequency values isω˜=0, in (b)ω˜=1.0.

Figure 18. Voltage v versus the d. c. component A of the applied flux for iB=2.5 and ω˜=0.5 (orange); ω˜=1.0(brown); ω˜=1.5(cyan). In (a) the a. c. component of the applied magnetic flux is B=0.50, in (b) B=1.0.

Effective Models of Superconducting Quantum Interference Devices

R. De Luca

1. Introduction

It is well known that the electrodynamic properties of SQUIDs (Superconducting Quantum Interference Devices) are obtained by means of the dynamics of the Josephson junctions in these superconducting system (Barone & Paternò, 1982; Likharev, 1986; Clarke & Braginsky, 2004). Due to the intrinsic macroscopic coherence of superconductors, r. f. SQUIDs have been proposed as basic units (qubits) in quantum computing (Bocko et al., 1997). In the realm of quantum computing non-dissipative quantum systems with small (or null) inductance parameter and finite capacitance of the Josephson junctions (JJs) are usually considered (Crankshaw & Orlando, 2001). The mesoscopic non-simply connected classical devices, on the other hand, are generally operated and studied in the overdamped limit with negligible capacitance of the JJs and small (or null) values of the inductance parameter. Nonetheless, r. f. SQUIDs find application in a large variety of fields, from biomedicine to aircraft maintenance (Clarke & Braginsky, 2004), justifying actual scientific interest in them.

As for d. c. SQUIDs, these systems can be analytically described by means of a single junction model (Romeo & De Luca, 2004). The elementary version of the single-junction model for a d. c. SQUID takes the inductance L of a single branch of the device to be negligible, so that β = LIJ0 ≈0, where Φ0 is the elementary flux quantum and IJ is the average value of the maximum Josephson currents of the junctions. In this way, the Josephson junction dynamics is described by means of a nonlinear first-order ordinary differential equation (ODE) written in terms of the phase variable φ, which represents the average of the two gauge-invariant superconducting phase differences, φ1 and φ2, across the junctions in the d. c. SQUID. By considering a device with equal Josephsons junction in each of the two symmetric branches, the dynamical equation of the variable φ can be written as follows (Barone & Paternò, 1982):

 dϕdτ+(−1)ncosπψexsinϕ=iB2 (1)

where n is an integer denoting the number of fluxons initially trapped in the superconducting interference loop, τ=2πRIJt0=t/τϕ, R being the intrinsic resistive junction parameter, ψexex0 is the externally applied flux normalized to Φ0 and iB= IB/IJ, is the bias current normalized to IJ. In what follows we shall consider zero-field cooling conditions, thus taking n=0. Eq. (1) is similar to the non-linear first-order ODE describing the dynamics of the gauge-invariant superconducting phase difference across a single overdamped JJ with field-modulated maximum current IJF (IJF=|cosπψex|) in which a normalized bias current iB/2 flows. This strict equivalence comes from the hypothesis that the total normalized flux ψ=Φ/Φ0 linked to the interferometer loop can be taken to be equal to ψex. However, being

 ψ=ψex+β(i1−i2), (2)

we may say that the above hypothesis may be stated merely by means of the following identity: β=0. Therefore, for finite values of the parameter β, Eq. (1) is not anymore valid and the device behaves as if the equivalent Josephson junction possessed a non-conventional current-phase relation (CPR). In fact, for small finite values of β, one can see that the following model may be adopted (Romeo & De Luca, 2004):

 dϕdτ+Xexsinϕ+πβYex2sin2ϕ=iB2 (3)

where Xex= cosπψex and Yex= sinπψex. A second-order harmonic in ϕ thus appears in addition to the usual sinϕ term. The sin2ϕ addendum, however, arises solely from electromagnetic coupling between the externally applied flux and the system, as described by Eq. (2), when β0. Therefore, the non-conventional CPR of the equivalent JJ in the SQUID model cannot be considered as a strict consequence of an intrinsic non-conventional CPR of the single JJs. The Josephson junctions in the device, in fact, could behave in the most classical way, obeying strictly to the Josephson current-phase relation; the interferometer, however, would still show the additional sin2ϕ term for finite values of β. In order to understand how the reduction in the dimensional order of the dynamical equations is possible, it is noted that the quantities τϕ = Φ0/2πRIJ and τψ = L/R, denoting the characteristic time scales of the variables ϕ and of the number of fluxons ψ in the superconducting SQUID loop, respectively, are intimately linked to the parameter β, since τψϕ =β. In this way, for constant applied magnetic fields, the flux dynamics for small values of β can be considered very fast with respect to the equivalent junction dynamics given in Eq. (3). As a consequence, the superconducting phase ϕ can be assumed to be adiabatic and the equation of motion for ψ in terms of the quasi-static variable ϕ can be solved by perturbation analysis. When the information for ψ is substituted back into the effective dynamical equation for ϕ, Eq. (3) is finally obtained.

The single-junction model can be adopted also when dealing with more complex systems, as one-dimensional Josephson junction arrays. In fact, by assuming small inductance values in the N current loops of a one-dimensional array containing N+1 identical overdamped Josephson junctions, the dynamical equations for the gauge-invariant superconducting phase differences can be reduced to a single non linear differential equation (Romeo & De Luca, 2005). The resulting time-evolution equation is seen to be similar to the single-junction dynamical equation with an appropriately defined current-phase relation. As specified before, the critical current, the I-V characteristics and the flux-voltage curves of the array can be determined analytically by means of the effective model. Furthermore, a one dimensional array of N cells of 0- and π-junctions in parallel can be considered (De Luca, 2011). In this case, by assuming that junctions parameters and effective loop areas alternate as one moves along the longitudinal direction of the array, going from 0- to π-junctions, an effective single junction model for the system can be derived. It can be shown that, by this model, interference patterns of the critical current as a function of the applied magnetic flux can be analytically found and compared with existing experiments (Scharinger et al., 2010).

Finally, a single-junction model for a d.c. SQUID is derived when we consider the effect of rapidly varying applied fields whose frequency ω is comparable withτψ1. By letting the applied magnetic flux have, in addition to a constant term A, an a. c. component, we can find, by similar reasoning as in the case of a constant applied magnetic flux, the effective reduced single-junction model for the system. In particular, for β = 0, the critical current of the device is seen to depend on A, and on the frequency ω and the amplitude B of the a. c. component of the applied magnetic flux in a closed analytic form. From the analysis of the voltage vs. applied flux curves it can be argued that the quantities ω and B can play the role of additional control parameters in the device.

The work will thus be organized as follows. In Section 2 the derivation of an effective single-junction model for a symmetric d. c. SQUID containing two equal junctions will be briefly reviewed. In Section 3 the extension to this model to Josephson junction arrays with equal junctions in all branches will be considered. In Section 4 the case of the alternate presence of 0- and π-junctions in the array is considered, the system being similar to multifacets Josepshon junctions. In Section 5 the effective single-junction model for a d. c. SQUID in the presence of rapidly varying field is derived. Finally, in Section 6 conclusions are drawn and further investigations are suggested.

2. Two-junction quantum interference devices

Let us consider a symmetric two junction interferometer with equal junctions of negligible capacitance, as shown in fig. 1. The dynamical equations for the variables ϕ and ψ characterizing this system, can be written in the following form (Romeo & De Luca, 2004):

 dϕdτ+(−1)ncos(πψ)sinϕ=iB2;                 (a)πdψdτ+(−1)nsin(πψ)cosϕ+ψ2β=ψex2β.  (b) (4)

where n is the number of initially trapped fluxons in the superconducting loop. Let us consider a new time variable θ=τ2πβ=RLt and write the solutions for ϕ a and ψ in the following form:ϕ(β,τ)ϕ0,β(τ)+βϕ1,β(τ);ψ(β,τ)ψ0(θ)+βψ1(θ). For simplicity, set n=0.

Figure 1.

Schematic representation of a two-junction quantum interferometer

This approach allows us not only to account for the regular part of the solution, as seen in (Grønbech-Jensen et al., 2003) and in (Romeo & De Luca, 2004), but also to consider its singular part. Moreover, as we shall see, the role of the two time variables will become evident in what follows, since one time scale is defined for Eq. (4a) and one for Eq. (4b). Consider then ϕ(β,τ) and ψ(β,τ) to be bounded, differentiable functions, and expand the sine and cosine functions appearing in Eqs. (4a-b) to first order inβ. By then collecting all coefficients of identical power of β, we can obtain a system of equations for the functions ϕk,β(τ) andψk(θ), with k = 0, 1, describing the k-th order solutions for ϕ and ψ, respectively. These approximate solutions are determined according to the following sequential scheme. As a first step, we use Eq. (4b) to determineψ0(θ). We adopt the solution found and substitute it in Eq. (4a) to determineϕ0,β(τ). The latter solution, on its turn, is substituted in Eq. (4b) to find ψ1(θ) and, finally, this solution is used in Eq. (4a) to findϕ1,β(τ). Note, therefore, that for defining first order solutions, knowledge of zero-th order solutions is required. Furthermore, we assume that the initial conditions are the following:

 ϕ(β,0)=ϕ0,β(0)+βϕ1,β(0); ψ(β,0)=ψ0(0)+βψ1(0) (5)
.

As for initial conditions, from Eq. (4b) we may notice that ψ0(τ)=ψex forβ=0, in which case we cannot even define the time variable θ. This condition, however, is inherited by the functionψ0(θ), since the following equalities are satisfied:

 ψ0(τβ)β=0=ψ0(θ)θ→∞=ψex. (6)

Furthermore, we may also notice thatψk(θ)θ=0=ψk(τβ)τ=0=ψk(0), for k = 0, 1.

By the general procedure described above we get the following differential equations for the superconducting phase variables:

 dϕ0,βdτ+cos(πψ0(θ))sinϕ0,β(τ)=iB2,                                                            (a)dϕ1,βdτ+ϕ1,β(τ)cos(πψ0(θ))cosϕ0,β(τ)=πψ1(θ)sin(πψ0(θ))sinϕ0,β(τ); (b) (7)

and the following for the flux number variables:

 dψ0dθ+ψ0(θ)=ψex,                                            (a)dψ1dθ+ψ1(θ)=−2sin(πψ0(θ))cosϕ0,β(2πβθ). (b) (8)

In Eqs. (7a-b) and (8a-b) we may notice the appearance of two different time scales the first, τψ=LR, linked to flux motion in and out the superconducting ring, the second, τϕ=Φ02πRIJ, pertaining to the dynamics of the superconducting phase difference value ϕ. We have already noticed thatτψτϕ=2πβ, so that, for negligible values of this ratio, the system behaves effectively as if the adiabatic time evolution of the superconducting phase difference variable ϕ could be studied by taking asymptotic solutions of ψ. In this case, therefore, we may first let the flux variable evolve, so that a stationary magnetic state is reached, and then solve for the superconducting phase difference time evolution of the system. This is exactly what is done, under the assumptions of negligible value of the ratior=τψτϕ, in (Grønbech-Jensen et al., 2003) and in (Romeo & De Luca, 2004). However, if one were to acquire the regular solution for the system dynamics, even when considering the approximate solution for the variables ϕ and ψ, one would follow the more general perturbation analysis described above, where the ratio r might not a priori be considered as negligible. Furthermore, considering that this ratio is proportional to the perturbation parameter β, one might wish to generalize the procedure described above to higher order in β to acquire a wider range of validity of the analysis.

Despite the fact that the more general approach allows extension to higher order approximations of the perturbation solutions, we wish to limit our analysis to the study of the electrodynamic properties of a two-junction or a multi-junction quantum interferometer with very small parameter β. Therefore, while in the present section we shall only be concerned with a single time scale, namelyτϕ, by assuming that the transient of the flux variable rapidly vanishes (τψτϕ<<1), in Section 5 we shall consider both evolution time, by introducing a rapidly varying external magnetic field.

By considering, for the time being, only the time scaleτϕ, the dynamical equations for the flux variable (Eqs. (8a-b)) give the following steady-state solution for ψ0 andψ1:

 ψ0=ψex,                                             (a)ψ1=−2sin(πψex)cosϕ0−2πdψexdτ, (b) (9)

where the term 2πdψexdτ has been inserted in Eq. (9b), in order to correctly take account of first order contributions in β, and the subscript β in ϕ0,β(τ) has been elided, as it will be done for ϕ1,β(τ) from this point on, since these functions will not depend on β in this limit.

In order to obtain some preliminary results, we start by considering ψex as constant. By the general procedure schematized above we get the following differential equations for the superconducting phase variables:

 dϕ0dτ+cos(πψex)sinϕ0(τ)=iB2,                                                   (a)dϕ1dτ+ϕ1(τ)cos(πψex)cosϕ0(τ)=πψ1(τ)sin(πψex)sinϕ0(τ),  (b) (10)

and the following for the flux number variables:

 ψ0=ψex,                                     (a)ψ1(τ)=−2sin(πψex)cosϕ0(τ). (b) (11)

According to the scheme described above, by having already setψ0=ψex, we may now solve for ϕ0(τ) in Eq. (10a). Let us therefore briefly discuss how to obtain this solution. In the case|iB2cos(πψex)|=1|a|>1, which characterizes the running state of the junctions, we have

 ϕ0(τ)=2tan−1[a+1−a2tan(γτ+tan−1ξ0)].  +2kπ (12)

where k is an integer, γ=12iB24cos2(πψex)andξ0=tan(ϕ0(0)2)a1a2.

On the other hand, in the case|iB2cos(πψex)|=1|a|<1, which characterizes the superconducting state of the junctions, we have

 ϕ0(τ)=2tan−1[a+a2−1tanh(γ˜τ+tanh−1(χ0))], (13)

where γ˜=12cos2(πψex)iB24 andχ0=tan(ϕ0(0)2)aa21.

Finally, in the case|a|=1, we have

 ϕ0(τ)=sgn(a)[2tan−1(iBτ2+ω0(±))−π2], (14)

whereω0(±)=tan(sgn(a)ϕ0(0)2+π4). Having found the time dependence of the variableϕ0, ψ1can be found by Eq. (11b) by substitution. Finally, by knowledge ofψ0, ϕ0andψ1, the function ϕ1 can be found by Eq. (10b), which is a standard first order linear differential equation. Solutions for ϕ0 are shown for cos(πψex)=0.3 in figs. 2a, 2b, and 2c foriB=1.6, iB=0.4, andiB=0.6, respectively, along with the solution obtained by numerically integrating Eqs. (4a-b). In figs. 2b-c the first order approximation of the solution is shown as a dotted line.

Solutions for ϕ1 are shown in figs. 3a-c for cos(πψex)=0.3 andiB=1.6, iB=0.4andiB=0.6, respectively. The above analysis thus leads to a solution in a closed form, to first order in the parameter β. Notice that in the case of time-dependent bias currents one should adopt a more general procedure.

As a simple application, let us calculate, to first order in the parameter β, the circulating current iS in the circuit, normalized toIJ, given by (Barone & Paternò, 1982):

 iS=ψ−ψexβ. (15)

For an arbitrary value n, which represents the number of fluxons initially trapped in the superconducting ring, we have

 iS=ψ1(τ)=−2sin(πψex)cosϕ0(τ). (16)

Graphs of circulating currents are shown in figs. 4a, 4b, 4c for n even, iB=2.2and forψex=0.1and0.3, ψex=0.5, ψex=0.7and0.9, respectively. The period T of these curves is equal to the pseudo-period of ϕ0,β which is given by the following expression in terms of ψex andiB:

 T=2π(iB2)2−cos2(πψex). (17)

Notice that the lowest value of the period is obtained for ψex=0.5 and that the curves for ψex=0.1and0.9 and forψex=0.3and0.7, although having the same period, as it can be argued from Eq. (17), are not equal.

Figure 2.

Average phase difference for cos(πψex)=0.3 and: a)iB=1.6; b)iB=0.4; c)iB=0.6. Dotted blue lines represent ϕ0(τ) for β = 0.02, full red lines represent ϕ(τ) as calculated to first order for β = 0.02, and the dashed black lines represents the numerical solution of the complete system. In a) the first order approximation of the solution is not shown, for clarity reasons.

Figure 3.

First order correction to ϕ0(τ) as calculated by the procedure described in the text by setting β=0.02, cos(πψex)=0.3and: a)iB=1.6; b)iB=0.4; c)iB=0.6.

Figure 4.

Circulating current iS as a function of the normalized time for null values of the initially trapped flux and for iB=2.2 and: a) ψex=0.1(orange), ψex=0.3(cyan) ; b)ψex=0.5; c) ψex=0.7(cyan), ψex=0.9(orange).

The above results have been obtained for the magnetic response of the system in the presence of a constant applied flux. In the following we shall analyze the electrodynamic response of the two junction quantum interferometer in the presence of a time-dependent external flux. For this purpose, we shall take a sinusoidal forcing term, in such a way thatψex(t)=A+Bcosωt, where A is the normalized d. c. component of the applied flux and B is the normalized amplitude of the a. c. signal.

Figure 5.

Circulating current iS in the presence of an oscillating applied magnetic field shown as a function of the normalized time for A=0, B=0.1, iB=2.5, β=0.01 and for normalized frequencies equal to: a)ω˜=0.03; b)ω˜=0.06; c)ω˜=0.09. No initially trapped flux in the superconducting loop is present.

Now, sincet=Φ02πRIJτ, we can write

 ψex(τ)=A+Bcosω˜τ, (18)

whereω˜=Φ02πRIJω. By considering normalized frequencies of the order ofτϕ1, we may still use the analysis described above. In Section 5 we shall relax the latter hypothesis. As also specified above, a different approach, which will be developed in Section 5, needs to be used when very rapidly varying applied fields are applied to the system. We shall assume that the normalized amplitude B of the oscillating signal is much less than one (B<<1). The perturbation analysis is then carried out in a way at all similar as done above.

Figure 6.

a) Time average value ic of the critical current ic as a function of the amplitude B of the oscillating magnetic flux for A=0.1 (blue) and A=0.5 (orange). (b) Time average value ic of the critical current ic as a function of A for B=0.1 (blue) and B=0.5 (orange). Both figures are obtained by taking null values of the initially trapped flux.

We start by setting, by Eq. (9a) and (9b), ψ0(τ)=ψex(τ)=A+Bcos(ω˜τ)and ψ1=2sin(πψex(τ))cosϕ02πdψexdτ and solve the equations for the phase differences

 dϕ0dτ+cos(πψex(τ))sinϕ0(τ)=iB2,                                                         (a)dϕ1dτ+ϕ1(τ)cos(πψex(τ))cosϕ0(τ)=−πsin2(πψex(τ))sin(2ϕ0(τ))  (b) (19)

By noticing, however, thatcos(πψex(τ))=cos(πA+πBcosω˜τ)cos(πA)πBsin, (πA)cosω˜τEq. (19a) can be written in the following form:

 dϕ0dτ+(a−bcosω˜τ)sinϕ0(τ)=iB2, (20)

where a=cos(πA) andb=πBsin(πA). In Eq. (20) we find a perturbed solution in terms of the parameter b, so that, by settingϕ0(τ)=η0(τ)+bη1(τ), we can write:

 dη0dτ+asinη0(τ)=iB2,                                     (a)dη1dτ+acosη0(τ)η1(τ)=sinη0(τ)cos(ω˜τ)   (b) (21)

Notice then that the solutions to the above equations can be found by exactly the same procedure described for the case of constant applied fields. Once the solution for φ0(τ) is found, by substituting in Eq. (19b), the solution for ϕ1(τ) can be determined by solving a first order linear differential equation with time-dependent coefficients. Assuming thus ϕ0(τ)=η0(τ)+bη1(τ) to be a known expression, we can then write:

 iS=−2sin(πψex(τ))cos(η0(τ)+bη1(τ))+2πω˜Bsinω˜τ (22)

As before, the above expression is equal to ψ1(τ) and represents the circulating current iS in the circuit. In figs. 5a, 5b, and 5c the time dependence of the current iS for normalized frequency valuesω˜=0.03,0.06,0.09, respectively, and for A=0, β=0.01 and iB=2.5 is represented. In these graphs we notice that the oscillating patterns, which we have already detected in the constant applied field case, are modulated by the externally applied oscillating signal.

Another important quantity to be measured in these systems is the critical currentic, defined as the maximum value of the current bias iB which can be injected in the two junction interferometer without giving rise to dissipation. By considering the stationary case of Eq. (21a) we write:

 iB=2cos(πψex(τ))sinϕ0(τ). (23)

Therefore, we have

 ic=2|cos(πA+πBcos(ω˜τ))|. (24)

Noticing that the time-averaged value <ic> of the critical current does not depend on the normalized frequency, it can be calculated in terms of solely A and B, the results being shown in fig. 6a and fig. 6b for null values of the initially trapped flux. In particular, in fig. 6a <ic> is shown as a function of the applied magnetic field amplitude B, for A=0.1 and A=0.4, while in fig. 6b, <ic> vs. A curves are shown for B=0.1 and B=0.2. In the curves in fig. 6a we notice Fraunhofer-like oscillations, while ordinary cosinusoidal oscillations are present in fig. 6b.

For what seen above, the electrodynamic properties of a symmetric quantum interferometer containing two identical junctions with negligible capacitance can be studied by means of a perturbation approach in the parameter β, whose value gives the strength of the electromagnetic coupling between the two junction in the system. The analysis is rather similar to what done in other works in the literature (Grønbech-Jensen et al., 2003; Romeo & De Luca, 2004). However, in the present section we have presented a rather general procedure to obtain the solution to the problem to first order in the parameter β. Considering at first transient solutions, we have noticed that the function ψ(β,θ) governs fluxon dynamics, where θ is the ordinary time t, normalized to the characteristic circuital time constantτΨ=LR. By this more general approach it becomes evident that the characteristic time constant τϕ=Φ02πRIJ of the dynamics of the average superconducting phase difference ϕ is different from the fluxon dynamics characteristic timeτΨ, so that the asymptotic solution for the system, proposed in the analyses carried out in (Grønbech-Jensen et al., 2003) and in (Romeo & De Luca, 2004), acquires a more precise meaning in this context. Indeed, when the parameter β is sufficiently small to allow, for finite values of τ, an asymptotic evaluation of ψ0(τ2πβ) andψ1(τ2πβ), the general solution given in the present work coincides with the asymptotic perturbed solution proposed in (Grønbech-Jensen et al., 2003) and in (Romeo & De Luca, 2004) in the limit of negligible junction capacitance.

The perturbation analysis has been first carried out for a constant applied magnetic flux. Successively, since it could be experimentally possible to force the system with a time-dependent magnetic field, it is noted that the perturbed solution for the flux number ψ, obtained for a sinusoidal magnetic flux, needs careful evaluation. In order to exhibit experimentally detectable quantities, the circulating current iS is evaluated as a function of time, for different values of the frequency of the forcing field, whose a. c. component is assumed to be small. Finally, the time average <ic> of the critical current of the device has been studied both as a function of the d. c. component A and of the small amplitude B of the oscillating part of the applied flux. In these curves two characteristic behaviors have been detected: A Fraunhofer-like pattern in <ic> vs. B curves; independence of <ic> fromω˜. We shall see in Section 5 of the present Chapter that dependence from ω˜ appears in a more general case, i. e., when ω˜1 and the coefficient B is not assumed to be small.

3. Multi-junction quantum interference devices

In this section we shall consider the one-dimensional Josephson junction array (1D-JJA) represented in fig. 7, consisting of N+1 identical overdamped junctions connected in parallel. In this figure we notice that the bias current IB is evenly applied to the two external branches of the array. By assuming perfectly identical overdamped Josephson junctions with resistive parameter R and maximum Josephson currentIJ, we take the inductance L of the horizontal upper branches to be such thatβ=LIJΦ0<<1, where Φ0 is the elementary flux quantum. This parameter has been defined as each single loop could be compared to a two-junction quantum interference device. By fluxoid quantization, the normalized magnetic flux Ψk=ΦkΦ0 linked to the k-th cell of the array is related to the superconducting phase differences ϕk1 and ϕk across the two junctions in the cell, so that:

 ϕk−1−ϕk+2πΨk=2πnk, (25)

where nk is an integer andk=1,2,...,N.

Figure 7.

Schematic representation of a one-dimensional array of equal Josephson junctions.

When an external magnetic fieldH, orthogonal to the plane of the array, is applied to the system so that the normalized geometric applied flux through each cell isΨex=ΦexΦ0=μ0HS0Φ0, where S0 is the cell area, we may set

 Ψk=β∑m=1kim−βiB2+Ψex, (26)

fork=1,2,...,N, where ik=IkIJ is the normalized current flowing in the k-th branch andiB=IBIJ. The dynamical equation for each Josephson junction in the array is written by means of the resistively shunted junction (RSJ) model (Barone & Paternò, 1982) as follows:

 dϕkdτ+sinϕk=ik, (27)

where k=0,1,2,...,N andτ=2πRIJΦ0t. Eqs. (25-27) can be used to define the dynamics of the gauge-invariant superconducting phase difference ϕk in terms of the forcing parameters Ψex andiB=IBIJ=k=0Nik. In addition, the instantaneous voltage v(τ) of the system can be obtained by setting:

 v(τ)=1N+1∑k=0Ndϕkdτ. (28)

Define now the partial sum sn (1nN) of the normalized fluxes as follows:

 sn=∑k=1nΨk. (29)

By fluxoid quantization (Eq. (25)), setting allnk’s to zero (under the hypothesis of zero initially trapped flux in the array), we can write:

 ϕk=ϕ0+2πsk, (30)

fork=1,2,...,N, so that the dynamical equations (Eq. (27)) can be rewritten as follows:

 dϕ0dτ+sinϕ0=i0,                                        (a)2πdsndτ+sin(ϕ0+2πsn)=in−i0+sinϕ0  (b) (31)

where n=1,2,...,N in the second equation. Expressing now, by means of Eq. (26), the currents i0 and in in terms of the forcing terms Ψex and iB and of the partial sums of the flux variablessn, we may finally rewrite Eqs. (31a-b) as follows:

 dϕ0dτ+sinϕ0=iB2+s1−Ψexβ,                                                                  (a)2πdsndτ+sin(ϕ0+2πsn)−sinϕ0−sn+1−2sn+sn−1β=Ψex−s1β−iB2,    (b)2πdsNdτ+sin(ϕ0+2πsN)−sinϕ0+sN−sN−1β=Ψexβ−s1−Ψexβ           (c) (32)

where n=1,2,...,N1 in Eq. (32b).

The above analysis has been carried out essentially to write the dynamical equations in terms of the effective superconducting phase difference ϕ0 and of the N partial sumssn.

We shall now develop a reduction of these variables to one, by assuming small values of the parameter β. Therefore, start by considering the dynamical equations of the system as written in Eqs. (32a-b). For small values of the parameter β we can set:

 sn≈sn(0)+β sn(1). (33)

By substituting the above expression in Eq. (32b) and, after having multiplied both members by β, by setting the coefficients of βm (m=0,1) equal to zero, we obtain, forn=1,2,...,N:

 sn(0)=nΨex. (34)

For the first order corrections, on the other hand, we need to solve the following set of equations:

 {y1−y0=s2(1)−3s1(1)−iB2y2−y0=s3(1)−2s2(1)−iB2y3−y0=s4(1)−2s3(1)+s2(1)−s1(1)−iB2....yN−1−y0=sN(1)−2sN−1(1)+sN−2(1)−s1(1)−iB2yN−y0=sN−1(1)−sN(1)−s1(1). (35)

whereyn=sin(ϕ0+2πnΨex),n=0,1,2,...,N. By solving fors1(1), which is the quantity required in Eq. (32a), we have:

 s1(1)=−N−1N+1iB2−1N+1∑k=1N(yk−y0), (36)

Substitution of the above results into Eq. (32a) gives:

 dϕ0dτ+1(N+1)∑n=0Nyn=iB(N+1), (37)

The above differential equation represents the effective model for the 1D-JJA represented in Fig. 7, containing N+1 identical over-damped junctions connected in parallel. We notice that the reduced model in Eq. (37), even being of first order in β, does not explicitly contain a β1 term, given that Eq. (32a) contains β1 terms.

We can now explicitly perform the sum in Eq. (37), so that we write:

 dϕ0dτ+BN+1(N+1)sin(ϕ0+NπΨex)=iB(N+1), (38)

where the absolute value of BN+1=sin((N+1)πΨex)sin(πΨex) is the normalized critical current of the device in this approximation, as already known from the literature (Likharev, 1986). We start by finding the voltage, as defined in Eq. (28), under the assumption that Eqs. (32b-c) can be replaced by the effective model given by first-order perturbation analysis above. Therefore, we may set:

 v(τ)=dϕ0dτ=iB(N+1)−BN+1(N+1)sin(ϕ0+NπΨex). (39)

In this way, we can find the I-V characteristics by simply integrating Eq. (39), recalling the well-known procedure for a single overdamped junction. Indeed, noticing that

 〈v〉=1T∫0Tdϕ0dτdτ=2πT, (40)
(40)

where T is the period for the instantaneous voltage curve and a pseudo-period for the superconducting phase difference ϕ0(i.e.,ϕ0(τ+T)=ϕ0(τ)+2π), we first find the expression for φ0 from Eq. (39), for  iB>|BN+1|and k integer, so that

 ϕ0(τ)+NπΨex=2tan−1[BN+1iB−iB2−BN+12iBtan(12(N+1)iB2−BN+12τ)] +2kπ . (41)

The pseudo-period T of the above solution can be found by inspection, so that:

 T=2π(N+1)iB2−BN+12. (42)

Therefore, by Eqs. (40) and (42), the I-V characteristics are given by the following expression:

 iB=(N+1)2〈v〉2+BN+12, (43)

where only the positive branch has been chosen.

In Figs. 8a-b BN+1vs.Ψex curves are shown. In particular, in Fig. 8a the number of JJ’s in the system (N+1) is taken to be equal to 10, while in Fig. 8b it is set equal to 15. We notice that, when the normalized flux approaches the first zero in the BN+1 vs. Ψexcurve for(N+1)=10, namelyΨex=1/10, the system goes toward a purely resistive behaviour. In Fig. 8b, the same resistive behaviour is expected forΨex=1/15, given that the BN+1 vs. Ψexcurve for (N+1)=15 reaches its first zero exactly atΨex=1/15.

Figure 8.

Critical current iC=|BN+1| as a function of the applied magnetic flux for an array of 10 (a) and 15 (b) Josephson junctions.

In Figs. 9a-b I-V characteristics for different externally applied flux values are shown. In Fig. 9a the number of JJ’s is taken to be equal to 10, while in Fig. 9b it is set equal to 15. Starting from Fig. 9a, we notice that, as the normalized flux approaches the first zero in the BN+1 vs. Ψexcurve for(N+1)=10, namelyΨex=1/10, the behaviour of the system indeed becomes purely resistive. In Fig. 8b, the same resistive behaviour is obtained atΨex=1/15, as predicted.

Figure 9.

a) I-V characteristics of an array of 10 Josephson junctions for the following values of the applied magnetic fluxΨex: 0.05 (blue); 0.075 (cyan); 0.10 (orange); 0.15 (red). (b) I-V characteristics of an array of 15 Josephson junctions for the following values of the applied magnetic fluxΨex: 0.025 (blue); 0.05 (cyan); 0.075 (orange); 0.10 (red).

For fixed bias current values, the voltage versus flux curves can be obtained by Eq. (43) and is given by the following:

 〈v〉=iB2−[BN+1(Ψex)]2(N+1). (44)

The above expression is similar to the homologous d. c. SQUID v vs. Ψexcurves.

In Figs. 10a and 10b we report the v vs. Ψexcurve for (N+1)=10 and for(N+1)=15, respectively. The normalized bias current values are iB=4.0, 10.0, 16.0 (from bottom to top) for Fig. 10a and iB=4.0, 10.0, 15.0 (from bottom to top) for Fig. 10b. Notice that the zero-voltage state regions on the Ψex-axis disappear at values of the bias current greater or equal to the critical current of the system. In Fig. 10a, indeed, we see that disappearance of zero-voltage states occurs, for increasing normalized flux values, exactly atiB=10.0. On the other hand, in Fig. 10b, the same behaviour is detected atiB=15.0. Notice also that, for relatively high bias currents, the Ψex intervals in which finite-voltage states are present tend to be flatter than in the case of a low number of Josephson junctions. This feature is more evident in Fig. 10b, where the region of zero-voltage states is narrower than in Fig. 10a. Therefore, for high enough values of the the number of Josephson junction in the array and of the bias current, away from integer values of the normalized applied flux, the interval in which finite-voltage states spread can be approximated by a horizontal segment.

Figure 10.

Average voltage vs. applied magnetic flux for an array of 10 (a) and 15 (b) Josephson junctions for three different values of the normalized bias currentiB. In (a) iB=4.0, 10.0, 16.0(from bottom to top). In (b) iB=4.0, 10.0, 15.0(from bottom to top).

In conclusion, by considering the dynamical equations of one-dimensional arrays containing N+1 identical overdamped Josephson junctions, the system of N+1 nonlinear first-order ordinary differential equation equations can be broken into two coupled subsystems, one consisting of only one equation for the superconducting phase of one junction in the array (arbitrarily chosen to be the first), the second describing the time evolution of N opportunely defined normalized flux variables.

When a solution of the latter N equations is found, by means of a perturbative approach to first order in the parameter β, the dynamical properties of the system are described by a single time-evolution equation. In this way, we may affirm that, for small values of β, the system may be described by an equivalent single-junction model, where the maximum Josephson current is appropriately defined in terms of the normalized applied fluxΨex. When we compare our present analysis to equivalent studies carried out for a two-junction interferometer, we realize that the degree of approximation of the present model in β is one order less than the first-order perturbation analysis carried out for the simplest two-junction system. This is a consequence of the approach followed in the present work, where we had to appropriately define partial sums of flux variables in order to separate the dynamical equations into two subsystems. When we refer to the SQUID case, then, we might state that the present analysis corresponds exactly to the β=0 limit. Further work is therefore needed to carry out more detailed information on the system behaviour for finite values of β.

Even though part of the present analysis reproduces known results, as, for example, the expression for the maximum Josephson current, it still represents a simple way of approaching the problem of the electrodynamic response of one-dimensional arrays of overdamped Josephson junctions by an equivalent single-junction model. In fact, by starting with the simple representation of the dynamics of the system given in Eq. (38), all the results obtained for a single Josephson junction can be reproduced for an array of N+1 equal overdamped Josephson junctions. In addition, in case the solutions of the normalized flux variable would be extended to second order in the parameter β, following the same analysis as in the present section, effects due to finite β values in the electrodynamic properties of the system would be detected. Finally, considering that the present analysis has been carried out in the absence of flux fluctuations, its extension to noise effects can be obtained by means of already known results obtained for a single overdamped Josephson junction (Ambegaokar & Halperin, 1969; Bishop & Trullinger, 1978). In this case, however, care must be taken in considering the stochastic terms on all branches of the array.

4. Parallel connections of N × (0-π) overdamped Josephson junctions

In the previous section we have considered an array of N+1 overdamped 0-junctions, without considering the possibility of inserting π-junctions in the system. We briefly recall that π-junctions (Bulaevskii et al., 1977; Geshkenbein et al., 1987; Baselman et al., 1999; Ryazanov et al., 2001 M. Weides et al., 2006), when compared to 0-junctions, possess an intrinsic phase difference exactly equal to π. By inserting a 0-junction and a π-junction in the same superconducting loop, π-SQUIDs may be realized. These non-conventional SQUIDs can be fabricated either by exploiting the symmetry properties of d-wave superconductors (Chesca, 1999; Schultz et al., 2000) or by utilizing both s-wave and d-wave superconductors (Wollman et al., 1993; Smilde et al., 2004). A π-SQUID can thus be viewed as an elementary cell of a N×(0-π) one-dimensional array of overdamped Josephson junctions shown in fig. 11.

Figure 11.

Schematic representation of a Josephson junction array with alternate parameters. In the branches of the parallel connections overdamped 0-junctions (green circles) and π-junctions (red diamonds) are alternately present.

Therefore, π-SQUIDs can be viewed as the building block of discretized models of multifacets Josephson junctions (MJJs) (Scharinger et al., 2010), in which the critical current density alternates between two opposite values along the junction length. However, even though some characteristic features of MJJs can be qualitatively reproduced by N×(0-π) one-dimensional arrays, one should bear in mind that the latter are, in general, less complex systems than MJJs.

For conventional arrays of overdamped Josephson junctions we have already shown that, for small enough values of the characteristic parameter β a series solution for the magnetic flux variable can be found by perturbation analysis. In this way, the multi-junction interferometer model reduces to a single non-linear ordinary differential equation. The same perturbation approach will be proposed again in the present section to derive the equivalent single-junction model of N×(0-π) one-dimensional arrays of overdamped Josephson junctions. Therefore, we start by considering the model system represented in fig. 11, consisting of identical overdamped junctions connected in parallel. In this system one half of the bias current IB is evenly applied to the two external branches of the array. This condition is obtained, for example, by injecting the current by means of a superconducting bar of width w equal to the length of the array. In order to have well focused bias, the penetration length of the superconducting bar would be much smaller than its width w.

4.1. The homogeneous case

Consider, as a first approach to the problem, the loop areasSk, the junctions resistive parameters Rk and maximum Josephson currentsIJk, and the inductances Lk in the horizontal upper branches to be all equal. In this way one may writeSk=S0, Rk=R0, IJk=IJ0, Lk=L0, for all allowed values of k. Define, as in the previous section in this case,β=L0IJ0Φ0. Consider, finally, as also represented in fig. 11, that 0-junctions (green circles) occupy even positions (k=0,2,4,..,M1) and π-junctions (diamonds) occupy odd positions (k=1,3,..,M). The reader will notice very many similarities with the analysis in the previous section. However, due to the higher degree of complexity of the problem, we need to proceed step by step.

By fluxoid quantization, the normalized magnetic flux Ψk=ΦkΦ0 linked to the k-th cell of the array is related to the superconducting phase differences ϕk1 and ϕk across the two junctions in the cell as follows:

 ϕk−1−ϕk+2πΨk=2π[nk+(−1)k−14], (45)

where nk is an integer and k=1,2,...,M(notice that the number of loops is one less the number of junctions). If an external magnetic fieldH, orthogonal to the plane of the array, is applied to the system, the normalized geometric applied flux through each cell now isΨex=ΦexΦ0=μ0HS0Φ0. By considering Φk as the sum of the applied and induced flux, we may write:

 Ψk=β(∑m=0k−1im−iB2)+Ψex, (46)

fork=1,2,...,M, where ik=IkIJ0 is the normalized current flowing in the k-th branch andiB=IBIJ0. The dynamical equation for each Josephson junction in the array can be written by means of the RSJ model (Barone & Paternò, 1982), so that:

 dϕkdτ+sinϕk=ik, (47)

where k=0,1,2,...,M andτ=2πR0IJ0Φ0t. Eqs. (45-47) can be used to define the dynamics of the gauge-invariant superconducting phase difference ϕk in terms of the forcing parameters Ψex andiB=IBIJ=k=0Mik. Define now the partial sum sn (1nM) of the normalized fluxes as follows:

 sn=∑k=1nΨk. (48)

By fluxoid quantization, setting allnk’s to zero in Eq. (45) under the hypothesis of zero initially trapped flux in the array, we can write:

 ϕk=ϕ0+2πsk+π[1−(−1)k2], (49)

fork=1,2,...,M, so that the dynamical equations can be rewritten as follows:

 dϕ0dτ+sinϕ0=i0,                                                                              (a)2πds2ndτ+sin(ϕ0+2πs2n)=i2n−i0+sinϕ0, (k=2n)                  (b)2πds2n−1dτ−sin(ϕ0+2πs2n−1)=i2n−1−i0+sinϕ0, (k=2n−1)    (c) (50)

where n=1,2,...,M12 in (50b), and n=1,2,...,M+12 in (50c). Expressing now, by means of Eq. (46), the currentsi0, i2n, and i2n1 in terms of the forcing terms Ψex and iB and of the partial sums of the flux variables sk defined in Eq. (48), we may finally rewrite Eqs. (50a-c) as follows:

 dϕ0dτ+sinϕ0=iB2+s1−Ψexβ,                                                                               (a)2πds2ndτ+sin(ϕ0+2πsn)−sinϕ0−s2n+1−2s2n+s2n−1β=Ψex−s1β−iB2,          (b)2πds2n−1dτ−sin(ϕ0+2πs2n−1)−sinϕ0−s2n−2s2n−1+s2n−2β=Ψex−s1β−iB2,  (c)2πdsMdτ−sin(ϕ0+2πsM)−sinϕ0+sM−sM−1β=Ψexβ−s1−Ψexβ,                   (d) (51)

where now n=1,2,...,M12 in Eqs. (51b) and (51c) and where it has been sets0=0.

Considering the dynamical equations of the system as written in Eqs. (51a-d), for small values of the parameter β we may assume that the solution, to first order in this perturbation parameter, can be written as follows:

 sn≈sn(0)+β sn(1). (52)

By substituting the above expression in Eqs. (51b-d) and, after having multiplied both members by β, by setting to zero the coefficients of βm (m=0,1), we obtain:

 sn(0)=nΨex. (53)

For the first order corrections, on the other hand, we need to solve the following set of equations:

 {−y1−y0=s2(1)−3s1(1)−iB2y2−y0=s3(1)−2s2(1)−iB2−y3−y0=s4(1)−2s3(1)+s2(1)−s1(1)−iB2....−yM−2−y0=sM−1(1)−2sM−2(1)+sM−1(1)−s1(1)−iB2yM−1−y0=sM(1)−2sM−1(1)+sM−2(1)−s1(1)−iB2−yM−y0=sM−1(1)−sM(1)−s1(1). (54)

whereyn=sin(ϕ0+2πnΨex),n=0,1,2,...,M. By solving fors1(1), which is the required quantity in Eq. (51a), we have:

 s1(1)=−M−1M+1iB2−1M+1∑k=1M[(−1)kyk−y0], (55)

Substitution of the above result into Eq. (51a) gives:

 dϕ0dτ+1(M+1)∑k=0M(−1)ksin(ϕ0+2πkΨex)=iB(M+1). (56)

It is now possible to explicitly calculate the finite sum in Eq. (56) to get, in terms of the number N=M+12 of the individual (0-π) cells:

 dϕ0dτ−A2N2Ncos[ϕ0+(2N−1)πΨex]=iB2N, (57)

whereA2N=sin(2NπΨex)cos(πΨex). The above equation represents the equivalent single-junction model for the homogeneous array consisting of N cells, each one containing one 0-junction and one π-junction.

4.2 The non-homogeneous case

Consider, next, a non-homogeneous array with alternating 0-π Josephson junctions. We shall take the parameters of all 0-junctions equal. The parameters of π-junctions, even being equal among them, are assumed to be different from those of the 0-junctions. In this case we can omit some of the calculations, having already treated the problem in detail in the previous subsection.

Considering again the 1D-JJA represented in fig. 11, we now assume that the loop areasSk, the junctions resistive parameters Rk and maximum Josephson currentsIJk, and the inductances Lk of the horizontal upper branches alternate in their values, as we go along the array. In this way we write

 Sk={S0                       k  evenS1=σS0             k  odd,   (a)Rk={R0                       k  evenR1=αS0             k  odd,  (b)IJk={IJ0                      k  evenIJ1=εIJ0            k  odd,   (c)Lk={L0                        k  evenL1=σS0              k  odd,  (d) (58)

for all allowed values of k. In this way, the additional parameters α, ε, and σ are implicitly defined. As before, we set τ=2πR0IJ0Φ0t andβ=L0IJ0Φ0; therefore, we may defineβ1=L1IJ1Φ0=εαβ.

Fluxoid quantization give the same relation as in Eq. (45) between the superconducting phases and the normalized magnetic flux Ψk linked to the k-th cell of the array. However, Eq. (46) is modified as follows

 Ψk={β(∑m=0k−1im−iB2)+Ψex k evenβσε(∑m=0k−1im−iB2)+σΨex k odd (59)

fork=1,2,...,M. The equations of the motion for the superconducting phases can now be written as follows:

 dϕ2ndτ+sinϕ2n=i2n,                   (a)dϕ2n−1dτ+εαsinϕ2n−1=αi2n−1,  (b) (60)

where n runs over all allowed k-values also in all following equations.

By adopting a first-order perturbation analysis in the parameter β, we set:

 s2n≈s2n(0)+β s2n(1)=n(σ+1)Ψex+β s2n(1),                      (a)s2n−1≈s2n−1(0)+β s2n−1(1)=[n(σ+1)−1]Ψex+β s2n−1(1),  (b) (61)

By substituting the above expression in Eqs. (49) to obtain the superconducting phases ϕk we then write the differential equations (60a-b) in terms of the sole phase variable ϕ0 and of the 2N normalized flux variablessk. By following the same steps as in the previous section, we obtain the following equivalent single-junction model for the superconducting array:

 dϕ˜0dτ+αA˜N(α+1)N[sinϕ˜0−εsin(ϕ˜0+2πσΨex)]=αiB(α+1)N, (62)

Where here A˜N=sin(πN(σ+1)Ψex)sin(π(σ+1)Ψex) andϕ˜0=ϕ0+π(N1)(σ+1)Ψex. Naturally, for α, ε, and σ all equal to one, the ordinary differential equation (62) reduces to Eq. (57).

4.3. Critical current

In order to find the critical current of arrays with alternating 0-π Josephson junctions, we proceed as follows. First, consider the homogeneous array described in Section 4.1. We look for the maximum value of the bias current iB which can be injected in the system at zero voltage (dϕ0dτ=0in Eq. (57)). Therefore, by maximizing with respect to ϕ0 the following expression

 iB=−A2Ncos[ϕ0+(2N−1)πΨex], (63)

we can express the critical current of the homogeneous device as follows:

 ic=|A2N|=|sin(2NπΨex)cos(πΨex)|. (64)

In order to understand the origin of the patterns we are going to show for the non homogeneous case, let us consider the result in Eq. (64) as the product of the envelop function

 E(Ψex)=1+ε2−2εcos(2πσΨex), (65)

(ε=1,σ=1in this case) and of the rapidly oscillating function

 F(Ψex)=|sin(Nπ(σ+1)Ψex)sin(π(σ+1)Ψex)|, (66)

withσ=1, giving rise to secondary peaks in the overall curves. The functions NE(Ψex) and 2F(Ψex) are separately shown in fig. 12a for ε=1,σ=1 and N=10. The factors in front of these functions are chosen in such a way that they can attain the same maximum value of the productE(Ψex)F(Ψex).

Figure 12.

a) Envelope function NE(Ψex) as a function of2F(Ψex). (b) Critical current as a function of the normalized applied flux for a homogeneous parallel array with N=12, ε=1, σ=1.

One notices that regular primary peaks of F(Ψex) appear at integer and half integer values ofΨex. By multiplying the functions E(Ψex) and F(Ψex) as reported in fig. 12b, we notice that primary peaks positioned at half-integer Ψex values are left unchanged, while those appearing at integer Ψex values are depressed to zero.

By proceeding in the same way for the non-homogeneous case, starting from Eq. (62) we find that the critical current of the non homogeneous array described in Section 3 can be written as follows:

 ic=E(Ψex)F(Ψex). (67)

Notice that the resulting pattern does not depend on α, which can be absorbed, in the effective dynamical equation, by a rescaling of the normalized time τ. Notice also that the periodicity of the envelop function E(Ψex) may now not be commensurable with that of the rapidly oscillating term|A˜N|. Consider, in fact, the case reported in figs. 13a and 13b, where the pattern is shown, as a full line, for a non homogeneous array with N=12, ε=1.3andσ=1.5.

Figure 13.

Critical current as a function of Ψex for a non-homogeneous parallel array with N=12, ε=1.3, σ=1.5 in the following ranges: (a) [-1.25, 1.25]; (b) [-2.5, 2.5]. In (a) the envelope function NE(Ψex) is shown as a dotted line.

In fig. 13a the envelop curve NE(Ψex) is shown as a dotted line. In the interference pattern in fig. 13a we may notice the finiteness of the small peak at Ψex=0 due to the non-vanishing value ofE(0)=|ε1|. Also, notice the reduction in height of the primary peaks close to half-integer values as compared to those in fig. 12b. Finally, notice the appearance of two extra peaks of equal height in between two successive primary peaks. In fig. 13b, where the range of Ψex is increased, we notice that this feature repeats over a periodΔΨex=2, as it can be directly confirmed by inspection of Eqs. (65-67). Differently from fig. 13a and 13b, when we consider an array with N=12, ε=1.3andσ=2, as in figs. 14a-c, the fundamental periods of the functions E(Ψex) and F(Ψex) are incommensurable, so that the overall pattern is not periodic. This feature can be detected by gradually increasing the range of the pattern, as it is done in fig. 14a and fig. 14b. In this way, when looking for the conditions giving periodicity, one realizes that the ratio σ+1σ needs to be a rational number. Therefore, when the parameter σ is a rational number, one has periodicity in the ic vs. Ψexcurves; otherwise, irregular patterns, like those shown in figs. 14a-b, are found.

Figure 14.

Critical current as a function of Ψex for a non-homogeneous parallel array with N=12, ε=1.3, σ=2in the following ranges: (a) [-1.75, 1.75]; (b) [-5.5, 5.5]. In (a) the envelope function NE(Ψex) is shown as a dotted line.

We have seen that a reduced single-junction model can be adopted to describe the overall dynamics of a one-dimensional arrays with alternating parameters of N×(0-π) Josephson junctions. This effective model is very useful, since it allows to obtain the critical current vs. normalized magnetic flux curves in closed analytic form. The interference patterns are seen to be qualitatively similar to recently obtained experimental results on multifacets Josephson junctions (Scharinger et al., 2010) in which the critical current density alternates many times between two opposite values along the junction length. Discrete Josephson junction arrays, even presenting some analogies with the latter devices, are much too simple systems to describe the complete behaviour of MJJs. As a matter of fact, shielding current effects is not taken into account by analysis carried out in this section in the lowest order approximation in β. Nevertheless, the analytic results obtained for the interference patterns shed some light on the causes of the presence or absence of periodicity and on the nature of primary and secondary peaks in the ic vs. Ψexcurves. We finally notice that the analysis relies on the choice of an even number of Josephson junction in the array. Different patterns are expected for an odd number of alternating junctions in the array, depending also on which type (0 or π) of junctions is predominant. Further work is therefore necessary to address this problem in its fullest extent.

5. Quantum interferometers in the presence of rapidly varying fields

In Section 2 we have analyzed the effective model for a two-junction quantum interferometer in the presence of an oscillating magnetic flux under the hypothesis that the frequency of oscillation is comparable with the inverse of the characteristic time evolution τϕ of the superconducting time variable ϕ. In this way, the quasi-static approach described in Section 2 has been proven to be applicable. In the present section, on the other hand, we shall consider rapidly varying externally applied fluxes, whose frequency ω is comparable withτψ1, so thatω>>τϕ1.

Let us therefore consider an externally applied flux having d. c. component A and a. c. amplitude B, so that

 ψex(t)=A+Bsinωt (68)

where ωτψ1>>τϕ1 is the frequency of the sinusoidal term. Let us again consider Eqs. (4a-b) rewritten, for n=0, as follows:

 dϕdτ+cosπψsinϕ=iB2,               (a)πdψdτ+sinπψcosϕ=ψex−ψ2β   (b) (69)

where the normalization τ=2πRIJΦ0t=tτϕ prescribes a τ-dependence of the externally applied flux as written in Eq. (18) with a normalized frequencyω˜=Φ02πRIJω=ωτϕ. We therefore need to consider again all steps in Section 2, having care to integrate opportunely the right hand side term of (69b), in order to obtain a solution for Ψ in terms of the superconducting phase by perturbation analysis on β for arbitrary values of A and B. In this way, the effective dynamics for φ can be found when the solution for Ψ is substituted in the cosine term of (69a).

Start by considering the cosine term in (69b) as a quasi-static quantity (i. e., it does not vary appreciably over an interval of time of the order ofτψ). This hypothesis is confirmed by what already stated in the previous section; i. e., while the variables φ varies on a characteristic time interval Δτϕ, the variable Ψ varies within a time interval ΔτΨ=βΔτφ «Δτφ. Within the former time interval Δτφ it is then possible to choose a subinterval, of the order of Δτψ, in which the variable φ does not vary appreciably. We can thus solve (69b), by perturbation analysis, by first setting τ=βθ and by rewriting it as follows:

 dψdθ+2βsinπψcosϕ+ψ(θ)=ψex(θ). (70)

By now setting

 ψ(θ)=ψ0(θ)+βψ1(θ) (71)

we again find the ODEs for ψ0 and ψ1 in (8a) and (8b). Recall that an ODE of the type

 dfdθ+f(θ)=g(θ), (72)

has solution f(θ)=e−θ θg(x)exdx. By now considering (68), by taking the non-decaying solutions of the system of ordinary differential equations (8a-b) and by considering the non-vanishing solution of (8a) for ψ0 at large values of θ, we have

 ψ0=A+B1+ω˜2(cosω˜θ+ω˜sinω˜θ). (73)

Having found the solution to (8a), we can find the solution to (8b) by the same type of reasoning. After some rather long calculations one finds

 ψ1=−2h(θ)cosϕ, (74)

where

 h(θ)=∑n=−∞+∞∑k=−∞+∞gn,k[sinαn,k(θ)−ω˜cosαn,k(θ)], (75)

withgn,k=Jn(γ)Jn(γω˜)1+(n+k)2ω˜2, Jn(x)being the Bessel function of order n, andαn,k(θ)=πA+(n+k)ω˜θ+nπ2. Therefore, the SQUID dynamics can be described, to first order in β, by the following ODE:

 dϕdτ+X0(θ)sinϕ+πβh(θ)Y0(θ)sin2ϕ=iB/2m, (76)

whereX0(θ)=cosπψ0 andY0(θ)=sinπψ0.

Equation (76) thus represents the differential equation describing the dynamics of the superconducting phase difference ϕ in a d. c. SQUID in the presence of a time-varying externally applied flux, whose frequency ω is considered to be comparable with τψ-1, in such a way thatω˜=ωτψ1. For slowly varying fields (ω˜<<1) one can readily verify from (8a) and (8b) that ψ0(θ)ψex(θ) andψ12sinψexcosϕ, respectively. In this way, the dynamics described by the quasi-static d. c. SQUID in (3) holds.

Forω˜=ωτψ1, time evolution of the two variables, ϕ and ψ, still occurs with two completely different time scales. In fact, as already stated, one has τψ=βτϕ «τϕ. Therefore, flux motion is very fast with respect to the dynamics of the phase variable ϕ. The only difference, here, is that the externally applied flux ψex(θ) is able to follow this fast dynamics. Having carefully solved (8a) and (8b), and having found the effective single-junction dynamical equation for a d. c. SQUID, we can determine the effective time-averaged equation, by taking the time average over the fast variable ψ. Therefore, we may write

 dϕdτ+〈X0(θ)〉sinϕ+πβ〈h(θ)Y0(θ)〉sin2ϕ=iB/2, (77)

where the symbol <x> stands for the time average of the variable x. Equation (77) can thus be considered an effective single-junction model for a d. c. SQUID in the presence of a rapidly varying magnetic field (ω˜=ωtL1). The average values X0(θ) and h(θ)Y0(θ) can be calculated as follows. First of all, set

 X0(θ)=Re{exp(iπψ0)}. (78)

Let us next express the exponential of a cosine and a sine terms in exp(iπψ0) through the following Bessel function identities

 eiacosx=∑n=−∞+∞inJn(a)einx, (a)eiasinx=∑n=−∞+∞Jn(a)einx.     (b) (79)

In this way, (78) becomes

 〈X0(θ)〉=Re{∑n,kJn(γ)Jk(γω˜)〈eiαn,k(θ)〉}. (80)

It is now easy to show that

 〈eiαn,k(θ)〉=eiπAinδn,−k, (81)

where the symbol δn,m is the Kronecker delta. By inserting (81) in (80), we have

 〈X0(θ)〉=cosπAρ(ω˜,B), (82)

with

 ρ(ω˜,B)=J0(γ)J0(γω˜)+2∑n=1+∞(−1)nJ2n(γ)J2n(γω˜). (83)

Proceeding in a similar way in finding the effective coefficient of the sin2ϕ term, we find:

 〈h(θ)Y0(θ)〉=12∑n,m,l=−∞+∞ξn,m,lgn,mJl(γ)Jn+m−l(γω˜). (84)

where

 ξn,m,l={(−1)n−l                             for (n-l) even−(n+k)(−1)n−l−12             for (n-l) odd. (85)

Having expressed the effective time-averaged terms in (77) in a closed analytic form, we can understand the effect of a high-frequency field on the electrodynamic behaviour a d. c. SQUID with extremely small value of the parameter β, for instance. The critical current ic of the device in this case (β = 0) can be expressed as follows:

 iC(ω˜,A,B)=2|cosπA||ρ(ω˜,B)|, (86)

where the quantity |ρ(ω˜,B)| is the extra-factor modifying the usual form of a d. c. SQUID, expressed, in terms of a constant applied flux A, as2|cosπA|. The above expression can be derived by inspection from Eq. (77), by setting β = 0 and getting the maximum stationary value for iB with respect to ϕ.

In fig. 15a-b we thus show the critical current ic in terms of the a. c. component B of the applied flux. In particular, in fig. 15a, we report the ic vs. B curves for various values of the d. c. component A of the applied flux. normalized frequencyω˜. In fig. 15b, on the other hand, the ic vs. B curves are shown for various values of the normalized frequencyω˜. From figs. 15a-b we notice Fraunhofer-like patterns of the critical current as shown as a function of the a. c. amplitude B.

The particular shape of these patterns in figs. 15a-b depends on the value of the normalized frequencyω˜. When the d. c. component A and the time-varying portion of the magnetic flux attain fixed values, and we let the normalized frequency vary with continuity, the ω˜ critical current of the device is the same as that of the quantity|ρ(ω˜,B)|. The ic vs. ω˜curves are represented in fig. 16 for A=0.1 and for various values of B. In this respect, we notice that, while for small fixed B values the quantity |ρ(ω˜,B)| is always increasing for increasing values of ω˜(see orange line in fig. 16) for values of B approaching 1.0, this character is lost (see brown and cyan lines in fig. 16).

Figure 15.

Critical current ic as a function of the oscillating amplitude B. (a) The normalized frequency is fixed at ω˜=1.0 and various values of the d. c. component A are represented: A=0.1 (orange); A=0.2 (brown); A=0.3 (cyan). (b) The d. c. component A is fixed at A=0.1 and the normalized frequency attains the following values: ω˜=0.5(orange); ω˜=1.0(brown); ω˜=1.5(cyan).

Figure 16.

Critical current ic as a function of the normalized frequency ω˜ for A=0.1 and for various values of the oscillating amplitude B: B=0.4 (orange); B=0.8 (brown); B=1.2 (cyan).

We can now calculate, for β=0, the flux-voltage curves (v vs. A) by the following well known expression (Barone & Paternò, 1982):

 v=〈dϕdτ〉=iB24−cos2(πA)ρ2(ω˜,B), (87)

foriB>2|cos(πA)ρ(ω˜,B)|. The important feature in this expression is that these curves depend both on B and ω˜ though the extra-factor|ρ(ω˜,B)|. Because of this dependence, the amplitude of the v vs. A curves can be varied and B and ω˜ can be viewed as control parameters. In figs. 17a-b we thus report the v vs. A curves foriB=2.5, for ω˜=0 (a) and ω˜=1.0 (b), and for various values of B. In particular, in fig. 17a we notice that the amplitude of the v vs. A curves obtained at B=0 and ω˜=0 (orange line) decreases as we let B increase to 0.15 first (brown line) and to 0.30 next (cyan line). The same decreasing behavior is detected in fig. 17b for ω˜=0.1 when B increases from 0.0 (orange line) to 0.15 (brown line) and to 0.30 (cyan line).

Figure 17.

Voltage v versus the d. c. component A of the applied flux for iB=2.5 and B=0.0 (orange); B=0.15 (brown); B=0.30 (cyan). In (a) the normalized frequency values isω˜=0, in (b)ω˜=1.0.

In fig. 18a-b, finally, by fixing the value of B first to 0.5 (a) and then to 1.0 (b), we notice that the amplitude of the v vs. A curves increases for increasing values ofω˜, as shown for ω˜=0.5 (orange line), ω˜=1.0(brown line), and ω˜=1.5 (cyan line).

Figure 18.

Voltage v versus the d. c. component A of the applied flux for iB=2.5 and ω˜=0.5 (orange); ω˜=1.0(brown); ω˜=1.5(cyan). In (a) the a. c. component of the applied magnetic flux is B=0.50, in (b) B=1.0.

6. Conclusion

We have studied the dynamical properties of quantum interferometers consisting of single or multiple superconducting loops, each containing two Josephson junctions.

A symmetric quantum interferometer containing two identical junctions with negligible capacitance has been considered first. The analysis of the system has been carried out by means of a perturbation approach in the parameter β, whose value gives the strength of the electromagnetic coupling between the two junction in the system. We have noticed that the flux-number function ψ(β,θ) governs fluxon dynamics, where θ is the laboratory time t normalized to the characteristic circuital time constantτψ=LR. By this general approach it becomes evident that the characteristic time constant τϕ=Φ02πRIJ of the dynamics of the average superconducting phase difference φ is different from the characteristic timeτψ. The perturbation analysis has been carried out for both a constant and a time-dependent applied magnetic flux. The circulating current iS and the time average of the critical current <ic> as a function of the d. c. and a. c. components of the applied flux are evaluated in the adiabatic limit, assuming that the oscillation frequency ω of the applied flux is much less thenτψ1. In this limit the Fraunhofer-like pattern in <ic> vs. B curves are shown to be independent from the normalized frequancyω˜=ωτψ.

Next, the dynamical equations of one-dimensional arrays containing N+1 identical overdamped Josephson junctions are considered. It has been noticed that the system of N+1 nonlinear first-order ordinary differential equations can be broken into two coupled subsystems, one consisting of only one equation for the superconducting phase of one junction in the array (arbitrarily chosen to be φ0), the second describing the time evolution of N opportunely defined normalized flux variables. When a solution of the latter N equations is found, by means of a perturbative approach to first order in the parameter β, the dynamical properties of the system are described by a single time-evolution equation for φ0. In this way, we may affirm that, for small values of β, the system may be described by an equivalent single-junction model, where the maximum Josephson current is appropriately defined. The analysis represents a simple way of approaching the problem of the electrodynamic response of one-dimensional arrays of overdamped Josephson junctions by an equivalent single-junction model.

The same approach is followed for one-dimensional arrays with alternating parameters of N×(0-π) Josephson junction. Even in this case an effective single-junction model can be adopted to describe the overall dynamics of the system. This model is very useful, since it allows us to obtain the critical current vs. normalized magnetic flux curves in closed analytic form. The interference patterns are seen to be qualitatively similar to recently obtained experimental results on multifacets Josephson junctions (Scharinger et al., 2010) in which the critical current density alternates many times between two opposite values along the junction length. The analytic results for the interference patterns clarify the presence or absence of periodicity and the nature of primary and secondary peaks in these curves. Further investigation on the dependence of ic on different distributions of the JJs in the array can be of interest.

Finally, by allowing the magnetic flux, applied to a two-junction superconducting quantum interference device, to have an a. c. component in addition to a constant term A, we derive the effective reduced single-junction model describing the dynamics of the average superconducting phase difference ϕ of the two junctions in the device. The difference between this case and the one previously treated is that the alternating flux now varies with a frequency ω of the same order of magnitude ofτψ1, so that the adiabatic approach does not apply. The single-junction model for the system is again obtained by perturbation analysis to first order in the parameter β. Averaging of the rapidly varying quantities in the differential equation for ϕ gives the effective dynamics of the two junctions in the system. In particular, for β = 0, the critical current of the device is seen to depend on A, on the frequency ω˜ and the amplitude B of the a. c. component of the applied magnetic flux in a closed analytic form. From the analysis of the voltage vs. applied flux curves it can be argued that the quantities ω˜ and B can play the role of additional control parameters in the device. Further work in extending the present analysis to finite values of β is necessary. Experimental work confirming the predictions of the present analysis needs to be performed. As far as non-normalized quantities are concerned, for direct experimental confirmation of the present results, we finally notice that the junction dynamics evolves with characteristic frequencies the order of 1 THz. Therefore one needs to run the experiment with very rapidly oscillating signals (10 THz or more) in such a way that normalized frequencies of ω˜=1.0 can be achieved.

References

1 - V. Ambegaokar, B. I. Halperin, 1969Voltage Due to Thermal Noise in the dc Josephson Effect. Phys. Rev. Lett., 2225June 1969), 136413660031-9007
2 - A. Barone, G. Paternò, 1982Physics and Applications of the Josephson Effect, John Wiley & Sons, 0-47101-469-9York, USA
3 - J. J. Baselmans, A. F. Morpurgo, B. J. Van Wees, M. Klapwijk, 1999Vortices with Half Magnetic Flux Quanta in “Heavy Fermion” Superconductors. Nature, 3976714January 1999), 43450028-0836
4 - A. R. Bishop, S. E. Trullinger, 1978Josephson Junction Threshold Viewed as a Critical Point. Phys. Rev. B, 175March 1978), 217521821098-0121
5 - M. F. Bocko, A. M. Herr, M. J. Feldman, 1997Prospect for Quantum Coherence Computation Using Superconducting Electronics. IEEE Transactions on Applied Superconductivity, 72June 1997), 363836410000-1051
6 - L. N. Bulaevkii, V. V. Kuzii, A. A. Sobyanin, 1977Superconducting System with Weak Coupling to the Current in the Ground State. JETP Lett., 257April 1997), 2902940021-3640
7 - B. Chesca, 1999Magnetic field dependencies of the critical current and of the resonant modes of dc SQUIDs fabricated from superconductors with s+idx2-y2 order parameter symmetries. Ann. Phys. (Leipzig), 86September 1999), 5115221521-3889
8 - J. Clarke, A. I. Braginsky, 2004The SQUID Handbook, Vol. I, Wiley-VCH, 3-52740-229-2Germany
9 - D. S. Crankshaw, T. P. Orlando, 2001Inductance Effect in the Persistent Current Qubit. IEEE Transactions on Applied Superconductivity, 111March 2001), 100610091051-8223
10 - R. De Luca, 2011Quantum interference in parallel connections of N×(0-π) overdamped Josephson junctions. Supercon. Sci. Technol., 246June 2011), 0650260953-2048
11 - V. B. Geshkenbein, A. I. Larkin, A. Barone, 1987Vortices with Half Magnetic Flux Quanta in “Heavy Fermion” Superconductors. Phys. Rev. B, 361July 1987), 2352381098-0121
12 - N. Grønbech-Jensen, D. B. Thompson, M. Cirillo, C. Cosmelli, 2003Thermal escape from zero-voltage states in hysteretic superconducting interferometers. Phys. Rev. B, 6722June 2003), 2245051098-0121
13 - K. K. Likharev, 1986Dynamics of Josephson Junctions and Circuits, Gordon and Breach, 2-88124-042-9The Netherlands
14 - F. Romeo, R. De Luca, 2004Effective Non-Sinusoidal Current-Phase Relations in Conventional d. c. SQUIDs. Physics Letters A, 3284-5August 2004), 3303340375-9601
15 - F. Romeo, R. De Luca, 2005A reduced model for one-dimensional arrays of overdamped Josephson junctions. Physica C, 4323-4November 2005), 1591660921-4534
16 - V. V. Ryazanov, V. A. Oboznov, A. Rusanov, Yu, A. V. Veretennikov, A. A. Golubov, J. Aarts, 2001Coupling of Two Superconductors Through a Ferromagnet: Evidence for a π Junction. Phys. Rev. Lett., 8611March 2001), 242724300031-9007
17 - S. Scharinger, C. Gürlich, R. G. Mints, M. Weides, H. Kohlstedt, E. Goldobin, D. Koelle, R. Kleiner, 2010Interference patterns of multifacet 20x(0-π) Josephson junctions with ferromagnetic barrier. Phys. Rev. B, 8117May 2010), 1745351098-0121
18 - R. R. Schultz, B. Chesca, B. Goetz, C. W. Schneider, A. Schmehl, H. Bielefeldt, H. Hilgenkamp, J. Mannhart, C. C. Tsuei, 2000Design and Realization of an all d-wave dc π Superconducting Quantum Interference Device. Appl. Phys. Lett., 767February 2000), 9129140003-6951
19 - H. J. H. Smilde, Ariando, D. H. A. Blank, H. Hilgenkamp, H. Rogalla, 2004SQUIDs based on Josephson Contacts Between High-Tc and Low-Tc Superconductors. Phys. Rev. B, 702July 2004), 0245191098-0121
20 - M. Weides, M. Kemmler, E. Goldobin, D. Koelle, R. Kleiner, H. Kohlstedt, . Buzdin, A. , 2006High quality ferromagnetic 0 and π Josephson tunnel junctions, Appl. Phys. Lett., 8912September 2006), 1225110000-0003
21 - D. A. Wollman, D. J. Van Harlingen, W. C. Lee, D. M. Ginsberg, A. J. Leggett, 1993Experimental Determination of the Superconducting Pairing State in YBCO from the phase coherence of YBCO-Pb dc SQUIDs. Phys. Rev. Lett., 7113September 1993), 213421370031-9007