Modern Aspects of Superconductivity - Theory of Superconductivity
Modern Aspects of Superconductivity - Theory of Superconductivity
Superconductivity
Theory of Superconductivity
Sergei Kruchinin
Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine
For photocopying of material in this volume, please pay a copying fee through the Copyright Clearance Center,
Inc., 222 Rosewood Drive, Danvers, MA 01923, USA. In this case permission to photocopy is not required from
the publisher.
ISBN-13 978-981-4261-60-9
ISBN-10 981-4261-60-2
Printed in Singapore.
“The single reason for our inability to treat the problems of
superconductivity consists in the absence of a sufficient
imagination.”
— Richard P. Feynman
Preface
Studies of superconductivity theory are among the most fruitful and promis-
ing trends in the theoretical physics of condensed matter, since supercon-
ductivity remains one of the most interesting research areas in physics.
The goal of this book is to give a representation of certain modern
aspects of superconductivity. We discuss important aspects of the theory of
superconductivity, such as the nature of high-Tc superconductivity, two-gap
superconductivity, room-temperature superconductivity, mesoscopic super-
conductivity, the pairing state and the mechanism of cuprate high-Tc super-
conductivity.
In Chap. 1, we consider the field-theoretical method of superconductivity
and discuss the basic idea of superconductivity and the elaboration of the
Ginzburg–Landau and Bardeen–Cooper–Schrieffer theories in the frame of
many-particle quantum field theory.
In Chap. 2, we consider the structures of high-Tc superconductors, phase
diagrams and the problem of pseudogaps, and analyze the mechanisms of
superconductivity. We present general arguments regarding the pairing sym-
metry in cuprate superconductors and investigate their thermodynamical
properties within the spin-fluctuation mechanism of superconductivity, by
using the method of functional integrals.
Chapter 3 concentrates on two-band and multiband superconductiv-
ity. We consider the physical properties of the superconductor MgB2 and
use our two-band model to explain the two coupled superconductor gaps
of MgB2 . To study the effect of the increasing Tc in MgB2 , we use the
renormalization-group approach and phase diagrams. In the field of super-
conductivity we meet the problem-maximum, which consists in the creation
of room-temperature superconductors. We consider this problem in our book
and make some recommendations on the search for these superconductors.
vii
viii Modern Aspects of Superconductivity: Theory of Superconductivity
Preface vii
1. Theory of Superconductivity 1
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Spinors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.2.1 Spinor . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2.2 Noether theorem and Nambu–Goldstone
theorem . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.3 Propagator . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.3.1 Hamiltonian . . . . . . . . . . . . . . . . . . . . . . 11
1.4 Noninteracting . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.5 Interacting . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.5.1 Unrestricted Hartree–Fock (HF) . . . . . . . . . . . 19
1.5.2 Gap equation for superconductivity . . . . . . . . . 20
1.6 Illustrative Example, Critical Temperature . . . . . . . . . 21
1.6.1 Bond alternation . . . . . . . . . . . . . . . . . . . . 21
1.6.2 Deformation energy . . . . . . . . . . . . . . . . . . 23
1.6.3 Polyacene, gap equation, critical temperature . . . . 24
1.6.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . 27
1.7 Linear Response Magnetic Resonance
in Normal and Superconducting Species;
Spin–lattice Relaxation Time . . . . . . . . . . . . . . . . . 28
1.7.1 Introduction . . . . . . . . . . . . . . . . . . . . . . 28
1.7.2 T1 in NMR . . . . . . . . . . . . . . . . . . . . . . . 31
1.7.3 Theory with Green’s function . . . . . . . . . . . . 32
1.7.4 Noninteracting . . . . . . . . . . . . . . . . . . . . . 35
1.7.5 Interacting; normal . . . . . . . . . . . . . . . . . . 36
ix
x Modern Aspects of Superconductivity: Theory of Superconductivity
Bibliography 203
Index 215
CHAPTER 1
Theory of Superconductivity
1.1 Introduction
In the preface to the book Superconductivity, edited by Parks,1 one of the
editors says: “During the preparation of this treatise one of the authors com-
mented that it would be the last nail in the coffin [of superconductivity].”
Some specialists in superconductivity told us that we would hardly be able
to find articles useful for our future investigations, except for Anderson’s
comments at the end of that book. Further, we learned that Anderson was
pessimistic about further advance of superconductivity; for example, high-
temperature or room-temperature superconductivity was most unlikely.
However, against his expectation, the discovery of high-temperature super-
conductivity, due to Bednorz and Müller,2 is astonishing. The traditional
Bardeen–Cooper–Schrieffer (BCS) theory3 has failed to explain the mecha-
nism of such superconductivity. Anderson5 proposed a new idea called the
resonating valence bond (RVB) theory or the t−J model (the term t implies
the transfer integral, and J the electron correlation). We have never known
his theory to be successful.
A comment by Feynman, found in his book Statistical Mechanics,4 says
that it will take almost 50 years for the problem of superconductivity to be
reduced to that of explaining the gap. Following the BCS theory, we will
explain the gap, and the theory is essentially correct, but we believe that
it needs to be made obviously correct. As it stands now, there are a few
seemingly loose ends to tie up.
The theory of superconductivity seems to be founded on the London pos-
tulate.7 Associated with the gauge transformation, the conserved current is
i
j = − (φ∗ ∇φ − φ∇φ∗ ) − e|φ|2 A. (1.1)
2
1
2 Modern Aspects of Superconductivity: Theory of Superconductivity
The current due to the first term is called the paramagnetic current; and that
due to the second, the diamagnetic current. In the superconducting state,
the first term on the right-hand side changes very slightly, and sometimes
the wave function is quite rigid, that only the diamagnetic current survives.
In this respect, the superconductor is the perfect diamagnetic substance.
The current is dominated by
j = −k2 A, (1.2)
∇ × B = j. (1.3)
∇2 B = k 2 B (1.4)
or
Bx = B0 e−kx . (1.5)
It is important to note that |φ|2 in Eq. (1.1) is very large, i.e. it is the classic
scale quantity, so that the magnetic field in Eq. (1.5) damps very rapidly.
This is the Meissner effect. On the book cover, we see the very spectacular
experiment demonstrating the Meissner effect — a permanent magnet hov-
ers above a superconducting plate. We know of similar phenomena — the
screening of the Coulomb interaction, or quark confinement.
The boson model of Cooper pairs is considerably successful. Equa-
tion (1.3) yields ∇2 A = k2 A, which is
∂ μ Aμ = −k2 Aμ (1.6)
The theory is handled as the phase transition. The Lagrangian for the
superconducting state is postulated as
1 1 e∗ A 2 h2
2
Fs = F0 + a|Ψ| + b|Ψ| + 4
−i∇ + + 8π , (1.7)
2 2m∗ c
where ∗ indicates the quantities in question referred to the superconductor
(×2). Note that there are |Ψ|2 and |Ψ|4 terms in the potential parts. These
drive the system to spontaneous symmetry breaking and lead to a phase
transition for the suitable choice of constants a and b. This is now called
the Higgs mechanism.8 Certainly, the GL treatment is a few years ahead of
the Nambu–Goldstone suggestion.
The GL theory is known as the macroscopic quantum mechanism, and
in this sense, the big Ψ is sometimes amusingly called the cat’s Ψ.
A very instructive presentation of the macroscopic quantum theory is
found in The Feynman Lectures on Physics. Vol. III, Chap. 21. Various
topics there, such as the Josephson junction, are quite readable.
The microscopic theory was prepared by Bardeen, Cooper and
Schrieffer.3 However, as a preliminary discussion, we present the Bogoliubov
treatment. The superconductivity is a kind of many-electron problems. The
most general Hamiltonian should be
1
H = dxψ+ (x)h(x)ψ(x) + dx dx ψ+ (x)ψ(x)v(x, x )ψ+ (x )ψ(x ).
2
(1.8)
Since the algebra of electrons is the spinor, the terms other than the
above identically vanish. In other words, the three- or four-body interac-
tions are useless. First, we specify the spin indices, and next the plane-wave
representation for the spatial parts. Equation (1.8) is simply written as
∗ ∗
H = k (a+
k + ak ) + vk,−k ak a−k a−k ak . (1.9)
The simplification or the mean-field approximation for the two-body part is
twofold — say,
Δk,−k akα akα a+
−kβ a−kβ
+
(1.10)
Δk,−k a+
kα akβ akα a−kβ + c.c.
+
The latter looks curious, since such an expectation value, akα a−kβ , vanishes
identically. Against this common sense, Bogoliubov put the Hamiltonian
∗ ∗ ∗
k + ak ) + Δk ak a−k + Δk a−k ak .
H = k (a+ (1.11)
4 Modern Aspects of Superconductivity: Theory of Superconductivity
with
u2k − vk2 = 1. (1.13)
The spirit of the Bogoliubov transformation is to mix operators ck↑ and
c+
−k↓ , which are different in spin. The quasiparticle yields the new ground
state, so that the particle pair or the hole pair arises near the chemical
potential. The stabilization energy thus obtained is called the gap energy
Δk .3, 10
We now follow the BCS microscopic treatment. The Green function has
been effectively used by employing the Nambu spinor. This makes the uni-
fied treatment of normal and superconducting states possible. However,
the temperature Green function (Matsubara function) is used from the
beginning.
1.2 Spinors
We start with the general many-electron Hamiltonian not restricted to the
BCS Hamiltonian. The BCS state responsible for the superconducting state
is easily recognized from the formal description. The simplest method of the
quantum chemistry should be the Hückel theory. This consists of the single
energy matrix,
βrs = dx h(x)ρrs (x), (1.14)
Any 2 × 2 matrix is expanded in the Pauli spin matrices together with the
unit matrix:
1 0 1 0 0 1 0 −i
σ =
0
, σ = 3
, σ =
1
, σ =
2
.
0 1 0 −1 1 0 i 0
(1.16)
However, we employ other combinations:
↑1 1 0
σ = (σ0 + σ3 ) = ,
2 0 0
↓ 1 0 0 0
σ = (σ − σ ) =
3
,
2 0 1
(1.17)
1 1 0 1
σ = (σ + iσ ) =
+ 2
,
2 1 0
1 0 0
σ− = (σ1 − iσ2 ) = .
2i 1 0
In Eq. (1.24), we then have
βrs = βμrs σμ ,
(1.18)
βμrs = Tr(σμ βrs ),
in detail:
β↑rs = β↑↑
rs β↓rs = β↓↓
rs ,
↓↑
(1.19)
rs = βrs
β+ β− ↑↓
rs = βrs .
1.2.1 Spinor
The algebra representing electrons is the spinor. The Dirac relativistic (spe-
cial relativity) function describes this property well. However, the relativity
seems not so important for the present problem. We now concentrate on the
spinor character of the electron. The field operator has two components in
the spin space:
φ↑ (x)
φ(x) = ,
φ+
↓ (x)
1 0
φ̄(x) = φ+ (x)σ3 = (φ+
↑ (x) φ↓ (x))
0 −1
= (φ+
↑ (x) −φ↓ (x)). (1.25)
This is called the Nambu representation.11 The negative sign in front of φ↓
is seen in the Dirac conjugate φ+ → φ̄.
The field operators satisfy, of course, the anticommutators
[φα (x), φ+
β (s )]+ = δα,β δ(x − x ), (1.26)
Theory of Superconductivity 7
where (α, β) = (↑, ↓) and x = (r, t). Then, for the spinors (1.25), the matrix
commutator holds:
φ↑ (x) +
[φ(x), φ̄(x )]+ = , φ↑ (x ) −φ↓ (x )
φ+↓ (x) +
+
φ↑ (x), φ↑ (x ) +
φ↓ (x ), φ↑ (x) +
= +
φ↓ (x), φ+ ↑ (x ) + φ↓ (x ), φ+
↓ (x) +
δ(x − x ) 0
= . (1.27)
0 δ(x − x )
φα (x) → φα (x)eiΛ ,
−iΛ
(1.28)
α (x) → φα (x)e
φ+ .
It is recognized that the Hamiltonian and the equation of motion are invari-
ant under this transformation. We can see that this transformation is a
rotation around the σ3 axis with Λ,
3
φ(x) → eiσ Λ φ(x), (1.29)
since
3
eiσ Λ = cos(σ3Λ) + i sin(σ3Λ)
(σ3Λ)2 (σ3Λ)4 (σ3Λ)3 (σ3Λ)5
=1− + + ··· + i 1 − + + ···
2! 4! 3! 5!
Λ2 Λ4 Λ3 Λ5
=1− + + · · · + iσ 1 −
3
+ + ···
2! 4! 3! 5!
iΛ
e 0
= σ cos Λ + iσ sin Λ =
0 3
.
0 e−iΛ
the density and the current are, in terms of the Nambu spinor,
j 0 = φ+ (x)σ3 φ(x),
(1.30)
j = −i (φ+ (x)∇φ(x) + ∇φ+ (x)φ(x)).
2m
Then the continuity relation is
∂t j0 (x) + ∇ · j = 0. (1.31)
If the system is static, the density must be conserved:
∂t j 0 = 0. (1.32)
Put
G= dxj 0 (x), (1.33)
and if it is found that
[G, Ψ (x)] = Ψ(x), (1.34)
and the expectation value of Ψ(x) over the ground state does not vanish,
0|Ψ|0 = 0, (1.35)
i.e. if the ground state satisfying the relation (1.33) does not vanish, we can
expect the appearance of a boson Ψ(x), whose mass is zero. This boson is
called a Goldstone boson, and the symmetry breaking takes place in the
system. This is what the Goldstone theorem insists on. The details are in
the standard book on the field theory.8
Now we apply this theorem to the superconductivity. The invariant
charge is
Q = d3 xφ+ α (x)φ α (x) = d3 xφ+ (x)σ3 φ(x). (1.36)
Notice that here the same symbol φ is used for the ordinary field with spin
and the Nambu spinor. The above are nothing but the Cooper pairs, and
we now find the Goldstone bosons Ψ∗ and Ψ. In literature, it is noted that
Δ±
0|φ+ (x)σ± φ(x)|0 = ± = 0, (1.39)
g
where Δ and g are the gap and the coupling parameter, respectively. We
expect the estimate
Δ+ ∼ Δ− = Δ
to be reasonable.
The Cooper pairs are now the Goldstone bosons. A comment about
the Goldstone boson or the massless elementary excitation with k = 0 is
required. Using Eq. (1.33), we write the Goldstone commutator (1.34) as
d3 y[0|j0 (y)|nn|φ (x)|0
which is met with the requirement (1.34). The excitation with Mn = 0 needs
no excitation energy, suggesting the Goldstone boson.
We further note that the mass-zero excitation is the imaginary quantity.
This suggests the current to be the phase current, as is seen in the Josephson
effect.
Before closing this preliminary discussion, we want to make a few
remarks. At the beginning, we mentioned London’s postulate that the super-
conducting state is characterized by a statement that the wave function is
rigid, so that the current is entirely the diamagnetic current due to only the
vector potential A. “Rigid” is not really rigid, but it is understood that the
spatial derivative is vanishing, or the current flows along the entirely flat
path, which is described in a textbook as the path going around the top of
a Mexican hat. Boldly speaking, the electron in the superconducting state
is massless. Also, we have pointed out that the vector potential A, which
leads to the Meissner effect, satisfies the covariant relation (1.6),
∂μ Aμ = −k2 , (1.44)
so that a photon is massive in the superconductor.
In the following chapters, we develop a substantial microscopic explana-
tion of the above assertions.
1.3 Propagator
In the previous section, we have found that, as is seen in Eq. (1.41), the
superconducting state strongly concerns the gap function or the anomalous
Green function. We want to deal with the solid-state substances. However,
the infinite crystals described by the single band have already been fully
investigated in literature, and the recent investigations were carried out on
objects with multiband structure.13, 14 The infinite system with many bands
is constructed from the unit cell, which is really a chemical molecule. The
atoms in this molecule give the band index of the real crystal. In this respect,
we first investigate the Green function of a unit cell. The Green function is
now shown in the site representation.
Corresponding to the spinors (1.24), we define the spinor in the site
representation as
ar↑
ar = , ār = (a+r↑ −ar↓ ). (1.45)
a+
r↓
Theory of Superconductivity 11
= −θ(τ1 − τ2 )ar (τ1 )ās (τ2 ) + θ(τ2 − τ1 )ās (τ2 )ar (τ1 ), (1.47)
where τ is the imaginary time, so that the propagator is the temperature
Green function or the Matsubara function. In what follows, we put
τ = τ1 − τ2 , (1.48)
1.3.1 Hamiltonian
Various Hamiltonians can be written by using the charge density matrix,
ρr↑s↑ (x) ρr↑s↓ (x)
ρrs (x) = = σμ ρμrs , (1.50)
ρr↓s↑ (x) ρr↓s↓ (x)
where the basis orbitals are put to be real, so that we do not need the
conjugation procedure for field operators.
The Hückel Hamiltonian or the single-particle Hamiltonian has the struc-
ture
we have17
1
H dir = (ār σa as )vrs;tu
ab
(āt σb au ) (a, b =↑, ↓),
2
1
H ex = − (ār σa as )vrs;ut
aa
(āt σa au ) (a =↑, ↓),
2 (1.54)
1
H super = {(ār σ+ as )vrs;tu
+−
(āt σ− au )
2
−+
+ (ār σ− as )vrs;tu (āt σ+ au )}.
For the direct interaction, the quantum-mechanical Hamiltonian is
H = dx dx χ∗r↑ (r)χs↑ (r)(x)v(x − x )χ∗↓ (x )χ∗u↓ (x ).
dir
1.4 Noninteracting
“Noninteracting” implies that the Hamiltonian is bilinear with respect to
operators so that diagonalization is always possible. It should be instruc-
tive to begin with the single-particle case, since even if we manipulate the
complicated two-particle case, the procedures are almost the same when the
mean-field approximation is employed.
The energy for the Hamiltonian (1.52) is
E 0 = Tr(ρ̂H 0 ) = hrs ār as , (1.55)
where ρ̂ is the statistical operator,
0 −Ω)
ρ̂ = eβ(H , (1.56)
Theory of Superconductivity 13
Then we have
1 1
as ār = ⎛ ⎞
β n (s|ii|r)↑
0
⎝ iωn −i↑ ⎠
(s|ii|r)↓
0 iωn −i↓
⎛ ⎞
(s|ii|r)↑
0
1 ⎜⎜ iωn − i↑
⎟
⎟
=
β n ⎝ (s|ii|r)↓ ⎠
0
iωn − i↓
(s|in(i )i|r)↑ 0
= , (1.66)
0 (s|in(i )i|r)↓
where
1
n(i↑ ) = . (1.67)
1+ ei↑ /kB T
Another sophisticated way starts from the decomposed Hamiltonian
H 0 = ār σμ hμrs as . (1.68)
The commutator is evaluated as
[as , ār σμ hμr s as ] = δsr σ0 σμ hμr s as = σμ hμss as . (1.69)
The equation of motion
(iωn δss − σμ hμss )as , ār ; ωn = σ0 δsr . (1.70)
In the matrix notation,
(iωn − σμ hμ )a, ā; ωn = σ0
or
σ0
a, ā; ωn = . (1.71)
iωn − σμ hμ
In the representation where h is diagonal,
r|h|s = r|ii|h|ii|s = i r|ii|s, (1.72)
the relation (1.71) becomes
(|ii|)(iωn + σμ μi )
a, ā; ωn =
(iωn − σμ μi )(iωn + σμ μi )
(|ii|)(iωn + σμ μi )
= . (1.73)
(iωn )2 − (μi )2
Theory of Superconductivity 15
Now
iωn + σμ μi iωn σμ μi
μ = +
(iωn )2 − (i )2 (iωn )2 − (i )2 (iωn )2 − (μi )2
μ
1 1 1 σμ 1 1
= + + −
2 iωn − μi iωn + μi 2 iωn − μi iωn + μi
1 1 1
= + (summing ωn )
2 iωn + μi iωn + μi
1 σμ
→ (n(μi ) + n(−μi )) + (n(μi ) − n(−μi )). (1.74)
2 2
We evaluate
σ0 μ μ σμ μ μ
G↑rs ↑
= Trσ r|ii|s (n(i ) + n(−i )) + (n(i ) − n(−i ))
2 2
= [r|ii|s]↑ n(↑i ). (1.75)
We then have
1 ↑ − 1 r|ii|s
Grs (0 ) = = r|in(↑i )i|s. (1.76)
β n β n iωn − ↑i
This relation holds for both ↑ and ↑, and we recover the result (1.66).
Along the way, we give the energy expression for Eq. (1.51):
E 0 = Tr(hrs Gsr )
h↑rs 0 s|in(↑i )i|r 0
= Tr ↓ . (1.77)
0 hrs 0 s|in(↓i )i|r
It may be needless to present another illustration:
E 0 = Tr(hrs Gsr ) = Tr(σμ hμrs )(συ Gυsr ). (1.78)
The result is meaningful if σμ συ = σ0 , which leads to, in the present case,
σμ = συ = σ↑ , or σμ = συ = σ↓ .
Such manipulations will be used in the later investigation.
1.5 Interacting
In this chapter, the electron–electron interactions are taken into account,
and we will discuss how they lead to the superconducting state. The
16 Modern Aspects of Superconductivity: Theory of Superconductivity
Taking the (r, s) matrix elements and the mode c which is achieved by the
operation,
Trσc Grs = Gc rs, (1.90)
we select the terms on the right-hand side with the mode c satisfying
Trσc σa = 1. (1.91)
Otherwise, the Tr operation leads to the vanishing result.
Carrying out the summation over ωn , we get
1 c
Gcrs (0− ) = G (ωn )
β n rs
1 c a 1 1 1
= r |i +
β n 2 iωn − (ηnorm
i + ηsup
i ) iωn + (ηnorm
i + ηsup
i )
1 c a Diia
× ia |sc + r |i
β n 2(ηnorm
i + ηsup
i )
1 1
× − ia |sc
iωn − (ηnorm
i + ηsup
i ) iωn + (ηnorm
i + ηsup
i )
1
= r c |ia (n(ηnorm
i + ηsup norm
i ) + n(−ηi − ηsup
i ))i |s
a c
2
Diia
+ r c |ia
2(ηnorm
i + ηsup
i )
× (n(ηnorm
i + ηsup
i ) − n(−ηi
norm
− ηsup
i )) i |s
a c
1 Diia ηi
= δrs − r |i
c a
sup tanh ia |sc , (1.92)
2 2(ηnorm
i + ηi ) 2kB T
where
1
n(ηai ) = a (1.93)
1+ eηi /kBT
and
n(η) + n(−η) = 1,
η
n(η) − n(−η) = −tanh .
2kB T
In the estimation of matrix elements, the chemical potential which has been
disregarded up to now is taken into account. Namely, the Hamiltonian has
an additive term −μār ar , which causes
ηai → ηai − μ < 0, while −ηai → −(ηai − μ) > 0.
Theory of Superconductivity 19
polyacene
t1
n n+1
t2
unit cell
been said that this discontinuity prevents the superconducting phase from
arising. However, this is an old-fashioned assertion and now seems sceptical.
In the Hückel theory, the interaction matrix elements are put as t1 for
the shorter bond and t2 for the longer bond:
H 0 = −t1 (a+ 2n a1n + a1n a2n ) − t2 (a1n+1 a2n + a2n an+1 ),
+ + +
(1.105)
where we consider the neighboring unit cells numbered n and n + 1, and
each cell has two kinds of bonds. The transfer integrals are parametrized as
δ
t1 = t − δ, t2 = t + δ δ > 0, 1 , (1.106)
t
and then the Hamiltonian is easily diagonalized as
k = ±[t2 (1 + cos k) + δ2 (1 − cos k)]1/2 , (1.107)
where + and − correspond to the conduction and valence bands, respec-
tively. We are interested in the features at the zone boundary, k = π:
cπ = 2δ, vπ = −2δ. (1.108)
Here, the superscripts c and v indicate the conduction and valence bands,
respectively. When δ = 0, certainly we have the gap, and if δ = 0, the two
bands continuously join into a single band called the half-filled band.
Next, we turn to the polyacene, whose unit cell is the butadiene molecule.
In this case, we obtain four bands:
1
ck = [t3 + (t23 + 4|t̃k |2 )1/2 ] = −vk ,
2
(1.109)
1
k = [−t3 + (t3 + 4|t̃k | ) ] = −k ,
c 2 2 1/2 v
2
where
t̃k = t1 + t2 eik . (1.110)
The usual pairing property in alternant hydrocarbons is also seen in this
case. Employing the parametrization (1.106) gives, at k = π,
2
1 δ
vπ = [t3 + (t23 + 16δ2 )1/2 ] ≈ 4t3 . (1.111)
2 t1
When δ = 0 (without bond alternation), we have
1
vk = {t3 − [t23 + 8t2 (1 + cos k)]1/2 }. (1.112)
2
Consider single-particle states. If δ = 0 and k = π we have, from Eq. (1.110),
t̃ = 0, so that the amplitudes at sites 2 and 5 vanish. Therefore, when δ = 0,
we get at k = π
1
|ψvk = (a+ − a+ 4k )|0,
2 1k
(1.113)
1 +
|ψk = (a1k + a4k )|0.
c +
2
Theory of Superconductivity 23
It is seen that |ψvk is antisymmetric about the C2v symmetry axis, and |ψvk
is symmetric. Thus, v and c bands are not continuous at k = π.
+ t(a+ + +
n+1,1 an2 + an+1,4 an3 + H.c.), (1.118)
where the second line connects the unit cells n and n + 1.21
The band structure of levels and the linear combination of atomic orbitals
(LCAO) coefficients U are
t
1 (k) = {1 + s(k)},
2
t
2 (k) = − {1 − s(k)},
2
(1.119)
t
3 (k) = {1 − s(k)},
2
t
4 (k) = − {1 + s(k)},
2
with
√
s(k) = 9 + 8 cos k, (1.120)
⎛ ⎞
N4 N3 N2 N1
⎜−N4 4 /t̃k −N3 3 /t̃k −N2 2 /t̃k −N1 1 /t̃k ⎟
U =⎜
⎝ N4 4 /t̃k
⎟, (1.121)
−N3 3 /t̃k N2 2 /t̃k −N1 1 /t̃k ⎠
−N4 N3 −N2 N1
where, for example,
|t̃k |2 1 eik/2
N12 = , = . (1.122)
2(|t̃k |2 + 21 ) t̃k 2t cos(k/2)
At this stage, we have completed, in principle, the usual many-electron
problem. The mean-field approximation makes the interaction problem a
one-particle problem even though the SCF treatment is required at each
step. In other words, from the viewpoint of the Hückel theory, the spin-
diagonal parts provide the answer. On the other hand, the spin-off diagonal
part, which means less in the case without electron–electron interactions, is
responsible for the superconductivity.
Up to the previous chapter, the problem had been investigated in the
site representation. That is to say, the system is considered to be composed
Theory of Superconductivity 25
of N sites. However, the real substance is formed from unit cells, so that the
system is a repetition of the unit cell. The usual band theory of polyacene
has thus been completed at this stage.
We turn to the onset of superconductivity. In this case, the single-particle
approach is almost meaningless, but the pair state — say, the wave function
of a Cooper pair — should be investigated. For this purpose, the Green
function of a Cooper pair is most preferable. The gap equation (1.98) is
nothing but the SCF equation for a Cooper pair.
The electronic structure of the single butadiene molecule referring to
k = 0 is suggestive. Let the total number of sites of high polymer polyacene
be N . The number of sites in the unit cell is four. Then N = 4n, with
the number of unit cells n. The numbers 1–4 are the band indices; then
we have the chemical potential between the 2n level and 3n level. The
Cooper pair should be the hole pair of the 2n level indicated by mode (−)
or the particle pair of the 3n level. The discussion is confined to these levels
in solid-state physics, even if the interaction with other bands is taken into
account.
The 2n and 2n + 1 levels are called highest occupied molecular
orbital (HOMO) and lowest unoccupied molecular orbital (LUMO). These
features and behaviors are not so far or qualitatively the same as those
for k = 0.
The electronic structure of a single butadiene molecule is suggestive. The
levels and the bond orders are (the unit = t)
1 −1.618
2 −0.618
(1.123)
3 0.618
4 1.618.
Here, 1 and 2 are occupied (valence) levels, while 3 and 4 are unoccupied
(conduction) ones. Let the probability amplitude, with which the electron
on level i is found at the site r, be r|i, then the bond order qrs is defined
as (at the zero temperature)
occ
qrs = r|ii|s, (1.124)
i
where the summation includes the spin state. In the determination of
the attractive electron–electron interaction, the bond order is of crucial
importance.
26 Modern Aspects of Superconductivity: Theory of Superconductivity
(1) The denominator of Eq. (1.98) is the sum of the ordinary Hartree–Fock
energy and that of the superconducting state, and the latter is
+− −
ηsup
i = −vii;tu Gtu . (1.127)
ηi = ηnorm
i + ηsup
i (1.128)
(4) The energy interval establishing the superconductivity in the BCS the-
ory is related to the electron–phonon interaction, ωD (Debye fre-
quency). In the present theory, it is replaced by the band width nearly
equal to |g|.
(5) Assuming that
− −
G+ +
s↑,r↓ (0 ) = Gt↑,u↓ (0 ) (1.130)
in Eq. (1.98), we have
|g|
dξ ξ + ηsup
1 = gN (0) tanh , (1.131)
0 ξ + ηsup 2kB T
where ξ = ηnorm . The critical temperature TC is determined by the
condition that ηsup vanishes at this temperature. The integration in
Eq. (1.131) is carried out as usual. Approximating
d3 k
= N (0) dξ
(2π)3
gives
|g| Z
dξ ξ dz |g|
tanh = tanh z, z=
0 ξ 2kB T 0 z 2kB T
∞
= [ln z tanh z]Z
0 − dz ln z sech2 z, (1.132)
0
where the upper limit in the second integration is replaced by ∞, which
makes it integrable15 :
∞
4eγ
dz ln z sec h2 z = − ln ; γ is the Euler constant.
0 π
A simple rearrangement of the result gives
2eγ
kB TC = |g|e−1/N (0)g ∼ 1.13|g|e−/N (0)g. (1.133)
π
The result is entirely the same as the current one. However, since it is
probable that
|g|
∼ 100, (1.134)
kB T
the critical temperature is, at most, enhanced by this value, even though it
is considerably reduced by the factor e−1/N (0)g .
1.6.4 Conclusion
As has been presented, superconductivity is not a too-complicated phe-
nomenon. If we employ the spinor representation, superconductivity is
28 Modern Aspects of Superconductivity: Theory of Superconductivity
kσ ck σ ± c−k −σ c−k−σ ).
Bkσ,k σ (c+ +
As seen above, the theory implicitly assumes that the system is a perfect
crystal and is described by the single wave vector k. The positive and nega-
tive signs of k indicate waves propagating from the vertex or off the vertex.
Before entering into the discussion about the relaxation time of the magnetic
resonance, we briefly review the Bogoliubov theory.3, 10, 15 The Bogoliubov
transformation defines a quasiparticle responsible for the superconductiv-
ity as
γk↑ uk −vk ck↑
= (1.136)
γ+
−k↓
vk uk c+−k↓
with
u2k − vk2 = 1.
The spirit of the Bogoliubov transformation is to mix the operators ck↑ and
c+
−k↓ , which are different in spin (as to the wave number, this mixing is not so
serious) and not mixed in the normal situation. The quasiparticle yields the
new ground state near the chemical potential. The stabilization energy thus
obtained is called the gap energy, Δk . What we have done in Sec. 1.3 is a sub-
stantial understanding of this reason. However, if we want to make the ±k
distinction meaningful, it is natural to adopt the four-component spinor —
say, the extended Nambu spinor (perhaps spurious).26 For ordinary states,
⎛ ⎞
ck↑
⎜ + ⎟
⎜c−k↓ ⎟
ck = ck↑ c−k↓ c−k↑ ck↓ , ck = ⎜
+ + +
⎜c+ ⎟,
⎟ (1.137)
⎝ −k↑ ⎠
ck↓
and, for the superconducting state,
⎛ ⎞
γk↑
⎜ + ⎟
⎜γ−k↓ ⎟
γ+
k = γ+
k↑ γ−k↓ γ−k↑ γ+
k↓ , γk = ⎜
⎜γ+ ⎟.
⎟ (1.138)
⎝ −k↑ ⎠
γk↓
These are connected with each other by the Bogoliubov transformation as
γk = Uk ck , γ+ +
k = ck Uk , (1.139)
where
uk 0 uk vk
Uk = , with uk = . (1.140)
0 u+
k −vk u+
k
Careful manipulation is instructive.
30 Modern Aspects of Superconductivity: Theory of Superconductivity
Case 1+ :
kσ ck σ + c−k −σ c−k−σ ) = ck Σ ck = γk Uk Σ Uk γk ,
(c+ + + 3 + 3 +
(1.141)
where
σ3
Σ =3
. (1.142)
−σ3
This is the 4 × 4 matrix manipulation; however, it is enough to note the
upper half of the result.
The upper half of Eq. (1.139) is
+ 1 0 ck ↑
(ck↑ c−k↓ )
0 −1 c+
−k ↓
+ u v 1 0 u −v γk ↑
= (γk↑ γ−k↓ )
−v u 0 −1 v u γ+−k ↓
uu − vv −uv − vu γk ↑
= (γ+ γ−k↓ ) , (1.143)
k↑ −vu − uv
vv − uu
γ+−k ↓
where u and v are the abbreviations of uk and vk , respectively, while u and
v are those of uk and vk .
Case 1− :
(c+ + + + +
kσ ck σ − c−k −σ c−k−σ ) = ck 1ck = γk Uk 1Uk γk . (1.144)
σ
The upper half of the above is
uu + vv
+ −uv + vu γk ↑
γk↑ γ−k↓ . (1.145)
−vu + uv vv + uu γ+
−k ↓
As is shown clearly, the original term is transformed in the quasiparticle
representation into a combination of the scattering term with the diagonal
element in the uk matrix and the creation or annihilation of a pair with the
off-diagonal element of u. These matrix elements are called the coherent
factors.
Let us turn to the case where the spin flip-flop is allowed.
Case 2+ :
kσ ck −σ + c−k σ c−k−σ = ck Σ ck = γk Uk Σ Uk γk
c+ + + J + J +
(1.146)
with
⎛ ⎞
1
⎜ −1 ⎟
ΣJ = ⎜
⎝ −1
⎟.
⎠ (1.147)
1
Theory of Superconductivity 31
uu + vv
+ uv − vu γ−k ↑
Eq. (1.147) = γk↑ γ−k↓ . (1.148)
−vu + uv −vv − uu γ+
k ↓
Case 2− :
with
⎛ ⎞
1
⎜ 1 ⎟
J=⎜
⎝ 1
⎟.
⎠ (1.150)
1
Note that, at present, the scattering terms are off-diagonal, and the cre-
ation and annihilation terms are diagonal.
By the use of the relations3, 10, 15
1 k 1 k
2
uk = 1+ 2
vk = 1− ,
2 Ek 2 Ek (1.152)
Ek = k + k
2 2 2
(k is the gap energy), the coherent factors are expressed substantially as
1 k k k k
1± : (uu ∓ vv )2 = 1+ ∓ ;
2 Ek Ek Ek Ek
(1.153)
2 1 k k
2± : (uv ∓ vu ) = 1∓ .
2 Ek Ek
In case 1+ , we have the ultrasonic attenuation, while the electromagnetic
interaction is in case 2+ and the magnetic resonances are in case 2− .
1.7.2 T1 in NMR
A detailed analysis of the spin–lattice relaxation time T1 in the nuclear
magnetic resonance will be presented in the next section. Here, the results
32 Modern Aspects of Superconductivity: Theory of Superconductivity
spectrum of the propagator GKM (τ) = aKM (τ)a+ KM gives the line shape
of the magnetic resonance.
The nuclear propagator GKM (τ) sees the electron sea, followed by the
electron excitation in the spin space. This gives the additional line width
of the nuclear magnetic resonance. The phenomenon looks like the vacuum
polarization in quantum electrodynamics. The self-energy part that has thus
arisen in the nuclear energy is the source of the line shape of the nuclear
magnetic resonance.23
The spin–lattice relaxation time of the nuclear spin I z is given by the
imaginary part of the magnetic susceptibility χzz , which is equal to (χ+− +
χ−+ )/2 in the spatially homogeneous system. Here, ± correspond to (I x ±
iI y )/2, respectively. The ensemble average of a change, δI + (t), is given by
the linear response theory as
t
δI + (t) = i dt Tr{ρG [H ex (t ), I + (t)]− }, (1.157)
−∞
where ρG is the grand canonical statistical operator. However, the chemical
potential is not given explicitly, unless otherwise stated. The rotating mag-
netic field causing the magnetic transition is, assuming a single mode for
simplicity,
H ex (t) = HR (I + (t)eiωt + I − (t)e−iωt ). (1.158)
As has been said, the spin–lattice relaxation arises from the interaction
between the nuclear spin and the electron spin. In other words, the electron
spins play the role of a lattice system.
H = γgβB F (R, r)I · S, (1.159)
where F (R, r) is a function of spatial coordinates of the nucleus, R, and
that of the electron, r. The term γ is the gyromagnetic ratio of the nucleus,
and g and βB are the g factor and the Bohr magneton of the electron,
respectively.
Now the second quantization of the above is carried out. First of all, the
orthonormalized wave function describing the nuclear behavior, |ξK (R)M ,
is introduced as
(HN + HM )|ξK (R)M = (K + M )|ξK (R), M
= KM |ξK (R, )M , (1.160)
with
KM = K + M,
34 Modern Aspects of Superconductivity: Theory of Superconductivity
where HN is the spatial part and HM is the Zeeman part. Then we have
I α → ξK (R)M |I α |ξK (R)M a+
KM aK M
= M |I α |M a+
KM aK M δKK . (1.161)
A similar equation is given for electrons,
(HS + Hm )φk (r)|m = (k + m)|φk (r), m = km |φk (r), m, (1.162)
where
km = k + m,
so that
H → hγgβB M |I β |M m|S α |m ξK (R)φk (r)|F (R, r)|ξK (R)φk (r)
× (a+ +
KM aK M )(ckm ck m ). (1.163)
If the nuclear motion is assumed to be that of a harmonic oscillator, a+KM
and aKM are the creation and annihilation operators of vibrational excita-
tions. When the nuclei carry noninteger spins, these are considered to obey
the Fermi statistics or satisfy the anticommutation relation
[a+
KM , aK M ]+ = δKK δM M . (1.164)
However, as will be seen in the following, this selection of the statistics is
not fatal for the theory. Needless to say, the operators for electrons satisfy
the anticommutation relations.
The change of I + in Eq. (1.157) can be written, in the interaction rep-
resentation, as (we retain the I − term in Eq. (1.158))
t
δI (t) = iγHR
+
dt e−iωt
−∞
t
= −γHR dt e−iωt DK,M −1,M (t − t)
−∞
∞
1
= − γHR e−iωt dse−iωs DK,M −1,M (s)
2 −∞ KM
1
= − γHR e−iωt DK,M −1,M (ω). (1.165)
2
From this result, the magnetic susceptibility of the present system is
γ
χ+− (ω) = − DK,M −1,M (ω). (1.166)
2
Theory of Superconductivity 35
In the course of derivation, the matrix elements of spin operators are put to
be equal to 1, and then the Tr operation is carried out. Here, DK,M −1,M (s)
is a retarded Green function,
DK,M −1,M (s) = −iθ(s)Tr{ρG [a+ +
KM −1 (s)aKM (s), aKM aKM −1 ]− }, (1.167)
and DK,M −1,M (ω) is its Fourier transform. Now our problem is to estimate
this retarded function.
The retarded Green function is easily obtained by analytical continua-
tion from the Matsubara function (or the temperature Green function with
imaginary time τ) which is causal in τ,
' (!
DK,M −1,M (τ) = −Tr ρG T τ a+ K,M −1 (τ)aK,M (τ)a+
a
K,M K,M −1 , (1.168)
for which the Feynman diagram analysis is available.15
1.7.4 Noninteracting
Here, we deal with the case without the spin–lattice interaction. It might
be trivial; however, it seems instructive for the later investigation. By the
use of the simplified notation, · · · = Tr(ρG · · · ), the Green function in this
case is written as
+ +
DK,M
0
−1,M (τ) = T τ[aKM −1 (τ)aK,M (τ)aKM aKM −1 ]
= GKM
0 0
(τ)GKM −1 (−τ), (1.169)
where
GKM
0
(τ) = −T τ[aKM (τ)a+
KM ], (1.170)
and the corresponding Fourier transform of Eq. (1.169) is
1 0
DK,M
0
−1,M (ωn ) = GKM (υn )GKM
0
−1 (υn − ωn ), (1.171)
β υ
n
where
1
GKM
0
(υn ) = . (1.172)
iυn − KM
Therefore,
1 1 1
DK,M
0
−1,M (ωn ) = ·
β υ iυn − iωn − KM −1 iυn − KM
n
1 1 1 1
= −
β υ iυn − iωn − KM −1 iυn − KM iωn + KM −1 − KM
n
1
= [n(KM −1 ) − n(KM )] , (1.173)
iωn − KI
36 Modern Aspects of Superconductivity: Theory of Superconductivity
where
KI = KM − KM −1 ,
and it is noted that ωn is even, and that it does not matter in
obtaining the particle number. The retarded Green function is obtained
simply by replacing iωn with ω + iη (η is a positive infinitesimal). Thus, we
obtain
1
δI + (t) = e−iωt HR (n(K− ) − n(K+ )) , (1.174)
ω − KI + iη
where the radio frequency stimulating the resonance is rewritten by ω0 . The
magnetic susceptibility χ+− thus becomes
1
χ+− (ω) = γ(n(K− ) − n(K+ )) ,
ω − KI + iη
γ = e−iωt HR , (1.175)
whose imaginary part is
χ+− (ω) = −πγ(n(K− ) − n(K+ ))δ(ω − KI ). (1.176)
This gives the sharp δ function-type energy spectrum. We have no line width
or the relaxation time, and states are stationary.
ω
m+1
λ ρ+ω
M-1 ρ
m
M ω
1
= (n(km ) − n(km+1 )) . (1.181)
ω − km+1 + km
The lattice (electron) gets ω from a nucleus to lift the electron spin from m
to m + 1, and at the other vertex, the inverse process occurs. This looks like
the radiation process in the photochemistry. We thus have the self-energy
part of GK;M −1,M ,
1
ΣK;M −1,M (υn ) = |K|2 (n(km ) − n(km+1 )) , (1.182)
ω − ks
where
ks = km+1 − km .
38 Modern Aspects of Superconductivity: Theory of Superconductivity
Another propagator GK;M , has the self-energy part ΣK;M,M −1 (υn ). This
is built by replacing ρn + ω with ρn − ω in Eq. (1.181).
1
ΣK;M,M −1 (υn ) = |K|2 (n(km ) − n(km−1 ))
ω − km−1 + km
1
= |K|2 (n(km ) − n(km−1 ))
,
ω − ks
1
≈ |K|2 (n(km+1 ) − n(km )) . (1.183)
ω − ks
We now turn to the evaluation of the propagator DK;M −1,M (τ):
DK;M −1,M (ωn )
1
= GKM (υn )GKM −1 (υn − ωn )
β υ
n
1
= (iυn − KM − ΣKM − )−1 (iυn − iωn − KM −1 − ΣK;M,M −1 )−1
β υn
1 1 1
= −
β υn
iυn − iωn − KM −1 − ΣKM −1 iυn − KM − ΣKM
1
×
iωn − KM −1 − ΣKM −1 + KM + ΣKM
1
= (N (KM ) − N (KM −1 )) . (1.184)
iωn − KI + ΣKM − ΣKM −1
Here,
KI = KM − KM −1 ,
and the self-energy parts are disregarded in obtaining the particle density.
We can see that the additional terms in the denominator modify the line
shape.
If we put iωn → ω + iη, we can obtain the retarded Green function,
whose imaginary part gives the line shape:
1
D(ω) ∼
ω + iη − KI + ΣKM − ΣKM −1
1
= P − iπδ(ω − KI + ΣKM − ΣKM −1 )
ω − −ΣKM −1 KI
i(Im)
= (at the resonance point). (1.185)
(Re)2 + (Im)2
The line shape is now changed from the δ function type to the Lorentz type,
as expected.
Theory of Superconductivity 39
with
1 0 0 1 0 −i
σ =
3
, σ =
1
, σ =
2
,
0 −1 1 0 i 0
1 1
σ+ = (σ1 + iσ2 ), σ− = (σ1 − iσ2 ). (1.192)
2 2
As for the electron–electron interaction, only those for the Cooper pairs are
selected; namely, Σ+ selects a Cooper pair in the particle state, and Σ− in
the hole state.
is invariant under the rotation about the Σ3 axis in the space spanned by
Σ3 , Σ+ and Σ− . Observing that
[Σ3 , Σ± ] = ±2Σ± (1.194)
suggests the phase transitions along the Σ± directions.
If we define
c̄ = c+ Σ3 ,
as has been done in Sec. 1.43, the discussions parallel to those there will be
possible in the following. However, this is not employed in this case.
The phase transition cannot be achieved by the perturbational approach,
but the effective Hamiltonian giving the phase transition should be included
at the beginning. For example, the modified Hamiltonian
−
k (k Σ + ρΣ + ηΣ )ck ,
H 0 = c+ 3 +
1 + + − −
(1.195)
H int = g {(ck Σ ck ) · (c+
k Σ ck ) − ck (ρΣ + ηΣ )ck },
+ +
2 kk
Theory of Superconductivity 41
where ρ and η are the so-called gap energies which are assumed to be
independent of k for simplicity. Note that the Hamiltonian H 0 is already
symmetry-broken. Here, we adopt a conventional method. In what follows,
H int is neglected, i.e. we start from the Green function due to the effec-
tive Hamiltonian and then take the interaction (1.158) into account. At this
stage, we may expect the result that will be obtained in such a manner that
the normal and superconducting states contribute additively.
The temperature Green function with the imaginary time τ is defined as
G0kk (τ) = −Tτ [ck (τ)c+
k ], (1.196)
where · · · = Tr{ρG · · · }. The equation of motion of G0kk is
∂τ G0kk (τ) = −∂τ [θ(τ)ck (τ)c+ +
k (0) − θ(−τ)ck (0)ck (τ)]
= δ(τ)[ck (τ), c+ +
k (0)]+ + T τ[ck (τ), H ]− , ck (0)
0
The first-order term has already been evaluated in Eq. (1.181). Then we
have to carry out the complicated manipulation due to the second term of
the last line in Eq. (1.203). However, the combinations other than those
satisfying the resonance condition ω ∼ ks [see Eq. (1.182)] should give a
small effect which is to be neglected.
We then consider that the procedures of the previous section need not
be repeated, and we are allowed to multiply the results there by the factor
(ρ + η)/2k↑ as the perturbing correction or the superconducting effect.
1.8.4 Conclusion
If we review the present investigations from the viewpoints of the line shape
problem in the magnetic resonance, three cases are clearly distinguished.
In the noninteracting case, the line shape is written in terms of the delta
function. In the interacting case of the normal phase, it is presented by
the Lorentz-like function, and in the superconducting phase it is further
multiplied by a coherent factor; whereas the statistical factors referring to
the nuclear states never change throughout.
Theory of Superconductivity 43
Fig. 1.3 In the BCS model, superconductivity arises due to the attractive electron–
electron interaction which appears owing to the scattering by phonons.
Fig. 1.4 Excitation of the ground state of a Cooper pair demonstrates the presence of
the energy gap Δ, which is (very approximately) isotropic in the momentum space.
46 Modern Aspects of Superconductivity: Theory of Superconductivity
a
Actually, dφR dφI = 1
2i
dφ∗ dφ. The constants are adsorbed in the normalization factor.
48 Modern Aspects of Superconductivity: Theory of Superconductivity
√
+ iβ πgα,k −k φα (υn )a∗k (ωn + υn )σ+ ak (ωn )
√
+ iβ πgα ,k−k φ∗α (υn )a∗−k (ωn )σ− a−k (ωn + υn )}. (1.222)
We can integrate with respect to fermion variables: note that31
∗
dz ∗ dze(z Az) = det A for a fermion,
∗
dz ∗ dze(z Az) = (det A)−1 for a boson. (1.223)
where we have made use of a symmetry property, gα,k −k = gα ,k−k , and
· · · stands for a matrix. This can be also rewritten as
Z = Dφ∗ Dφ exp−βW (φ), (1.225)
= 0. (1.228)
We thus obtain
1 (uδ,k uα,k ) · φ̄α
πφ̄δ (υn ) = , (1.229)
β (iωn )2 − Ek2
iωn
where
+
Ek2 (υ) = 2k + (uα,k uα ,k )φ̄∗α (υn )φ̄α (υn ). (1.230)
Note that the field strength φ∗α (υn )φα (υn ) has the dimension of energy.
In obtaining Eq. (1.229), we first rationalize the denominator, and then
the Tr operation on σ’s is carried out. The numerator that remains is the
coefficient of Tr σ+ σ− = 1. A similar equation is obtained for φ∗ . Equation
(1.229) corresponds to the gap equation.
Now we employ an approximation in the rest of this study where, as was
mentioned at the beginning, uα,k = u (a positive constant). Thus, to the
first approximation, φ is nearly constant:
π 1 1
= . (1.231)
u β ω (iωn )2 − Ek2
n
Theory of Superconductivity 51
For the convergent integrals which are the second and third terms, the
integration limit is extended to infinity, and then we obtain (putting =
ω tan θ)
ωD 1
2π d β 1 π
= tanh + u|φα |2
uN (0) 0 2 β4 ω3n
1 3π 1
+ (u|φα |2 )2 + ··· . (1.240)
β 16 ω5n
The integration of the first term on the right-hand side, which is similar
to (1.235), has β instead of βc . Then, using (1.236), we have, in the lowest
approximation,
βc 1π 1
log = u|φα |2 + ··· . (1.241)
β β4 ω3n
Using
∞
1 2p − 1
= ζ(p) (1.242)
(2n + 1)p 2p
0
1 1
even
√
− Tr [(ωn )i u(Φα (υn ) · τ)]n , (1.246)
βω υ n n
n n
where
√
iωn + k τ3 + i u(φ̄α · τ)
k (ωn ) = . (1.247)
(iωn )2 − Ek2
Note that now the propagator includes the mean-field effect. What we are
interested in is the second line of Eq. (1.246).
The term with n = 2 is precisely written as
−1
(a) ∗(a )
Tr uk (ωn ) Φα (υn )τ(a) k+α (ωn + υn ) Φα (υn )τ(a )
2β ω υ
n n
= u|Φα (υn )|2 Pα (υn ), (1.248)
υn
where careful manipulations about τ’s are required. Except for the explicit
ones combined with Φ’s, we have τ’s inside ’s. Let us call the terms with
a = a = 1 and a = a = 2 the direct interactions, and the terms with
a = 1, a = 2 and a = 2, a = 1 the cross interactions. We then have
−2 iωn (iωn + iυn ) − k k+α + uφ̄∗ · φ̄
Pα (υn ) = 2 ], (1.249)
β ω [(iωn )2 − Ek2 ][(iωn + iυn )2 − Ek+α
n
where the first two terms in the numerator arise from the direct interaction
and the third from the cross interaction. Here, use has been made of the
54 Modern Aspects of Superconductivity: Theory of Superconductivity
11 1 1
=β −
β υ 2 iυn − (Ek+α − Ek ) iυn + (Ek+α − Ek )
n
β
= [−nB (Ek+α − Ek ) + nB (−Ek+α + Ek )],
2
−β β(Ek+α − Ek )
= cot h ,
2 2
where nB is the Bose function, and use of Eq. (1.232) has been made for
the even frequency υn . Then Eq. (1.250) becomes
β β(Eα+γ − Eα )
Pγ (υn ) = n(Eα )(1 − n(Eα+γ )) cot h
υn
2 2
Eα+γ Eα − α+γ α + u(φ̄∗ · φ)
× . (1.251)
Eα+γ Eα
We then obtain the effective action for φ up to the second order in Φ, as
follows:
1
Wα (φ∗ , φ) = π φ∗α (υn )φα (υn ) − Tr
υn
βω υ
n n
√
× log{−β(iωn − k τ − uφ̄α (υn ) · τ)}
3
Ek+α − Ek
− n(Eα )(1 − n(Eα+γ ))
υn
(iυ) − (Ek+α − Ek )2
2
Ek+α Ek − k+α k + u(φ̄∗ · φ̄)
× . (1.252)
Ek+α Ek
Theory of Superconductivity 55
1.9.8 Instability
In order to discuss the dynamics of Φα (υn ), we must obtain an expression
for the action of Φα (υn ) similar to that for aα (ωn ) given in the exponent of
Eq. (1.208). Since we are interested in the energy region much lower than
that for the electronic excitation (Eα+γ − Eα ), the denominator of Pα (υn )
in Eq. (1.250) is expanded as
1 −1 (iυn )2
= 1+ + ··· .
(iυ)2 − (Ek+α − Ek )2 (Ek+α − Ek )2 (Ek+α − Ek )2
(1.253)
Substituting this into Pα (υn ) in Eq. (1.250) and taking up the terms of the
order of (iυ)2 in the above, we have
(Ek+α − Eγ )2
(iυn )2 + u(φ̄∗ φ̄) + ··· . (1.254)
Ek+α Ek
The second term is a small and complicated, but positive quantity. We thus
obtain the effective Lagrangian for Φ up to the second order of Φ∗ Φ:
LB ∼ Φ∗ (iυn )2 Φ + ηΦ∗ Φ + · · · , η > 0. (1.255)
In the above, iυn is replaced by ∂/∂τ in the future. At this stage, the details
of η are immaterial except that this is positive.
To put forward the problem, we need the term proportional to Φ∗ Φ∗ ΦΦ.
This is obtainable from the term with n = 4 in Eq. (1.246). Without any
detailed calculation, we may presume a positive quantity for this, noted
by Λ2 . Then the total Lagrangian of the boson field becomes
LB ∼ Φ∗ (iυn )2 Φ + ηΦ∗ Φ − Λ2 Φ∗ Φ∗ ΦΦ. (1.256)
This looks like the Lagrangian that Ginzburg and Landau have used for
investigating the superconductivity. We thus expect a kind of phase transi-
tion. Let us optimize, with respect to Φ∗ , the potential part of Eq. (1.256):
V = −ηΦ∗ Φ + Λ2 Φ∗ Φ∗ ΦΦ, (1.257)
∂V
= −ηΦ + 2Λ2 Φ∗ ΦΦ = (−η + 2Λ2 Φ∗ Φ)Φ = 0. (1.258)
∂Φ∗
If Φ = 0, nothing happens; while, in the case of Φ = 0, it is possible that
η
Φ∗ Φ = ≡ ā2 . (1.259)
2Λ2
We now put
Φ = (ρ + ā)eiθ . (1.260)
56 Modern Aspects of Superconductivity: Theory of Superconductivity
Namely, only the radial part is affected. This choice is allowed in the spirit
of the unitary gauge. If we insert this into (1.257), it follows that
Observing the final result, we realize that the potential becomes mini-
mum at |Φ| = ā; namely, at this position apart from the origin, |Φ| = 0,
which corresponds to the HF state. The mass of the field Φ is given by
the coefficient of Φ∗ Φ, which is 4Λ2 ā2 for a particle of the ρ field, while
the term involving θ is completely lost. That is to say, the particle in the
θ field is massless, and it is called the Goldstone boson. This is due to the
occurrence of the infinite degeneracies along the direction perpendicular to
the ρ coordinate.
1.9.9 Supersymmetry
The total Lagrangian of the system is
where
In the above, the second line is obtained from Eqs. (1.256) and (1.261),
and then the third line from Eq. (1.245) supplied by necessary terms, and a
compact notation is used, and the immaterial constant terms are omitted.
The L in the above is, strictly speaking, the Lagrangian density, and the
action which we are interested in is the space–time integral of the Lagrangian
density, so that the partial integration with respect to τ (after the replace-
ment iυn → ∂τ ) gives the first term with the minus sign. If we define the
momentum conjugate to Φ,
∂Ltot ∂Ltot
Π= = −Φ̇∗ , Π∗ = = −Φ̇, (1.264)
∂ Φ̇ ∂ Φ̇∗
the total Hamiltonian for the boson field is obtained as
Htot = Π∗ Φ̇ + ΠΦ̇∗ − L
√
= Π∗ Π + i u(a∗ Φτ+ a + a∗ Φ∗ τ− a) + Λ2 [Φ∗ Φ − ā2 ]2 . (1.265)
Theory of Superconductivity 57
Λ2 = 2u, (1.273)
37
which is reasonable. However, as mentioned before, this had to be obtained
analytically from the term with n = 4 in Eq. (1.252). If it is assumed that
the present system is the case of supersymmetry, we can avoid this tedious
calculation.
The nilpotency of the charge operators,
Q2 = 0, (Q∗ )2 = 0, (1.274)
is clearly preserved from a2 = 0 and (a∗ )2 = 0. We have thus completed a
supersymmetric analysis.
1.9.11 Discussion
Some people seemed to be puzzled by the curious look of the GL equa-
tion when it was published: the superconductivity being a macroscopic,
thermodynamically stable phase — why is it predicted by a Schrödinger-
like wave function, to which the microscopic phenomena are subjected?
Several years later, the microscopic theory was established by BCS, and
soon Abrikosov had successfully correlated these two treatments. He made
clear that the GL function is deeply concerned with the gap function char-
acteristic of the superconductivity and is a wave function describing the
condensed Cooper pairs.
The superconducting state is certainly a stable thermodynamical state,
and a phase transition from the normal state to the superconducting state is
interpreted as a long-distance correlation between the Cooper pairs. Yang38
developed a unified treatment of the phase transition in terms of the density
matrix. According to his treatment, the onset of the superconductivity is
understood in such a way that the off-diagonal long-range order (ODLRO)
of the second-order density matrix has a nonvanishing value. This concept
is clearly related to London’s rigid function, the quasiboson condensation
being widely seen as a powerful model of superconductivity and the vari-
ational wave function tried by BCS.3 Recently, Dunne et al.39 applied the
concept of ODLRO to argue the high-temperature superconductivity in cop-
per oxide, where the attractive interaction between electrons is assumed to
have originated from Friedel oscillations in the screened potential.
However, the previous presentation of Abrikosov looks to be a detour;
a complicated and tedious procedure in which the condensed pair is char-
acterized by an anomalous temperature Green function which is difficult to
manipulate for beginners. Therefore, the direct way to reach the GL theory
from the BCS Hamiltonian should be preferable. What we have done in the
present investigation is the following: the auxiliary boson field driven by
the Hubbard–Stratonovitch transformation to eliminate the quartic term of
electron operators is just the GL function.
In conclusion, from the treatments used so far, the conditions which are
necessary for the occurrence of superconductivity are:
(1) The wave function which means the ground-state average of operators
of a Cooper pair must be complex or two-dimensional. If one of these
two degrees of freedom gets a new stable structure, the other degree of
freedom, whose direction is perpendicular to the former, offers the infi-
nite degeneracy. In the Nambu theory, the first refers to the σ1 direction,
and the second to σ2 in the fictitious spin space.
60 Modern Aspects of Superconductivity: Theory of Superconductivity
2.1 Introduction
In 1986, K. A. Müller and J. G. Bednorz discovered the phenomenon of
superconductivity in Cu-oxide compounds of lanthanum and barium at
the temperature Tc = 35 K,2 and they were awarded the Nobel Prize
in 1987. This discovery gave an additional push to the intensification of
the scientific activity in the field of superconductivity. Over the last 23
years, the superconducting transition temperature has been increased to
140 K. Moreover, one may expect the discovery of new superconductors with
higher critical temperatures. Recently, the high critical temperatures of new,
Mg-containing superconductors (MgBr2 ) (Tc = 39 K) and superconducting
oxypnictides LnO(1−x) Fx FeAs (Tc = 50 K) were registered.
The main purpose of solid-state physics is the development and fabri-
cation of substances which possess superconductivity at room temperature.
At present, there is an intensive search for such high-temperature supercon-
ductors (HTSCs). Especially promising is the method of producing super-
conducting materials with the help of laser spraying of layers.
After the discovery of HTSCs, the following problems became urgent:
(1) Conceptual–theoretic: the problem of clarifying mechanisms of high-
temperature superconductivity.
(2) Engineering and technical: the problem of practical applications of
HTSCs.
(3) Research: the problem of the search for materials with higher Tc .
The great interest in the phenomenon of superconductivity is caused by the
basic possibility of using it in the future for electric power transfer with-
out losses and for the construction of quantum high-power generators. This
61
62 Modern Aspects of Superconductivity: Theory of Superconductivity
superconductors with pairing of the BCS type, the energy gap is isotropic
(s pairing), and the temperature dependence of the heat capacity has an
exponential form ∼ exp−Δ/kB T , where Δ is the superconductor’s gap. In
superconductors with anisotropic pairing, the temperature dependence of
the heat capacity has a power character, namely T n . The appearance of
such temperature dependences is related to the fact that the superconduc-
tor’s gap is zero on the Fermi surface. The development of the thermody-
namics of this model is the actual problem of HTSCs.
In the following section, we consider the history of the development of
studies on the phenomenon of superconductivity. The structure of HTSCs
and their physical properties are analyzed. We make a survey of the mech-
anisms of superconductivity and discuss the problem on the symmetry of
the order parameter.
We draw the conclusion that, on the whole, the most probable is the “syn-
ergetic” mechanism of pairing which includes, as components, the electron–
phonon and spin-fluctuation interactions in cuprate planes.
Only if the interaction of all the degrees of freedom (lattice, electronic
and spin ones) is taken into account can a number of contradictory prop-
erties observed in the superconducting and normal phases be explained.
Moreover, one needs to consider the complicated structure of HTSCs. It
will be emphasized that new experiments should be executed in order to
explain the available experimental data.
In this section, we also consider the thermodynamics of the spin-
fluctuation mechanism of pairing, present the method of functional inte-
gration for the calculation of thermodynamical properties on the basis of
the Pines spin-fluctuation Hamiltonian, and deduce the Schwinger–Dyson
equations for Green functions and equations for the thermodynamic poten-
tial. Based on the Schwinger–Dyson equation, we obtain equations for the
superconductor gap, which are used in numerical calculations of the thermo-
dynamics of HTSCs. We present analytic formulas for the thermodynamic
potential and its jump. The numerical calculations of the temperature
dependence of the electron heat capacity indicate that the temperature
dependence of the heat capacity is proportional to the square of the tem-
perature. We emphasize that such temperature dependence is related to the
d pairing. It is shown that the measurement of the temperature dependence
of the heat capacity can be a supplementary test for the establishment of
a type of the symmetry of pairing in HTSCs. The jump of the heat capac-
ity of cuprate superconductors near the critical temperature is evaluated
as well.
64 Modern Aspects of Superconductivity: Theory of Superconductivity
accepted. Therefore, the BCS theory was needed in substantiation with the
help of stricter arguments. This was made by academician N. N. Bogoliubov
in Ref. 42. He developed a microscopic theory of the phenomena of super-
conductivity and superfluidity. By using the Fröhlich effective Hamiltonian,
Bogoliubov calculated the spectrum of excitations of a superconductor
within the method of canonical transformations proposed by him in 1947.
In 1962, B. D. Josephson advanced the theory of tunnel effects in super-
conductors (the Josephson effect), which was marked by the Nobel Prize in
1973.43 This discovery strongly intensified the experimental studies of super-
conductivity and was applied to electronics and computer technologies.
In 1986, K. A. Müller and J. G. Bednorz discovered HTSCs,2 which ini-
tiated the huge “boom” manifested in numerous publications. The principal
thought of Müller was the following: by selecting the suitable chemical com-
position, one can enhance the electron–phonon interaction and thus increase
the critical temperature Tc . Müller and Bednorz found a new class of HTSCs,
the so-called cuprate superconductors.
HTSCs have been studied for 23 years by making significant efforts, but
the pattern is still not completely clear. This is related to the complicated
structure of cuprates, difficulties in the production of perfect single crystals,
and difficult control over the degree of doping. The comprehension of HTSCs
will be attained if our knowledge about them approaches a critical level
sufficient for the understanding of a great amount of experimental data from
a single viewpoint. In Fig. 2.1, we present a plot of the critical temperatures
of superconductivity over years.44, 45
Fig. 2.1 The evolution of critical temperatures since the discovery of superconductivity.
Fig. 2.2 Elementary cell of yttrium ceramics YBa2 CuO7 : 1–yttrium; 2–barium; 3–
oxygen; 4–copper.
N 1 2 3 4 5 6
Fig. 2.5 Structure of energy bands near the Fermi energy for (a) high-temperature super-
conductor Bi2 Sr2 CaCu2 O8 and (b) low-temperature superconductor Bi2 Sr2 CaCuO6 .
We note the importance of the structure of zones near the Fermi level.
For example, the structures of the bands of two superconductors (high-
temperature Bi2 Sr2 CaCu2 O8 and low-temperature Bi2 Sr2 CaCuO6 ) are sig-
nificantly different at point K. The energy interval from the band bottom
to the Fermi surface is equal to 0.7 eV at point K for the high-temperature
superconductor and 0.1 eV for the low-temperature one. The explanation
for the influence of the multiband structure on Tc is given in Chap. 3.
It is worth noting that all copper-oxide superconductors have a block
character in their structure. The main block defining the metallic and super-
conductive properties of a compound is the plane with CuO2 which con-
tains the square lattice of copper ions coupled with one another through
oxygen ions. Depending on the composition, the elementary cell of a high-
temperature compound can have one, two, three or more cuprate layers.
Moreover, the critical temperature of the superconducting transition has a
nonmonotonous dependence on the number of cuprate layers.46
and spin correlations in the metallic phase became possible only with the
help of the neutron scattering.
“transfer system” in metals which is distinct from the phonon system. The
general scheme for the interaction of electrons via a transfer system X can
be presented as
e1 + X → e11 + X ∗ ,
(2.2)
e2 + X ∗ → e12 + X,
where ei corresponds to an electron with momentum pi , X is the ground
state and X ∗ is the excited state of the transfer system.
As a result of this reaction, the system returns to the initial state, and
the electrons make an exchange by momenta.
74 Modern Aspects of Superconductivity: Theory of Superconductivity
It can be shown that such an interaction leads to the attraction, and the
critical temperature is given by the formula
−1
Tc ∼ ΔE exp ,
λ
where ΔE is the difference in energies of the states X and X ∗ , and λ depends
on the interaction of electrons with the system X.51
the doping of HTSC compounds, there arise holes which can form the
complexes with spinons — holons. Superconductivity is explained by the
pairing of holons, i.e. by the creation of spinless bosons with a double
charge. That is to say, the pairing of carriers in the RVB model is realized
due to the exchange by magnons. At low temperatures, the paired holons
form a superconducting condensate. The RVB model has played a positive
role, by attracting the attention of researchers to the study of antiferromag-
netism in HTSCs, though spinons and holons have not been experimentally
identified.
The effective mass m̃ can be very high and can exceed the mass of a free
quasiparticle in the conduction band by several orders.
To calculate the superconducting transition temperature Tc which cor-
responds to the Bose condensation, the authors applied the formula for the
ideal Bose gas:
3.312 N 2/3
kB Tc = . (2.12)
m̃
At m̃ = 100 and N = 1021 cm−3 , this formula yields Tc ∼ 28 K.
Analogously to the bipolaronic model, the bisoliton mechanism of high-
temperature superconductivity was proposed in Ref. 45 and 46.
where Δo is the maximum value of the gap, and a is the lattice constant.
The gap is strongly anisotropic along the direction (110) in the k space. In
this case, the order parameter sign is changed in the directions along kx
and ky .
Apart from the d symmetry, it is worth considering also the s symmetry,
for which we can choose two collections of basis functions:
Δ(k) = Δo . (2.15)
They have eight parts with alternating signs, and eight nodes which are split
by ±γπ/2 along the direction (110). G. Kotliar and coworkers used mixed
s + id states63 :
[s + id] Δ(k) = Δo [ε + i(1 − ε)[cos(kx a) − cos(ky a)]]. (2.18)
R. B. Laughlin62 analyzed the mixed states dx2 +y2 + idxy on the basis of
the anyonic mechanism of pairing:
[d + id] Δ(k) = Δo [(1 − ε)[cos(kx a) − cos(ky a)]
+ iε[2 sin(kx a) sin(ky a)]], (2.19)
where ε is the share of s or dxy states mixed with the dx2 +y2 states, and εΔo
is the minimum value of the energy gap. These mixed states are of interest
because they are not invariant with respect to the inversion in the time.
The value and phase of the superconductor’s order parameter, as functions
of the direction in cuprate planes CuO2 , are given in Fig. 2.8 for various
kinds of pairing symmetry.
Fig. 2.8 Types of pairing symmetries which are considered for HTSCs.
The symmetric states of the system can be presented by the indication of all
possible subgroups of the full group, relative to which the order parameter is
invariant.66 The full symmetry group of a crystal includes a point symmetry
group G, operations of inversion of the time R, and the group of calibration
84 Modern Aspects of Superconductivity: Theory of Superconductivity
Fig. 2.9 Critical current in a Josephson junction vs the applied magnetic field: (a) stan-
dard tunnel contact, (b) corner tunnel junction for a superconductor with the s symmetry
of the order parameter, (c) corner tunnel junction for a superconductor with the d sym-
metry of the order parameter.67
Fig. 2.10 Critical current as a function of the magnetic field in a Josephson junction
YBCO-Au-Pb in two geometries: (a) standard and (b) corner.67
Here, Aij is the constant which characterizes the junction of ij; θi and
θj are the angles of the crystallographic axes with the boundary plane;
and Φij is the difference of phases of the order parameters on both sides
88 Modern Aspects of Superconductivity: Theory of Superconductivity
where Δo is the maximum value of the gap and a is the lattice constant.
The gap is strongly anisotropic in the direction (110) in the k space, and
the sign of the order parameter is changed along the directions kx and ky .
The main results of this section have been published in Refs. 78–81.
Fig. 2.12 Distribution of signs of the gap function in the limits of the Brillouin first band.
H = H0 + Hint , (2.29)
1 +
s(q) = Ψk+q σαβ Ψkβ (2.31)
2
α,β,k
where s(q) is the operator of spin density, σαβ is the Pauli matrix, Ψ+ k+q,α
is the operator of creation of an electron with the momentum k + q and
the spin projection α, Ψkβ is the operator of creation of a hole with the
momentum k and the spin projection β, and S(−q) is the operator of spin
fluctuations, whose properties are set by the correlator χ(q, ω)69, 70 ;
χij (n, m) the spin susceptibility which is modeled by
χQ
χ(q, ω) = , (2.32)
1+ ξ2 (q − Q)2 − iω/ωSF
qx > 0, qy > 0,
1
H0 = −t ψ+
α (n)ψα (n + p) + Si (n)χ−1
ij (n, m)Sj (m), (2.37)
n,p
2 n,m
σ
i
Hint = g ψ+
α (n) ψβ (n)Si (n), (2.38)
n
2 αβ
where the sum is taken over all the sites of the infinite lattice (the lattice
constant is equal to a), ρ is the unit vector joining the neighboring sites, Si
.
is the spin operator, N = n ψ+ α (n)ψα (n) is the operator of the number of
particles, t is the half-width of the conduction band69, 70 and χij (n, m) is
the spin correlation function.
It is necessary to calculate the grand partition function:
exp{−βΩ(μ, β, g)} ≡ Tr exp{−β(H − μN )}, (2.39)
1
β= ,
kT
where μ is the chemical potential, g is the coupling constant and Ω(μ, β, q)
is the thermodynamic potential.
It is convenient to use the formalism of continual integration for a sys-
tem of Fermi particles. The method of functional integration, i.e. the inte-
gration in the space of functions, was proposed by N. Wiener in 1925, but
the physicists paid no attention to this method. Continual integrals were
introduced into physics by R. Feynman4 in the 1940s and were used for the
reformulation of quantum mechanics. Continual integration is one of the
most powerful methods of contemporary theoretical physics, allowing one
to simplify, accelerate and clarify the process of analytic calculations. The
94 Modern Aspects of Superconductivity: Theory of Superconductivity
with
∂
L(τ) = ψ+
α (n, τ)− μ ψα (n, τ) − t ψ+
α (n, τ)ψα (n + p, τ)
n
∂τ n,p
σ
i
+g ψ+
α (n, τ) ψβ (n, τ)Si (n, τ)
n
2 α,β
1
+ Si (n, τ)χ−1
ij (n, m, τ)Sj (m, τ), (2.41)
2 n,m
where L(τ) is the Lagrangian of the system, and N is the normalizing factor:
) β
−1
N = dSi (n) dψα (n, τ) dψα (n, τ) exp −
+
dτL(τ, μ = g = 0) .
n 0
(2.42)
We will use the matrix formalism to construct the theory of perturbations
for the Green functions. It is convenient to introduce a four-component
bispinor (Majorana):
ψ
Ψ= , (2.43)
−σ2 ψ∗
ψ1
ψ= ,
ψ2
where σ1 are the Pauli spin matrices. The Majorana spinor is a Weyl spinor
written in the four-component form.
Now, we can write L in the form
1 0 ∂
L= Ψ(n, τ) Γ − Γ Γ5 μ δ̂ − t(τ − δ̂)Γ Γ5
0 0
(2.44)
2 n,p ∂τ
g
× ψ(n + p, τ) + Ψ(n, τ)Γi Γ5 Ψ(n, τ)S i (n, τ)
4 n
1
+ Si (n, τ)χ−1
ij (n, m, τ)Sj (m, τ), (2.45)
2 n,m
Physics of High-Tc Superconductors 95
with
0 I
Γ =0
, (2.46)
I 0
0 −σi
Γ =
i
,
σi 0
1 0
I= ,
0 1
First, we consider the case of a free system. By integrating over the fermion
fields in the functional integral, we obtain
1 1
F0 = βΩ(J) = − Tr Ln(G−1 0 − Γ Γ5 μ + J) + Tr Ln G0 ,
0 −1
(2.52)
2 2
−1
0 ∂
G0 = Γ − t(1 − δ̂)Γ Γ5
0
. (2.53)
∂τ
For the interacting system, the free energy takes the form
F = F0 + Fint . (2.54)
.
We represent Fint as a series in g: Fint = n=1 (g2 )Fn . In the lowest order
in g, we get the relation
g2
F1 = − Tr{Γi G T rΓj G − 2GΓi GΓj }χij , (2.55)
32
where the Green function for fermions and spins is as follows:
Ψ(x, τ)Ψ̄(y, τ ) = G(x, τ; y, τ ),
Si (x, τ)Sj (y, τ ) = χij (x, τ; y, τ ). (2.56)
Taking the condition δF/δG = 0 into account, we obtain the Schwinger–
Dyson equation:
g2 i
G−1 = G−1 0 − Γ 0
Γ5 μ + {Γ T rΓj Gχij − 2Γi GΓj χij }, (2.57)
12
where Γ0 and Γi are Dirac matrices (Fig. 2.13).
Here, the straight line corresponds to the Green function for fermions G,
and the wavy line to the Green function for spins χi,j (= δij χ). The equation
for free energy
g2
F1 = − Tr{Γi G TrΓj G − 2GΓi GΓj }χij , (2.58)
32
where G and χij are, respectively, the Green fermion and spin functions,
corresponds to the contribution of two vacuum diagrams (Fig. 2.14).
0
G−1 (k, iωn ) = G−1
0 (k, iωn ) − Σ(k, iωn ), (2.62)
and Eq. (2.57) is the equation for Σ. We write the solution for Σ in the form
Σ = Γ0 A + Γ0 Γ5 B + Δ + Γ0 Γ5 Γi Δi , (2.63)
Γ0 (ωn − A) − Γ0 Γ5 (ε(k) − μ + B) + Δ + Γ0 Γ5
G(k, iωn ) = . (2.64)
(ωn − A)2 − (ε(k) − μ + B)2 − Δ2 − Δ2i
98 Modern Aspects of Superconductivity: Theory of Superconductivity
where
Taking into account that Im χ(q, −ω) = −Im χ(q, ω), Eq. (2.70) can be
reduced to
1 ∞ 1 1
χ(q, iωn ) = − dω Im χ(q, ω ) − . (2.73)
π 0 iωn − ω iωn + ω
With the help of the formula (2.72), we can rewrite (2.67) in the form
∞
g2 d2 p
Δ(k, iωn ) = dωF (p, ω) dω Im χ(k − p, ω )
8π(2πi) (2π)2 C 0
1 1 βω
×
−
th , (2.74)
iωn − ω − ω iωn − ω + ω 2
where
Δ(k, iωn)
F (p, ω) = .
ω2 − (ε(p) − μ)2
100 Modern Aspects of Superconductivity: Theory of Superconductivity
ω2 − (ε(p) − μ)2 − δ2
Re F (p, ω + iδ) = ReΔ(p, ω) . (2.76)
[ω2 − (ε(p) − μ)2 − δ2 ]2 + ω2 δ2
If the superconducting gap depends weakly on the frequency, then
Δ(k, ω) Δ(k, 0). In this case, we set ω = 0 in Eq. (2.75). Considering
the relation Reχ(q, ω) = Re χ(q, −ω) and making the simple transforma-
tions, we obtain an equation similar to that in Ref. 70:
g2 π/a d2 p th β(ε(p)−μ)
Δ(k) = Re χ(k − p, ε(p) − μ) 2
8 −π/a (2π)2 ε(p) − μ
∞
dυ βυ
+ 2 cth Im χ(k − p, υ)
0 π 2
(ε(p) − μ)2 − υ2 + δ2
× Δ(p). (2.77)
[(ε(p) − μ)2 − υ2 + δ2 ]2 + 4υ2 δ2
It is important that our method of calculations yields a more general
equation for the superconducting gap than that obtained by D. Pines. Our
equation coincides with the Pines equations after some simplification, which
is the test for our calculations.
Physics of High-Tc Superconductors 101
where G−1
o (k, iωn ) = r iωn − r r5 ε(k), and G(k, iωn ) is given by the for-
o o
β
d2 k 2 ch 2 (ε − μ)2 + δ2
Ω(Δ) = V − ln
(2π)2 β ch βε
2
Δ2 β (ε − μ)2 + Δ2
+ th . (2.85)
2 (ε − μ)2 + Δ2 2
By making some transformations and calculating the traces of Γ matrices,
we arrive at the equation
V dw d2 k 4
Ω(Δ) − Ω(0) = Δ (k, w)
2 C 2πi (2π)2 i
1 1
× . (2.86)
eBw + 1 [w − (ε(k − μ)2 ]2
2
F = βΩ; Ω(Δ) and Ω(0) are the thermodynamic potentials at T < Tc and
T > Tc , respectively, and V is the two-dimensional volume (the area of a
cuprate layer).
It is easy to verify that despite the presence of the factors [ε(k) − μ] in
the denominator of the formula (2.87), no singularity on the Fermi surface
ε(k) − μ is present.
Equations (2.85) and (2.87) (at T ∼ Tc ) can be a basis for calculations of
various thermodynamic quantities, including a jump of the heat capacity.
A similar method of calculations was developed by S. Weinberg.85
β
d2 k 2 ch 2 ( (ε(k) − μ)2 + Δ(k)2 )
Ω(Δ) = V − ln
(2π)2 β ch βε
2
Δ(k)2 (ε(k) − μ)2 + Δ(k)2
+ th . (2.88)
2 (ε(k) − μ)2 + Δ(k)2 2
Here, V is the two-dimensional volume (the area of a cuprate layer), Δ(k)
is the superconductor gap, k is the momentum of an electron, and ε(k) =
−2t[cos(kx a) + cos(ky a)] gives the spectrum of two-dimensional electrons.
The heat capacity is calculated by the formula
∂2Ω
C = −T . (2.89)
∂T 2
The results of calculations of the heat capacity C are given in Ref. 45.
We have carried out computer-based calculations of the heat capacity of
the compound YBa2 Cu3 O6.63 . The equation for the superconductor gap Δ
was obtained in Sec. 2.7.5. The task of solving the integral equation was
reduced to that of an algebraic equation which was solved by the method
of iterations.
The equation for the superconductor gap depends on the spin correla-
tion function. For it, the necessary data were taken from Refs. 69 and 70.
In calculations, we used the following values of parameters of the correlator:
ωSF (TC ) ≈ 7.7 meV, χS (Q) = 44 eV, ξ/a ≈ 2.5, t = 2 eV, μ = 0.25 and
Tc = 95 K. We chose a value of the constant g which satisfies the relation
2Δ/kTc = 3.4, similar to that in Refs. 66 and 67. The results of numerical
104 Modern Aspects of Superconductivity: Theory of Superconductivity
calculations are presented in Fig. 2.15, where we can see the temperature
dependence of the electron heat capacity: 1 — a curve which is an approxi-
mation of the results of computer-based calculations (crosses); 2 — a curve
which describes the exponential BCS dependence. These results are in good
agreement with experiments made by A. Kapitulnik and co-workers86 at low
temperatures.
The calculations indicated a square dependence of the heat capacity on
the temperature in the temperature range from zero to the superconducting
temperature. It was shown45,46,79,80 that the dependence has the form
C = AT 2 , (2.90)
where A = 0.126 J/K3 mol.
Fig. 2.16 Temperature dependence of the electron heat capacity: 1 — a curve which
is an approximation of the results of computer-based calculations (crosses); 2 — a curve
which describes the exponential BCS dependence, with the points giving the experimental
data.86
Fig. 2.17 The parameter R(Tc ), which characterizes the heat capacity jump as a function
of the doping for Y0.8 Ca0.2 Ba2 Cu3 O7−x . The solid curve presents the results of theoretical
calculations; the points mark the experimental data.89
106 Modern Aspects of Superconductivity: Theory of Superconductivity
urgent at the present time. Many experiments have confirmed the d sym-
metry of the pairing.67, 68 In particular, we mention new experiments,68 in
which the researchers have studied the quantization of magnetic flows in a
ring which includes three Josephson junctions. These works supporting the
existence of d pairing in HTSCs were marked by the Barkley Prize. It was
also shown that the sign of the order parameter in YBa2 Cu3 O7−δ depends
on the direction. This dependence corresponds to the dx2 −y2 symmetry.
The temperature dependence of the electron heat capacity obtained in
our calculations79, 80 is related to the d pairing. These thermodynamic cal-
culations can be a supplementary test in the determination of the pairing
symmetry in HTSCs.
It should be noted that thermodynamic calculations clarifying the behav-
ior of HTSCs have also been performed in Refs. 92–95 in the frame of the
other mechanisms of pairing.
2.8 Summary
High-temperature superconductivity is a dynamical field of solid-state
physics which has been intensively developed by theorists and experi-
menters. By summarizing the physical properties and the mechanisms of
superconductivity in new HTSCs , we wish to separate the main properties
and the theoretical problems arising in the studies of HTSCs.
In order to comprehend the nature of the superconducting state, it is nec-
essary to construct a consistent microscopic theory which would be able to
describe superconductive and normal properties of HTSCs. It is seen from
the above-presented survey that many mechanisms of pairing in HTSCs,
proposed as the explanations for this phenomenon, have been advanced. On
the whole, the most probable seems to be the “synergetic” mechanism of
pairing, whose constituents are the electron–phonon interaction, spin fluc-
tuation, and other types of interaction in cuprate planes.
We believe that the known challenging properties — which are contra-
dictory to a certain extent — of many chemical compounds in the super-
conducting and normal phases can be explained only by considering the
interaction of all the degrees of freedom, such as lattice-, electron- and spin-
related ones. In this case, it is also necessary to take into account the com-
plicated structure of HTSCs. Further development of the theory will require
not only the execution of bulky numerical calculations, but also the solu-
tion of a number of fundamental problems concerning the strong electron
Physics of High-Tc Superconductors 107
Multiband Superconductivity
3.1 Introduction
In this chapter, we present the theory of superconductivity with regard
to a complicated multiband structure of superconductors. Calculations of
the band structure of cuprate superconductors indicate that several energy
bands intersect one another on the Fermi surface in these compounds,96 and
the Fermi surface passes through high-symmetry points which correspond
to the Lifshitz electronic topological transition. In addition, the discovery of
two-gap superconductivity in two-band superconductors MgB2 allows one
to consider the possibility of using a multiband theory of superconductivity.
In Sec. 3.2, we analyze the problems of multiband superconductivity and
superconductivity at room temperature. In Sec. 3.3, we study the physical
properties of superconductors MgB2 and problems of two-gap supercon-
ductivity, as well as the phase diagram for superconductors MgB2 on the
basis of the renormalization group approach. In the field of superconduc-
tivity, we meet the problem-maximum — the creation of room-temperature
superconductors. We consider this problem in our book and give some rec-
ommendations on the search for these superconductors. The results of this
chapter were obtained by the authors and published in Refs. 96–103.
The theory of superconductivity arising from the electron–phonon inter-
action mechanism by Bardeen, Cooper and Schrieffer (BCS)104 has been
well established, and it is now the standard theory for superconductiv-
ity.2, 3, 6, 15, 16, 104 In Chap. 2, we have considered non-BCS mechanisms via
spin fluctuations, charge fluctuations (plasmons) and electron excitations
(excitons), which have attracted great interest, for example in relation to the
possibility of high-Tc superconductivity. These mechanisms have the com-
mon characteristic that the electron–electron (e–e) interaction is an origin
108
Multiband Superconductivity 109
3.2.1 Hamiltonian
We start from the Hamiltonian for two bands i and j :
H = H0 + Hint , (3.1)
with
H0 = [εi − μ]a+
ikσ aikσ + [εj − μ]ajkσ ajkσ ,
+
(3.2)
k,σ
1
Hint = αβγδ aip1 α aip2 β aip3 γ aip4 δ + (i → j)
Γiiii + +
4
δ(p1 +p2 ,p3 +p4 )αβγδ
+ Γiijj
αβγδ aip1 α aip2 β ajp3 γ ajp4 δ + (i → j)
+ +
+ Γijij + +
αβγδ aip1 α ajp2 β aip3 γ ajp4 δ + (i → j) , (3.3)
with
ip1 αjp2 β|kp3 γlp4 = dr1 dr2 φ∗ip1 α (r1 )φ∗jp2 β (r2 )
Γiijj jjii
αβγδ = Γαβγδ = g2 (δαδ δβγ − δαγ δβδ ),
Γijij jiji
αβγδ = Γαβγδ = g3 δαδ δβγ − g4 δαγ δβδ , (3.10)
where
∂aυkσ (t)
i = [aυkσ , H]. (3.13)
∂t
The equation for the Green function (3.12) is rewritten after inserting
Eq. (3.13) as
∂1 2
i −iT aυkσ (t)a+ υkσ (t )
∂t
1 2
= δ(t − t ) + [ευ − μ] −iT aυkσ (t)a+
υkσ (t )
1
+ g1 (δσδ δβγ − δσγ δβδ )Gγδβσ
2υυυυ (p3 , p4 , p2 , k; t, t )
2
βγδ δ(k+p2 ,p3 +p4 )
1
+ g2 (δσδ δβγ − δσγ δβδ )Gγδβσ
2υ υ υυ (p3 , p4 , p2 , k; t, t )
2
υ βγδ δ(k+p2 ,p3 +p4 )
1
+ (g3 δσδ δβγ − g4 δσγ δβδ )Gγδβσ
2υυ υ υ (p3 , p4 , p2 , k; t, t ),
2
υ βγδ δ(k+p2 ,p3 +p4 )
(3.14)
where
1 2
Gγδβσ
2υυυυ (p3 , p4 , p2 , k; t, t ) = −iT aυp3 γ (t)aυp4 δ (t)aυp2 β (t − 0)aυkσ (t ) ,
+ +
(3.15)
1 2
Gγδβσ
2υ υ υυ (p3 , p4 , p2 , k; t, t ) = −iT aυ p3 γ (t)aυ p4 δ (t)aυp2 β (t − 0)aυkσ (t ) ,
+ +
(3.16)
1 2
Gγδβσ
2υυ υ υ (p3 , p4 , p2 , k; t, t ) = −iT aυp3 γ (t)aυ p4 δ (t)aυ p2 β (t − 0)aυkσ (t ) .
+ +
(3.17)
where
2 − nυk+q − nυk −q
K(ω, k, k , q) = . (3.20)
2ω − ευk+q − ευk −q
where
√
N (εf ) = 2πm∗υ m∗υ ε 2− nυk+q − nυk −q , (3.22)
ε=εf
1
nυk = . (3.23)
exp ευk − εf /T + 1
Therefore, we obtain the equation for coupled states in the electron system:
υ υ 2Δ 2Δ
2
1 − g2 N εf N εf ln 1 − 1− = 0. (3.27)
aυ aυ
From Eq. (3.27), we find the possibility of the existence of a solution for
the coupled states, if g2 = 0. Thus, the effective e–e interaction g 2 with a
positive value contributes to the superconductivity.
Δ Δ
1 − (g1 + g2 )N (εf ) ln1 − 1 − (g1 − g2 )N (εf ) ln1 − = 0.
a a
(3.29)
When g1 + g2 < 0 or g1 − g2 < 0, we can find solutions to this equation.
expect that the transition temperature is higher (300 K) than that for copper
oxides by the cooperative mechanism.
Let us discuss the problem of superconductors operating at room tem-
perature (RTSCs). It is obvious that the main task in the field of super-
conductivity is the fabrication of materials with superconducting properties
at room temperature. The study of HTSCs is only a stage on the way to
the main purpose, namely to the development of RTSCs. At the present
time, the highest known Tc = 135 K (at the atmospheric pressure). Now
there occurs a wide-scale search for such superconductors. Since none of the
known physical laws allows one to exclude the possibility for RTSCs to exist,
the future discovery of RTSCs seems to be without doubt. New materials
are created by means of physical and chemical modifications of the known
compounds, including the application of nanotechnological methods. The
creation of new materials with preassigned physical properties is one of the
actual problems of modern science. At present, various approaches are being
developed in order to avoid the labor-consuming sorting of different chem-
ical compounds and conditions of synthesis, i.e. to optimize the solution of
this problem. One of the most efficient and promising solution is the method
of structural design. In Ref. 121, the structural design was used in solving
the problems related to the search for new HTSCs on the basis of complex
copper oxides.
We now present the recommendations concerning the search for new
HTSCs with higher critical temperatures which were advanced by the Nobel
118 Modern Aspects of Superconductivity: Theory of Superconductivity
has a hexagonal structure118 ; see Fig. 3.3.23 The results of calculations of the
temperature dependence of the specific heat of MgB2 given in Fig. 3.5 follow
the corresponding experimental data. The band structure of MgB2 has been
calculated in several works since the discovery of superconductivity.128 The
electronic properties of MgB2 are plotted in Fig. 3.4. The band structure of
MgB2 is similar to that of graphite and is formed by the σ and π zones.
Fig. 3.3 The structure of MgB2 containing graphite-type B layers separated by hexagonal
close-packed layers of Mg.
T (K)
Fig. 3.5 The MgB2 specific heat vs. temperature. The line indicates theory, and the
circles represent the experimental data.
e–p interaction mechanism for the superconductivity. In this case, the pos-
sibility of two-band superconductivity has also been discussed in relation
to two-gap functions experimentally and theoretically. Very recently, two-
band or multiband superconductivity has been theoretically investigated
in relation to the superconductivity arising from Coulomb repulsive inter-
actions. The two-band model was introduced by Kondo.129 We have also
investigated anomalous phases in the two-band model by using the Green
function techniques.96–99, 101, 103 Recently, we have pointed out the impor-
tance of multiband effects in high-Tc superconductivity.96–98 The expres-
sions for the transition temperature for several phases have been derived,
and the approach has been applied to the superconductivity in molecular
crystals by charge injection and the field-induced superconductivity.13 In
previous papers,96–99, 101, 103 we have investigated the superconductivity by
using the two-band model and the two-particle Green function techniques.
We have applied the model to the e–p mechanism for the traditional BCS
method, the e–e interaction mechanism for high-Tc superconductivity,2 and
the cooperative mechanism. In the framework of the two-particle Green
function techniques,101 it has been shown that the temperature dependence
of the superconductivity gap for high-Tc superconductors is more compli-
cated than that predicted in the BCS approach. In Ref. 99, phase diagrams
for the two-band model superconductivity have been investigated, by using
the renormalization group approach. Below, we will discuss the possibil-
ity of the cooperative mechanism of two-band superconductivity in relation
to high-Tc superconductivity and study the effect of the increase of Tc in
MgB2 due to the enhanced interband pairing scattering. In this section,
we will investigate our two-band model for the explanation of the multi-
gap superconductivity of MgB2 . We apply the model to the e–p mechanism
for the traditional BCS method, the e–e interaction mechanism for high-Tc
superconductivity, and the cooperative mechanism in relation to multiband
superconductivity.
Tc (K)
Fig. 3.6 The influence of doping on Tc in Mg1−x Alx B2 . The line indicates theory, and
the circles represent the experimental data.
effects: the first one is a coupling effect related to the changes of the car-
rier concentration, and the second one depends on the introduction of new
scattering centers leading to a modification of the interband and intraband
scattering.
H = H0 + Hint , (3.30)
with
H0 = [i (k) − μ]a†ikσ aikσ + [j (k) − μ]a†jkσ ajkσ , (3.31)
k,σ
1 † †
Hint = αβγδ aip1 α aip2 β aip3 γ aip4 δ + (i ↔ j)
Γiiii
4
δ(p1 +p2 ,p3 +p4 ) α,β,γ,δ
† †
+ Γiijj
αβγδ aip1 α aip2 β ajp3 γ ajp4 δ + (i ↔ j)
† †
+ Γijij
αβγδ aip1 α ajp2 β aip3 γ ajp4 δ + (i ↔ j) , (3.32)
Γijkl
αβγδ = ip1 αjp2 β|kp3 γlp4 δδαδ δβγ − ip1 αjp2 β|lp4 δkp3 γδαγ δβδ , (3.33)
Multiband Superconductivity 123
with
where gi1 and gj1 represent the ith and jth intraband two-particle nor-
mal scattering processes, respectively, and g2 indicates the intraband two-
particle umklapp scattering (see Fig. 3.7). For simplicity, we consider the
three cases in Refs. 99 and 100: (i) g1 = 0 and others = 0; (ii) g2 = 0; and
others = 0; (iii) gi1 and gj1 , and others = 0, using the two-particle Green
function techniques (see Fig. 3.7).
It was shown that, possibly, the two-gap superconductivity arises in
case (iii). The superconductivity arising from the e–p mechanism (g1 < 0
and g1 < g2 ) such as for MgB2 is in the two-gap region. On the other hand,
the superconductivity of copper oxides (g1 > g2 ) is outside the two-gap
region. These results predict that we may observe two-gap functions for
MgB2 and only a single-gap function for copper oxides.
124 Modern Aspects of Superconductivity: Theory of Superconductivity
Fig. 3.7 Electron–electron interactions. Solid and dashed lines indicate π and σ bands,
respectively; gi1 , gj1 , g2 , g3 and g4 contribute to superconductivity.
H = H0 + Hint , (3.39)
where
H0 = [i − μ]a†ikσ aikσ + [j − μ]a†jkσ ajkσ , (3.40)
k,σ
Hint = g1i a†ik a†i−k ai−k aik + i→j+ g2 a†ik a†i−k aj−k ajk . (3.41)
We now define the order parameters which are helpful in constructing the
mean-field Hamiltonian:
1 † † 2
Δi = aip↑ ai−p↓ , (3.42)
p
1 2
Δj = a†jp↑ a†j−p↓ . (3.43)
p
where
μ−Ec
dξ (ξ2 + Δ2i )1/2
fi = tan h ,
μ (ξ2 + Δ2i )1/2 2T
(3.45)
μ−Ej (ξ2 + Δ2j )1/2
dξ
fj = tan h ,
μ−Ec (ξ2 + Δ2j )1/2 2T
We have tried to estimate the coupling constant of the pair electron scat-
tering process between the π and σ bands of MgB2 and have calculated
the parameters by using a rough numerical approximation. We focus on one
π band and σ band of MgB2 and consider electrons near Fermi surfaces.
We found that the parameter g1 = −0.4 eV, by using the transfer integral
between the π and σ bands. We estimate the coupling parameter g2 of the
pair electron scattering process by the expression
1,2
g2 = Vk1,k2 , (3.47)
k1,k2
1,2
Vk1,k2 = u∗1,r (k1)u∗1,s (k1)vrs u2,t (k2)u2,u (k2), (3.48)
r,s,t,u
where the labels 1 and 2 are the π band and the σ band, respectively, ui,r (ξ)
is the LCAO coefficient for the ith band and ξ is the moment.115, 116 The
indices k1 and k2 are summed over each Fermi surface. However, it is diffi-
cult to perform the sum exactly. In this case, we used a few points near the
Fermi surface. The coupling constant of the pair electron scattering between
the π band and the σ band is g2 = 0.025 eV. From numerical calculations
of Eqs. (3.44)–(3.46), we can also obtain the temperature dependence of
the two-gap parameters (see Fig. 3.8). We have used the density of states
of the π and σ bands (ρi = 0.2 eV −1 ; ρj = 0.14 eV−1 ), chemical potential
μ = −2.0, the top energy of the σ band Ej = −1.0 and the fitting param-
eters (gi1 = −0.4 eV; gj1 = −0.6 eV; g2 = 0.02 eV). These calculations have
qualitative agreement with experiments.118, 130, 131 The expression for the
transition temperature of superconductivity derived in a simple approxima-
tion is:
−1
Tc+ = 1.13(ζ − Ej ) · exp , (3.49)
g+ ρ
126 Modern Aspects of Superconductivity: Theory of Superconductivity
where
1
g+ = (B + B 2 − 4A) (3.50)
24
and
A = g1i g1j − g22 , (3.51)
B = g1i + g1j , (3.52)
ζ = −μ. (3.53)
From the expressions for Tc+ , we can see the effect of increase of Tc+ due to
the enhanced interband pairing scattering (g2 ).
Figure 3.9 shows a schematic diagram of the mechanism of pairing for
two gaps. The scenario is as follows. Electrons from the π and σ bands
make up the subsystems. For g2 , we have two independent subsystems with
the different transition temperatures of superconductivity Tcπ and Tcσ and
two independent superconducting gaps. In our model, we have two coupled
superconducting gaps with the relation (3.44) and one transition tempera-
ture of superconductivity Tc+ , which is in agreement with experiments. In
this model, we have two channels of superconductivity: conventional chan-
nel (intraband g1 ) and unconventional channel (interband g2 ). Two gaps
appear simultaneously in different bands which are like BCS gaps. The gap
Multiband Superconductivity 127
Fig. 3.9 Schematic diagram of the mechanism of pairing for two gaps.
in the π band is bigger than that for the σ band, because the density of
state is 0.25 eV in the π band and 0.14 eV in the σ band. The current of
Cooper pairs flows from the π band into the σ band, because the density
of Cooper pairs in the π band is much higher. The tunneling of Cooper
pairs also stabilizes the order parameter in the σ band. In this way, we can
predict the physical properties of the multigap superconductivity if we have
the superconductors with a multiband structure, as shown in Fig. 3.9.
Thus, we have presented our two-band model with the intraband two-
particle scattering and interband pairing scattering processes to describe
two-gap superconductivity in MgB2 . We defined the parameters of our model
and made numerical calculations of the temperature dependence of two
gaps in qualitative agreement with experiments. We have proposed a two-
channel scenario of superconductivity: the conventional channel (intraband
g1 ), which is connected with the BCS mechanism in different bands; and
the unconventional channel (interband g2 ), which describes the tunneling of
Cooper pairs between two bands. The tunneling of Cooper pairs also sta-
bilizes the order parameters of superconductivity and increases the critical
temperature of superconductivity.
role of the e–p interaction has been pointed out for the superconductiv-
ity of MgB2 . Recent band calculations of the transition temperature for
MgB2 127, 128 with the McMillan formula132 have supported the e–p interac-
tion mechanism for the superconductivity. In this case, the possibility of the
two-band superconductivity in relation to two-gap functions has also been
considered experimentally and theoretically. Very recently, the two-band or
multiband superconductivity has been theoretically investigated in relation
to the superconductivity arising from Coulomb repulsive interactions.129
The two-band model was introduced by Kondo.129 Recently, we have pointed
out the importance of multiband effects in high-Tc superconductivity.99, 138
We have also investigated anomalous phases in the two-band model by using
the Green function techniques.101, 136 The expressions of the transition tem-
perature for several phases have been derived, and the approach has been
applied to the superconductivity in several crystals by charge injection and
the field-induced superconductivity.136, 139 In the previous section, we have
investigated superconductivity by using the two-band model and the two-
particle Green function techniques.99, 101 We have applied the model to the
e–p mechanism for the traditional BSC method, the e–e interaction mecha-
nism for high-Tc superconductivity and the cooperative mechanism. In the
framework of the two-particle Green function techniques,103 it has been
shown that, in the e–p system, a class of new so-called coupled states arises.
In this section, we investigate two- or multiband effects in superconductivity
by using the two-band model within the renormalization group approach.
Renormalization equations for the two-band superconductivity are derived
from the response function and the vertex correction of the model. Phase
diagrams numerically obtained from the renormalization equations are pre-
sented. We also discuss the superconductivity arising from the e–e repulsive
interaction in relation to the two-band superconductivity.
1 † †
Hint = αβγδ aip1 α aip2 β aip3 γ aip4 δ + (i ↔ j)
Γiiii
4
δ(p1 +p2 ,p3 +p4 ) α,β,γ,δ
† †
+ Γiijj
αβγδ aip1 α aip2 β aj p3 γ aj p4 δ + (i ↔ j)
† †
+ Γijij
αβγδ aip1 α aj p2 β aip3 γ aj p4 δ + (i ↔ j) , (3.56)
Γijkl
αβγδ = ip1 αjp2 β|kp3 γlp4 δ δαδ δβγ − ip1 αjp2 β|lp4 δkp3 γ δαγ δβδ , (3.57)
with
where gi1 and gj1 represent the ith and jth intraband two-particle nor-
mal scattering processes, respectively, and g2 indicates the intraband two-
particle umklapp scattering. Note that the Γ’s are given by
Γjjjj
αβγδ = gj1 (δαδ δβγ − δαγ δβδ ),
Γiijj jjii
αβγδ = Γαβγδ = g2 (δαδ δβγ − δαγ δβδ ), (3.62)
i i j j
=gi1 =gj1
i i j j
i j j i
=g2 =g2
i j j i
g̃2 g2 −gi1 g2 −g2 gj1 Li
= + , (3.66)
g̃2 g2 −gj1 g2 −g2 gi1 Lj
with
T
Li = Πi (k, ω) = Gi (q, ω )Gi (k − q, ω − ω ),
N
q,ω
T
Lj = Πj (k, ω) = Gj (q, ω )Gj (k − q, ω − ω ), (3.67)
N
q,ω
Fig. 3.11 Diagrams of the first-order vertex correction. (a) and (b) contribute to g̃i1 , (c)
and (d) are diagrams for g̃j1 , and (e)–(h) for g̃2 .
132 Modern Aspects of Superconductivity: Theory of Superconductivity
where G and T are the temperature Green function and temperature, respec-
tively. For the special case ω = 0, k = 0, the functions Li and Lj become
ui /y u /y
Li = − tan h + tan h i ln ξ − 2A, (3.68)
2ξ 2ξ
uj /y uj /y
Lj = − tan h + tan h ln ξ − 2A, (3.69)
2ξ 2ξ
where
A= dx ln x sech2 x, (3.70)
ui (uj ) and ui (uj ) are dimensionless functions expressed by the chemical
potential, the cutoff energy, the top energy of the jth band, and the density
of state for the ith (jth) band.
Next, we consider a first-order response function for a singlet Cooper
pair, as shown in Fig. 3.12. Then the first-order vertex function Λ for the
ith and jth bands can be written as
Λi + Λj 2 −g1 −g2 Li + Lj
= + . (3.71)
Λi − Λj 0 −g1 g2 Li − Lj
Fig. 3.12 Diagrams of the first-order response function. (a) and (b) contribute to Λi ,
and (c) and (d) show diagrams for Λj .
Multiband Superconductivity 133
In a similar way, using Eqs. (3.64) and (3.71), we obtain the differential
equations written as
∂
ln Λ+ = −g̃1 − g̃2 , (3.74)
∂x
∂
ln Λ− = −g̃1 + g̃2 , (3.75)
∂x
where Λ+ = Λi + Λj and Λ− = Λi − Λj .
(a) (b)
4 4
2 2
2
0
2
g
g
-2 -2
-4 -4
-5 0 5 -5 0 5
g g
1 1
g2 < −g1 with g2 > 0 (two-gap region), we can expect that two-gap func-
tions are observed. In the former region, those superconducting gaps may
be expressed by |Λ+ | > |Λ− |, and the latter may be |Λ+ | < |Λ− |. On
the other hand, we expect only a single-gap function in the other region.
These results agree with the previous solutions96, 99 derived by using the
two-particle Green function techniques. The superconductivity arising from
the e–p mechanism (g1 < 0 and |g2 | < |g1 |) such as that in MgB2 is in the
two-gap region. On the other hand, the superconductivity such as in copper
oxides (|g2 | > |g1 |) is outside the two-gap region. These results predict that
we may observe two-gap functions for MgB2 and only a single-gap function
for copper oxides. In the two-band model for negative g1 with transferring
or tunneling of Cooper pairs between two bands, we can expect that the
transition temperature becomes higher than that derived from the single-
band model. The tunneling of Cooper pairs causes stabilization of the order
parameter of superconductivity.119, 120, 129 We can also expect higher Tc of
the superconductivity than that for copper oxides in two regions (g1 < 0,
g2 < 0 and g1 < 0, g2 > 0) by the cooperative mechanism. Phase diagrams
for CDW, SDW, and singlet superconductivity derived from a more general
Hamiltonian will be presented elsewhere.
Thus, we have derived the renormalization equations and presented the
phase diagrams for the two-band superconductivity. In the framework of the
two-band model, the present results predict that superconductivity appears
even if the e–e interaction is positive. We can expect that the transition
temperature becomes higher than that of copper oxides by the cooperative
mechanism.
CHAPTER 4
Mesoscopic Superconductivity
4.1 Introduction
Recent advances in nanoscience have demonstrated that fundamentally new
physical phenomena are found when systems are reduced in size to dimen-
sions which become comparable to the fundamental microscopic lengths of
the investigated material. Superconductivity is a macroscopic quantum phe-
nomenon, and it is therefore especially interesting to see how this quantum
state is influenced when the samples are reduced to nanometer sizes. Recent
developments in nanotechnologies and measurement techniques allow exper-
imental investigation of the magnetic and thermodynamic superconducting
properties of mesoscopic samples in this regime.
In this book, we develop some theoretical models for describing such
nanoscale superconducting systems and explore possible new experimental
phenomena which we can predict based upon these models. In bulk sam-
ples, the standard BCS theory gives a good description of the phenomenon
of superconductivity.104 However, it was noticed by Anderson in 1959 that,
as the size of a superconductor becomes smaller, and the quantum energy
level spacing of the electrons in the sample approaches the superconduct-
ing gap, the BCS theory will fail.6 The exact solution to the reduced BCS
Hamiltonian for finite-size systems was developed by Richardson in the con-
text of nuclear physics a long time ago.149 This shows that, while the grand
canonical BCS wave function gives a very accurate solution to the BCS
Hamiltonian in the limit where the number of electrons is very large N 1,
for small values of N one has to use exact analytical methods to obtain
reliable results. The recent experimental advances in fabricating and mea-
suring superconductivity in ultrasmall mesoscopic and nanoscale grains has
renewed theoretical interest in the Richardson solution.150–153 In this book,
135
136 Modern Aspects of Superconductivity: Theory of Superconductivity
Here, a†jσ (ajσ ) and b†jσ (bj,σ ) are the creation (annihilation) operator in sub-
levels 1 and 2 with spin σ and energies ε1j and ε2j , respectively, the operators
for each sublevel satisfy the anticommutation relations, and the operators
between sublevels are independent, μ is the chemical potential, the second
term of Eq. (4.1) is the interaction Hamiltonian, g1 and g2 are the effective
interaction constant for sublevels 1 and 2, and g12 is an effective interac-
tion constant which corresponds to the pair scattering process between two
bands. The sums of j and k in Eq. (4.3) are over the set I of N1I states
which correspond to the half-filled band 1 with fixed width 2ω1D and the
set J of N2J states for band 2, respectively.
In this study, we assume that the Debye energies for two sublevels are the
same: ω1D = ω2D = ωD . Within this assumption, N1I and N2J are relatively
estimated by the density of state (DOS) for two bands as N1I /N2J = ρ1 /ρ2 ,
where ρ1 and ρ2 are the DOS for two bands. The interaction constants
g1 and g2 can be written as d1 λ1 and d2 λ2 , respectively. d1 = 2ωD /N1I
and d2 = 2ωD /N2J are the mean energy level spacings, and λ1 and λ2 are
the dimensionless parameters for two sublevels. We take the intersublevel
Mesoscopic Superconductivity 139
√
interaction constant g12 = d1 d2 λ12 . We then obtain the relation ρ1 /ρ2 =
N1I /N2J = d2 /d1 .
1 ∗ ∗
2 [g2 Δ1 + g1 Δ2 + g12 (Δ1 Δ2 + Δ1 Δ2 )],
2 2
+ (4.9)
g1 g2 − g12
where 1j = (ξ21j + Δ21 )1/2 and 2k = (ξ22k + Δ22 )1/2 . In Eq. (4.9), the values
of Δ1 and Δ2 must be chosen in such a way as to minimize Ω. From the
minimization of Ω, we obtain a coupled gap equation at zero temperature
for the two-gap system:
g . 1 .
−g12 k 212k Δ
Δ1 1 j 21j 1
= . 1 . 1 . (4.10)
Δ2 −g12 j g2 k Δ2
21j 22k
From the coupled gap equation (4.10), we formally obtain an expression for
the bulk gap for two-gap superconductivity at zero temperature:
−1 1
Δ̃1 = ω sinh , (4.11)
η1
−1 1
Δ̃2 = ω sinh , (4.12)
η2
where
1 λ2 + α± [η1 , η2 ]λ12
= , (4.13)
η1 λ1 λ2 − λ212
1 λ1 + α−1
± [η1 , η2 ]λ12
= , (4.14)
η2 λ1 λ2 − λ212
with
sinh η11
α± [η1 , η2 ] = ± . (4.15)
sinh η12
For the two-band superconductivity, we can consider two cases for the phase
of the gaps: sgn(Δ̃1 ) = sgn(Δ̃2 ) and sgn(Δ̃1 ) = −sgn(Δ̃2 ). For the same
phase, α+ is used in Eqs. (4.13) and (4.14), and we use α− for the opposite
phase. Note that Δ̃1 = −Δ̃2 in the limit of strong intersublevel coupling
λ12 , i.e. the opposite phase. For λ12 = 0, we find the same results for two
bulk gaps derived from the conventional BCS theory for two independent
sublevels.
Mesoscopic Superconductivity 141
C
EN 1 ,b1 ;N2 ,b2
(λ1 , λ2 , λ12 ) = EN
G
1 ,b1 ;N2 ,b2
(0, 0, 0) − EN
G
1 ,b1 ;N2 ,b2
(λ1 , λ2 , λ12 )
− n1 λ1 d1 − n2 λ2 d2 , (4.16)
where EN G
1 ,b1 ;N2 ,b2
(λ1 , λ2 , λ12 ) is the ground state energy of the (N1 + N2 )-
electron system. From Eqs. (4.4) and (4.9), the condensation energy of the
two-sublevel system can be expressed by the condensation energy of inde-
pendent single-level systems:
λ212
C
EN 1 ,b1 ;N2 ,b2
(λ1 , λ2 , λ12 ) = EN
C
1 ,b1
C
(λ1 ) + EN 2 ,b2
(λ2 ) −
λ1 λ2 − λ212
2
Δ1 Δ22 2(Δ∗1 Δ2 + Δ1 Δ∗2 )
× + + √ , (4.17)
d1 λ1 d2 λ2 d1 d2 λ12
where EN C C
(λ1 ) and EN (λ2 ) correspond to the condensation energy for
1 ,b1 2 ,b2
the single-band case. In the same phases of Δ1 and Δ2 , the condensation
energy (4.17) decreases, i.e. there appears the instability by the coupling
constant λ12 . On the other hand, in the opposite phases, the condensation
energy becomes larger, because Δ∗1 Δ2 + Δ1 Δ∗2 < 0. We can expect that
the condensation energy for two-gap superconductivity results in a higher
stability than that of two independent systems, due to the intersublevel
coupling λ12 and the opposite phases.
(d) (e)
n1+2 n1+2 μ
n1+1 n2+1 μ n1+1
n2+1
n1 n1 n2
n2
n1-1 n1-1
n2-1 n2-1
Fig. 4.4 Position of the chemical potential relative the electronic energy levels in a two-
gap superconducting grain. Solid and dotted lines mean two sublevels. (a) Half-filled
system with 2n1 + 2n2 electrons; (b) (2n1 + 1 + 2n2 )-electron system; (c) [2(n1 + 1) + 2n2 ]-
electron system; (d) [2(n1 +1)+2n2 +1]-electron system; (e) [2(n1 +1)+2(n2 +1)]-electron
system.
144 Modern Aspects of Superconductivity: Theory of Superconductivity
4.3.2 Hamiltonian
We consider a Hamiltonian for two bands 1 and 2 written as
H = H1 + H2 + Hint , (4.32)
where
H1 = ε1j a†jσ ajσ − g1 a†j↑ a†j↓ ak↓ ak↑ , (4.33)
jσ jk
H2 = ε2j b†jσ bjσ − g2 b†j↑ b†j↓ bk↓ bk↑ , (4.34)
jσ jk
Hint = g12 a†j↑ a†j↓ bk↓ bk↑ + g12 b†j↑ b†j↓ ak↓ ak↑ . (4.35)
jk jk
The first and second terms of Eq. (4.32) correspond to the reduced BCS
Hamiltonian for bands 1 and 2, respectively. The third term means a cou-
pling between them and corresponds to the pair scattering process between
these two bands (see Fig. 4.5). a†jσ (ajσ ) and b†jσ (bjσ ) are the creation
(annihilation) operator in bands 1 and 2 with spin σ and the single-particle
levels ε1j and ε2j , respectively. The sums of j and k are taken over a set of
Fig. 4.5 Two-band system. The dotted line means the chemical potential; εnj is the
single-particle energy for band n and level j; −g1 and −g2 are the intraband pair inter-
action coupling constants and g12 are the interband pair interaction coupling constant.
Mesoscopic Superconductivity 147
N1 states for band 1 with fixed width 2ω1D and a set of N2 states for band
2 with fixed width 2ω2D , respectively.
In this study, we assume that the Debye energies for two bands coin-
cide, i.e.
ω1D = ω2D = ωD . (4.36)
Within this assumption, N1 and N2 are relatively estimated by the DOS for
two bands as
N1 ρ1
= , (4.37)
N2 ρ2
where ρ1 and ρ2 are the DOS for two bands, respectively. The interaction
constants g1 and g2 can be written as
g1 = d1 λ1 , g2 = d2 λ2 , (4.38)
where d1 and d2 are the mean single-particle level spacings,
2ωD 2ωD
d1 = , d2 = , (4.39)
N1 − 1 N2 − 1
and λ1 and λ2 are the dimensionality interaction parameters for two bands.
We define the interband interaction constant as
g12 = d1 d2 λ12 . (4.40)
In summary, we obtain the relation
ρ1 N1 − 1 d2
≈ = . (4.41)
ρ2 N2 − 1 d1
The system we are considering consists of two half-filled bands, each of
which has equally spaced Nn single-particle levels and Mn (= Nn /2) doubly
occupied pair levels (n = 1, 2). We take the single-particle level spacing as
our energy unity. Thus, the single-particle spectrum is given by
εnj = dn j − ωD , j = 1, 2, . . . , Nn (n = 1, 2). (4.42)
Richardson has obtained his solution within the single-band model for
an arbitrary set of single-particle levels. For simplicity, we assume that there
are no singly occupied single-particle levels. As can be seen from Eqs. (4.33)–
(4.35), these levels are decoupled from the rest of the system. They are said
to be blocked and contribute to the total energy with their single-particle
energies. The above simplification implies that every single-particle level j
is either empty (i.e. |vac) or occupied by a pair of electrons (i.e. a†j↑ a†j↓ |vac
and b†j↑ b†j↓ |vac). These are called the unblocked level.
148 Modern Aspects of Superconductivity: Theory of Superconductivity
d†2
j = 0, [dj , d†k ] = δjk (1 − 2d†j dj ), [d†j dj , d†k ] = δjk d†j , (4.46)
which reflects the Pauli principle for the fermions they were constructed
from.
The Hamiltonian (4.32) for the unblocked levels can then be written as
N1
N1
N2
N2
HU = 2 ε1j c†j cj − g1 c†j ck + 2 ε2j d†j dj − g2 d†j dk
j jk j jk
N1
N2
N2
N1
+ g12 c†j dk + g12 d†j ck . (4.47)
j k j k
N1
N2
HU = 2 ε1j c†j cj − g1 C0† C0 +2 ε2j d†j dj − g2 D0† D0
j j
c†j d†j
[c†j cj , CJ† ] = , [d†j dj , DJ† ] = , (4.53)
2ε1j − E1J 2ε2j − E2J
N1
1 − c†j cj
N2
1 − d†j dj
[C0 , CJ† ] = , [D0 , DJ† ] = , (4.54)
2ε1j − E1J 2ε2j − E2J
j j
N1
1 − c†j cj
[HU , CJ† ] = E1J CJ† + C0† + g1 C0†
2ε1j − E1J
j
N1
1 − c†j cj
+ g12 D0† , (4.55)
2ε1j − E1J
j
N2
1 − d†j dj
[HU , DJ† ] = E2J DJ† + D0† + g1 D0†
2ε2j − E2J
j
N2
1 − d†j dj
+ g12 C0† . (4.56)
2ε2j − E2J
j
⎛ ⎞
†
M2
N2
g12
M2
2g12
+ C0 ⎝ − ⎠|M1 ; M2 (K)U
2ε2j − E2K E2K − E2K
K=1 j K =K
⎛ ⎞
M2
N2
g2
M2
2g2
+ D0† ⎝1 − + ⎠|M1 ; M2 (K)U ,
2ε2j − E2K
E2K − E2K
K=1 j K =K
(4.57)
where
)
L−1 )
M1 )
M2
|M1 (L); M2 U = CJ† CJ† †
DK |vac, (4.58)
J=1 J =L+1 K=1
)
M1 )
L−1 )
M2
|M1 ; M2 (L)U = CJ† †
DK †
DK |vac. (4.59)
J=1 K=1 K =L+1
Comparing Eq. (4.57) with Eq. (4.48), we obtain, for arbitrary J and K,
1 + g1 A1J −g12 A2K |M1 (J); M2 U
C0† †
D0 = 0, (4.60)
−g12 A1J 1 + g2 A2K |M1 ; M2 (K)U
where
Nn
1
Mn
2
AnL = − + . (4.61)
2εnj − EnL
EnL − EnL
j L =L
where we assume that the number of pairs is even. Since the complex
pair energies appear in complex conjugate pairs, the total energy will be
expressed as the real component.
A further transformation is necessary in order to remove the divergences
from the first term of AnL . We define new variables xnλ and ynλ as
where
Then we can rewrite FJK , by using the new variables, and define the result
as Fαβ . We extract the real and imaginary parts of Fαβ as
1 ∗
+
Fαβ = (Fαβ + Fαβ ) = 1 + g1 R1α + g2 R2β + (g1 g2 − g12
2
)R1α R2β
2
− (g1 g2 − g12
2
)I1α I2β cos(φ1α + φ2β )
− {g1 + (g1 g2 − g12
2
)R2β } I1α sin φ1α
− {g2 + (g1 g2 − g12
2
)R1α } I2β sin φ2β , (4.68)
− 1 ∗
Fαβ = (Fαβ − Fαβ ) = −(g1 g2 − g12
2
)I1α I2β sin(φ1α + φ2β )
2i
+ {g1 + (g1 g2 − g12
2
)R2β }I1α cos φ1α
+ {g2 + (g1 g2 − g12
2
)R1α }I2β cos φ2β , (4.69)
152 Modern Aspects of Superconductivity: Theory of Superconductivity
where
Mn
ξnμλ (ξ2nμλ + η2nμ + η2nλ )
+4
(ξ2nμλ + η2nμ + η2nλ )2 − 4η2nμ η2nλ
μ=λ
Nn
2εnλ − ξnλ
− , (4.70)
(2εnλ − ξnλ )2 + η2nλ
j=2λ−1,2λ
1 − ynλ
2
Inλ =
(1 − ynλ )2 Δε2n2λ − (1 + ynλ )2 d2n x2nλ
Mn
ξ2nμλ − η2nμ + η2nλ
− 4ynλ
(ξ2nμλ + η2nμ + η2nλ )2 − 4η2nμ η2nλ
μ=λ
⎫ ,
⎬ 2
Nn
1 Δεn2λ − d2n x2nλ
+ ynλ × , (4.71)
(2εnλ − ξnλ )2 + η2nλ ⎭ ynλ
n=2λ−1,2λ
with
+
Fαβ = 0,
−
Fαβ = 0 (α = 1, 2, . . . , M1 /2; β = 1, 2, . . . , M2 /2). (4.73)
Fig. 4.6 Single-particle levels near the Fermi level in the case of two-band superconduc-
tivity. The dotted lines represent the chemical potential. The left and right bands are
band 1 and band 2, respectively. d1 and d2 are the mean level spacings. (a) 2M1 + 2M2
electron system, where Mn is the number of pair levels; (b) (2M1 + 1) + 2M2 electron
system; (c) (2M1 + 2) + 2M2 electron system.
Numerical calculations are carried out under the condition that N1 :N2 =
3:2, ωD = 50 and λ1 = λ2 = λ.
Fig. 4.7 Typical behavior of pair energy levels of two bands for the ground state. The
parameters used in calculation are N1 = 12, N2 = 8, M1 = 6, M2 = 4, ωD = 50,
0 ≤ λ1 = λ2 ≤ 1.5 and 0 ≤ λ12 ≤ 0.3. The solid and broken lines correspond to the pair
energy levels of band 1 and band 2, respectively.
(a) (b)
Fig. 4.8 Condensation energy. The parameters used in calculation are ω “ D = 50 ”and
λ1 = λ2 = 0.5. Values are normalized by the bulk gap, Δ1 = ωD sinh−1 λ λ λ−λ 2
2 or
“ ” 1 2 12
−1 λ1
Δ2 = ωD sinh λ1 λ2 −λ2
. The solid and broken lines correspond to the condensation
12
energy for band 1 and band 2, respectively. The lines plotted with squares, triangles
and circles are for the 2M1 + 2M2 electron system, (2M1 + 1) + 2M2 electron system
and (2M1 + 2) + 2M2 electron system, respectively. (a) The condensation energy for the
interband coupling parameter λ12 = 0.01; (b) the condensation energy for λ12 = 0.1.
of band 2 for the (2M1 +2)+2M2 electron system is, however, different from
the others. We can also see that the condensation energy is affected by the
interband interaction λ12 . This was mentioned in our previous work.156
Fig. 4.9 Parity gap. The parameters used in calculation are ωD = 50 and λ1 = λ2 = 0.5.
The solid and broken lines correspond to the parity gap for band 1 and band 2, respectively.
The lines plotted with triangles and squares are for the interband coupling parameter
λ12 = 0.01 and λ12 = 0.1, respectively. Values are normalized by the bulk gap.
point is, however, much less than that for the case of calculations for the sin-
gle band. The parity gap is almost independent of the interband interaction
λ12 . This was also mentioned in our previous work.156
We have extended Richardson’s exact solution to the two-band system,
and have derived a new coupled equation. To investigate the properties of
the two-band superconductivity, we have solved the equation numerically
and determined the behavior of pair energy levels, the condensation energy
and the parity gap.
The band whose mean level spacing is larger than that of the other band
degenerates and condenses faster. The behavior of the condensing band is
qualitatively the same as that in the case of calculations for the single band.
The coexistence of the normal band and the condensed one may be reflected
in the opposite phases of the gaps of these bands. This phase character
appears in all the results of numerical calculations. Therefore, the phase of
a gap is important for stabilizing the two-band superconductivity.
We have also calculated the condensation energy and the parity gap for
two-band superconductivity. The results suggest that the interband inter-
action λ12 affects the condensation energy, but not the parity gap.
In summary, the expression of Richardson’s exact solution for two-band
superconductivity has been presented, by solving numerically a new cou-
pled equation. Then, the behaviors of pair energy levels, the condensation
Mesoscopic Superconductivity 157
energy and the parity gap have been determined. The results for the con-
densed band are almost qualitatively the same as those for the single-band
calculation, and the coexistence of the normal band and the condensed one
may have originated from the opposite phases of the gaps of these bands.
moments from the impurities tend to align the spins of the electron pairs
in the superconductor, which often results in a strongly reduced transition
temperature. Buitelaar et al. have experimentally investigated the Kondo
effect in a carbon nanotube quantum dot coupled to superconducting Au/Al
leads.183 They found that the superconductivity of the leads does not destroy
the Kondo correlations in the quantum dot at the Kondo temperature. A
more subtle interplay has been proposed for exotic and not well-understood
materials such as heavy-fermion superconductors, in which the two effects
might actually coexist.189
In this paper, we investigate the Kondo effect and the superconductivity
in ultrasmall grains by using a model which consists of the sd and reduced
BCS Hamiltonians with the introduction of a pseudofermion. The mean-
field approximation for the model is introduced, and we calculate physical
properties of the critical level spacing and the condensation energy. These
physical properties are discussed in relation to the coexistence of the super-
conductivity and the Kondo regime. Finally, we derive the exact equation
for the Kondo regime in a nanosystem and discuss the condensation energy
from the viewpoint of the correlation energy.
4.4.2 Model
We consider a model coupled to the superconductivity for quantum dots
to investigate the Kondo effect in normal metals, which can be expressed
by the effective low-energy Hamiltonian obtained by the Schrieffer–Wolff
Mesoscopic Superconductivity 159
transformation190 :
H = H0 + H1 + H2 , (4.77)
where
H0 = εk a†kσ akσ + Eσ d†σ dσ , (4.78)
k,σ σ
' (
H1 = J S+ a†k ↓ ak↑ + S− a†k ↑ ak↓ + Sz a†k ↑ ak↑ − a†k ↓ ak↓ , (4.79)
k,k
H2 = −g a†k↑ a†k↓ ak ↓ ak ↑ . (4.80)
k,k
a†kσ (akσ ) and d†σ (dσ ) are the creation (annihilation) operator corresponding
to conduction electrons and the effective magnetic particle as an impurity,
respectively. In this study, we assume that the magnetic particle is a fermion
with S = 1/2, for simplicity. E represents an extraction energy given by
E↑,↓ = −E0 ±Ez , including the Zeeman effect. The second term of Eq. (4.77)
is the interaction between conduction electrons and the spin in a quantum
dot. S is the spin operator as S+ = d†↑ d↓ , S− = d†↓ d↑ and Sz = (d†↑ d↑ −
d†↓ d↓ )/2. The third term corresponds to the interaction between conduction
electrons included in the pairing-force Hamiltonian.
Here, we introduce a pseudofermion for the magnetic particle opera-
tor 191 as
d†↑ = f↓ , d↑ = f↓† ,
and we have |σ = fσ† |0. The spin operator S can be presented as S+ = f↑† f↓ ,
S− = f↓† f↑ and Sz = (f↑† f↑ − f↓† f↓ )/2. The Hamiltonian can be rewritten as
†
H0 = ε̃k ckσ ckσ + Efσ† fσ , (4.83)
k,σ σ
H1 = J fσ† fσ c†k σ ckσ , (4.84)
k,k ,σ,σ
160 Modern Aspects of Superconductivity: Theory of Superconductivity
H2 = −g c†k↑ c†k↓ ck ↓ ck ↑ , (4.85)
k,k
. . †
where ckσ = i Uik aiσ and ε̃k = i,j Uki [εi δij − J/2]Ujk . For the sake of
simplicity, we focus only on Ez = 0 without an external magnetic field:
E = E0 .
Using these order parameters in Eqs. (4.84) and (4.85), we obtain the mean-
field Hamiltonian
† √
HMF = ε̃k ckσ ckσ + Ẽfσ†fσ + 2J [Ξfσ c†kσ + Ξ∗ ckσ fσ† ]
k,σ σ k,σ
−g [Δ∗ ck↓ ck↑ + Δc†k↑ c†k↓ ]. (4.88)
k
The constraint (4.82) is taken into account by the second term with a
Lagrange multiplier λ. In this study, we assume a constant DOS with
the energy region of the Debye energy, and the coupling constants can be
expressed as J = dJ˜ and g = dλ.
Fig. 4.10 Gap function and spin-singlet order. (a) The gap function. This vanishes in the
region of λ values less than λc . (b) Spin singlet order parameter. In the case of J˜ < J˜c ,
the singlet order vanishes. The system consists of eight energy levels and eight electrons
with the level spacing d = 1.0 and ωD = 1.0.
162 Modern Aspects of Superconductivity: Theory of Superconductivity
expressed as
Ξ(ξk − x)
Ξ= , (4.89)
(ξk − x)2 + Ξ2
k
' + (
where ξk = ε̃k −μ, x = ε̃k + Ẽ ± (ε̃k − Ẽ)2 + 4Ξ2 /2 and μ is the chemical
potential. For the case of the critical level spacing, the solution shows that
the spin singlet order parameter vanishes. From Eq. (4.89), we can find the
critical level spacing dKondo
c for the Kondo regime.
γ 1
dcKondo
= 4ωD e exp − √ . (4.90)
2 2J˜
√
When the coupling parameter J˜ is smaller than J˜c = [2 2(ln(4ωD /d) +
γ)]−1 , the spin singlet order parameter vanishes.
Figure 4.10(b) presents the spin singlet order parameter given by
Eq. (4.86) in the case g = 0. In the region of J˜ < J˜c , the order param-
eter vanishes. This result suggests the critical level spacing in the Kondo
effect.
Fig. 4.11 Physical properties in a coupled system. (a) Gap function and spin singlet order
parameter. (b) Condensation energy. J/˜ J˜c = 0, 0.94741, 1.8948, 2.8422. Other parameters
are the same.
Figure 4.11(b) shows the condensation energy for several λ and J˜ values.
We have found that the condensation energy of the coupled system between
the superconductivity and the Kondo regime becomes lower than that for the
pure superconductivity. In the coexistence region, the highest value of the
condensation energy appears in all cases.
N
λ
n
2λ
2− + = 0, (4.91)
ε̃k − Ei El − Ei
k=1 l=1, l=i
164 Modern Aspects of Superconductivity: Theory of Superconductivity
Fig. 4.12 Exact solution for the superconductivity. (a) Condensation energy of the exact
solution and that obtained in the mean-field approximation. (b) Pairing energy level with
energy level obtained in the mean-field approximation. Eight energy levels, eight electrons,
d = 1.0, ωD = 4.0.
where N and n are the number of orbitals and the number of the occupied
orbitals, respectively, and Ei corresponds to the exact orbital. Figure 4.12
shows the condensation energy and the pairing energy level for the nanosize
superconductivity. Note that the physical properties obtained in the mean-
field approximation give a good description of the high DOS (d → ∞).
We have found the different behavior of the condensation energy from
that obtained in the mean-field approximation, as shown in Fig. 4.12(a).
Figure 4.12(b) presents the qualitative behavior of the pairing energy level in
the ground state. At λ about 1.6, more than two energy levels in Fig. 4.12(b)
are completely paired. The pairing behavior has already been reported by
many groups.138, 172
Let us derive the exact equation for the Kondo regime in ultrasmall
grains. We can consider the Hamiltonian H = H0 + H1 in Eq. (4.77). We
introduce a creation operator describing all excited states at the spin singlet
coupling between a conduction electron and a pseudofermion:
†
c† fσ
kσ
Bj = , (4.92)
ε̃k − Ej
k,σ
where Ej represents the exact eigenenergies in the Kondo regime. The exact
eigenstate |Ψn for the Kondo regime can be written as |Ψn = Πnυ=1 Bυ† |0.
Other electrons, which are not related to the spin singlet order, contribute
.
Esingle = nk=1 ε̃k to the eigenenergy. The ground-state energy EGS can be
.
written as EGS = nk=1 [Ek + ε̃k ].
Mesoscopic Superconductivity 165
Fig. 4.13 Condensation energy for the Kondo regime: all parameters used in the system
are as follows: eight energy levels, eight electrons, d = 1.0 and ωD = 4.0.
166 Modern Aspects of Superconductivity: Theory of Superconductivity
approximation. From the results, we have found that strong local magnetic
moments from the impurities reduce the transition temperature for super-
conductivity. However, weak couplings λ of the superconductivity do not
destroy the spin singlet order parameter at all. These results are in good
agreement with the experimental results.183 We have found that there is a
coexistence region for the superconductivity and the Kondo regime.
Finally, we have derived the exact equation for the Kondo regime in
a nanosystem, which was not an easy task, and have discussed the con-
densation energy from the viewpoint of energy levels. Further study of the
properties in the Kondo regime with the use of the exact equation will be
presented elsewhere.
In summary, we have investigated the Kondo effect and the supercon-
ductivity in ultrasmall grains by using a model which involves the sd and
reduced BCS Hamiltonians with the introduction of a pseudofermion. The
mean-field approximation for the model has been introduced, and we have
calculated physical properties of the critical level spacing and the conden-
sation energy. These physical properties have been discussed in relation to
the coexistence of the superconductivity and the Kondo regime. Finally, we
have derived the exact equation for the Kondo regime in a nanosystem and
discussed the condensation energy from the viewpoint of energy levels.
Let us introduce a coordinate system with the origin at the first magnetic
moment m1 and the polar axis along the line connecting magnetic moments.
In this coordinate system, the magnetic moments have the components mi =
mi (cos ϕi sin θi , sin ϕi sin θi , cos θi ), and their interaction energy (4.110) is
written in the form
μ0 exp(−R12 /λ)
U= m1 m2 · 3 · f (θi , ϕi , R12 /λ), (4.111)
4π R12
where
and ϕ = ϕ2 − ϕ1 .
Differentiating the function f (θi , ϕi , x) with respect to angular variables
and equating the results to zero, we find that there are four possible stable
172 Modern Aspects of Superconductivity: Theory of Superconductivity
energy configurations:
θ1 = θ2 = 0,
θ1 = θ2 = π, (4.113)
0 ≤ ϕ < 2π,
π
θ1 = θ2 = ,
2 (4.114)
ϕ = π,
π
θ1 = θ2 = ,
2 (4.115)
ϕ = 0,
θ1 = 0, θ2 = π,
(4.116)
0 ≤ ϕ < 2π,
which are illustrated in Fig. 4.15.
Fig. 4.15 Four energy saddle points of a pair of ferromagnetic quantum dots, as defined
in Eq. (4.113). Of the four, 1 and 2 correspond to the ground state, depending on the
condition (4.118). States 3 and 4 are never stable.
Mesoscopic Superconductivity 173
Further analysis shows that configurations (4.115) and (4.116) are saddle
points, not energy minima. Evaluating the second derivatives of (4.112), we
obtain the stability condition of the configuration (4.113):
2
R R
− − 1 ≤ 0. (4.117)
λ λ
This implies that the ferromagnetic ordering of the pair of magnetic
moments is possible if
R 1 √
≤ (1 + 5). (4.118)
λ 2
It turns out that if the condition (4.118) is violated, then the alternative
configuration (4.114) is a stable energy minimum. We can conclude that if
the temperature changes, and the penetration depth parameter λ(T ) varies
in such a way that the condition (4.118) is not satisfied, then the ground-
state orientation will change from (4.113) to (4.114).
This result is readily generalized to ordered arrays of ferromagnetic gran-
ules, and so we conclude that orientational phase transitions are possible in
systems of quantum dots in a superconducting matrix. For example, it is
clear that the condition (4.118) can be applied to linear chains of quantum
dots. On the other hand, the results for square or cubic lattices remain to
be determined.
We have considered the interactions between nanoscale magnetic dots
embedded in a bulk superconducting material. Our approach is valid for
materials which are well described in the London limit r0 λ, since RKKY
interactions are negligible. We have shown that, depending on the dimen-
sionless parameter R/λ, different stable ground states occur. So, as the tem-
perature varies, orientational phase transitions will take place for periodic
arrays of such quantum dots. Of course, our calculation does not include all
types of interactions which define the orientation of magnetic moments in
space. In particular, we neglect the energy of a magnetic anisotropy of gran-
ules which is determined by the shape of granules or the type of their crystal
lattice. However, when the shape of granules is close to the spherical one
and the lattice of a ferromagnet has the cubic symmetry, then Eq. (4.109)
will essentially be exact.
In the experimental systems studied by V. V. Moshchalkov,201 a square
of Pt/Co magnetic nanodots was deposited on the surface on the Pb super-
conductor, which is of type I(k = 0.48). The dots were about 0.26 μm in
diameter, and they were deposited on the grid with a spacing of 0.6 μm.
For Pb, the penetration depth is 39 nm at low temperatures. So, this array
174 Modern Aspects of Superconductivity: Theory of Superconductivity
was in the limit R > λ and the dot–dot interaction would be expected to
correspond to the antiferromagnetic alignment shown in the second state
in Fig. 4.15. With increase in the temperature, the transition to the ferro-
magnetic alignment would occur according to Eq. (4.118) at λ = 0.36R, i.e.
219 nm.
−1/2
According to the Casimir formula λ(T ) = λ(0)(1 − t4 ) , with t =
T /Tc , this would occur at T = 7.14 K, compared with Tc = 7.2 K. Therefore,
the experimental conditions for the transition to be observed are certainly
feasible. Of course, for an exact comparison with theory in this case, our
theory should be generalized to deal with magnetic particles near the surface
rather than with those embedded in the bulk of a superconductor.
Of course, it would be interesting in the future to generalize our results
to superconductors in the Pippard limit, where the RKKY interactions
between quantum dots will dominate over dipolar forces.203
The expression (4.110) can be used to study, by means of numerical
methods, the magnet configurations and the orientational phase transitions
in an ensemble of nanogranules. It is possible to determine the conditions
of orientational transformations in the analytic form for the ordered struc-
tures (a chain of granules, plane and volume lattices). Inasmuch as the state
of the magnetic subsystem of a specimen at phase transitions is changed,
this phenomenon can be experimentally observed under the change in the
magnetic susceptibility in the region of low fields.
Fig. 4.16 Two-dimensional lattice of magnetic points. Magnetic moments (black) are
directed up or down (white); a is the lattice constant.
Fig. 4.17 Two-dimensional lattice of magnetic points. Magnetic moments are aligned
in the plane in the form of chains; a is the lattice constant.
Fig. 4.18 Part of the lattice with a modulated planar distribution of the magnetization; a
is the lattice constant. Circles represent magnetic points, and the arrows on them indicate
the directions of magnetic moments mn k in the base plane (n and k are the spatial indices
of magnetic points). The angle ϕ defines the deviation of the moments of magnetic points
from the principal direction of the lattice. The states of all points in the given configuration
are equivalent. The net magnetic moment is zero. The separated circles schematically
denote magnetic vortices. The dotted lines are tangents to the directions of magnetic
moments at sites of the lattice.
simpler:
N α · r (r · m )(r · m )
μ0 N ∂ ∂2 i i i i 0,0
U =− 3 1− + 2 exp −
4π 2 ∂α ∂α δ 5
ri
i
α · r m · m
∂2
N
N ∂ i i 0,0
+ 1− + 2 exp − . (4.119)
2 ∂α ∂α δ ri3
i
By writing the formula (4.119), we used the method of differentiation with
respect to the parameter α, which should be set equal to 1 after the calcu-
lations. The index i stands for the summation over all sites of the lattice,
δ — the penetration depth for 2D case. Performing the summation in the
relation (4.119), it is convenient to introduce the pair of indices (n, k) defin-
ing the position of a site in the lattice (Fig. 4.18) instead of the running
index of magnetic points i. It is easy to see that the system represented in
Fig. 4.18 possesses the translational invariance with a period of 2a so that
mn,k = mn+2l, k+2p ,
l, p = ±1, ±2, . . . . (4.120)
178 Modern Aspects of Superconductivity: Theory of Superconductivity
The lattice has only four types of magnetic points differ which from
one another by a spatial orientation of magnetic moments. Their vector
components depend on the angle ϕ in the following manner:
⎛ ⎞
cos ϕ
mn,k = mn+2l,k+2p , m2l,2p = m0,0 = m⎝sin ϕ ⎠,
0
⎛ ⎞
cos ϕ
m2l+1,2p = m1,0 = m⎝− sin ϕ⎠,
0
⎛ ⎞
− cos ϕ
m2l,2p+1 = m0,1 = m⎝ sin ϕ ⎠,
0
⎛ ⎞
− cos ϕ
m2l+1,2p+1 = m1,1 = m ⎝ − sin ϕ ⎠,
0
l, p = ±1, ±2, . . . ,
l, p = ±1, ±2, . . . , (4.121)
where m is the modulus of the magnetic moment of a site.
After the substitution of Eq. (4.121) in Eq. (4.119) and the summation
over sites of the unbounded lattice, we get the following interesting result.
It turned out that the interaction energy of the system of magnetic points
(see Fig. 4.18) does not depend on the angle ϕ and is determined by the
relation
U μ0 m2
= · F (a/δ), (4.122)
N 4π a3
where N is the number of sites of the lattice, m is the magnetic moment of
a granule and F (a/δ) is the energy characteristic of a magnetic state which
is a universal function of a single parameter and determines the dependence
of the energy on both the period and the field penetration depth in the case
where the magnetic moments are distributed in the base plane of the lattice.
It can be represented in the form of a sum:
∞ ∞
μ0 ∂ ∂2 1 exp(−2α(a/δ) l2 + p2 )
F(a/δ) = −1 + + ·
4π ∂α ∂α2 8 (l2 + p2 )3/2
l=1 p=0
∞ ∞
∂ ∂2 1
− −1 + + 2 ·
∂α ∂α 4
l=−∞ p=−∞
Mesoscopic Superconductivity 179
(4.125)
The calculation of the functions F (a/δ) and Φ(a/δ) on the basis of the
relations (4.123) and (4.125) is not a difficult task and can be realized with
180 Modern Aspects of Superconductivity: Theory of Superconductivity
Fig. 4.19 Plots of the energy characteristics F (a/δ) and Φ(a/δ) of states of the lattice
with normal and planar orientations of magnetic moments, respectively. Values of F (0)
and Φ(0) correspond to the transition of the superconducting matrix to the normal state.
|f = a|0 + b|1,
182 Modern Aspects of Superconductivity: Theory of Superconductivity
where
a2 + b2 = 1. (4.126)
Since the devices can register only classical quantities, the measurement of
a quantum bit will give the states |0 or |1 with probabilities determined
by the squares of the coefficients a and b.
The use of the principle of superposition allows one to increase the infor-
mational space exponentially with linear growth of the size of a physical
system, because the register including n qubits can be in a superposition of
at once 2n states. In addition, quantum mechanics admits the existence of
the so-called entangled states possessing no analogs in classical physics. An
ensemble of qubits is a collection of qubits in different but given states.
It is worth noting that, even on the level of mathematics, a quantum
computer operates in a basically different way than a classical computer.
Input data are coded in “quantum cells of memory.” In this case, the col-
lection of qubits becomes a single quantum system. This system undergoes
the sequence of elementary quantum operations. Quantum computations
are a realization of the most astonishing idea of applying the principles of
quantum mechanics to the world of computers. Ordinary computers, despite
their complexity, use classical laws of mechanics. In recent years, the the-
ory of classical computations has been developed on the basis of works
by A. Turing. With the appearance of quantum computations, new possi-
bilities have arisen, and the situation has been radically changed. Quantum
methods can be successfully used in the solution of mathematical prob-
lems, though the time consumed for the solution of mathematical problems
increases exponentially with the complexity of the problem.
For the theory of quantum calculations, the physical nature of qubits is
not of crucial importance; the basically important point is that the system
in the course of calculations obeys the laws of quantum mechanics.
(atoms and electrons) but also for large electronic schemes, which include
superconducting bits (qubits).
Many years ago, attempts to construct a Josephson computer on the
basis of the Josephson tunnel logic failed hopelessly. The principal rea-
son for the failure was the huge technological dispersion of parameters of
tunnel junctions, which did not allow one to produce large microcircuits.
A Josephson computer can be created only on the way, having nothing in
common with that based on semiconductors; namely, it can be just a quan-
tum computer.
The hope is related to two circumstances. First, the fabrication of super-
conducting qubits is quite possible in the framework of up-to-date tech-
nology. Second, the presence of a gap in the spectrum of excitations of
a superconductor allows one to expect the suppression of generation in a
system.
Let us consider the first steps on the way toward construction of a super-
conducting quantum computer. First, we consider, in brief, the Josephson
effect.
superconducting current
Jc = J0 sin ϕ, (4.127)
dϕ 2eV
= , (4.128)
dt
where ϕ is the phase difference on the interface of superconductors, V the
applied voltage, and J0 the critical current through the junction.
The Josephson effect indicates the existence of the electron ordering in
superconductors, namely the phase coherence: in the ground state, all elec-
tron pairs have the same phase ϕ which characterizes their wave function
√
Ψ1 = ns1 eiϕ1 . According to quantum mechanics, the presence of a phase
difference must cause a current through the junction. The discovery of such
a current in experiment proves the existence of macroscopic phenomena in
nature which are directly determined by the phase of a wave function:
√
Ψ = neiψ . (4.129)
The Josephson current in the scheme drawn in Fig. 4.21 can be easily deter-
mined with the help of the elementary Feynman approach to a Josephson
junction as a two-level quantum-mechanical system.4 By introducing an
intermediate object, some island, with the wave function Ψ0 into the two-
level system characterized by the wave functions Ψ1 = ns1 /2eiϕ1 and
Ψ2 = ns2 /2eiϕ2 , we write the Schrödinger equation in the form
dΨ1 eV
i = Ψ1 + KΨ0 , (4.130)
dt 2
dΨ0
i = KΨ1 + KΨ2 + E0 Ψ0 , (4.131)
dt
dΨ2 eV
i = KΨ0 − Ψ2 . (4.132)
dt 2
The formula for a constant current running through two series-connected
Josephson junctions at the zero external potential difference takes the form
√
∂ nS1 K2 √
J = =− nS2 sin(ϕ2 − ϕ1 ). (4.133)
∂t E0
The Josephson effect is still one of the phenomena that make superconduc-
tors such a fascinating area of study. Despite more than 40 years of intensive
studies and numerous applications, it remains an important field of research
in connection with the use of small superconducting grains.
4.7.3.2 SQUIDs
where
Δϕ1 = ϕ2A − ϕ1A ,
Δϕ2 = ϕ2B − ϕ1B
are the phase differences on the first and second Josephson junctions. There
occurs a distinctive interference of the superconducting currents running
through these connections.
Inside the superconductor, the current is zero: j ≡ 0. We will use the
following formula for the current j:
2e
j ∼ h∇ϕ − A. (4.137)
c
We can write
2e
ϕ1B − ϕ1A = A · dl, (4.138)
hc C1
2e
ϕ2B − ϕ2A = A · dl. (4.139)
hc C2
Summing up these two equations, we get 7
2e Φ
ϕ1B − ϕ2B + ϕ2A − ϕ1A = A · dl = 2π . (4.140)
hc Φ0
Thus, we have
Φ
Δϕ1 − Δϕ2 = 2π , (4.141)
Φ0
where Φ is the total quantum flux.
The flux quantum is defined as
h
Φ0 = . (4.142)
2e
For a balance SQUID ring system, we can write
Φ
Δϕ1 = ϕ0 + π ,
Φ0
Φ
Δϕ2 = ϕ0 − π . (4.143)
Φ0
The total current in the SQUID is
I = I0 sin(Δϕ1 ) + I0 sin(Δϕ2 )
Φ Φ
= I0 sin ϕ0 + π + I0 sin ϕ0 − π
Φ0 Φ0
Φ Φ
= 2I0 sin(ϕ0 ) cos π
= Imax cos π , (4.144)
Φ0 Φ0
with Imax = 2I0 sin(ϕ0 ).
In this case, the critical current turns out to be periodically dependent
on the flow of an external magnetic field, which allows one to use such a
unit for exact measurement of the magnetic field.
Mesoscopic Superconductivity 187
Let us consider the first steps on the way towards creation of a super-
conducting computer. The simplest superconducting system demonstrating
the coherence is the SQUID, which is a superconducting ring including a
Josephson junction at one point. The energy of this system contains two
terms, the Josephson transition energy (cos Φ) and the energy related to
the ring L:
Φ (Φ − Φx )2
H = −EJ cos 2π + . (4.145)
Φ0 2L
Here, Φ is the difference of superconductive phases at the junction. The
superconductive phase in the ring is proportional to a magnetic flow applied
to the ring (quantization of the magnetic flow). If Φx is equal exactly to a
half of the magnetic flow, the potential of a SQUID becomes doubly degen-
erate.
Two minima of the well correspond to the currents in the ring passing
in the clockwise and counterclockwise directions, respectively.208
A superposition of these states in SQUIDs, was observed experimentally
in Refs. 208, 211. These experiments demonstrated clearly the possibility of
creating a superposition of states in a system with a macroscopic number of
particles. In the given case, the circular current including 1013 electrons was
registered in a loop. The states participating in a superposition were macro-
scopically distinguishable, by differing from one another by the currents,
whose difference was several microamperes. Recently, an important notion
has been introduced in the course of studies of structures with Josephson
junctions. It is the notion of macroscopic quantum coherence. In such sys-
tems, the Josephson energy can have two almost-degenerate minima at val-
ues of the phase which are separated by a potential barrier (see Fig. 4.23).
It is possible that the phase passes from one minimum to another due to
the quantum-mechanical tunneling, and the eigenstates of the system are
superpositions of the states localized in the first and second minima.
Fig. 4.24 Cooper-pair box [× denotes a Josephson junction (JJ)]. Cg and Vg are the
capacitor and the controlling voltage.
Mesoscopic Superconductivity 189
There exists one more possibility of realizing SQUIDs with the use of HTSCs
possessing lattice, d pairing (π-loop ones).210
The physics of d pairing was considered in Chap. 2. Quantum processors
on the base of these SQUIDs are being developed at the Canadian company
D-wave.
We note also that, in addition to superconducting qubits, quantum com-
puters use qubits possessing other physical properties. Scientists at Yale
University used a very fine aluminum plate in the fabrication of a quantum
chip. A single qubit consists of one billion aluminum atoms, which never-
theless behaves as a single unit that can be in two energy states, denoted as
0 and 1. Such quantum-mechanical states of a qubit cannot be long-term —
their lifetime is about one microsecond. But it is sufficient for a chip to
solve the so-called algorithm. We have considered the technologies of super-
conducting computers which represent a new type of quantum computers.
These computers are based on the other mechanism obeying the laws of
quantum mechanics. We recall that, till recently, the principle of devices
was invariable, and the archetype of such devices is mechanical clocks. In
such devices, all stages of their relative motion can be observed; therefore, it
is quite simple to understand their structure. On the contrary, the principle
190 Modern Aspects of Superconductivity: Theory of Superconductivity
In this book we have concentrated our attention on the most urgent prob-
lems of superconductivity, such as the nature of high-temperature supercon-
ductivity, mechanisms and symmetry of pairing, two-gap superconductivity
in magnesium diborades, mesoscopic superconductivity and the problems of
room-temperature superconductivity.
We considered some questions concerning the application of quantum
field theory to the problems of superconductivity. Quantum field theory pro-
vides the original and powerful means for the solution of certain problems of
superconductivity. In the field of superconductivity, we meet the problem-
maximum — it consists the creation of room-temperature superconductors.
We considered this problem in our book, gave some recommendations on the
search for these superconductors and analyzed the possibility of the fabri-
cation of artificial materials possessing the property of superconductivity at
room temperature.
We also touched on the questions of the application of superconducting
qubits to quantum computers. It was shown that the description of super-
conducting qubits is based on the laws of quantum mechanics.
It is worth noting that HTSCs have been investigated for more than
two decades with great efforts but the whole pattern of the phenomenon
is not yet available. We are sure that the comprehension of HTSCs will be
attained when our knowledge of HTSCs reaches the critical level which will
be sufficient for understanding the huge amount of experimental data from
a single viewpoint. The study is in progress. In particular, we mention the
191
192 Modern Aspects of Superconductivity: Theory of Superconductivity
1
υ
+√ χsq a+υ s
k,σ ak−q,σ Qq + Ωqs b+s s
q bq , (A.1)
N k,σ,υ,q,s q,s
where a+vkσ are creation operators for the electron with the momentum k in
band υ, avkσ are annihilation operators for the electron with the momentum
k in band υ, b+s q and bsq are creation and annihilation operators for the
phonons with the momentum q, εvk is the energy of the electron with the
impulse k, Ωqs is the energy of the phonons, Qsq = b+s s
q + bq , s is the number
of phonon branches, V −q = Vq , Xq∗s = Xqs are the Fourier components
of the Coulomb interactions of the electron and their coupling constant
with the lattice phonons, respectively. We can see that both constants are
independent of the spin index of an electron. Let us introduce new operators
for electron and phonon systems by the rules
193
194 Modern Aspects of Superconductivity: Theory of Superconductivity
where we have united two indices and so υ = (υ, σ) and μ = (μ, σ ) are the
complex indices which characterize the number of crystals in the electron
zone and the spin of the electron. The Hamiltonian describing the system
of interacting electrons and phonons of the crystal after transformation by
a unitary operator is written in the form98–102
2
1 χsq
H= εk −
v
A+k,v Ak,v
N s,q Ωs,q
k,v
1 χsq 2
+ Vq − 2 A+ +
k,v Ak ,v Ak −q,v Ak−q,v
2N Ωs,q
q,k,v,k ,v s
+ Ωs,q β+
s,q βs,q + higher-order term. (A.4)
s,q
The equation for the Green function can be found from the equation of
motion:
3 4
∂8
9
∂
i −iT Ak,v (t)Ak,v (t ) = δ(t − t ) + −iT
+ +
Ak,v (t) Ak,v (t ) ,
∂t ∂t
(A.6)
Appendix A: Two-particle Green Function for Multiband Superconductors 195
where
∂Ak,v 8 9 8 9
i = [Ak,v · H] = −iT Ak,v (t)A+ k,v (t ) = −iT Ak,v (t)A+
k,v (t )
∂t
1
= ε̃vk Ak,v + Ṽq A+
k2 ,v2 Ak2 +q,v2 Ak−q,v
N
q,k2 ,v2
1
= ε̃k,v Ak,v − Ṽq Ak+q,v Ak2 −q,v2 A+
k2 ,v2 , (A.7)
N
q,k2 ,v2
with
1 χsq 2
ε̃vk = εvk − , (A.8)
N s,q Ωs,q
1
ε̃k,v = ε̃k +
v
Ṽq=0 − n Ṽq . (A.9)
v
N q
2
Therefore, the equation for the Green function (A.6) is changed, after we
insert (A.7):
∂ 8
9
i −iT Ak,v (t)A+ k,v (t )
∂t
8 9
= δ(t − t ) + ε̃k,v −iT Ak,v (t)A+
k,v (t
)
1 8 9
− Ṽq −iT Ak+q,v (t)Ak2 −q,v2 (t)A+
k2 ,v2 (t − 0)A+
k,v (t ) .
N
q,k2 ,v2
(A.10)
Such a two-particle Green function satisfies an equation of the Bethe–
Salpeter type (we do not split this function into two one-particle Green
functions of the Gorkov type). The solution of this equation according to
the Bogoliubov–Tyablikov method gives rise to the following expression for
the Fourier component of the two-particle Green function:
.
k2 ,υ;k1 ,μ f (k1 , μ; k2 , υ; ω) σ,σ ϕ(μ, υ; σ, σ )
G2 ω ∼ , (A.11)
k1 + q, μ; k2 − q, υ 1 − V K(k1 , μ; k2 , υ; ω)
with
μ υ
1 1 − nk1 +q − nk2 +q
K(k1 , μ; k2 , υ; ω) = , (A.12)
N q ω − εk1 +q,μ − εk2 +q,υ
196 Modern Aspects of Superconductivity: Theory of Superconductivity
Fig. A.1 The effect of temperature on the K(k, μ; k, υ; ω) function (E = k2 /2m = 0) for
different structures of the energy zones: (a) Δ1 = −1, m∗1 = 1; (b) Δ2 = 0.2, m∗2 = −2;
(c) Δ1 = −1, m∗1 = 1, Δ2 = 0.2, m∗2 = −2; curves 4, T = 2 K; curves 3, T = 10 K; curves
2, T = 50 K; curves 1, T = 100 K; arrows 5–7 correspond to different 1/V values (arrow
5, V = −0.07; arrow 6, V = −0.06; arrow 7, V = −0.05). All energy values are taken as
arbitrary values, i.e. V = V /M , where M is a scale factor (M = 1 eV, for convenience;
m∗i = mm
i
, where m is the free electron mass).
Fig. A.2 The temperature dependence of the superconducting gap. Dashed line —
theory; continuous line — BCS; a axis, b axis — experiment.
203
204 Modern Aspects of Superconductivity: Theory of Superconductivity
21. Kimura, M., Kawabe, H., Nakajima, A., Nishikawa, K. and Aono, S. (1988).
Possibility of superconducting and other phases in organic high polymers of poly-
acene carbon skeletons, Bull. Chem. Soc. Jpn. 61, pp. 4239–4224.
22. Kubo, R. and Tomita, J. (1954). A general theory of magnetic resonance absorption,
Phys. Soc. Jpn. 9, pp. 888–919.
23. Doi, Y. and Aono, S. (1974). Theory of line shape in nuclear magnetic resonance,
Prog. Theor. Phys. 51, pp. 1019–1029.
24. Hebel, L. C. and Slichter, C. P. (1959). Nuclear spin relaxation in normal and super-
conducting aluminum, Phys. Rev. 113, pp. 1504–1519.
25. Fibich, M. (1965). Phonon effects on nuclear spin relaxation in superconductors,
Phys. Rev. Lett. 14, pp. 561–564.
26. Maki, K. (1967). Gapless Superconductors in Superconductivity, ed. Parks, R. D.
(Marcel Dekker, New York).
27. Uchino, K., Maeda, A. and Terasaki, I. (1959). Nature of the High Temperature
Superconductor (in Japanese) (Baifuukann).
28. Aono, S. (1998). Superconductivity as symmetry breaking classical treatment, Bull.
Chem. Soc. Jpn. 71, pp. 49–56.
29. Scalapino, D. J. (1995). The case for dx2 −y 2 pairing in the cuprate superconductors,
Phys. Rep. 250, pp. 329–360.
30. Gorkov, L. P. (1959). Microscopic derivation of the Ginzburg–Landau equations in
the theory of superconductivity, Sov. Phys. JETP 9, pp. 1364–1375.
31. Negele, J. W. and Orland, H. (1987). Quantum Many-Particle Systems (Addison-
Wesley, Reading, Masscchusetts).
32. Kawabe, H., Nishikawa, K. and Aono, S. (1994). Path integral approach to correla-
tion energy, Int. J. Quantum Chem. 51, pp. 265–283.
33. Feynman, R. P. (1948). Space–time approach to non-relativistic quantum mechanics,
Rev. Mod. Phys. 20, pp. 367–387.
34. Hubbard, J. (1959). Calculation of partition functions, Phys. Rev. Lett. 3, pp. 77–78.
35. Das, A. (2006). Field Theory: A Path Integral Approach (World Scientific,
Singapore).
36. DeWitt, B. (1984) Supermanifolds (Cambridge University Press).
37. Ferraz, A. and Nguyen, A. V. (1995). Supersymmetry and electron–hole excitations
in semiconductors, Phys. Rev. B 51, pp. 10548–10555.
38. Yang, C. N. (1962). Concept of off-diagonal long-range order and the quantum phases
of liquid He and of superconductors, Rev. Mod. Phys. 34, pp. 694–704.
39. Dunne, L. J., Murrell, J. N. and Brändas, E. J. (1997). Off-diagonal long-range
order, eta-pairing, and Friedel oscillations in high Tc cuprate superconductors and
the ground state of the extended Hubbard model, Int. J. Quantum Chem. 63,
pp. 675–684.
40. Ginzburg, V. L. and Kirzhnitz, D. A. (1982) (eds). High-Temperature Super-
conductivity (Consultants Bureau, New York).
41. Fröhlich H. (1954). On the theory of superconductivity, Proc. R. Soc. London A
223, p. 296–305.
42. Bogolyubov, N. N. (1958). A new method in the theory of superconductivity, Zh.
Eksp. Teor. Fiz. 34, pp. 58–65.
43. Josephson, B. D. (1974). The discovery of tunneling superconductivity, Rev. Mod.
Phys. 1074, pp. 251–254.
44. Maroushkine, A. (2004). Room-Temperature Superconductivity (Cambridge
International Science, Cambridge).
Bibliography 205
89. Loram, J. W., Mirza, K. A., Cooper, J. R. and Tallon, J. L. (1998). Specific
heat evidence on the normal state pseudogap, J. Phys. Chem. Solids 59, 10–12,
pp. 2091–2094.
90. Revaz, B., Genond, J. Y., Junod, A., Naumdier, K., Erb, A. and Walker. E. (1998).
d-wave scaling relations in the mixed-state specific heat YBa2 Cu3 O7 , Phys. Rev.
Lett. 80, 15, pp. 3364–3367.
91. Hasselbach, K., Kirtley, J. R. and Flouquet, J. (1994). Symmetry of the gap in
superconducting UReSi, Phys. Rev. B 47, pp. 509–512.
92. Muhrer, G. and Schachinger, E. (2000). Free energy of a d-wave superconductor in
the nearly antiferromagnetic Fermi liquid model, J. Low Temp. Phys. 88–120, 1–2,
pp. 65–88.
93. Dorbolo, S., Ausloos, M. and Houssa, M. (1998). Electronic specific heat of super-
conductors with van Hove singularities: effects of a magnetic field and thermal fluc-
tuations, Phys. Rev. B 57, pp. 5401–5411.
94. Dobrosavljevic-Crujic, L. and Miranovic, P. (2003). Low temperature specific heat
in anisotropic superconductors, Physica C 397, pp. 117–122.
95. Bogolyubov, Jr., N. N. and Kruchinin, S. P. (2003). Modern approach to the calcu-
lation of the correlation function for superconductivity model, Mod. Phys. Lett. B
17, 10–12, pp. 709–724.
96. Nagao, H., Kruchinin, S. P., Yaremko, A. M. and Yamaguchi, K. (2002) Multiband
superconductivity, Int. J. Mod. Phys. B 16, pp. 3419–3428.
97. Kruchinin, S. P. and Nagao, H. (2005) Two-gap superconductivity in MgB, Phys.
Particles Nuclei 36, Suppl. 1, pp. 127–130.
98. Kruchinin, S. P. and Yaremko, A. M. (1998). Many zone effect in cuprate supercon-
ductors, Supercond. Sci. Technol. 11, pp. 4–8.
99. Nagao, H., Kawabe, H., Kruchinin, S. P., Manske, D. and Yamaguchi, K. (2003).
Theoretical studies on many band effects superconductivity by using renormalization
group approach, Mod. Phys. Lett. B 17, 10–12, pp. 423–431.
100. Nagao, H., Kawabe, H., Kruchinin, S. P. and Yamaguchi, K. (2003). Supercon-
ductivity in two-band model by renormalization group approach, Int. J. Mod. Phys.
B 17, pp. 3373–3376.
101. Yaremko, A. M., Mozdor, E. V. and Kruchinin, S. P. (1996). Coupled states of
electron–phonon system, superconductivity and HTSC of crystals, Int. J. Mod. Phys.
B 10, pp. 2665–2674.
102. Kruchinin, S. P. and Nagao, H. (2006). Multi-gap superconductivity in MgB2 , in
Proc. NATO Advanced Research Workshop “Symmetry and Heterogeneity in High-
Tc Superconductors”, ed. Bianconi, A. (Kluwer, Dorderecht), pp. 43–53.
103. Nagao, H., Yaremko, A., Kruchinin, S. P. and Yamaguchi, K. (2002). Many band
effects in superconductivity, in Proc. NATO Advanced Research Workshop “New
Trends in Superconductivity”, (eds.) Annett, J. and Kruchinin, S. P. (Kluwer,
Dordrecht), pp. 155–167.
104. Bardeen, J., Cooper, L. N. and Schrieffer, J. R. (1957). Microscopic theory of super-
conductivity, Phys. Rev. 108, pp. 1175.
105. Anderson, P. W. and Zou, Z. (1988). “Normal” tunneling and “normal” transport:
diagnostics for the resonating-valence-bond state, Phys. Rev. Lett. 60, p. 132.
106. Kampf, A. P. (1994). Magnetic correlations in high temperature superconductivity,
Phys. Rep. 249, p. 219.
107. Emery, V. J. (1987). Theory of high-Tc superconductivity in oxides, Phys. Rev. Lett.
58, p. 2794.
208 Modern Aspects of Superconductivity: Theory of Superconductivity
108. Hirsch, J. E. (1985). Attractive interaction and pairing in fermion systems with
strong on-site repulsion, Phys. Rev. Lett. 25, p. 1317.
109. Zhang, F. C. and Rice, T. M. (1988). Effective Hamiltonian for the superconducting
Cu oxides, Phys. Rev. B 37, p. 3759.
110. Fradkin, E. (1991). Field Theories of Condensed Matter Systems (Addison Wesley).
111. Nagaosa, N. and Lee, P. A. (1990). Normal-state properties of the uniform
resonating-valence-bond state, Phys. Rev. Lett. 64, p. 2450.
112. Baskaran, G., Zou, Z. and Anderson, P. W. (1987). The resonating valence bond
state and high-Tc superconductivity — a mean field theory, Solid State Commun.
63, p. 973.
113. Fukuyama, H. and Yoshida, K. (1987). Critical temperature of superconductivity
caused by strong correlation, Jpn. J. Appl. Phys. 26, pp. 371–373.
114. Yamamoto, S., Yamaguchi, K. and Nasu, K. (1990). Ab initio molecular-orbital study
on electron correlation effects in CuO6 clusters relating to high-Tc superconductivity,
Phys. Rev. B 42, pp. 266–272.
115. Nagao, H., Kitagawa, Y., Kawakami, T., Yoshimoto, T., Saito, H. and Yamaguchi, K.
(2001). Theoretical studies on field-induced superconductivity in molecular crystals,
Int. J. Quantum Chem. 85, pp. 608–618.
116. Kimura, M., Kawabe, H., Nakajima, A., Nishikawa, K. and Aono, S. (1988).
Possibility of superconducting and other phases in organic high polymers polyacene
carbon skeletons, Bull. Chem. Soc. Jpn. 61, p. 4239.
117. Torrance, J. B., Bezinge, A., Nazzal, A. I. and Parkin, S. S. (1989). Disappearance of
high temperature superconductivity induced by high carrier concentrations, Physica
C 162–164, p. 291.
118. Ord, T., Kristoffel, N. and Rago, K. (2003). MgB2 superconductivity properties in
a two-gap model, Mod. Phys. Lett. B 17, pp. 667–673.
119. Suhl H., Matthias, B. T. and Walker, R. (1959). Bardeen–Cooper–Schrieffer theory
of superconductivity in the case of overlapping bands, Phys. Rev. Lett. 3, p. 552.
120. Moskalenko, V. A. (1959). Superconductivity metals with overlapping energetic
bands, Fiz. Metalloved 8, p. 503.
121. Antipov, E. and Abakumov, A. (2008). Structural design of superconductors based
on complex copper oxides, Phys. Usp. 51, p. 180.
122. Müller, K. A. (2006). The search for new high temperature superconductors,
Supercond. Sci. Technol. 58, pp. 1–3.
123. Bennemann, K. H. and Ketterson, J. B. (2008). Superconductivity (Springer,
Heidelberg), vol. 1.
124. Nagamatsu, J., Nakamura, N., Muranaka, T., Zentani, Y. and Akimitsu, J. (2001).
Superconductivity at 39 K in magnezium deboride, Nature 410, p. 63.
125. Canfield, J. P. C., Bud’ko, S. L. and Finnemore, D. K. (2003). An overview of the
basic physical properties of M gB2 , Physica C 385, 1–2, pp. 1–3.
126. Binnig, G., Baratoff, A., Hoenig, H. E. and Bednorz, J. G. (1980). Two-band super-
conductivity in Nb-doped SrT iO3 , Phys. Rev. Lett. 45, p. 1352.
127. An, J. M. and Picket, W. E. (2001). Superconductivity of MgB2 : covalent bonds
driven metallic, Phys. Rev. Lett. 86, p. 4366.
128. Kortus, J., Mazin, I. I., Belashenko, K. D., Antropov, V. P. and Boyer, I. L. (2001).
Superconductivity of metallic boron in MgB2 , Phys. Rev. Lett. 86, p. 4656.
129. Kondo, J. (1963). Superconductivity in transition metals, Prog. Theor. Phys. 29,
p. 1.
130. Kristoffel, N., Ord, T. and Rago, K. (2003). MgB2 two-gap superconductivity with
intra- and interband couplings, Europhys. Lett. 61, pp. 109–115.
Bibliography 209
131. Szabo, P., Samuely, P., Klein, T., Marcus, J., Fruchart, D., Miraglia, S., Marcenat C.
and Jansen, A. (2001). Evidence for two superconducting energy gaps in M gB2 by
point-contact spectroscopy, Phys. Rev. Lett. 87, p. 137005.
132. McMillan, W. L. (1968). Transition temperature of strong-coupled superconductors,
Phys. Rev. 167, p. 331.
133. Konsin, P., Kristoffel, N. and Örd, T. (1988). The interband interaction as a possible
cause of high-temperature superconductivity, Phys. Lett. A 129, pp. 339–342.
134. Combescot, R. and Leyronas, X. (1995). Coupling between planes and chains in
YBa2 Cu3 O7 : a possible solution for the order parameter controversy, Phys. Rev.
Lett. 75, pp. 3732–3735.
135. Konsin, P. and Sorkin, B. (1998). Electric field effects in high-Tc cuprates, Phys.
Rev. B 58, pp. 5795–5802.
136. Nagao, H., Nishino, M., Shigeta, Y., Yoshioka, Y. and Yamaguchi, K. (2000).
Theoretical studies on superconducting and other phases: triplet superconductiv-
ity, ferromagnetism, and ferromagnetic metal, J. Chem. Phys. 113, pp. 721–732.
137. Kondo, J. (2001). Superconductivity of the two-dimensional Hubbard models with
a small U., J. Phys. Soc. Jpn. 70, p. 808.
138. Nagao, H., Mitani, M., Nishino, M., Shigeta, Y., Yoshioka, Y. and Yamaguchi, K.
(1999). Theoretical studies on anomalous phases in molecular systems with external
field: possibility of photo-induced superconductivity, Int. J. Quantum Chem. 75,
pp. 549–561.
139. Nagao, H., Nishino, M., Mitani, M., Yoshioka, Y. and Yamaguchi, K. (1997).
Possibilities of charge- and/or spin-mediated superconductors and photo-induced
superconductors in the intermediate region of metal-insulator transitions, Int. J.
Quantum Chem. 65, pp. 947–957.
140. Nagao, H., Mitani, M., Nishino, M., Shigeta, Y., Yoshioka, Y. and Yamaguchi, K.
(1988). Possibility of charge-mediated superconductors in the intermediate region of
metal–insulator transitions, Int. J. Quantum Chem. 70, pp. 1075–1084.
141. Bogoliubov, N. N. and Bogoliubov, N. N., Jr. (1986). An Introduction to Quantum
Statistical Mechanics (Nauka, Moscow) (transl. into English 1994).
142. Bogoliubov, N. N. and Shirkov, D. V. (1959). Introduction to the Theory of Quantized
Field (Interscience, New York).
143. Mattuck, R. D. (1976). A Guide to Feynman Diagrams in the Many-Body Problem
2nd edn. (McGraw-Hill, New York).
144. Kimura, M., Kawabe, H., Nishikawa, K. and S. Aono (1988). Possibility of supercon-
ducting and other phases in organic high polymers of polyacene carbon skeletons II:
screened electron–electron interaction, Bull. Chem. Soc. Jpn. 61, pp. 4245–4252.
145. Kimura, M., Kawabe, H., Nishikawa, K. and Aono, S. (1986). Superconducting and
other phases in organic high polymers of polyacenic carbon skeletons I: the method
of sum of divergent perturbation series, J. Chem. Phys. 85, pp. 3090–3096.
146. Kimura, M., Kawabe, H., Nishikawa, K. and Aono, S. (1985). Superconducting and
other phases in organic high polymers of polyacenic carbon skeletons II: the mean
field method, J. Chem. Phys. 85, pp. 3097–3100.
147. Kimura, M., Kawabe, H., Nakajima, A., Nishikawa, K. and Aono, S. (1988).
Possibility of superconducting and other phases in organic high polymers of poly-
acene carbon skeletons II: screened electron-electron interaction. Bull. Chem. Soc.
Jpn. 61, pp. 4239–4244.
148. Bogoliubov, N. N. and Tyablikov, S. (1959). Retarded and advanced Green func-
tions in statistical physics, Dokl. Acad. Sci. USSR 126, p. 53; Sov. Phys. Dokl. 4,
p. 589.
210 Modern Aspects of Superconductivity: Theory of Superconductivity
169. Matveev, K. A. and Larkin, A. L. (1997). Parity effect in ground state energies of
ultrasmall superconducting grains, Phys. Rev. Lett. 78, pp. 3749–3752.
170. Richardson, R. W. and Sherman, N. (1964). Pairing models of Pb206, Pb204 and
Pb202, Nucl. Phys. 523, pp. 253–268.
171. Richardson, R. W. (1977). Pairing in the limit of a large number of particles,
J. Math. Phys. 6, pp. 1802–1811.
172. Richardson, R. W. (1966). Numerical study of the 8–32-particle eigenstate of the
pairing Hamiltonian, Phys. Rev. 141, pp. 949–955.
173. Sierra, G., Dukelsky, J., Dussel, G. G., von Delft, J. and Braun, F. (2000). Exact
study of the effect of level statistics in ultrasmall superconducting grains, Phys. Rev.
B 61, pp. 11890–11893.
174. Kouwenhoven, L. and Glazman, L. (2001). The revival of the Kondo effect, Phys.
World 14, 1, pp. 33–38.
175. van der Wiel, G., De Franceschi, S., Fujisawa, T., Elzerman, J. M., Tarucha, S. and
Kouwenhoven, L. P. (2000). The Kondo effect in the unitary limit, Science 289,
pp. 2105–2108.
176. Inoshita, T. (1998). Kondo effect in quantum dots, Science 281, pp. 526–527.
177. Sasaki, S., De Franceschi, S., Elzerman, J. M., van der Wiel, W. G., Eto, M.,
Tarucha, S. and Kouwenhoven, L. P. (2000). Kondo effect in an integer-spin quantum
dot, Nature (London) 405, pp. 764–767.
178. Inoshita, T., Shimizu, A., Kuramoto, Y. and Sakaki, H. (1993). Correlated elec-
tron transport through a quantum dot: the multiple-level effect, Phys. Rev. B 48,
pp. 114725–114728.
179. Izumida, W., Sakai, O. and Shimizu, Y. (1998). Kondo effect in single quantum dot
systems study with numerical renormalization group method, J. Phys. Soc. Jpn. 67,
p. 2444.
180. Yeyati, A. L., Flores, F. and Martin-Rodero, A. (1999). Transport in multilevel
quantum dots: from the Kondo effect to the Coulomb blockade regime, Phys. Rev.
Lett. 83, pp. 600–603.
181. Pasupathy, A. N., Bialczak, R. C., Martinek, J., Grose, J. E., Donev, L. A. K.,
McEuen, P. L. and Ralph, D. C. (2004). The Kondo effect in the presence of ferro-
magnetism, Science 306, pp. 86–90.
182. Matsubayashi, D. and Eto, M. (2007). Spin splitting and Kondo effect in quantum
dots coupled to noncollinear ferromagnetic leads, Phys. Rev. B 75, p. 165319.
183. Buitelaar, M. R., Nussbaumer, T. and Schönenberger, C. (2002). Quantum dot in
the Kondo regime coupled to superconductors, Phys. Rev. Lett. 89, p. 256801.
184. Ueda, A. and Eto, M. (2006). Resonant tunneling and Fano resonance in quantum
dots with electron–phonon interaction, Phys. Rev. B 73, pp. 235353–235365.
185. Eto, M. and Nazarov, Y. V. (2000). Enhancement of Kondo effect in quantum dots
with an even number of electrons, Phys. Rev. Lett. 85, pp. 1306–1309.
186. Reich, S., Leitus, G., Popovitz-Biro, R. and Schechter, M. (2003). Magnetization of
small lead particles, Phys. Rev. Lett. 91, p. 147001.
187. Yamaji, K. and Shimoi, Y. (1994). Superconducting transition of the two-
chain Hubbard model indicated by diagonalization calculations, Physica C 222,
pp. 349–360.
188. Combescot, R. and Leyronas, X. (1995). Coupling between planes and chains in
YBa2 Cu3 O7 : a possible solution for the order parameter controversy, Phys. Rev.
Lett. 75, pp. 3732–3735.
212 Modern Aspects of Superconductivity: Theory of Superconductivity
209. Nakamura, Y., Pashkin, Y. A. and Tsai, J. S. (1999). Quantum oscillations in two
coupled charge qubits, Nature 398, pp. 786–788.
210. Makhlin, Y., Schon, G. and Shnirman, A. (2001). Quantum-state engineering with
Josephson-junction devices, Rev. Mod. Phys. 73, pp. 357–400.
211. Chiorescu, I., Nakamura, Y., Harmans, C. J. P. M. and Mooij, J. E. (2003). Coherent
quantum dynamics of a superconducting flux qubit, Science 299, pp. 1869–1871.
INDEX
π junction 88 167
acoustic phonons 74 79
anharmonic model 72 76
antiferromagnetic 62 64 70 71
88 93 101 107
111 138 193
anyonic superconductivity 78
basis functions 81 82
BCS theory 27 32 44 51
64 65 72 78
79 88 89 104
114 116 135 137
140 142 145 161
163
bipolaronic mechanism 72 78
Bose condensation 80
Bose particles 78
boson field 44 47 48 53
55 59
gap equation 20 21 24 25
50 125 140 141
Gauss theorem 170 201
Ginzburg–Landau equation 58 59
Green functions 67 94 113 114
139 197
instability 55 90
interacting 15 36 39 42
168
interaction constant 138 142 144
interferometer 185
Jahn–Teller effect 77
Josephson junction 3 84 106 182
184
jump of the heat 89 104 107
Knight shift 84
Kondo effect 68 115 157 158
160 165 166 197
Lagrangian 3 7 44 53
55 56 94
length scale 168 169
level spacing 135 141 161 165
166
London equation 168
London superconductors 166 168 199
LUMO 25
magnetic resonance 28 31 32
magnons 74 90
mechanism of pairing 62 63 72 74
77 78 80 93
106 126 127
mesoscopic superconductivity 135
Mg1-xAlxB2 121 122
mixed state 81 82
multiband superconductivity 61 108 135 190
Néel temperature 70 71
Nambu spinor 4 7 29 39
Nambu–Goldstone theorem 7
nanosize two-gap superconductivity 137 138 141 143
NMR 31 84 91 92
Noether theorem 7
propagator 10 11 16 17
19 33 37 41
42 53
pseudogap 70 71
self-energy 33 37 38 42
170
size effect 136
spin dynamics 41 92
theory of superconductivity 1 2 32 43
62 80 88 108
thermodynamics 62 88 89 93
107 111 138 193
Tl-based 67 68
two-dimensional 59 76 77 81
84 90 101 118
119 174
two-gap superconductivity 65 66 108 112
113 119 127 137
143 157 194 196
This page has been reformatted by Knovel to provide easier navigation.
Index Terms Links
yttrium ceramics 66 67