Chapter 10: Superconductivity: Bardeen, Cooper, & Schrieffer April 26, 2017
Chapter 10: Superconductivity: Bardeen, Cooper, & Schrieffer April 26, 2017
Chapter 10: Superconductivity: Bardeen, Cooper, & Schrieffer April 26, 2017
Contents
1 Introduction 3
1.1 Evidence of a Phase Transition . . . . . . . . . . . . . . . . . . . . . 3
1.2 Meissner Effect . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
3 Cooper Pairing 11
3.1 The Retarded Pairing Potential . . . . . . . . . . . . . . . . . . . . . 11
3.2 Scattering of Cooper Pairs . . . . . . . . . . . . . . . . . . . . . . . . 13
3.3 The Cooper Instability of the Fermi Sea . . . . . . . . . . . . . . . . 15
1
6 BCS ⇒ Superconducting Phenomenology 33
9 Tunnel Junctions 43
10 Unconventional Superconductors 48
10.1 D-wave Superconductors . . . . . . . . . . . . . . . . . . . . . . . . . 49
10.2 Triplet Superconductors . . . . . . . . . . . . . . . . . . . . . . . . . 53
10.3 Odd-frequency Superconductors . . . . . . . . . . . . . . . . . . . . . 54
2
1 Introduction
3
C (J/mol°K)
Cn ∼ γT
CS
T T
c
Figure 1: The specific heat of a superconductor CS and and normal metal Cn . Below
the transition, the superconductor specific heat shows activated behavior, as if there is
a minimum energy for thermal excitations.
Cs ∼ e−β∆ (1)
4
within the bulk of a superconductor. This is fundamentally dif-
ferent than an ideal conductor, for which Ḃ = 0 since for any
closed path
Superconductor
I Z Z
1 ∂B
0 = IR = V = E ·dl = ∇×E ·dS = − ·dS , (2)
S c S ∂t
or, since S and C are arbitrary
1
0 = − Ḃ · S ⇒ Ḃ = 0 (3)
c
Thus, for an ideal conductor, it matters if it is field cooled or
zero field cooled. Where as for a superconductor, regardless
of the external field and its history, if T < Tc, then B = 0
inside the bulk. This effect, which uniquely distinguishes an
ideal conductor from a superconductor, is called the Meissner
effect.
5
Ideal Conductor
Zero-Field Cooled Field Cooled
T > Tc T > Tc
B=0 B≠0
T < Tc T < Tc
B=0 B≠0
T<T T < Tc
c
B≠0 B=0
Figure 3: For an ideal conductor, flux penetration in the ground state depends on
whether the sample was cooled in a field through the transition.
µ−1
B = µH = 0 ⇒ µ = 0 M = χH = H (4)
4π
1
χSC = − (5)
4π
Ie., the measured χ, Fig. 4, in a superconducting metal is very
large and negative (diamagnetic). This can also be interpreted
6
χ ∝ D(E )
Pauli F
Tc
0
T
js
χ M
∼
∼ ∼
∼
H
-1
4π
7
H
Normal
Hc
S.C.
Tc T
∂js e2ns
= E (First London Eqn.) (8)
∂t m
8
Then, using the Maxwell equation
1 ∂B m ∂js 1 ∂B
∇×E =− ⇒ ∇ × + =0 (9)
c ∂t nse2 ∂t c ∂t
or
∂ m 1
∇ × js + B = 0 (10)
∂t nse2 c
This described the behavior of an ideal conductor (for which
ρ = 0), but not the Meissner effect. To describe this, the
constant of integration must be chosen to be zero. Then
nse2
∇ × js = − B (Second London Eqn.) (11)
mc
m
or defining λL = ns e2
, the London Equations become
B ∂js
= −λL∇ × js E = λL (12)
c ∂t
4π
If we now apply the Maxwell equation ∇×H = c j ⇒ ∇×B =
4π
c µj then we get
4π 4πµ
∇ × (∇ × B) = µ∇ × j = − 2 B (13)
c c λL
and
1 4πµ
∇ × (∇ × j) = − ∇×B=− 2 j (14)
λLc c λL
or since ∇ · B = 0, ∇ · j = 1c ∂ρ
∂t = 0 and ∇ × (∇ × a) =
SC
^ ^
∂Bx
B j ∝∇×B∝z×x
s ∂z
y
j z
Figure 6: A superconducting slab in an external field. The field penetrates into the
q
mc2
slab a distance ΛL = 4πne 2µ .
10
3 Cooper Pairing
region of + +
8 + + + +
vF ∼ 10 cm/s positive charge
attracts a second
+
electron
Figure 7: Origin of the retarded attractive potential. Electrons at the Fermi surface
travel with a high velocity vF . As they pass through the lattice (left), the positive ions
respond slowly. By the time they have reached their maximum excursion, the first
electron is far away, leaving behind a region of positive charge which attracts a second
electron.
11
tice deforms slowly in the time scale of the electron. It reaches
its maximum deformation at a time τ ∼ 2π
ωD ∼ 10−13s after the
electron has passed. In this time the first electron has traveled
◦
8 cm −13
∼ vF τ ∼ 10 s · 10 s ∼ 1000 A. The positive charge of
the lattice deformation can then attract another electron with-
out feeling the Coulomb repulsion of the first electron. Due
to retardation, the electron-electron Coulomb repulsion may be
neglected!
The net effect of the phonons is then to create an attrac-
tive interaction which tends to pair time-reversed quasiparticle
states. They form an antisymmetric spin singlet so that the
k↑
e
ξ ∼ 1000Α°
e
- k↓
Figure 8: To take full advantage of the attractive potential illustrated in Fig. 7, the
spatial part of the electronic pair wave function is symmetric and hence nodeless. To
obey the Pauli principle, the spin part must then be antisymmetric or a singlet.
2
Ek ∼ k
k’1
k’ k’
1
2
ω
D
k’ k2
2 k2 k1
Figure 9: Pair states scattered by the exchange of phonons are restricted to a narrow
scattering shell of width h̄ωD around the Fermi surface.
13
this scattering process is also conserved
scattering shell
k1 -k
2
K
Figure 10: If the pair has a finite center of mass momentum, so that k1 + k2 = K,
then there are few states which it can scatter into through the exchange of a phonon.
14
3.3 The Cooper Instability of the Fermi Sea
Now consider these two electrons above the Fermi surface. They
will obey the Schroedinger equation.
h̄2 2
− (∇1 + ∇22)ψ(r1r2) + V (r1r2)ψ(r1r2) = ( + 2EF )ψ(r1r2)
2m
(19)
If V = 0, then = 0, and
1 ik1·r1 1 ik2·r2 1 ik(r1−r2)
ψV =0 = e e = e , (20)
L3/2 L3/2 L3
where we assume that k1 = −k2 = k. For small V, we will
perturb around the V = 0 state, so that
1 X
ψ(r1r2) = 3 g(k)eik·(r1−r2) (21)
L
k
15
The Schroedinger equations may be converted to a k-space
equation by multiplying it by
Z
1 3 −ik0 · r
d r e ⇒ S.E. (24)
L3
so that
h̄2k 2 1 X
g(k) + 3 g(k0)Vkk0 = ( + 2EF )g(k) (25)
m L 0
k
where Z
−k0 )·r 3
Vkk0 = V (r)e−i(k dr (26)
Z EF +h̄ωD
dE
1 = V0 Z(EF ) (32)
EF 2E − − 2EF
1 − 2h̄ωD
1 = V0Z(EF ) ln (33)
2
2h̄ωD V0
= 2/(V Z(E ))
' −2h̄ωD e−2/(V0Z(EF )) < 0, as → 0
1−e 0 F EF
(34)
18
pair states labeled by k
or
|ψk i = uk |0ik + vk |1ik (39)
Where σk+ and σk−, describe the creation and anhialation of the
19
state (k ↑, −k ↓)
0 1
σk+ = 12 (σk1 + iσk2 ) = (42)
0 0
0 0
σk− = 21 (σk1 − iσk2 ) = (43)
1 0
0 1
Of course σk+ =
1 0
k
20
Thus the reduction of the potential energy is given by hφBCS |V | φBCS i
(
V0 Y X
hφBCS |V | φBCS i = − 3 (up h0| + vp h1|) σk+σk−0
L p 0
kk
Y
up0 |0ip0 + vp0 |1ip0 (47)
0
p
21
∂WBCS V0 X
= 0 = −4ξk cos θk sin θk − 3 cos 2θk sin 2θk0 (52)
∂θk L 0
k
1 V0 X
ξk tan 2θk = − 3 sin 2θk0 (53)
2L 0
k
p
Conventionally, one introduces the parameters Ek = ξk2 + ∆2, ∆ =
V0 P V0 P
L3
u v
k k k = L3 k cos θk sin θk . Then we get
∆
ξk tan 2θk = −∆ ⇒ 2uk vk = sin 2θk = (54)
Ek
−ξk
cos 2θk = = cos2 θk − sin2 θk = vk2 − u2k = 2vk2 − 1 (55)
Ek
!
1 −ξk 1 ξk
wk = vk2 = 1− = 1− p 2 (56)
2 Ek 2 ξk + ∆2
∆ 2 1 ξk
If we now make these substitutions 2uk vk = E , vk = 2 1 − E
k k
wk = v 2 T=0
k clearly kinetic
energy increases
1
2 2
h k
0 ξ k = -E +
F 2m
22
into WBCS , then we get
L3 2
X ξk
WBCS = ξk 1 − − ∆. (57)
Ek V0
k
23
see this consider
L3∆2
X 1 ξk
WBCS = 2ξk 1− − (61)
2 Ek V0
k
↓ Lots of algebra (See I&L)
X
WBCS = − 2Ek vk4 (62)
e
-k′↓
2
w = v2 = 1 vk′ =0
k′ k′
p
Figure 12: Breaking a pair requires an energy 2 ξk2 + ∆2 ≥ 2∆
X X q
∆E = − 2vk4 Ek + 2vk4 Ek = 2Ek0 = 2 ξk20 + ∆2 (63)
k6=k 0 k
h̄2 k 02
Then since ξk0 = 2m − EF , the smallest such excitation is just
∆Emin = 2∆ (64)
24
This is the minimum energy required to break a pair, or create
an excitation in the BCS ground state. It is what is measured
by the specific heat C ∼ e−β2∆ for T < Tc.
Now consider some experiment which adds a single electron,
or perhaps a few unpaired electrons, to a superconductor (ie
tunneling). This additional electron cannot find a partner for
superconductor normal
metal
Figure 13:
h̄2 k 02
For ξk2 ∆, Ek = ξk = 2m − EF , which is just the energy of
a normal metal state. Thus for energies well above the gap, the
normal metal continuum is recovered for unpaired electrons.
25
To calculate the density of unpaired electron states, recall
that the density of states was determined by counting k-states.
These are unaffected by any phase transition. Thus it must be
that the number of states in d3k is equal.
3
kz d k
ky
3
π
L
k
x
26
E
Density of additional
∆ electron states only!
1 Ds Dn
Figure 15:
should be calculated.
V0 X V0 X V0 X ∆
∆ = 3 sin θk cos θk = 3 uk vk = 3 (68)
L L L 2Ek
k k k
1 V0 X ∆
∆ = (69)
2 L3
p
k
ξk2 + ∆2
Convert this to sum over energy states (at T = 0 all states with
2 2
ξ < 0 are occupied since ξk = h̄2mk − EF ).
Z h̄ωD
V0 Z(EF + ξ)dξ
∆ = ∆ p (70)
2 −h̄ωD ξ 2 + ∆2
Z h̄ωD
1 dξ
= p (71)
V0Z(EF ) 0 ξ 2 + ∆2
1 h̄ω D
= sinh−1 (72)
V0Z(EF ) ∆
27
For small ∆,
h̄ωD 1
∼ e V0 Z(EF )
(73)
∆
1
− V Z(E
∆ ' h̄ωD e 0 F) (74)
sinh x
∼ ex
Figure 16:
∆E ∼ 2∆e−β2∆ T Tc (75)
28
∂∆E ∂β ∆2 −β2∆
C∼ ∼ 2e (76)
∂β ∂T T
10
hω
microwave
hω
hω = 2∆
B
Figure 17: If B > Bc or h̄ω > 2∆, then absorption reduces the intensity to the
normal-state value I = In . For B = 0 the microwave intensity within the cavity is
large so long as h̄ω < 2∆
29
tensity is high I = Is. On the other hand, if h̄ω > 2∆ ,or
B > Bc, then the intensity falls in the cavity I = In due to
absorbs ion by the walls.
Note that this also allows us to measure ∆ as a function of
T. At T = Tc, ∆ = 0, since thermal excitations reduce the
number of Cooper pairs and increase the number of unpaired
electrons, which obey Fermi-statistics. The size of (Eqn. 71) is
k′↑ e
kT ∼ 2∆
e
-k′↓
Figure 18:
30
p
1 − 2f ξ2 + ∆2 + EF , T . Thus for T 6= 0
Z h̄ωD
1 dξ n p o
= p 1 − 2f ξ 2 + ∆ 2 + EF , T
V0Z(EF ) 0 ξ 2 + ∆2
(77)
p
Note that as ξ 2 + ∆2 ≥ 0, when β → ∞ we recover the
T = 0 result.
This equation may be solved for ∆(T ) and for Tc. To find Tc
In Pb
∆(T) Sn
∆(0)
Real SC data (reflectivity)
T/Tc
1
Figure 19: The evolution of the gap (as measured by reflectivity) as a function of tem-
perature. The BCS approximation is in reasonably good agreement with experiment.
T
consider this equation as Tc → 1, the first solution to the gap
equation, with ∆ = 0+, occurs at T = Tc. Here
Z h̄ωD
1 dξ ξ
= tanh (78)
V0Z(EF ) 0 ξ 2kB Tc
which may be solved numerically to yield
1.14h̄ωD
1 = V0Z(EF ) ln (79)
kB Tc
31
kB Tc = 1.14h̄ωD e−1/{V0Z(EF )} (80)
Table 1: Note that the value 2.15 for ∆(0)/kB Tc for Pb is higher than BCS predicts.
Such systems are labeled strong coupling superconductors and are better described by
the Eliashberg-Migdal theory.
32
1
Thus Tc ∼ M − 2 . This has been confirmed for most normal
superconductors, and is considered a ”smoking gun” for phonon
mediated superconductivity.
(r - r0 )/ΛL
j = j0 e
H H
• ⊗
S r0 j
0
dl
Λ
L
Figure 20: Integration contour within a long thick superconducting wire perpendicular
to a circulating magnetic field. The field only penetrates into the wire a distance ΛL .
4π
∇×H= j (83)
c
33
along the contour shown in Fig. 20.
Z Z Z
4π
∇ × HdS = H · dl = j · ds (84)
c
4π
2πr0H = 2πr0ΛLj0 (85)
c
If j0 = jc (jc is the critical current), then
4π
Hc =ΛLjc (86)
c
Since both Hc and jc ∝ ∆, they will share the temperature-
dependence of ∆.
At T = 0, we could also get an expression for Hc by noting
that, since the superconducting state excludes all flux,
1 1 2
(W n − W BCS ) = H (87)
L3 8π c
However, since we have earlier
1 1
3
(W n − W BCS ) = N (0)∆2, (88)
L 2
we get
p
Hc = 2∆ πN (0) (89)
c
We can use this, and the relation derived above jc = 4πΛL Hc ,
34
to get a (properly derived) relationship for jc.
c p
jc = 2∆ πN (0) (90)
4πΛL
However, for most metals
n
N (0) ' (91)
EF
s
mc2
ΛL = (92)
4πne2µ
taking µ = 1
r s
c 4πne2 πn2m √ ne
jc = 2∆ = 2∆ (93)
4π mc2 h̄2kF2 h̄kF
This gives a similar result to what Ibach and Lüth get, but
for a completely different reason. Their argument is similar to
one originally proposed by Landau. Imagine that you have a
fluid which must flow around an obstacle of mass M . From the
perspective of the fluid, this is the same as an obstacle moving
in it. Suppose the obstacle makes an excitation of energy and
momentum p in the fluid, then
E0 = E − P0 = P − p (94)
Figure 21: A superconducting fluid which must flow around an obstacle of mass M .
From the perspective of the fluid, this is the same as an obstacle, with a velocity equal
and opposite the fluids, moving in it.
E′
(a) (b)
P P′
M M
E p
Figure 22: A large mass M moving with momentum P in a superfluid (a), creates an
excitation (b) of the fluid of energy and momentum p
P 02 P2 P·p p2
− =− + = E0 − E = (95)
2M 2M M 2M
P θ
v = P/M
P′
Figure 23:
pP cos θ p2
= − (96)
M 2M
p2
= pv cos θ − (97)
2M
36
If M → ∞ (a defect in the tube which carries the fluid could
have essentially an infinite mass) then
= v cos θ (98)
p
Then since cos θ ≤ 1
v ≥ (99)
p
Thus, if there is some minimum ,then there is also a mini-
mum velocity below which such excitations of the fluid cannot
happen. For the superconductor
min 2∆
vc = = (100)
p 2h̄kF
Or
ne
jc = envc = ∆ (101)
h̄kF
This is the same relation as we obtained with the previous
√
thermodynamic argument (within a factor 2). However, the
former argument is more proper, since it would apply even for
gapless superconductors, and it takes into account the fact that
the S.C. state is a collective phenomena ie., a minuet, not a
waltz of electric pairs.
37
7 Coherence of the Superconductor ⇒ Meisner
effects
3
4πn ξcp
∼ 108 (105)
3 2
other pairs have their center of mass.
38
Figure 24: Many electron pairs fall within the volume of a Cooper wavefunction.
This leads to a degree of correlation between the pairs and to rigidity of the pair
wavefunction.
ξ coh > ξ cp
Figure 25:
2e
j=− {ψp∗ψ ∗ + ψ ∗pψ} (106)
4m
39
where pair mass = 2m and pair charge = −2e.
2e
p = −ih̄∇ − A (107)
c
A current, or a CM momentum K, modifies the single pair state
1 X
ψ(r1, r2) = 3 g(k)eiK· (r1+r2)/2eik· (r1−r2) (108)
L
k
φ = K · (R1 + R2 + · · · ) (111)
∇ = ∇R + ∇r ≈ ∇R (112)
Thus
2e X 2eA
js ≈ Φ∗BCS −ih̄∇Rν + ΦBCS
4m ν c
∗
2eA
+ΦBCS ih̄∇Rν + Φ∗BCS (113)
c
40
or
( )
2e 2 4eA 2
X
js = − |Φ(0)| + 2h̄ |Φ(0)| ∇R ν φ (114)
2m c ν
ne2
∇×j=− B (116)
mc
which is the second London equation which as we saw in Sec. 2
leads to the Meissner effect. Thus the second London equation
can only be derived from the BCS theory by assuming that the
BCS state is spatially homogeneous.
41
distance ΛL from the surface).
e2ns eh̄ns X
js = − A− ∇Rν φ (117)
mc 2m ν
e2ns
Z Z Z
eh̄ns X
◦js · dl = − ◦A · dl − ◦∇Rν φ · dl (118)
ms 2m ν
Presumably the phase of the BCS state ΦBCS = eiφΦ(0) is
superconducting loop
C
X
X
X
X
X
X X X
B
X X X X ΛL
X X
X
X X
Figure 26: Magnetic flux penetrating a superconducting loop is quantized. This may
be seen by integrating Eq. 114 along a contour within the superconducting bulk (a
distance ΛL from the surface).
single valued, so
XZ
∇Rν φ · dl = 2πN N ∈Z (119)
ν
Also since the path l may be taken inside the superconductor
by a depth of more than ΛL, where js = 0, we have that
Z
js · dl = 0 (120)
42
so
e2ns e2ns
Z Z
eh̄ns
− A · dl = − B · ds = 2N π (121)
ms ms 2m
Ie., the flux in the loop is quantized.
9 Tunnel Junctions
V0
Figure 27:
0 0
ψa = A1eikx + B1e−ikx ψb = A2eik x + B2e−ik x
ψc = B3e−ikx (122)
43
These are solutions to the S.E. if
√
2mE
k = in a & c (123)
h̄
p
2m(E − V0)
k0 = in b (124)
h̄
The coefficients are determined by the BC of continuity of ψ
and ψ 0 at the barriers x = 0 and x = d. If we take B3 = 1
and E < V0, so that
p
2m(E − V0)
k 0 = iκ = (125)
h̄
then, the probability of having a particle tunnel from left to
right is
( 2 2 )−1
2
|B3| 1 1 1 k κ 1 k κ
Pl→r ∝ = = − − + + cosh 2κd
|B1|2 |B1|2 2 8 κ k 8 κ k
(126)
For large κd
44
−2
k κ
Pl→r ∝ 8 + e−2κd (127)
κ k
−2 ( p )
k κ 2d 2m(V0 − E)
∝ 8 + exp − (128)
κ k h̄
eV ∼ ∆ (129)
4kB Tc
We know that 2∆
kB Tc ∼ 4, ∆ ∼ 2 ∼ 10◦K. However typical
metallic densities of states have features on the scale of electron-
volts ∼ 104◦K. Thus, on this energy scale we may approximate
45
S I N
E
eV
X N(E)
Figure 28: Electrons cannot tunnel accross the barrier since no unoccupied states are
available on the left with correspond in energy to occupied states on the right (and
vice-versa). However, the application of an appropriate bias voltage will promote the
state on the right in energy, inducing a current.
46
EF
Figure 29: If eV= 0, but there is a small overlap of occupied and unoccupied states on
the left and right sides, then there still will be no current due to a balance of particle
hopping.
suppressed. Then
Z
I ∼ P Nr (EF ) df ( − eV )Nl () (132)
and
∂f ( − eV )
Z
dI
∼ P Nr (EF ) d Nl () (133)
dV ∂V
∂f
∼ eδ( − eV − EF ) (T EF ) (134)
∂V
dI
' P Nr (EF )Nl (eV + EF ) (135)
dV
dI
Thus the low temperature differential conductance dV is a mea-
sure of the superconducting density of states.
47
dI
I
dV
V
∆/e V ∆/e
10 Unconventional Superconductors
48
supererconductors, to very unusual ideas, that have not yet
been observed, such as odd-frequency superconductors. In each
case, I will assume that the superconducting state is formed via
Cooper pairing of fermions, so that the order parameter must
remain odd under a product of symmetry operations, such as
parity, time reversal, spin, etc. as summarized in Tab. 2.
Type spin symmetry inversion symmetry time reversal symmetry
50
We assume that the pairing interaction is strongly peaked at
the antiferromagnetic ordering wavenumbers k = (±π, ±π), or
V (k) = −V0δ(k − (±π, ±π)) where V0 > 0 and is only finite
at energies near the fermi surface. Since the pairing is due to
magnetic correlations, the width of the scattering shell is now
roughly J which is assumed to be J EF . In this case the
gap equation becomes
0 1 X ∆(k)V (k − k0)
∆(k ) = 3 p (136)
2L ξk2 + ∆(k)2
k
52
10.2 Triplet Superconductors
53
This triplet state should have a number of experimental con-
sequences. Perhaps the most obvious is that, like the d-wave
superconductors, the pairing should be very sensitive to diorder,
at least disorder with a mean-free path that is shorter than the
pairing length. I.e., the stronger the pairing, the less sensitive
the state is do disorder.
54