Solid State Physics

Solid State Physics
Chetan Nayak
Physics 140a
Franz 1354; M, W 11:00-12:15
Office Hour: TBA; Knudsen 6-130J
Section: MS 7608; F 11:00-11:50
TA: Sumanta Tewari
University of California,
Los Angeles
September 2000
Contents
1 What is Condensed Matter Physics? 1

1.1 Length, time, energy scales . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Microscopic Equations vs. States of Matter, Phase Transitions, Critical
Points . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Broken Symmetries . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.4 Experimental probes: X-ray scattering, neutron scattering, NMR, ther-
modynamic, transport . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.5 The Solid State: metals, insulators, magnets, superconductors . . . . 4
1.6 Other phases: liquid crystals, quasicrystals, polymers, glasses . . . . . 5
2 Review of Quantum Mechanics 7

2.1 States and Operators . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.2 Density and Current . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.3 δ-function scatterer . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.4 Particle in a Box . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.5 Harmonic Oscillator . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.6 Double Well . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.7 Spin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.8 Many-Particle Hilbert Spaces: Bosons, Fermions . . . . . . . . . . . . 15
3 Review of Statistical Mechanics 18
ii
3.1 Microcanonical, Canonical, Grand Canonical Ensembles . . . . . . . . 18
3.2 Bose-Einstein and Planck Distributions . . . . . . . . . . . . . . . . . 21
3.2.1 Bose-Einstein Statistics . . . . . . . . . . . . . . . . . . . . . . 21
3.2.2 The Planck Distribution . . . . . . . . . . . . . . . . . . . . . 22
3.3 Fermi-Dirac Distribution . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.4 Thermodynamics of the Free Fermion Gas . . . . . . . . . . . . . . . 24
3.5 Ising Model, Mean Field Theory, Phases . . . . . . . . . . . . . . . . 27
4 Broken Translational Invariance in the Solid State 30

4.1 Simple Energetics of Solids . . . . . . . . . . . . . . . . . . . . . . . . 30
4.2 Phonons: Linear Chain . . . . . . . . . . . . . . . . . . . . . . . . . . 31
4.3 Quantum Mechanics of a Linear Chain . . . . . . . . . . . . . . . . . 31
4.3.1 Statistical Mechnics of a Linear Chain . . . . . . . . . . . . . 36
4.4 Lessons from the 1D chain . . . . . . . . . . . . . . . . . . . . . . . . 37
4.5 Discrete Translational Invariance: the Reciprocal Lattice, Brillouin
Zones, Crystal Momentum . . . . . . . . . . . . . . . . . . . . . . . . 38
4.6 Phonons: Continuum Elastic Theory . . . . . . . . . . . . . . . . . . 40
4.7 Debye theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
4.8 More Realistic Phonon Spectra: Optical Phonons, van Hove Singularities 46
4.9 Lattice Structures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
4.9.1 Bravais Lattices . . . . . . . . . . . . . . . . . . . . . . . . . . 48
4.9.2 Reciprocal Lattices . . . . . . . . . . . . . . . . . . . . . . . . 50
4.9.3 Bravais Lattices with a Basis . . . . . . . . . . . . . . . . . . 51
4.10 Bragg Scattering . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
5 Electronic Bands 57
5.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
5.2 Independent Electrons in a Periodic Potential: Bloch’s theorem . . . 57
iii
5.3 Tight-Binding Models . . . . . . . . . . . . . . . . . . . . . . . . . . 59
5.4 The δ-function Array . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
5.5 Nearly Free Electron Approximation . . . . . . . . . . . . . . . . . . 66
5.6 Some General Properties of Electronic Band Structure . . . . . . . . 68
5.7 The Fermi Surface . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
5.8 Metals, Insulators, and Semiconductors . . . . . . . . . . . . . . . . . 71
5.9 Electrons in a Magnetic Field: Landau Bands . . . . . . . . . . . . . 74
5.9.1 The Integer Quantum Hall Effect . . . . . . . . . . . . . . . . 79
iv
Chapter 1
What is Condensed Matter

Physics?
1.1 Length, time, energy scales

We will be concerned with:
• ω, T 1eV
• |xi − xj |, 1q 1Å
as compared to energies in the M eV for nuclear matter, and GeV or even T eV , in

particle physics.
The properties of matter at these scales is determined by the behavior of collections
of many (∼ 1023 ) atoms.
In general, we will be concerned with scales much smaller than those at which
gravity becomes very important, which is the domain of astrophysics and cosmology.
1
Chapter 1: What is Condensed Matter Physics? 2
1.2 Microscopic Equations vs. States of Matter,

Phase Transitions, Critical Points
Systems containing many particles exhibit properties which are special to such sys-
tems. Many of these properties are fairly insensitive to the details at length scales
shorter than 1Å and energy scales higher than 1eV – which are quite adequately
described by the equations of non-relativistic quantum mechanics. Such properties
are emergent. For example, precisely the same microscopic equations of motion –
Newton’s equations – can describe two different systems of 1023 H2 O molecules.
d2~xi X
~ i V (~xi − ~xj )
m 2 =− ∇ (1.1)
dt j6=i
Or, perhaps, the Schrödinger equation:

 
2
h̄ X
∇2i +
X
− V (~xi − ~xj ) ψ (~x1 , . . . , ~xN ) = E ψ (~x1 , . . . , ~xN ) (1.2)
2m i i,j
However, one of these systems might be water and the other ice, in which case the
properties of the two systems are completely different, and the similarity between
their microscopic descriptions is of no practical consequence. As this example shows,
many-particle systems exhibit various phases – such as ice and water – which are
not, for the most part, usefully described by the microscopic equations. Instead, new
low-energy, long-wavelength physics emerges as a result of the interactions among
large numbers of particles. Different phases are separated by phase transitions, at
which the low-energy, long-wavelength description becomes non-analytic and exhibits
singularities. In the above example, this occurs at the freezing point of water, where
its entropy jumps discontinuously.
1.3 Broken Symmetries

As we will see, different phases of matter are distinguished on the basis of symmetry.
The microscopic equations are often highly symmetrical – for instance, Newton’s
laws are translationally and rotationally invariant – but a given phase may exhibit
much less symmetry. Water exhibits the full translational and rotational symmetry of
of Newton’s laws; ice, however, is only invariant under the discrete translational and
rotational group of its crystalline lattice. We say that the translational and rotational
symmetries of the microscopic equations have been spontaneously broken.
1.4 Experimental probes: X-ray scattering, neu-

tron scattering, NMR, thermodynamic, trans-
port
There are various experimental probes which can allow an experimentalist to deter-
mine in what phase a system is and to determine its quantitative properties:
• Scattering: send neutrons or X-rays into the system with prescribed energy,
momentum and measure the energy, momentum of the outgoing neutrons or
X-rays.
• NMR: apply a static magnetic field, B, and measure the absorption and emission
by the system of magnetic radiation at frequencies of the order of ωc = geB/m.
Essentially the scattering of magnetic radiation at low frequency by a system
in a uniform B field.
• Thermodynamics: measure the response of macroscopic variables such as the

energy and volume to variations of the temperature, pressure, etc.
• Transport: set up a potential or thermal gradient, ∇ϕ, ∇T and measure the

electrical or heat current ~j, ~jQ . The gradients ∇ϕ, ∇T can be held constant or
made to oscillate at finite frequency.
1.5 The Solid State: metals, insulators, magnets,

superconductors
In the solid state, translational and rotational symmetries are broken by the arrange-
ment of the positive ions. It is precisely as a result of these broken symmetries that
solids are solid, i.e. that they are rigid. It is energetically favorable to break the
symmetry in the same way in different parts of the system. Hence, the system resists
attempts to create regions where the residual translational and rotational symmetry
groups are different from those in the bulk of the system. The broken symmetry can
be detected using X-ray or neutron scattering: the X-rays or neutrons are scattered
by the ions; if the ions form a lattice, the X-rays or neutrons are scattered coherently,
forming a diffraction pattern with peaks. In a crystalline solid, discrete subgroups of
the translational and rotational groups are preserved. For instance, in a cubic lattice,
rotations by π/2 about any of the crystal axes are symmetries of the lattice (as well
as all rotations generated by products of these). Translations by one lattice spacing
along a crystal axis generate the discrete group of translations.
In this course, we will be focussing on crystalline solids. Some examples of non-
crystalline solids, such as plastics and glasses will be discussed below. Crystalline
solids fall into three general categories: metals, insulators, and superconductors. In
addition, all three of these phases can be further subdivided into various magnetic
phases. Metals are characterized by a non-zero conductivity at T = 0. Insulators
have vanishing conductivity at T = 0. Superconductors have infinite conductivity for
T < Tc and, furthermore, exhibit the Meissner effect: they expel magnetic fields.
In a magnetic material, the electron spins can order, thereby breaking the spin-
rotational invariance. In a ferromagnet, all of the spins line up in the same direction,
thereby breaking the spin-rotational invariance to the subgroup of rotations about this
direction while preserving the discrete translational symmetry of the lattice. (This
can occur in a metal, an insulator, or a superconductor.) In an antiferromagnet,
neighboring spins are oppositely directed, thereby breaking spin-rotational invariance
to the subgroup of rotations about the preferred direction and breaking the lattice
translational symmetry to the subgroup of translations by an even number of lattice
sites.
Recently, new states of matter – the fractional quantum Hall states – have been
discovered in effectively two-dimensional systems in a strong magnetic field at very
low T . Tomorrow’s experiments will undoubtedly uncover new phases of matter.
1.6 Other phases: liquid crystals, quasicrystals,

polymers, glasses
The liquid – with full translational and rotational symmetry – and the solid – which
only preserves a discrete subgroup – are but two examples of possible realizations
of translational symmetry. In a liquid crystalline phase, translational and rotational
symmetry is broken to a combination of discrete and continuous subgroups. For
instance, a nematic liquid crystal is made up of elementary units which are line seg-
ments. In the nematic phase, these line segments point, on average, in the same
direction, but their positional distribution is as in a liquid. Hence, a nematic phase
breaks rotational invariance to the subgroup of rotations about the preferred direction
and preserves the full translational invariance. Nematics are used in LCD displays.
In a smectic phase, on the other hand, the line segments arrange themselves into
layers, thereby partially breaking the translational symmetry so that discrete transla-
tions perpendicular to the layers and continuous translations along the layers remain
unbroken. In the smectic-A phase, the preferred orientational direction is the same
as the direction perpendicular to the layers; in the smectic-C phase, these directions
are different. In a hexatic phase, a two-dimensional system has broken orientational
order, but unbroken translational order; locally, it looks like a triangular lattice. A
quasicrystal has rotational symmetry which is broken to a 5-fold discrete subgroup.
Translational order is completely broken (locally, it has discrete translational order).
Polymers are extremely long molecules. They can exist in solution or a chemical re-
action can take place which cross-links them, thereby forming a gel. A glass is a rigid,
‘random’ arrangement of atoms. Glasses are somewhat like ‘snapshots’ of liquids, and
are probably non-equilibrium phases, in a sense.
Chapter 2
Review of Quantum Mechanics
2.1 States and Operators

E
A quantum mechanical system is defined by a Hilbert space, H, whose vectors, ψ

are associated with the states of the system. A state of the system is represented by
E
the set of vectors eiα ψ . There are linear operators, Oi which act on this Hilbert

space. These operators correspond to physical observables. Finally, there is an inner

D E E E
product, which assigns a complex number, χψ , to any pair of states, ψ , χ . A
E
state vector, ψ gives a complete description of a system through the expectation

D E E D E
values, ψ Oi ψ (assuming that ψ is normalized so that ψ ψ = 1), which would

be the average values of the corresponding physical observables if we could measure

E
them on an infinite collection of identical systems each in the state ψ .

The adjoint, O† , of an operator is defined according to

D E D E
χ Oψ = χO† ψ (2.1)
E E
In other words, the inner product between χ and Oψ is the same as that between
E E
O† χ and ψ . An Hermitian operator satisfies

O = O† (2.2)
7
Chapter 2: Review of Quantum Mechanics 8
while a unitary operator satisfies
OO† = O† O = 1 (2.3)
If O is Hermitian, then
eiO (2.4)
is unitary. Given an Hermitian operator, O, its eigenstates are orthogonal,

D E D E D E
λ0 Oλ = λ λ0 λ = λ0 λ0 λ (2.5)
For λ 6= λ0 ,
D E
λ0 λ = 0 (2.6)

If there are n states with the same eigenvalue, then, within the subspace spanned by
these states, we can pick a set of n mutually orthogonal states. Hence, we can use
E
E
the eigenstates λ as a basis for Hilbert space. Any state ψ can be expanded in
the basis given by the eigenstates of O:
E X E
ψ = c λ λ (2.7)

λ
with
D E
cλ = λψ (2.8)
A particularly important operator is the Hamiltonian, or the total energy, which

we will denote by H. Schrödinger’s equation tells us that H determines how a state
of the system will evolve in time.
∂ E E
ih̄ ψ = H ψ (2.9)

∂t
If the Hamiltonian is independent of time, then we can define energy eigenstates,

E E
H E = E E (2.10)

which evolve in time according to:

E Et
E
E(t) = e−i h̄ E(0) (2.11)

An arbitrary state can be expanded in the basis of energy eigenstates:

E X E
ψ = ci Ei (2.12)

i
It will evolve according to:

E Ej t E
cj e−i h̄ Ej
X
ψ(t) = (2.13)

j
For example, consider a particle in 1D. The Hilbert space consists of all continuous
complex-valued functions, ψ(x). The position operator, x̂, and momentum operator,
p̂ are defined by:
x̂ · ψ(x) ≡ x ψ(x)
∂
p̂ · ψ(x) ≡ −ih̄ ψ(x) (2.14)
∂x
The position eigenfunctions,
x δ(x − a) = a δ(x − a) (2.15)
are Dirac delta functions, which are not continuous functions, but can be defined as
the limit of continuous functions:
1 x2
δ(x) = lim √ e− a2 (2.16)
a→0 a π
The momentum eigenfunctions are plane waves:

∂ ikx
−ih̄ e = h̄k eikx (2.17)
∂x
Expanding a state in the basis of momentum eigenstates is the same as taking its
Fourier transform:
Z ∞ 1
ψ(x) = dk ψ̃(k) √ eikx (2.18)
−∞ 2π
where the Fourier coefficients are given by:

1 Z∞
ψ̃(k) = √ dx ψ(x) e−ikx (2.19)
2π −∞
If the particle is free,

h̄2 ∂ 2
H=− (2.20)
2m ∂x2
then momentum eigenstates are also energy eigenstates:
h̄2 k 2 ikx
Ĥeikx = e (2.21)
2m
If a particle is in a Gaussian wavepacket at the origin at time t = 0,
1 x2
ψ(x, 0) = √ e− a2 (2.22)
a π
Then, at time t, it will be in the state:
1 Z∞ a h̄k2 t 1 2 2
ψ(x, t) = √ dk √ e−i 2m e− 2 k a eikx (2.23)
2π −∞ π
2.2 Density and Current

Multiplying the free-particle Schrödinger equation by ψ ∗ ,
∂ h̄2 ∗ ∂ 2
ψ ∗ ih̄ ψ=− ψ ψ (2.24)
∂t 2m ∇2
and subtracting the complex conjugate of this equation, we find
∂ ih̄ ~ ∗ ~
~ ∗ ψ
(ψ ∗ ψ) = ∇ · ψ ∇ψ − ∇ψ (2.25)
∂t 2m
This is in the form of a continuity equation,
∂ρ ~ · ~j
=∇ (2.26)
∂t
The density and current are given by:
ρ = ψ∗ψ
~j = ih̄ ∗ ~
~ ∗ ψ
ψ ∇ψ − ∇ψ (2.27)
2m
The current carried by a plane-wave state is:
~j = h̄ ~k 1 (2.28)
2m (2π)3
2.3 δ-function scatterer

h̄2 ∂ 2
H=− + V δ(x) (2.29)
2m ∂x2

 eikx + Re−ikx

if x < 0
ψ(x) = (2.30)
 T eikx

if x > 0
1
T = mV
1− h̄2 k
i
mV
h̄2 k
i
R = mV (2.31)
1 − h̄2 k i
There is a bound state at:

mV
ik = (2.32)
h̄2
2.4 Particle in a Box

Particle in a 1D region of length L:
h̄2 ∂ 2
H=− (2.33)
2m ∂x2
ψ(x) = Aeikx + Be−ikx (2.34)
has energy E = h̄2 k 2 /2m. ψ(0) = ψ(L) = 0. Therefore,
nπ

ψ(x) = A sin x (2.35)
L
for integer n. Allowed energies

h̄2 π 2 n2
En = (2.36)
2mL2
In a 3D box of side L, the energy eigenfunctions are:
nx π ny π nz π

ψ(x) = A sin x sin y sin z (2.37)
L L L
and the allowed energies are:
h̄2 π 2 2 2 2

En = n x + n y + n z (2.38)
2mL2
2.5 Harmonic Oscillator

h̄2 ∂ 2 1
H=− 2
+ kx2 (2.39)
2m ∂x 2
q
Writing ω = k/m, p̃ = p/(km)1/4 , x̃ = x(km)1/4 ,
1 2
H= ω p̃ + x̃2 (2.40)
2
[p̃, x̃] = −ih̄ (2.41)
Raising and lowering operators:

√
a = (x̃ + ip̃) / 2h̄
√
a† = (x̃ − ip̃) / 2h̄
(2.42)
Hamiltonian and commutation relations:

1

†
H = h̄ω a a +
2
[a, a† ] = 1 (2.43)
The commutation relations,
[H, a† ] = h̄ωa†
[H, a] = −h̄ωa (2.44)
imply that there is a ladder of states,
Ha† |Ei = (E + h̄ω) a† |Ei

Ha|Ei = (E − h̄ω) a|Ei (2.45)
This ladder will continue down to negative energies (which it can’t since the Hamil-
tonian is manifestly positive definite) unless there is an E0 ≥ 0 such that
a|E0 i = 0 (2.46)
Such a state has E0 = h̄ω/2.

We label the states by their a† a eigenvalues. We have a complete set of H eigen-
states, |ni, such that
1

H|ni = h̄ω n + |ni (2.47)
2
and (a† )n |0i ∝ |ni. To get the normalization, we write a† |ni = cn |n + 1i. Then,
|cn |2 = hn|aa† |ni

= n+1 (2.48)
Hence,
√
a† |ni = n + 1|n + 1i
√
a|ni = n|n − 1i (2.49)
2.6 Double Well

h̄2 ∂ 2
H=− + V (x) (2.50)
2m ∂x2
where 
∞ if |x| > 2a + 2b






V (x) = 0 if b < |x| < a + b



if |x| < b

 V0

Symmetrical solutions:

 A cos k 0 x

if |x| < b
ψ(x) =  (2.51)
 cos(k|x| − φ) if b < |x| < a + b
with s
0 2mV0
k = k2 − (2.52)
h̄2
The allowed k’s are determined by the condition that ψ(a + b) = 0:
1

φ= n+ π − k(a + b) (2.53)
2
the continuity of ψ(x) at |x| = b:
cos (kb − φ)
A= (2.54)
cos k 0 b
and the continuity of ψ 0 (x) at |x| = b:
1

k tan n+ π − ka = k 0 tan k 0 b (2.55)
2
If k 0 is imaginary, cos → cosh and tan → i tanh in the above equations.
Antisymmetrical solutions:

 A sin k 0 x

if |x| < b
ψ(x) = (2.56)
 sgn(x) cos(k|x| − φ) if b < |x| < a + b

The allowed k’s are now determined by

1

φ= n+ π − k(a + b) (2.57)
2
cos (kb − φ)
A= (2.58)
sin k 0 b
1

k tan n+ π − ka = − k 0 cot k 0 b (2.59)
2
Suppose we have n wells? Sequences of eigenstates, classified according to their
eigenvalues under translations between the wells.
2.7 Spin
The electron carries spin-1/2. The spin is described by a state in the Hilbert space:
α|+i + β|−i (2.60)
spanned by the basis vectors |±i. Spin operators:

 
1  0 1 
sx =
2
 
1 0
 
1  0 −i 
sy =
2
 
i 0
 
1  1 0 
sz = (2.61)
2
 
0 −1
Coupling to an external magnetic field:
~
Hint = −gµB ~s · B (2.62)
States of a spin in a magnetic field in the ẑ direction:

g
H|+i = − µB |+i
g2
H|−i = µB |−i (2.63)
2
2.8 Many-Particle Hilbert Spaces: Bosons, Fermions

When we have a system with many particles, we must now specify the states of all
of the particles. If we have two distinguishable particles whose Hilbert spaces are
spanned by the bases

E
i, 1 (2.64)

and
E
α, 2 (2.65)

Then the two-particle Hilbert space is spanned by the set:

E E E
i, 1; α, 2 ≡ i, 1 ⊗ α, 2 (2.66)

Suppose that the two single-particle Hilbert spaces are identical, e.g. the two particles
are in the same box. Then the two-particle Hilbert space is:
E E E
i, j ≡ i, 1 ⊗ j, 2 (2.67)

E E
If the particles are identical, however, we must be more careful. i, j and j, i must
be physically the same state, i.e.
E E
i, j = eiα j, i (2.68)

Applying this relation twice implies that

E E
i, j = e2iα i, j (2.69)

so eiα = ±1. The former corresponds to bosons, while the latter corresponds to
fermions. The two-particle Hilbert spaces of bosons and fermions are respectively
spanned by:
E E
i, j + j, i (2.70)

and
E E
i, j − j, i (2.71)

The n-particle Hilbert spaces of bosons and fermions are respectively spanned by:
X E
iπ(1) , . . . , iπ(n) (2.72)
π
and
E
(−1)π iπ(1) , . . . , iπ(n)
X
(2.73)

π
In position space, this means that a bosonic wavefunction must be completely sym-
metric:
ψ(x1 , . . . , xi , . . . , xj , . . . , xn ) = ψ(x1 , . . . , xj , . . . , xi , . . . , xn ) (2.74)
while a fermionic wavefunction must be completely antisymmetric:
ψ(x1 , . . . , xi , . . . , xj , . . . , xn ) = −ψ(x1 , . . . , xj , . . . , xi , . . . , xn ) (2.75)

Chapter 3
Review of Statistical Mechanics
3.1 Microcanonical, Canonical, Grand Canonical

Ensembles
In statistical mechanics, we deal with a situation in which even the quantum state
of the system is unknown. The expectation value of an observable must be averaged
over:
X
hOi = wi hi |O| ii (3.1)
i
where the states |ii form an orthonormal basis of H and wi is the probability of being
in state |ii. The wi ’s must satisfy
P
wi = 1. The expectation value can be written in
a basis-independent form:
hOi = T r {ρO} (3.2)
i wi |iihi|.
P
where ρ is the density matrix. In the above example, ρ = The condition,
P
wi = 1, i.e. that the probabilities add to 1, is:
T r {ρ} = 1 (3.3)
We usually deal with one of three ensembles: the microcanonical emsemble, the
canonical ensemble, or the grand canonical ensemble. In the microcanonical ensemble,
18
Chapter 3: Review of Statistical Mechanics 19
we assume that our system is isolated, so the energy is fixed to be E, but all states
with energy E are taken with equal probability:
ρ = C δ(H − E) (3.4)
C is a normalization constant which is determined by (3.3). The entropy is given by,
S = − ln C (3.5)
In other words,

S(E) = ln # of states with energy E (3.6)
Inverse temperature, β = 1/(kB T ):

!
∂S
β≡ (3.7)
∂E V
Pressure, P : !
P ∂S
≡ (3.8)
kB T ∂V E
where V is the volume.
First law of thermodynamics:
∂S ∂S
dS = dE + dV (3.9)
∂E ∂V
dE = kB T dS − P dV (3.10)
Free energy:
F = E − kB T S (3.11)
Differential relation:
dF = −kB S dT − P dV (3.12)
or, !
1 ∂F
S=− (3.13)
kB ∂T V
!
∂F
P =− (3.14)
∂V T
while
E = F + kB T S !
∂F
= F −T
∂T V
∂ F
= −T2 (3.15)
∂T T
In the canonical ensemble, we assume that our system is in contact with a heat
reservoir so that the temperature is constant. Then,
ρ = C e−βH (3.16)
It is useful to drop the normalization constant, C, and work with an unnormalized

density matrix so that we can define the partition function:
Z = T r {ρ} (3.17)
or,
e−βEa
X
Z= (3.18)
a
The average energy is:

1X
E = Ea e−βEa
Z a
∂
= − ln Z
∂β
∂
= − kB T 2 ln Z (3.19)
∂T
Hence,
F = −kB T ln Z (3.20)
The chemical potential, µ, is defined by

∂F
µ= (3.21)
∂N
where N is the particle number.

In the grand canonical ensemble, the system is in contact with a reservoir of heat
and particles. Thus, the temperature and chemical potential are held fixed and
ρ = C e−β(H−µN ) (3.22)
We can again work with an unnormalized density matrix and construct the grand
canonical partition function:
e−β(Ea −µN )
X
Z= (3.23)
N,a
The average number is:

∂
N = −kB T ln Z (3.24)
∂µ
while the average energy is:
∂ ∂
E=− ln Z + µkB T ln Z (3.25)
∂β ∂µ
3.2 Bose-Einstein and Planck Distributions
3.2.1 Bose-Einstein Statistics
For a system of free bosons, the partition function
e−β(Ea −µN )
X
Z= (3.26)
Ea ,N
can be rewritten in terms of the single-particle eigenstates and the single-particle

energies i :
Ea = n0 0 + n1 1 + . . . (3.27)
P P
e−β ( n −µ n )
X
Z = i i i i i
{ni }
!
−β(ni i −µni )
Y X
= e
i ni
Y 1
= (3.28)
i 1 − e−β(i −µ)
1
hni i = (3.29)
eβ(i −µ) − 1
The chemical potential is chosen so that
X
N = hni i
i
X 1
= (3.30)
i eβ(i −µ) −1
The energy is given by
X
E = hni i i
i
X i
= β( −µ)
(3.31)
i e i −1
N is increased by increasing µ (µ ≤ 0 always). Bose-Einstein condensation occurs

when
X
N> hni i (3.32)
i6=0
In such a case, hn0 i must become large. This occurs when µ = 0.
3.2.2 The Planck Distribution
Suppose N is not fixed, but is arbitrary, e.g. the numbers of photons and neutrinos
are not fixed. Then there is no Lagrange multiplier µ and
1
hni i = (3.33)
eβi − 1
Consider photons (two polarizations) in a cavity of side L with k = h̄ωk = h̄ck and
2π
k= (mx , my , mz ) (3.34)
L
X D E
E = 2 ωmx ,my ,mz nmx ,my ,mz (3.35)
mx ,my ,mz
We can take the thermodynamic limit, L → ∞, and convert the sum into an
integral. Since the allowed ~k’s are 2π
L
(mx , my , mz ), the ~k-space volume per allowed ~k
is (2π)3 /L3 . Hence, we can take the infinite-volume limit by making the replacement:
1 X ~
f (~k) = f (k) (∆~k)3
X
k
~ 3
(∆k) k
L3 Z 3~ ~
= d k f (k) (3.36)
(2π)3
Hence,
Z
d3 k
kmax h̄ωk
E = 2V
0 (2π)3 eβh̄ωk − 1
Z kmax 3
dk h̄ck
= 2V 3 βh̄ck
0 (2π) e −1
4 Z βh̄ckmax 3
V kB 4 x dx
= T (3.37)
π 2 (h̄c)3 0 ex − 1
For βh̄ckmax 1,
4 Z ∞ 3
V kB 4 x dx
E= 2 3
T (3.38)
π (h̄c) 0 ex − 1
and
4 Z ∞ 3
4V kB 3 x dx
CV = 2 3
T (3.39)
π (h̄c) 0 ex − 1
For βh̄ckmax 1 ,
3
V kmax
E= kB T (3.40)
3π 2
and
3
V kmax kB
CV = (3.41)
3π 2
3.3 Fermi-Dirac Distribution

For a system of free fermions, the partition function
e−β(Ea −µN )
X
Z= (3.42)
Ea ,N
can again be rewritten in terms of the single-particle eigenstates and the single-particle
energies i :
Ea = n0 0 + n1 1 + . . . (3.43)
but now
ni = 0, 1 (3.44)
so that
P P
e−β ( n −µ n )
X
Z = i i i i i
{ni }
 
1
e−β(ni i −µni ) 
Y X
= 
i ni =0
Y
= 1 + e−β(i −µ) (3.45)
i
1
hni i = (3.46)
eβ(i +µ) +1
The chemical potential is chosen so that
X 1
N= (3.47)
i eβ(i +µ) + 1
The energy is given by

X i
E= β( +µ)
(3.48)
i e i +1
3.4 Thermodynamics of the Free Fermion Gas

Free electron gas in a box of side L:
h̄2 k 2
k = (3.49)
2m
with
2π
k= (mx , my , mz ) (3.50)
L
Then, taking into account the 2 spin states,

X D E
E = 2 mx ,my ,mz nmx ,my ,mz
mx ,my ,mz
h̄2 k2
Z kmax d3 k 2m
= 2V
0 (2π)3 β h̄2 k2
2m
−µ
e +1
(3.51)
Z kmax d3 k 1
N = 2V
0 (2π)3 β h̄2 k2
2m
−µ
e +1
(3.52)
At T = 0,
h̄2 k 2
!
1
=θ µ− (3.53)
β h̄2 k2
2m
−µ 2m
e +1
All states with energies less than µ are filled; all states with higher energies are empty.
We write √
2mµT =0
kF = , F = µT =0 (3.54)
h̄
Z kF 3
N dk kF3
=2 = (3.55)
V 0 (2π)3 3π 2
E Z kF 3
d k h̄2 k 2
= 2
V 0 (2π)3 2m
2 5
1 h̄ kF
=
π 2 10m
3N
= F (3.56)
5 V
3 1
Z
d3 k m2 22 Z 1
2 3
= 2 3 d 2 (3.57)
(2π) π h̄
For kB T eF ,
3 1
N m2 22 Z ∞ 1 1
= 3 d 2
V π 23h̄ 1 0 eβ(−µ) +1
3 1 Z 3 1
m 2 Z µ
m 2 µ 1 m2 22 Z ∞ 1

2 2 1 2 2 1 1
= 2 3 d + 2 3
2 d 2
β(−µ)
−1 + 2 3 d 2 β(−µ)
π h̄ 0 π h̄ 0 e +1 π h̄ µ e +1
3 3 1 3 1
(2m) 2 3 m 2 2 2 Z µ 1 1 m2 22 Z ∞ 1 1
= µ 2 − d 2 + d 2
3 3 3
3π 2 h̄ π 2 h̄ 0 e−β(−µ) + 1 π 2 h̄ µ eβ(−µ) + 1
3 3 1 Z
(2m) 2 3 m 2 2 2 ∞ kB T dx 1 1
2 − (µ − k T x) 2 −βµ
= µ 2 + (µ + k B T x) B + O e
3π 2 h̄3 π 2 h̄3 0 ex + 1  
3 3 ∞ 3
(2m) 2 3 (2m) 2 X 3 1 Γ 2
Z ∞
x 2n−1
= µ2 + 2 3 (kB T )2n µ 2 −2n dx x 
3π 2 h̄3

π h̄ n=1 (2n − 1)! Γ 5 − 2n 0 e +1
 2
3 !2
(2m) 2 3  3 kB T
= 2 3 µ
2 1+ I1 + O T 4  (3.58)
3π h̄ 2 µ
with
Z ∞ xk
Ik = dx (3.59)
0 ex + 1
We will only need
π2
I1 = (3.60)
12
Hence,
 !2 
3 3 3 kB T
(F ) = µ 1 +
2 2 I1 + O T 4  (3.61)
2 µ
To lowest order in T , this gives:

 !2 
kB T
µ = F 1 − I1 + O(T 4 )
F
 !2 
π2 kB T
= F 1 − + O(T 4 ) (3.62)
12 F
3 1
E m2 22 Z ∞ 3 1
= 3 d 2
β(−µ)
V π 23h̄ 1 0 e 3 1+ 1 3 1
m 2 Z µ
m2 22 Z µ 1 m2 22 Z ∞ 1

2 2 3 3 3
= 2 3 d + 2 3
2 d 2
β(−µ)
−1 + 2 3 d 2 β(−µ)
π h̄ 0 π h̄ 0 e +1 π h̄ µ e +1
3 3 1 3 1
(2m) 2 5 m 2 2 2 Z µ 3 1 m2 22 Z ∞ 3 1
= µ2 − 2 3 d 2 −β(−µ) + 2 3 d 2 β(−µ)
5π 2 h̄3 π h̄ 0 e +1 π h̄ µ e +1
3 3 1 Z
(2m) 2 5 m 2 2 2 ∞ kB T dx 3 3
2 − (µ − k T x) 2 −βµ
= µ 2 + (µ + k B T x) B + O e
5π 2 h̄3 π 2 h̄3 0 ex + 1  
3 3 ∞ 5
(2m) 2 5 (2m) 2 X 5 1 Γ 2
Z ∞
x 2n−1
= µ2 + 2 3 (kB T )2n µ 2 −2n dx x 
5π 2 h̄3

π h̄ n=1 (2n − 1)! Γ 7 − 2n 0 e +1
2
3
 !2 
(2m) 2 5  15 kB T 4
= 2 3µ
2 1+ I1 + O(T )
5π h̄ 2 µ
 !2 
3N 5π 2 kB T
= F 1 + + O(T 4 ) (3.63)
5 V 12 F
Hence, the specific heat of a gas of free fermions is:

π2 kB T
CV = N kB (3.64)
2 F
Note that this can be written in the more general form:
CV = (const.) · kB · g (F ) kB T (3.65)
The number of electrons which are thermally excited above the ground state is ∼
g (F ) kB T ; each such electron contributes energy ∼ kB T and, hence, gives a specific
heat contribution of kB . Electrons give such a contribution to the specific heat of a
metal.
3.5 Ising Model, Mean Field Theory, Phases

Consider a model of spins on a lattice in a magnetic field:
Siz ≡ 2h Siz
X X
H = −gµB B (3.66)
i i
with Siz = ±1/2. The partition function for such a system is:
!N
h
Z = 2 cosh (3.67)
kB T
The average magnetization is:
1 h
Siz = tanh (3.68)
2 kB T
The susceptibility, χ, is defined by
!
∂ X z
χ= S (3.69)
∂h i i h=0
For free spins on a lattice,

1 1
χ= N (3.70)
2 kB T
A susceptibility which is inversely proportional to temperature is called a Curie suc-
septibility. In problem set 3, you will show that the susceptibility is much smaller for
a system of electrons.
Now consider a model of spins on a lattice such that each spin interacts with its
neighbors according to:
1X z z
H=− JS S (3.71)
2 hi,ji i j
This Hamiltonian has a symmetry
Siz → −Siz (3.72)
For kB T J, the interaction between the spins will not be important and the
susceptibility will be of the Curie form. For kB T < J, however, the behavior will be
much different. We can understand this qualitatively using mean field theory.
Let us approximate the interaction of each spin with its neighbors by an interaction
with a mean-field, h:
hSiz
X
H=− (3.73)
i
with h given by
JhSiz i = JzhSiz i
X
h= (3.74)
i
where z is the coordination number. In this field, the partition function is just
2 cosh kBhT and
h
hS z i = tanh (3.75)
kB T
Using the self-consistency condition, this is:
JzhS z i
hS z i = tanh (3.76)
kB T
For kB T < Jz, this has non-zero solutions, S z 6= 0 which break the symmetry
Siz → −Siz . In this phase, there is a spontaneous magnetization. For kB T > Jz,
there is only the solution S z = 0. In this phase the symmetry is unbroken and the
is no spontaneous magnetization. At kB T = Jz, there is a critical point at which a
phase transition occurs.
Chapter 4
Broken Translational Invariance in

the Solid State
4.1 Simple Energetics of Solids

Why do solids form? The Hamiltonian of the electrons and ions is:
X p2i X Pa2 X e2 X Z 2 e2 X Ze2
H= + + + − (4.1)
i 2me a 2M i>j |ri − rj | a>b |Ra − Rb | i,a |ri − Rb |
~a → R
It is invariant under the symmetry, ~ri → ~ri + ~a, R ~ a + ~a. However, the energy
can usually be minimized by forming a crystal. At low enough temperature, this will
win out over any possible entropy gain in a competing state, so crystallization will
occur. Why is the crystalline state energetically favorable? This depends on the type
of crystal. Different types of crystals minimize different terms in the Hamiltonian.
In molecular and ionic crystals, the potential energy is minimized. In a molecular
crystal, such as solid nitrogen, there is a van der Waals attraction between molecules
caused by polarization of one by the other. The van der Waals attraction is balanced
by short-range repulsion, thereby leading to a crystalline ground state. In an ionic
crystal, such as NaCl, the electrostatic energy of the ions is minimized (one must
be careful to take into account charge neutrality, without which the electrostatic
30
Chapter 4: Broken Translational Invariance in the Solid State 31
energy diverges). In covalent and metallic crystals, crystallization is driven by the

minimization of the electronic kinetic energy. In a metal, such as sodium, the kinetic
energy of the electrons is lowered by their ability to move throughout the metal. In
a covalent solid, such as diamond, the same is true. The kinetic energy gain is high
enough that such a bond can even occur between just two molecules (as in organic
chemistry). The energetic gain of a solid is called the cohesive energy.
4.2 Phonons: Linear Chain
4.3 Quantum Mechanics of a Linear Chain

As a toy model of a solid, let us consider a linear chain of masses m connected by
springs with spring constant B. Suppose that the equilibrium spacing between the
~ j,
masses is a. The equilibrium positions define a 1D lattice. The lattice ‘vectors’, R
are defined by:
~ j = ja
R (4.2)
~ and R
They connect the origin to all points of the lattice. If R ~ 0 are lattice vectors,
~ +R
then R ~ 0 are also lattice vectors. A set of basis vectors is a minimal set of vectors
which generate the full set of lattice vectors by taking linear combinations of the basis
vectors. In our 1D lattice, a is the basis vector.
Let ui be the displacement of the ith mass from its equilibrium position and let
pi be the corresponding momentum. Let us assume that there are N masses, and
let’s impose a periodic boundary condition, ui = ui+N . The Hamiltonian for such a
system is:
1 X 2 1 X
H= pi + B (ui − ui+1 )2 (4.3)
2m i 2 i
Let us use the Fourier transform representation:
1 X
uj = √ uk eikja
N k
1 X
pj = √ pk eikja (4.4)
N k
As a result of the periodic boundary condition, the allowed k’s are:
2πn
k= (4.5)
Na
We can invert (4.4):
1 X ik0 ja 1 XX 0
√ uj e = uk ei(k−k )ja
N j N k j
1 X ik0 ja
√ uj e = uk 0 (4.6)
N j
Note that u†k = u− k, p†k = p− k since u†j = uj , p†j = pj . They satisfy the commutation
relations:
1 X ikja ik0 j 0 a
[pk , uk0 ] = e e [pj , uj 0 ]
N j,j 0
1 X ikja ik0 j 0 a
= e e − ih̄δjj 0
N j,j 0
1 X i(k−k0 )ja
= −ih̄ e
N j
= −ih̄δkk0 (4.7)
Hence, pk and uk0 commute unless k = k 0 .

The displacements described by uk are the same as those described by uk+ 2πm for
a
any integer n:
1 X 2πm
uk+ 2πm = √ uj e−i(k+ a )ja
a
N j
1 X
= √ uj e−ikja
N j
= uk (4.8)
Hence,
2πm
k≡k+ (4.9)
a
Hence, we can restrict attention to n = 0, 1, . . . , N − 1 in (4.5). The Hamiltonian can
be rewritten: !2
X 1 ka
H= pk p−k + 2B sin uk u−k (4.10)
k 2m 2
Recalling the solution of the harmonic oscillator problem, we define:
 
 !2  14
1  ka  i 
ak = √  4mB sin

uk + p 

2h̄  2 2 14 k 
4mB sin ka
2
 
 !2  14
1  ka i 
a†k = √  4mB sin  u−k − p 
2 14 −k 
 (4.11)
2h̄  2
ka
4mB sin 2
(Recall that u†k = u−k , p†k = p−k .) which satisfy:
[ak , a†k ] = 1 (4.12)
Then:
1

a†k ak
X
H= h̄ωk + (4.13)
k 2
with
12
B ka

ωk = 2 sin (4.14)

m 2
Hence, the linear chain is equivalent to a system of N independent harmonic oscilla-
tors. Its thermodynamics can be described by the Planck distribution.
The operators a†k , ak are said to create and annihilate phonons. We say that a

E
state ψ with
E E
a†k ak ψ = Nk ψ (4.15)

has Nk phonons of momentum k. Phonons are the quanta of lattice vibrations,

analogous to photons, which are the quanta of oscillations of the electromagnetic
field.
Observe that, as k → 0, ωk → 0:
!1
Ba 2
ωk→0 = k
m/a
!1
Ba 2
= k (4.16)
ρ
The physical reason for this is simple: an oscillation with k = 0 is a uniform transla-
tion of the linear chain, which costs no energy.
Note that the reason for this is that the Hamiltonian is invariant under translations
ui → ui + λ. However, the ground state is not: the masses are located at the points
xj = ja. Translational invariance has been spontaneously broken. Of course, it could
just as well be broken with xj = ja + λ and, for this reason, ωk → 0 as k → 0. An
oscillatory mode of this type is called a Goldstone mode.
Consider now the case in which the masses are not equivalent. Let the masses
alternate between m and M . As we will see, the phonon spectra will be more com-
plicated. Let a be the distance between one m and the next m. The Hamiltonian
is:
X 1 1 2 1 1

H= p21,i + p2,i + B (u1,i − u2,i )2 + B (u2,i − u1,i+1 )2 (4.17)
i 2m 2M 2 2
The equations of motion are:
d2
m u1,i = −B [(u1,i − u2,i ) − (u2,i−1 − u1,i )]
dt2
d2
M 2 u2,i = −B [(u2,i − u1,i ) + (u2,i − u1,i+1 )] (4.18)
dt
Going again to the Fourier representation, α = 1, 2
1 X
uα,j = √ uα,k eikja (4.19)
N k
Where the allowed k’s are:
2π
k= n (4.20)
N
if there are 2N masses. As before,
uα,k = uα,k+ 2πn (4.21)

a
     
2
m dtd 2 0 u1,k  2B −B 1+e ika
u1,k 
= (4.22)
   
   
d2
0 M dt2 u2,k − B 1 + e−ika 2B u2,k
Fourier transforming in time:

     
2
 − mωk 0 u1,k   2B − B 1 + eika   u1,k 
=  (4.23)

  
0 − M ωk2 u2,k − B 1 + e−ika 2B u2,k
This eigenvalue equation has the solutions:

 v 
u 2 !2
 1 1 u 1 1 4 ka
2
ω± =B + ± + − sin (4.24)
t 
M m M m nM 2

Observe that !1
− Ba 2
ωk→0 = k (4.25)
2(m + M )/a
This is the acoustic branch of the phonon spectrum in which m and M move in phase.
As k → 0, this is a translation, so ωk− → 0. Acoustic phonons are responsible for
sound. Also note that
12
2B

−
ωk=π/a = (4.26)
M
Meanwhile, ω + is the optical branch of the spectrum (these phonons scatter light), in
which m and M move in opposite directions.
s
1 1

+
ωk→0 = 2B + (4.27)
m M
12
2B

+
ωk=π/a = (4.28)
m
so there is a gap in the spectrum of width
12 12
2B 2B

ωgap = − (4.29)
m M
Note that if we take m = M , we recover the previous results, with a → a/2.

This is an example of what is called a lattice with a basis. Not every site on the
chain is equivalent. We can think of the chain of 2N masses as a lattice with N sites.
Each lattice site has a two-site basis: one of these sites has a mass m and the other
has a mass M . Sodium Chloride is a simple subic lattice with a two-site basis: the
sodium ions are at the vertices of an FCC lattice and the chlorine ions are displaced
from them.
4.3.1 Statistical Mechnics of a Linear Chain
Let us return to the case of a linear chain of masses m separated by springs of force
constant B, at equilibrium distance a. The excitations of this system are phonons
2πm
which can have momenta k ∈ [−π/a, π/a] (since k ≡ k + a
), corresponding to
energies
12
B ka

h̄ωk = 2 sin

(4.30)
m 2
Phonons are bosons whose number is not conserved, so they obey the Planck distri-
bution. Hence, the energy of a linear chain at finite temperature is given by:
X h̄ωk
E =
k eβh̄ωk
−1
Z π
a dk h̄ωk
= L −π βh̄ω
(4.31)
a
(2π) e k − 1
Changing variables from k to ω,

√
L Z 4B/m 2 dω h̄ω
E = 2· q
βh̄ω
2π 0 a 4B − ω 2 e −1
√ m
2N Z 4B/m dω h̄ω
= q
βh̄ω
π 0 4B
−ω e
2 −1
m√
2N (kB T )2 βh̄ 4B/m
Z
1 x dx
= r (4.32)
π h̄ 0
2 ex − 1
4B x
m
− βh̄
q
In the limit kB T h̄ 4B/m, we can take the upper limit of integration to ∞ and
drop the x-dependent term in the square root in the denominator of the integrand:
N m (kB T )2 Z ∞ x dx
r
E =
π 4B h̄ 0 ex − 1
(4.33)
Hence, Cv ∼ T at low temperatures.

q
In the limit kB T h̄ 4B/m, we can approximate ex ≈ 1 + x:
√
2N (kB T )2 Z βh̄ 4B/m dx
E = r
π h̄ 0 4B

x
2
m
− βh̄
= N kB (4.34)
In the intermediate temperature regime, a more careful analysis is necessary. In

q q
4B
particular, note that the density of states, 1/ m
− ω 2 diverges at ω = 4B/m; this
is an example of a van Hove singularity. If we had alternating masses on springs,
then the expression for the energy would have two integrals, one over the acoustic
modes and one over the optical modes.
4.4 Lessons from the 1D chain

In the course of our analysis of the 1D chain, we developed the following strategy,
which we will apply to crystals more generally in subsequent sections.
• Expand the Hamiltonian to Quadratic Order
• Fourier transform the Hamiltonian into momentum space
• Identify the Brillouin zone (range of distinct ~ks)
• Rewrite the Hamiltonian in terms of creation and annihilation operators

• Obtain the Spectrum
• Compute the Density of States
• Use the Planck distribution to obtain the thermodynamics of the vibrational

modes of the crystal.
4.5 Discrete Translational Invariance: the Recip-

rocal Lattice, Brillouin Zones, Crystal Momen-
tum
2πn 2πn
Note that, in the above, momenta were only defined up to a
. The momenta a
form a lattice in k-space, called the reciprocal lattice. This is true of any function
which, like the ionic discplacements, is a function defined at the lattice sites. For

~ , defined on an arbitrary lattice, the Fourier transform
such a function, f R
~ ~

f˜ ~k = ~
eik·R f R
X
(4.35)
R
satisfies

f˜ ~k = f˜ ~k + G
~ (4.36)
~ is in the set of reciprocal lattice vectors, defined by:

if G
~ ~ ~
eiG·R = 1 , for all R (4.37)
The reciprocal lattice vectors also form a lattice since the sum of two reciprocal lattice
vectors is also a reciprocal lattice vector. This lattice is called the reciprocal lattice
or dual lattice.
In the analysis of the linear chain, we restricted momenta to |k| < π/a to avoid
double-counting of degrees of freedom. This restricted region in k-space is an example
of a Brillouin zone (or a first Brillouin zone). All of k-space can be obtained by
translating the Brillouin zone through all possible reciprocal lattice vectors. We
could have chosen our Brillouin zone differently by taking 0 < k < 2π/a. Physically,
there is no difference; the choice is a matter of convenience. What we need is a set of
points in k space such that no two of these points are connected by a reciprocal lattice
vector and such that all of k space can be obatained by translating the Brillouin zone
through all possible reciprocal lattice vectors. We could even choose a Brillouin zone
which is not connected, e.g. 0 < k < π/a. or 3π/a < k < 4π/a.
Later, we will consider solids with a more complicated lattice structure than our
linear chain. Once again, phonon spectra will be defined in the Brillouin zone. Since

f˜ ~k = f˜ ~k + G
~ , the phonon modes outside of the Brillouin zone are not physically
distinct from those inside. One way of defining the Brillouin zone for an arbitrary
lattice is to take all points in k space which are closer to the origin than to any
other point of the reciprocal lattice. Such a choice of Brillouin zone is also called the
Wigner-Seitz cell of the reciprocal lattice. We will discuss this in some detail later
but, for now, let us consider the case of a cubic lattice. The lattice vectors of a cubic
lattice of side a are:
~ n1 ,n2 ,n3 = a (n1 x̂ + n2 ŷ + n3 ẑ)
R (4.38)
The reciprocal lattice vectors are:
~ m1 ,m2 ,m3 = 2π (m1 x̂ + m2 ŷ + m3 ẑ)

G (4.39)
a
The reciprocal lattice vectors also form a cubic lattice. The first Brillouin zone
2π
(Wigner-Seitz cell of the reciprocal lattice) is given by the cube of side a
centered
at the origin. The volume of this cube is related to the density according to:
Z
d3 k 1 Nions
3
= 3 = (4.40)
B.Z. (2π) a V
As we have noted before, the ground state (and the Hamiltonian) of a crystal is in-
variant under the discrete group of translations through all lattice vectors. Whereas
full translational invariance leads to momentum conservation, lattice translational

symmetry leads to the conservation of crystal momentum – momentum up to a re-
ciprocal lattice vector. (See Ashcroft and Mermin, appendix M.) For instance, in a
collision between phonons, the difference between the incoming and outgoing phonon
momenta can be any reciprocal lattice vector, G. Physically, one may think of the
missing momentum as being taken by the lattice as a whole. This concept will also
be important when we condsider the problem of electrons moving in the background
of a lattice of ions.
4.6 Phonons: Continuum Elastic Theory

Consider the lattice of ions in a solid. Suppose the equilibrium positions of the ions are
~ i . Let us describe small displacements from these sites by a displacement
the sites R
~ i ). We will imagine that the crystal is just a system of masses connected by
field ~u(R
springs of equilibrium length a.
Before considering the details of the possible lattice structures of 2D and 3D
crystals, let us consider the properties of a crystal at length scales which are much
larger than the lattice spacing; this regime should be insensitive to many details of
the lattice. At length scales much longer than its lattice spacing, a crystalline solid
~ i ) by ~u(~r) (i.e. we replace the
can be modelled as an elastic medium. We replace ~u(R
~ i , by a continuous variable, ~r). Such an approximation is valid at
lattice vectors, R
length scales much larger than the lattice spacing, a, or, equivalently, at wavevectors
q 2π/a.
In 1D, we can take the continuum limit of our model of masses and springs:
!2
1 X dui 1 X
H = m + B (ui − ui+1 )2
2 i dt 2 i
!2
1 m X dui 1 X ui − ui+1 2
= a + Ba a
2 a i dt 2 i a
 !2 !2 
Z
1 du 1 du
→ dx  ρ + B  (4.41)
2 dt 2 dx
where ρ is the mass density and B is the bulk modulus. The equation of motion,
d2 u d2 u
= B (4.42)
dt2 dx2
has solutions
uk ei(kx−ωt)
X
u(x, t) = (4.43)
k
with v
u
uB
ω=t k (4.44)
ρ
The generalization to a 3D continuum elastic medium is:

~ ∇
ρ∂t2~u = (µ + λ) ∇ ~ · ~u + µ∇2~u (4.45)
where ρ is the mass density of the solid and µ and λ are the Lamé coefficients. Under
a dilatation, ~u(~r) = α~r, the change in the energy density of the elastic medium is
α2 (λ + 2µ/3)/2; under a shear stress, ux = αy, uy = uz = 0, it is α2 µ/2. In a crystal
– which has only a discrete rotational symmetry – there may be more parameters
than just µ and λ, depending on the symmetry of the lattice. In a crystal with cubic
symmetry, for instance, there are, in general, three independent parameters. We will
make life simple, however, and make the approximation of full rotational invariance.
The solutions are,
~
~u(~r, t) = ~ ei(k·~r−ωt) (4.46)
where is ~ a unit polarization vector, satisfy

−ρω 2~ = − (µ + λ) ~k ~k · ~ − µk 2~ (4.47)
For longitudinally polarized waves, ~k = k~,

s
2µ + λ
ωkl =± k ≡ ±vl k (4.48)
ρ
while transverse waves, ~k · ~ = 0 have

s
µ
ωkt =± k ≡ ±vs k (4.49)
ρ
Above, we introduced the concept of the polarization of a phonon. In 3D, the
displacements of the ions can be in any direction. The two directions perpendicular
to ~k are called transverse. Displacements along ~k are called longitudinal.
The Hamiltonian of this system,
1 ˙ 2 1 ~ · ~u 2 − 1 µ ~u · ∇2~u
Z
H= d3 r ρ ~u + (µ + λ) ∇ (4.50)
2 2 2
can be rewritten in terms of creation and annihilation operators,
" s #
1 √ ρ
ak,s = √ ρωk,s ~s · ~uk + i ~s · ~u˙ k
2h̄ " ωk,s
s #
1 √ ρ
a†k,s = √ ρωk,s ~s · ~u−k − i ~s · ~u˙ −k (4.51)
2h̄ ωk,s
as
1

a†k,s ak,s
X
H= h̄ωk,s + (4.52)
k,s 2
Inverting the above definitions, we can express the displacement ~u(r) in terms of
the creation and annihilation operators:
s
h̄ ~
~ a~k,s + a†−~k,s eik·~r
X
~u(~r) = s s
(4.53)
k,s 2ρV ωk
s = 1, 2, 3 corresponds to the longitudinal and two transverse polarizations. Acting
with ~uk either annihilates a phonon of momentum k or creates a phonon of momentum
−k.
The allowed ~k values are determined by the boundary conditions in a finite system.
For periodic boundary conditions in a cubic system of size V = L3 , the allowed ~k’s
are 2π
L
(n1 , n2 , n3 ). Hence, the ~k-space volume per allowed ~k is (2π)3 /V . Hence, we
can take the infinite-volume limit by making the replacement:
1 X ~
f (~k) = f (k) (∆~k)3
X
k
~ 3
(∆k) k
V Z 3~ ~
= d k f (k) (4.54)
(2π)3
4.7 Debye theory

Since a solid can be modelled as a collection of independent oscillators, we can obtain
the energy in thermal equilibrium using the Planck distribution function:
XZ d3 k h̄ωs (k)
E=V
s B.Z. (2π)3 eβh̄ωs (k) − 1
(4.55)
where s = 1, 2, 3 are the three polarizations of the phonons and the integral is over
the Brillouin zone.
This can be rewritten in terms of the phonon density of states, g(ω) as:
Z ∞ h̄ω
E=V dω g(ω)
0 eβh̄ω − 1
(4.56)
where
XZ d3 k
g(ω) = δ (ω − ωs (k)) (4.57)
s B.Z. (2π)3
The total number of states is given by:

Z ∞ Z ∞ XZ d3 k
dω g(ω) = dω δ (ω − ωs (k))
0 0 s B.Z. (2π)3
d3 k Z
= 3
B.Z. (2π)3
Nions
= 3 (4.58)
V
The total number of normal modes is equal to the total number of ion degrees of
freedom.
For a continuum elastic medium, there are two transverse modes with velocity vt
and one longitudinal mode with velocity vl . In the limit that the lattice spacing is
very small, a → 0, we expect this theory to be valid. In this limit, the Brillouin zone
is all of momentum space, so
Z
d3 k
gCEM (ω) = (2δ (ω − vt k)) + δ (ω − vl k)))
(2π)3 !
1 2 1
= + ω2 (4.59)
2π 2 vt3 vl3
In a crystalline solid, this will be a reasonable approximation to g(ω) for kB T
h̄vt /a where the only phonons present will be at low energies, far from the Brillouin
zone boundary. At high temperatures, there will be thermally excited phonons near
the Brillouin zone boundary, where the spectrum is definitely not linear, so we cannot
use the continuum approximation. In particular, this g(ω) does not have a finite
integral, which violates the condition that the integral should be the total number of
degrees of freedom.
A simple approximation, due to Debye, is to replace the Brillouin zone by a sphere
of radius kD and assume that the spectrum is linear up to kD . In other words, Debye
assumed that: 
3


2π 2 v 3
ω2 if ω < ωD
gD (ω) = 
 0 if ω > ωD
Here, we have assumed, for simplicity, that vl = vt and we have written ωD = vkD .
ωD is chosen so that
Z ∞
Nions
3 = dω g(ω)
V Z0ωD
3
= dω 2 3 ω 2
0 2π v
3
ωD
= (4.60)
2π 2 v 3
i.e.
1
ωD = 6π 2 v 3 Nions /V 3
(4.61)
With this choice,

Z ωD 3 h̄ω
E = V dω 2 3
ω 2 βh̄ω
0 2π v e −1
3(kB T )4 Z βh̄ωD x3
= V dx (4.62)
2π 2 v 3 h̄3 0 ex − 1
In the low temperature limit, βh̄ωD → ∞, we can take the upper limit of the
integral to ∞ and:
3(kB T )4 Z ∞ x3
E≈V dx (4.63)
2π 2 v 3 h̄3 0 ex − 1
The specific heat is:
4
x3
" Z ∞ #
3 12kB
CV ≈ T V dx x (4.64)
2π 2 v 3 h̄3 0 e −1
The T 3 contribution to the specific heat of a solid is often the most important con-
tribution to the measured specific heat.
For T → ∞,
3(kB T )4 Z βh̄ωD
E ≈ V dx x2
2π 2 v 3 h̄3 0
3
ωD
= V kB T
2π 2 v 3
= 3Nions kB T (4.65)
so
CV ≈ 3Nions kB (4.66)
The high-temperature specific heat is just kB /2 times the number of degrees of free-
dom, as in classical statistical mechanics.
At high-temperature, we were guaranteed the right result since the density of
states was normalized to give the correct total number of degrees of freedom. At
low-temperature, we obtain a qualitatively correct result since the spectrum is linear.
To obtain the exact result, we need to allow for longitudinal and transverse velocities

which depend on the direction, vt k̂ , vl k̂ , since rotational invariance is not present.
Debye’s formula interpolates between these well-understood limits.
We can define θD by kB θD = h̄ωD . For lead, which is soft, θD ≈ 88K, while for
diamond, which is hard, θD ≈ 1280K.
4.8 More Realistic Phonon Spectra: Optical Phonons,

van Hove Singularities
Although Debye’s theory is reasonable, it clearly oversimplifies certain aspects of the
physics. For instance, consider a crystal with a two-site basis. Half of the phonon
modes will be optical modes. A crude approximation for the optical modes is an
Einstein spectrum:
Nions
gE (ω) = δ(ω − ωE ) (4.67)
2
In such a case, the energy will be:
3(kB T )4 Z βh̄ωmax x3 Nions h̄ωE

E=V 3 dx x
+V βh̄ω
(4.68)
2 3
2π v h̄ 0 e −1 2 e E −1
with ωmax chosen so that
Nions ω3
3 = max (4.69)
2 2π 2 v 3
Another feature missed by Debye’s approximation is the existence of singularities

in the phonon density of states. Consider the spectrum of the linear chain:
12
B ka

ω(k) = 2 sin (4.70)

m 2
The minimum of this spectrum is at k = 0. Here, the density of states is well described
by Debye theory which, for a 1D chain predicts g(ω) ∼ const.. The maximum is at
k = π/a. Near the maximum, Debye theory breaks down; the density of states is
singular:
2 1
g(ω) = q (4.71)
πa ωmax
2 − ω2
In 3D, the singularity will be milder, but still present. Consider a cubic lattice.
The spectrum can be expanded about a maximum as:
2
ω(k) = ωmax − αx (kxmax − kx )2 − αy kymax − ky − αz (kzmax − kz )2 (4.72)
Then (6 maxima; 1/2 of each ellipsoid is in the B.Z.)

Z ωmax
G(ω) ≡ dω g(ω)
ω
1 V

= 6· · vol. of ellipsoid
2 (2π)3
3
V 4 (ωmax − ω) 2
= 3 π (4.73)
(2π)3 3 (αx αy αz ) 12
Differentiating:
1
3V (ωmax − ω) 2
g(ω) = 2 (4.74)
4π (αx αy αz ) 12
~ k ω(k) = 0, but the
In 2D and 3D, there can also be saddle points, where ∇
eigenvalues of the second derivative matrix have different signs. At a saddle point,
the phonon spectrum again has a square root singularity. van Hove proved that every
3D phonon spectrum has at least one maximum and two saddle points (one with one
negative eigenvalue, one with two negative eigenvalues). To see why this might be
true, draw the spectrum in the full k-space, repeating the Brillouin zone. Imagine
drawing lines connecting the minima of the spectrum to the nearest neighboring
minima (i.e. from each copy of the B.Z. to its neighbors). Imagine doing the same
with the maxima. These lines intersect; at these intersections, we expect saddle
points.
4.9 Lattice Structures

Thus far, we have focussed on general properties of the vibrational physics of crys-
talline solids. Real crystals come in a variety of different lattice structures, to which
we now turn our attention.
4.9.1 Bravais Lattices
Bravais lattices are the underlying structure of a crystal. A 3D Bravais lattice is

~
defined by the set of vectors R
n o
~ R
R ~ = n1~a1 + n2~a2 + n3~a3 ; ni ∈ Z (4.75)
where the vectors ~ai are the basis vectors of the Bravais lattice. (Do not confuse with
a lattice with a basis.) Every point of a Bravais lattice is equivalent to every other
point. In an elemental crystal, it is possible that the elemental ions are located at the
vertices of a Bravais lattice. In general, a crystal structure will be a Bravais lattice
with a basis.
The symmetry group of a Bravais lattice is the group of translations through the
lattice vectors together with some discrete rotation group about (any) one of the
lattice points. In the problem set (Ashcroft and Mermin, problem 7.6) you will show
that this rotation group can only have 2-fold, 3-fold, 4-fold, and 6-fold rotation axes.
There are 5 different types of Bravais lattice in 2D: square, rectangular, hexago-
nal, oblique, and body-centered rectangular. There are 14 different types of Bravais
lattices in 3D. The 3D Bravais lattices are discussed in are described in Ashcroft and
Mermin, chapter 7 (pp. 115-119). We will content ourselves with listing the Bravais
lattices and discussing some important examples.
Bravais lattices can be grouped according to their symmetries. All but one can
be obtained by deforming the cubic lattices to lower the symmetry.
• Cubic symmetry: cubic, FCC, BCC
• Tetragonal: stretched in one direction, a × a × c; tetragonal, centered tetragonal
• Orthorhombic: sides of 3 different lengths a × b × c, at right angles to each

other; orthorhombic, base-centered, face-centered, body-centered.
• Monoclinic: One face is a parallelogram, the other two are rectangular; mono-
clinic, centered monoclinic.
• Triclinic: All faces are parallelograms.
• Trigonal: Each face is an a × a rhombus.
• Hexagonal: 2D hexagonal lattices of side a, stacked directly above one another,

with spacing c.
Examples:
• Simple cubic lattice: ~ai = a x̂i .
• Body-centered cubic (BCC) lattice: points of a cubic lattice, together with the
centers of the cubes ∼
= interpenetrating cubic lattices offset by 1/2 the body-
diagonal.
a
~a1 = a x̂1 , ~a2 = a x̂2 , ~a3 = (x̂1 + x̂2 + x̂3 ) (4.76)
2
Examples: Ba, Li, Na, Fe, K, Tl
• Face-centered cubic (FCC) lattice: points of a cubic lattice, together with the
centers of the sides of the cubes, ∼
= interpenetrating cubic lattices offset by 1/2
a face-diagonal.
a a a
~a1 = (x̂2 + x̂3 ) , ~a2 = (x̂1 + x̂3 ) , ~a3 = (x̂1 + x̂2 ) (4.77)
2 2 2
Examples: Al, Au, Cu, Pb, Pt, Ca, Ce, Ar.
• Hexagonal Lattice: Parallel planes of triangular lattices.
a √
~a1 = a x̂1 , ~a2 = x̂1 + 3 x̂2 , ~a3 = c x̂3 (4.78)
2
Bravais lattices can be broken up into unit cells such that all of space can be
recovered by translating a unit cell through all possible lattice vectors. A primitive
unit cell is a unit cell of minimal volume. There are many possible choices of primitive
unit cells. Given a basis, ~a1 , ~a2 , ~a3 , a simple choice of unit cell is the region:
n o
~r ~r = x1~a1 + x2~a2 + x3~a3 ; xi ∈ [0, 1] (4.79)

The volume of this primitive unit cell and, thus, any primitive unit cell is:
~a1 · ~a2 × ~a3 (4.80)
An alternate, symmetrical choice is the Wigner-Seitz cell: the set of all points
which are closer to the origin than to any other point of the lattice. Examples:
Wigner-Seitz for square=square, hexagonal=hexagon (not parallelogram), oblique=distorted
hexagon, BCC=octohedron with each vertex cut off to give an extra square face (A+M
p.74).
4.9.2 Reciprocal Lattices
If ~a1 , ~a2 , ~a3 span a Bravais lattice, then
~b1 = 2π ~a2 × ~a3

~a1 · ~a2 × ~a3
~b2 ~a3 × ~a1
= 2π
~a1 · ~a2 × ~a3
~b3 ~a1 × ~a2
= 2π (4.81)
~a1 · ~a2 × ~a3
span the reciprocal lattice, which is also a bravais lattice.

~
The reciprocal of the reciprocal lattice is the set of all vectors ~r satisfying eiG·~r = 1
~ i.e. it is the original lattice.
for any recprocal lattice vector G,
As we discussed above, a simple cubic lattice spanned by
ax̂1 , ax̂2 , ax̂3 (4.82)

has the simple cubic reciprocal lattice spanned by:

2π 2π 2π
x̂1 , x̂2 , x̂3 (4.83)
a a a
An FCC lattice spanned by:
a a a
(x̂2 + x̂3 ) , (x̂1 + x̂3 ) , (x̂1 + x̂2 ) (4.84)
2 2 2
has a BCC reciprocal lattice spanned by:
4π 1 4π 1 4π 1
(x̂2 + x̂3 − x̂1 ) , (x̂1 + x̂3 − x̂2 ) , (x̂1 + x̂2 − x̂3 ) (4.85)
a 2 a 2 a 2
Conversely, a BCC lattice has an FCC reciprocal lattice.
The Wigner-Seitz primitive unit cell of the reciprocal lattice is the first Brillouin
zone. In the problem set (Ashcroft and Mermin, problem 5.1), you will show that the
Brillouin zone has volume (2π)3 /v if the volume of the unit cell of the original lattice
is v. The first Brillouin zone is enclosed in the planes which are the perpendicular
bisectors of the reciprocal lattice vectors. These planes are called Bragg planes for
reasons which will become clear below.
4.9.3 Bravais Lattices with a Basis
Most crystalline solids are not Bravais lattices: not every ionic site is equivalent to
every other. In a compound this is necessarily true; even in elemental crystals it is
often the case that there are inequivalent sites in the crystal structure. These crystal
structures are lattices with a basis. The classification of such structures is discused
in Ashcroft and Mermin, chapter 7 (pp. 119-126). Again, we will content ourselves
with discussing some important examples.
• Honeycomb Lattice (2D): A triangular lattice with a two-site basis. The trian-
gular lattice is spanned by:
a √ √
~a1 = 3 x̂1 + 3 x̂2 , ~a2 = a 3 x̂1 (4.86)
2
The two-site basis is:

0, a x̂2 (4.87)
Example: Graphite
• Diamond Lattice: FCC lattice with a two-site basis: The two-site basis is:
a
0, (x̂1 + x̂2 + x̂3 ) (4.88)
4
Example: Diamond, Si, Ge
• Hexagonal Close-Packed (HCP): Hexagonal lattice with a two-site basis:
a a c
0, x̂1 + √ x̂2 + x̂3 (4.89)
2 2 3 2
Examples: Be,Mg,Zn, . . . .
• Sodium Chloride: Cubic lattice with Na and Cl at alternate sites ∼

= FCC lattice
with a two-site basis:
a
0, (x̂1 + x̂2 + x̂3 ) (4.90)
2
Examples: NaCl,NaF,KCl
4.10 Bragg Scattering

One way of experimentally probing a condensed matter system involves scattering a
photon or neutron off the system and studying the energy and angular dependence
of the resulting cross-section. Crystal structure experiments have usually been done
with X-rays.
Let us first examine this problem intuitively and then in a more systematic fashion.
Consider, first, elastic scattering of X-rays. Think of the X-rays as photons which
can take different paths through the crystal. Consider the case in which ~k is the
wavevector of an incoming photon and ~k 0 is the wavevector of an outgoing photon.

Let us, furthermore, assume that the photon only scatters off one of the atoms in
the crystal (the probability of multiple scattering is very low). This atom can be
any one of the atoms in the crystal. These different scattring events will interfere
constructively if the path lengths differ by an integer number of wavelengths. The
~ as compared to an atom at the
extra path length for a scattering off an atom at R,
origin is:

~ · ~k + −R
R ~ · ~k 0 = R
~ · ~k − ~k 0 (4.91)
~ then scattering interferes

If this is an integral multiple of 2π for all lattice vectors R,
constructively. By definition, this implies that ~k − ~k 0 must be a reciprocal lattice
vector. For elastic scatering, |~k| = |~k 0 |, so this implies that there is a reciprocal lattic
~ of magnitude
vector G

G = 2|~
~
k| sin θ (4.92)
where θ is the angle between the incoming and outgoing X-rays.

To rederive this result more formally, let us assume that our crystal is in thermal
equilibrium at inverse temperature β and that photons interact with our crystal via
the Hamiltonian H 0 . Suppose that photons of momentum ~ki , and energy ωi are
scattered by our system. The differential cross-section for the photons to be scattered
into a solid angle dΩ centered about ~kf and into the energy range between ωf ± dω
is:
d2 σ X kf D E2
= ~kf ; m |H 0 | ~ki ; n e−βEn δ (ω + En − Em ) (4.93)

dΩ dω k
m,n i
where ω = ωi − ωf and n and m label the initial and final states of our crystal. Let
~q = ~kf − ~ki .
Let us assume that the interactions between the photon and the ions in our system
is of the form:

H0 = ~ + ~u R
~
X
U ~x − R (4.94)
R
Then
* +
D
~kf ; m |H 0 | ~ki ; n
E Z
1 X
~ + ~u R~
= d3~x ei~q·~x m U ~x − R n

V
R

* +
1X ~ ~
e−i~q·(R+~u(R)) n Ũ (~q)

= m
V
R
1 X −i~q·R~

~
m e−i~q·~u(R) n Ũ (~q)

= e
V
" R #
1 X −i~q·R~ D −i~q·~u(0) E
= e m e n Ũ (~
q) (4.95)
V R
Let us consider, first, the case of elastic scattering, in which the state of the crystal
E
does not change. Then |ni = |m , |ki | = |kf | ≡ k, |q| = 2k sin 2θ , and:
" #
D
~kf ; n |H | ~ki ; n
0
E 1 X −i~q·R~ D −i~q·~u(0) E
= e n e n Ũ (~
q) (4.96)
V R
Let us focus on the sum over the lattice:
1 X −i~q·R~ X
e = δq~,G~ (4.97)
V R R.L.V. G
The sum is 1 if ~q is a reciprocal lattice vector and vanishes otherwise. The scattering
cross-section is given by:
d2 σ
" #
D E2 2
e−βEn n e−i~q·~u(0) n Ũ (~q)
X X
= δq~,G~ (4.98)
dΩdω n R.L.V. G
In other words, the scattering cross-section is peaked when the photon is scattered
through a reciprocal lattice vector ~kf = ~ki + G.
~ For elastic scattering, this requires
2 2
~ki = ~ki + G
~ (4.99)
or,
2
~
G = −2 ~ki · G
~ (4.100)
This is called the Bragg condition. It is satisfied when the endpoint of ~k is on a Bragg
plane. When it is satisfied, Bragg scattering occurs.
When there is structure within the unit cell, as in a lattice with a basis, the
formula is slightly more complicated. We can replace the photon-ion interaction by:

H0 = ~ + ~b + ~u(R
~ + ~b)
XX
Ub ~x − R (4.101)
R b
Then,
1 X −i~q·R~
e Ũ (~q) (4.102)
V R
is replaced by
1 X −i~q·R~
fq e (4.103)
V R
where
Ub (~x) e−i~q·~x
X
fq = (4.104)
b
As a result of the structure factor, fq , the scattering amplitude need not have a peak
at every reciprocal lattic vector, ~q.
Of course, the probability that the detector is set up at precisely the right angle
to receive ~kf = ~ki + G
~ is very low. Hence, these experiments are usually done with
a powder so that there will be Bragg scattering whenever 2k sin 2θ = |G|. By varying
θ, a series of peaks are seen at, e.g. π/6, π/4, etc., from which the reciprocal lattice
vectors are reconstructed.
−1
Since |k| ∼ |G| ∼ 1Å , the energy of the incoming photons is ∼ h̄ck ∼ 104 eV
which is definitely in the X-ray range.
Thus far, we have not looked closely at the factor:
D E2
m e−i~q·~u(0) n (4.105)

This factor results from the vibration of the lattice due to phonons. In elastic scat-
tering, the amplitude of the peak will be reduced by this factor since the probability
of the ions forming a perfect lattice is less than 1. The inelastic amplitude will con-
tain contributions from processes in which the incoming photon or neutron creates a
phonon, thereby losing some energy. By measuring inelastic neutron scattering (for
which the energy resolution is better than for X-rays), we can learn a great deal about
the phonon spectrum.
Chapter 5
Electronic Bands
5.1 Introduction
Thus far, we have ignored the dynamics of the elctrons and focussed on the ionic
vibrations. However, the electrons are important for many properties of solids. In
metals, the specific heat is actually CV = γT + αT 3 . The linear term is due to the
electrons. Electrical conduction is almost always due to the electrons, so we will need
to understand the dynamics of electrons in solids in order to compute, for instance,
the conductivity σ(T, ω).
In order to do this, we will need to understand the quantum mechanics of electrons
in the periodic potential due to the ions. Such an analysis will enable us to understand
some broad features of the electronic properties of crystalline solids, such as the
distinction between metals and insulators.
5.2 Independent Electrons in a Periodic Potential:

Bloch’s theorem
Let us first neglect all interactions between the electrons and focus on the interactions
between each electron and the ions. This may seem crazy since the inter-electron
57
Chapter 5: Electronic Bands 58
interaction isn’t small, but let us make this approximation and proceed. At some
level, we can say that we will include the electronic contribution to the potential
in some average sense so that the electrons move in the potential created by the
ions and by the average electron density (of course, we shoudld actually do this self-
consistently). Later, we will see why this is sensible.
When the electrons do not interact with each other, the many-electron wavefunc-
tion can be constructed as a Slater determinant of single-electron wavefunctions.
Hence, we have reduced the problem to that of a single electron moving in a lattice
of ions. The Hamiltonian for such a problem is:
h̄2 2 X ~ − ~u(R))
~
H=− ∇ + V (~x − R (5.1)
2m R
~
expanding in powers of ~u(R),
h̄2 2 X ~ − ∇V
X
~ (~x − R)
~ · ~u(R)
~ + ...
H=− ∇ + V (~x − R) (5.2)
2m R R
The third term and the . . . are electron-phonon interaction terms. They can be treated
as a perturbation. We will focus on the first two terms, which describe an electron
moving in a periodic potential. This highly simplified problem already contains much
of the qualitative physics of a solid.
Let us begin by proving an important theorem due to Bloch.
~ = V (~r) for all lattice vectors R
Bloch’s Theorem: If V (~r + R) ~ of some given
lattice, then for any solution of the Schrödinger equation in this potential,
h̄2 2
− ∇ ψ(~r) + V (~r)ψ(~r) = Eψ(~r) (5.3)
2m
there exists a ~k such that
~ = ei~k·R~ ψ(~r)
ψ(~r + R) (5.4)
Proof: Consider the lattice translation operator TR which acts according to
~
TR χ(~r) = χ(~r + R) (5.5)
Then
~
TR Hχ(~r) = HTR χ(~r + R) (5.6)
i.e. [TR , H] = 0. Hence, we can take our energy eigenstates to be eigenstates of TR .

Hence, for any energy eigenstate ψ(~r),
~
TR ψ(~r) = c(R)ψ(~
r) (5.7)
The additivity of the translation group implies that,
~ R
c(R)c( ~ 0 ) = c(R
~ +R
~ 0) (5.8)
Hence, there is some k such that

~ ~
~ = eik·R
c(R) (5.9)
~ ~ ~ is a reciprocal lattice vector, we can always take ~k to be in the
Since eiG·R = 1 if G
first Brillouin zone.
5.3 Tight-Binding Models

Let’s consider a very simple model of a 1D solid in which we imagine that the atomic
nuclei lie along a chain of spacing a. Consider a single ion and focus on two of its
electronic energy levels. In real systems, we will probably consider s and d orbitals,
but this is not important here; in our toy model, these are simply two electronic
states which are localized about the atomic nucleus. We’ll call them |1i and |2i, with
energies 01 and 02 . Let’s further imagine that the splitting 02 − 01 between these levels
is large. Now, when we put this atom in the linear chain, there will be some overlap
between these levels and the corresponding energy levels on neighboring atoms. We
can model such a system by the Hamiltonian:
X
10 |R, 1ihR, 1| + 20 |R, 2ihR, 2| − (t1 |R, 1ihR0 , 1| + t2 |R, 2ihR0 , 2|)
X
H=
R R,R0 n.n
(5.10)
We have assumed that t1 is the amplitude for an electron at |R, 1i to hop to |R0 , 1i,
and similarly for t2 . For simplicity, we have ignored the possibility of hopping from
|R, 1i to |R0 , 2i, which is unimportant anyway when 0 is large. The eigenstates of
this Hamiltonian are:
eikR |R, 1i
X
|k, 1i = (5.11)
R
with energy
1 (k) = 01 − 2t1 cos ka (5.12)
and
eikR |R, 2i
X
|k, 2i = (5.13)
R
with energy
2 (k) = 02 − 2t2 cos ka (5.14)
Note, first, that k lives in the first Brillouin zone since

2πn E
|k, ii ≡ k + ,i (5.15)

a
Now, observe that the two atomic energy levels have broadened into two energy
bands. There is a band gap between these bands of magnitude 02 − 01 − 2t1 − 2t2 . This
is a characteristic feature of electronic states in a periodic potential: the states break
up into bands with energy gaps separating the bands.
How many states are there in each band? As we discussed in the context of
phonons, there are as many allowed k’s in the Brillouin zone as there are ions in the
crystal. Let’s repeat the argument. The Brillouin zone has k-space extent 2π/a. In
a finite-size system of length L with periodic boundary conditions, allowed k’s are
of the form 2πn/L where n is an integer. Hence, there are L/a = Nions allowed k’s
in the Brillouin zone. (This argument generalizes to arbitrary lattices in arbitrary
dimension.) Hence, there are as many states as lattice sites. Each state can be filled
by one up-spin electon and one down-spin electron. Hence, if the atom is monovalent
– i.e. if there is one electron per site – then, in the ground state, the lower band,
|k, 1i is half-filled and the upper band is empty. The Fermi energy is at 01 . The
Fermi momentum (or, more properly, Fermi crystal momentum) is at ±π/2a. At
low temperature, the fact that there is a gap far away from the Fermi momentum is
unimportant, and the Fermi sea will behave just like the Fermi sea of a free Fermi
gas. In particular, there is no energy gap in the many-electron spectrum since we
can always excite an electron from a filled state just below the Fermi surface to one
of the unfilled states just above the gap. For instance, the electronic contribution to
the specific heat will be CV ∼ T . The difference is that the density of states will
be different from that of a free Fermi gas. In situations such as this, when a band
is partially filled, the crystal is (almost always) a metal. (Sometimes inter-electron
interactions can make such a system an insulator.)
If there are two electrons per lattice site, then the lower band is filled and the
upper band is empty in the ground state. In such a case, there is an energy gap
Eg = 02 − 01 − 2t1 − 2t2 between the ground state and the lowest excited state
which necessarily involves exciting an electron from the lower band to the upper
band. Crystals of this type, which have no partially filled bands, are insulators. The
electronic contribution to the specific heat will be suppressed by a factor of e−Eg /T .
Note that the above tight-binding model can be generalized to arbitrary dimension
of lattice. For instance, a cubic lattice with one orbital per site has tight-binding
spectrum:
(k) = −2t (cos kx a + cos ky a + cos kz a) (5.16)
Again, if there is one electron per site, the band will be half-filled (and metallic); if
there are two electrons per site the band will be filled (and insulating).
The model which we have just examined is grossly oversimplified but can, never-
theless, be justified, to an extent. Let us reconsider our lattice of atoms.
Electronic orbitals of an isolated atom:
ϕn (~r) (5.17)
with energies n :
h̄2 2
− ∇ ϕn (~r) + V (~r)ϕn (~r) = n ϕn (~r) (5.18)
2m
We now want to solve:
h̄2 2 X
~ ψk (~r) = (k)ψk (~r)
− ∇ ψk (~r) + V (~r + R) (5.19)
2m R
Let’s try the ansatz:

~ ~
~
cn eik·R ϕn (~r + R)
X
ψk (~r) = (5.20)
R,n
which satisfies Bloch’s theorem. Substituting into Schrödinger’s equation and taking
the matrix element with ϕm , we get:
~ 0 ) PR,n cn ei~k·R~ ϕn (~r + R)

2

R
d3 r ϕ∗m (~r) − 2m
h̄
∇2 +
P
R0 V (~r + R ~ =
~ ~ ~ (5.21)
(k) d3 r ϕ∗m (~r) eik·R ϕn (~r + R)
R P
R,n cn
Let’s write
h̄2 2 X ~ 0 h̄2 2 ~ +
X
~ 0)
− ∇ + V (~r + R ) = − ∇ + V (~r + R) V (~r + R
2m R0
2m R0 6=R
= Hat,R + ∆VR (r) (5.22)
Then, we have
Z Z
~ ~ ~ ~
d3 r ϕ∗m (~r) ~ = (k)
cn eik·R (Hat,R + ∆VR (r)) ϕn (~r + R) d3 r ϕ∗m (~r) ~
cn eik·R ϕn (~r + R)
X X
R,n R,n
(5.23)
Z Z
~ ~ ~ ~
cn n eik·R ~ +
d3 r ϕ∗m (~r)ϕn (~r + R) cn eik·R ~ =
d3 r ϕ∗m (~r)∆VR (r)ϕn (~r + R)
X X
R,n R,n Z
i~k·R
~ ~ (5.24)
d3 r ϕ∗m (~r)ϕn (~r + R)
X
(k)cm + (k) e cn
R6=0,n
i~k·R
~ ~ +
d3 r ϕ∗m (~r)ϕn (~r + R)
P R
c m m + R6=0,n cn n e
i~k·R
~ ~ =
d3 r ϕ∗m (~r) [∆VR (r)] ϕn (~r + R)
P R
R,n cn e
i~k·R
~ ~ (5.25)
d3 r ϕ∗m (~r)ϕn (~r + R)
P R
(k)cm + (k) R6=0,n e cn
Writing:
Z
~ =
αmn (R) ~
d3 r ϕ∗m (~r)ϕn (~r + R)
Z
~ = −
γmn (R) ~
d3 r ϕ∗m (~r) [∆VR (r)] ϕn (~r + R) (5.26)
We have:
~ ~ ~ ~
~ =
cn (n − (k)) eik·R αmn (R) ~
cn eik·R γmn (R)
X X
cm (m − (k)) + (5.27)
R6=0,n R,n
~ and γmn (R)

Both αmn (R) ~ are exponentially small, ∼ e−R/a0 . In particular, αmn (R)
~
~ are much larger for nearest neighbors than for any other sites, so let’s
and γmn (R)
~ n.n. ), γmn = γmn (R
neglect the other matrix elements and write αmn = αmn (R ~ n.n. ),
vmn = γmn (0). In problem 2 of problem set 7, so may make these approximations.
Suppose that we make the approximation that the lth orbital is well separated in
~ and γln (R)
energy from the others. Then we can neglect αln (R) ~ for n 6= l. We write
β = vll . Focusing on the m = l equation, we have:

~ ~ ~ ~
eik·R = β + γll eik·R
X X
(l − (k)) + αll (l − (k)) (5.28)
Rn.n. Rn.n.
Hence,
~ ~
β + γll Rn.n. eik·R
P
(k) = l − (5.29)
1 + αll Rn.n. ei~k·R~
P
If we neglect the α’s and retain only the γ’s, then we recover the result of our
phenomenological model. For instance, for the cubic lattice, we have:
(k) = [l − β] − 2γll [cos kx a + cos ky a + cos kz a] (5.30)
Tight-binding models give electronic wavefunctions which are a coherent super-

position of localized atomic orbitals. Such wavefunctions have very small amplitude
in the interstitial regions between the ions. Such models are valid, as we have seen,
when there is very little overlap between atomic wavefunctions on neighboring atoms.
In other words, a tight-binding model will be valid when the size of an atomic orbital
is smaller than the interatomic distance, i.e. a0 R. In the case of core electrons,
e.g. 1s, 2s, 2p, this is the case. However, this is often not the case for valence elec-
trons, e.g. 3s electrons. Nevertheless, the tight-binding method is a simple method
which gives many qualitative features of electronic bands. In the study of high-Tc
superconductivity, it has proven useful for this reason.
5.4 The δ-function Array

Let us now consider another simple toy-model of a solid, a 1D array of δ-functions:
∞
h̄2 d2
!
X
− + V δ(x − na) ψ(x) = E ψ(x) (5.31)
2m dx2 n=−∞
Between the peaks of the δ functions, ψ(x) must be a superposition of the plane waves
eiqx and e−iqx with energy E(q) = h̄2 q 2 /2m. Between x = 0 and x = a,
ψ(x) = eiqx+iα + e−iqx−iα (5.32)
with α complex. According to Bloch’s theorem,
ψ(x + a) = eika ψ(x) (5.33)
Hence, in the region between x = a and x = 2a,

ψ(x) = eika eiq(x−a)+iα + e−iq(x−a)−iα (5.34)
Note that k which determines the transformation property under a translation x →

x + a is not the same as q, which is the ‘local’ momentum of the electron, which
determines the energy. Continuity at x = a implies that
cos(qa + α) = eika cos α (5.35)

or,
cos qa − eika
tan α = (5.36)
sin qa
Integrating Schrödinger’s equation from x = a − tox = a + , we have,
2mV ika
sin(qa + α) − eika sin α = e cos α (5.37)
h̄2 q
or,
2mV
h̄2 q
eika − sin qa
tan α = (5.38)
cos qa − eika
Combining these equations,
!
2ika mV
e − 2 cos qa + 2 sin qa eika + 1 = 0 (5.39)
h̄ q
The sum of the two roots is cos ka:
cos(qa − δ)
cos ka = (5.40)
cos δ
where
mV
tan δ = (5.41)
h̄2 q
h i
For each k ∈ − πa , πa , there are infinitely many roots q of this equation, qn (k). The
energy spectrum of the nth band is:
h̄2
En (k) = [qn (k)]2 (5.42)
2m
±k have the same root qn (k) = qn (−k). Not all q’s are allowed. For instance, the
values qa − δ = nπ are not allowed. These regions are the energy gaps between bands.
Consider, for instance, k = π/a.
cos(qa − δ) = cos δ (5.43)
This has the solutions

qa = π , π + 2δ (5.44)
2mV a
For V small, the latter solution occurs at qa = π + πh̄2
. The energy gap is:
π π π 2mV a π

E2 − E1 ≈ E + 2 −E
a a a πh̄ a
≈ 2V /a (5.45)
5.5 Nearly Free Electron Approximation

According to Bloch’s theorem, electronic wavefunctions can be expanded as:
~ ~
ck−G ei(k−G)·~x
X
ψ(x) = (5.46)
G
In the nearly free electron approximation, we assume that electronic wavefunctions are
given by the superposition of a small number of plane waves. This approximation is
valid, for instance, when the periodic potential is weak and contains a limited number
of reciprocal lattice vectors.
Let’s see how this works. Schrödinger’s equation in momentum space reads:
h̄2 k 2
!
− (~k) ck + ck−G VG = 0
X
(5.47)
2m G
Second-order perturbation theory tells us that (let’s assume that Vk = 0)

|VG |2
(~k) = 0 (~k) +
X
(5.48)
~ ~ ~
G6=0 0 (k) − 0 (k − G)
where
h̄2 k 2
0 (~k) = (5.49)
2m
Perturbation theory will be valid so long as the second term is small, i.e. so long as
|VG | 0 (~k) − 0 (~k − G)

~ (5.50)
For generic ~k, this will be valid if VG is small. The correction to the energy will

be O |VG |2 . However, no matter how small VG is, perturbation theory fails for
degenerate states,
h̄2 k 2 h̄2 (k − G)2
= (5.51)
2m 2m
or, when the Bragg condition is satisfied,
G2 = 2~k · G
~ (5.52)
In other words, perturbation theory fails when ~k is near a Brillouin zone boundary.
Suppose that VG is very small so that we can neglect it away from the Brillouin
zone boundaries. Near a zone boundary, we can focus on the reciprocal lattice vector
~ and ignore VG0 for G
which it bisects, G ~ = ~ 0 . We keep only ck and ck−G , where
6 G
0 (k) ≈ 0 (k − G). We can thereby reduce Schrödinger’s equation to a 2 × 2 equation:

0 (~k) − (~k) ck + ck−G VG = 0

0 (~k − G)
~ − (~k − G)
~ ck−G + ck VG∗ = 0 (5.53)
~ =
VG0 for G ~ 0 can be handled by perturbation theory and, therefore, neglected in
6 G
the small VG0 limit. In this approximation, the eigenvalues are:
" r #
1 2
± (~k) = 0 (~k) + 0 (~k − G)
~ ± 0 (~k) − 0 (~k − G)
~ 2
+ 4|VG | (5.54)
2
At the zone boundary, the bands have been split by
+ (~k) − − (~k) = 2 |VG | (5.55)
~ =
The effects of VG0 for G ~ 0 are now handled perturbatively.
6 G
To summarize, the nearly free electron approximation gives energy bands which
are essentially free electron bands away from the Brillouin zone boundaries; near the
Brillouin zone boundaries, where the electronic crystal momenta satisfy the Bragg
condition, gaps are opened.
Though intuitively appealing, the nearly free electron approximation is not very
reasonable for real solids. Since
4πZe2
VG ≈ ∼ 13.6 eV (5.56)
G2
while F ∼ 10eV ,
|VG | ∼ 0 (~k) − 0 (~k − G)
~ (5.57)
and the nearly free electron approximation is not valid.
5.6 Some General Properties of Electronic Band

Structure
Much, much more can be said about electronic band structure. There are many
approximate methods of obtaining energy spectra for more realistic potentials. We
will content ourselves with two observations.
Band Overlap. In 2D and 3D, bands can overlap in energy. As a result, both the
first and second bands can be partially filled when there are two electrons per site.
Consider, for instance, a weak periodic potential of rectangular symmetry:
2π 2π
V (x, y) = Vx cos x + Vy cos x (5.58)
a1 a2
with Vx,y very small and a1 a2 . Using the nearly free electron approximation,
we have a spectrum which is essentially a free-electron parabola, with small gaps
opening at the zone boundary. Since the Brillouin zone is much shorter in the kx -
direction, the Fermi sea will cross the zone boundary in this direction, but not in
the ky -direction. Hence, there will be empty states in the first Brillouin zone, near
(0, ±πa2 ) and occupied states in the second Brillouin zone, near (±πa1 , 0). This is a
general feature of 2D and 3D bands. As a result, a solid can be metallic even when
it has two electrons per unit cell.
van Hove singularities. A second feature of electronic energy spectra is the ex-
istence of van Hove singularities. They are singularities in the electronic density of
states, g()
Z Z
d3 k
d g() f () = f ((k)) (5.59)
(2π)3
They occur for precisely the same reason as in the case of phonon spectra – as a result
of the lattice periodicity.
Consider the case of a tight-binding model on the square lattice with nearest-
neighbor hopping only.
(k) = −2t (cos kx a + cos ky a) (5.60)
~ k (k) = 2ta (sin kx a + sin ky a)

∇ (5.61)
The density of states is given by:

Z
d2 k
g() = 2 δ ( − (k)) (5.62)
(2π)2
Let’s change variables in the integral on the right to E and S which is the arc length
around an equal energy contour = (k):
1 Z dE
g() = dS δ ( − E)
2π 2
∇~ k (k)
1 Z 1
= dS (5.63)
2π 2
∇~ k (k)
The denominator on the right-hand-side vanishes at the minimum of the band, ~k =

(0, 0), the maxima ~k = (±π/a, ±π/a) and the saddle points ~k = (±π/a, 0), (0, ±π/a).
At the latter points, the density of states will have divergent slope.
5.7 The Fermi Surface

The Fermi surface is defined by
n (k) = µ (5.64)
By the Pauli principle, it is the surface in the Brillouin zones which separates the
occupied states, n (k) < µ, inside the Fermi surface from the unoccupied states
n (k) > µ outside the Fermi surface. All low-energy electronic excitations involve
holes just below the Fermi surface or electrons just above it. Metals have a Fermi
surface and, therefore, low-energy excitations. Insulators have no Fermi surface: µ

lies in a band gap, so there is no solution to (5.64).
In the low-density limit the Fermi surface is approximately circular (in 2D) or
spherical (in 3D). Consider the 2D tight-binding model
(k) = −2t (cos kx a + cos ky a) (5.65)
For ~k → 0,

(k) ≈ −4t + ta2 kx2 + kx2 (5.66)
Hence, for µ + 4t t, the Fermi surface is given by the circle:
µ + 4t
kx2 + kx2 = (5.67)
ta2
Similarly, in the nearly free electron approximation,
h̄2 k 2 X |VG |2
(~k) = + h̄2 (k−G)2
(5.68)
2m h̄2 k2
G6=0 2m − 2m
For µ → 0 and VG small, we can neglect the scond term, and, as in the free electron
case, the Fermi surface is given by
1q
k= 2mµ (5.69)
h̄
Away from the bottom of a band, however, the Fermi surface can look quite
different. In the tight-binding model, for instance, for µ = 0, the Fermi surface is the
diamond kx ± ky = ±π/a.
The chemical potential at zero temperature is usually called the Fermi energy, F .
The key measure of the number of low-lying states which are available to an electronic
system is the density of states at the Fermi energy, g(F ). When g(F ) is large, the
CV , σ, etc. are large; when g(F ) is small, these quantities are small.
5.8 Metals, Insulators, and Semiconductors

Earlier we saw that, in order to compute the vibrational properties of a solid, we
needed to determine the phonon spectra of the crystal. A characteristic feature of
these phonon spectra is that there is always an acoustic mode with ω(k) ∼ k for k
small. This mode is responsible for carrying sound in a solid, and it always gives a
CVph ∼ T 3 contribution to the specific heat.
In order to compute the electronic properties of a solid, we must similarly de-
termine the electronic spectra. If we ignore the interactions between electrons, the
electronic spectra are determined by the single-electron energy levels in the periodic
potential due to the ions. These energy spectra break up into bands. When there is a
partially filled band, there are low energy excitations, and the solid is a metal. There
will be a CVel ∼ T electronic contribution to the specific heat, as in a free fermion
gas. When all bands are either filled or completely empty, there is a gap between the
many-electron ground state and the first excited state; the solid is an insulator and
there is a negligible contribution to the low-temperature specific heat. Let us recall
how this works. Once we have determined the electronic band structure, n (k), we
can determine the electronic density-of-states:
XZ d2 k
g() = 2 δ ( − n (k)) (5.70)
n B.Z. (2π)2
With the density-of-states in hand, we can compute the thermodynamics. In the limit
kB T F ,
N Z ∞ 1
= d g()
V Z0µ eβ(−µ)+1
Z µ
1
Z ∞ 1
= d g() + d g() −1 + d g()
Z0F Z0 µ eZβ(−µ) +1 µ eβ(−µ) +1
µ 1 Z ∞ 1
= d g() + d g() − d g() + d g()
0 F
Z ∞0
+1 µe−β(−µ) +1 eβ(−µ)
N kB T dx
−βµ

≈ + (µ − F ) g(F ) + (g (µ + k B T x) − g (µ − k B T x)) + O e
V 0 ex + 1
N X∞
(kB T )n+1 (n) Z ∞ x2n−1
= + (µ − F ) g(F ) + g (µ) dx x
V n=1 n! 0 e +1
N
≈ + (µ − F ) g(F ) + (kB T )2 g 0 (F ) I1 (5.71)
V
with
Z ∞ xk
Ik = dx x (5.72)
0 e +1
We will only need
π2
I1 = (5.73)
6
Hence, to lowest order in T ,
(µ − F ) g(F ) ≈ −(kB T )2 g 0 (F ) I1 (5.74)
Meanwhile,
E Z ∞ 1
= d g()
V Z0F eβ(−µ) +1
Z µ Z µ
1
Z ∞ 1
= d g() + d g() + d g() −1 + d g()
0 F
Z µ 0
+ 1Z eβ(−µ)
µ eβ(−µ) +1
E0 1 ∞ 1
≈ + (µ − F ) F g(F ) − d g() −β(−µ) + d g() β(−µ)
V 0 e +1 µ e +1
Z ∞
E0 kB T dx
= + (µ − F ) F g(F ) + (µ + kB T x) g (µ + kB T x) −
V 0 ex + 1 !

(µ − kB T x) g (µ − kB T x) + O e−βµ
E0
≈ + (µ − F ) F g(F ) + (kB T )2 [g (F ) + F g 0 (F )] I1 (5.75)
V
Substituting (5.74) into the final line of (5.75), we have:
E E0
= + (kB T )2 g (F ) I1 (5.76)
V V
Hence, the electronic contribution to the low-temperature specific heat of a crys-

talline solid is:
CV π2 2
= k T g (F ) (5.77)
V 3 B
In a metal, the Fermi energy lies in some band; hence g (F ) is non-zero. In an
insulator, all bands are either completely full or completely empty. Hence, the Fermi
energy lies between two bands, and g (F ) = 0.
Each band contains twice as many single-electron levels (the factor of 2 comes
from the spin) as there are lattice sites in the solid. Hence, an insulator must have an
even number of electrons per unit cell. A metal will result if there is an odd number
of electrons per unit cell (unless the electron-electron interactions, which we have
neglected, are strong); as a result of band overlap, a metal can also result if there is
an even number of electrons per unit cell.
A semiconductor is an insulator with a small band gap. A good insulator will
have a band gap of Eg ∼ 4eV. At room temperature, the number of electrons which
will be excited out of the highest filled band and into the lowest empty band will
be ∼ e−Eg /2kB T ∼ 10−35 which is negligible. Hence, the filled and empty bands will
remain filled and empty despite thermal excitation. A semiconductor can have a band
gap of order Eg ∼ 0.25 − 1eV. As a result, the thermal excitation of electrons can be
as high as ∼ e−Eg /2kB T ∼ 10−2 . Hence, there will be a small number of carriers excited
into the empty band, and some conduction can occur. Doping a semiconductor with
impurities can enhance this.
The basic property of a metal is that it conducts electricity. Some insight into
electrical conduction can be gained from the classical equations of motion of a electron,
i.e. Drude theory:
d 1
~r = p~
dt m
d e
p~ = −eE(~r, t) − p~ × B(~r, t) (5.78)
dt m
If we continue to treat the electric and magnetic fields classically, but treat the elec-
trons in a periodic potential quantum mechanically, this is replaced by:
d 1~
~r = vn (k) = ∇ k n (k)
dt h̄
d~ h̄
h̄ k = −eE(~r, t) − e vn (k) × B(~r, t) − ~k (5.79)
dt τ
The final term in the second equation is the scattering rate. It is caused by effects
which we have neglected in our analysis thus far: impurities, phonons, and electron-
electron interactions. Without these effects, electrons would accelerate forever in a
constant electric field, and the conductivity would be infinite. As a result of scattering,
σ is finite. Hence, a finite electric field leads to a finite current:
~j =
XZ d3 k 1 ~
∇k n (k) (5.80)
n (2π)3 h̄
Filled bands give zero contribution to the current since they vanish by integration
by parts. Since an insulator has only filled or empty bands, it cannot carry cur-
rent. Hence, it is not characterized by its conductivity but, instead, by its dielectric
constant, .
5.9 Electrons in a Magnetic Field: Landau Bands

In 1879, E.H. Hall performed an experiment designed to determine the sign of the
current-carrying particles in metals. If we suppose that these particles have charge
e (with a sign to be determined) and mass m, the classical equations of motion of
charged particles in an electric field, E = Ex x̂ + Ey ŷ, and a magnetic field, B = Bẑ
are:
dpx
= eEx − ωc py − px /τ
dt
dpy
= eEy + ωc px − py /τ (5.81)
dt
where ωc = eB/m and τ is a relaxation rate determined by collisions with impurities,

other electrons, etc. These are the equations which we would expect for free particles.
In a crystalline solid, the momentum p~ must be replaced by the crystal momentum
and the velocity of an electron is no longer p~/m, but is, instead,
~ p (p)
~v (p) = ∇ (5.82)
We won’t worry about these subtleties for now. In the systems which we will be
considering, the electron density will be very small. Hence, the electrons will be close
to the bottom of the band, where we can approximate:
h̄2 k 2
(k) = 0 + + ... (5.83)
2mb
where mb is called the band mass. For instance, in the square lattice nearest-neighbor
tight-binding model,
(k) = −2t (cos kx a + cos ky a)

≈ −4t + ta2 k 2 + . . . (5.84)
Hence,
h̄2
mb = (5.85)
2ta2
In GaAs, mb ≈ 0.07me . Once we replace the mass of the electron by the band mass,
we can approximate our electrons by free electrons.
Let us, following Hall, place a wire along the x̂ direction in the above magnetic
fields and run a current, jx , through it. In the steady state, dpx /dt = dpy /dt = jy = 0,
m
we must have Ex = ne2 τ
jx and
B −e h Φ/Φ0
Ey = − jx = jx (5.86)
ne |e| e2 N
where n and N are the density and number of electrons in the wire, Φ is the magnetic
flux penetrating the wire, and Φ0 = h/e is the flux quantum. Hence, the sign of the
charge carriers can be determined from a measurement of the transverse voltage in
a magnetic field. Furthermore, according to (5.86), the density of charge carriers –
Figure 5.1: ρxx and ρxy vs. magnetic field, B, in the quantum Hall regime. A
number of integer and fractional plateaus can be clearly seen. This data was taken
at Princeton on a GaAs-AlGaAs heterostructure.
i.e. electrons – can be determined from the slope of the ρxy = Ey /jx vs B. At high
temperatures, this is roughly what is observed.
In the quantum Hall regime, namely at low-temperatures and high magnetic fields,
very different behavior is found in two-dimensional electron systems. ρxy passes
1 h
through a series of plateaus, ρxy = ν e2
, where ν is a rational number, at which
ρxx vanishes, as may be seen in Figure 5.1. The quantization is accurate to a few
parts in 108 , making this one of the most precise measurements of the fine structure
e2
constant, α = h̄c
, and, in fact, one of the highest precision experiments of any kind.
Some insight into this phenomenon can be gained by considering the quantum
mechanics of a single electron in a magnetic field. Let us suppose that the electron’s
motion is planar and that the magnetic field is perpendicular to the plane. For now,
we will assume that the electron is spin-polarized by the magnetic field and ignore
the spin degree of freedom. The Hamiltonian,
1
H= (−ih̄∇ + e A)2 (5.87)
2m
takes the form of a harmonic oscillator Hamiltonian in the gauge Ax = −By, Ay = 0.

(Here, and in what follows, I will take e = |e|; the charge of the electron is −e.) If
we write the wavefunction φ(x, y) = eikx x φ(y), then:
1 1

Hψ = (eB y + h̄kx )2 + (−ih̄∂y )2 φ(y) eikx x (5.88)
2m 2m
The energy levels En = (n+ 21 )h̄ωc , called Landau levels, are highly degenerate because
the energy is independent of k. To analyze this degeneracy, let us consider a system
of size Lx × Ly . If we assume periodic boundary conditions, then the allowed kx
values are 2πn/Lx for integer n. The harmonic oscillator wavefunctions are centered
at y = h̄k/(eB), i.e. they have spacing yn − yn−1 = h/(eBLx ). The number of
these which will fit in Ly is eBLx Ly /h = BA/Φ0 . In other words, there are as many
degenerate states in a Landau level as there are flux quanta.
It is often more convenient to work in symmetric gauge, A = 12 B × r Writing
z = x + iy, we have:
h̄2
" ! ! #
z̄ z 1
H= −2 ∂ − 2 ∂¯ + 2 + 2 (5.89)
m 4`0 4`0 2`0
with (unnormalized) energy eigenfunctions:

|z|2
− 2
m
ψn,m (z, z̄) = z Lm
n (z, z̄)e
4`
0 (5.90)
at energies En = (n + 21 )h̄ωc , where Lm

n (z, z̄) are the Laguerre polynomials and `0 =
q
h̄/(eB) is the magnetic length.
Let’s concentrate on the lowest Landau level, n = 0. The wavefunctions in the
lowest Landau level,
|z|2
−
m 4`2
ψn=0,m (z, z̄) = z e 0 (5.91)
are analytic functions of z multiplied by a Gaussian factor. The general lowest Landau
level wavefunction can be written:
|z|2
−
4`2
ψn=0,m (z, z̄) = f (z) e 0 (5.92)
The state ψn=0,m is concentrated on a narrow ring about the origin at radius rm =
q
`0 2(m + 1). Suppose the electron is confined to a disc in the plane of area A. Then
the highest m for which ψn=0,m lies within the disc is given by A = π rmmax , or,
simply, mmax + 1 = Φ/Φ0 , where Φ = BA is the total flux. Hence, we see that in the
thermodynamic limit, there are Φ/Φ0 degenerate single-electron states in the lowest
Landau level of a two-dimensional electron system penetrated by a uniform magnetic
flux Φ. The higher Landau levels have the same degeneracy. Higher Landau levels
can, at a qualitative level, be thought of as copies of the lowest Landau level. The
detailed structure of states in higher Landau levels is different, however.
Let us now imagine that we have not one, but many, electrons and let us ignore
the interactions between these electrons. To completely fill p Landau levels, we need
Ne = p(Φ/Φ0 ) electrons. Lorentz invariance tells us that if
e2
n=p B (5.93)
h
then
e2
jx = p Ey (5.94)
h
i.e.
e2
σxy = p (5.95)
h
The same result can be found by inverting the semi-classical resistivity matrix, and
substituting this electron number.
Suppose that we fix the chemical potential, µ. As the magnetic field is varied, the
energies of the Landau levels will shift relative to the chemical potential. However,
so long as the chemical potential lies between two Landau levels (see figure 5.2), an
integer number of Landau levels will be filled, and we expect to find the quantized
Hall conductance, (5.95).
These simple considerations neglected two factors which are crucial to the obser-
vation of the quantum Hall effect, namely the effects of impurities and inter-electron
interactions.1 The integer quantum Hall effect occurs in the regime in which impuri-
2
ties dominate; in the fractional quantum Hall effect, interactions dominate.
1
We also ignored the effects of the ions on the electrons. The periodic potential due to the lattice
has very little effect at the low densities relevant for the quantum Hall effect, except to replace
the bare electron mass by the band mass. This can be quantitatively important. For instance,
mb ' 0.07 me in GaAs.
2
The conventional measure of the purity of a quantum Hall device is the zero-field mobility, µ,
which is defined by µ = σ/ne, where σ is the zero-field conductivity. The integer quantum Hall effect
was first observed by von Klitzing, Pepper, and Dorda in Si mosfets with mobility ≈ 104 cm2 /Vs
Figure 5.2: (a) The density of states in a pure system. So long as the chemical poten-
tial lies between Landau levels, a quantized conductance is observed. (b) Hypothetical
density of states in a system with impurities. The Landau levels are broadened into
bands and some of the states are localized. The shaded regions denote extended
states. (c) As we mention later, numerical studies indicate that the extended state(s)
occur only at the center of the band.
5.9.1 The Integer Quantum Hall Effect
Let us model the effects of impurities by a random potential in which non-interacting

electrons move. Clearly, such a potential will break the degeneracy of the different
states in a Landau level. More worrisome, still, is the possibility that some of the
states might be localized by the random potential and therefore unable to carry any
current at all. As a result of impurities, the Landau levels are broadened into bands
and some of the states are localized. The possible effects of impurities are summarized
in the hypothetical density of states depicted in Figure 5.2.
e2
Hence, we would be led to naively expect that the Hall conductance is less than h
p
when p Landau levels are filled. In fact, this conclusion, though intuitive, is completely
wrong. In a very instructive calculation (at least from a pedagogical standpoint),
Prange analyzed the exactly solvable model of electrons in the lowest Landau level
interacting with a single δ-function impurity. In this case, a single localized state,
which carries no current, is formed. The current carried by each of the extended states
is increased so as to exactly compensate for the localized state, and the conductance
e2
remains at the quantized value, σxy = h
. This calculation gives an important hint of
the robustness of the quantization, but cannot be easily generalized to the physically
relevant situation in which there is a random distribution of impurities. To understand
while the fractional quantum Hall effect was first observed by Tsui, Störmer, and Gossard in GaAs-
AlGaAs heterostructures with mobility ≈ 105 cm2 /Vs. Today, the highest quality GaAs-AlGaAs
samples have mobilities of ≈ 107 cm2 /Vs.
Figure 5.3: (a) The Corbino annular geometry. (b) Hypothetical distribution of
energy levels as a function of radial distance.
the quantization of the Hall conductance in this more general setting, we will turn to
the beautiful arguments of Laughlin (and their refinement by Halperin), which relate
it to gauge invariance.
Let us consider a two-dimensional electron gas confined to an annulus such that
all of the impurities are confined to a smaller annulus, as shown in Figure 5.3. Since,
as an experimental fact, the quantum Hall effect is independent of the shape of the
sample, we can choose any geometry that we like. This one, the Corbino geometry,
is particularly convenient. States at radius r will have energies similar to to those
depicted in Figure 5.3.
Outside the impurity region, there will simply be a Landau level, with energies
that are pushed up at the edges of the sample by the walls (or a smooth confining
potential). In the impurity region, the Landau level will broaden into a band. Let us
suppose that the chemical potential, µ, is above the lowest Landau level, µ > h̄ωc /2.
Then the only states at the chemical potential are at the inner and outer edges of the
annulus and, possibly, in the impurity region. Let us further assume that the states
at the chemical potential in the impurity region – if there are any – are all localized.
Now, let us slowly thread a time-dependent flux Φ(t) through the center of the
annulus. Locally, the associated vector potential is pure gauge. Hence, localized
states, which do not wind around the annulus, are completely unaffected by the flux.
Only extended states can be affected by the flux.
When an integer number of flux quanta thread the annulus, Φ(t) = pΦ0 , the
flux can be gauged away everywhere in the annulus. As a result, the Hamiltonian
in the annulus is gauge equivalent to the zero-flux Hamiltonian. Then, according
to the adiabatic theorem, the system will be in some eigenstate of the Φ(t) = 0
Hamiltonian. In other words, the single-electron states will be unchanged. The
only possible difference will be in the occupancies of the extended states near the
chemical potential. Localized states are unaffected by the flux; states far from the
chemical potential will be unable to make transitions to unoccupied states because
the excitation energies associated with a slowly-varying flux will be too small. Hence,
the only states that will be affected are the gapless states at the inner and outer edges.
Since, by construction, these states are unaffected by impurities, we know how they
are affected by the flux: each flux quantum removes an electron from the inner edge
R R R dΦ
and adds an electron to the outer edge. Then, I dt = e and V dt = dt
= h/e,
so:
e2
I= V (5.96)
h
Clearly, the key assumption is that there are no extended states at the chemical
potential in the impurity region. If there were – and there probably are in samples
that are too dirty to exhibit the quantum Hall effect – then the above arguments
break down. Numerical studies indicate that, so long as the strength of the impurity
potential is small compared to h̄ωc , extended states exist only at the center of the
Landau band (see Figure 5.2). Hence, if the chemical potential is above the center of
the band, the conditions of our discussion are satisfied.
The other crucial assumption, emphasized by Halperin, is that there are gapless
states at the edges of the system. In the special setup which we assumed, this was
guaranteed because there were no impurities at the edges. In the integer quantum
Hall effect, these gapless states are a one-dimensional chiral Fermi liquid. Impurities
are not expected to affect this because there can be no backscattering in a totally
chiral system. More general arguments, which we will mention in the context of
the fractional quantum Hall effect, relate the existence of gapless edge excitations to
gauge invariance.
One might, at first, be left with the uneasy feeling that these gauge invariance
arguments are somehow too ‘slick.’ To allay these worries, consider the annulus with
a wedge cut out, which is topologically equivalent to a rectangle. In such a case,
some of the Hall current will be carried by the edge states at the two cuts (i.e. the
edges which run radially at fixed azimuthal angle). However, probes which measure
the Hall voltage between the two cuts will effectively couple these two edges leading,
once again, to annular topology.
Laughlin’s argument for exact quantization will apply to the fractional quantum
Hall effect if we can show that the clean system has a gap. Then, we can argue that
for an annular setup similar to the above there are no extended states at the chemical
potential except at the edge. Then, if threading q flux quanta removes p electrons
from the inner edge and adds p to the outer edge, as we would expect at ν = p/q, we
p e2
would have σxy = q h
.

Solid State Physics

Uploaded by

Document Informationclick to expand document information

Copyright:

Available Formats

Solid State Physics

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Solid State Physics

Uploaded by

Copyright:

Available Formats

Solid State Physics

1 What is Condensed Matter Physics? 1

2 Review of Quantum Mechanics 7

3 Review of Statistical Mechanics 18

4 Broken Translational Invariance in the Solid State 30

What is Condensed Matter

1.1 Length, time, energy scales

as compared to energies in the M eV for nuclear matter, and GeV or even T eV , in

1.2 Microscopic Equations vs. States of Matter,

Or, perhaps, the Schrödinger equation:

1.3 Broken Symmetries

1.4 Experimental probes: X-ray scattering, neu-

• Thermodynamics: measure the response of macroscopic variables such as the

• Transport: set up a potential or thermal gradient, ∇ϕ, ∇T and measure the

1.5 The Solid State: metals, insulators, magnets,

1.6 Other phases: liquid crystals, quasicrystals,

Review of Quantum Mechanics

2.1 States and Operators

space. These operators correspond to physical observables. Finally, there is an inner

be the average values of the corresponding physical observables if we could measure

The adjoint, O† , of an operator is defined according to

while a unitary operator satisfies

is unitary. Given an Hermitian operator, O, its eigenstates are orthogonal,

A particularly important operator is the Hamiltonian, or the total energy, which

If the Hamiltonian is independent of time, then we can define energy eigenstates,

which evolve in time according to:

An arbitrary state can be expanded in the basis of energy eigenstates:

It will evolve according to:

x δ(x − a) = a δ(x − a) (2.15)

The momentum eigenfunctions are plane waves:

where the Fourier coefficients are given by:

If the particle is free,

2.2 Density and Current

The current carried by a plane-wave state is:

2.3 δ-function scatterer

There is a bound state at:

2.4 Particle in a Box

ψ(x) = Aeikx + Be−ikx (2.34)

has energy E = h̄2 k 2 /2m. ψ(0) = ψ(L) = 0. Therefore,

for integer n. Allowed energies

2.5 Harmonic Oscillator

Raising and lowering operators:

Hamiltonian and commutation relations:

The commutation relations,

[H, a] = −h̄ωa (2.44)

imply that there is a ladder of states,

Ha† |Ei = (E + h̄ω) a† |Ei

Such a state has E0 = h̄ω/2.

|cn |2 = hn|aa† |ni

2.6 Double Well

The allowed k’s are now determined by

α|+i + β|−i (2.60)

spanned by the basis vectors |±i. Spin operators:

States of a spin in a magnetic field in the ẑ direction:

2.8 Many-Particle Hilbert Spaces: Bosons, Fermions

spanned by the bases

Then the two-particle Hilbert space is spanned by the set:

Applying this relation twice implies that

while a fermionic wavefunction must be completely antisymmetric:

CV = (const.) · kB · g (F ) kB T (3.65)

where is ~ a unit polarization vector, satisfy

For longitudinally polarized waves, ~k = k~,

while transverse waves, ~k · ~ = 0 have