Skip to main content

Full text of "Andreev Reflection In Ferromagnet-Superconductor Junctions"

See other formats

Andreev reflection in ferromagnet— superconductor junctions 

M. J. M. de Jong'^'* and C. W. J. Beenakker^ 
(a) Philips Research Laboratories, 5656 A A Eindhoven, The Netherlands 
(h) Instituut-Lorentz, University of Leiden, 2300 RA Leiden, The Netherlands 
(Submitted September 12, 1994 — |cond-mat/9410014| ) 


The transport properties of a ferromagnet-super conductor (FS) junction are 
studied in a scattering formulation. Andreev reflection at the FS interface is 
strongly affected by the exchange interaction in the ferromagnet. The con- 
ductance Gps of a ballistic point contact between F and S can be both larger 
or smaller than the value Gfn with the superconductor in the normal state, 
depending on the ratio of the exchange and Fermi energies. If the ferromagnet 
contains a tunnel barrier (I) , the conductance Gfifs exhibits resonances which 
do not vanish in linear response — in contrast to the Tomasch oscillations for 
non- ferromagnetic materials. 

PACS numbers: 74.80.Fp, 72.10.Bg, 74.50.+r 


Typeset using REVT^ 

Electrons in a metal can not penetrate into a superconductor if their excitation energy 
with respect to the Fermi level is below the superconducting gap A. Still, a current may flow 
through a normal-metal-superconductor (NS) junction in response to a small applied voltage 
V < A/e, by means of a scattering process known as Andreev reflection |l|]: An electron in 
the normal metal is retroreflected at the NS interface as a hole and a Cooper pair is carried 
away in the superconductor. Andreev reflection near the Fermi level conserves energy and 
momentum but does not conserve spin — in the sense that the incoming electron and the 
Andreev reflected hole occupy opposite spin bands. This is irrelevant for materials with 
spin-rotation symmetry, as is the case for normal metals. However, the change in spin band 
associated with Andreev reflection may cause an anomaly in the conductance of (metallic) 
ferromagnet-superconductor (FS) junctions, because the spin-up and the spin-down band in 
the ferromagnet are different. This paper contains a theoretical study of Andreev reflection 
in FS junctions. We use a scattering approach based on the Bogoliubov-de Gennes equation 
to study the transport properties for zero temperature and small V {eV -C A). We will 
concentrate on two distinct effects, which we think are experimentally observable. First, 
due to the change in spin band there is no complete Andreev reflection at the FS interface. 
This has a clear influence on the conductance and the shot-noise power of clean FS point 
contacts. Second, the different spin-up and spin-down wavevector at the Fermi level may 
lead to quantum-interference effects. This shows up in the linear-response conductance of 
FIFS junctions, where the ferromagnet contains an insulating tunnel barrier (I). 

In the past, FS junctions with an insulating layer between the ferromagnet and the 
superconductor have been used in spin-dependent tunneling experiments 0. There the 
emphasis was on the voltage scale eV ^ A and Andreev reflection did not play a role. 
Tunneling through S-Fi-S junctions, where Fi is a magnetic insulator, has been studied 
both experimentally and theoretically [^Q. In addition, there has been theoretical work 
on the Josephson effect in SFS junctions An experimental investigation of the boundary 
resistance of sputtered SFS sandwiches has also been reported j^. The importance of phase 
coherence was demonstrated in a recent experiment ||^, in which the effect of a remote 
superconducting island on the conductance of a ferromagnet was observed. We do not 
know of any previous theoretical work on the influence of Andreev reflection on the sub-gap 
conductance of an FS junction. 

In order to clarify the effects we are aiming at, let us flrst give an intuitive and sim- 
ple description of the conductance through a ballistic FS point contact. A ferromagnet is 
contacted through a small area with a superconductor. The transverse dimensions of the 
contact area are much smaller than the mean free path and the interface is clean, so that 
the conductance is completely determined by the scattering processes that are intrinsic to 
the FS interface. In a semiclassical approximation all scattering channels (transverse modes 
in the point contact at the Fermi level) are fully transmitted, when the superconductor is in 
the normal state. Let N^{Nf^) be the number of up(down)-spin channels, so that N-^ > N^. 
At zero temperature, the spin channels do not mix and the conductance is given by the 
Landauer formula 

Gfn = jiNi + iVr) • (1) 

In the superconducting state, the spin-down electrons of all the A^^ channels are Andreev 
reflected into spin-up holes. They give a double contribution to the conductance since 2e is 


transferred at each Andreev reflection. However, only a fraction Ni/N^ of the channels 
can be Andreev reflected, because the density of states in the spin-down band is smaller 
than in the spin-up band. Therefore, the resulting conductance is 

G^s = j{2N^ + 2j;^N^)=4-N^. (2) 

Comparison of Eqs. (|1|) and shows that Gps may be either larger or smaller than Gfn 
depending on the ratio N^/N^. If A^|/A^| < 1/3 then G-ps < Gfn, and vice versa. This 
qualitative argument can be substantiated by an explicit calculation, as we now show. 

For the conduction electrons inside the ferromagnet we apply the Stoner model, using 
an effective one-electron Hamiltonian with an exchange interaction. The effect of the fer- 
romagnet on the superconductor is twofold. First, there is the influence of the exchange 
interaction on states near the interface. This will be fully taken into account. Second, there 
is the effect of the magnetic field due to the magnetization of the ferromagnet. Since this 
field — which is typically a factor thousand smaller than the exchange field — does not 
break spin-rotation symmetry it will be neglected for simplicity. Note that in typical layered 
structures the magnetization is parallel to the FS interface, so that it has no influence on 
the superconductor at all. 

Transport through NS junctions has successfully been investigated through the 
Bogoliubov-de Gennes equation |TD|-0. Here, we adopt this approach for an FS junc- 
tion. In the absence of spin-flip scattering in the ferromagnet, the Bogoliubov-de Gennes 
equation breaks up into two independent matrix equations, one for the up-electron, down- 
hole quasiparticle wavefunction {u^,vi) and another one for {ui,v^). Each matrix equation 

has the form 14 

Ho-h A 

A* -{Ho + h) 


Here, e is the quasiparticle energy measured from the Fermi energy 

p^/2m + V — Ep is the single-particle Hamiltonian, with V^(r) the potential energy, /i(r) is 
the exchange energy, and A(r) is the pair potential. For simplicity, it is assumed that the 
ferromagnet and the superconductor have identical T^o- For comparison with experiment, 
our model can easily be extended to include differences in effective mass and band bottom. 
We adopt the usual step-function model for the pair potential |1T0|-[13| and do the same for 

the exchange energy Defining the FS interface at x = with S at a; > 0, we have 

A(r) = /\Q{x) and /i(r) = hQQ{—x), with Q{x) the unit step-function. 

A scattering formula for the linear-response conductance of an NS junction is given by 
Takane and Ebisawa [|12[. Application to the FS case is straightforward. 

Gfs = 2- ^ Trr^^^r,,^,,, (4) 

where the matrix v^^^f.^ contains the reflection amplitudes from incoming electron modes 
with spin a to outgoing hole modes with spin a (opposite to a) evaluated at the Fermi 
level (e = 0). We first consider a ballistic point contact. We assume that the dimensions 
of the contact are much greater than the Fermi wavelength, as is appropriate for a metal. 


so that quantization effects can be neglected. The number Ni of minority spin modes in 
the point contact (with area fl) is = Nq{1 — ho/Ep), with Nq = kpfl/An the number 
of modes per spin for a non-ferromagnetic {ho = 0) contact of equal area. The reflection 
matrices for this case can be evaluated by matching the bulk solutions for the ferromagnet 
and for the superconductor at the interface. An incoming electron from the ferromagnet is 
either normally reflected as an electron of the same spin or Andreev reflected as a hole with 
the opposite spin. (Transmission into the superconductor is not possible at e = 0.) The 
reflection matrices are diagonal, with elements 

^ _ _ k^ki - 

k'\ki + q^ 

where the longitudinal wavevectors km) in the ferromagnet and q in the superconductor are 
defined in terms of the energy En of the n-th transverse mode by 

q=^{2m/h^){EF-En), (6a) 
= ^{2m/h^){EF-En + ho) , (6b) 
ki = ^{2m/n^){EF-En-ho) . (6c) 

In the above expressions terms of order A/Ep are neglected ||T5[. Note that |reeP+|r/ieP = 1 

as required from quasiparticle conservation. It follows from Eq. (|g) that a clean FS junction 
does not exhibit complete Andreev reflection, in contrast to the NS case. This is due to the 
potential step the particle passes when being Andreev reflected to the opposite spin band. 

Because of the large number of modes the trace in Eq. (^) can be replaced by an inte- 
gration, which can be evaluated analytically. The result is 


Gfs = 4: — No 

h 15?7 


where t] = ho/Ep. The conductance is plotted in Fig. |l], and compared with the semiclassical 
estimate from Eq. (|), which turns out to be quite accurate. Since + N'^ = 2Nq one has 
from Eq. (|I|) Gfn > Gps if ho > 0A7Ep, or equivalently NjNi^ < 0.36. 

Further information on the Andreev reflection at the FS interface can be obtained from 
the shot-noise power P of the junction. Shot noise is the time-dependent fluctuation in the 
current due to the discreteness of the charges. For uncorrelated electron transmission, one 
has the maximal noise power of a Poisson process Ppoisson = 2eJ, with / the mean current. 

On the one hand, correlations due to the Pauli principle reduce P below Ppoisson [16,17]. On 

the other hand. Cooper-pair transport across an NS junction has been shown to manifest 

itself as a doubling of the maximal noise power [16,18]. We apply the general result of Ref. 
]|T8| to the FS junction 



Substitution of Eq. (|5bD into Eq. (H) yields the shot-noise power of a balhstic point contact, 
plotted in Fig. 0. The shot noise increases from complete suppression for a non-ferromagnetic 
{Hq = 0) junction to twice the Poisson noise for a half-metallic ferromagnet {Hq = Ep). The 
initial increase is slow, indicating that the modes undergo nearly complete Andreev re- 
flection. However, for higher exchange energies the Andreev reflection probability decreases 
in favour of the normal reflection probability. This is manifested by the increase in the 
shot-noise power. 

The second system we consider is an FIFS junction which contains a planar tunnel barrier 
(I) at X = —L. The barrier is modeled by a channel- and spin-independent transmission 
probability F G [0, 1]. The matrix r[- ^^r^- in Eq. (||) is diagonal, with elements 



F^lr^el' { l + 2p2cos(xT-Xi)+p' 

+ 2reep(l + p^) (cos XT + COS Xl) 

+ 2r,V[l + cos(xT + Xi)]}"\ (9) 

where p = a/1 — F and Xu = '^k^^L. Eq. describes resonant Andreev reflection: Due 
to the different wavevector of up electrons and down holes, | function 
of Xt S'lid Xl between F^, the value for a two-particle tunneling process, and 1 for full 
resonance. The conductance Gfifs is evaluated by substitution of Eq. into Eq. (^). It 
is depicted in Fig. |^ as a function of L for = 0.2Ep and F = 0.1. The resonances have 
a dominant period 5L = 7ihvF/2hQ{= biikp^ in Fig. 0), which is caused by the simplest 
round-trip containing two Andreev reflections and two barrier reflections. Superimposed 
one sees oscillations with smaller period, caused by longer trajectories in which also normal 
reflections at the FS interface occur. This becomes clear when we calculate Gps with ree set 
to zero, which is also shown in Fig. |^. For large L, Gfifs approaches the classical (i.e. all 
interferences are neglected) value 4(e^//i)A^|F/(2 — F). The oscillations in Fig. |^ are distinct 
from the Tomasch oscillations known to occur in the non-linear differential conductance of 

NINS junctions |T^. There, quasi-bound states arise because electron and hole wavevectors 
disperse if e > 0. However, in linear response Gnins = 4(e^//i)Ai'oF^/(2 — F)^, independent 
of L In the ferromagnetic junction the resonances do not vanish in linear response, in 
contrast to the Tomasch oscillations. The quasi-bound states at the Fermi level are a direct 
consequence of the change in spin band upon Andreev reflection. 

We believe that both phenomena are experimentally accessible. The FS point contact can 
be constructed according to the nanofabrication technique of Ref. [^. The FIFS junction 
can be made by growing a wedge-shaped layer of ferromagnet on a superconducting substrate 
and then depositing a thin oxide layer. This allows a measurement of Gfifs for different 
values of L. It is not necessary for the contact area to be small, so that no nanofabrication 
techniques are needed. (Note, that in order to observe the resonances due to the quasi- 
bound states it is not essential that the contact on top of the barrier is a ferromagnet.) 
To estimate the effect of disorder (growth imperfections and impurities) on the resonances, 
we have numerically calculated Gfifs for a disordered ferromagnet between the barrier and 

the FS interface. The computations are similar to the NS case treated in Ref. I^ll. The 
disordered region is modeled by a tight-binding Hamiltonian on a square lattice with a 
random impurity potential at each site. (For computational efficiency the geometry is two- 
dimensional, but this makes no qualitative difference.) The matrix r/i^eo- is obtained by 


combining the scattering matrix of the disordered region with the reflection coefficients for 
the FS interface (|]). We then calculate Gps through Eq. (^. The result for various disorder 
strengths is shown in Fig. ^ For the clean case we recognize a behavior similar to Fig. 
H. Adding some disorder removes the small-period oscillations but preserves the dominant 
oscillations. Only quite a strong disorder (for the top curve kp x bulk mean free path ^ 9) 
is able to smooth away the resonances. 

In summary, we have shown that the transport properties of ferromagnet-superconductor 
junctions are qualitatively different from the non-ferromagnetic case, because the Andreev 
reflection is modified by the exchange interaction in the ferromagnet. Two illustrative exam- 
ples have been given: For a ballistic FS point contact it is found that the conductance can 
be both larger or smaller than the normal-state value and for an FIFS junction containing 
a tunnel barrier conductance resonances are predicted to occur in linear response. 

We are especially grateful to H. van Houten for suggesting the problem treated in this 
paper. Furthermore, we thank P. J. Kelly and C. M. Schep for useful discussions. This 
research was supported by the Dutch Science Foundation NWO/FOM. 



[1] A. F. Andreev, Zh. Eksp. Teor. Fiz. 46, 1823 (1964) [Sov. Phys. JETP 19, 1228 (1964)]. 
[2] A recent review is: R. Meservey and P. M. Tedrow, Phys. Rep. 238, 173 (1994). 
[3] F. Stageberg, R. Cantor, A. M. Goldman, and G. B. Arnold, Phys. Rev. B 32, 3292 

[4] M. J. DeWeert and G. B. Arnold, Phys. Rev. Lett. 55, 1522 (1985); Phys. Rev. B 39, 
11307 (1989). 

[5] S. V. Kuplevakhskii and I. I. Fal'ko, Fiz. Met. Metalloved. 71, 68 (1991) [Phys. Met. 

Metallogr. 71, 65 (1991)]. 
[6] L. N. Bulaevskii, A. I. Buzdin, and S. V. Panjukov, Solid State Commun. 44, 539 (1982). 
[7] S. V. Kuplevakhskii and I. I. Fal'ko, Fiz. Met. Metalloved. 62, 13 (1986) [Phys. Met. 

Metallogr. 62, 8 (1986)]; Pis'ma Zh. Eksp. Teor. Fiz. 52, 957 (1990) [JETP Lett. 52, 

340 (1990)]. 

[8] C. Fierz, S.-F. Lee, J. Bass, W. P. Pratt, Jr., and P. A. Schroeder, J. Phys. Condens. 

Matter 2, 9701 (1990). 
[9] V. T. Petrashov, V. N. Antonov, S. V. Maksimov, and R. Sh. Shaikhaidarov, Pis'ma 
Zh. Eksp. Teor. Fiz. 59, 523 (1994) [JETP Lett. 59, 551 (1994)]. 
[10] G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. Rev. B 25, 4515 (1982). 
[11] C. J. Lambert, J. Phys. Condens. Matter 3, 6579 (1991). 
[12] Y. Takane and H. Ebisawa, J. Phys. Soc. Jpn. 61, 1685 (1992). 
[13] C. W. J. Beenakker, Phys. Rev. B 46, 12841 (1992). 

[14] P. G. de Gennes, Superconductivity of Metals and Alloys (Benjamin, New York, 1966). 
[15] This is the Andreev approximation One can easily go beyond it by including terms of 

order A/Ep in Eqs. (|^) and (^. We have checked that this has only a small influence on 

our final results. In fact, the larger ho, the more accurate is the Andreev approximation. 
[16] V. A. Khlus, Zh. Eksp. Teor. Fiz. 93, 2179 (1987) [Sov. Phys. JETP 66, 1243 (1987)]. 
[17] G. B. Lesovik, Pis'ma Zh. Eksp. Teor. Fiz. 49, 513 (1989) [JETP Lett. 49, 592 (1989)]; 

M. Biittiker, Phys. Rev. Lett. 65, 2901 (1990). 
[18] M. J. M. de Jong and C. W. J. Beenakker, Phys. Rev. B 49, 16070 (1994). 
[19] W. J. Tomasch, Phys. Rev. Lett. 15, 672 (1965); 16, 16 (1966); W. L. McMillan and P. 

W. Anderson, Phys. Rev. Lett. 16, 85 (1966); A. Hahn, Phys. Rev. B 31, 2816 (1985). 
[20] P. A. M. Holweg, J. A. Kokkedee, J. Caro, A. H. Verbruggen, S. Radelaar, A. G. M. 

Jansen, and P. Wyder, Phys. Rev. Lett. 67, 2549 (1991). 
[21] I. K. Marmorkos, C. W. J. Beenakker, and R. A. Jalabert, Phys. Rev. B 48, 2811 (1993). 



FIG. 1. The conductance Gps (full curves) and the shot-noise power Pfs (dashed) of a ballistic 
point contact in a ferromagnet-superconductor junction (see inset), as a function of the exchange 
energy /iq. The thick line represents the exact result (^) for Gpsi the thin line the estimation (^). 

FIG. 2. The conductance Gfifs of a clean FIFS junction containing a planar tunnel barrier 
(transparency F) on the ferromagnetic side, as a function of the separation L from the interface 
(see inset). The thick solid line is computed from Eq. (^) for F = 0.1, /iq = Q.2Ep. For the thin 
line normal reflection at the FS interface is neglected {r^e = 0). The dashed line is the classical 
large- -L limit. 

FIG. 3. Numerical calculation of the effect of disorder in the ferromagnet on the oscillations 
shown in Fig. |2| for a clean junction. The disordered region is modeled by a L x IF square lattice 
(lattice constant o) with random on-site disorder (uniformly distributed between ±[7/2). The 
width W = 101a is fixed and the length L is varied on the horizontal axis. The results shown are 
for Ep = fi^ /2'mo? , /iq = 0.2Ef, F = 0.1, and for various U. For each disorder strength U the bulk 
mean free path i is given. Thick lines belong to one realization of disorder, thin to an average over 
20 realizations.