Aug 6, 2010 - Î¾F = âD/h, while the oscillation period is given by. 2ÏÎ¾F . Here D is ..... cases of diffusive SNS junction with h = 011â14 and S...

2 downloads 17 Views 262KB Size

arXiv:1008.1251v1 [cond-mat.supr-con] 6 Aug 2010

Institute of Solid State Physics, Chernogolovka, Moscow reg., 142432 Russia (Dated: July 3, 2018) The Josephson current through a long s-wave superconductor/weak ferromagnet/s-wave superconductor weak link is studied theoretically in the regime of nonequilibrium spin-dependent occupation of electron states in the ferromagnetic intelayer. While under the considered nonequilibrium conditions the standard supercurrent, carried by the singlet part of current-carrying density of states, practically is not modified, the additional supercurrent flowing via the triplet part of the current-carrying density of states appears. Depending on voltage, controlling the particular form of spin-dependent nonequilibrium in the interlayer, this additional current can enhance or reduce the usual current of the singlet component and also switch the junction between 0- and π-states. PACS numbers: 74.45.+c, 74.50.+r

New states have been predicted and observed in Josephson weak links in the last years. One of them has common origin with the famous LOFF-state1,2 . This mesoscopic LOFF-state, which is induced in superconductor/weak ferromagnet/superconductor (SFS) Josephson junction, was predicted theoretically3,4 and observed experimentally5–8 . In this state Cooper pair acquires the total momentum 2Q or −2Q inside the ferromagnet as a response to the energy difference between the two spin directions. Here Q ∝ h/vf , where h is an exchange energy and vf is the Fermi velocity. Combination of the two possibilities results in the spatial oscillations of the condensate wave function Ψ(x) in the ferromagnet along the direction normal to the SF interface. Ψs (x) ∝ cos(2Qx) for the singlet Cooper pair9 . The same picture is also valid in the diffusive limit. The only thing we need to add is an extra decay of the condensate wave function due to scattering. In the regime h ≫ |∆|, where ∆ is a superconducting order parameter (OP) in the leads, the decay length is equal to the magnetic coherence length p ξF = D/h, while the oscillation period is given by 2πξF . Here D is the diffusion constant in the ferromagnet, ~ = 1 throughout the paper. The presence of an exchange field also leads to the formation of the triplet component of the condensate wave function in the interlayer. In the case of a homogeneous exchange field only the component with zero spin projection on the field direction Sz = 0 is induced. Combining the two pairs with the total momenta 2Q and −2Q into the triplet combination we see that the inplane (Sz = 0) triplet condensate wave function component Ψt (x) ∝ sin(2Qx), that is oscillates in space with the same period as the singlet one, but is shifted by π/2 with respect to it. Here we do not discuss the other triplet components with Sz = ±1, which are typically induced in case of inhomogeneous magnetization10 . The energy spectrum of the superconducting correlations is expressed in a so-called supercurrent- carrying density of states (SCDOS)11–14 . This quantity represents the density of states weighted by a factor proportional to the current that each state carries in a certain direction.

Under equilibrium conditions the supercurrent can be expressed via the SCDOS as13 Z j ∝ dεNj (ε) tanh ε/2T , (1) where ε stands for the quasiparticle energy, ϕ(ε) = tanh ε/2T is the equilibrium distribution function and Nj (ε) is SCDOS. In the presence of spin effects SCDOS becomes a matrix 2 × 2 in spin space and can be repreˆj = Nj,s + Nj,t σ, where σi are Pauli matrices sented as N in spin space. The scalar in spin space part of SCDOS Nj,s is referred to as the singlet part of SCDOS in the paper and the vector part Nj,t is referred to as the triplet part. Nj,t is directly proportional to the triplet part of the condensate wave function. It is well known that the spin supercurrent can not flow through the singlet superconducting leads. Therefore, Nj,t does not contribute to the supercurrent in equilibrium. Having in mind that the triplet part of SCDOS is even function of quasiparticle energy, one can directly see that this is indeed the case. Under nonequilibrium conditions one can change the value of the critical Josephson current through the junction and even realize the π-state by manipulating the quasiparticle distribution in an interlayer region. This effect was predicted theoretically11,12 and observed experimentally15–17 for a diffusive SNS junction. The point is that positive and negative parts of SCDOS give energy-dependent contributions to the supercurrent in the positive and negative direction. The size and direction of the total supercurrent depends therefore on the occupied fraction of these states, which is analogous to the occupation of the discrete Andreev bound states in a ballistic system. That is one can obtain negative Josephson current response to small phase differences and, hence, switch the system into π-state by creating an appropriate nonequilibrium quasiparticle distribution in the weak link region. The combination of the exchange field h ≪ |∆| and the spin independent nonequilibrium distribution function has been considered as well18 . Under these conditions the influence of the nonequilibrium distribution function is also consists in the redistribution

2 of the quasiparticles between the energy levels. However, it was shown that in this limit of small exchange fields the combined effect of the exhange field and the nonequilibrium distribution function is nontrivial. For instance, part of the field-suppressed supercurrent can be recovered by adjusting a voltage between additional electrodes, which controls the ditribution function. In the present paper we investigate the effect of nonequilibrium occupation of the supercurrent-carrying states on the Josephson current in SFS junction in the parameter range |∆| ≪ h ≪ εf , where εf is the Fermi energy of the ferromagnet. This regime is relevant to weak ferromagnetic alloys, which were used for the experimental realization of magnetic π-junctions. It is shown that if the distribution function becomes nonequilibrium and spin-dependent, the supercurrent carried by the SCDOS triplet component Nj,t in the ferromagnet is non-zero. The magnitude of this current cotribution jt can be of the same order or even larger than the current contribution js carried by the singlet component Nj,s . Due to the fact that the singlet and triplet components of the anomalous Green’s function have the same oscillation period but shifted in phase by π/2, jt can increase the usual supercurrent, carried by the singlet component of SCDOS, or weaken it, or even reverse the sign of the total supercurrent, thus switching between 0 and π-state. Experimentally the most probable way to realize the spindependent nonequilibrium in the interlayer is to apply a voltage to (or to pass the dissipative current through) a spin-active material. Then under appropriate conditions even quite small voltages should be enough to switch the system from 0 to π-state and vice versa. That is we consider another mechanism of supercurrent manipulation by creating a nonequilibrium quasiparticle distribution in the interlayer, which cannot be reduced to the redistribution of the quasiparticles between the energy levels. In principle, the both mechanisms can be realized in a junction simultaneously. However, in this particular study we assume h ≫ ∆ and the Thouless energy

j=

−d eRF

+∞ Z

−∞

εT h = D/d2 ≪ h, that is the interlayer length d ≫ ξF . As it is shown below, in this regime (and for the case of low-transparency SF interfaces) a spin-independent nonequilibrium distribution of quasiparticles in the ferromagnet practically does not affect the Josephson current, that is the described above mechanism of the critical current reversal by the spin-independent redistribution of supercurrent-carrying states population is irrelevant in this case. For high-transparency SF interfaces the both mechanisms contribute to supercurrent. For a quantitative analysis we use the formalism of quasiclassical Green-Keldysh functions in the diffusive limit19 . The fundamental quantity for diffusive transport is the momentum average of the quasiclassical Green’s function gˇ(x, ε) = hˇ g(pf , x, ε)ipf . It is a 8 × 8 matrix form in the product space of Keldysh, particle-hole and spin variables. Here x is the coordinate measured along the normal to the junction. The electric current should be calculated via Keldysh part of the quasiclassical Green’s function. For the plane diffusive junction the corresponding expression reads as follows −d j= eRF

+∞ Z

τ0 + τ3 dε K , (2) Tr (ˇ g (x, ε)∂ g ˇ (x, ε)) 4 x 8π 2 2

−∞

where e is the electron charge and RF stands for the reK sistance of the ferromagnetic region. (ˇ g(x, ε)∂x gˇ(x, ε)) is 4 × 4 Keldysh part of the corresponding combination of full Green’s functions. τi are Pauli matrices in particlehole space. It is convenient to express Keldysh part of the full Green’s function via the retarded and advanced components and the distribution function: gˇK = gˇR ϕˇ − ϕˇ ˇg A . Here the argument (x, ε) of all the functions is omitted for brevity. The distribution function is diagonal in particleˆ˜ 2 (τ0 − τ3 )/2. All the hole space: ϕˇ = ϕ(τ ˆ 0 + τ3 )/2 + σ2 ϕσ matrices denoted by ˆ are 2 × 2 matrices in spin space throughout the paper. In terms of the distribution function current (2) takes the form

i h dε ˆ ˆ ˆ ˆ g A ∂x gˆA + fˆA ∂x f˜A ) . ˆg A − fˆR ∂x ϕˆ˜f˜A + (ˆ g R ∂x gˆR + fˆR ∂x f˜R )ϕˆ − ϕ(ˆ Tr2 −π 2 ∂x ϕˆ − gˆR ∂x ϕˆ 2 8π

We assume that the direction of the exchange field h is spatially homogeneous and choose the quantization axis along the field. In this case the distribution function and the normal part gˆR,A of the Green’s function are diagonal matrices in spin space. The anomalous Green’s functions can be represented as fˆR,A = fˆdR,A iσ2 and ˆ ˆ ˆ f˜R,A = −iσ2 f˜dR,A , where fˆdR,A and f˜dR,A are diagonal in spin space.

(3)

The retarded and advanced Green’s functions are obtained by solving the Usadel equations19 supplemented with Kupriyanov-Lukichev boundary conditions at SF interfaces20 . It is worth to note that we can safely apply these boundary conditions to the problem of plane diffusive junction even for high enough dimensionless conductance g of the SF interface. This can be done in spite of the fact that they are the linear in transparency approximation of more general Nazarov boundary conditions21 .

3 The reason is that the effective number of interface channels N ∼ dy /l is large and a separate channel transparency T ∼ g(l/dy ) is usually considerably less than unity. Here dy is the junction width and l is the mean free path. Further, the condition d ≫ ξF allows us to find the solution analytically even for an arbitrary SF interface transparency and low temperature, that is in the parameter region, where the equations cannot be linearized with respect to the anomalous Green’s function. We start from the completely incoherent junction (that is, consider the left and right SF interfaces separately) and then calculate the corrections up to the first order of the small parameter exp [−d/ξF ] to the Green’s functions. Within this accuracy the anomalous Green’s functions in the vicinity of left and right SF interfaces (at x = ∓d/2) take the form i h R,A −iαχ/2 R,A iαχ/2 , fdσ = κiπ sinh ΘR,A e + 4Σ (x)e σ σ R,A R,A f˜dσ = −fdσ (χ → −χ) .

(4)

Here σ =↑, ↓ (or +1(−1) within Equations) is the electron spin projection on the quantization axis, κ = +1(−1) corresponds to the retarded (advanced) functions, α = +1(−1) in the vicinity of the left (right) SF interface and χ is the order parameter phase difference between the superconducting leads. The first term represents the anomalous Green’s function at the ferromagnetic side of the isolated SF boundary and does not enter the following results. So, we do not give it explicitly. The second term is the first order correction, originated from the anomalous Green’s function extended from the other SF interface: d −( − αx)(1 + iκσ)/ξF R,A R,A Σσ (x) = Kσ e 2 , (5) where KσR,A is determined by the equation RF ξF 2 (1 + iκσ)KσR,A (1 − KσR,A ) = sinh ΘR,A (1+ s 4Rg d i 2 4 R,A R,A 2 +KσR,A + KσR,A ) − cosh ΘR,A 4K (1 + K ) (6) . s σ σ

Rg stands for the resistance of each SF interface, which are supposed to be identical. sinh ΘR,A and cosh ΘR,A s s originate from the Green’s functions at the superconducting side of SF interfaces. We assume that p the parameter (RF ξs /Rg d)(σF /σs ) ≪ 1, where ξs = D/∆ is the superconducting coherence length in the leads, σF and σs stand for conductivities of ferromagnetic and superconducting materials, respectively. It allows us to neglect the suppression of the superconducting order parameter in the S leads near the interface and take the Green’s functions at the superconducting side of the boundaries to be equal to their bulk values. In this case −κiε cosh ΘR,A =p s 2 |∆| − (ε + κiδ)2 −κi|∆| sinh ΘR,A = p , (7) s 2 |∆| − (ε + κiδ)2

where δ is a positive infinitesimal. The SCDOS, entering the current (3) takes the form ˆ ˆ Njσ (ε) = gˆR ∂x gˆR + fˆR ∂x f˜R − gˆA ∂x gˆA − fˆA ∂x f˜A = σ 8π 2 sin χ d R = Im ΣR . (8) σ (x = −α ) sinh Θs σF Rg 2 The nonequilibrium distribution function in the interlayer is proposed to be created by applying a voltage along the y direction between two additional electrodes Nb and Nt , which are attached to the central part of the interlayer. It is supposed that the conductances of Nb F and Nt F interfaces are spin-dependent and equal to gb σ and gt σ , respectively. The voltage Vt σ (Vb σ ) between the superconducting leads and Nt (Nb ) electrode can also be spin-dependent. It can be realized, for example, by attaching one (or both) of the electrodes Nb or Nt to a strong ferromagnet and applying a voltage between the other one and the ferromagnet. In order to obtain the distribution function up to the first order of the parameter exp [−d/ξF ] we solve the kinetic equation for the distribution function, which is derived from the Keldysh part of the Usadel equation. The boundary conditions to the kinetic equation are also obtained from the Keldysh part of the general KypriyanovLukichev boundary conditions. Further it is assumed that |eVt,b ↑,↓ | < |∆|. Under this condition the part of the current associated with the first three terms in Eq. (3) is to the only determined by the first order correction ΣR,A σ anomalous Green’s function and the distribution function (0) ϕσ , calculated up to the zero order of the parameter exp [−d/ξF ]: (0) (0) ϕ˜σ − ϕσ ˆ −π 2 ∂x ϕˆ − gˆR ∂x ϕˆ ˆg A − fˆR ∂x ϕˆ˜f˜A = Njσ , 2 σ (9) where Njσ (ε) is expressed by Eq. (8). The distribution function ϕ(0) does not depend on x. So, it is convenient to calculate it in the middle of the interlayer, where disregarding the parameter exp [−d/ξF ] means disregarding the superconducting proximity effect. Then, by considering Nb /F/Nt junction and applying Kypriyanov-Lukichev boundary conditions Nb /F and F/Nt interfaces, we come to the following expression for ϕ(0) (inelastic scattering processes are not taken into account):

ϕ(0) σ =

tσ tanh ε−eV gt σ (σF + dy gb σ ) 2T + σF (gt σ + gb σ ) + 2dy gt σ gb σ

+ (0)

bσ tanh ε−eV gb σ (σF + dy gt σ ) 2T , σF (gt σ + gb σ ) + 2dy gt σ gb σ

(0)

(10)

ϕ˜σ is connected to ϕσ by the symmetry relation (0) (0) ϕ˜↑,↓ (ε) = −ϕ↓,↑ (−ε). We focus on the case gtσ ≪ σF /dy or gbσ ≪ σF /dy , when the y-dependence of the distribution function ϕ(0) can be disregarded.

4 Substituting Eqs. (8)-(9) into Eq. (3) we find that the Josephson current takes the form Z i dε X h (0) −d (ϕσ + ϕ˜(0) , (11) j= σ )Njσ 2 2eRF 8π σ (0)

where Njσ (ε) is expressed by Eq. (8) and ϕσ should be taken from Eq. (10). Elecrtical current (11) can be divided into two parts: j = js + jt , Z i dε h (0) −d (0) (ϕ + ϕ ˜ )N , js = js,c sin χ = j,s 0 0 eRF 8π 2 Z i −d dε h (0) (0) jt = jt,c sin χ = (12) , (ϕ + ϕ ˜ )N j,t z z eRF 8π 2 where ϕ0 = (ϕ↑ + ϕ↓ )/2, ϕz = (ϕ↑ − ϕ↓ )/2, ϕ˜0,z = (ϕ˜↑ ± ϕ˜↓ )/2 = ∓ϕ0,z (−ε), Nj,s = (Nj↑ + Nj↓ )/2 is the singlet part of SCDOS and Nj,t = (Nj↑ − Nj↓ )/2 is the z-component of the SCDOS triplet part (the other components equal to zero for the considered case of homogeneous magnetization). It is seen from Eq. (12) that Nj,t gives rise to the additional contribution to the spinless electrical current if the quasiparticle distribution is spin-dependent. While Eqs. (11)-(12) are valid for arbitrary SF interface transparency, at first we concentrate on the discussion of the tunnel limit g˜ ≡ RF ξF /Rg d ≪ 1, where Eq. (6) can be easily solved and the integral over energy can be calculated analytically. For the analytical analysis we choose the most simple form for the distribution function Eq. (10) by setting gt σ → 0 and Vb ↓ = 0. Then (0) (0) ϕ↑ = tanh[(ε − eVb ↑ )/2T ], while ϕ↓ = tanh[ε/2T ]. As it is demonstrated below, the results corresponding to another set of parameters, which determines the particular form of the distribution function satisfying Eq. (10), qualitatively do not differ from that ones represented here. At low temperature T ≪ |eVb↑ | we obtain for the Josephson current the following result √ π RF ξF sin χ d −d/ξF + )+ |∆|e j= 2π cos( 2 4eRg d ξF 4 |∆| + eVb↑ d π 1 √ log sin( + ) . (13) |∆| − eVb↑ ξF 4 2 The first term in Eq. (13) represents the part of the supercurrent js carried by the singlet component of SCDOS. Under the conditions T ≪ |eVt,b | and |eVt,b ↑,↓ | < |∆| it is not affected by the fact that the distribution function is nonequilibrium, as can be seen from Eq. (13). The reason is that for a long junction and h ≫ |∆| in the tunnel limit the singlet part of SCDOS is concentrated in the narrow energy intervals around ε = ±|∆|, as it is illustrated in panel (a) of Fig. 1. This is opposed to the cases of diffusive SNS junction with h = 011–14 and SFS junction with h ≪ ∆18 , where the singlet part of SCDOS is finite and exhibits nontrivial energy dependence in the subgap energy region |ε| < |∆|. Under the conditions

T ≪ |eVt,b | and |eVt,b ↑,↓ | < |∆| the distribution function (0)

(0)

ϕ0 + ϕ˜0 ≈ (sgn[ε − eVb↑ ] + sgn[ε + eVb↑ ] + 2sgn[ε])/2, which enters the expression for js [Eq.(12)], practically does not differ from its equilibrium value for energy intervals around ε = ±|∆|. That is, under the considered conditions the widely discussed in the literature mechanism of supercurrent manipulation by nonequilibrium redistribution of quasiparticles between energy levels11,12,15–18 does not contribute. The second term jt is caused by the triplet component of SCDOS and vanishes in the equilibrium Vb↑ = 0. As it is seen in panel (b) of Fig. 1, the triplet part of SCDOS Nj,t is an even function of energy and has finite value in the subgap energy region. So, multiplying it by the distribution function ϕ(0) ˜(0) z +ϕ z ≈ (sgn[ε − eVb↑ ] − sgn[ε + eVb↑ ])/2 , (14) one obtains current contribution jt . The absolute value of this contribution is roughly proportional to Vb↑ for small enough values of this parameter and increases sharply when Vb↑ approaches |∆|. This behavior is a consequence of two facts: (i) distribution function (14) is a constant within energy interval [−eVb↑ , eVb↑ ] and vanishes outside it, (ii) the triplet part of SCDOS has a particular shape shown in panel (b) of Fig. 1. To calculate the Josephson current for the case of arbitrary transparency of SF interfaces one needs to solve Eq. (6) numerically and make use of Eq. (11). The resulting curves as functions of eVb↑ are plotted in Fig. 2. Panel (a) shows the current for low enough dimensionless conductance of SF interface g˜ = 0.1, while panel (b) represents the case of highly transparent interface g˜ = 3.0. Different curves correspond to different lengths d of the junction. In dependence on Vb ↑ the current can be enhanced or reduced with respect to its value at Vb ↑ = 0. If the length of the equilibrium junction is not far from the 0-π transition, then small enough voltage can switch between the states. Separate plots of the current contributions js and jt are represented in panel (a) of Fig. 3 for the low-transparency junction with g˜ = 0.1 and in panel (b) for the high-transparency junction with g˜ = 3.0. It is seen, that for low-transparency junction the current and its separate contributions js and jt behave just as described by the tunnel limit discussed above. While the tunnel limit qualitatively captures the essential physics for the high-transparency junction as well, there are some new features, which are discussed below. The singlet part of SCDOS for the case of hightransparency junction is plotted in panel (c) of Fig. 1. It is clearly seen that as distinct from the limit of tunnel junction, it is not only concentrated around superconducting gap edges, but is finite in the whole subgap region. It results in the sensitivity of js to the fact that (0) (0) the distribution function ϕ0 + ϕ˜0 is nonequilibrium. In other words, upon increasing of the SF-interface conductance the mechanism of the current controlling based on the redistribution of the quasiparticles between energy

0.02

g˜=0.1 d=3.4ξF

-1

1

0

0.01 0 -0.01 -2

-1

1

0

0.1 0.05 0 -0.05 -0.1

ε/∆ 2

g˜=3.0 d=3.6ξF

(c)

ε/∆ 2 g˜=0.1 d=3.4ξF

(b)

Nj,s (ε)

0.04 0.03 (a) 0.02 0.01 0 -0.01 -0.02 -0.03 -2

-1

-2

Nj,t(ε)

Nj,t(ε)

Nj,s (ε)

5

0.08 0.06 0.04 0.02 0 -0.02

1

0

ε/∆ 2 g˜=3.0 d=3.6ξF

(d)

-1

-2

1

0

ε/∆ 2

FIG. 1. The singlet and triplet parts of SCDOS in dependence on the quasiparticle energy. The left column corresponds to the low-transparency limit g˜ = 0.1, while the right column represents the high-transparency case g˜ = 3.0. 0.004

0.003

g˜ = 0.1 3.4 3.7

0.001

3.96

0

4.3

0.002

0.001 0

(b)

-0.5

0

g˜ = 3.0

0.5 eVb↑ /∆ 1 3.6

0.02

0.01

4.46 4.8 5.2

-0.5

0

0.5 eVb↑ /∆ 1

-0.5

0

1

0.5

eVb,↑ /∆

0.03

4.1

0

-0.002 -1

jc |e|Rg /∆

5.0

-1 jc |e|Rg /∆

g˜=0.1 d=3.4 ξF

-0.001

-0.001

-1

(a)

0.003

0.002

0.03

jc |e|Rg /∆

jc |e|Rg /∆

0.004

(a)

(b)

0.02

0.01 0 g˜=3.0 d=3.6 ξF

-0.01 -1

-0.5

0

0.5

1 eVb,↑ /∆

FIG. 2. The critical Josephson current as a function of eVb↑ /∆. The different curves correspond to different lengths d of the junction, which are measured in units of ξF .

levels starts to contribute. It is seen from panel (b) of Fig. 3 that the difference js,c (Vb↑ ) − js,c (Vb↑ = 0) considerably grows when |eVb↑ | approaches |∆|. It is worth to note that the absolute value of js,c is always reduced by this mechanism. This fact leads to the reduction of the absolute value of the high-transparency total supercurrent (see Fig. 2(b)) at eVb↑ → ∆. This behavior should be compared to the tunnel limit, where the dependence of js,c on Vb↑ is negligible and absolute value of the total supercurrent only increases for eVb↑ → ∆ due to growing contribution of jt,c . Now we discuss the dependence of the obtained results

FIG. 3. Separate plots of the current contributions js,c (dashed line) and jt,c (dotted line) together with the total critical current js,c + jt,c (solid line) as a function of eVb↑ /∆. Panel (a) represents the low-transparency case and panel (b) corresponds to high-transparency case.

to the particular shape of the distribution function. It is illustrated in Fig. 4. Four different examples of the (0) (0) distribution function z-component ϕz + ϕ˜z , satisfying Eq. (10) are shown in the right panel. Due to symmetry relations this combination is always an even function of energy. The scalar part of the distribution function is not plotted in the figure because it does not influence the supercurrent in the tunnel limit. The corresponding plots of the critical current in the low-transparency limit are represented in the left panel. It can be con-

6

0.004

0.003

g˜=0.1 d=3.4ξF 3

0.002

0.001

cluded that, while the quantitative value of the critical current depends on the particular choice of the distribution function, this choice has no qualitative effect on the supercurrent behavior, as it was already pointed before.

1

ϕz + ϕ˜z

jc |e|Rg /∆

0.005

2

1 3

4 2

0

4

-0.001 -1

-0.5

0

0.5

1 eVb,↑ /∆

-2

-1

0

1

ε/∆

2

FIG. 4. Dependence of the critical current on the particular shape of the distribution function. The right panel represents four different examples of the distribution function vector part. The curves marked by the numbers from 1 to 4 correspond to the following sets of parameters in Eq. (10): 1. gtσ → 0, Vb↓ = 0; 2. gt↑ , gb↓ → 0, eVt↓ = 0.8∆; 3. gt↑ , gb↓ → 0, eVt↓ = −0.8∆ and 4. gt↑ = 0.1σF /dy , gt↓ = 4σF /dy , gb↑ = 3σF /dy , gb↓ = 0.2σF /dy , Vb↓ = 0.5Vb↑ , eVt↑ = 0.35∆, eVt↓ = 0.7∆. The respective critical current plots in dependence on eVb↑ , calculated for the same sets 1-4 of the distribution function parameters are shown in the left panel. The dashed vertical line marks the position eVb↑ = 0.2∆, for which the distribution functions are calculated.

1

2 3

4

5

6

7

8

9

A.I. Larkin and Yu.N. Ovchinnikov, Sov. Phys. JETP 20, 762 (1965) [Zh. Eksp. Teor. Fiz. 47, 1136 (1964)]. P. Fulde and R.A. Ferrel, Phys.Rev. 135, A550 (1964). A. I. Buzdin, L. N. Bulaevsky, and S. V. Panyukov, JETP Lett. 35, 178 (1982) [Pis’ma Zh. Eksp. Teor. Fiz. 35, 147 (1982)]. A. I. Buzdin, B. Bujicic, and M. Yu. Kupriyanov, Sov. Phys. JETP 74, 124 (1992) [Zh. Eksp. Teor. Fiz. 101, 231 (1992)]. V. V. Ryazanov, V. A. Oboznov, A. Yu. Rusanov, A. V. Veretennikov, A. A. Golubov, and J. Aarts, Phys. Rev. Lett. 86, 2427 (2001). T. Kontos, M. Aprili, J. Lesueur, F. Genet, B. Stephanidis, R. Boursier, Phys. Rev. Lett. 89, 137007 (2002). Y. Blum, A. Tsukernik, M. Karpovski, and A. Palevski, Phys. Rev. Lett. 89, 187004 (2002). W. Guichard, M.Aprili, O. Bourgeois, T. Kontos, J. Lesueur, and P. Gandit, Phys. Rev. Lett. 90, 167001 (2003). E. A. Demler, G. B. Arnold, M. R. Beasley, Phys. Rev. B 55, 15174 (1997).

In conclusion, we have studied the effects of nonequilibrium spin-dependent electron distribution in a weakly ferromagnetic interlayer on the Josephson current through SFS junction. It is shown that the nonequilibrium spindependent electron distribution gives rise to the supercurrent carried by the triplet component of SCDOS. Depending on voltage, controlling the particular form of spin-dependent nonequilibrium in the interlayer, this additional current can enhance or reduce the usual current of the singlet component and also switch the junction between 0- and π-states. We gratefully acknowledge discussions with V.V. Ryazanov, A.S. Melnikov and M.A. Silaev. The support by RFBR Grant 09-02-00779 and the programs of Physical Science Division of RAS is acknowledged. A.M.B. was also supported by the Russian Science Support Foundation.

10

11 12

13 14

15

16

17

18

19 20

21

F.S. Bergeret, A.F. Volkov, and K.B. Efetov, Rev. Mod. Phys. 77, 1321 (2005). A.F. Volkov, Phys. Rev. Lett. 74, 4730 (1995). F.K. Wilhelm, G. Sch¨ on, and A.D. Zaikin, Phys. Rev. Lett. 81, 1682 (1998). S.-K. Yip, Phys. Rev. B 58, 5803 (1998). T.T. Heikkil¨ a, J. S¨ arkk¨ a, and F.K. Wilhelm, Phys. Rev. B 66, 184513 (2002). J.J.A. Baselmans, A.F. Morpurgo, B.J. van Wees, and T.M. Klapwijk, Nature (London) 397, 43 (1999). J. Huang, F. Pierre, T.T. Heikkil¨ a, F.K. Wilhelm, and N.O. Birge, Phys. Rev. B 66, 020507 (2002). M.S. Crosser, J. Huang, F. Pierre, P. Virtanen, T.T. Heikkil¨ a, F.K. Wilhelm, and N.O. Birge, Phys. Rev. B 77, 014528 (2008). T.T. Heikkil¨ a, F.K. Wilhelm, and G. Sch¨ on, Europhys. Lett. 51, 434 (2000). K.D. Usadel, Phys.Rev.Lett. 25, 507 (1970). M.Yu. Kupriyanov and V.F. Lukichev, Sov. Phys. JETP 67, 1163 (1988). Yu.V. Nazarov, Superlattices Microstruct. 25, 1221 (1999).