arxiv: v1 [cond-mat.supr-con] 4 Oct 2014

Size: px

Start display at page:

Download "arxiv: v1 [cond-mat.supr-con] 4 Oct 2014"

Arnold Gibson
5 years ago
Views:

1 Effect of current injection into thin-film Josephson junctions V. G. Kogan, and R. G. Mints, Ames Laborator, US Department of Energ, Ames, Iowa 5, USA The Ramond and Beverl Sackler School of Phsics and Astronom, Tel Aviv Universit, Tel Aviv 69978, Israel arxiv:.88v [cond-mat.supr-con] Oct New thin-film Josephson junctions have recentl been tested in which the current injected into one of the junction banks governs Josephson phenomena. One thus can continuousl manage the phase distribution at the junction b changing the injected current. A method of calculating the distribution of injected currents is proposed for a half-infinite thin-film strip with source-sink points at arbitrar positions at the film edges. The strip width W is assumed small relative to Λ = λ /d, λ is the bulk London penetration depth of the film material, d is the film thickness. PACS numbers: v Ec, w, 85.5.Cp I. INTRODUCTION In recent ears, the phsics of Josephson phenomena enjoed a number of important developments. Introduction of π and -π junctions, various was to have a different from π phase shift,,3 effect of vortices in the junction vicinit, 6 to name a few. Striking improvements in managing Josephson phenomena came after injection of currents into one of the thin-film banks was introduced that allowed for continuous control of the phase difference on the junction 7 9 and, in particular, to imitate the -π behavior. This development necessitates evaluation of the injected supercurrent distribution in one of the junction banks since this determines the distribution of the superconducting phase. S= S=I S= FIG.. (a) A sketch of two semi-infinite thin-film strips forming the Josephson junction in the overlapping shaded region. (b) The top film; arrows show positions of injection point contacts. (c) The stream function S at the film edge. In one of common realizations, the junction is formed b two half-infinite thin-film strips with overlapped edges. Etra current injectors are attached to the edges (a) (b) (c) of one of the films, shown schematicall in Fig.. The injected current affects the phase distribution in the thinfilm bank where it flows and thus the phase difference on the junction. We show in this communication that for sufficientl thin films with the size W smaller than the Pearl length Λ = λ /d the problem of the injected currents can be solved under ver general assumptions, so that the design of junctions with needed properties becomes possible. II. STREAM FUNCTION Consider a half-infinite thin-film strip of a width W Λ = λ /d where λ is the London penetration depth of the film material and d is the film thickness. Choose along the strip, < <, and across so that < < W, Fig. b. Let the injection points be at (, ) and (, ) at the film edge. The London equation integrated over the film thickness reads: h z + πλ c curl z g =. () Here, g is the sheet current densit and h is the self-field of the current g. The Biot-Savart integral for h z in terms of g shows that h z is of the order g/c, whereas the second term on the left-hand side of Eq. () is of the order gλ/cw g/c. Hence, in narrow strips with W Λ, the self-field can be disregarded. Introducing the scalar stream function S via g = curl[s(, )ẑ], we obtain instead of Eq. (): S =. () Phsicall, this simplification comes about since in narrow films the major contribution to the sstem energ is the kinetic energ of supercurrents, while their magnetic energ can be disregarded. The boundar condition of zero current component normal to edges, e.g., g = S = at the edge =, translates to S = constant along the edges. This constant, however, is not necessaril the same everwhere,

2 in particular, it should eperience a finite jump at injection points. Consider a point contact at the edge as illustrated at Fig., take its position as the origin of polar coordinates, and integrate the r component of the current g along the small half-circle centered at the injection point. The total injected current is: I = π g r r dφ = π S dφ = S(π) S(). (3) φ Hence, on two sides of the injection point the stream function eperiences a jump equal to the total injected current. Clearl, at the current sink S-jump has the opposite sign. Thus, we can choose S = everwhere at the edges of the half-infinite strip, ecept the segment between the injection and sink points, where S = I. To solve the Laplace equation () we first emplo the conformal mapping of the half-strip to a half-plane: 6, u + iv = i cosh π( + i). () It is seen that the half-plane u > is transformed to the half-strip of a width (hereafter we use W as a unit length). Eplicitl, this transformation reads: u =, v = cosh π cos π. (5) Hence, we have to solve the Laplace equation on a halfplane u > subject to boundar conditions S = I at the edge u = in the interval v < v < v with v = cosh π cos π, v = cosh π cos π, (6) and S = otherwise. To proceed, we first write the step-function S(, v) at the edge u = as a Fourier integral: S(, v) = S(, k) = I v v dk π S(, k)eikv, (7) dv e ikv = I ik ( e ikv e ikv). (8) Since S(, v) is a linear superposition of plane waves e ikv, we first consider the solution of the Laplace equation ( u + v)s(u, v) = subject to the boundar condition s(, v) = e ikv. Separating variables we obtain s(u, v) = e ikv e k u. Hence, the solution for the actual boundar condition is: S(u, v) = dk π S(, k)eikv k u. (9) Substituting here S(, k) of Eq. (8) one obtains: S(u, v) = I ( tan v v tan v v ). () π u u It is seen that at u, S = I if v < v < v and S = othewise, as it should be. Thus, we have the steam function at the half-plane u > for arbitrar positions v and v of the current contacts at the edge u =. We now can go back to the (, ) plane and specif the injection positions. A. Injectors at the edge = Let the injector and the sink be at (, ) and (, ). We obtain: S(, ) I/π = cos π tan cosh π cos π tan cos π cosh π cos π. () The lines of the current g are given b g dr = or b S d + S d = ds =, in other words, b contours of S = const. An eample is shown in Fig FIG.. The current distribution for the injection points (,.3) and (,.7). B. The injector at = and the sink at = This situation corresponds to positions (, ) and (, ): S(, ) I/π = cosh π tan cosh π cos π tan cosh π cosh π cos π. () An eample is shown in Fig. 3. It is worth noting that the same method can be emploed for currents injected to thin-film samples of an polgonal shape. According to the Schwartz-Christoffel theorem an polgon can be mapped onto a halfplane. The general solution () on the (u, v) plane will hold. Therefore, knowing the function which realizes the needed transformation, one obtains S(, ). The onl phsical precondition for this method to work is the requirement of a small sample size on the scale of Pearl length λ /d, that allows one to reduce the problem to the Laplace equation for the stream function S. The method can be applied for more than two injection points or to etended injections, which require though different boundar conditions imposed on S.

3 φ j FIG. 3. The current distribution for the asmmetric injection points (., ) and (,.7). III. PHASE FIG.. The phase ϕ(, ) for j = at fied distances from the edge = as a function of the transverse coordinate for smmetric contacts at =.3 and =.7. The solid line is for =.5, the dashed line is for =.. We now note that the sheet current is epressed either in terms of the gauge invariant phase ϕ or via the stream function S: g = cφ π ϕ = curl Sz Λ This relation written in components shows that S(r) and (cφ /π Λ)ϕ(r) are the real and imaginar parts of an analtic function. It is eas to construct the phase on the (u, v) plane since i ln(u + iv) = tan (v/u) i ln u + v. Hence, the phase corresponding to the stream function of Eq. () obes: or cφ π Λ ϕ = I π ln u + (v v ) u + (v v ), (3) ϕ = I I ln u + (v v ) u + (v v ), I = cφ πλ. () The characteristic current I depends on the Pearl length. Thus, the phase is proportional to the reduced injected current j = I/I and a factor ϕ/j depending on the film geometr and injection positions. Substituting here u(, ) and v(, ) one obtains the phase as a function of (, ). Eamples of the phase near the edge = are shown in Fig.. Clearl, one can make the phase jump steeper b putting injection contacts closer. It is worth noting that the phase change due to injected currents can be used, e.g., to imitate properties of -π junctions, or in fact to have an phase shift b choosing properl the injected current. IV. JOSEPHSON CRITICAL CURRENT The total Josephson current through a rectangular patch (the shaded region in Fig. a) with the size along the ais is J(I/I ) J c = d d sin[ϕ(,, I/I ) + ϕ ], (5) where J c is the critical Josephson current densit (which in the following is set equal to ) and ϕ is an overall phase imposed b the transport current through the junction. To find the critical current, we maimize this relative to ϕ to obtain J c = A + B where A = B = d d d sin[ϕ(,, ; j)], A. Smmetric injection d cos[ϕ(, ; j)]. (6) Consider J c (j) for the rectangular junction of the width =. W, similar to the eperimental set up. The injection contacts are at the edge = and =.3 and =.7, i.e. the are smmetric relative to the strip middle = /. The current distribution for this case is given in Fig.. J c (j) evaluated with the help of Eq. (6) is shown in Fig. 5. We note that J c () is proportional to the junction width, since the Josephson critical current densit is constant in the absence of injected currents. Fig. 6 shows J c (j) in the same junction with contacts at =.8 and =.5 so that the are separated b =., ten times closer than in the previous eample. Comparing these plots we see that the first zero of J c at j.5 in the first graph whereas it is at j 5 in the second. We then conclude that zeros roughl scale as the

4 the last epression in powers of : πj ϕ(, ) = + O( )3. (8) cos π Since cos π is odd relative to the strip middle, and so is sin ϕ(, ), we have A =, and = B, see Eq. (6): Z Z d d cos[ϕ(, )] = Z d cos =. 3 5 j FIG. 5. (j) for the junction width =. and smmetric contacts at the edge =. =.3 and =.7 so that the distance between contacts is = j FIG. 6. (I/I ) for the junction width =. and smmetric contacts at the edge = and =.8 and =.5 so that the distance between contacts is =.. inverse of the contact separation, /. We also observe that maima of seem to be independent of contacts separation. It is shown below that these properties of (j) can be traced back to general epressions (5) and (6) for narrow junctions,, and small contacts separations. To this end, we note that being a solution of the Laplace equation for a half-infinite strip, the phase ϕ(, ) changes considerabl on distances of the order W = from the edge = where the contacts are placed. Hence, for narrow junctions with, one can set in the first approimation: (v v ) (v v ) (cos π cos π ) = j ln. (cos π cos π ) πj. cos π (9) Clearl, when no current is injected, () = as is seen in Figs. 5 and 6. The integral over can be done: 3 η, () = η p Fq,,, where p Fq is the generalized hpergeometric function and η = πj. This function is epressed in terms of Bessel and Struve functions so that it oscillates when η changes: π J (η) [J (η)h (η) J (η)h (η)]. () It is worth noting that of Eqs. (9) and () depends on the injected current j and the contacts separation onl via the product η = πj. In other words, the curves (j, ) for different contact separations can be rescaled to a universal curve (η). In particular, the zeros of (j) and the positions of its maima should scale as /. Unlike their positions, the absolute value of its n-th maimum is independent of the separation. The two first roots of (η) found numericall are. and.6. For =., we then obtain the two first zeros of (j): j = η /π. and j.6. For =., the two first zeros are at. and 6.. These are close to what we have at Figs. 5 and 6 although here we stud an approimate solution for whereas the figures are results of eact numerical integration. Let us now consider the magnetic field H applied parallel to the junction plane (, ). In general, the critical current should depend on both H and j, = (j, H). In the absence of injected currents, (, H) has a standard shape with maimum at H =. The presence of zeros of (j) for smmetric injection in zero field, has an important consequence. If the injected current jn is such that (jn ) =, application of the magnetic field will result in the pattern (jn, H) such that (jn, ) =, i.e., the curve (jn, H) will have zero at H = instead of the standard maimum. This situation is similar to the famous case of -π junction, however, here the injected currents cause a necessar phase shift. Precisel this situation has been seen in eperiment7,8 with smmetric injectors. ϕ(, ) ϕ(, ) = j ln B. (7) Here,, = ( ± )/, and for one can epand Asmmetric injection If the injection contacts are arranged asmmetricall as, e.g., at Fig. 3, the minima of (j) do not reach zeros as shown in Fig. 7. Without discussing a variet of

5 J c j FIG. 7. J c(j) corresponding to the current distribution of Fig. 3 for the injection points (., ) and (,.7). asmmetric injections, we note that J c (j) of all of them have the propert that their minima do not reach zeros. V. DISCUSSION In summar, we have shown that the current distribution in thin film samples small on the scale of the Pearl length Λ = λ /d can be found b solving the Laplace equation for the stream function under boundar conditions specified for injection sources at arbitrar points at sample edges. When this film constitutes one of the Josephson junction banks, the contribution of the phase associated with injected currents to the junction phase difference is proportional to the injected current j = πλi/cφ. Hence, the thinner the film (or the larger the Pearl Λ) the smaller injected currents are needed for the same effect upon the junction properties. The critical Josephson current J c (j), for certain (smmetric) infection geometries, has zeros, the position of which scales as the inverse distance between the injection points. If j n is one of these zeros, application of a field H parallel to the junction plane results in a pattern J c (H, j n ) with zero at H = instead of a standard maimum, the propert seen in eperiments. 7,8 The authors are grateful to E. Goldobin for sharing eperimental information and man helpful discussions. The Ames Laborator is supported b the Department of Energ, Office of Basic Energ Sciences, Division of Materials Sciences and Engineering under Contract No. DE-AC-7CH358. kogan@ameslab.gov mints@post.tau.ac.il L. N. Bulaevskii, V. V. Kuzii and A. A. Sobanin, Sol. St. Com, 5, 53 (978). C. C. Tsuei and J.R. Kirtle, Rev. Mod. Phs. 7, 969 (). 3 R. G. Mints, Phs. Rev. B57, R3 (998). T. Golod, A. Rdh, and V. M. Krasnov, Phs. Rev. Lett., 73 (). 5 V. G. Kogan, R. G. Mints, Phs. Rev. B89, 56 (). 6 V. G. Kogan, R. G. Mints, Phs. C Superc. 5, 58 (). 7 A. V. Ustinov, Appl. Phs. Lett. 8, 353 (). 8 E. Goldobin, A. Sterck, T. Gaber, D. Koelle, and R. Kleiner, Phs. Rev. Lett.9, 575 (). 9 A. Dewes, T. Gaber, D. Koelle, R. Kleiner, and E. Goldobin, Phs. Rev. Lett., 7 (8). P. M. Morse and H. Feshbach Methods of Theoretical Phsics, McGraw-Hill, 953. E. Goldobin, private communication.

Shapiro steps in Josephson junctions with alternating critical current density

PHYSICAL REVIEW B 76, 5458 27 Shapiro steps in Josephson junctions with alternating critical current density M. Moshe and R. G. Mints* School of Physics and Astronomy, Raymond and Beverly Sackler Faculty