Statistical Physics Overview and Concepts
Statistical Physics Overview and Concepts
Gianni Blatter
2020
ii
Thanks: I wish to thank Pascal Steger for his help in writing this elec-
tronic form of the script and numerous students for their valuable comments.
G. Blatter
0 Literature 1
1 Thermodynamics 3
1.1 Basic concepts and definitions . . . . . . . . . . . . . . . . . . 4
1.1.1 State variables, equation of state . . . . . . . . . . . . 4
1.1.2 Functions of many variables, differentials, integrability 5
1.1.3 Process types . . . . . . . . . . . . . . . . . . . . . . . 7
1.1.4 Units and coefficients . . . . . . . . . . . . . . . . . . 9
1.2 Three basic Laws of Thermodynamics . . . . . . . . . . . . . 11
1.2.1 The First Law of Thermodynamics, energy . . . . . . 11
1.2.2 The Second Law of Thermodynamics, entropy . . . . . 14
1.2.3 The third law of thermodynamics, zero temperature . 17
1.3 Multicomponent systems . . . . . . . . . . . . . . . . . . . . . 20
1.4 Thermodynamic potentials . . . . . . . . . . . . . . . . . . . 21
1.4.1 Thermodynamic potentials with extensive variables . . 21
1.4.2 Constraining parameters . . . . . . . . . . . . . . . . . 22
1.4.3 Legendre transformation . . . . . . . . . . . . . . . . . 24
1.4.4 Thermodynamic potentials with intensive variables . . 26
1.4.5 Equation of Gibbs-Duhem . . . . . . . . . . . . . . . . 28
1.5 Phenomenology of phase transitions . . . . . . . . . . . . . . 30
1.5.1 Gibbs’ phase rule . . . . . . . . . . . . . . . . . . . . . 30
1.5.2 Gibbs surfaces . . . . . . . . . . . . . . . . . . . . . . 32
1.5.3 Ehrenfest classification of phase transitions . . . . . . 33
1.5.4 Phase diagrams . . . . . . . . . . . . . . . . . . . . . . 34
iii
iv CONTENTS
2 Statistical Description 45
6 Fermions 99
6.0.1 Dilute Fermi gas, δ = nλ3 1 . . . . . . . . . . . . . 101
6.0.2 Degenerate Fermi gas, δ = nλ3 1 . . . . . . . . . . . 101
6.1 Electrons in a magnetic field . . . . . . . . . . . . . . . . . . . 106
6.1.1 Landau diamagnetism . . . . . . . . . . . . . . . . . . 108
6.1.2 Pauli spin paramagnetism . . . . . . . . . . . . . . . . 113
7 Bosons 117
7.1 Photons . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117
7.2 Phonons . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 122
7.3 Ideal Bose-gas:
Bose-Einstein-condensation . . . . . . . . . . . . . . . . . . . 125
16 k-space RG 259
16.1 Gauss model . . . . . . . . . . . . . . . . . . . . . . . . . . . 261
16.2 Φ4 theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 264
16.3 Dangerously irrelevant variables . . . . . . . . . . . . . . . . . 274
16.4 β-function and ‘charge’ . . . . . . . . . . . . . . . . . . . . . . 275
16.5 Two-loop expansion for G and finite η . . . . . . . . . . . . . 277
Literature
K. Huang, Statistical Mechanics (John Wiley & Sons, New York, 1987),
good survey.
A. Sommerfeld, Thermodynamik und Statistik (Vorlesungen über theo-
retische Physik. Band V, Harri Deutsch, 1977), good survey, old style.
L.E. Reichl, A Modern Course in Statistical Physics (E. Arnold (Publish-
ers), 1980, 1987), good survey.
D. Chandler, Introduction to Modern Statistical Physics (Oxford Univer-
sity Press, New York, 1987), good survey.
H. Smith, H. Jensen, Transport Phenomena, (Claredon Press, Oxford,
1989), Transportphenomena.
N. Straumann, Thermodynamik (Lecture Notes in Physics, Springer, Berlin,
1986), Focus on Thermodynamics.
S.K. Ma, Statistical Mechanics (World Scientific, Singapore, 1985), Focus
on Statistical Mechanics.
R. P. Feynman, Statistical Mechanics, a set of lectures (Frontiers in Physics,
Benjamin/Cummings, Reading, Massachusetts, 1982), Path Integrals.
N. Goldenfeld, Lectures on Phase Transitions and the Renormalization
Group (Frontiers in Physics, Addison-Wesley, Reading, Massachusetts,
1994), Phase Transitions, Renormalization Group Theory.
M. Plischke, B. Bergersen, Equilibrium Statistical Mechanics (Prentice-
Hall, 1989), Phase Transitions.
S.K. Ma, Modern Theory of Critical Phenomena (Frontiers in Physics,
Benjamin/Cummings, Reading, Massachusetts, 1976), Phase Transi-
tions, Renormalization Group Theory.
1
2 CHAPTER 0. LITERATURE
Thermodynamics
3
4 CHAPTER 1. THERMODYNAMICS
famous examples being those of the ideal gas, its equation of state pV = nRT
(n = moles of substance, R = 8.314 J/mol K the gas constant) and its caloric
equation U = (3/2)nRT . Note that in statistical mechanics these laws
read pV = N kB T and U = (3/2)N kB T with the Boltzmann constant kB =
1.381 10−23 J/K and n = N/NA , R = kB NA with NA = 6.022 1023 /mol the
Avogadro number. Furthermore, note the R is a phenomenological constant
while kB is a microscopic quantity.
B
A f ( p,V,T ) = 0
γ1 T
B
A γ2
p
Fig. 1.1: States A and B in state space described by the thermodynamic vari-
ables p, V and T . Thermodynamic equilibrium states are defined through
the state function f (p, V, T ) = 0; the condition f (p, V, T ) = 0 defines the
surface of allowed equilibrium states. The state B can be reached from the
state A via many different paths. A state function Z has the characteristic
property that its value ZA in A can be obtained from the value ZB in B via
arbitrary paths γ connecting A and B.
integral condition on Z
I
dZ = 0 (1.5)
that guarantees its integrability. When dealing with two independent vari-
ables x and y, any 1-form Xdx + Y dy can be made exact with the help of an
integrating factor. An example in thermodynamics is the relation between
entropy and heat; the heat δQ(V, T ) is not an exact differential (hence we
write δ) but the product T −1 δQ = dS is. As a result, S is a state func-
tion but Q is not. The statement that dS is an exact differential is the
Second Law of Thermodynamics in mathematical form. In the same vain,
we have that δW and δQ are not exact differentials but the combination
dU = δQ − δW is—this this just the statement given by the First Law of
Thermodynamics in mathematical form.
Integrability becomes more cumbersome when dealing with more than
two variables.
Pn The integrability of a 1-form with n variables (or in n di-
mensions) k Xk dxk requires, that the n(n − 1)/2 conditions ∂Xi /∂xk =
∂Xk /∂xi are satisfied. This can be written in compact form as the condition
~ x ) = X,
for the vanishing of the rotation of the n-dimensional vector field Z(~ ~
~ ~ ~ ~
∇ ∧ Z = 0. For n = 3 we can write ∇ ∧ Z as a vector; for n = 4 we ob-
tain an anti-symmetric tensor. Integrability of a 1-form in three dimensions
1.1. BASIC CONCEPTS AND DEFINITIONS 7
that follows from z(x, y = const) via dz = ∂x z|y dx, 1 = ∂x z|y ∂z x|y . Let
us include the additional function w(x, y) in our considerations. We can
construct the pair of functions x(y, w) und y(z, w) and their differentials
∂x ∂x ∂y ∂y
dx = dy + dw, dy = dz + dw.
∂y w ∂w y ∂z w ∂w z
Fixing w and inserting one in the other, we obtain dx = ∂y x|w ∂z y|w dz and
hence the chain rule at fixed w,
∂x ∂y ∂x
= . (1.9)
∂y w ∂z w ∂z w
irreversible
T ? ?
A f ( p,V,T ) = 0
reversible
B
A
p
Fig. 1.2: While a reversible process proceeds within the equilibrium manifold
defined through f (p, V, T ) = 0, an irreversible one leaves the surface and has
to return to it via relaxation. While the reversible process can be described
by a minimal number/evolution of thermodynamic variables, the irreversible
process cannot be easily handled.
0 VA VB 0 VA VB
T T
Fig. 1.3: Isothermal expansion of a gas: in the reversible process on the left,
the cylinder picks up heat from the reservoir and can convert it into work
extracted from the displacement of the piston, hence quite a useful process.
The irreversible process on the right where the gas just flows over in the
second compartment misses this opportunity—there is no heat transferred
from the reservoir to the gas and there is no work done by its expansion.
When attempting to calculate the entropy change of the right process, we
substitute it by the left one and obtain the result ∆S = nR ln(VB /VA ).
We measure
temperature T
the
pressure p
the
10 CHAPTER 1. THERMODYNAMICS
heat δQ
the
amount n of a substance
in mol, g,
X
1 mol substance ≈ atomic weight · g, [gram-mol].
atoms in the substance
(1.14)
E.g., one mole of water
ˆ 2 · 1(H) + 1 · 16(O) = 18 g
1 mol H2 O = (1.15)
1 mol ˆ 6.022 ·
= 1023 molecules
1 mol 12
C = 6.022 · 1023 · 12 · 1.66044 · 10−24 g = 12 g, (1.16)
Atom.
Material coefficients
dU = δQ − δW, (1.18)
Alternatively, the first law states that any loop integral (or cyclic process)
over the combination δQ − δW vanishes,
I
δQ − δW = 0, (1.20)
This finding is accurate for an ideal gas; for a real gas (e.g., a van der Waals
gas, see later) the internal energy turns out to depend on both T and V .
This is clear if we think in microscopic terms (which we should do in stat
mech but not in thdyn, but we do it anyway), as the atoms/molecules in
the Van der Waals gas attract, and hence expanding the gas enlarges the
distance between the particles and kinetic energy is transformed to potential
energy, thus cooling the gas internally. For a monoatomic ideal gas, CV is a
constant, hence U = CV T . This is a typical application of thermodynamic
principles—do an experiment (phenomenology) and use the thermodynamic
laws to obtain new insights, here, a simplification of a general state function
U (V, T ) to a simpler one U (T ).
In extending the above exercise, we compare a process that changes the
gas’ state at fixed volume V versus one where the pressure p is fixed,
where we have assumed the same change in temperature dT for the two
processes. Hence, the First Law gives us a new insight on the relation
between specific heats,
Cp − CV = nR > 0, (1.25)
thus relating experiments that probe different caloric properties of the sys-
tem. Again, thinking about the result, it is clear that Cp > CV as in the
second experiment the system produces work (on a piston that holds p)
while changing its temperature by the same amount, so the system has to
pick up more heat in the second process.
Finally, we discuss the adjabatic reversible expansion of an ideal gas
which is characterized by δQ = 0 and δW = p dV , see Fig. 1.4. Making use
of the First Law and the caloric equation dU = CV dT , we find that
CV dT = dU = δQ − δW = −p dV, (1.26)
1.2. THREE BASIC LAWS OF THERMODYNAMICS 13
cylinder
Fig. 1.4: Geometry of the
p cylinder-piston configuration in
T
the adiabatic-reversible expansion.
V piston
δQ = 0
Clausius: heat can never pass from a colder to a warmer body without
some other change, connected therewith and occurring at the same time.
T1 T1
Q1 Q1
W W
Q2 Q2
T2 T2
Fig. 1.6: Carnot’s cyclic processes: the engine absorbs heat Q1 from the hot
reservoir and transforms part of it into work W = Q1 − Q2 , while the rest
flows as waste heat Q2 into the cold reservoir. The pump or refrigerator
consumes work W = Q1 − Q2 to transfer the heat Q2 to the hot reservoir
(all quantities W , Q1 , and Q2 are positive in this accounting).
Q1 T1
= . (1.30)
Q2 T2
4 3
T2
V
The differential of S is
∂S ∂S
T dS = T dT + T dV
∂T V ∂V T
and comparing the two expressions for T dS, we obtain the following relations
between partial derivatives,
∂S 1 ∂U ∂S 1 ∂U
= , = +p .
∂T V T ∂T V ∂V T T ∂V T
The Third Law of Thermodynamics (due to Planck and Nernst) tells, that
T
Cp (p, T 0 )
Z
S(p, T ) = dT 0 + S(p, 0), (1.47)
0 T0
and set S(V, 0) = 0, S(p, 0) = 0, since S(T → 0) becomes independent on V
and p. In order for the integrals not to diverge as T → 0, the specific heats
Cp , CV (T → 0) → 0 must vanish. From condensed matter courses, we know
that indeed CV ∝ T 3 for an elastic system (phonons) and CV ∝ T in the
electron system of a metal.
1.2. THREE BASIC LAWS OF THERMODYNAMICS 19
Checking matters for the ideal gas, we run into a problem—the entropy
S(V, T ) in Eq. (1.40) is ill defined as T → 0 and the specific heat CV does
not approach 0. Again, quantum mechanics comes to our rescue, telling us
that the ideal classical gas is a high temperature approximation of physical
ideal gases which are made from bosons or fermions. Ideal fermions then
become degenerate as T is lowered and CV ∝ T , while bosons condense and
CV ∝ T 3/2 as T → 0.
An important consequence of the Third Law is the statement that the ab-
solute zero in temperature cannot be reached. We show this for a mechanical
process—the translation to other processes, e.g., adjabatic demagnetization,
is straightforward. Using (1.43) and (1.47) we can write
T
Cp (p, T 0 ) T
dT 0
Z Z
∂S ∂ ∂Cp
Vα=− =− dT 0 =− .
∂p T ∂p 0 T0 0 ∂p T0 T0
Vα
⇒ → const 6= 0 für T → 0 (1.48)
Cp
T dS = Cp dT − T V α dp,
Then use that the prefactor of dT is Cp and the prefactor of dp is T V α, with α the thermal
expansion coefficient.
20 CHAPTER 1. THERMODYNAMICS
When dealing with more complex systems we have to account for many
components (say r, components i ∈ {1, . . . , r}) that can occur in various
phases or aggregates (say ν; possible phases α ∈ {1, . . . , ν} may be solid,
liquid, or gas). A component i can show up in any of these phases α, hence
the total amount of component i is given by
ν
(α)
X
ni = ni , (1.51)
α=1
(α)
where we count the amount of substance in moles and ni quantifies the
amount of substance i in the phase α ∈ {1, . . . , ν}. Similarly, the total
internal energy, entropy, and volume of the overall system results from the
sum over its phases
X X X
U= U (α) , S = S (α) , V = V (α) . (1.52)
α α α
(α)
We define the chemical potential of the i-th component in phase α, µi , as
(α)
the change in energy under a reversible change in the amount of moles ni ,
(α) ∂U
µi = (α)
. (1.53)
(β) (α)
∂ni S,V,nj without ni
(α)
The chemical potentials µi are equilibrium parameters, the same way as
the other intensive variables p and T . Indeed, a system in thermal equi-
librium has a uniform distribution of temperature, pressure, and chemical
potentials throughout its phases,4
(µ) (ν)
T (µ) = T (ν) = T, p(µ) = p(ν) = p, µi = µi = µi . (1.54)
A gradient (or difference) in the chemical potential between two phases will
lead to a phase transformation between the phases of that component such
as to equalize the chemical potentials.
Finally, the differentials of the entropy and the internal energy pick up
additional terms accounting for a change in mole (n) or particle (N ) number,
r
X
T dS = dU + pdV − µi dni , (1.55)
i=1
r
X
dU = T dS − p dV + µi dni .
i=1
4
The intensive variables are also homogeneous throughtout the individual phases, see
Sec. 1.4.2
1.4. THERMODYNAMIC POTENTIALS 21
Here, ai are the intensive variables dual to Ai , obvious pairs being ai and
Ai = −p and V , µ and N , σ = surface tension and A = area, φ = electric po-
tential and Q = charge, E = electric field and P = electric polarization, H =
magnetic field and M = magnetic moment, etc. The set U, A1 , . . . , An , with
U the caloric variable, is assumed to provide a complete characterization of
the system’s state through extensive variables. The equation (1.56) is known
as the Gibbs fundamental equation. Its knowledge allows to determine the
caloric and all state equations of the system, that is why S(U, A1 , . . . , An )
is called a thermodynamic potential.
In order to find the equations of state, we exploit the integrability of
dS, guaranteeing the existence of the state function S(U, A1 , . . . , An ). Its
expansion
n
∂S X ∂S
dS = dU + dAk (1.57)
∂U A1 ,...,An ∂Ak U,A1 ,...,An without Ak
k
22 CHAPTER 1. THERMODYNAMICS
1 ∂S
= → T = T (U, A1 , . . . , An ) → U = U (T, A1 , . . . , An ),
T ∂U A1 ,...,An
∂S
ak = −T → ak = ak (T, A1 , . . . , An ).
∂Ak U,A1 ,...,An no Ak
These are nothing but the caloric and thermal equations of state. Hence,
knowing a thermodynamic potential of a system allows for the immediate
derivation of its equations of state and thus its full thermodynamic charac-
terization.
The second thermodynamic potential follows immediately from the en-
tropy. Indeed, the inner energy U is a state function as well and expressing it
through the extensive variables S, A1 , . . . , An , we obtain the thermodynamic
potential U (S, A1 , . . . , An ). Its differential reads
n
X
dU = T dS + ak dAk (1.58)
k
P P’
p1, V1 p2 , V2
0 x
Fig. 1.9: Piston system with a constraint. The system volume V is parti-
tioned into two parts V1 and V2 such that V = V1 + V2 . The inner piston
P’ can be moved in order to change the volumes V1 and V2 , the constraint
demanding that ∆V1 = −∆V2 .
Letting P’ free, the piston will move such as to reach equilibrium and po-
tentially produce (useful) work. Thereby, the entropy will increase (see Fig.
1.10)
δS|U,V > 0. (1.62)
0 V1
T1 = T2 in equilibrium. (1.63)
S, V fixed
Fig. 1.11: The behavior of
U the internal energy U under a
(virtual) constraint Y .
0 Y
yx
f (x) Fig. 1.12: Legendre trans-
formation Lf (y) of the func-
y tion f (x).
x
L f ( y)
to introduce an additional operation, the Legendre transformation. The
Legendre transformation is a tool that allows to substitute variables without
loss of information.
We start from the fact that S is concave in U and V , i.e.,
f Lf
C
A
B BC
B
x y
AB
∂F ∂F ∂F
= −p, = −S, = µ, (1.72)
∂V T,N ∂T V,N ∂N T,V
and the condition for the integrability results in one of the famous Maxwell
relations
∂p ∂S
= . (1.73)
∂T V ∂V T
Replacing the entropy dS by dS = δQrev /T , we find the relation
∂p 1 δQ
= , (1.74)
∂T V T dV T
1.4. THERMODYNAMIC POTENTIALS 27
that will appear in a similar form as the Clapeyron equation in the context
of phase transitions.
By construction, F is convex in the extensive variables V aund N und
concave (note the different signs in (1.68) und (1.70)) in the intensive vari-
able T and the following minimum principle applies,
Enthalpy H
dH = T dS + V dp + µ dN,
Gibbs potential G
The Gibbs potential G(T, p, N ) derives from the internal energy U (S, V, N )
by replacing S and V ,
dG = −S dT + V dp + µ dN,
Grand potential Ω
Maxwell-relations
dU = T dS + S dT − p dV + V dp + µ dN + N dµ,
and exploiting the usual form of the differential dU in its natural variables
S, V , and N ,
dU = T dS − p dV + µ dN, (1.85)
0 = S dT − V dp + N dµ, (1.86)
stating, that the set of intensive variables is not independent. Indeed, the
function p(S, V, Nk ) is homogeneous of order zero and we have
p = p(S/N, V /N, x1 , . . . , xr )
= p(S/N, V /N, x1 , . . . , xr−1 , 1 − x1 − x2 − . . . − xr−1 ). (1.88)
U − TS = F
F + pV = G
G − µn = 0
⇒ 0 = SdT − V dp + ndµ
in the variables T, p, µ.
30 CHAPTER 1. THERMODYNAMICS
Indeed, all these phases, gas, liquid, and solid do exist under appropri-
ate conditions and transform into one another in so-called phase transitions.
These phases and transitions then are conveniently represented in so-called
phase diagrams. While we discuss here the phenomenological aspects of
phase transitions, their ‘microscopic’ properties and understanding will be
one of the central topics of statistical mechanics. To fix ideas, we will con-
sider a system that can be characterized by the parameters p and V , T and
S, and µ and N (or mole number n)—other systems may be characterized
by alternative pairs, e.g., a magnetic system is described by H and M rather
than p and V (although a piece of magnetic material requires both pairs for
its complete description, e.g., in magnetostriction phenomena the volume V
is non-trivially related to M ).
order into this zoo of possibilities by outlining the basic structure a system’s
phase diagram can assume. Here, we discuss a specific situation where
the pressure p, temperature T , and amount of substances (in moles) ni ,
i ∈ {1, . . . , r} are prescribed—the situation where other ‘external’ parame-
ters are implosed can be handled in a similar way. The appropriate thermo-
(α)
dynamic potential for these variables is G and the chemical potentials µi
of the i-th substance in the α-th phase, α ∈ {1, . . . , ν},
(α) (α) (α)
µi p, T, x1 , . . . , xr−1 , (1.89)
f = 2 + r − ν. (1.91)
cross in a single point, the triple point pαβγ , Tαβγ . According to Gibbs’
phase rule there is no possibility for the coexistence of four phases in the
one-component system with only two intensive variable p and T . A phase
diagram may have a form as sketched in Fig. 1.15. The critical endpoint of
the α-β phase line can be understood when studying the condition (1.92) in
more detail.
Every phase α defines a surface µ(α) (p, T ) (the Gibbs surface) in µ–p–T
space, see Fig. 1.16. Consider the situation (e.g., at appropriate pressure p
µ µ
β β
C
p p
1.
pαβ(T) pαβ(T)
α α
T T 2.
Fig. 1.16: Left: Gibbs surface for two phases α and β. According to the
thermodynamic minimization principle, it is the phase with a lower energy
that will be realised at given p and T . The phases transform into one
another at the crossing line between the surfaces that defines the transition
line pαβ (T ) when projected into the p–T plane. The sudden change of slope
at the crossing of the minimal surfaces implies jumps in the conjugated
variables of p and T , the volume V = ∂p G and the entropy S = −∂T G.
These jumps are related by Clapeyron’s equation. Right: Gibbs surface
describing the situation with a second-order transition: the line where the
two Gibbs surfaces cut at an angle, thus defining a first-order transition,
terminates in a critical point C describing a second-order phase transition.
and temperature T near a transition line pαβ (T )) where two phases α and
β are relevant, see Fig. 1.16. In equilibrium, the Gibbs energy
is always minimal, i.e., for p > pαβ phase β is realized and we have n(α) = 0,
while for p < pαβ the phase α has lower energy and hence n(β) = 0. In
general, the Gibbs surfaces cut under a finite angle such that the extensive
variables V (conjugate to p) und S (conjugate to T ; here, we use molar
quantities s and v) undergo a jump at the phase transformation,
∂µ(α) ∂µ(β)
= − 6= 0, (1.95)
∂p T ∂p T
∆s(T ) = s(α) pαβ (T ), T − s(β) pαβ (T ), T
∂µ(α) ∂µ(β)
= − + 6= 0.
∂T pαβ ∂T pαβ
T ∆s = ` (1.97)
has to be provided. This latent heat T dS modifies both the system’s internal
energy dU as well as its volume dV (as the system provides work under
expansion), T dS = dU + p dV . The CC (Clausius-Clapeyron)-equation
then reads
dpαβ `
T = . (1.98)
dT ∆v
Measuring ` and ∆v at fixed p and T for different values of T allows us to
find the transition line pαβ through simple integration.
namic limit (i.e., in a system with infinitely many degrees of freedom5 ) and
usually manifests itself through jumps in derivatives. A system described
with a Gibbs potential that has jumps in the first derivatives at the tran-
sition is undergoing a first-order phase transition, see the above example.
A phase transition is called to be of n-th order if the jumps appear in the
n-th derivative. Other (modern) terminology distinguishes between transi-
tions with a latent heat (first-order, with a jump in the order parameter)
and second-order transitions (smooth development of the order parameter,
with divergent susceptibiliy and correlation length), also sometimes called
continuous transitions. But beware, that there are also infinite-order tran-
sitions (that do not break a symmetry and have no order parameter), the
most well known example being the Kosterlitz-Thouless transition; some-
times, ‘continuous’ is referring to an infinite-order transition. One way to
arrive at a second-order transition is via the termination of a first-order
line where two Gibbs surfaces cut, see Fig. 1.16. The point C is called a
critical point or critical end-point. Note that the typical characteristics of a
second-order transition (e.g., the critical exponents) near a critical end-point
manifest themselves along the continuation of the first-order phase line, in
our example the continuation of pαβ (T ).
Often, it is convenient to study other diagrams, different from the p–T phase-
diagram, in order to understand and characterize a system. Alternative
diagrams are those providing the pressure p or temperature T as a function
of the volume V , see Fig. 1.17.
p p
C T
α p
pαβ
pαβ T pαβ(T)
β
T v v (α ) v (β ) v v (α ) v (β ) v
features associated with the jump in volume at the transition. Hence, while
the p–T diagram tells everything about the existing phases and transitions,
replacing p by V to arrive at a T –V diagram adds information on the jump
in volume at fixed p, but no longer provides a survey of phases.
p p
T T
C l C
l s K g
K
g T 3−phase
coexistence
v v
The three diagrams in Fig. 1.17 are nothing but the projections of the
surface of thermodynamic states in p-v-T space describing a one-component
system with two phases α and β (e.g., gas and liquid). A typical example is
shown in Fig. 1.18, together with the corresponding surface for a three-phase
system including a third phase γ (e.g., a solid phase). The ruled surfaces K
are made from straight lines (projecting onto individual points in the p–T
phase diagram) and describe regions of phase coexistence; their projection
into the p, T plane defines the transition lines, e.g., pαβ (T ) in Fig. 1.17. The
phases α and β cannot be distiguished from each other beyond the critical
point C. The phase diagram for the three-phase system is shown in Fig.
1.19 and includes a triple-point T .
namic potentials that actually generate the latter. We have started from
the ‘good’ variables p and T for which the Gibbs potential generates the
p–T phase diagram. When going to the variables v and T one may ask what
characteristic form the associated thermodynamic potential f (V, T ) actually
has in the presence of a transition—this is the subject of the next section.
Let us consider the free energy (per mole) f (v, T ). The phases α and β define
the free energies f (α) and f (β) . The (convex) isotherm f (v, T = const) must
have the form as shown in Fig. 1.20, with the tangents t(α) = t(β) = t in v (α)
and in v (β) connecting the states A and B following from the relation
∂f
p=− (1.99)
∂v T
v (α) v (β) v
Along the tangent, the free energy assumes the linear form
v (β) (T ) − v
f (t) (T, v) = f (α) T, v (α) (t)
(1.102)
v (β) (T ) − v (α) (T )
v − v (α) (T )
+f (β) T, v (β) (T ) (β)
v (T ) − v (α) (T )
1.5. PHENOMENOLOGY OF PHASE TRANSITIONS 37
f
T = const Fig. 1.21: Two-phase coexistence
α
α & β, metastable states α0 and β 0 ,
α’ A’ unstable and unstable regime with negative
A
B’ compressibility.
α&β β’
β
B
v (α) v (β) v
It is interesting to return from the free energy f (v, T ) in Fig. 1.20 back
to the Gibbs energy g(p, T ) in the slightly different formulation involving
38 CHAPTER 1. THERMODYNAMICS
T
Fig. 1.23: Free energy at fixed tem-
f perature T . The transition occurs
at the pressure p = plg (T ) and cor-
responds, via ∂v f = −plg (T ), to the
straight segment in f .
0 vli vgas v
In returning to the Gibbs free energy µ = f + pv, we keep the volume (or
density) as a constraint parameter and plot µ(p, T ; v) = f (v, T )+p v at fixed
p and T as a function of v in Fig. 1.24. The minima of the constraint Gibbs
metastable
µ p > plg
gas phase
stable
liquid
phase
p = plg (T )
first−order
p < p lg
phase transition
at p = p lg
constraining
0 −1 parameter
ρ−1
liquid ρ−1
gas ρ=v
for µ
metastable liquid stable
phase gas phase
s→u→f →µ
l C
s
µ
g
T s
l
C
s C
T g
l C
u T g p s l
s g
T
v T
the surface tension (note that such an interface has an intermediate density
ρgas < ρint < ρli and hence costs more energy, see Fig. 1.24). Hence, the
droplet growth involves the competition between a volume gain and a surface
cost, totalling to a Gibbs free energy
4π 3
G(r) = r ρli (µli − µgas ) + 4πr2 σlg . (1.105)
3NA
Depending on the pressure p, the function G(r) changes its behavior, see
Fig. 1.27: At small radii r, the positive quadratic term always dominates.
When increasing the pressure across plg (T ), the cubic term changes sign and
produces a maximum at the critical radius rc as obtained from the condition
δG = 0, 7
and using the liquid’s properties of small volume vli vgas = RT /p and
poor compressibility ∂p ρli ≈ 0, we arrive at the differential equation
Making use of the boundary condition p(rc = ∞) = plg , we obtain the final
results (see Fig. 1.28)
2σ v 1 2σ vli
li
p = plg exp and rc (p, T ) = . (1.110)
RT rc RT ln(p/plg )
p (T ) p
lg
3
He −A Fig. 1.29: Nucleation of a 3 He-B
droplet in 3 He-A as an example of
a quantum phase transition. Simi-
3
lar ideas of ‘false vacuum’ have been
He −B formulated about our universe (see
S. Coleman)
44 CHAPTER 1. THERMODYNAMICS
Chapter 2
Statistical Description
We will now proceed with the microscopic derivation and extension of the
results from thermodynamics. To this end, we will consider a many-body
system which can be described by the coordinates
∂H ∂H
q̇i = , ṗi = − .
∂pi ∂qi
45
46 CHAPTER 2. STATISTICAL DESCRIPTION
While, these computations still fall short of describing systems with 1023
particles which can be experimentally observed over milliseconds to days,
such huge simulations are often not required and smaller systems/times will
do. Another important aspect, however, is that such computations produce
vast quantities of data that do not really interest us in their great detail,
but from which we will only extract certain specific characteristics of the
system.
In the statistical mechanics approach, the time-average (2.2) is replaced
by the ensemble-average
R 3N 3N
d p d q M(p, q) ρ(p, q)
hM i = R (2.3)
d3N p d3N q ρ(p, q)
as the quantity of interest. Here the density-function
ρ(p, q)d3N p d3N q (2.4)
describes the probability that the system is in the volume d3N p d3N q around
the point (p, q) in phase-space Γ.
While the phase-space density ρ(p, q, t) is a time-dependent function in
general, its total derivative dρ/dt vanishes by Liouville’s theorem: A change
of density ∂ρ/∂t in the volume ω ⊂ Γ can only be due to a flow of density
through the boundary ∂ω,
Z Z
∂ρ
dω = − d~σ · (ṗ, q̇)ρ
ω⊂Γ ∂t
Z∂ω
Gauss
= − dω (∇p , ∇q ) · [(ṗ, q̇)ρ], ∀ω
ω
∂ρ
→ = −(∇p , ∇q ) · [(ṗ, q̇)ρ] (2.5)
∂t
3N
X ∂(ṗi ρ) ∂(q̇i ρ)
= − −
∂pi ∂qi
i=1
3N
X ∂ρ ∂H ∂ρ ∂H
= − = {H, ρ}, (2.6)
∂pi ∂qi ∂qi ∂pi
i=1
where we have used (above, {f, g} denotes the Poisson bracket of f and g)
∂ ṗi ∂ q̇i
= 0, = 0,
∂pi ∂qi
∂H ∂H
= −ṗi , = q̇i .
∂qi ∂pi
According to Eq. (2.5), the material derivative of ρ
dρ
[∂t + (ṗ, q̇) · (∇p , ∇q )]ρ = =0 (2.7)
dt
47
vanishes; this is Liouville’s theorem. However, note that the partial deriva-
tive ∂t ρ 6= 0 does not vanish in general2 .
In statistical mechanics, we are not interested in an arbitrary state of
the system but specifically in the thermodynamic equilibrium state. For an
isolated system in equilibrium with the volume V , the particle number N ,
and the energy E fixed, the principle of equal a-priori probabilities provides
us with the density-function ρ(p, q) in the form
(
const, E < H(p, q) < E + ∆,
ρ(p, q) = (2.8)
0, otherwise,
i.e., ρ is uniform over all states (p, q) which are compatible with the con-
straints. As ρ(p, q) = ρ(H(p, q)), we immediately have ∂t ρ = 0 from (2.6),
i.e., ρ is time-independent in thermodynamic equilibrium.
A set of systems prepared with identical macroscopic parameters (for
example V , N , E fixed) is called an ensemble. Such a (thermodynamically
equilibrated) ensemble is associated with a corresponding density function
ρ(p, q); the specific ensemble characterized by fixed V , N , E and described
by the density ρ(p, q) as given in (2.8) is called the microcanonical ensemble.
The fundamental idea of statistical mechanics then is that the time-
average and ensemble-average of an observable M coincide,
T
M = hMi, (2.9)
time-average = ensemble-average.
where wi are the probabilities to find the system in the quantum state
|Ψi i. The average then is the combination of a classical average with
probabilities wi and a quantum average over the wave functions or
amplitudes Ψi .
In this lecture, we will discuss the aspects 1 and 2 in detail, while the
non-equilibrium aspects described in 3 will be left for another advanced
course.
Chapter 3
Classical statistical
mechanics
49
50 CHAPTER 3. CLASSICAL STATISTICAL MECHANICS
The constant h has the dimension of an action and vanishes from the ex-
pectation hMi in Eq. (3.1). The factor N ! resolves the Gibbs paradox aris-
ing in the context of the mixing of indistinguishable particles. The factor
N ! makes sure that all permutations of the state (pπ(1) , . . . , pπ(N ) , qπ(1) , . . . ,
qπ(N ) ), π ∈ S N , are counted only once, resulting in the correct ‘Boltzmann-
counting’. That microscopic particles are truly indistinguishable is a postu-
late of quantum statistics and carries over (in its principle) to the classical
description.
We define the (dimensionless) volume
Z
Γ(E) ≡ d3N p d3N q ρ(p, q)
d3N p d3N q
Z
= 3N
(3.3)
E<H(p,q)<E+∆ N ! h
Making use of the volume Γ defined in Eq. (3.3), we introduce the expression
S(E, V, N ),
S(E, V, N ) = kB log Γ(E), (3.6)
where kB = 1.380·10−16 erg/K = (1/11605) eV/K and 1 eV = 1.602·10−19 J.
The expression S is (see the proof below)
− extensive in E, V, N
− maximal in an isolated system, i.e., (3.7)
S(E, V, N ) > S(E, V, N ; constrained).
3.1. MICROCANONICAL ENSEMBLE 51
Given the statements in (3.7), we can identify S(E, V, N ) with the system’s
entropy. This provides us with a microscopic interpretation of entropy: it
counts the number of available states in Γ-space and its maximization directs
the isolated system towards a homogeneous distribution of states compatible
with the fixed variables E, V, N .
E E2
E p2
p1 E 1+ Δ
E1
0
Δ << E
E 2 +Δ 0
q1 q2
E/∆
X
Γ(E) = Γ(E1i )Γ(E − E1i ) and thus
i
E/∆
X
S(E, V, N ) = kB log Γ(E1i )Γ(E − E1i ). (3.11)
i
52 CHAPTER 3. CLASSICAL STATISTICAL MECHANICS
We then have
E
Γ(Ē1 )Γ(Ē2 ) ≤ Γ(E) ≤ Γ(Ē1 )Γ(Ē2 ). (3.12)
∆
As log Γ(Ēi ) ∝ Ni , log E/∆ ∼ log N , we immediately find that
Let us next investigate how we split the system when choosing the divi-
sion E = Ē1 + Ē2 . As the energies Ē1 and Ē2 are the main contributors to
the entropy, the variation
δΓ1 δΓ2
= − , and using δE1 = −δE2
Γ1 Γ2
∂ ∂
log Γ1 = log Γ2
∂E1 ∂E2
1 ∂S1 ∂S2 1
→ ≡ = ≡ . (3.15)
T1 ∂E1 Ē1 ∂E2 Ē2 T2
1 ∂S
= . (3.16)
T ∂E V,N
Combining the results (3.13) and (3.15), we find that separating our
system into two subsystems while keeping them in equilibrium, our function
S is additive and hence extensive, up to small log corrections. Note that the
result in (3.13) also implies that the most likely (equlibrium) distribution of
energies is E = Ē1 + Ē2 : nearly all configurations have the optimal energy
distribution, i.e., the one belonging to the largest phase space volume.
other equilibrium parameters) and that will entail a rise in entropy. Hence,
a division of the system into two parts that remain in equilibrum implies
extensiveness, S = S1 + S2 , while division into two parts out of equilibrium
implies a rise in entropy, S > S1 + S2 . 2
Repeating the above steps and relaxing the constraints in V and N leads
us to the additional equilibrium parameters p (pressure) and µ (chemical
potential),
∂S ∂S
p=T , µ = −T . (3.17)
∂V E,N ∂N E,V
With the obvious identification of U ≡ E and using the definition of entropy
through Eq. (3.6), we can identify S expressed as a function of the extensive
variables E, V, N as a thermodynamic potential
S(U ≡ E, V, N ) = kB log Γ(E). (3.18)
Hence, the calculation of the phase-space volume Γ(E) (or Σ(E), ω(E))
provides us with the entire thermodynamics of an isolated system that is
characterised by the extensive variables E, V, N . Note that log Γ ∼ log Σ ∼
log ω. We can go over to open systems by using the Legendre transformations
F = U − T S, G = U − T S + pV .
S(E, V, N ) = kB log Σ
h V 4πmE 3/2 5 i
= N kB log + log + . (3.24)
N 3N h2 2
Here, the prefactor is extensive ∝ N and the terms in parenthesis are inten-
sive ∝ V /N (from N !) and ∝ E/N (from Γ(3N/2 + 1)). The internal energy
3h2 N 2/3 h2 S 5i
U (S, V, N ) = N exp − (3.25)
4πm V 3 N kB 3
is obtained by inversion and the definition of the temperature T via the
derivative of U , providing us with the caloric equation of state,
∂U 2 U 3
T = = , U = N kB T = E. (3.26)
∂S V 3 N kB 2
Making use of the specific heat
∂U 3
CV = = N kB , (3.27)
∂T 2
we recover the thermodynamic expression for S(V, T )
V
S = N kB log + CV log T + const. (3.28)
N
1 z=n
The integer and half-integer values of the Γ-function derive from Γ(z + 1) = zΓ(z) =
√
n! and Γ(1/2) = π and thus
√ √
π π ππ 4π 3
S1 = √ R = 2R, S 2 = R2 , S 3 = √ R3 = R .
π/2 1 (3/4) π 3
Exercise: Calculate the formula for the surface and the volume of an n-dimensional sphere
by using the integration of a Gaussian in Cartesian and polar coordinates. Show that the
volume of a high-dimensional R ∞sphere concentrates itself on the surface. Calculate the
Stirling formula from N ! = 0 dx eS(x) , S(x) = N log x − x and evaluate the integral by
the method of steepest descent.
3.1. MICROCANONICAL ENSEMBLE 55
Equipartition theorem
XD ∂H ∂H E ai ,bi 6=0 X
pi + qi = h2Hi = 2 kB T. (3.33)
∂pi ∂qi
i i
As a result, every kinetic and potential degree of freedom carries the energy
kB T /2.
Proof:
we find
D ∂H E Z Z
1 ∂H ∆ ∂ ∂H
xi = dp dq xi = dp dq xi
∂xj Γ E<H<E+∆ ∂xj Γ ∂E H<E ∂xj
Z Z
δij ∂ δij δij
= dp dq (E − H) = dp dq = Σ
ω ∂E H<E ω H<E ω
1 kB
= δij = δij = δij kB T.
∂E log Σ ∂E S
3N D 3N
X ∂H E DX E
− qi = qi ṗi = −3N kB T. (3.34)
∂qi
i=1
| i=1{z }
virial from classical mechanics
3
The virial theorem from mechanics tells us that the average value of the total kinetic
energy hT P
i in a bound system is linked to the total potential energy hV i, more precisely,
2hT i = − N ~ ri i, with F
~i the force acting on the particles.
i hFi · ~
3.2. CANONICAL ENSEMBLE 57
The quantity F is extensive5 and relates to the internal energy U and the
entropy S via (proof follows)
F = U − T S. (3.40)
The quantity F (T, V, N ) thus obeys the differential equation (3.41) which
is the one defining the energy potential U (S, V, N ) from the free energy
potential F (T, V, N ) via the Legendre transformation; using −S = ∂T F , we
have U = F + T S = F − T ∂T F .
For a canonical ensemble, we then can recover all of the system’s ther-
modynamic behaviour by calculating its partition function
e−β H(p,q)
Z
ZN (T, V ) = d3N p d3N q
N ! h3N
∂F
S = − ,
∂T V,N
U = F + T S, caloric equation,
∂F
p = − , thermodynamic equation,
∂V T,N
∂F
µ = , thermodynamic equation.
∂N T,V
Using the equations of state for the pressure ∂V F = −p and the chemical
potential ∂N F = µ that characterize the reservoir T, V2 , N2 , it follows from
(3.42) that
ZN2 (T, V2 )
= e−β[F (T,V −V1 ,N −N1 )−F (T,V,N )]
ZN (T, V )
≈ e−β[pV1 −µN1 ] . (3.48)
and using the density function (3.51) and its normalization (3.47), we find
the grand potential 7
−Ω = pV = kB T log Z(T, V, z). (3.53)
The grand canonical partition function directly gives us the thermal equation
of state p = p(T, V, µ). We find the caloric equation of state by calculating
the average energy8
XZ
U = hHi = d3N p d3N q H(p, q)ρ(p, q, N )
N
d3N p d3N q
X Z
= e|−βpV zN H(p, q)e−βH(p,q) ,
{z } N ! h3N
1/Z N | {z }
−∂β ZN (T,V )
∂
U = − log Z(β, V, z) . (3.54)
∂β V,z
8
Note that here, we view z as an independent variable for µ and not as a function of
β and µ, i.e., ∂β z = 0.
3.4. FLUCTUATIONS* 61
3.4 Fluctuations*
In the canonical ensemble the energy E is not fixed, and similar is the case for
N in the grand canonical ensemble. Let us then determine the fluctuations in
energy (E) in the canonical ensemble and of the particle number (N ) in the
grand canonical ensemble. We show that √ the fluctuations in the extensive
quantities E and N are small, of order N ,
hH 2 i − hHi2 = kB T 2 CV , hHi ∼ N, CV ∼ N,
N2 1 ∂V
hN 2 i − hN i2 = kB T κT , κT = − . (3.57)
V V ∂p T
d3N p d3N q
Z
1
U = hHi = He−βH
ZN N ! h3N
d3N p d3N q
Z
[U − H]e−β(H−F ) = 0.
N ! h3N
d3N p d3N q
Z
∂U h i
+ [U − H] F + β∂β F −H e−β(H−F ) = 0,
∂β N ! h3N | {z }
=hHi=U
from which we obtain the desired result for the energy fluctuations
∂U
h(U − H)2 i = − . (3.58)
∂β
∂U ∂U
= −kB T 2 = −kB T 2 CV
∂β ∂T
62 CHAPTER 3. CLASSICAL STATISTICAL MECHANICS
1 ∂2F
F(E) ' F(Ē) + (E − Ē)2 ,
2 ∂E 2 Ē
∂F ∂S 1
= 1−T = 0 → ∂E S|E=Ē = ,
∂E Ē ∂E Ē T
we write ∂E S = T̃ (E)−1 , T̃ (Ē) = T,
∂2F ∂2S
= −T (3.61)
∂E 2
Ē ∂E 2 Ē
∂ ∂S 1 ∂ T̃ 1
= −T = = .
∂E ∂E T ∂E T CV
Ē=U 1
→ F(E) ' [U − T S(U )] + (E − U )2 . (3.62)
2T CV
9
For the ideal gas, we can write F(E) = N [e − (3/2)kB T log e] + rest with e = E/N , see
Eq. (3.24). This spectral free energy has a minimum at ē = (3/2)kB T as trivially obtained
from ∂E F(E) = [1 − (3/2)kB T /e]ē = 0.
3.4. FLUCTUATIONS* 63
Z ∞ h (E − U )2 i
ZN ≈ exp −β(U − T S(U )) dE exp − (3.63)
| {z } 0 2kB T 2 CV
F̄ √
| {z }
2πkB T 2 CV
p
≈ 2πkB T 2 CV e−β(U −T S) = e−βF .
U E
1
F = U − T Smicrocan − kB T log CV ,
2
1
→ Scan = Smicrocan + log CV . (3.64)
2 | {z }
O(log N )
0 U E
The particle number fluctuations ∆N appear from (3.52) (also see (3.55))
and are given by
∂ ∂
hN 2 i − hN i2 = z z log Z(T, V, z)
∂z ∂z | {z }
pV /kT
∂ 1 ∂
∂z
= βz ∂µ V 1 ∂2p ∂2p
= = V k B T . (3.65)
kB T β 2 ∂µ2 ∂µ2
We rewrite ∂µ2 p by using the homogeneity of F
With
∂F ∂f
p = − =− ,
∂V ∂v
∂F ∂f
µ = =f −v ,
∂N ∂v
we find that ∂µ2 p is proportional to the compressibility κT ,10
∂2p ∂ ∂p
=
∂µ2 ∂µ ∂µ
1 ∂v 1 1 1 κT
= − 2 = 3 2 =− 3 = 2
v ∂µ v ∂v f v ∂v p v
and we obtain the result (3.57b)
Note that both fluctuations ∆E and ∆N are given by their respective linear
response coefficients11 CV and κT ; this property is related to the fluctuation–
dissipation theorem, see Landau-Lifschitz.
Since the fluctuations in N are small, we can identify N = hN i = N̄ in
(3.52) as the dominant term and obtain
The chemical potential in the grand canonical ensemble plays the role of
a Lagrange multiplier H → H − µN that guarantees the conservation
of particle numbers. If particles may spontaneously be created then the
difference between particle and antiparticle numbers is usually conserved
NT − NAT = const. Correspondingly the expression µN is replaced by
µ(NT − NAT ) in the formula for the grand canonical partition function
X
Z= ZNT ZNAT eβµ(NT −NAT ) . (3.71)
NT ,NAT
W gas
liquid & gas
liquid
κT = κT <
8
V/ v g V /v li N
Quantum statistical
mechanics
HΦn = En Φn . (4.2)
67
68 CHAPTER 4. QUANTUM STATISTICAL MECHANICS
where cn are the amplitudes of the reservoir1 and Φn are the states of the
subsystem (we choose a time-independent basis {Φn } of the subsystem). As
we do not know the state of the reservoir, we replace the coefficients c∗n cm
by expectation values hcn |cm ires ; similarly |cn |2 becomes hcn |cn ires . The
quantity measured in an experiment
P is the time average of M. With a time-
independent norm hΨ|Ψi = n hcn |cn ires in the isolated system = system +
reservoir, we have
P T
T m,n hcn |cm ires hΦn |M|Φm i
hMi = . (4.6)
hΨ|Ψi
T
As a first step, we consider the non-diagonal elements hcn |cm ires , n 6= m.
The reservoir dynamics generates amplitudes cn with random phases, such
that the time-averages of matrix elements of differing amplitudes vanish
T
hcn |cm ires = 0 if n 6= m. (4.7)
This is the postulate of random phases. Note that inelastic collisions be-
tween particles are sufficient to randomize the phases.2 As a result of phase
T
randomization, only the diagonal weights hcn |cn ires remain in the average
(4.6). In a next step, we have to make a statistical assumption about the
system in equilibrium, as is usually done with the assumption of molecular
chaos in kinetic theory or the equal a priori probabilities in phase-space
Γ for the microcanonical ensemble (at fixed extensive variables E, V , N )
in classical statistical mechanics. The latter is needed to replace the time-
average by the ensemble-average and brings us to the ergodic hypothesis:
The trajectories of particles explore equally often all parts of phase-space
that are compatible with fixed E, V , and N . Thus, in the present case,
we also define the microcanonical ensemble by giving all states Φn that are
compatible with the constraints of fixed energy E, volume V , and particle
1
Here the reservoir acts as a generator of random phases as opposed to one exchanging
energy or particles. The latter types of reservoirs will be added in the context of the
canonical and grand canonical ensembles.
2
The assumption (4.7) breaks down when the reservoir does not properly randomize
phases. This can happen in small systems when the phase is preserved over distances
that are of similar length or longer than the sample. Also, phase information may be
kept in a small reservoir as is often discussed in the context of quantum thermodynamics.
The important element entering the discussion then is the phenomenon of entanglement
between the system and the reservoir, preventing the system from being factorized.
4.1. EXPECTATION VALUES IN QUANTUM STATISTICS 69
with weights (
1, E < En < E + ∆,
wn = (4.10)
0, otherwise.
If wn = δn0 the state is called pure; if more than one weight wn isP
non-zero it
is mixed. Usually, the weights wn are normalized according to n wn = 1.
Here, we consider the normalisation explicitly and write
T Tr(ρM) X
hMi = , Tr(A) = Trace(A) = hΦn |A|Φn i, (4.12)
Tr(ρ) n
for the expectation values. The equivalence of (4.9) and (4.12) is trivial.
The time evolution of ρ is given by the quantum mechanical version of
Liouville’s theorem, the Neumann differential equation
[H, ρ] = 0 (4.14)
Using the concept of density matrices ρ from the beginning, one may
proceed in a more straightforward way. Consider a system that is iso-
lated but still coupled to the rest of the world (a reservoir) with the total
Hilbert space Hsys ⊗ Hres . Note that this reservoir neither exchanges energy
nor particles with our isolated system, hence E, V , N are fixed (alterna-
tively, we look at a subsystem in a larger system). In quantum mechanics,
the measurement of a system observable M provides us with the expec-
tation value hMi = Tr(ρM). Since we do not know anything about the
rest of the world (or reservoir), we trace it out, hMi = Trsys Trres (ρM) =
Trsys [Trres (ρ)M] ≡ Trsys (ρsys M), since M only acts in Hsys . In tracing ρ
over the reservoir, we have eliminated the entanglement with the rest of the
world, what corresponds to the phase randomization above. Requiring ρsys
to describe a time-independent equilibrium state, we have ∂t ρsys = 0 and
As a result, ρsys is diagonal in the eigenstates |Φn i of
hence [Hsys , ρsys ] = 0. P
Hsys and hence ρsys = i wi |Φn ihΦn | with real coefficients wi (since ρ = ρ† ).
Finally, the a priori assumption of equal probabilities wi = 1/Γ(E) provides
the final result for the properly normalized (i.e., Tr(ρsys ) = 1, cf. Eq. (4.18))
density matrix of the microcanonical ensemble, ρsys = 1l/Γ(E).
With S = kB log Γ the entropy counts the number of states that are avail-
able. According to ρ ∼ e−H/T only the ground state is available when T = 0
and thus S = 0. If the ground state is g ∼ N -times degenerate, we still have
s = S/N ∼ (log N )/N → 0, and the entropy per particle vanishes. A more
relevant problem is the continuity of the spectrum—the argument above
then is still correct but mainly academic as the ground state is only realized
for T < ∆E ∼ ~2 /mL2 ∼ 5 · 10−15 K (m = mnukleon , L = 1 cm). Thus the
density of states ω(E) is relevant in the low temperature phases of matter.
Usual matter transform into one of the following phases for low temperatures
T → 0: crystals (nearly all) or superfluids/superconductors (3 He, 4 He, met-
als). The density of states in these systems then has a gap (e.g., fermionic
superfluids and superconductors) or vanishes as ω(E) ∼ E d−1 for E → 0
(Goldstone modes) and thus S → 0 for T → 0 (see Section 7.2). The notoric
exceptions are frustrated systems, glasses, two-level systems, etc., showing
a residual entropy as T → 0.
Chapter 5
The classical partition function for a free particle follows from (3.38)
Z 3 3
d pd q −p2 /2mkB T 2πmk T 3/2
B
Z1 = e = V
h3 h2
Z
where we have used d3 q = V,
Z ∞
2 √
and dx e−x /2σ = 2πσ. (5.1)
−∞
73
74 CHAPTER 5. APPLICATIONS AND TECHNIQUES
where |~
p i is a basis and using
2
H = −~2 ∇ ~ /2m, (5.3)
~ √
p i = eik·~x / V ,
h~x |~
2π
p~ = ~~k, ~k = ~n, ni ∈ Z,
L
and V = L3 the quantization volume
X Z L 3 Z
3
= d n= d3 k
2π
~
n
d3 p −p2 /2mkB T
Z mk T 3/2
−p2 /2mkB T
X B
→ Z1 = e =V e = V . (5.4)
(2π~)3 2π~2
~
n
The quantum mechanical calculation gives the classical result and fixes
the quantity h in phase-space: In accord with the uncertainty principle,
we should divide µ-space (the Γ space for one particle) into boxes of size
d3 p d3 q = (2π~)3 and count one state per box,
→ class. Z1 = q.m. Z1 if h = 2π~,
~ = the reduced Planck’s constant.
With the thermal de Broglie wavelength
s
2π~2
λ= (5.5)
mkB T
the partition function is
V L 3
Z1 = = . (5.6)
λ3 λ
The thermal de Broglie wavelength is the wavelength of a particle with
energy E = πkB T ,
~2 k 2 ~2 (2π)2 2π~2 1/2
E= = = πkB T → λ = . (5.7)
2m 2mλ2 mkB T
∂ h ~2 i
− p |e−βH |~q i = −
h~ p |e−βH |~q i.
∆ + V h~ (5.12)
∂β 2m
1
For independent particles Z1 → Z1N and the factor 1/N ! takes care about the Gibbs
paradoxon.
76 CHAPTER 5. APPLICATIONS AND TECHNIQUES
1 ~ p2
∂β S 1 (∇S0 )2 + V =
= +V
2m 2m
→ S1 = βH(~ p, ~q ),
1 ~ ~ 1 ) − 1 ∆S0 = p~ β ∇V ~
∂β S 2 = (∇S0 )(∇S
m 2m m
β2 ~
→ S2 = p~ · ∇V
2m
1 ~ 1~ ~ 1
∂β S 3 = (∇S1 )2 + ∇S 0 · ∇S2 − ∆S1
2m m 2m
β2 ~ 2 β2 ~ )2 V − β ∆V,
= (∇V ) + 2
p·∇
(~
2m 2m 2m
β3 ~ 2 β3 2
~ )2 V − β ∆V,
→ S3 = (∇V ) + (~
p · ∇
6m 6m2 4m
plus further orders Sn≥4 ∝ ~ in the exponent.2 The partition function
3
For the classical expectation values one easily finds (we only point out some
critical steps)
Z
hS2 icl ∼ d3 p p~ · (. . .) = 0, (integral odd in pi )
Z 3 3
1 d p d q −βH β3
hS3 icl = cl 3
e [. . . + pi pj ∂i ∂j V ] → h∆V i,
Z1 ~ R
6m2 |{z}
pi pj e−βH →δij mkB T
d3 p d3 q −βH β 4
Z
1
hS22 icl = e pi pj (∂i V )(∂j V ) → h(∇V )2 i,
Z1cl ~3 R
4m 2 |{z}
pi pj e−βH →δij mkB T
and using
Z Z Z Z
−βH −βH ~ −βH ~
e ~
∆V = − ∇e ~
∇V = β ∇V e ∇V = e−βH β(∇V
~ )2 ,
with the quantum correction of order (λ/l)2 V with l ≈ V /|∇V~ | the typical
3
scale of variations in the potential V . Indeed, the condition λ/l 1 of
smooth potential is the standard requirement for the validity of the quasi-
classical (WKB) approximation. Hence, we see that the configurational
integral is modified due to an additional quantum smearing of the potential.
For the free energy F1 = −kB T ln Z1 , we obtain the quantum corrected
version
~2 β 2 ~ 2
F1 = F1cl + h(∇V ) icl + O(~3 ),
24m
where we have again used h∆V i = βh(∇V ~ )2 i. Note that we could perform
2
the expansion in ~ systematically. The expansion in β (for example via
3
We assume that there is only one scale in V (q) such that ∂q2 V /V ∼ (∂q V /V )2 ∼ 1/`2 .
78 CHAPTER 5. APPLICATIONS AND TECHNIQUES
HΦn = En Φn
Elastic manifolds
with the mass density ρ and the elasticity C. This Hamiltonian derives from
the Lagrangian
Z h ρ ∂u 2 C i
L = dd ~q ~ 2 − V [u(~q )]
− (∇u) (5.16)
2 ∂t 2
∂L ∂u
Π = =ρ ,
∂u,t ∂t
h = Πu̇ − `
comprising kinetic energy density, elastic energy density, and the potential.4
The displacement field ~u(~q ) can be vectorial in general. For d = 1 at V = 0
(V 6= 0) (5.16) describes a free (bound) string.5
4
The variation produces the equations of motion in Lagrangian or Hamiltonian form,
∂` ∂` ∂`
0 = ∂t + ∂i − ,
∂(∂t u) ∂(∂i u) ∂u
∂h ∂h
−Π̇ = − ∂i ,
∂u ∂(∂i u)
∂h
u̇ = .
∂Π
5
Note the symmetry in the dynamic and elastic
p terms for the “Euclidean” Lagrangian
LE = T + V (compare to (5.60)): t → t0 = t c/ρ gives a d + 1-dimensional isotropic
problem.
5.2. GASES 79
5.2 Gases
The topic of quantum gases has attracted huge interest during recent years.
The ability to laser- and evaporatively cool atoms to amazingly low tem-
peratures (µK and nK) using optical traps and the development of suitable
(magneto-optic) traps allowed Eric Cornell and Carl Wiemann at JILA and
Wolfgang Ketterle, Kendall David, and Marc-Oliver Mewes at MIT to pro-
duce a Bose-Einstein condensate in 1995; in 2001 Cornell, Wiemann, and
Ketterle received the physics Nobel prize for this achievement. At the be-
ginning of the millenium, the cooling of fermionic atoms into a quantum-
degenerate state with kB T εF followed. Today quantum optics is an ever
blooming field. The focus is on interesting systems described by compli-
cated Hamiltonians that generate non-trivial physics. The quantum gases
are then used as ‘quantum simulators’ for such Hamiltonians. Elements
of such Hamiltonians include non-trivial periodic lattices (→ simulations
of graphene, topological insulators assisted by optical lattices), interactions
(tunable by Feshbach resonances), interactions mediated by gauge fields
(generated by structured laser fields), disorder, Bose-Fermi mixtures, mul-
ticomponent gases, etc.
Note the difference between quantum gases and strongly correlated quan-
tum gases, i.e., gases with strong interactions, and quantum liquids that
commonly interact strongly and have a high density. Typical examples
of quantum liquids are fermionic 3 He and bosonic 4 He. The subject of
fermionic quantum liquids brings one to the Landau theory of Fermi-liquids
(encompassing the electron gas/liquid in a metal), that is renormalized due
to the strong interactions but otherwise behaves phenomenologically like a
fermionic quantum gas.
In what follows, we first discuss the statistics of classical and quantum
gases before we study Fermi and Bose gases more thouroughly in the follow-
ing two chapters, respectively.
For later comparison, we first start with the classical or so called “Boltz-
mann” gas. We calculate6 the partition function in the canonical ensemble
6
We discretize phase-space and work with discrete momenta p
~i .
80 CHAPTER 5. APPLICATIONS AND TECHNIQUES
P
−β np~i εpi
X
ZN = e p
~i
· gc , (5.17)
|{z}
{np~i }N configurations
X
np~i = # particles with momentum p~i , np~i = N,
p
~i
εp i = p2i /2m.
The correct Boltzmann counting accounts for the number of possible ways
to distribute the N particles among the momenta p~0 , . . . , p~i , . . .,7
N!
g̃c = .
np~0 !np~1 ! . . . np~i ! . . .
g̃c → gc = g̃c /N !.
and then replace the sum over momenta by integrals, counting one state per
phase-space volume h3 → (2π~)3 (alternatively replace d3 n → (L/2π)3 d3 k →
V /(2π~)3 d3 p)
d3 p p2
Z
X
−βεpi
V
e =V exp − = 3,
(2π~)3 2mkB T λ
p
~i
where s
2π~2
λ= (5.19)
mkB T
7
In the correct ‘Boltzmann counting’, we first treat count configurations with the num-
ber g̃c > 1 of their appearance and then account for the ‘indistinguishability’ of particles
with the Gibbs factor 1/N !. In quantum mechanics, indistinguishability is accounted for
by the fermionic or bosonic statistics.
8
The multinomial theorem says that
!
X k
(x1 + x2 + · · · + xn )k = xk1 xk2 · · · xknn ,
k1 , . . . , k n 1 2
k1 +···+kn =k
where !
k k!
= .
k1 , . . . , k n k1 ! · · · kn !
5.2. GASES 81
1 V N
ZN = , (5.20)
N ! λ3
V h
3
i
log ZN = N log + N = N − log(nλ ) + 1 ,
N λ3 h i
F = −kB T log ZN = N kB T log(nλ3 ) − 1 ,
∂F
S = − = −F/T + (3/2)N kB = (−F + U )/T,
∂T V
∂F
p = − = N kB T /V,
∂V T
with the gas parameter δ ≡ nλ3 . For an ideal classical gas δ 1, i.e.,
the gas is dilute (small density n) or hot (large temperature T , small λ).
Using U = (3/2)N kB T , the above results produce the Sackur-Tedrode result
(3.24)
for the entropywith E replaced by T and h by 2π~, S(T, V, N ) =
N kB − log(nλ3 ) + 5/2 .
The quantum gases are easiest to describe in the grand canonical ensemble
(this removes the constraint on the particle number N , i.e., all values of N
are admitted)
∞
X
Z(V, T, z) = z N ZN (V, T )
N =0
X∞ X X
= z N exp −β np~ εp
N =0 {np~ }N p
~
| {z }
np
~
Q −βεp
p
~ ze
XX X h np~ np~ i
= ... ... ze−βεp0 0 . . . ze−βεpk k ...
np~0 np~1 np~k
YhX n i
= ze−βεp . (5.21)
p
~ n
we obtain (
p = (kB T /λ3 )f5/2 (z),
fermions (5.26)
N/V = n = (1/λ3 )f3/2 (z),
where the functions f5/2 and f3/2 derive from [we use the expansion log(1 +
2
s) = s − s2 /2 + s3 /3 − s4 /4 + · · · = ` (−1)`+1 s` /` and s = ze−x ]
P
Z ∞
4 2
f5/2 (z) = √ dx x2 log(1 + ze−x ) (5.27)
π 0
Z ∞ ∞ ∞
4 X (−1)`+1 z ` −`x2 X (−1)`+1 z `
= √ dx x2 e = ,
π 0
`=1
`
`=1
`5/2
∞
X (−1)`+1 z `
∂
f3/2 (z) = z f5/2 (z) = ,
∂z
`=1
`3/2
5.2. GASES 83
with
∞ ∞
z`
Z
4 2 −x2
X
g5/2 (z) = − √ dx x log(1 − ze )= , (5.29)
π 0 `=1
`5/2
X z` ∞
∂
g3/2 (z) = z g5/2 (z) = .
∂z `3/2 `=1
For z → 1, the term (1/V ) z/(1 − z) = hn0 i/V will single out a fraction of
the total density as a k-space condensate, i.e., hn0 i ∝ N ; we will discuss
in detail later how the occupation number hn0 i can become macroscopically
occupied. The consequence is the appearance of quantum mechanical effects
on a macroscopic scale: MQP = macroscopic quantum phenomena. The
internal energy U of the quantum gas is obtained from (5.23),
1 X N X X P
U (V, T, z) = z np~ εp e−β p~ np~ εp
Z
N {np~ }N p
~
(
∂ ∂ h V f5/2 (z) i
= − log Z(V, T, z) = −
∂β ∂β λ3 g5/2 (z)
(
3 V kB T f5/2 (z), fermions,
= (5.30)
2 λ3 g5/2 (z), bosons,
3
= p(V, T, z)V.
2
In the last expression, valid for the pure Bose gas, we have ignored the con-
densation terms for the bosons. Free Fermi and Bose gases have a multitude
of interesting properties and applications which we will come back to later.
9
use again the expansion of the logarithm, log(1 − s) = −(s + s2 /2 + s3 /3 + s4 /4 +
WeP
· · · ) = ` s` /`.
84 CHAPTER 5. APPLICATIONS AND TECHNIQUES
H0 |Φn i = En |Φn i,
N
X p2i 2π
En = , p~i = ~~ni , ~ni ∈ Z3 ,
2m L
i=1
1 X
h~q1 , ~q2 , . . . ~qN |Φn i = √ (±1)π φp~1 (~qπ1 ) . . . φp~N (~qπN )
N ! π∈S
N
ei~p·~q/~
φp~ (~q ) = √ .
V
VN
X Z
→ d3N p
n
N ! (2π~)3N
VN
Z
ZN = d3N p d3N q |Φn (~q1 , . . . , ~qN )|2 e−βEn . (5.32)
N !(2π~)3N
The factor 1/N ! takes into account that permutations of the vector p~i in
|Φn i does not generate any new states. We simplify the expression for the
probability |Φn |2 according to
1 X 0
|Φn |2 = (±1)π+π · φp∗~1 (~qπ1 )φp~1 (~qπ10 ) . . . φp∗~N (~qπN )φp~N (~qπN
0 )
N! 0
π,π
1 X
= (±)π ei~p1 ·(~q1 −~qπ1 )/~ . . . ,
VN π
where we have used that every permutation π 0 under the integral d3N p
R
λ2
P
For T → ∞, → 0, the only term in the sum π that survives is the
identity 1lPand thus we obtain the classical result (5.20). For T < ∞ the
terms in π become smaller the more particles are permuted by π. As long
as h(~qi − ~qj )2 i > λ2 (i.e., at low densities), we can expand the sum in the
number of permutations and obtain
X X h 2π i
(±1)π [. . .] ≈ 1 ± exp − 2 (~qi − ~qj )2 + 3 particle − terms.
π
λ
i<j
VS
Fig. 5.1: Statistical inter-
action for bosons (B, attrac-
F tive) and fermions (F, repul-
sive with a hard-core).
0 q /λ
1
B
−k B T log 2
Comparing (5.35) with the classical result (5.10), we see that the quantum
statistics of the particles expresses itself as an effective potential under a
configurational integral, the so called statistical interaction
2
V S (q) = −kB T log 1 ± e−2π(q/λ) . (5.36)
86 CHAPTER 5. APPLICATIONS AND TECHNIQUES
This interaction is short ranged (λ) and depends on the particle statistics,
attractive for bosons and repulsive for fermions, see Fig. 5.1.
1 r0
k cot δ0 ≈ − + k 2 + . . . . (5.38)
a 2
The scattering length is a > 0 for a repulsive interaction and a < 0 for
an attractive interaction. The parameter r0 gives the effective range of the
r Ψ( r )
r Ψ (r )
V
0 a r a 0 r
V
Fig. 5.2: Shape of the wave function Ψ(r) for a repulsive (left) and an
attractive (right) interaction V (r).
that of hard spheres of radius a > 0. The s-wave scattering of hard spheres
can conveniently described by the pseudopotential
a 4πa~2 ∂
V`=0 (~r ) = δ(~r ) r (5.39)
m ∂r
i.e., the two Hamiltonians
p2
H = + Vhard sphere, R=a ,
2m
p2 a
H = + V`=0 (~r ), (5.40)
2m
generate the same scattering solutions Ψ`=0,k (~r ) in the exterior space r >
a for ` = 0, k → 0; correspondingly they generate the same eigenvalues
E`=0 (k).
a (~
A further simplification arises when the pseudopotential V`=0 r ) is treated
perturbatively. Within a perturbative analysis, one uses plane waves as the
unperturbed solutions. As these solutions, say ϕ(~r ), are regular at the ori-
gin ~r = 0, we have δ(~r )(∂r rϕ) = δ(~r )ϕ(0) and one may replace11 (5.39) by
a (~
V`=0 r ) = (4πa~2 /m)δ(~r ). To treat repulsively interacting quantum gases
(either fermions or bosons) in the limit of small temperatures T → 0 the
following strategy is suitable:
X p2 4πa~2 X
i
H= + δ(~qi − ~qj ) (5.41)
2m m
i i<j
for the real space representation of the density matrix. The idea then is that,
while the calculation of ρ(x, x0 ; τ ) is difficult, the determination of ρ(x, x0 ; ε)
is simplified by the presence of the small parameter ε. This will become
apparant when we discuss the particle in a potential V (x) further down—
before going into this, however, we discuss the simplest situation, which is
free particle.
For a free particle, the imaginary time propagator over an interval ε has
the form m 1/2 h m i
ρ0 (x, x0 ; ε) = exp − (x − x0 )2 , (5.50)
2π~ε 2~ε
and inserting this expression into (5.49), we arrive at the full density matrix
in the form
Z Z
0 1 dx1 dxn−1
ρ0 (x, x ; τ ) = p ... p ... p
2π~ε/m 2π~ε/m 2π~ε/m
h mε x − x 2 i h mε x − x0 2 i
n−1 1
exp − . . . exp − .
2~ ε 2~ ε
90 CHAPTER 5. APPLICATIONS AND TECHNIQUES
xk − xk−1 dx h mε i h 1Z τ m dx 2 i
→ , Π exp − 2
(. . .) → exp − dτ 0
ε dτ 2~ ~ 0 2 dτ
| {z }
SE0 [x(τ )]
with SE0 [x(τ )] the Euklidean action (and the index 0 standing for the free
particle). Furthermore, we abbreviate the product of integrals over interme-
diate positions xi as an integral over paths,
Z n−1 Z x(τ )=x
1 Y dxi
p p → D[x(τ )].
2π~ε/m i=1
2π~ε/m x(0)=x0
i~∂t G0 = H0 G0 ,
G0 (x, x0 ; 0) = δ(x − x0 ),
the action Z
i i m dx 2
S0 [x(τ )] = dt
~ ~ 2 dt
transforms to the Euclidean action by a Wick rotation
Z Z
i m dx 2 1 m dx 2 SE0 [x(τ )]
→ −idτ =− dτ =− , (5.54)
~ 2 −idτ ~ 2 dτ ~
thus reproducing the result (5.51).
In a second step, we focus on a particle in the potential V (x). We can
analyze the imaginary time Schrödinger equation (5.43) by splitting the
Hamilonian into a free part H0 and a potential part V that we treat sepa-
rately, as usually done in a perturbatve approach. 12 Given that
−~∂τ ρ0 = H0 ρ0 (5.55)
has been solved, the task is to find the solution to −~∂τ ρ = Hρ. As usual,
we rewrite the correction eH0 τ /~ ρ = ρD in the Dirac representation. The
dynamical equation for ρD is easily found,
−Hρ
z }| {
H0 τ /~ H0 τ /~
−~∂τ ρD = −e H0 ρ − e ~∂τ ρ = eH0 τ /~ V ρ
= VD ρD , where VD = eH0 τ /~ V e−H0 τ /~ .
where in the second equation we use the smallness of τ (not of V ). Using this
perturbative result, we can generalize the expression (5.50) for Rρ0 (x, x0 ; ε) to
the situation of a particle in a potential V (x). Inserting 1 = dx0 |x00 ihx00 |
into (5.56), we obtain
ρ(x, x0 ; ε) ≈ ρ0 (x, x0 ; ε) (5.57)
Z ∞ Z ε
1
− dx00 dε0 ρ0 (x, x00 ; ε − ε0 )V (x00 )ρ0 (x00 , x0 ; ε0 ).
~ −∞ 0
For small ε the main contribution to the integral dx00 comes from the
R
region
(ε − ε0 )x0 + ε0 x if x ≈ x0
x̄00 = ≈ x ≈ x0 ,
ε
which can be verified by minimizing the exponent of ρ0 · ρ0 ,
∂ (x − x00 )2 (x00 − x0 )2
+ = 0.
∂x00 ε − ε0 ε0 x̄00
hβ
imaginary time
q ( τ)
τ
q2
q
1
0
Z ~β
SE dτ h m i
= (∂τ ~q )2 + V (~q )
~ 0 ~ 2
SE0
≈ + βV (~q ), (5.62)
~
94 CHAPTER 5. APPLICATIONS AND TECHNIQUES
Z
→ρ ≈ D[~q (τ )]e−SE0 /~ e−V (~q )/kB T = ρ0 e−V (~q )/kB T ,
Z 3 Z 3
d q −V (~q )/kB T cf. (5.8) V d q −V (~q )/kB T
→ Z1 ≈ Z10 e = 3
e .
V λ V
hβ
T
8
top view
τ
hβ
τ
qy
0
q
x
Fig. 5.4: World lines in the limit β → ∞ (left) and the crossover into
classical world lines (right) for β → 0, T → ∞. The top perspective in the
center illustrates the quantum fluctuations (“smearing”) of the particle at
low temperatures.
Z
d3N q ρ(~q1 , . . . , ~qN ; ~q10 , . . . , ~qN0 ; β~)|q~i0 =~qi , (5.63)
we have to include all world lines with permuted endpoints ~qi0 = ~qP i , P ∈
S N , in the partition function. For bosons, we must add all the contributions
5.3. PATH INTEGRALS 95
disentangled entangled
hβ
Fig. 5.5: Disentangled (left) and entangled (right) world lines. Bosons at
low temperatures entangle easily and thereby attract one another, as sep-
arating entangled bosons stretches their worldlines leading to an increase
of the action; the same result has been obtained earlier on, see Eq. (5.36).
Furthermore, there is an interesting insight into the Bose condensation phe-
nomenon, as Bosons exchange their identity through the periodic boundary
conditions. Bosons entangling in a “ring”within so-called “cooperative ring
exchanges”, assume the “same identity”. Once these permutation “rings”
extend throughout the entire sample, the system undergoes condensation.
Note that there is only one disentangled state but there are N ! states that
involve entanglement.
Monte Carlo simulation one calculates and adds all configurations as those
sketched in Fig. 5.5. For fermions the problem is more complicated: The sum
P
P
p in (5.64) must be weighted by a (−1) prefactor, there is a sign problem
in the PIQMC-technique for simulating fermions. The corresponding path-
integrals involve Grassmann variables.
A cute application is the study of two-boson/two-fermion systems and
the resulting understanding of the statistical interaction in the language of
96 CHAPTER 5. APPLICATIONS AND TECHNIQUES
world lines.
F ≤ F0 + hH − H0 i0 . (5.67)
We assume that S and S0 are real and use the convexity of e−x ,
f1 f f2 f
C0 = 0,
C1 = M1 = hV (x)i,
C2 = M2 − M12 = hV 2 (x)i − hV (x)i2 ,
C3 = M3 − 3M2 M1 − 2M13 ,
...
Fermions
We summarize the results for the Fermi gas in the previous chapters:
Partition function:
Y
Z(V, T, z) = (1 + ze−βεp ). (6.1)
p
~
Equation of state:
kB T kB T
p= log Z = 3 f5/2 (z).
V λ
Particle density:
z 1
n= ∂z log Z = 3 f3/2 (z).
V λ
Occupation number:
1
hnp~ i = . (6.2)
eβ(εp −µ) +1
∞ ∞
(−1)`+1 z `
Z
4 2
X
f5/2 (z) = √ dx x2 log(1 + ze−x ) = ,
π 0 `=1
`5/2
∞ ∞
x2 (−1)`+1 z `
Z
4 X
f3/2 (z) = √ dx = .
π 0 1 + ex2 /z `=1
`3/2
We carry out the following program: We first solve the equation relating the
gas parameter nλ3 and the fugacity z = eβµ ,
99
100 CHAPTER 6. FERMIONS
we can then find the thermal equation of state via (see 8.41)
2U
p= (6.6)
3V
and we are done.
The dimensionless density parameter
δ = nλ3 (6.7)
characterizes the quantumness of the gas, see Fig. 6.1,
δ 1 → classical gas,
δ 1 → quantum gas. (6.8)
We expand f3/2 (z) for large and small z using the series expansion for
z 1,
z2 z3
f3/2 (z) ≈ z − 3/2 + 3/2 − . . . , (6.9)
2 3
and the integral representation for z 1,
Z ∞
4 x2
f3/2 (z) = √ dx (6.10)
π 0 1 + ex2 −βµ
Z ∞ √ Z ∞
2 y 4 y 3/2 ey−βµ
= √ dy = √ dy
π 0 1 + ey−βµ 3 π 0 (1 + ey−βµ )2
where z = eβµ . The function
ey−ν d 1
y−ν 2
=−
(1 + e ) dy 1 + ey−ν
evaluated at fixed ν = log z = βµ is sharply peaked around ν motivating an
expansion of y 3/2 around this point (Sommerfeld expansion),
3 3
y 3/2 ≈ ν 3/2 + ν 1/2 (y − ν) + 1/2 (y − ν)2 + . . . .
2 8ν
101
Starting from δ = f3/2 (z), (6.3), and using (6.9), we find (fMB (~
p ) is the
1
classical Maxwell-Boltzmann distribution )
z2 1
δ ≈ z− √ z ≈ δ + √ δ 2 > 0,
→ (6.14)
2 2 2 2
2π~2 3/2 h p2 i
hnp~ i ≈ ze−βεp = n exp −
mkB T 2mkB T
3
≈ h fMB (~
p ).
The series expansion of f5/2 (z) provides us with the thermal equation of
state
V z2
pV ≈ kB T z − 5/2 + . . . (6.15)
λ3 2
nλ3
≈ N kB T 1 + √ + . . . , (6.16)
4 2
describing the dilute Fermi Gas as a classical ideal gas with quantum cor-
rections ∝ δ. Note that the chemical potential µ and its dependence on
n appears in the prefactor of the Maxwell-Boltzmann distribution (via the
fugacity z = exp(βµ)).
between the density n and the chemical potential2 µ, and thus to lowest
order
~2
µ(T = 0) = (6π 2 n)2/3 = εF . (6.18)
2m
h π 2 kB T 2 i
µ(T ) ≈ εF 1 − . (6.20)
12 εF
1
p ) = hnp i ≈
fFD (~
eβ(εp −µ) +1
For energies εp~ that are g-times degenerate, the corresponding result is
6π 2
kF3 = n.
g
For example, for an electron gas with s = 1/2, g = 2, and one obtains
kF3 = 3π 2 n. (6.22)
2
Note that εF is one of the characterising parameters of a Fermi gas, while µ describes
the thermodynamic chemical potential.
3
First, recast (6.17) into the form
h π 2 kB T i2/3 ~2
µ 1+ ≈ (6π 2 n)2/3 (6.19)
8 µ 2m
3 h 5π 2 kB T 2 i
≈ N εF 1 + + ... . (6.23)
5 12 εF
With the Sommerfeld expansion to order T 2 , we can obtain the specific heat
π 2 kB T T →0
C V = N kB → 0, (6.24)
2 εF
in agreement with the third law.
Finally, we find the thermal equation of state (or the partition function)
to be
2U 2 h 5π 2 kB T 2 i
p= = nεF 1 + + ... . (6.25)
3V 5 12 εF
Note that the pressure, different from the classical ideal gas, does not vanish
at T → 0 but stays large. It corresponds to the pressure of a system at the
degeneracy temperature
~2 kF2 ~2 6π 2 2/3
kB TF = εF = = n . (6.26)
2m 2m g
Particle-hole asymmetry
The ∝ T 2 correction in the expression for the chemical potential µ(T ), see
Eq. (6.20), is a consequence of the particle-hole asymmetry: The density
of states ρ(ε) of the free Fermi gas in 3D is (we use ρ(k) = k 2 /π 2 for fermions
with spin and ∂k ε = ~2 k/m)
√ r
dk mk m 2mε ε
ρ(ε) = ρ(k) = 2 2
= 2 3
= ρ(εF ) , (6.28)
dε π ~ π ~ εF
104 CHAPTER 6. FERMIONS
with ρ(εF ) = 3n/2εF the density of states for an electron gas of total density
√
n including both spin directions.4 We see that ρ(ε) ∝ ε, i.e., there are
more states available for the electrons at energies ε > εF than for the holes
at ε < εF . Thus, for T > 0 the chemical potential µ has to decrease such
as to conserve the particle number, see Fig. 6.2. Indeed, in the integral
providing the fixed number of particles N ,
Z ∞
N =V dε fFD (ε)ρ(ε), (6.29)
0
the distribution function fFD (ε) is symmetric with respect to µ, but ρ(ε)
is not. For a free massive particle, we find a constant density of states
ρ(ε) (and hence symmetric around εF ) in dimension d = 2, as ρ(k) ∝ k
and ∂k ε ∝ k. In general, with a dispersion ε(k) ∝ k n and ρ(k) ∝ k d−1
in d dimensions, we have a constant ρ(ε) if d = n. Note that electrons
moving in a periodic potential have complicated dispersions and particle-
hole symmetry may appear in specific regions of the Brioullin zone. In
1D mesoscopic physics, the dispersion around the Fermi energy is often
linearized, εk − εF ≈ vF ~(k − kF ), implying particle-hole symmetry within
this approximation.
0 µ εF ε
The most prominent example of a Fermi gas is the electron gas in a metal.
However, other interesting examples are know from astronomy, i.e., the mat-
ter from which stars are made. Depending on the properties of such matter,
gravitational properties (i.e., mass density) versus pressure (i.e., equation of
state), a star may be stable or not.
In white dwarves (made from a 4 He2+ , e− plasma), the pressure is gener-
ated by the Fermi gas of the electrons. At a mass density of ρ ∼ 106 g/cm3 ,
the electron density is of the order of n ∼ 1029 cm−3 (m4 He ≈ 7 · 10−24 g)
4
Sometimes one may encounter the formula ρσ = 3n/4εF for the density of states per
spin σ in an electron gas.
105
much larger than the effective temperature T ∼ 107 K in the core of the
star. We may notice a neat division of tasks between the nuclei (the 4 He-
cores) and the electrons: The nucleons generate the attractive gravitational
force while the electrons, through their degeneracy pressure, generate the
stabilizing pressure. In the non-relativistic case, we obtain the polytrope (or
equation of state)
p = Kρ1+2/3 = Kρ5/3 . (6.31)
Metallic sound
5 2
κ−1 = −V ∂V p = n∂n p = p = nεF
3 3
5
At low temperatures, adjabatic and isothermal compressibilities are the same.
106 CHAPTER 6. FERMIONS
1
2ZnεF 1 mZ p2F Zm 2
c2 = =
= = v , (6.33)
ρκ
nM 3 3 M m 2 3M F
r
Zm
→c = vF , (6.34)
3M
where we have used that ρ = nM/Z with n the electron density, M the mass
of the ions, and Z is the ion’s valency. Using
pF ~(3π 2 n)1/3
vF ∼ ∼ ∼ 108 cm/s,
m m
we obtain a typical sound (or phonon) velocity cmetal ∼ 106 cm/s in a metal.
Note that in these collective excitations both the electrons and the ions are
oscillating. What happens if only the electrons oscillate?
X [~ ~ qi )]2
pi + (e/c)A(~ X
H = − ~
~ i · H,
µ (6.35)
2m
i i
~i ,
~ i = −µB g S
µ
e~
µB = = 0.579 · 10−8 eV/G = 1.165 · 10−19 G cm3 ,
2mc
g = 2[1 + α/2π + . . .] ≈ 2,
A(~ ~ r ) − Φ0 ∇χ(~
~ r ) → A(~ ~ r ),
2π
Φ0
φ(~r ) → φ(~r ) + ∂t χ(~r ),
2πc
Ψ(~r ) → eiχ(~r ) Ψ(~r ), (6.36)
6
√
In 3 He-fluids and in neutron stars c = vF / 3.
6.1. ELECTRONS IN A MAGNETIC FIELD 107
hc
Φ0 = ≈ 4 · 10−7 G cm2 (6.37)
e
is the quantum of magnetic flux.
The statistical mechanics of this system is given by calculating the par-
tition function ZN or Z. A new quantity that we must define in the present
context is the magnetization density 7 (we choose H~ = (0, 0, H) without loss
of generality) (
kB T ∂ ZN ,
m= log (6.38)
V ∂H Z.
It characterizes the response of the system to an applied magnetic field
H. The external field H induces currents or changes the direction of the
magnetic moments µ ~i associated with the spins S
~ i = −g µB S ~i (= ~σi /2) in
the system. Alternatively, we may define m as the change in energy induced
by an applied magnetic field H,
1 D ∂H E
m = −
V ∂H
1 DX E 1 D ∂H E
orb
= µz,i − (6.39)
V V ∂H
i
B = H + 4π mSpin > H,
requires us to compute the spectrum of (6.35) (we ignore the spin compo-
nent). The single particle problem
[~ ~ q )]2
p + (e/c)A(~
H= , HΨλ = ελ Ψλ (6.43)
2m
can be reduced to the displaced harmonic oscillator. We choose the gauge
~ = (−Hy, 0, 0)
A (6.44)
~ = (0, 0, H) according to H
that generates the correct field H ~ =∇
~ ∧ A.
~ The
eigenvalue problem then has the form
~2 h 2πi 2 i
− ∂x − Hy + ∂y2 + ∂z2 Ψλ = ελ Ψλ . (6.45)
2m Φ0
with the cyclotron frequency ωc and the magnetic length `. The latter intro-
duces a new length scale in the problem (besides the interparticle distance
or kF ). The eigenvalues and eigenfunctions are those of a harmonic oscillator
and we find
1 ~2 kz2
ελ = ε(n, kz ) = ~ωc n + + , (6.47)
2 2m
2 2
Ψλ = eikx x eikz z Nn Hn [(y − y0 )/`] e−(y−y0 ) /2` ,
we have kz = (2π/Lz )nz , kx = (2π/Lx )nx , 0 < y0 = `2 (2π/Lx )nx < Ly , and
hence
Lx Ly HA ΦA
nx ≤ = = , (6.48)
2π`2 Φ0 Φ0
with the flux ΦA through the surface determining the degeneracy of the
Landau levels.
ρ ρ
0 H= 0 ε 0 H>0 ε
Fig. 6.4: Sketch p of the density of states for the free electron gas at H = 0
(left) ∝ (3εF /2n) ε/εF and at a finite field H > 0 (right) with the first
peak at ~ωc /2 followed by further equidistant peaks at (n + 1/2)~ωc .
√
discrete spectrum (for the harmonic oscillator) plus a 1D DOS ∝ 1/ ε
DOS for the motion along z, see Fig. 6.4. Going to 2D, one then expects
the quantization into Landau levels to fully manifest itself in a discrete
spectrum, see Fig. 6.5.
ρ(ε ) ΦA Φ states
H= 0
ε
δ−peaks
ρ(ε )
h ωc
H> 0
h ωc / 2 ε
With the degeneracy ΦA /Φ0 per Landau level and the inclusion of all
Landau levels, one obtains a complete set of orthogonal states. Indeed, in
2D, the orbits Ψλ=[n,kx (y0 )] with fixed n and varying kx cover the entire plane
with a resolution given by the Hermite function Hn [(y − y0 )/`], see Fig. 6.6.
The n = 0 orbitals have the extent of the magnetic length, hy 2 i1/2 ∼ `.
The size of the orbital at εF , n = εF /~ωc , is given by the Larmor radius
rL = vF /ωc = `(kF `) (note the increasing powers of y in Hn ).
n= 0 n= 1 n= 2
l
y0 y0 y0
Fig. 6.6: Coverage of the plane with orbits of increasing resolution associ-
ated with higher Landau levels.
At the boundaries, the orbitals are deformed by the box potential and
their energy is raised, see Fig. 6.7. Within a semi-classical picture, the orbits
transform from circular-localized in the bulk to extended skipping orbits at
the boundary. These skipping orbits are the simplest semi-classical picture
of the edge-currents in the quantum Hall effect (QHE).
PΦA /Φ0
By superposition i=1 ani Ψn (y −y0i ) of degenerate states within each
6.1. ELECTRONS IN A MAGNETIC FIELD 111
0 y0 y
Landau level, we can generate a new set of eigenfunctions for the n-th level.
For example, with the gauge A~ = (B/2)(y, −x, 0), one obtains the circular
solutions in Fig. 6.8
2∂ m ∂ n −z z̄/2`2
Ψm,n (z = x + iy) ∼ ez z̄/4` e , (6.49)
∂ z̄ m ∂z n
with ∂z = ∂x − i∂y and ∂z̄ = ∂x + i∂y ,
1
HΨmn = ~ωc n + Ψmn ,
2
n = energy q# and m = angular momentum q#. (6.50)
y y
circular
x x
lattice
Fig. 6.8: Different gauges and superpositions produce other sets of Lan-
dau level functions, e.g., the circular solutions (left) or the lattice solutions
(right).
Using the spectrum (6.47), we can calculate the grand canonical partition
112 CHAPTER 6. FERMIONS
with the three sums over the level degeneracy, the Landau levels, and the
kinetic energy in the z direction. In addition, we have assumed a high
temperature approximation when expanding the logarithm. In the high
temperature limit, the density relates to the fugacity via z = nλ3 , see (6.14).
For small fields, x = ~ωc /2kB T = λ2 /4π`2 1, we can expand and obtain
the following results,
nλ3 V H 1 − x + x2 /2
log Z '
Φ0 λ 1 − (1 − 2x + 4x2 /2 − 8x3 /6)
nλ2 H 1
= V (1 − x + x2 /2) · (1 + x + x2 /3)
Φ0 2x
nλ2 H Φ0 mkB T x2
= V 1 −
Φ0 2π~2 H 6
2
π ~ λ H 2 2 1 λ4
' nV 1 − = N 1 − , (6.51)
12 mkB T Φ20 96π 2 `4
kB T ∂ π n~2 λ2 H
m = log Z = −
V ∂H 6 m Φ20
nµ2B H nµB λ2
= − =− , (6.52)
3kB T 12π `2
∂m nµ2B
χL = =− , (6.53)
∂H 3kB T
1 n 2
χL = − ρ(εF )µ2B = − µ . (6.54)
3 2εF B
0 εF T
|δH| 1 e2
= 4π|χL | = kF 2 ∼ 10−5 (6.55)
H 3π mc
with kF ∼ 1/Å and the classical electron radius e2 /mc2 ∼ 3 · 10−5 Å.
An interesting question is the one about the extent of the regions near the
surface where the currents flow.
We again consider (6.35), but now ignore the gauge field in the kinetic energy
(we write µ~ i = −gµB S~i ≈ −µB~σi with ~σ the Pauli matrices, g ≈ 2, and the
~ of the spin is included in µB ),
X p2
i
H= + µ B σi H , σi = ±1,
2m
i
p2
εp~,σ = + σµB H, (6.56)
2m
with the occupation of these states given by np~,σ ∈ {0, 1}. DefiningP
the total
number of electrons with up (+) and down (−) moments µi , N± = p~ np~,∓1 ,
the total energy of a state {np~,σ } is given by the expression
X p2
E[{np~,σ }] = np~,σ − µB H(N+ − N− ). (6.57)
2m
p
~,σ
114 CHAPTER 6. FERMIONS
X N
X h i
ZN = e−βE[{np~,σ }] = exp βµB H(2N+ − N ) (6.58)
{np~,σ }N N+ =0
X h X p2 i X h X p2 i
× exp −β np~,−1 exp −β np~,x−1 .
2m 2m
{np~,−1 }N+ p
~ {np~,1 }N− p
~
As usual, it is the largest term in the sum that plays a dominant role: Using
f (N+ ) = 2µB HN+ − F0 (N+ ) − F0 (N − N+ ), we have to find the value
N+ = N̄+ maximizing f ,
∂f
= 0. (6.61)
∂N+ N̄+
With ∂N F0 = µ the chemical potential of the ideal spinless Fermi gas, see
(6.14) and (6.18),
kB T log nλ3 < 0,
µ(n) =
π 2
k T
2 (6.62)
εF (n) 1 − B
0,
12 εF
the condition (6.61) tells us that the chemical potentials of the spin ↑ and
spin ↓ systems differ by the magnetic energy 2µB H
Using N̄± = (1/2 ± ε)N and (6.63), we find the relative disbalance
µB HV
ε= . (6.64)
∂n µ N
in the spin population. For the degenerate Fermi gas
∂n 1 3n
= ρ(εF ) = ,
∂µ 2 4εF
6.1. ELECTRONS IN A MAGNETIC FIELD 115
where the additional factor 1/2 appears since we consider the individual spin
configuration separately. The spin disbalance then becomes
2µB H 3 n 3 µB H
N̄+ − N̄− = 2εN = N= N
n 4 εF 2 εF
and the magnetization density takes on the value
χP = −3χL . (6.67)
For the non-degenerate Fermi gas we have the Curie law behavior
∂n 1 ∂ 1 µ/kB T n µB H
= 3
e = , N̄+ − N̄− = N,
∂µ 2 ∂µ λ 2kB T kB T
µ2 H n 2
m = B n, χP = µ . (6.68)
kB T kB T B
These results are summarized in the sketch of Fig. 6.11.
3 n µ2
Fig. 6.11: Temperature de-
2 εF B pendence of the paramagnetic
Pauli spin-susceptibility in a
χP Curie
Fermi gas.
~ 1/ T
0 εF T
116 CHAPTER 6. FERMIONS
8
The finite magnetization density |m(H = 0)| > 0 is a consequence of spontaneous
symmetry breaking where the spins arrange in parallel along some direction in the absence
of an external field.
Chapter 7
Bosons
We study two types of bosonic systems: Photons are the quantum par-
ticles constituting the classical waves of electromagnetic radiation—when
discussing photons in a general context, their particle number is not con-
served and their chemical potential is zero, µ = 0, meaning the cost of
adding a photon is always its energy. Similar is the case for the phonons,
the quantum particle associated with sound waves in a material. Quantum
engineered photonic systems involving cavities and waveguides are an area of
current research, in particular, the non-equilibrium aspect that is associated
with the omnipresent drive and dissipation in this type of systems—steady
states then replace the equilibrium states. The second class of bosons are
massive particles, in particular atoms, that are experimentally studied in
the form of trapped cold gases; below, we will discuss the properties of an
ideal Bose-gas and its Bose-Einstein condensation.
7.1 Photons
in Fock space, where the ~k, λ-mode is occupied n~k,λ ∈ N0 times. The pa-
rameters ~k and λ describe the wavevector and the polarization of the mode.
In a quantization volume L3 , the wavevectors are given by
~k = 2π ~n, ~n ∈ Z3 . (7.2)
L
117
118 CHAPTER 7. BOSONS
For every ~k there exist two polarizations (either linear or circular), specified
by the polarization vectors
Energy, momentum and spin momentum of the photon field are given by
X
E[{n~k,λ }] = ~ωk n~k,λ ,
~k,λ
X
P~ [{n~k,λ }] = ~~k n~k,λ ,
~k,λ
X
~ ~ }] =
S[{n ~k̂ (n~k,+ − n~k,− ). (7.5)
k,λ
~k
1~ ~ = 0 in the Coulomb gauge; the physical fields are transverse, with Coulomb’s
∇·A
~ ~
law ∇ · E = 0 and Gauss’ law ∇ ~ ·B ~ = 0.
2
Alternatively consider a (classical) harmonic oscillator and calculate its (canonical)
partition function by a gaussian integral over phase space (p, q). Then go over to the
quantum mechanical oscillator with energies En = ~ωn and calculate the corresponding
result Z = Tr[exp(−βH)] = 1/(1 − exp[−β~ω)] (without quanta/particle interpretation);
for many independent oscillators with frequencies ωk , we find the canonical partition
function Z = Πk Zk , which is exactly the result found in (7.6). The equality Z = Z (or
F = Ω) follows from (4.20) with z = 1, µ = 0. Thus, we can interpret the quanta of
the harmonic oscillator as bosonic particles with chemical potential µ = 0. The linear
spectrum in n of the harmonic oscillator is an essential part of this argumentation.
7.1. PHOTONS 119
1 ∂ log Z 1
hn~k,λ i = − = β~ω ,
β ∂~ωk e k −1
(7.7)
∂ X
U = − log Z = ~ωk hn~k,λ i (7.8)
∂β
~k,λ
1 ∂ 1 ∂ U
p= log Z = − log Z = . (7.9)
β ∂V 3V ∂β 3V
Here we used that Z depends on V through x = βωk : With ω~k = 2πc|~n |V −1/3 ,
we find ∂V = (∂V x)∂x and ∂V x = −x/3V , ∂x = ωk−1 ∂β , thus ∂V =
−(β/3V )∂β . Note that (7.9) (m = 0) is different from (8.41) (m > 0)
(for massive particles, we have εp ∝ V −2/3 and ∂V = −(2β/3V )∂β ). To
calculate U one substitutes
Z ∞ Z ∞
d3 k
Z
X V →∞ V 2 V
→ 2V = 2 dk k = 2 3 dω ω 2 (7.10)
(2π)3 π 0 π c 0
~k,λ
π 2 (kB T )4 π2 V
= V = kB T, (7.11)
15 (~c)3 15 λ3
4π 2 kB4 T 3 T →∞
cV = → ∞ (7.12)
15 (~c)3
~ω 3 1
u(ω, T ) = 2 3
(7.17)
π c exp ~ω/kB T − 1
dS 1 1 d2 S 1
= =− ln U, 2
=− . (7.20)
dU T αω dU αωU
dS 1 kB d2 S kB
→ = = , 2
= − 2. (7.21)
dU T U dU U
d2 S 1
=− (7.22)
dU 2 αωU + U 2 /kB
7.2 Phonons
Quite often, the vacuum state of a condensed matter system can be viewed
as an elastic deformable medium. The excitations of this medium, the
phonons, then play the same role as the photons in the electromagnetic
vacuum. The important difference is that condensed matter systems are
composed of atoms: The continuum nature of the vacuum in electromag-
netism is replaced by the discrete atomic lattice of condensed matter. In-
stead of infinitely many modes the condensed matter system has 3N modes
corresponding to the 3N degrees of freedom of the N atoms. The spectrum
in momentum space is reduced to the Brillouin zone (BZ) due to the fact that
large wave numbers k > Kn (with Kn the reciprocal lattice vectors) do not
make sense because such modes resolve the (inexistent) structure between
the lattice sites. Furthermore, the spectrum at the BZ edge is deformed as
the phonons “scatter on the lattice” and are thus reflected ∂ω/∂k|BZ = 0.
Thirdly, there are now “3 polarisations”, two transversal and one longitu-
dinal, for every ~k-vektor, transversal shearing modes and longitudinal com-
pression modes. The latter correspond to the hydrodynamic sound in gas
(longitudinal compression modes)—the shearing modes are a characteristic
of the hard condensed matter system (they appear neither in liquids nor in
gases). Note that the classification of longitudinal (k) and transversal (⊥)
(to ~k ) directions only make sense for special directions of ~k in the BZ.
ω ω
+/−
ck
2x
0 k − π/ a 0 k [1,0,0] π/a
For a simple isotropic elastic medium, we can model the phonon disper-
sion by
ωλ = cλ |~k | (7.24)
7.2. PHONONS 123
and the density of states has the form (with d3 k ∼ 4πk 2 dk)
X L 3
ρ(ω)dω ' 4πk 2 dk
2π
λ
X V 1
2
= ω dω
2π 2 c3λ
λ
3V ω 2
= dω, (7.25)
2π 2 c̄3
3 2 1
= + .
c̄3 c3⊥ c3k
We find the size of the isotropic Brillouin zone through the constraint
Z ωmax 3
V ωmax V 3
ρ(ω)dω = 2 3
= 2 kBZ = 3N
0 2π c̄ 2π
3
→ kBZ = 6π 2 n, cf. (6.22)
ωmax = c̄kBZ . (7.26)
kB Θ = ~ωmax . (7.30)
At temperatures T above Θ, all modes are excited and we find the classical
result of Dulong-Petit for the specific heat of a solid
CV = 3N kB , T > Θ, (7.31)
124 CHAPTER 7. BOSONS
12π 4 T 3
CV = N kB , T < Θ. (7.32)
5 Θ
Here is a collection of Debye temperatures
Ag 215 K
Al 428 K
Au 170 K
Cu 343 K
He 25 K
C (diamond) 2250 K
C (graphite) 413 K
Si 645 K
MgO 750 K
SiO2 (quartz) 470 K
NaCl 310 K
Small values of Θ describe “soft”, while large values refer to “hard” materi-
als.
In general, every excitation/degree of freedom contributes to CV ; it is
their spectral density as it follows from the dispersion ω ∼ k n and the
dimensionality of the system (→ ρ ∼ k d−1 ) that determines the temperature
dependence of CV . Simple dimensional considerations, valid for small T such
that the cutoff in the integral is on T and not on ρ, then provide the scaling
behavior
1 ω∼T
Z Z
d ~ω
U ∼ d k β~ω ∼ dk ρ,
e −1 β
Z ω∼T Z
CV ∼ dk ρ ∼ dk k d−1 ∼ k d ∼ ω d/n |ω∼T ∼ T d/n .
7.3. IDEAL BOSE-GAS: BOSE-EINSTEIN-CONDENSATION 125
When all modes are exhausted (with the integral cutoff by ρ) we have U ∝ T
and CV is constant. For phonons in 3D, we have n = 1, d = 3,
CV ∼ T 3 . (7.33)
Equation of state:
kB T kB T kB T
p= log Z = 3 g5/2 (z) − log(1 − z). (7.35)
V λ V
Particle density:
z ∂ 1 1 z
n= log Z = 3 g3/2 (z) + . (7.36)
V ∂z λ V 1−z
Occupation number:
1
hnp~ i = . (7.37)
eβ(εp −µ) −1
Z ∞ ∞
4 2
X z`
g5/2 (z) = − √ dx x2 log(1 − ze−x ) = 5/2
,
π 0 `
`=1
Z ∞ 2 ∞ `
4 x X z
g3/2 (z) = √ dx x2 = . (7.38)
π 0 e /z − 1 `3/2 `=1
( √ √
z + z 2 /2 2 + z 3 /3 3 + . . . , z 1,
g3/2 (z) = (7.40)
ζ(3/2) = 2.612, z = 1.
126 CHAPTER 7. BOSONS
The slope of g3/2 (z) diverges in z = 1. We see that for nλ3 < 2.612, equation
(7.39) can be satisfied with a value z < 1 deriving from nλ3 = g3/2 (z) with
the second term playing no role in the limit V → ∞. In contrast when
nλ3 > 2.612, we find a solution to (7.39) requires that z = 1 − O(1/V ) and
the second term has to compensate the contribution nλ3 − 2.612. In the
limit V → ∞, we find z → 1 and using (7.39), the ground state occupation
hn0 i = hnp~=0 i
z V
hn0 i = = 3 (nλ3 − nBE λ3 ) ∝ V (7.41)
1−z λ | {z }
=2.612
~2 n2/3 4π
kB TBE = · 2/3
(7.43)
2m
| {z } |2.612
{z }
Energie 6.63
2 n 0 /N
λ α1 /T BE−condensation
1
2/3
1/ TBE α 1 / n
0 n 0 TBE T
Fig. 7.5: Left: particles Bose condense when sufficiently many particles
occupy the de Broglie volume λ3 . Right: the fraction of particles in the
condensate grows continuously with lowering temperature T .
One may ask why only the p~ = 0 state condenses, as the momentum
p1 = 2π~/L also goes to 0 for V = L3 → ∞. Thus, we would have
N 1 hn0 i hn1 i
= 3 g3/2 (z) + + (7.46)
V λ V V
P
with the last term the next in the sum p~ hnp~ i. We fix T below TBE and
0 p
determine hn1 i from Eq. (7.37) for a large volume and βε1 1 (we use
hn0 i/V = 1/α > 0),
we can evaluate the equation of state (see 7.35) of the ideal Bose-gas
(
kB T nλ3 T > TBE ,
p = g (z) → k T n 1 − √ + . . . , (7.48)
λ3 5/2
B
4 2 n < nBE ,
(
kB T kB T T→0 T < TBE ,
p = 3
g5/2 (1) = 1.342 3 ∝ T 5/2 → 0, (7.49)
λ λ n > nBE .
This result defines the transition line pBE (T ) in the p–T diagram, see Fig.
7.7,
m 3/2
pBE (T ) = g5/2 (1) (kB T )5/2 .
2π~2
gas
0 T
In a further step, we determine the p–v diagram for the ideal Bose gas.
For a volume v < vBE (T ) (v = 1/n is the volume per particle, n > nBE (T )),
we have p = pBE = const. This is the vapour pressure above the conden-
sate with the gas-phase assuming a fraction nBE (T )/n of the density (the
remaining particles are in the condensate and generate no pressure). At
the transition, nλ3 = λ3 /v = g3/2 (1) and hence TBE (v) ∝ v −2/3 . Using the
expression (7.49) for the pressure at the transition, we find the transition
line in the p–v diagram,
(7.50)
7.3. IDEAL BOSE-GAS: BOSE-EINSTEIN-CONDENSATION 129
The phase diagrams in the Figs. 7.7 and 7.8 resemble those for the real
gas with the liquid–gas transition. We can make an analogy between Bose
gases and real classical gases with the Bose condensed (C) phase correspond-
ing to the liquid and the Bose gas phase (G) with p > 0 corresponding to
the usual gas phase (G). The condensate replacing the liquid, however, has
different, atypical properties. In particular, it does not sustain a pressure,
and thus cannot oppose any compression. One then may observe the fol-
lowing similarities and differences: First, the transition between G and C in
Fig. 7.8 obeys the Clausius-Clapeyron relation with the volume jump from
vG = vBE (T ) to vC = 0, ∆v = vBE (T ). Correspondingly, one obtains the
latent heat (per particle joining the condensate)
where we have used vBE /λ3 = 1/g3/2 (1) for the volume jump. On the other
hand, we cannot reach the p > pBE (T ) region of the p-T diagram. If we
want to cool further below TBE , there is no state that can resist a pressure
and thus the system follows the transition line pBE (T ) until it reaches p = 0,
T = 0 (up to the impossibility to reach T = 0 exactly), see Fig. 7.9. The full
condensation can then be realized by i) ‘pressing’ the gas into the condensate
at p und T (reduce the volume), or ii) reducing the temperature at constant
volume (stretching the de Broglie wavelength beyond the system volume);
naturally, variations of these two methods are possible.
Further, we show (see below) that CV (T ) only shows a kink. The inter-
pretation of the BEC as a phase transition of first or second order according
to the Ehrenfest classification is not really applicable, but note that L&L
classify the transition as being second order.5
For the internal energy U , the entropy S, and the specific heat CV , one
5
There is no finite free energy surface that defines a finite slope ∂p g or ∂T g below the
pBE (T ) transition line.
130 CHAPTER 7. BOSONS
8
sure or by lowering the tem-
perature T → 0 at fixed vol-
condensation ume and following the transi-
via T 0.
tion line pBE (T ).
0 T
finds
U 3 3 kB T3 g5/2 (z),
2 nλ
= pv = (7.52)
N 2 3 kB T g (1),
3 2 nλ 5/2
5 g5/2 (z)
S 2 nλ3 − log z,
= (7.53)
N kB 5 g5/2 (1) ,
2 nλ3
g (z)
15 5/2 9 g3/2 (z)
CV 4 nλ3 − 4 g1/2 (z) ,
= (7.54)
N kB 15 g5/2 (1)
4 nλ3
.
The internal energy appears from the pressure and U = (3/2)pV . The
entropy S follows from
G = N µ = N kB T log z,
∂G
S = − , (7.55)
∂T
∂z 5z g5/2 (z)
= − , (7.56)
∂T 2T g3/2 (z)
where the last equation is derived from ∂ T g5/2 (z) = pλ3 /kB T (using the
and making use of ∂ T g3/2 (z) = nλ3 . The dependence CV (T ), see Fig.
C/Nk B
Fig. 7.10: Specific heat of a
BEC with a cusp at TBE .
3/2
T 3/2
0 TBE T
g5/2 (1) 5
s= kB , (7.58)
g3/2 (1) 2
i.e., the entropy S (and the pressure p) are generated purely by the gas phase.
For T → 0 the gas phase vanishes along with the entropy and the specific
heat, in accordance with the third law of thermodynamics. The latent heat
` = T ∆s = T s originating from (7.58) agrees with (7.51). Note that the
latent heat and volume jump are realized when hitting the transition point
pBE (T ) and either condensing the gas by further decrease of p and T along
the transition line or by squeezing the gas into the condensate at fixed p and
T . Thereafter, there is no new thermodynamic phase which can be further
cooled with a new thermodynamic Gibbs function, i.e., we cannot calculate
an entropy/volume jump as a jump in derivatives of two free energy surfaces
across the transition line as is the case for a usual first-order transition.
132 CHAPTER 7. BOSONS
Chapter 8
Magnetic order
1
I.e., the critical exponents do not depend on the dimensionality d and n of space and
of the order parameter.
133
134 CHAPTER 8. MAGNETIC ORDER
J X~ ~ X
~i · ~h = HH + Hext ,
H = − Si · Sj − µS (8.1)
2
hi,ji i
J > 0, ferromagnetic coupling, ↑↑: E = −J; ↑↓: E = J,
~
|Si | = 1, classical spins,
µ = magnetic moment.
The energies Es,t = hΨs,t |HH2 |Ψs,t i of the singlet and triplet states are
separated by
i.e., the effective Hamiltonian for the spin problem is of the Heisenberg type.
The exchange interaction J = ∆E/2 (1/2 from the definition |S ~i | = 1) is
purely electronic (non magnetic). Typical values for J can go up to ∼ 103 K,
hence, the temperature scale of electronic magnetism can be rather high.
An important parameter in the study of magnetic systems and phase
transitions in general is the range of the interaction. Restricting the latter
to nearest neighbors, it is the number of nearest neighbours z that plays an
important role. For a cubic lattice, z = 2d where d denotes the dimension-
ality of the crystal.
The interaction J can be negative—in this situation one talks about
anti-ferromagnetic coupling. If the lattice is bipartite, an anti-ferromagnetic
order can appear with opposite feromagnetic orders on the two sub-lattices.
If the lattice is not bipartite the system will be frustrated, for example on a
triangular lattice in d = 2, see Fig. 8.1. In a magnetically doped solid (for
?
J<0
example Mn in Cu, Ag, Au) with different distances between the magnetic
atoms, the coupling Jij can be +/− randomly4 and the result is a spin glas
with maximal frustration.
The Heisenberg model has O(3)-Symmetry. Generally one talks about
O(n)-symmetric spin models with the continuity of the symmetry playing
~i in the Heisenberg model is
a particularly relevant role—the spin-field S
4
E.g, generated via the RKKY interaction (Ruderman-Kittel-Kasuya-Yosida) and ran-
dom distances between dopants.
136 CHAPTER 8. MAGNETIC ORDER
8.1.2 XY model
If one limits the degree of freedom of the spin to one axis (an “easy axis”
anisotropy in solids) one obtains the Ising model,
JX z z X
H=− Si Sj − µh Siz = HI + Hext . (8.4)
2
hi,ji i
Its symmetry is given by the discrete group Z2 = {1l, R with R2 = 1l} (Ising
symmetry). The Ising model of magnetism is equivalent to the lattice gas
model where a lattice site can be occupied or not and to the problem of
binary alloys, for example Zn-Cu in a bcc lattice. Furthermore, the Ising
model and the van der Waal gas-liquid transition are in the same universality
class.
~i denotes the
in magnetic systems, see also Eqs. (6.38) and (6.39). Here, S
(classical) spin on site i and µ the magnetic moment. In the homogeneous
case hm(~
~ r )i = hmi
~ and thus
~
M
hmi
~ = (8.7)
V
~ being the (extensive) magnetic mo-
defines the magnetisation density, M
ment. As for the density n = 1/V , the magnetisation density m
~ takes on
the role of a thermodynamic variable that interdepends on quantities such
as entropy or energy.
Another quantity that characterizes the magnetic system and that can be
conveniently measured is the magnetic stiffness or (magnetic) susceptibility
∂m
χ = . (8.8)
∂h
This quantity should be viewed as the magnetic analog of the compressibility
κ = V −1 ∂p V in the mechanical context with variables p and V or the specific
heat CV = T ∂T S in the caloric context with thermal variables T und S. In
the sections 6.1.1 and 6.1.2, we have already seen the Landau- and Pauli-
susceptibilites of the electron gas.
The spin system described by the Hamiltonian (8.1) has a fixed number
N of spins or degrees of freedom and thus the canonical ensemble is the
appropriate one for its statistical description. We then have to adapt the
canonical partition sum (3.38)
Z
ZN (V, T ) = d3N p d3N q ρ(p, q)
Z 3N 3N
d pd q H(p, q)
= exp − ,
h3N N ! kB T
to the spin lattice problem. The continuous variables (p1 , . . . , q3N ) are
replaced by the set (S~1 , . . . , S
~N ) and the integral over phase space goes
5
Be careful about units; here, m(~
~ r ) assumes the unit of moment per volume since it
is a field, while m ~i is a moment. We go from one to the other by multiplying with
~ i = µS
the unit cell volume Vuc .
138 CHAPTER 8. MAGNETIC ORDER
over to the sum over all spin configurations {S~i }N with an integration
1
over angles for n > 1 on each site. The situation is particularly sim-
ple for the Ising model, where the partition sum is over all realizations
of {S1 = ±1, S2 = ±1, . . . , SN = ±1}.
Good thermodynamic variables for the spin system are the temperature
T (and its conjugate entropy S) and the (intensive) external magnetic field
~h (with the magnetization M ~ assuming the role of the conjugate extensive
variable). The variables p and V , pressure and volume, are usually ignored
when discussing generic magnetic phenomena. However, in real compressible
material they may have to be included into the description as well, e.g., in
analyzing the phenomenon of magnetostriction. The partition function of a
magnetic system then is usually written in the form
X ~i }]
H[{S
ZN (T, ~h) = exp − (8.9)
kB T
~i }
{S
with angular integrations per site included in the sum when n > 1. The
factor
R h−3N can be dropped as we have replaced the phase-space integral
dpdq by a dimensionless sum. Similarly, the Gibbs factor 1/N ! becomes
irrelevant since the moments/spins are fixed on the lattice sites and there is
no permutation of variables.
The associated thermodynamic potential is a function of the intensive
variables T and ~h, as well as the number N playing the role of the extensive
variable. Hence, the good potential is the Gibbs fee energy,
∂mi 1 ∂2G
χij = =− . (8.12)
∂hj V ∂hi ∂hj
F = G + ~h · M
~ = F (T, M
~ , N ).
8.4. BROKEN SYMMETRY AND ORDER 139
~h = ∂ ~ F, (8.13)
M
∂2F
χ−1
ij = V . (8.14)
∂Mi ∂Mj
hm(T,
~ ~h, N )i = 1 ∂ ln ZN (T, ~h ),
βV ∂~h
hm(T
~ )i = lim lim hm(T,
~ ~h, N )i. (8.16)
~h→0 N →∞
We then obtain the phase diagram as shown in figure 8.3 (compare to the
p–T phase diagram in Fig. 1.17). The transition at Tc in an h = 0 field
is of second order: The magnetisation grows steadily from 0 when T drops
below Tc . The transition at T = const < Tc and h = 0 is first-order: The
magnetisation jumps from −m ê to m ê with m = |m|. ~ A sketch of the
equation of state hm(T,
~ h)i is shown in Fig. 8.4. For ~h = 0 the spontaneous
magnetisation hm(T
~ )i plays the role of an order parameter that specifies
how, i.e., in which direction, the symmetry is spontaneously broken.
As in the crystal and the superfluid, the spontaneous symmetry breaking
goes hand in hand with the appearance of long-range ordering (LRO) in
the system.6 The latter is quantified by the spin-spin or magnetisation-
magnetisation correlator
C(~r, ~r 0 ) = hm(~ ~ r 0 )i = µ2 hS(~
~ r ) · m(~ ~ r 0 )i .
~ r ) · S(~ (8.17)
6
A special case may appear in 2D, see later.
8.4. BROKEN SYMMETRY AND ORDER 141
|~ r 0 |→∞
r−~
C(~r, ~r 0 ) → ~ )i2 6= 0
hm(T (LRO). (8.18)
δ 2 ln Z m(~r 0 )e−βH
R
δ
= (8.21)
∂βh(~r )δβh(~r 0 )
R
δβh(~r ) e−βH
m(~r 0 )m(~r )e−βH m(~r 0 )e−βH m(~r )e−βH
R R R
= R −
e−βH
R R
e−βH e−βH
= hm(~r )m(~r 0 )i−hm(~r )ihm(~r 0 )i ≡ hhm(~r )m(~r 0 )ii.
from which the magnetic information derives from (again, be careful with
units, specially volumes)
δF
hi (~r ) = , (8.22)
δhmi (~r )i
δ2F δhi (~r )
χ−1 r, ~r 0 ) =
ij (~ 0
= . (8.23)
δhmi (~r )iδhmj (~r )i δhmj (~r 0 )i
The last line is found from (where cr denotes the chain rule and fd stands
142 CHAPTER 8. MAGNETIC ORDER
We begin with the simplest case of the Ising model and then proceed with
the mean-field theory of the O(n) model involving the XY and Heisenberg
spins.
JX z z X
H=− Si Sj − µh Siz , (8.27)
2
hi,ji i
where each pairing hi, ji appears twice. Every spin is exposed to the external
field h and the effective or “Weiss” field of its neighbours,
X
µh0i = µh + JSjz ≈ µh + zJhSi. (8.28)
j
7
The functional derivative of F [m(~r )] with respect to m(~r 0 ) is defined as
~
e.g., with F [m(~r )] = dd r f [m(~r ), ∇m(~
R
r )], we find that
δF
Z
∂f δm(~r )
h ∂f ~
δ ∇m(~ r )i
= dd r +
δm(~r 0 ) ∂m(~r ) δm(~r )0 ~
∂ ∇m(~ r ) δm(~r 0 )
Z h ∂f ∂f i
= dd r δ(~r − ~r 0 ) − ∇
~ δ(~r − ~r 0 )
∂m(~r ) ~
∂ ∇m(~ r)
∂f ~ ∂f
= −∇ .
∂m(~r 0 ) ~
∂ ∇m(~ r 0)
8.5. MEAN-FIELD THEORY 143
J X
HMF = zN hSi2 − µh0 Siz . (8.29)
2
i
The following alternative derivation makes the type of approximation, the
dropping of quadratic fluctuation terms, more transparent. We rewrite the
Hamiltonian (8.27)
JX z X
H = − [(Si − hSi) + hSi][(Sjz − hSi) + hSi] − µh Siz
2 | {z }
hi,ji i
| {z }
δSiz δSjz
J X X X
= − N zhSi2 − JzhSi δSiz − µh Siz − J δSiz δSjz
2
i i hi,ji
J X
≈ zN hSi2 − µh0 Siz ,
2
i
where we counted each pairing twice and ignored the quadratic terms in the
fluctuations; this derivation again leads us to the mean-field Hamiltonian
(8.29).
Yet another approach is via Bragg-Williams theory,Pwhere one writes
an Ansatz for a single-particle Hamiltonian HMF = −α i Siz + const and
minimises the Bogoliubov functional (5.67) with respect to α, see later.
The Hamiltonian (8.29) is (up to the constant) a sum of single parti-
cle Hamiltonians H1 = −αS z , α = µh0 and the corresponding statistical
mechanics can be easily calculated, ZN = Z1N with Z1 the single spin parti-
tion function. The average spin hSi is obtained through the self consistency
equation derived from statistical averaging,
∂
hSi = ln Z1 , Z1 = eβα + e−βα = 2 cosh βα,
∂βα
∂ µh + zJhSi
= ln 2 cosh βα = tanh .
∂βα kB T
Inverting the above equation provides us with a transzendental equation for
hSi,
µh = kB T artanh hSi − zJhSi (8.30)
which can be solved graphically as sketched in Fig. 8.5. The critical tem-
perature Tc separates regimes with positive and negative slopes ∂h/∂hSi at
the origin,
∂µh ∂
' (kB T hSi − zJhSi) = 0
∂hSi ∂hSi
→ kB Tc = zJ. (8.31)
144 CHAPTER 8. MAGNETIC ORDER
the dimensionality d of space appears only in the value of the critical tem-
perature kB Tc = 2dJ, indicating that higher dimensions favour the ordered
phase more strongly.
The self-consistency equation (8.32) provides us with further details of
the magnetic transition. At Tc , Eq. (8.32) reduces to µh = kB T hSi3 /3,
telling us that the order parameter hmi = (µ/Vuc )(3µh/kB T )1/3 on the crit-
ical isotherm scales as
see Fig. 8.4. The critical exponent δ for the order parameter as a function
of h on the critical isotherm T = Tc assumes the value δ = 3 in mean-field
theory. Furthermore, we can calculate the magnetic suceptibility
χ = ∂h hmi. (8.38)
Above Tc , hmi = 0 and we find (from (8.32)) the derivative ∂hSi h = (kB Tc /µ)
(T /Tc − 1), from which we obtain χ = (µ2 /Vuc kB Tc )/τ and hence the scaling
of the susceptibility at h = 0 near Tc is
The same exponent γ but with a different prefactor is obtained when evalu-
ating χ below Tc where hmi =6 0. The critical exponent γ for the magnetic
susceptibility χ assumes the value γ = 1 in mean-field theory.
The last exponent characterizing the phase transition is the one for the
specific heat cV . Evaluating this quantity from the free energy f0 (see
below) one finds a jump, telling us that the critical exponent vanishes within
mean-field theory,
cV ∝ τ α with α = 0. (8.40)
The critical exponent α for the specific heat cV assumes the value α = 0
in mean-field theory.
A different formulation of the mean-field theory for the Ising model has
been given by Bragg and Williams. Using the Bogoliubov variational princi-
ple, we construct the free energy F = U −T S by calculating S and estimating
146 CHAPTER 8. MAGNETIC ORDER
This expression vanishes for the ordered phases with hSi = ±1 and amounts
to S = N kB ln 2 for an unpolarized (disordered) state with hSi = 0. For the
internal energy one uses the mean-field approximation
1
Uh = − N zJhSi2 − N µhhSi. (8.42)
2
T >Tc
Fig. 8.7: Free energy
f0 f0 (T, hSi) at h = 0. As
T drops below Tc the
minimal free energy is
Tc T < Tc attained at a finite value
of hSi.
0 S
∂f0 ∂fh
h= or = 0, (8.44)
∂µhSi ∂µhSi
8.5. MEAN-FIELD THEORY 147
Starting from the Hamiltonian (8.1) of the Heisenberg model, one again
makes a mean-field ansatz (see 8.28)
Using ~h = hêx , hS
~ i = hS x iêx , one finds
XY model: S x = cos ϕ
0
dϕ eβµhx cos ϕ cos ϕ
R
x
hS i = , (8.48)
dϕ eβµh0x cos ϕ
R
Evaluating the integrals over angles ϑ and ϕ, one arrives at the self-consistency
equations for the XY model,
I1 [βµhx + βzJhS x i]
XY model: hS x i = , (8.50)
I0 [βµhx + βzJhS x i]
0 x
The expansion for small hS x i provides us with the expression for Tc and the
T -dependence of the order parameter around Tc . Using x ≡ βzJhSi we find,
kB T x2
Ising: =1− → kB Tc = zJ,
zJ 3 √ √
hSi = 3 1 − t. (8.52)
kB T 1 x2
XY: = − → kB Tc = zJ/2,
zJ 2 16 √ √
hSi = 2 1 − t. (8.53)
kB T 1 x2
Heisenberg: = − → kB Tc = zJ/3,
zJ 3 45 p √
hSi = 5/3 1 − t. (8.54)
these three models, we have results for the critical temperature that depend
on d and n as
2d
kB Tc = J , (8.56)
n
i.e., larger values of d (more neighbors) increases Tc , while a larger n pro-
viding more degrees of freedom to the order parameter decreases Tc .
The mean field theory provides us with a simple compact and universal
picture of magnetic phase transitions. In fact, mean-field theory always gives
the same structural result: The continuous (second-order) phase transitions
are described by a self-consistency equation which around Tc is given by
T 4u 3
Φ=Φ+ Φ ,
Tc r0
Φ = order parameter, from where we obtain a free energy (by integration)
|r0 | T 2
f (T, Φ) = − 1− Φ + uΦ4 .
2 Tc
Unfortunately, mean-field theory misses an important point: Its super (i.e.,
independent on d and n) universal character predicts transitions that do not
exist and gives us exponents which are independent of the dimensionalities
d and n of the model. This deficiency will be cured when fluctuations are in-
cluded into the picture. In spite of these drawbacks, its simple and compact
form motivates the use of mean-field theory as a first starting point when
dealing with a phase transition. There are even situations where mean-field
theory is “correct”. This is, e.g., the case when the mean field is supported
by many (or even an infinite number of) neighbours. In this situation, the
averaged field acting on the order parameter is well defined, i.e., the devi-
ations (= fluctuations) from the mean are small. Many neighbours appear
when going to higher dimensions and long-range interactions, and mean-field
theory turns out correct for the magnetic transitions of the Ising-, XY-, and
Heisenberg models in d > 4. Furthermore, mean-field theory is very accurate
for BCS (Bardeen-Cooper-Schrieffer) superconductors with a large correla-
tion or coherence length ξBCS , nξBCS3 1 with n the electron density. In
the following, we study how the mean-field theory fails in magnetic systems
in low dimensions, where z is small and fluctuations in the mean field are
expected to be large.
A simple argument demonstrates that the (short range8 ) Ising model cannot
order at finite temperatures in one dimension. Consider the ordered Ising
8
See D.J. Thouless, Phys. Rev. 187, 732 (1969) for a discussion of long-range-order in
one-dimensional Ising systems with long-range (slowly decaying) interactions.
150 CHAPTER 8. MAGNETIC ORDER
chain in Fig. 8.9. The introduction of a domain wall costs an energy 2J.
For N → ∞, we have Fdw < 0 for all T > 0 and thus there is no order present
at any finite temperature. The only place where an ordered ferromagnetic
state can survive is at T = 0 and we can call this a T = 0 critical point,
Tc = 0.
Domain walls are topological objects or defects. The only way to generate
them is in pairs or by dragging them into the system across its boundaries,
see Fig. 8.10.
The 1D-Ising model can be solved exactly with the help of transfer
matrices. Consider N Spins Si = ±1 with periodic boundary conditions
S1 = SN +1 . The partition sum ZN (T, h) can be written in the form
N
XX X X µh
ZN (T, h) = ... exp β JSk Sk+1 + (Sk + Sk+1 ) . (8.58)
2
S1 S2 SN k=1
|Sk ihSk | = 1)
X X
(with = hS1 |T N |S1 i
k S1
= λN+ + λN
−. (8.60)
∂m µ2
χ= = e2J/kB T (8.65)
∂h Vuc kB T
diverges for T → 0, signalling a T = 0 phase transition. Indeed, going to
high temperatures, hm(T, h)i ≈ µ2 h/Vuc kB T and we find the Curie-Weiss
law χ ∝ (T − Tc )−1 with Tc = 0, again indicating a T = 0 phase transition.
9
Here, the bra and ket vectors denote classical spin states.
152 CHAPTER 8. MAGNETIC ORDER
We first give an argument (due to Peierls) that the Ising model in 2D shows
a phase transition at finite temperatures. We will avoid discussing the exact
solution here (this exact solution is due to Onsager and produces a critical
temperature given by kB Tc ≈ 2.269185 J). Next, we discuss the continuous-
spin XY model and its special quasi–long-range ordered low-temperature
phase, as well as the special type of transition between the short-range or-
dered phase at high temperatures and the quasi–long-range ordered state
(due to Berezinskii and Kosterlitz-Thouless). Finally, we make a few com-
ments on the non-ordering (in 2D) Heisenberg model.
We show that at h = 0 and positive but sufficiently small T > 0 there exists
a finite magnetisation hSi > 0 in the 2D Ising model.
Peierls’ argument:
Consider the T = 0 ground state of the 2D Ising model with all spins pointing
up, Siz = ↑, denoted by + in Fig. 8.11. As the temperature is raised, spin–
↓ domains (denoted by − in Fig. 8.11) appear. The associated domain
walls separating spin–↑ from spin–↓ regions cost an energy 2JL, where L
is the length of the domain’s circumference. The growth/extension of such
domains costs an energy ∝ L, the same scaling as the entropy associated
with the geometry and location of such a domain. As a result, the energy
dominates over the entropy at sufficiently low but finite temperatures T > 0
and the phase stays ordered.10
10
Note that a domain in 1D costs L0 , in 2D L1 , in d dimensions Ld−1 , hence additional
dimensions “stabilise” the ordered phase.
8.7. MAGNETIC ORDER IN TWO DIMENSIONS 153
The factor g(L) > 1 favors the appearance of an additional loop that flips
the spin, while the Boltzmann factor accounts for the loop’s cost in energy
154 CHAPTER 8. MAGNETIC ORDER
and reduces its probability to show up—obviously, this is the usual interplay
of entropy and energy associated with a loop of size L.
In our estimate of g(L), we start at an arbitrary point in the domain.
For the surface of the domain, we set the maximum value (L/4)2 . In the
first step, we choose one from four directions, in every following step from
three (we ignore crossings). Thus we obtain
L 2 1
g(L) ∼ · 4 · 3L−1 · . (8.66)
4 2L
The last factor takes into account that we can start at any point of the loop
and can go around it in two directions, clock- and anti-clockwise. We then
find
∞
X L L −2JL/kB T
p− < 3 e
24
L=4,6,...
∞
1 X −d −αL 2J
< e , α= − ln 3
24 dα kB T
L=4,5
1 d e−4α 5 e−4α
= − < . (8.67)
24 dα 1 − e−α 24 1 − e−α
For sufficiently large α or small temperatures T , the average spin hS ∗ i =
1 − 2p− is bounded from below,
5 e−4α
hS ∗ i > 1 − (8.68)
12 1 − e−α
and choosing kB T < J, e−α < 0.41, we find that hS ∗ i > 0.98. Thus below
kB T = J the phase is definitely ordered; in fact, the ordered phase survives
until kB Tc ≈ 2.269185 J (Onsager).
This type of proof for the existence of an ordered phase can be used
for other models. The central object entering the proof is the loop of
length L flipping a domain. One then may argue, that a finite-temperature
phase transition appears if there is a temperature 0 < T < ∞ where
loops start proliferating. The free energy of such a loop is Floop (L) ∼
2LJ − T ln(L3L /24) ∼ L(2J − T ln 3) and we find that loops flood the plane
when T > 2J/ ln 3.
8.7.2 XY model in 2D
~i = (cos ϑi , sin ϑi ),
The XY Hamiltonian can be written in polar coordinates S
X X
H = −J ~i · S
S ~j = −J cos(ϑi − ϑj ) (8.69)
hi,ji hi,ji
8.7. MAGNETIC ORDER IN TWO DIMENSIONS 155
Ja 2 Z d
d r ~ 2
≈ E0 + (∇ϑ) , (8.70)
2 ad
where E0 is the ground state energy and ϑ(~r ) is the (continuous) deforma-
tion field away from ϑ0 . Note that we kept the underlying space dimension
d arbitrary; we will focus on the planar case d = 2 later. The parameter J
describes the stiffness of the ordered phase; a is the lattice spacing.
We want to check whether the ordered phase as obtained from the mean-
field solution, is stable under fluctuations of the field ϑ(~r ) or whether it gets
destroyed as in the 1D Ising model. To this end, we calculate the expectation
value (we set ϑ0 = 0)
D[ϑ(~r )]e−βHel eiϑ(~r )
R
hSx i = hcos ϑ(~r )i = R (8.71)
D[ϑ(~r )]e−βHel
with the Gaussian Hamiltonian
Ja2
Z d
d r ~ 2
Hel = [∇ϑ]
2 ad
Ja2 X 2
= k ϑ~k ϑ−~k , (8.72)
2
~k
Ja2 k2 h|ϑ~k |2 i/2 = kB T /2. S1d is the surface of the unit sphere in d dimensions.
156 CHAPTER 8. MAGNETIC ORDER
~ r ) · S(0)i
C(~r, 0) = hS(~ ~ = < exp[i(ϑ(~r ) − ϑ(0))]
= exp(−h[ϑ(~r ) − ϑ(0)]2 i/2). (8.74)
Using
2
a d Z 2kB T (1 − cos ~k · ~r )
h[ϑ(~r ) − ϑ(0)] i = dd k
2π Ja2 k 2
d
S1 2 kB T
, d > 2,
(2π)d d−2 J
∼ kB T R (8.75)
πJ ln a , d = 2,
kB T x ,
d = 1,
J a
we find that in
r→∞
3D: C(~r, 0) → const. > 0, LRO, (8.76)
k T R
~ 0) R→∞ B
2D: C(R, → exp − ln
2πJ a
a kB T /2πJ R→∞, alg.
= → 0, QLRO,
R
k T x x→∞, exp.
x→∞ B
1D: C(x, 0) → exp − → 0, SRO.
2J a
In words: thermally excited, elastic, low-energy fluctuations tend to destroy
the T = 0 order at T > 0. In 3D, C(~r ) stays finite assymptotically and long-
range order survives for 0 < T < Tc . In 1D, C(x) goes to 0 exponentially
for every temperature T > 0 and the order is destroyed. The 2D case is
marginal: C(R) decays to 0, but only slowly, algebraically. An algebraically
decaying correlator is typical for a system at a phase transition Tc , where
no length scale is apparent in the system12 . In 2D, at all temperatures T ,
the system behaves like a critical one at Tc . As the long-range order is
only slowly destroyed one talks about quasi–long-range ordered = QLRO
phase. This result is suspicious—it is hardly possible that QLRO survives
to arbitrary high temperatures. Indeed, we will see below, a very special
phase transition (due to Bereszinskii-Kosterlitz-Thouless) will intervene and
replace the algebraic decay by an exponential one, such that the order is
short-ranged at high temperatures.13
12
Algebraic functions have “no scale”.
13
The mistake is in our use of the Gaussian approximation, which fails at high temper-
atures where another type of fluctuations, vortices, comes into play.
8.7. MAGNETIC ORDER IN TWO DIMENSIONS 157
Hohenberg-Mermin-Wagner Theorem
Thus, we obtain the Einstein relation between the diffusion coefficient D of the Brownian
motion hδx2 i(t) ∼ Dt and the dissipation η or mobility µ = 1/η. As an excercise one may
calculate the T → 0 limit of this “high temperature” result which will replace h[x(t) −
x(0)]2 i ∝ T t by h[x(t) − x(0)]2 i ∝ ~ ln ωUV t, with ωUV a UV cutoff.
158 CHAPTER 8. MAGNETIC ORDER
Crystal, phonon = Goldstone mode: Using the Hamiltonian and the ex-
pression for the displacement field u of the phonons (note the appearance of ωk = ck,
c2 = ε/µ, the ratio of elasticity and mass density)
X
H = ~ωk (a~†k a~k + 1/2), ωk = ck,
~
k
~ 1/2
u~k = (a~k + a†−~k ) = displacement,
2ωk µL
In the limit T → 0, all modes k > 0 remain unoccupied and the mode k = 0 leaves the
order unchanged; an ordered ferromagnetic state survives.
mc2 n0
Z
n0 ∼ n − dd k ha~†k a~k i, hn~k i =
[1 + 2n(Ek )],
Ek 2n
mc2 n0 T
Z Z
T%
∼ n − dd k → ··· − T dk kd−3 ,
nEk2 0
where we have used the Bogoliubov theory of interacting Bose gases to obtain the expres-
sion for hn~k i (with c2 = nV /m the speed of sound and n, n0 the total and condensate
densities). With Ek2 ∝ k2 at small k, we find that the normal density n − n0 diverges
in d = 1, 2 and there is no room for a condensate. Comparing the three cases, we see
that for sound waves (in crystals and in the BEC) the occupation number n(ωk ) and the
8.7. MAGNETIC ORDER IN TWO DIMENSIONS 159
Fig. 8.13: Spin configurations with ϑ(R)~ = kθ + ϕ defining the angle of the
XY spin with respect to the x-axis at the position R~ = R(sin θ, sin θ). From
left to right k = 1, ϕ = 0; k = −1, ϕ = 0; k = 1, ϕ = π/2; k = −1, ϕ = π/2.
~ R
Using ∇ϑ( ~ ) = ±R
~ ⊥ /R2 , one finds that
I
d~` · ∇ϑ
~ = ±2π; (8.81)
amplitudes (hu2 i ∝ ωk and hn~k i ∝ 1/Ek ) both contribute a factor k, producing the overall
factor k2 in the denominator, while for the ferromagnet, the dispersion is quadratic and
all k2 contribution is from n(εk ). Note that the ideal Bose gas behaves as the ferromagnet
with a quadratic dispersion and spectral factor unity. See HẆatanabe and H. Murayama
for more on this.
160 CHAPTER 8. MAGNETIC ORDER
E VV
Fig. 8.15: Attractive
potential between a vor-
tex and an anti-vortex
ln R /a
at a distance R.
0 1 R/a
−1
1 1
−1
With Eq. (8.87), we have found a Hamiltonian that can describe the
phase transition that leads to the loss of the quasi-LRO that is present in
the low-temperature phase of the XY -model. The need for the existence of
such a phase transition in 2D became clear from the work of Berezinskii, who
found that the low-temperature phase is characterized by quasi-LRO and
stiffness, and high-temperature expansions, telling that the order is short-
ranged, C(R) ∼ e−R/ξ . It was Kosterlitz and Thouless,16 who in 1973 (in J.
Phys. C 6, 1181(1973)) proposed the following scenario for the transition: In
the low temperature phase the vortices are pairwise bound and thus do not
destroy the quasi-LRO. More pairs appear with increasing temperature—
they screen the interaction between larger pairs which finally leads to the
dissociation of pairs at a temperature TBKT ≈ πJ/2kB . The appearance
of free vortices above TBKT destroys the quasi-LRO. The calculation is a
beautiful example of renormalisation in classical statistical mechanics, see
Sec. 17. They obtained the Nobel prize for this work in 2016 (together with
Haldane, in the more general setting of ‘topology in condensed matter’).
Berezinskii-Kosterlitz-Thouless transition
and thus it becomes favourable for the system to generate such excitations.
Above TBKT , vortices are bound at short distances R < ξ(T ) and therefore
the quasi-order survives on the scale of ξ. At distances R > ξ, free vortices
appear and C(R) goes to 0 exponentially. Below TBKT , we have ξ(T <
TBKT ) = ∞. Above the transition at TBKT , the coherence length diverges
upon approaching TBKT from above according to
b
ξ(T ) ∼ exp p .
(T − TBKT )/TBKT
Ja2
Z d Z d
d r d r ~
with H = d
(∂i~n ) · (∂i~n ) − µh · ~n,
2 a ad
~n = (~π , σ),
where we chose ~h = hêz and the symmetry broken along ~n0 = (~0, 1) = êz .
The analysis of the nlσ model is non-trivial and we cannot go into detail
here. The renormalisation group analysis delivers a ferromagnetic phase
transition at
d−2
kB Tc ≈ 2πJ , (8.92)
n−2
and hence Tc = 0 for the n = 3, d = 2 Heisenberg model in 2D (note
that in addition to spin-waves, also topological excitations in the form of
Skyrmions are relevant in the characterization of the system, e.g., the cor-
relation length, see Polyakov’s book). In 3D the Heisenberg model shows
a transition at Tc ∼ J. The overall behavior of the critical temperatures
in magnetic transitions as a function of n and d is sketched in Fig. 8.18,
where we include our knowledge gained so far from mean-field theory and
the analysis of fluctuations.
164 CHAPTER 8. MAGNETIC ORDER
MF
MF
MF
QLRO
Tc n BKT
3 Heis
2 XY
LRO
1 Ising
0 1 2 3 4 d
Gas–Liquid Transition
165
166 CHAPTER 9. GAS–LIQUID TRANSITION
with the parameters a and b describing the attractive and repulsive interac-
tions. The p-v-diagramm then exhibits the characterisitc isotherms shown
in Fig. 9.1.
Tc
vc v
vA vB v
We then need the expressions for the internal energy u and the entropy s:
the internal energy u of a real gase depends on the volume v, see (1.42),
∂u a
∂v u|T = T ∂T p|v − p ⇒ = , (9.7)
∂v T v2
a consequence of the interatomic attraction (i.e., the stored potentiel en-
ergy). We obtain the energy difference ugas − uli by straightforward integra-
tion
Z B
Tli =Tgas a a a
ugas − uli = dv 2 = − + . (9.8)
A v vgas vli
1 1 ∂u ∂u p
ds = (du + pdv) = dT + dv + dv (9.9)
T T ∂T V ∂v T T
168 CHAPTER 9. GAS–LIQUID TRANSITION
cv R
ds = dT + dv
T v−b
Z B
Tli =Tgas R vgas − b
⇒ sgas − sli = dv = R ln . (9.10)
A v−b vli − b
On the other hand, we can determine the work or area in the p-v diagram
Z B Z B Z B
dv a dv
dv p = RT −
A A v − b li v2
B a B
= RT ln(v − b) + (9.13)
A v A
and find that it agrees with the right hand side of (9.12), hence
Z B
dv p = p(vgas − vli ). (9.14)
A
We find that the equilibrium condition µli = µgas can be reduced to the
condition of equal shaded areas in Fig. 9.2, which is commonly known as
the Maxwell construction for the physical correction of the van der Waals
Isotherms. One easily shows, that this construction corresponds to the con-
vex completion of the free energy surface in section 1.5.5.
and B 0 ) and regions with metastable phases (α0 and β 0 ). Repeating this
construction for all pressures p < pc , we obtain the p-v diagram in Fig.
9.4. At the fog line first liquid droplets appear in the gas phase, while
at the boiling line first gas bubbles show up in the liquid. Going on a
path γ around the critical point we can go continuously from the gas- into
the liquid phase—this is the core statement made in the title of van der
Waals’ thesis. Note that one cannot liquefy the gas through pressurizing
when the temperature is high, T > Tc —gas- and liquid phases cannot be
distinguished at these high temperatures. Interesting further topics are the
vc v
boiling line
d3N q
Z
1 V N h X i
ZN (V, T ) = exp −β V (~
q i − ~
q j ) . (9.15)
N ! λ3 VN
i<j
The r−6 interaction between the neutral atoms is long-ranged, i.e., every
atom interacts with many neighbors, motivating a mean-field description
with a potential
1 V0 −r/R0
V (r) = − e (9.16)
8π R03
170 CHAPTER 9. GAS–LIQUID TRANSITION
d3N q Y
Z
sr
ZN = χσ (~qi − ~qj ) (9.19)
VN
i<j
with χσ (r) the characteristic function on r > σ. This is the partition func-
tion of a gas of hard spheres of radius σ which we approximate by
sr
ZN = (V − N b)N /V N , (9.20)
where b is the excluded volume per particle. Collecting factors, we find the
partition function of the van der Waals gas
1 V N βV0 N 2 /2V V − N b N
ZN (T, V ) = e (9.21)
N ! λ3 V
as well as its free energy
aN 2
F (T, V, N ) = Fid−gas − − kB T N log(1 − bN/V ) (9.22)
V
with a = V0 /2. The equation of state follows from the partial derivative
p = −∂V F and we obtain the van der Waals equation
kB T a kB T b kB T a
p= − 2+ = − , (9.23)
v v v v−b v − b v2
with v = V /N the volume per particle. The free energy per particle assumes
the form (see Eq. (8.23))
a
f = kB T log[λ3 /(v − b)] − 1 − ,
(9.24)
v
which coincides with the ideal gas result when a = 0 and b = 0. The entropy
per particle is given by −∂T f ,
Finally, the Gibbs free energy per particle assumes the form
2a kB T v
g = f + p v = kB T log[λ3 /(v − b)] − 1 −
+ . (9.27)
v v−b
It is interesting to plot
R p this quantity at fixed T by using the fact that
g(p, T ) = g(T, p0 ) + p0 v(p)dp (since v = ∂p g|T ) and use the form of v(p)
sketched in Fig. 9.5(a): we follow this curve from D to E to F to G, see
Fig. 9.5(b), and use that ∂p v is the curvature, then g first rises as a concave
curve, becomes convex at E, and concave again at F. √ Near E (F), we have
p−pE ∝ −(v −vE )2 (p−pF ∝ (v −vF )2 ) and since δv ∝ δp, we find that the
curvature diverges and g has cusps at E and F. The convex branch between
E and F is unstable.
E
p G
(a) g (b) C’’
G
F C C’
vli v gas v
E
C’ C
C’’
D D
F p
Fig. 9.5: p-v diagram (a) and Gibbs energy per particle g(p) at fixed T
(b) for the van der Waals gas with two stable (solid lines), two metastable
(dashed) and one unstable branch (dotted).
ν ∝ π 1/3 (9.31)
for the ‘order parameter’ ν(h) versus its conjugate field π at Tc . This
corresponds to the mean-field exponent δ = 3, ν ∝ π 1/δ . The corresponding
result in the magnetic setup is identical, m ∝ h1/3 . The correction ν 4 tells
that cubic curve is asymmetric, steeper at negative values of ν and flatter
on the positive side. This tells that the gas volume grows faster than the
liquid volume shrinks when the temperature T decreases below Tc , see Fig.
9.6.
which relates πgl to the volumes νgas and νli , see Fig. 9.6,
2
πgl = 4τ − 3τ (νgas + νli ) + 3τ (νgas + νli νgas + νli2 ) + . . . . (9.33)
We then have to solve the cubic problems π(νgas ) = πgl , π(νli ) = πgl , with
π(ν) given by (9.29), for νgas (πgl ) and νli (πgl ) and solve Eq. (9.33) for πgl (τ );
9.1. THE VAN DER WAALS GAS 173
the third solution νunst to the cubic equation defines the unstable volume
with maximal Gibbs energy. Fortunately, to lowest p order we can solve the
problem more easily. Assuming a solution ν = α |τ |, we can rewrite Eq.
(9.29) in the form
3
π − 4τ = − α3 − 6α |τ |3/2 + O(τ 2 ), (9.34)
2
with the solutions (using α(α2 − 4) = 0)
p p
πgl ≈ 4τ < 0, νgas ≈ 2 |τ |, νunst ≈ 0, νli ≈ −2 |τ | (9.35)
The coexistence region together with the diameter is shown in the v–T
diagram of Fig. 9.7; it corresponds to the ferromagnetic region in the m–T
diagram of magnetism, see Fig. 8.6.
In Fig. 9.8, we sketch the analog of the m–T –h diagram in Fig. 8.4 for
the gas–liquid transition; note that Tc in Fig. 9.8(b) is corresponds to the
tip of the ruled surface in Fig. 9.8(a) in rotated coordinates.
When transforming v to p, we obtain the p–T diagram with the first-
order transition line pgl (T ). We can obtain the shape of the transition
line by inserting the results (9.36) into (9.33) and obtain
48 2
πgl ≈ 4τ + τ . (9.38)
5
The insets in Fig. 9.7 show the first-order line in the gas–liquid transition
together with the analogous line in the magnetic transition.
Continuing the transition line πgl within the π–τ diagram beyond τ = 0,
we move on the critical isochore where ν = 0, see (9.29)
π h
gas
ν m
liquid 0 τ 0 τ
gas
0 τ 0 τ
liquid
Fig. 9.7: Gas–liquid coexistence region in the v–T diagram. When com-
pared to the ferromagnetic magnetization in the corresponding m–T diagram
of magnetism, we note the absence of symmetry m ↔ −m in the gas–liquid
coexistence region. Insets: Note the similar asymmetry in the p–T dia-
gram versus the h–T diagram with its symmetric first-order transition line
at h = 0.
m
p
T
C h
0
l
Tc
g T
v
Fig. 9.8: Surfaces of thermodynamic equilibrium states for the gas–liquid
transition in p–v–T space as compared to the h–m–T diagram for the ferro-
magnetic system.
The critical isochore in the vdW problem plays the same role as the h = 0,
T > Tc line in the magnetic problem. It is along this line where the critical
scaling behavior develops when approaching the second-order critical point
‘from above’. When the critical point is approached along another line in
the π–τ (or h-τ ) diagram, this scaling behavior is modified.
Next, we calculate the isothermal compressibility κT = −v −1 ∂p v—
its analogue in the magnetic system is the suceptibility χ = ∂h m, see Eq.
(8.38). Above Tc , τ > 0, the derivative of Eq. 9.29 on the critical isochore
ν = 0 provides us with the result
1 ∂ν 1
κT = − ≈ . (9.40)
pc ∂π 6pc τ
Below, we determine the analogous quantity along the coexistence bound-
9.1. THE VAN DER WAALS GAS 175
aries πgas (νgas ) and πli (νli ). Using again Eq. (9.29), we have ∂ν π = −6τ +
18τ ν − (9/2)ν 2 and inserting ν 2 ≈ −4τ > 0, we obtain the compressibility
below Tc ,
1 ∂ν 1
κT = − ≈ , (9.41)
pc (1 + ν) ∂π 12pc |τ |
half the value at temperatures T > Tc , τ > 0. The results Eqs. (9.40) and
(9.41) scale ∝ 1/|τ | as predicted by mean-field theory. The analogous result
for the magnetic system is χ ∝ 1/|τ |γ with γ = 1.
As a last quantity, we calculate the specific heat cv (T ) at constant
volume. On the critical isochore ν = 0, we use the expression Eq. (9.26) for
the internal energy and obtain the standard ideal gas result
3
cv = ∂T u|v = kB (9.42)
2
valid above Tc . Below Tc , we have to take the volume term into account,
since the system is in a mixed gas–liquid state on the critical isochore. The
volume ratios at a given ν are
3k
2 B
0 τ
Chapter 10
Lee-Yang Theorem∗
V
Fig. 10.1: Shape of the inter-
action potential describing (in-
finitely) hard spheres.
a
r
−V
0
For a hard sphere system, only a finite number of particles M (V ) can fit
in a volume V , ZN >M (V ) (V ) = 0, and the grandcanonical partition function
must be a polynomial of order M in the fugacity z,
Z(T, V, z) = 1 + zZ1 (T, V ) + z 2 Z2 (T, V ) + . . . + z M ZM (T, V ). (10.1)
The equations of state are obtained from
kB T
p = log Z(T, V, z),
V
1 z
= ∂z log Z(T, V, z). (10.2)
v V
As all the coefficients Zk in (10.1) are positive, the function Z(z ∈ C) has no
zeroes on R+ ⊂ C. Therefore, the functions p(z), v(z) and p(v) are analytic
177
178 CHAPTER 10. LEE-YANG THEOREM∗
on R+ and p, −∂p/∂v ≥ 0,
∂p ∂p ∂z kB T 1
= =
∂v ∂z ∂v zv ∂v/∂z
kB T V
= − 3 ≤ 0,
v (hN 2 i − hN i2 )
where hN 2 i − hN i2 = z∂z (z∂z log Z) = z∂z (V /v).
As there are no singularities, no phase transition appears in a finite system.
Only when going to the thermodynamic limit V → ∞, n = N/V = 1/v =
const, there can be singularities in the quantities
1
p = kB T lim log Z(T, V, z),
V →∞ V
1 z
= lim ∂z log Z(T, V, z). (10.3)
v V →∞ V
Lee and Yang have shown that phase transitions are controlled by the dis-
tribution of zeroes in Z(z ∈ C)—phase transitions appear when these zeroes
approach the real axis as sketched in Fig. 10.2, thereby cutting the axis into
separate segments.
exists for all z > 0 and is independent of the volume V if the surface
∂V ≤ V 2/3 .
Theorem 2: Let SR be a segment of R+ , U (SR ) ⊂ C a region around SR
with no zeroes of Z. Then (1/V ) log Z converges uniformly in U for
V → ∞ and its limit F∞ (z) is analytic.
1. Let SR = R+ , see Fig. 10.3. Then there exists only a free phase.
z z v
2. Let R+ = SR1 + SR2 . Then there exist two phases that transform into
each other via a phase transition.
The behavior of the function F∞ (z) ∝ p across the two segments SR1
and SR2 determines the type of transition: The transition is first-order if
the derivative ∂z F∞ exhibits a discontinuity or jump. A continuous first
derivative and a discontinuity or jump in the second derivative ∂z2 F∞ signals
a second-order transition, see Fig. 10.4
The Lee-Yang circle theorem tells that for the hard sphere lattice gas
with a potential as sketched in Fig. 10.1 the zeroes of Z(T, V, z) lie on a unit
circle in the complex plane. At any finite V , these zeroes appear in discrete
complex-conjugate pairs away from the real axis. In the thermodynamic
limit, the zeros accumulate and approach the real axis, ‘pinching it off’ at
z = 1, what results in two phases with a transition at z = 1.
180 CHAPTER 10. LEE-YANG THEOREM∗
Chapter 11
Landau Theory
In this chapter, we turn our focus to phenomena that appear in the proximity
of a phase transition. We start with a phenomenological mean-field theory
for the order parameter Φ (including also a gradient term), the Landau
theory of phase transitions, and then discuss the various types of transitions,
second-order, first-order, and multicritical points. We start with a short
overview.
11.1 Overview
181
182 CHAPTER 11. LANDAU THEORY
The original Landau expansion with powers of Φ has later been extended
to a field-theoretical framework (Landau-Ginzburg-Wilson) including also
gradient terms ∇Φ. The Landau free-energy functional F [Φ(~r )] (or effective
coarse grained Hamiltonian) over the order-parameter field Φ(~r ) can be
understood as the result from taking a partial trace (coarse-graining) over
microscopic degrees of freedom Φµ described by a microscopic Hamilonian
H[Φµ (~r )],
Z Y
Z = D[Φ(~r )]TrΦµ δ[Φ(~r ) − Σ(Φµ (~r ))]e−βH[Φµ (~r )]
~
r
Z
= D[Φ(~r )]e−βF [T,Φ(~r )] , (11.2)
Q
with the product ~r δ[Φ(~r ) − Σ(Φµ (~r ))] filtering out all Φµ configurations
which give the same Φ. Here, Σ denotes a function that takes the micro-
scopic degrees of freedom Φµ (~r ) near ~r as an input to generate a coarse-
grained variable Φ(~r ), e.g., see (14.26). Hence, while the Hamiltonian H
and field Φµ describe the microscopic problem, the free energy F and the
order-parameter field Φ describe an effective coarse-grained problem. The
underlying idea is that F describes the correct long wave length and low-
energy physics that is relevant for the analysis of the phase transition. As an
illustration of such a coarse graining procedure, we may refer to the discrete
microscopic spin degrees of freedom S~i in an Ising model that become con-
tinuous variables φ(~r ) = hS(~r )i under averaging over small local domains;
the resulting effective φ4 theory, see Eq. (8.43), then describes the Ising spin
problem as an effective continuous theory. Note that the coefficients in the
coarse grained free energy or effective Hamiltonian in general depend on T ,
which is due to the coarse-graining over microscopic degrees of freedom in
Eq. (11.2).
Note the differences and similarities between our mean-field approach in
chapters 8 and 9 and the present discussion: in Chaps. 8 and 9, we have
started from a microscopic Hamiltonian and have introduced a Weiss-field
by replacing the neighboring microscopic degrees of freedom with its local
average. This mean-field theory has 1/z as its small parameter, with 1/z
small implying small fluctuations. The original Landau free energy was just
an expansion of a phenomenological free energy f (T, Φ) near Tc in a small
order parameter Φ, so the small parameter here is the order parameter Φ
itself. Taking the saddle-point approximation to the path integral (11.2),
ZMF = e−βFMF (T ) .
then provides us with a mean-field theory of the type constructed in Chaps.
8 and 9, bringing us full circle back.
However, using the path-integral in Eq. (11.2), we can go further and
include fluctuations around the saddle-point. The simplest approximation is
11.2. LANDAU FREE ENERGY FUNCTIONAL 183
Within Landau theory, we consider the part F (T, Φi (~r )) of the free energy
that is associated with the phase transition, Ftotal = F + FRest . The field
Φi (~r ), i = 1, . . . , n describes an n components local order parameter. Ac-
cording to Landau, one constructs the free energy F under the following
provisions:
Note that Φ ~ does not describe the microscopic spin field S(~ ~ r ), but
the corresponding order parameter. In particular, in the Ising model
|S z | = 1, but 0 < Φ < 1; according to mean-field theory, we have
Φ = hS z i, Φ(~r ) = hS z (~r )i` , where h. . .i` is a local average over a few
cells.
3. A field ~h conjugate to Φ
~ breaks the symmetry of F and generates
~ is found by
a finite order parameter even above Tc . The value of Φ
minimising the Legendre transformation
Z
G = F − dd r ~h(~r ) · Φ(~
~ r)
~ r ) = − δG .
Φ(~ (11.4)
δ~h(~r )
Inversely, the conjugate field ~h follows from
~h(~r ) = δF δG
→ = 0, (11.5)
~ r)
δ Φ(~ ~ r)
δ Φ(~
~ as a constraining parameter in G.
where we have used Φ
4. The parameters r, u, u0 , w, c, . . . in f , Eq. (11.1), are found by compar-
ing experimental and theoretical results (fitting parameters). When
n > m experimentally determined quantities have fixed the m param-
eters of the theory the latter becomes predictive. In some cases the
parameters can also be calculated from a microscopic theory. This
is the case, e.g., for the derivation of Ginzburg-Landau theory from
BCS-theory (Gorkov).
For the Ising, O(n), gas–liquid,1 . . . problems, we use a scalar order param-
eter Φ and expand f in Φ2 ,
r
f (T, Φ) = Φ2 + uΦ4 . (11.6)
2
1
There is no broken symmetry in the gas–liquid transition; nevertheless, we can let the
volume v per particle assume the role of an order parameter
11.3. SECOND-ORDER PHASE TRANSITIONS 185
h = rΦ + 4uΦ3 .
This provides us with two results, the critical isotherm, see Fig. 11.9,
h 1/3
Φ= ∼ h1/δ , T = Tc , (11.9)
4u
186 CHAPTER 11. LANDAU THEORY
above and below the transition and taking the second derivative in T we
obtain the change in the specific heat associated with the transition,
(
∂2f 0, T > Tc ,
cV = −T = (11.12)
∂T 2 2
T r0 /8u, T < Tc .
We thus find that in the Landau mean-field approximation the specific heat
jumps. Denoting the exponent of the specific heat by α, see Fig. 11.9,
1
cV ∼ (11.13)
|T − Tc |α
this jump is consistent with a value α = 0. Note that the specific heat jump
described by (11.12) describes only the singular contribution associated with
the transition; other degrees of freedom contribute to cV , entropy s, and f
with a regular part.
Examples of second-order transitions are the magnetic transitions of the
Ising, XY, and Heisenberg models. As discussed in detail in Chap. 9, the
mean-field theory of the van der Waals gas exhibits a second-order transition
at the critical point pc , vc , Tc as well, although no symmetry is broken. The
role of the order parameter then is assumed by the volume (or density)
difference Φ ↔ vgas − vli .
Another classic example for a second-order phase transition is the normal-
metal–superconductor (or fermionic superfluid) transition. The microscopic
origin of the superconducting transition is found in a pairing of (fermionic)
electrons into (bosonic) Cooper pairs which then can condense in k-space.
2
The derivative ∂h of (11.9) gives 1 = (r + 12uΦ2 )∂h Φ and thus χ = 1/r at Φ = 0,
χ = 1/(r − 12ur/4u) at Φ2 = −r/4u.
11.3. SECOND-ORDER PHASE TRANSITIONS 187
linear/metallic
0 Tc T 0 1 2 3 4K T
Fig. 11.2: Left, sharp specific heat jump typical for a mean-field type transi-
tion as it appears in the fermionic normal-metal–superconductor transition.
The mean-field character is a consequence of the overlapping Cooper pairs
(→ many ‘neighbors’ defining a sharp ‘Weiss’ field). This is different for
a bosonic superfluid transition, right, where the coherence length is of the
order of the particle distance.
Different from the Weiss mean-field theory in Chap. 8, the Landau free
energy functional involves a gradient term (∇Φ)2 and we can define further
quantities (and their critical exponents) associated with the spacial behavior
of Φ. One such quantity is the non-local susceptibility χ(~q ) defined in
188 CHAPTER 11. LANDAU THEORY
For a (n) vector-valued order parameter Φi , one can define two suceptibilities
χk and χ⊥ , that are directed longitudinally and transverse to the direction
of order, respectively. The relevant susceptibility tensor has been defined
in (8.23) and proceeding as in (11.15) with (11.3), we obtain (see also Eq.
(11.16))
χ−1 q ) = T G−1
ij (~ q ) = (r + 4uΦ2 + cq 2 )δij + 8uΦi Φj .
ij (~ (11.21)
5
The superfluid transition in 4 He is called a λ-transition because of its broad smearing
(in the shape of a λ).
190 CHAPTER 11. LANDAU THEORY
A ξ B
ξ
ξ
K
K
Tc T
Fig. 11.4: Domains of alternative phases near the transition are of typical
size ξ. E.g., a typical critical (K) bubble/droplet in the disordered (B) /
ordered (A) phases above and below Tc measures ξ(T ). The critical domains
K represent the scale invariant phase at the transition. On approaching Tc ,
the correlation length ξ diverges from either side.
where6
1 1
χk (~q ) = , χ⊥ (~q ) = . (11.23)
r + 12uΦ2 + cq 2 r + 4uΦ2 + cq 2
While χk is identical to the scalar susceptibility χ in (11.16), we find a
different expression for the transversal susceptibility
(
1/(r + cq 2 ), T > Tc ,
χ⊥ (~q ) = (11.24)
1/cq 2 , T < Tc .
Thus, the transverse part of the order parameter correlator becomes mass-
less below Tc (above Tc , the parameter r > 0 plays the role of a ‘mass’;
alternatively, it generates a finite ‘length scale’ for the correlations),
T
G⊥ (~q ) = hδΦ⊥ (~q )δΦ⊥ (−~q )i = ,
cq 2
1
G⊥ (~r ) ∼ . (11.25)
rd−2
6
We write χ−1 2 2 2 2
ij = αδij + βφi φj /φ with α = r + 4uΦ + cq and β = 8uΦ . Using
k
δij δij = n, δij φi φj /φ2 = 1, (φi φj /φ2 )(φi φj /φ2 ) = 1 and χ−1
k = Pij χ−1 −1
ij , (n − 1)χ⊥ =
⊥ −1
Pij χij , 11.23 follows immediately.
11.4. FIRST-ORDER PHASE TRANSITIONS 191
mexican hat
f Fig. 11.5: Mexican hat illus-
trating longitudinal modes that
change the order parameter am-
plitude and cost quadratic (in
0 δΦ) energy and transverse modes
φ that cost ‘no’ energy in the limit
k → 0.
φ φ
The simplest generic form of the Landau free energy f (T, Φ) for a first order
phase transition is given by
r
f (T, Φ) = Φ2 − wΦ3 + uΦ4 , (11.27)
2
where
r = r0 (T − T ∗ ), w, u > 0, constant.
(r − 3wΦ + 4uΦ2 )Φ = 0,
(r/2 − wΦ + uΦ2 )Φ2 = 0,
T =Tc
f f
T >Tc
T = T **
T < Tc
0 st Φ 0 Φ
1
T =T *
that
9w2
r∗ = 0, r∗∗ = . (11.28)
16u
Note the different separations r∗∗ − rc = w2 /16u and rc − r∗ = 8w2 /16u.
The caloric properties of the transition can be calculated from the free
energy (
0, T > Tc ,
f= 2
(11.29)
(r − rc )Φmin /2, T < Tc ,
from which we obtain the entropy jump ∆S = −∂f /∂T ,
r0 w2
SΦmin − S0 = −
8u2
and the latent heat ` = −Tc ∆S,
r0 w2 Tc
`= . (11.30)
8u2
The jump in S vanishes with the cubic term w → 0. Note that on top of
the entropy associated with the transition, there is a positive regular part
11.4. FIRST-ORDER PHASE TRANSITIONS 193
T* Tc T ** T
ξ(T ) appearing in χ(~q ) is cut off at Tc . The first-order phase transition thus
shows no scaling laws and we cannot define any critical exponents.
Note that our analysis shows a jump in the order parameter from 0 to
Φmin > 0. The expansion of f as a power series in Φ requires that Φ > Φmin
remains small; in this situation, one talks about a weak first-order transition.
Otherwise, higher-order corrections to f become important.
Typical examples for the appearance of a first-order transition are the
nematic to isotropic transition in liquid-crystals (Maier-Saupe model), see
Fig. 11.8 and the liquid to crystal transition. The order parameter of the
I-N transition is not a vector but a director written in a matrix form.7 The
7
For a uniaxial nematic, Qij = S(ni nj −δij /3) with ~n a unit vector, the Frank director.
194 CHAPTER 11. LANDAU THEORY
molecules
Fig. 11.8: The transition be-
tween isotropic (left) and ne-
matic (right) phases of molecules
is first-order.
isotropic
nematic
order parameters in the liquid–crystal transition are the density waves ρG~
defining the density modulation
X ~
ρ = ρ0 + ρG~ eiG~r .
~
G
~ i add up to zero,
The cubic term in f is allowed if the three wavevectors G
~1 + G
G ~2 + G~ 3 = 0, and is given by
X
ρG~ 1 ρG~ 2 ρG~ 3 .
~ 1 ,G
G ~ 2 ,G
~3
A first order transition can also appear in a model without cubic terms, for
example f = (r/2)Φ2 + uΦ4 + u0 Φ6 with u < 0. This model can also describe
a second-order phase transition, namely for r < 0, u > 0. In a system with
many order-destroying variables, for example T and a mixing concentration
x, one can manipulate r and u such that both transitions are present in
the phase diagram and meet at the tricritical point. The appearance of
such multicritical points is the focus of the next section. Before, we briefly
summarize and compare the main results of our discussion of second- and
first-order transitions.
χ
Φ
cV jump
scaling
h
0
scaling Tc Tc T Tc T
T
Φ cV χ
h δ− function
0
jump Tc Tc T T* Tc T ** T
T
hysteresis
Fig. 11.9: Phase diagrams (left), specific heat (middle), and susceptibility
(right) of second-(top) and first-order (bottom) mean-field transitions. The
scaling behavior and divergencies in the second-order transition are replaced
by jumps and a knee in the first-order transition. The correlation length ξ(T )
exhibits the same behavior as the susceptibility χ(T ).
1
f (T, ∆, Φ) = r(T, ∆)Φ2 + u(T, ∆)Φ4 + u0 Φ6 . (11.31)
2
The parameter u0 > 0 guarantees the stability of f . For u > 0, the line
r(T, ∆) = 0 gives a line of second-order phase transitions. For u < 0, there
are first-order transitions when f = ∂Φ f = 0 at Φ 6= 0,
The resulting u-r phase-diagram with the tricritical point at the origin is
shown in Fig. 11.10.
The transition lines of both the first- and second-order transitions in
the u-r plane are easily transformed to the physical variables T and ∆, see
Fig. 11.11. Again, the two transition lines join smoothly with their slopes
196 CHAPTER 11. LANDAU THEORY
u f
f
nd Φ
Φ 2
tricritical point
r
f f
st
1
Φ
Φ
properly aligning,8
∂∆ ∂r/∂T rT
= − ≡− ,
∂T r=0 ∂r/∂∆ r∆
∂∆ 0
rT − uuT /u TC: r=0=u rT ∂∆
= − 0
= − = . (11.33)
∂T r=u2 /2u0 r∆ − uu∆ /u r∆ ∂T r=0
Note that the conjugate (to T and ∆) variables S and x jump on the first-
∆
order line. As usual, these jumps ∆S and ∆x are related through Clapeyron
equations. For example one finds the jump in s to be
∂f rT |u|
∆s = − ≈− → 0 at TC. (11.34)
∂T 4 u0
A new feature is the coexistence of three phases with Φ = 0, Φ = ±Φc on
the first-order phase boundary.
8
Along r(T, ∆) = 0, we have 0 = dr = rT dT + r∆ d∆.
11.4. FIRST-ORDER PHASE TRANSITIONS 197
The critical behaviour at the tricritical point is different from the normal
second-order critical behaviour. We approach the tricritical point TC along
the line u = 0 with r = r0 TTC τ , τ = 1 − T /TTC . Then, the order parameter
Φ increases as
−r 1/4
Φ=± ∼ τ β with β = 1/4. (11.35)
6u0
With h the field conjugate to Φ, we find the equation of state from h = ∂Φ f ,
Finally, to find the critical exponent α of the specific heat, we require the
free energy f ,
(
0, T > TTC ,
f= p (11.38)
(r/2) −r/6u , T < TTC → f ∼ |τ |3/2 ,
0
α β γ δ ν η (11.40)
C : 0 1/2 1 3 1/2 0
TC : 1/2 1/4 1 5 1/2 0
By comparing the terms uΦ4 and u0 Φ6 to (r/2)Φ2 , we can identify the regions
where critical and tricritical fluctuations dominate. Within these regions,
the free energies can be approximated by
r
C : f ≈ Φ2 + uΦ4 , u0 Φ6 small,
2
r
TC : f ≈ Φ2 + u0 Φ6 , uΦ4 small. (11.41)
2
198 CHAPTER 11. LANDAU THEORY
3 He acts as an order parameter for phase separation. The ± signs for the
4 He sites allow to define the (magnetic) average hSi i as an order parameter
that describes superfluidity. The spin Hamiltonian driving these magnetic
and phase-separation order parameters looks pretty complicated, see Blume,
Emery, Griffiths, PRA 4 1071 (1971). One finds the T −∆−h-diagram shown
in Fig. 11.14 (Griffiths, PRL 24, 715 (1970)) with a three-phase coexistence
(two magnetic phases that can get split by h and a paramagnetic phase)
on the first-order critical line D and three second-order lines meeting in a
point; the merging of these three second-order lines gives the name to the
tricritical point.
T
λ−line
2 nd
A TC
B’
B
∆
D
h
Fig. 11.14: T -h-∆ phase diagram for the 3 He–4 He mixture. The region (or
surface) A below the λ-line (a line of critical points corresponding to the
λ-line in Fig. 11.13) describes the superfluid phase in the He mixture and a
first-order or coexistence surface with two degenerate magnetic phases (with
OP ±hSi selected by ±|h| → ±0) in the spin lattice model. The surface A
is connected via D (corresponding to the split lines terminating the phase
separated region in Fig. 11.13) to two other first-order (coexistence) surfaces
B and B 0 extending into the finite h regions. In the magnetic analogon, the
two phases coexisting in the surface B involve different values of hS 2 i (and
similar for B 0 ). The three critical lines terminating A (the λ-line), B, and
B 0 meet in the tricritical point TC. Along the line D, three phases coexist.
We close with a brief overview of multicritical points, see Fig. 11.15. The
lesson to learn here is that different phase transition lines connect with
different constraints. E.g., a second-order line can branch off a first-order
line at a finite angle, while it extends smoothly across a tricritical point.
200 CHAPTER 11. LANDAU THEORY
h h λ−line h
tetracritical
critical CP point
point
critical
endpoint
T T T
λ tricritical λ
T T
point TC
x
konj.
to ∆ phase
separation
∆ x
bicritical
h m
point
m phase
separation BP
konj.
to h
T T
−c −c
0 0
Fig. 11.16: Phase diagram with a Lifshitz point involving a spatially mod-
~ is vector-valued, the transition
ulated phase (left). If the order parameter Φ
to the modulated phase is second-order (right).
9
with ferromagnetic in-plane coupling, ferromagnetic out-of-plane coupling, and anti-
ferromagnetic next-nearest-neighbor coupling.
202 CHAPTER 11. LANDAU THEORY
Chapter 12
203
204 CHAPTER 12. FIELD THEORY: GAUSS AND GINZBURG
LMFT Φmin
Φ
− dim Φ −space
8
theory is not exactly solvable. The uΦ4 term introduces interactions between
the fluctuations (mode-mode-coupling). We will make a first attempt to
treat effects of the Φ4 term in Sec. 12.3 on the level of the Hartree- or self-
consistent field approximation; the full treatment via the renormalization
group will be discussed later in Sec. 16.2. Next, we discuss how to deal with
the functional integral over the order parameter field in Eq. (11.2), both on
a lattice in real space and in ~k-space.
2
Note the similarity with the problem of an elastic manifold in a potential V (u) as
described by the configurational integral
Z
dd r {C[∇u(~
~ r )]2 /2+V [u(~
R
Z= D[u(~r )] e−β r )]}
(12.3)
and touched upon in Sec. 5.1.5 (Eq. (12.3) describes the classical limit of (5.15) after a
trivial integration over π). Hence, the order-parameter field Φ(~r ) can be understood as
the deformation field u(~r ) of an elastic manifold trapped in a square/quartic potential
V (u).
3
We will generalize this later by first expanding (12.2) near a homogeneous minimum
and only then setting u to zero.
12.1. LATTICE FIELD THEORY AND CONTINUUM LIMIT 205
Φmin
x f
Φ Φ
f
x
Φ Φ min Φ
Fig. 12.2: Minimal configuration (solid line) and fluctuations (dashed line)
for the potential shown to the right. Below Tc , the symmetry has been
broken by selecting one minimum.
identify the minimal configuration Φmin and the energetically less favourable
fluctuations δΦ. The continuum limit of (12.4) is obtained from the limit
v → 0 with V = N v = const and the substitutions
ν → ~r,
Z d
X d r
f (Φν ) → f [Φ(~r )],
ν
v
Z d
1X 2 1 d r ~ 2
Cνν (Φν − Φν ) →
0 0 c [∇Φ] (12.5)
2 0 2 v
ν,ν
1X
with c = Cν,0 Rν2 for an isotropic model.
d ν
4
We do not worry about the normalization of the functional integral over Φν , what
accounts for an undetermined constant in the entropy, which is irrelevant in the discussion
of critical phenomena.
206 CHAPTER 12. FIELD THEORY: GAUSS AND GINZBURG
Dropping the elastic and quartic terms in the Hamiltonian of Eq. (12.7), we
easily find the partition function
Z N Z ∞
r N
−βH[Φ(~
r )]
Y 2πkB T −(r/2kB T )Φ2ν
ZN = D[Φ(~r )]e = dΦν e = .
ν −∞ r
(12.8)
Its associated free energy scales linearly with N . In order to include the
elastic term as well, we analyze the model in Fourier space.6 The inclusion
5
The expression (12.7) strictly follows from the discrete version (12.4); later in the text,
we will sometimes use v = 1, implying that the coupling constants scale differently with
length. Such different notation is more convenient, e.g., in the discussion of scaling. Also
note that N in the Fourier transforms below goes over to V when we choose v = 1.
6
Alternatively, the discrete quadratic problem (12.4) without quartic term can be solved
directly using the following statements about Gauss integrals: We start from the simplest
expression r
Z ∞
2 2π λ2 /2C
dx e−Cx /2+λx = e (12.9)
−∞ C
that we reduce to a simple Gauss integral through quadratic completion. Its n × n matrix
generalisation reads
n
Z Y r
h 1 i (2π)n h1
−1
i
dxi exp − xi Cij xj + λi xi = exp λi Cij λj , (12.10)
i=1
2 det C 2
where the matrix Cij = hi|C|ji is taken to be real and symmetric. In this case, C can be
diagonalised, hp|C|p0 i = δpp0 Cp , with Cp real and the transformation matrix Mpi = hp|ii
can be chosen to be real and symmetric (with hi|pi the orthonormal eigenfunctions of C).
Using
X X
xi = hi|pixp , λj = hj|piλp ,
p p
| det M | = | det M 2
| 1/2
= | det 1|1/2 = 1,
we find
Z Yn n Z n
s
−xi Cij xj /2+λi xi
Y −Cp x2
p /2+λp xp
Y 2π λ2p /2Cp
dxi e = dxp e = e ,
i=1 p=1 p=1
Cp
12.1.1 ~k-space
dd r −i~k·~r
Z Z d
1 d r i(~k−~k 0 )·~r
Φ~k = e Φ(~r ), e = δ~k,~k 0 ,
v N v
1 X i~k·~r 1 X i~k·(~r−~r 0 )
Φ(~r ) = e Φ~k , e = δ(~r − ~r 0 ).
N V
~k ~k
dd k
Z
1 X
→ . (12.11)
V (2π)d
~k
where ~0k refers to the constraint of taking only half of all ~k values. Note
P
that with this counting, N degrees of freedom in real space go over into N
degrees of freedom in ~k-space, namely N/2 ~k-points and 2 real fields <Φ~k
and =Φ~k per ~k-point.
The path integral in ~k-space is written as
Z Z Y Z Y
0 0
D[Φ(~r )] → dΦ~k dΦ~∗k /N → (2/N ) d<Φ~k d=Φ~k , (12.14)
k<Λ k<Λ
7
Note that with this definition, Φ(~r ) and Φ~k have the same units. For v = 1, we have
a change of units and have to replace N by V .
208 CHAPTER 12. FIELD THEORY: GAUSS AND GINZBURG
where the factors N and 2 in the last two equations originate from the Ja-
cobians when going from real- to k-space and changing Φ~k , Φ~∗ to the real
k
fields <Φ~k and =Φ~k (note that the Jacobian for the Fourier transform is
unity only for a unitary, i.e., symmetric definition of the Fourier transform).
The (UV = ultra-violett) cutoff Λ in the product takes into account the un-
derlying lattice, and thus the coarse graining of the Hamiltonian, Λ ∼ 2π/a.
The integrals over real- and imaginary parts <Φ~k and =Φ~k are independent
but with constrained ~k vectors, hence either the product 0 is over half of
Q
~k-space or over all ~k vectors with the result be taken to the half power.
Please note that, here, we have been more careful than usual about a con-
sistent scaling between descriptions in real and k-space in order to provide
a good reference point; however, keep in mind that we have no control on
the absolute normalization in the (classical) functional integral, such that
keeping all factors may appear as an overkill as it goes beyond precision.
The Gaussian model with u = 0 in (12.4) and (12.12) can be solved exactly.
The Hamiltonian of the Gaussian model is defined as (see (12.7))
Z d
1 d r ~ r )|2 ]
HG = [rΦ2 (~r ) + c|∇Φ(~ (12.15)
2 v
1 X0 1 X0 ∗ −1 ~
= [r + ck 2 ] |Φ~k |2 ≡ Φ~k G0 (k )Φ~k
N Nβ
~k ~k
with8
kB T
G0 (~k) = . (12.16)
r + ck 2
8
The kernel G−1 0 in the Hamiltonian is the original definition of the Green function of a
free theory; it is then straightforward to show (via a Gauss integration) that the correlator
is the inverse of G−10 , hΦΦi ∼ G0 .
12.2. GAUSSIAN MODEL 209
∂ 2 h kB T dd k
Z
2πkB T i
cV = T log
∂T 2 2 (2π)d r(T ) + ck 2
dd k (T ∂T r)2
Z
∼ kB + further less singular terms,
(2π)d (r + ck 2 )2
where the dominant singular contribution is generated by the temperature
dependence in r(T ) = r0 (T − p Tc ) = r0 Tc τ . We write the denominator
r0 Tc τ + ck 2 as c(ξ −2 + k 2 ), ξ = c/r0 (T − Tc ) and switch to dimensionless
variables x = kξ to obtain
Z Λξ
4−d xd−1
cV ∼ kB ξ dx . (12.20)
0 (1 + x2 )2
For τ → 0+ the correlation length ξ → ∞ diverges, while the integral in
(12.20) remains finite for 1 < d < 4: for x → 0 we have xd−1 and for x → ∞
we obtain the scaling xd−5 . Thus, we find that cV diverges at the critical
point,
τ → 0+
cV ∼ const ξ 4−d → ∞ in 1 < d < 4. (12.21)
210 CHAPTER 12. FIELD THEORY: GAUSS AND GINZBURG
The result (12.21) defines a new value for the critical exponent α,
cV ∼ τ −(4−d)/2 ∼ τ −α ,
we find that within the Gaussian approximation
α = 2 − d/2. (12.22)
Tracing back the origin of the divergence in (12.20)
Z Z
d 1 1
cV ∼ d k −2 2 2
∼ dd k 4 ,
(ξ + k ) ξ −1 k
Next, we turn our interest to the susceptibility χ0 = ∂h hΦi and the cor-
relator G0 = hδΦ δΦi. According to (8.20) (which is generally valid for
any Hamiltonian since it is a linear response result), they are related by
χ0 = G0 /kB T . The inverse susceptibility χ−1
0 is directly obtained from the
functional derivative9
v 2 δ 2 HG
χ−1 r, ~r 0 ) =
0 (~ = v(r − c∇2 )δ(~r − ~r 0 )
δΦ(~r )δΦ(~r 0 )
1
↔ χ0 (~q ) = . (12.23)
r + cq 2
Next, let us explicitly calculate the order parameter correlator hδΦ δΦi, see
(8.20); starting from
hδΦ(~r )δΦ(~r 0 )i = hΦ(~r )Φ(~r 0 )i,
we calculate the order-parameter correlator,
1 X i~k·~r i~k 0 ·~r 0
hΦ(~r )Φ(~r 0 )i = 2 e e hΦ~k Φ~k 0 i
N
~k,~k 0
Z Y
1 0
hΦ~k Φ~k 0 i = (dΦq~ dΦ∗q~ /N ) Φ~k Φ~k 0 e−βHG . (12.24)
ZG
q~
9
Note that our measure is dd r/v rather then dd r, that generates the additional factor
v per functional derivative.
12.2. GAUSSIAN MODEL 211
For ~k 6= ~k 0 , the equation (12.24) reduces to hΦ~k ihΦ~k 0 i and hence vanishes.
For ~k = ~k 0 , we have to calculate hΦ~2 i, which will vanish by angular inte-
k
gration of dΦ~k . The only term surviving is the expectation value hΦ~k Φ−~k i,
R∞ R 2 −(β/N )(r+ck2 )|Φ~k |2
0 d|Φ~k | dϕ~k |Φ~k | |Φ~k | e N kB T
hΦ~k Φ−~k i = R ∞ 2 2 = , (12.25)
−(β/N )(r+ck )|Φ~k | r + ck 2
R
0 d|Φ ~k
| dϕ~
k
|Φ ~
k
| e
where we have used dϕ = 2π, d|Φ| |Φ| = dz 2 /2, and a standard inte-
R R R
δH
= 0,
δΦ(~r ) Φsaddle
Hsaddle = H(Φsaddle ). (12.29)
In an external field h(~r ), the saddle point10 Φmin (~r ) of H − dd r h(~r )Φ(~r )
R
δH
= h(~r ). (12.31)
δΦ(~r ) Φmin (~
r)
δ2H
Z
1
HG = H(hΦ(~r )i) + dd r dd r0 δΦ(~r ) δΦ(~r 0 ), (12.32)
2 δΦ(~r )δΦ(~r 0 ) hΦ(~r )i
δ2H
H2 ≡ , (12.33)
δΦ(~r )δΦ(~r 0 ) hΦ(~
r )i
Susceptibility (we use a general expectation value hΦ(~r )i) accounting for an
external field h(~r )):
χ−1
0 (~r, ~r 0 ) = [r + 12uhΦ(~r )i2 − c∇2 ]δ(~r − ~r 0 ) [= kB T G−1
0 ], (12.38)
−1 ~ 2 2
χ0 (k ) = r + 12uhΦi + ck ,
where the latter expression applies for a homogeneous situation. The ap-
proximation (12.37) to the Φ4 Hamiltonian can be used in both regimes
T > Tc (where hΦi = 0) and T < Tc (where hΦi = 6 0, resulting in a mass
= r + 12uhΦi2 = 2|r|). Close to Tc , we have r → 0, hΦi → 0, and the
fluctuations can no longer be treated in the Gaussian approximation.
Next, we analyse the statistical mechanics of the Φ4 model in the Gaus-
sian approximation for both regions T > Tc and T < Tc .
For large temperatures T > Tc , we use Φmin = hΦi = 0 and can make
straightforward use of the results obtained in our discussion of the Gaussian
model above, e.g., the free energy (12.19), the specific heat (12.21) with
the exponent α = (4 − d)/2, the susceptibility (12.23) and the correlator
(12.26). Analysing (12.23) as done before in Sec. 11.3 provides us with the
exponents γ = 1, ν = 1/2 and η = 0, all three unchanged when compared
to the mean-field result.11
11
One can take the point of view that the exponents for ν and η first appear in the
Gaussian model, depending on whether one wants to define the non-local χ(~ q ) within the
mean-field theory or within the Gauss theory, which naturally includes the derivative term
describing inhomogeneities in Φ.
214 CHAPTER 12. FIELD THEORY: GAUSS AND GINZBURG
At zero field, the only change then is the substitution r → r+12uhΦ2 i = 2|r|.
~ is vector-valued (with n components), there is an additional contri-
If Φ
bution to the Hamiltonian (12.39)
1 X0 ∗ p
Φ~k,i 2h u/|r| + ck 2 Φ~k,i ,
V
i,~k
dd k
Z
kB T 2πkB T
− (n − 1) d
log p (12.41)
2 (2π) 2h u/|r| + ck 2
exp α β γ δ ν η
MFA 0 1/2 1 3 1/2 0
∗ ∗
GA 2 − d/2 1/2 1 3 1/2 0
We thus see that a study of cV is best geared towards measuring the validity
divergence in GT
cV Fig. 12.3: Illustration of the var-
ious contributions to the specific
heat cV (T ) near Tc : mean-field
regular jump on top of the regular back-
background ground and the divergence due to
MF−jump Gaussian fluctuations producing
a critical exponent α = 2 − d/2.
Tc T
dd q G−1
Z
kB T 0 (q)
F = FMF + log , (12.42)
2 (2π)d 2π
with the inverse Green function G−10 (q), in a real space representation (see
−1
(12.38)) kB T G0 = (r + 12uhΦ(~r )i2 − c∇2 )δ(~r − ~r 0 ). By taking the func-
tional derivative12 of (12.42), we obtain a new expression for G−1 (after
taking derivatives, we assume temperatures T > Tc such that hΦi = 0 and
12
We use (with v = 1)
δ 2 (Φ(x))2 δΦ(x)
δ(x − y) = 2 δ(r − x)δ(x − y) = 2δ(r0 − x)δ(r − x)δ(x − y).
δΦ(r)δΦ(r0 ) δΦ(r0 )
216 CHAPTER 12. FIELD THEORY: GAUSS AND GINZBURG
kB T G−1 2
0 = r + ck ),
δ 2 βF
G−1 (~r, ~r 0 ) =
δhΦ(~r )iδhΦ(~r 0 )i
δ 2 G−1
Z
−1 0 1 0 (~
x, ~y )
= G0 (~r, ~r ) + dd x dd y G0 (~x, ~y )
2 δhΦ(~r )iδhΦ(~r 0 )i
= G−1 r, ~r 0 ) + 12u δ(~r − ~r 0 ) G0 (~r, ~r 0 ),
0 (~ (12.43)
and hence
dd q
Z
−1 ~ −1 ~
kB T G (k ) = kB T G0 (k ) + 12u G0 (~q ), (12.44)
(2π)d
dd q
Z
2 kB T
= r + ck + 12u .
(2π)d r + cq 2
dd q
Z
Σ1 = 12u G0 (~q ). (12.45)
(2π)d
The next term within a systematic expansion is given by Σ2 , see Fig. 12.4;
with the idea of resummation in the back, we only include diagrams that
cannot be divided into two by cutting a single internal line, that defines the
irreducible self-energy. The second-order term Σ2 plays an important role in
the RG of the Φ4 theory. After resummation, we can write G−1 = G−1 0 −Σ
or G = G0 /(1 − G0 Σ). Within the skeleton approximation, we replace G0
in this expression by G (and vertices by vertex functions).
The self-consistent field approximation (SCFA) picks up this idea by
replacing G0 in the loop integral Σ1 by a self-consistent expression for G. A
12.3. SELF-CONSISTENT FIELD, HARTREE APPROXIMATION 217
G G0
= + + + + ...
u
+ + ... + + ...
~ drop
~ +
1 2
= + + ...
lazy derivation goes like this: We focus on the Φ4 term of our Hamiltonian
H and contract two of the four Φ-factors H,
in one of the six possible ways (there are 4 · 3/2 = 6 ways to choose 2 of
4 factors Φ to contract). This is equivalent to replacing r by r + 12uhΦ2 i
(note the factor 1/2 in front of r) and thus we obtain the improved result
for the Green function,
Note that hΦi = 0 above Tc within the Gauss approximation (12.37), while
the square amplitude hΦ2 i does not vanish above Tc . Since
dd k
Z
hΦ2 i = G(~r − ~r 0 = 0) = G(~k )
(2π)d
dd k
Z
kB T
= , (12.48)
(2π) r + 12uhΦ2 i + ck 2
d
218 CHAPTER 12. FIELD THEORY: GAUSS AND GINZBURG
dd q
Z
kB T G−1
(~k ) = kB T G−1 ~
0 (k ) + 12u G(~q ), (12.49)
(2π)d
dd k
Z
1
χ−1 (~q = 0) ≡ r̃ = r + 12ukB T (12.50)
(2π)d r̃ + ck 2
Tc ξ d−2
0
1− ∼ Gi <0 (12.53)
Tc0 a
• Large fluctuations (at short scales) near Tc reduce the critical temper-
ature Tc0 (MF) → Tc (SCFA) < Tc0 . For d ≤ dL we have Tc = 0. For
the Φ4 -theory, we found dL = 2 (in SCFA).
as well as (1 − t)2 = (ξ0 /ξ)4 /4, we find (according to (12.58) fluctuations are
relevant when (ξ/ξ0 )d−4 < #)
ξ d−4 Ad kB 1 Ad kB 2/(4−d)
= = Gi or 1 − tG = .(12.61)
ξ0 T =TG ∆cV ξ0d 2 ∆cV ξ0d
bosons
cV fermions cV λ −transition
BCS
0 TBCS T 0 Tλ T
critical region
Above, we derived the criterion (12.61) for the width of the critical region
inside Landau mean-field theory for a second-order phase transition. More
12.4. THE GINZBURG CRITERION 223
G(~k = 0)
Z Z
1 d 1
d r hδΦ(~r )δΦ(~r )i ∼ dd r G(~r ) ∼
Vξ Vξ Vξ Vξ ξd
kB T χ(~k = 0)
∼ ∼ |1 − t|dν−γ ,
ξd
hΦi2 ∼ (1 − t)2β ∼ (1 − t)2β , (12.64)
2
to find the upper dimension dU where both quantities hδΦ i and hΦi2 scale
equally in 1 − t,
and hence
2β + γ
dU = . (12.66)
ν
For the tricritical point, we have β = 1/4, γ = 1, ν = 1/2 and dU = 3: a
tricritical point is already well described by mean-field theory in 3D.
224 CHAPTER 12. FIELD THEORY: GAUSS AND GINZBURG
Chapter 13
The elastic term tells us that the dimension of the field Φ̃, (L = length) is
225
226 CHAPTER 13. CRITICAL FLUCTUATIONS AND SCALING
R u0 6= 0, u0 1, we
For u0 = 0, this is a trivial task (Gaussian model). For
can try perturbation theory in the interaction Hint = dd x u0 ϕ4 ,
Z 1 2
Z ≈ D[ϕ]e−H0 [ϕ] 1 − Hint + Hint − ... . (13.9)
2!
Analyzing the “small” parameter u0 more accurately, we find that
(
u 4−d 0, d > 4,
u0 = 2
ξ → (13.10)
βc T →Tc ∞, d < 4,
Changing the length scale L → L0 ≡ εL, Ĝ transforms according to Ĝ0 (~k 0 ) L02 =
Ĝ(~k ) L2 (since increasing the unit leads to a reduction of the number at fixed
physical quantity) and hence
1
Ĝ0 (~k 0 ) = Ĝ(~k ). (13.15)
ε2
The dimensional analysis with just one length scale in the problem leads
to the conclusion that the mean-field theory provides us with the correct
exponents, contrary to experimental evidence. What went wrong?
The root of the problem lies in our assumption that only one length
scale determines the physics. In fact there is a second, microscopic length
228 CHAPTER 13. CRITICAL FLUCTUATIONS AND SCALING
scale Λ−1 , that influences the fluctuations in the critical region. In fact,
you might have been bothered by the appearance of the cutoff Λ in the var-
ious expressions of the thermodynamic quantities. Technically, the renor-
malization group theory will solve the problem by eliminating high-energy
short-wavelength fluctuations in a way that accounts for their presence via
simultaneously modifying the long-wavelength fluctuations by renormalising
the coupling constants r and u.
On a non-technical level, we may still ask how to set up consistent scaling
laws that fullfill (13.15) as well as (13.16) with η 6= 0? The idea is to consider
the second length scale a = Λ−1 explicitly. For example, replacing (13.16)
by
Ĝ(~k, Tc ) ∼ aη k −2+η , (13.17)
we can consistently modify the k-scaling of Ĝ under a change in length scale
L0 = εL, with the overall scaling Ĝ0 ∝ Ĝ/ε2 preserved,
a η
Ĝ0 (~k 0 , Tc ) ∼ a0η k 0−2+η ∼ (εk)−2+η
ε
1 1
∼ 2 aη k −2+η ∼ 2 Ĝ(~k, Tc ).
ε ε
By studying Ĝ at Tc and constant ~k as a function of a, we find the anomalous
exponent η from Ĝ ∼ aη .
Similarly, we can find the anomalous dimension of the correlation length
ξ: Using [ξ] = L, [a] = L, [r̃] = 1/L2 , we modify the scaling law for ξ by
introducing a scaling function f (x),
1
ξ = √ f (r̃a2 ), (13.18)
r̃
with a dimensionless argument x. If f (x → 0) is regular, we again obtain
ν = 1/2, the mean-field value. If, however,
f (x → 0) ∼ xθ , (13.19)
with the canonical dimension 1/2 and the anomalous dimension θ of ξ.1 In
the same way one says that the order-parameter field
Φ̃ ∼ L1−d/2−η/2 (13.21)
1
Alternatively, we may assume that the correlation length ξ has to be built from c/r
and from a, ξ 2 ∼ (c/r)1−2θ a2θ , from which we find ξ ∼ (c/r)1/2−θ ∼ |τ |−1/2+θ and hence
ν = 1/2 − θ.
13.1. SCALING LAWS 229
The scaling laws carry the names of their inventors Fisher, Rushbrooke,
Widom, and Josephson and read
(F ) : γ = ν(2 − η),
(R) : α + 2β + γ = 2,
(W ) : γ = β(δ − 1),
(J) : νd = 2 − α. (13.22)
dd r G(~r ) as well as β dd r HΦ = 1 or
R R
(note that we have used kB T χ =
Φ = −∂H f ) and the definitions
β(δ − 1) = γ. (13.27)
2
Note that we assume two anomalous scalings where the microscopic scale a has entered,
one for lengths, i.e., ξ (implying ν 6= 1/2 is non-trivial) and one for the field Φ (implying
η 6= 0 is non-trivial). Near the critical point, the microscopic length a, that was needed in
getting nontrivial exponents for distances r and fields Φ, only appears via the anomalous
values of ν and η; it does not appear again in other quantities due to the assumption of ξ
being the only relevant scale. Indeed, since all quantities can be built from lengths/volumes
and fields, assuming a scaling for lengths/volumes L ∼ ξ, V ∼ ξ d and fields Φ ∼ ξ 1−d/2−η/2
does the job. Still, one may ask whether this assumption of two anomalous exponents
(and the appearance of ξ alone as the relevant length near Tc ) is correct. E.g., one could
wonder, whether another quantity Q involves yet one more anomalous scaling exponent,
i.e., follows a scaling Q ∼ ξ y−κ aκ with anomalous exponent κ rather than Q ∼ ξ y . This
question has to be answered by calculation, e.g., using RG. The question is, whether
the integrals appearing in the perturbative analysis are determined by small wavelengths
(k ∼ 1/ξ) or whether large k ∼ 1/a-values play a role.
3
The term −2 in the exponent |τ |dν−2 originates from the derivative ∂T ξ ∝ τ −ν−1 .
13.1. SCALING LAWS 231
β(δ + 1) = νd = 2 − α → α + 2β + γ = 2.
The following consideration gives a deeper insight into the scaling laws and
is due to Widom: To fix our ideas, we assume a magnetic system and analyse
the set of curves m(τ, h) (Φ → m, τ = (T − Tc )/Tc , βH = h conjugate to
m). Widom then suggested combining the two relations
(
0, τ > 0,
m(τ, 0) ∼ β
(13.28)
±|τ | , τ < 0,
and
m(0, h) ∼ ±|h|1/δ , τ = 0, (13.29)
into a single scaling law
h
m(τ, h) ∼ |τ |β M± . (13.30)
|τ |∆
In doing so, the exponents β and ∆, the so-called “gap exponent”, are taken
to be universal and M± are two scaling functions, + for τ > 0, − for τ < 0.
The ansatz (13.30) tells us that by choosing the three parameters β, ∆ and
Tc correctly, we can reduce all the curves m(τ, h) to only two curves, one
describing the situation above and one below Tc . This is an elegant ansatz
that allows
• to reduce the experimental data to two curves and thus represent the
data in a consistent manner, and
A well known example is the data collapse at the magnetic transition in Ni,
see Fig. 13.1
The scaling functions M± (x) fullfill certain conditions, e.g., they are odd
under a sign change in x,
(this follows from the symmetry m(τ, −h) = −m(τ, h)) and
M+ (0) = 0, M− (0) 6= 0,
232 CHAPTER 13. CRITICAL FLUCTUATIONS AND SCALING
4 8 12 16 20 ·10 5
∂m 1 ∂m |τ |β
χ(h = 0) = = ∼ ∆ M0± (0) ∼ |τ |−γ
∂H kB T ∂h |τ |
→ ∆ = β + γ. (13.31)
M± (x → ∞) ∼ xλ (13.32)
Correlator:
So far, we have used only the scale ξ in the problem (but accounted for the
existence of the second, high energy cutoff scale Λ ∼ π/a (or lattice cutoff
a) by assuming an anomalous contribution θ to the canonical dimension
1/2 for the exponent ν, see Eq. (13.20)). When discussing the correlator at
T = Tc , we have to include the microscopic scale explicitly when analyzing
the scaling of the correlator G (since ξ → ∞ and we need the second scale
a to have an anomalous exponent η). Let a be the microscopic length; we
distribute the scaling factor L2 over a and ξ,
G(x → ∞) ∼ xε , (13.46)
we should have
(the factor a2−s only corrects the length scaling and is otherwise irrelevant
in this context) and we obtain Fisher’s law
(2 − η)ν = γ (F ). (13.48)
Thus, we have rederived all the scaling laws assuming a scaling ansatz for
various quantities. Note, that all the above discussion is based on the scaling
hypothesis that has an empirical root. A firm derivation of scaling is only
provided by the RG treatment introduced in the following chapter.
Chapter 14
Renormalization Group
Theory
Consider an Ising model, see Eq. (8.4), on a d-dimensional lattice ~ri with
lattice constant a (we set µ = 1 and drop the superscript on Siz ) and N
sites,
X X
H = −J Si Sj − h Si . (14.1)
hi,ji i
Near the critical point, the correlation length ξ diverges and spins are
strongly correlated, allowing us to go over to a coarse grained description
with averaged spins
X
S̄I = b−d Si , (14.2)
i∈BI
with the Block BI containing bd spins, see Fig. 14.1x (in order to keep
the spins normalized, we may replace S̄I → S I = S̄I /|S̄I |, that implies a
235
236 CHAPTER 14. RENORMALIZATION GROUP THEORY
2a
Fig. 14.1: Block-spin transformation: The spins on the original lattice ~ri
with lattice constant a are first combined into block spins S̄I . In a second
step, both the spins and the lattice are rescaled, S̄I → s0 and R~ I → ~r 0 .
somewhat different scaling analysis for G, see below). The new lattice with
positions R~ I now has a lattice constant ba, fewer particles b−d N , and a
correlation length that is reduced by b, ξ → ξ/b. Since the new system has a
shorter correlation length it has moved further away from the critical point.
The above transition to Blockspins goes under the idea of ‘integration’
over small-scale or high-energy degrees of freedom. Given the scale invari-
ance of the problem near a critical point, this ‘integration’ shall not change
the form of the theory, the Hamiltonian, or other quantities, such as the
correlator. However, having integrated over small scale degrees of freedom,
all our quantities now live on the scale ba. Instead of keeping track of the
scale, we simply rescale all quantities, which is the second idea or concept
in building the renormalization group.
Before doing that, let us agree on the following notation to denote coor-
dinates and fields after integration and rescaling,
~ I ↔ ~r −→ ~r 0 = ~r/b,
R
S̄I ↔ s(~r ) −→ s0 (~r 0 ) = bλs s(~r ), (14.3)
λs = (d − 2 + η)/2 (14.7)
as it allows us to include the field rescaling b2λs into the decay |~r1 −~r2 |−(d−2+η) ,
Note the difference to Eq. (13.4) providing the trivial field scaling [Φ̃] =
L1−d/2 —the Block spin s0 is already assumed to scale with an anomalous
correction η that is different from the original ‘geometric’ scaling of the field
s.2
Next, let us analyze the behavior of the conjugate field h under coarse
graining. Again, we find the scaling of h(~ri ) (assumed to be smooth
P on the
scale ba) from the requirement that the free energy term Hh = − i h(~ri )Si
remains form-invariant under the transformation to block spins,
X X X
−Hh = h(~ri )Si = ~ I )bd S̄I =
h(R h(~r )bd s(~r ) (14.9)
i ~I
R ~
r
!
X X X
= h(~r )bd−λs s0 (~r 0 ) = h(b~r 0 )bd−λs s0 (~r 0 ) = h0 (~r 0 )s0 (~r 0 ),
~
r r0
~ r0
~
Hence, the scaling exponent for the field h is (see Eq. (13.45))
d+2−η β+γ
λh = d − λs = = , (14.11)
2 ν
in agreement with the gap exponent ∆ = λh ν = β + γ found in Eqs. (13.44)
and (13.36). In (14.11), we have made use of the scaling relations (??),
(d + 2 − η)/2 = βδ/ν, and Widom’s scaling βδ = β + γ.
1
If we want to rescale the normalized field S I , we make the Ansatz that |S I | ∼ b−λs
and G0 (~r 01 , ~r 02 ) ∼ (1/|S I |)2 G(~r1 , ~r2 ).
2
Technically, a finite anomalous exponent η is produced as a result of renormalization
of the couplings in the Hamiltonian when changing the scale parameter b; e.g., in the
Φ4 -theory of chapter 16.2, a nonzero value for η is obtained at order 2 , with = 4 − d.
238 CHAPTER 14. RENORMALIZATION GROUP THEORY
λχ = d − λφ . (14.12)
This relation is helpful in the discussion of the scaling of the reduced Hamil-
tonian KJ = HJ /T . Expanding around the critical point, we write
X
KJ ≈ K c + τ K c s(~r )s(~r + d~ )
r,d~
~
X
= Kc + Kc τ (~r ), (14.13)
~
r
then assumes the role of the field conjugate to τ , and imposing the scaling
ansatz for and τ ,
λ = d − λτ . (14.16)
From (13.43), we know that λτ = 1/ν (since τ scales as ξ −1/ν ) and we obtain
the scaling exponent for the energy
dν − 1
λ = . (14.17)
ν
Note that we cannot find the scaling of via field rescaling since we cannot
disentangle the interactions between spins when going to a block spin. Also,
on a more abstract level, the scaling of energies is related to ν with its
own anomalous contribution θ, rather than field scaling with the anomalous
exponent η.
Finally, we can find the homogeneity relation for the free energy F . De-
termining the latter by calculating the partition sum over all configurations,
we expect to arrive at the same expression independent of the starting point,
the original or the block variables, hence F (τ, h, χ, . . . ) = N f (τ, h, χ, . . . ) =
N 0 f (τ 0 , h0 , χ0 , . . . ). Using N = bd N 0 as well as the scaling properties of the
various fields, we obtain the scaling or homogeneity relation for the free
energy density
This relation is true for any value of b; choosing b = |τ |−ν (i.e., rescaling by
the maximal value ξ/a) we recover the scaling relation (13.44)
f (τ, h, χ, . . . ) = ξ −d f0 (h/|τ |∆h , χ/|τ |∆χ , . . . ) (14.19)
with the gap exponents related to the scaling exponents via ∆ = λν.
This is a very deep result as it tells us about the relevance of terms in
the Hamiltonian: If the scaling exponent λ or gap exponent ∆ of a variable
v is negative, λ, ∆ < 0 (note that ν > 0) then the variable is reduced
under rescaling and vanishes ultimately, v → bλ v → 0 under rescaling or
v/|τ |∆ = v|τ ||∆| → 0 on approaching the critical point τ → 0. Accordingly,
variables are classified as relevant if λ > 0, irrelevant if λ < 0, or marginally
relevant if λ = 0. This is the origin of universality: although Hamiltonians
may look different with additional terms or terms of different forms, once
they approach the same form under rescaling (i.e., differences are reduced
and finally vanish under rescaling) they describe the same physics.
The above discussion was simple and illustrative but a bit naive. E.g., it is
rather clear that the transition to block spin variables cannot be accurately
done without generation of new terms in the Hamiltonian that differ from
the simple nearest-neighbor Ising form. In order to overcome this deficiency,
we consider a most general Hamiltonian H[{K}] that depends on a complete
set of coupling constants {K}. In fact, to start with, the Hamiltonian may
involve some subset of finite couplings, but new terms may be generated in
the step of renormalization, therefore the set {K} shall be complete. E.g.,
for an Ising spin system such a Hamiltonian has the form
X
H[{K}] = Kα Sα (14.20)
α
with
Y
Sα = Si (14.21)
i∈Iα
and Iα denoting arbitrary sets of site labels. The Ising Hamiltonian (14.1)
then contains only two types of sets, the single sites with K1 = −h and
the set of nearest neighbor pairs with K2 = −J. The general Hamiltonian
H[{K}], however, is supposed to involve all possible sets
X X
H[{K}] = K0 N + K1 Si + K2 Si Sj (14.22)
i hi,jinn
X X
+ K20 Si Sj + K3 Si Sj Sk + . . .
hi,jinnn hi,j,kinn
240 CHAPTER 14. RENORMALIZATION GROUP THEORY
The projection operator in (14.26) plays the role of the filter in Eq. (11.2)
when going from microscopic to coarse grained variables. Furthermore, we
require that P ({Si }, {SI0 }) respects the symmetries of the problem in order
0 0
P K[K , {SI }] involves
to make sure that only symmetry-allowed terms. Using,
the property {S 0 } P ({Si }, {SI0 }) = 1 (since a given set {Si } defines the set
I
{SI0 }, try it out with 6 spins and 2 blocks), the partition sum is conserved,
0 0
ZN 0 [K 0 ] = Tr{SI0 } eK[K ,{SI }] (14.27)
= Tr{SI0 } Tr{Si } P ({Si }, {SI0 })eK[K,{Si }]
= Tr{Si } eK[K,{Si }] = ZN [K],
K ∗ = Rb [K ∗ ]. (14.29)
In addition, we will look at the flow around these fixed points which is
determined by the linearized equation around each K ∗ , see below—knowing
the fixed points and the flow in their vicinity will give us a good picture
of the overall flow of couplings and will provide us with the major physical
insights into the system’s behavior.
Let us first consider the scaling of lengths. As explained above, the RG
takes the system further away from a critical point since lengths shrink by
b and thus
Hence, all the K points in the basin of attraction are critical; they define
the critical manifold associated with the fixed point. This is the basis for
universality—all the different Hamiltonians in the basin of attraction of a
given critical fixed point exhibit the same critical behavior characterized by
the same exponents.
Next, let us find these exponents. In order to do so, we have to go
slightly away from the critical point and determine the flow in its vicinity.
Hence, we expand K around K ∗ ,
with
∂Rb [K]n
Rnm [K ∗ ] = . (14.34)
∂Km K=K ∗
Equation (14.33) is the linearized RG equation describing the flow near the
fixed point K ∗ . Assuming that Rnm is symmetric (which is often, but not
3
A point K is in the basin of attraction of K ∗ if limn→∞ Rbn [K] → K ∗ .
14.3. FIXED POINTS AND EXPONENTS 243
always, the case, see later about strategies for non-symmetric situations),
we diagonalize the matrix
Rnm [K ∗ ]vm
α
= bλα vm
α
(14.35)
with coefficients
X
cα = vnα δKn = v α · δK, (14.38)
n
where in the last equation, we have to remind that δK and v α are vectors.
Equation (14.37) is not valid for non-symmetric matrices Rmn .
Applying now the RG transformation to δK expressed in the form (14.37),
we immediately see that the coefficients get multiplied by the eigenvalues
bλα ,
K3 critical
fixed point
critical manifold
K1
trajectories / flow K2
under RG
Fig. 14.2: Sketch of the critical point with its critical manifold and flow
under the renormalization group in parameter space {K}. Points on the
critical manifold flow into the fixed point and exhibit the same universal
behavior. Moving away from the critical manifold, trajectories first approach
the fixed point and then move away along an unstable direction. This flow
is governed by the same eigenvectors and eigenvalues for all points near
the critical manifold, thus generating the universal behavior with identical
exponents.
Fixed points then can be classified through their flow pattern: A fixed
point associated with a bulk phase is a pure sink, with the couplings K ∗
characterizing the bulk phase. A critical fixed point exhibits both inward
and outgoing flows. The inward flows define the critical manifold, while
the outgoing flows in the relevant directions typically take the system away
from the fixed point and drive it into a stable bulk phase. Hence, all relevant
couplings have to be carefully tuned in order to ultimately reach the desired
endpoint of the trajectory, be it a critical point or a bulk phase. A generic
situation of the flow in coupling space is shown in Fig. 14.2; we will encounter
specific flow patterns in our applications in sections 15–17.
A final remark concerns the notion of universality: All couplings in the
critical manifold will ultimately flow to the fixed point K ∗ . Going slightly
off this manifold will lead to trajectories that first drive the system towards
the fixed point but then will turn away from K ∗ as the relevant directions
14.3. FIXED POINTS AND EXPONENTS 245
start to kick in. This drive to and away from the fixed point is governed
by the same exponents λα characterizing the fixed point K ∗ , which tells us
that the behavior of all these systems is universal, i.e., governed by the same
exponents or fixed point.
where vLα is a row vector. It is a few steps to show that left- and right-vector
eigenvalues are the same and that the eigenvectors are orthogonal, vLα · vRα =
δα,α0 , assuming proper normalization. The only change then appears in the
expansion of δK which is done in the eigenbasis vRα ,
X
δK = cα vRα (14.41)
α
Expanding in ε, we obtain
dKb
= B[Kb ] = ∂ε R1+ε [Kb ]|ε=0 . (14.45)
d log b
In section 8.6, we have discussed the exact solution for the 1D-Ising model.
This solution provides an ideal example for the illustration of the Kadanoff
block spin transformation as it can be carried out easily (and exactly) with
the help of transfer matrices. Before doing so, we complete the exact solution
with a calculation of the spin correlations in this model. Furthermore, we
will discuss the results of section 8.6 in the light of our expanded knowledge
on critical phenomena.
The Hamiltonian for the Ising model is given in Eq. (14.1). Due to trans-
lation invariance, the spin correlator G(i, j) = G(n) = hSi Sj i only depends
on the distance n = |i − j| between sites. We define the dimensionless cou-
pling K = J/kB T , set h = 0, and use the transfer matrix technique to obtain
247
248 CHAPTER 15. REAL SPACE RG, SPIN MODELS
(choosing i = 1 and j = n + 1)
1 X
G(n) = S1 hS1 |T |S2 ihS2 |T |S3 i . . . Sn+1 hSn+1 |T |Sn+2 i . . . hSN |T |S1 i
ZN
{Si }
1 X
= S1 T n Sn+1 T N −n
ZN
S1 ,Sn+1
1 h i
= (2 cosh K)N (tanh K)n [1 + (tanh K)N −2n ] . (15.1)
ZN
Here, we have used that T (S, S 0 ) = cosh K (1 + SS 0 tanh K) and
With ZN = (2 cosh K)N and taking the limit N → ∞ (note that tanh K < 1
at any finite temperature, hence the term ∝ (tanh K)N vanishes), we find
that
with
e−K
K+L
eC
T (K, L, C) = e , (15.11)
e−K eK−L
T (K 0 , L0 , C 0 ) = T b (K, L, C) (15.12)
u(1 + v)2
u0 = , (15.15)
(u + v)(1 + uv)
v(u + v)
v0 = ,
1 + uv
w2 uv 2
w0 = .
(1 + v)2 (u + v)(1 + uv)
Let us analyze the resulting flow, first for L = 0, see Fig. 15.1. Eq. (15.13)
defines the fixed points
All values K < ∞ flow to K0∗ = 0 which thus is a stable or attractive fixed
point—it describes the high-temperature paramagnetic phase of the Ising
model. Only the point K∞∗ = ∞ has tanh K ∗ = 1 and hence maps to itself.
∞
∗ = 0 unstable fixed point describing the ferromagnetic phase
This is the T∞
of the Ising model. There is no finite temperature phase transition in the
model. Combining (15.13) and the definition of ξ in (15.4), we obtain the
scaling
ξ 0 = ξ/b, (15.17)
unstable stable
fixed point fixed point
F
1 h =0
h >0
v = e −2L
paramagnetic line
isolated
fixed point
0
0 u = e −2K 1
T 0 T 8
ξ ξ 0
8
1
gsing (ξ, L) = g0 (Lξ), (15.20)
ξ
For completeness, we may also study the flow near the other fixed points.
Near u∗ = v ∗ = 0, the small deviations u and v flow as u0 = u/(u + v) and
v 0 = v(u + v), hence a point (u, 0) ends up at u0 = 1, while a point (0, v)
goes to v 0 = v 2 , in agreement with the flow sketched in Fig. 15.1. Near
u∗ = 1, v ∗ = v, we find that a deviation δ < 0 from u∗ = 1 scales as
δ 0 = −δ 2 v/(1 + v)2 while v changes to v 0 = v[1 + δ(1 − v)/(1 + v)], hence
the flow is downwards towards the paramagnetic line.
a
1 2 Fig. 15.2: 2D Ising model
3 on a triangular lattice. Blocks
with three spins form a new
3a lattice √
with larger lattice con-
stant 3a. The interaction
1 2 between blocks involves two
bonds 13 and 23.
3
The 2D Ising model has been exactly solved by Onsager and provides
us with the exacts exponents α = 0 and δ = 15—we can use these results
and determine the other exponents, allowing us to check the quality of our
RG analysis. From Josephson’s law 2 − α = νd, see (13.22), we find ν = 1.
Using Widom’s law, we obtain δ = (β + γ)/β. Using Eq. (14.11), β + γ =
(ν/2)(2 + d − η) as well as (13.26) and Fisher’s law to find 2β = ν(d − 2 + η),
15.2. 2D ISING MODEL 253
we divide the lattice into blocks BI and define the new block spins SI through
the majority rule,
where the blocked variables (S)I = (SI,1 , SI,2 , SI,3 ) group the original spins
Si into blocks, SI,i denoting the i-th spin in the I-th block. The set {(S)I }
then covers all spin configurations of the original lattice. We can relate
the coarse-grained Hamiltonian for a specific block configuration {SI } by
summing the original Hamiltonian over all configurations {(S)I } compatible
with {SI },
X
eKb [{SI }] = eK[{SI },{(S)I }] . (15.25)
{(S)I }
In order to evaluate the partial sum in Eq. (15.25), we split the Hamiltonian
into a part K0 involving only bonds within blocks and a part K0 accounting
for bonds connecting different blocks,
XX
K0 = K SI,i SI,j , (15.26)
I (S)I
X X
0
K = K SI,i SJ,j . (15.27)
I6=J (S)I ,(S)J
with
0 0 02i 0 i2 ]/2+...
heK i0 = ehK i0 +[hK 0 −hK 0 (15.29)
and
K0 [{SI },{(S)I ]}
P
{(S)I } A[{SI }, {(S)I ]}]e
hA[{SI }]i0 ≡ P K0 [{SI },{(S)I ]}
. (15.30)
{(S)I } e
with
X
ZI (K) = eK(SI,1 SI,2 +SI,2 SI,3 +SI,3 SI,1 ) = 3e−K + e3K (15.32)
(S)I
the partition function of the block I which only depends on K and N/3 is
the number of blocks.
We approximate the first factor in Eq. (15.28) by its lowest-order term
hK0 i0 in the expansion (15.29) which involves just two bonds between adja-
cent blocks, see Fig. 15.2,
hK0 (SI , SJ )i0 = KhSI,1 SJ,3 + SI,2 SJ,3 i0 = 2KhSI,1 SJ,3 i0 . (15.33)
The Hamiltonian K0 does not couple spins in different blocks and hence the
correlator in (15.33) factorizes
e−K + e3K
hSI,1 i0 = SI ≡ s(K) SI (15.36)
3e−K + e3K
with the field scaling factor s(K) ≤ 1. Hence, the term hK0 i0 assumes the
simple form
with rescaled block spin variables and the rescaled Hamiltonian reads
X
Kb [{SI }] = (N/3) ln ZI (K) + K 0 SI SJ + O(K2 ). (15.38)
hI,Ji
Hence, we find that the RG has generated a constant term (trivial energy
shift) and has renormalized the couplings between spins
K 0 = 2[s(K)]2 K. (15.39)
into Eq. (15.41) and evaluating the right hand side to lowest (i.e., zeroth)
order in K0 , we obtain
X X X
Lb [{SI }] = L hSI,i i0 = 3LhSI,1 i0 . (15.43)
I (S)I I
Making use of Eq. (15.36), we obtain the recursion relation for the magnetic
field
L0 = 3[s(K)] L. (15.44)
256 CHAPTER 15. REAL SPACE RG, SPIN MODELS
The RG equations
tell us all universal features of the 2D Ising model. We have 5 fixed points,
√
K ∗ = 0, ∞, s−1 (1/ 2) at L∗ = 0 (15.46)
L∗ = ±∞ at K ∗ = 0,
where s−1 is the inverse function of s, i.e., the finite K (and non-trivial or
critical) fixed point is defined by [s(K ∗ )]2 = 1/2. We first focus on this
critical fixed point. Solving for K ∗ ≡ Kc , we obtain the location of critical
fixed point at1
√
Kc = (1/4) log(1 + 2 2) ≈ 0.3356, (15.47)
L critical
h >0
fixed point
paramagnetic
continuity T 0
P F
fixed point T
8
C K ferromagnetic
discontinuity
fixed point
h <0
first−order
transition
Fig. 15.3: Sketch of the flow diagram for the 2D Ising model. The 2D
Ising model has a finite temperature critical point C and two relevant flow
directions. At h = 0 (an invariant manifold), the flow is towards the low-
temperature ferromagnetic critical point when starting above Kc and to-
wards the high-temperature paramagnetic fixed point when starting below
Kc . The ferromagnetic fixed point is a discontinuity fixed point describing
the first-order phase transition at T < Tc , h = 0 of the generic magnetic
phase diagram in Fig. 8.3 (discontinuous order parameter jumping between
T = 0 bulk ferromagnetic phases with spontaneously broken symmetry with
spins/magnetization pointing up or down). The paramagnetic point is a con-
tinuity fixed point; both F and P are unstable towards the sinks S describing
field-aligned spin-ordered phases (without spontaneously broken symmetry).
actually attractive along K. Another speciality of this fixed point is the fact
that the scaling exponent λh = d. This result tells us that this fixed point
describes a first-order transition (Nienhuis-Nauenberg criterion) in h = 0 at
temperatures T < Tc —this is the usual first-order transition in the magnetic
phase diagram of a ferromagnetic system, see Fig. 8.3. Finally, the fixed
points at K ∗ = 0, L∗ = ±∞ describe high-temperature phases with spins
aligned by the strong magnetic field—they are sinks and hence describe bulk
(ferromagnetic) phases. A sketch of the overall flow is given in the phase
diagram of Fig. 15.3.
258 CHAPTER 15. REAL SPACE RG, SPIN MODELS
Chapter 16
k-space renormalization of
field theories
+ · · · − hΦ~k=0 . (16.1)
Alternatively, we replace
dd k
Z Z
1 X
→ ≡ (dd k) (16.2)
V (2π)d
~k
and cut the integral at Λ ∼ 2π/a with a the lattice constant. The central
quantity to be determined is the partition function
Y0 Z
Z(c, r, u, . . . ) = e−F (c,r,u,... ) = N dΦ~k dΦ~∗k e−K[{Φ~k }] , (16.3)
k<Λ
with the prime reminding us that the product is to be taken over half of
~k-space since Φ ~ = Φ~ for a real order parameter field Φ(~r ). The prefactor
−k k
N reminds us that we actually have an undefined normalization in the mea-
sure of the path integral; this normalization drops out in the calculation of
259
260 CHAPTER 16. K-SPACE RG
1 Λ/b d
Z
Kb [{Φ~k }] = (d k)[c0 k 2 + r] |Φ~k |2 + . . . (16.5)
2
depends on a reduced set of modes with k < Λ/b and in general picks up a
new elasticity c0 6= 1 which we will deal with in step 3).
2) rescaling: we restore the old cutoff Λ by increasing the unit of length
by a factor b (as we did before with the block spin transformation), r →
r0 = r/b and hence
k 0 = bk. (16.6)
3) Field rescaling: we restore the unit elasticity (and hence the original
form of the Hamiltonian) by absorbing the factor c0 /bd+2 in the field Φ,
Defining the scaling factor zk for the field Φ~k through Φ~0 0 = zk Φ~k , we find
k
that
zk = 1/b(d+2)/2 (16.9)
Eq. (16.8) corresponds to the introduction of new block spin variables in the
Kadanoff procedure used before. Here, we have used the simpler version of
16.1. GAUSS MODEL 261
1 Λ d 0 02
Z
0
Kb [{Φ~k 0 }] = (d k )[k + r0 ] |Φ~0k 0 |2 + . . . (16.10)
2
with the renormalized ‘mass’
b2
r0 = r. (16.11)
c0
Hence, we have absorbed the volume bd and have shifted the renormalization
of the elasticity to that of the mass by rescaling the field Φ → Φ0 in order to
arrive at a form invariant Hamiltonian. The partition function takes again
the form (16.3) but with a new (physically irrelevant) prefactor. In the
following, we will apply this scheme to two examples, the Gauss model and
the Φ4 theory.
in real space,
Z Λ/b Z Λ
d i~k·~
r ~
Φ(~r ) = (d k) Φ~k e + (dd k) Φ~k eik·~r (16.14)
0 Λ/b
= Φ (~r ) + Φ> (~r ).
<
262 CHAPTER 16. K-SPACE RG
The Hamiltonian then splits into slow and fast modes as well,
Z Z
1 d 2 < 2 1
K = (d k)[r + k ]|Φ~ | + (dd k)[r + k 2 ]|Φ~> |2 − hΦ~k=0
2 k 2 k
< >
= K +K , (16.15)
where we have included the magnetic field term −hΦ~k=0 into K< . The
partition function separates into factors describing slow and fast modes as
well (we abbreviate DΦ~k = dΦ~k dΦ~∗ ),
k
Z Y
0
ZG (r, Λ) = e−FG = DΦ~k e−K (16.16)
k<Λ
Z Y Z Y
0 < 0 >
= DΦ~k e−K DΦ~k e−K = ZG (r, Λ/b)ZG>
k<Λ/b Λ/b<k<Λ
with
h1 X 2π i
ZG> = exp ln 2 (16.17)
2 k
Λ/b<k<Λ
and where we have set r = 0 in the second, regular factor ZG> . Focus-
ing on the critical behavior of the Gauss model, we can ignore this factor.
The partition function Z G(r, Λ/b) can be brought into the form ZG (r, Λ) by
rescaling momenta k 0 = bk and fields Φ~0 0 = (1/bd+2 )1/2 Φ~< ~ 0 . These steps
k k=k /b
2) and 3) of the general RG scheme then provide us with the renormalized
Hamiltonian
1 Λ d
Z
0
Kb [{Φ~k }] = (d k)[k 2 + r0 ] |Φ~0k |2 − h0 Φ0 (0) (16.18)
2
and the renormalized couplings
where we have used the fact that the integration over modes Λ/b < K < Λ
does not renormalise the elasticity, c0 = 1, a consequence of the factorization
of the modes. Choosing b = 1 + ` close to unity, we can expand, r` =
(1 + `)2 r ≈ r(1 + 2`), to obtain the differential scaling laws
Note that ` ≈ log b; the differential equations (16.20) give us direct access
to the scaling exponents λ, in agreement with (14.45).
The recursion relations (16.19) have two fixed-points at r∗ = 0, ∞, of
which the Gaussian one at r∗ = 0, h∗ = 0 is the interesting one. This is
a critical fixed-point describing the critical behavior of the phase transition
16.1. GAUSS MODEL 263
We ignore the factor ZG> (it could be absorbed into the normalization factor
N of the path integral) and define
with M (Λ) the number of modes in the integration volume k < Λ. Compar-
ing the right hand sides of the above two equations, we obtain the scaling
relation
or
i.e., we can describe the physics originating from m0 and Λ (up to an ir-
relevant shift in energy ∝ ln zs ) through the physics described by a new
partition function (or free energy FG ) that lives on the scale λ with a cou-
pling constant
λ2
m(λ) = m0 (16.27)
Λ2
that depends on the scale (running coupling). All we know when looking at
the world with a resolution λ is the associated mass m(λ). From which scale
Λ the ‘microscopic’ mass m0 originated is neither relevant nor accessible.
The ensuing free energy then is perfectly regular, see Eq. (12.19),
16.2 Φ4 theory
We have already seen, that the Gauss model does not allow us to go across
the critical point at Tc of the Landau-Ginzburg free-energy functional and
we need to include the regularization property of the Φ4 term in order to do
so. Accounting for this term is non-trivial and quite some work, as modes
now start interacting.
Let us first quickly check the relevance of this term. In order to do
so, we consider the trivial scaling of the various terms in the real space
version of K, see (16.1). Rescaling lengths and fields, r → r0 = r/b and
φ(~r ) → Φ0 (~r 0 ) = zr (b)Φ(~r ), we find that the mass r, elasticity c, interaction
u, and field h energies scale as
Assuming that the elastic term does not rescale, we set zr2 = 1/b2−d (note
the difference to the field scaling (16.8) in k-space) and obtain the trivial
scaling relations (due to power counting)
in agreement with (16.19). Thus the coupling u is relevant for d R< 4 and ir-
relevant for d > 4. Higher order contact couplings of the form up dd rΦp (~r )
scale with
and hence powers p > 4 are less dangerous, i.e., they become irrelevant
already at smaller dimensions d, e.g., the Φ6 is irrelevant for d > 3. Hence,
in our physical dimension d = 3, we are well advised to account for the Φ4
term and may be more casual about the Φ6 term in the Landau-Ginzburg
expansion. As a general rule, we learn that every field Φ adds a power
b1−d/2 and every gradient comes with a factor b−1 in the power count, so
more fields Φ and more gradients ∇~ render a term less relevant. Below, we
will consider dimensions
=4−d (16.32)
in our -expansion with assumed to be small, hence we are safe from the
Φ6 term for < 1.
Next, we outline the strategy of the calculation. We start from the
Hamiltonian (16.1) at h = 0, written in the form
Z Λ
1 2
(dd k) ck + r |Φ~k |2
K(c, r, u) = (16.33)
2
Z Λ
+ Π4i (dd ki )u(2π)d δ(Σ4i ~ki )Φ~k1 Φ~k2 Φ~k3 Φ~k4 .
We split the modes into slow and fast ones as in (16.14) and divide the
Hamiltonian as in (16.15),
Note that K< and K> are simple Gauss Hamiltonians over the slow and fast
fields.
266 CHAPTER 16. K-SPACE RG
with Z > given in (16.17) and the expectation value h. . . i> to be taken with
the Gaussian K> describing fast modes. We drop Z > (as it does not interfere
with the critical behavior) and rescale with k 0 = bk and Φ~0 0 = zk Φ~<0 to
k k /b
obtain (see (16.23))
In lowest order (i.e., 1 ), the elasticity does not renormalize (see the more
detailed discussion later) and the field rescales with zk = 1/b(d+2)/2 such as
to keep c0 = c (see Eq. (16.9)). We derive the RG equations for r and u to
order = 4 − d, slightly away from the point d = 4 where effects from u
become irrelevant,
where the first terms follow from simple power counting. Including cor-
rections to order u2 or 2 in the equation for u, we will find a non-trivial
fixed-point u∗ ∼ , hence r∗ ∼ as well, and we do not need corrections ∼ u2
in the equation for r. In general, the RG will generate additional terms ∝ Φp
with p > 4 as well as dispersive couplings, e.g., a dispersive elasticity c0 (k)
is obtained in order 2 and will generate a finite η ∝ 2 ; we will not include
such higher-order in corrections in our lowest-order analysis below.
16.2. Φ4 THEORY 267
Next, we calculate the cumulants hVi> and [hV 2 i> − hVi> 2 ]. We start
out with the linear correction hVi> and insert the decomposition into slow
and fast modes Φ(~r ) = Φ< (~r ) + Φ> (~r ), Eq. (16.14), into the definition of
V. As a result, we obtain terms such as Φ< 4 , Φ< 3 Φ> , Φ< 2 Φ> 2 , Φ< Φ> 3 ,
Φ> 4 . The first term generates the trivial scaling term b−3d z −4 Φ0 4 = b4−d Φ0 4
with z −4 = b2d+4 and the factor b−3d deriving from the k-integrations (b−4d )
combined with the δ-function (bd ). The second and fourth terms are odd
in Φ> and hence vanish when averaging. The last term provides a constant
shift of the free energy which we ignore. Remains the terms Φ< 2 Φ> 2 which
has the form
Z Λ/b
>
hVi |Φ< 2 Φ> 2 = u (dd k1 )(dd k2 ) Φ~< Φ~< (16.40)
k1 k2
Z Λ
× (dd k3 )(dd k4 ) (2π)d δ(Σ4i ~ki )hΦ~> Φ~> i> .
k3 k4
Λ/b
and hence
Z Λ/b
>
hVi |Φ< 2 Φ> 2 = uI1 (dd k1 )(dd k2 )(2π)d δ(Σ2i ~ki ) Φ~< Φ~< (16.42)
k1 k2
with
Z Λ
I1 = (dd k) [r + ck 2 ]−1 . (16.43)
Λ/b
Finally, we carry out the integral over ~k2 and rescale to arrive at
Z Λ
hVi> |Φ< 2 Φ> 2 = uI1 b2 (dd k)|Φ~0k |2 (16.44)
This term has the identical shape as the usual mass term in the Landau-
Ginzburg (or Gauss) theory. There are 6 ways to choose the two wavevectors
of the fast modes out of the four wavevectors ~ki . To lowest order in V, we
find that these terms add to the renormalization of r (note that the mass
term appears with a factor 1/2),
r0 = b2 r + 12 b2 uI1 . (16.45)
Fig. 16.1: (a) four-leg vertex representing the term uφ4 in ~k-space. The
contact interaction translates into the condition of conserved momenta for
the fields on the legs. (b) Same with distinction between fast (dashed legs)
and slow (solid legs) modes. (c) one-loop diagram renormalizing r; the
paired up fields generate a Green’s function G> (~k ) and the momentum ~k
is integrated over. (d) Cluster diagram shifting the energy, of relevance in
the discussion of the free energy density, but irrelevant in the context of
renormalized couplings.
The Feynman diagrams are a mnemonic for properly organizing the per-
turbation theory. A perturbation typically involves a coupling constant g
with p fields, e.g., for our contact interaction we have g = u and p = 4. This
defines the vertex shown in Fig. 16.1(a), a point with four legs. The contact
nature in real space transforms into the conservation of momenta in ~k-space,
hence, the momenta ~ki attributed to the four legs add up to zero. When
separating modes into slow and fast ones, we draw solid and dashed lines for
these fields in order to distinct them, see Fig. 16.1(b). Taking averages over
the fast modes means joining up two dashed lines into a loop that gener-
ates correlators or propagaters, i.e., Green’s functions. Wick’s theorem tells
us, how to bind up lines/fields into Green’s functions: Taking the average
h. . . iG of a 2n-correlator/propagater over a Gaussian weight (or quadratic
Hamiltonian), we have to sum over all possible pairings,
X
hΦ1 Φ2 . . . Φk . . . Φ2n−1 Φ2n iG = Πn Gij , (16.46)
pairs ij
with Gij = hΦi Φj iG , what reduces the task of calculating a horrible 2n-
correlator to finding all allowed pairings and properly counting them (the
clue is that this works only for Gaussian/quadratic weights, this is why we
need perturbation theory at all). Correlators with odd number of fields
vanish.
In the calculation of the average hVi> , we have three different types
of diagrams: The one with four solid lines is the tree-level diagram and
generates the trivial scaling term b u in Eq. (16.39). The six diagrams with
two solid and two dashed lines generate the loop diagram in Fig. 16.1(c)
16.2. Φ4 THEORY 269
(actually, 1-loop) that give rise to the term 12b2 uI1 in Eq. (16.45)—we have
seen how the calculation of this type of diagram is done. The six versions of
this diagram arise from 4 possibilities to place the first dashed line on the
four momenta ~ki and the 3 possibilities to select the second moment, gives
12 possibilities, divided by two for dividing out trivial interchanges in what
is first and second, gives 6 (note the factor 1/2 in front of r, which takes the
factor 6 to 12/2). The last diagram is a cluster diagram, see Fig. 16.1(d).
It involves the integration over two Green’s functions, actually the square
of an integral over a Green’s function
Z Λ 2
>
hVi |Φ> 4 = u d ~
(d k) G(k ) (16.47)
Λ/b
and is a mere number producing a shift in the free energy (note that there
are 3 versions of this diagram).
[96]
[144]
[24]
[72]
[144]
(c)
[16]
[96]
Fig. 16.2: (a) Fully contracted (cluster) diagrams shifting the free energy.
(b) Two-loop diagrams with 2 external legs provide corrections to r in order
u2 (or 2 ) and are beyond the requested precision. (c) The six slow-fields
diagram generate a Φ6 term which is irrelevant at the Gaussian fixed point.
(d) Of the two diagrams with four slow fields the first is zero and the second
generates the renormalization of u to order . The numbers in square brack-
ets denote the symmetry factors (multiplicity of appearance of the term in
the perturbation series); note that a factor 1/2 has to be included in the
cumulant expansion of Eq. (16.36).
In going to the second order in V, we start with two 4-leg vertices; the
corrections will be of order u2 ∼ 2 and contribute only to the RG equation
270 CHAPTER 16. K-SPACE RG
for u. We divide the legs into two classes describing slow (solid lines) and
fast (dashed lines); in performing the average hV 2 i> we can pair dashed
lines (fast modes) on individual vertices or between vertices. The former
ones appear in the product hVi> 2 as well and exactly mutually cancel one
another. What remains are the connected diagrams and we write
2 2
hV 2 i> − hVi> = hV 2 i> >
c = hhVii , (16.48)
where the last bracket symbol hh. . . ii (double brackets) usually denotes a
cumulant, hence, the index c can be read as ‘connected’ or ‘cumulant’. So
the great news is that only diagrams with dashed lines contracted between
vertices matter. However, there are still many of them. A number of them is
not relevant in our discussion: the fully contracted (cluster) diagrams in Fig.
16.2(a) are mere numbers shifting the energy. The two-loop diagrams with
2 external legs, see Fig. 16.2(b), provide corrections to r in order u2 and are
beyond the requested precision. The six slow-fields diagram in Fig. 16.2(c)
generates a Φ6 term which is irrelevant at the Gaussian fixed point (see our
power counting argument above). Remain the four slow-fields diagram in
Fig. 16.2(d), of which the first is zero and the second (which is of a bow-tie
form), relevant one, is (cf. (16.38))
hV 2 i> 2 ~ ~
c |bt = u b I2 (k3 + k4 ) (16.49)
Z Λ
× Π4i (dd ki )(2π)d δ(Σ4i ~ki ) Φ~k1 Φ~k2 Φ~k3 Φ~k4
r0 = b2 r + 12 b2 uI1 , (16.52)
0 2
u = b u − 36 b u I2 .
16.2. Φ4 THEORY 271
3
r` − r = 2r + 2 u(Λ2 − r) `,
(16.55)
2π
h 9 i
u` − u = u − 2 u2 `.
2π
3
∂` r` = 2r` − 2 u` r` − 1 , (16.56)
2π
h 9 i
∂` u` = u` − 2 u` .
2π
The first terms on the right hand side of these equations, 2r` and u` , derive
from field rescaling, while the second terms account for the integration over
short-wave degrees of freedom. Furthermore, the term ∝ u` derives from
large k ∼ Λ; it is responsible for the shift in Tc , see also the discussion in
Sec. 12.3, and is the relevant term in the determination of the non-trivial
fixed point. The term ∝ u` r` acts only as a correction to this shift, but is
the relevant one in the determination of the critical exponent λτ , see below.
Besides the Gaussian fixed point at rG∗ = 0, u∗G = 0, these equations have
a second non-trivial (Wilson-Fisher) fixed point of order at
∗
rWF = −/6, (16.57)
∗ 2
uWF = (2π /9) .
272 CHAPTER 16. K-SPACE RG
The exponents associated with these fixed points follow from the linearized
equations
2 3/2π 2
B(r, u) = (16.59)
0
λτ = 2 and λu = (16.60)
τ u 2
~v = (1, 0) and ~v = (−3/4π , 1).
Hence, the flow directions at the two fixed points align and the Wilson-
Fisher fixed point lies on the ~v u direction of the Gaussian fixed point, see
Fig. 16.3.
Let us then have a look at the flow of the Landau-Ginzburg theory. This
depends crucially on the dimension. When d > 4 the Wilson-Fisher fixed
point resides in the unphysical region u < 0 and the only fixed point re-
maining is the Gaussian one. So the system’s behavior follows that dictated
by the Gaussian fixed point, where our results tell that u is an irrelevant
variable since λu = < 0 (although dangerously irrelevant, see below), while
r or τ is relevant.
When the dimension drops below 4, d < 4, = 4 − d > 0 and things
change quite appreciably. Now the Gaussian fixed point has become fully
repulsive and the critical behavior is governed by the Wilson-Fisher fixed
point which now resides in the physical region with u > 0. The ~v u di-
rection now is attractive with respect to the WF fixed-point and defines a
critical manifold—all Hamiltonians on this line flow to the WF fixed-point.
Deviating from this line, the parameters first flow towards the WF fixed-
point and then turn away to follow the (shifted) r or τ axis to the high- or
low-temperature phase. E.g., viewing the LG Hamiltonian as a continuum
16.2. Φ4 THEORY 273
cri
tica
u
d< 4
lm
flow
ani
under RG
fol
d
WF fixed point
8
T 0
G fixed point
r
unphysical region
u
d >4
G fixed point
unphysical r
region WF fixed point
description of the Ising model, the flow is towards the disordered paramag-
netic or ordered ferromagnetic phase. The transition then is governed by the
exponents ν = 1/λτ ≈ 1/(2 − /3) ≈ 1/2 + /12 > 1/2 and has an (anoma-
lous) correction θ ≈ −/6, see Eq. (13.20). Since elasticity is not normalized
at this order, the anomalous exponent η ≈ 0 vanishes—it will become finite
when going to order 2 , see below. Knowing the two exponents ν and η,
we can determine the remaining ones using the scaling laws: The exponent
γ = ν(2−η) ≈ 1+/6 follows from Fisher’s relation and α = 2−dν ≈ /6 fol-
lows from Josephson’s relation. The exponent β = (2 − γ − α)/2 ≈ 1/2 − /6
is derived with the help of Rushbrooke’s law and δ = 1 + γ/β ≈ 3 + is the
result of Widom’s law. In d = 3, = 1, we get the exponents α ≈ 0.166,
274 CHAPTER 16. K-SPACE RG
Above, we have carried out the renormalization on the level of the Hamil-
tonian—this is one of the great achievements of Ken Wilson. An important
element in this machinery, besides integration, is the process of rescaling.
But even before Wilson’s RG, renormalization was used in an alternative
manner without rescaling, the idea being, that a physical quantity is renor-
malized by the elimination of degrees of freedom through integration; the
elimination of these (short-scale) degrees of freedom then is compensated by
a change in the coupling or in the quantity under consideration. A famous
example of this type of ‘old-style’ RG is the renormalization of the charge
in QED, corresponding to the renormalization of u in the present context of
the φ4 -theory. In this RG, the integration is not over a small k-intervall, but
is usually carried out over all k-values from Λ until a large distance cutoff
R ≤ ξ(τ ) is reached; the flow of the quantity is then studied as a function
of the scale R, see also the discussion at the end of Sec. 16.1 above.
The RG equation (16.56) for the coupling u as a function of scale R
involves only u itself and we write it in the form
This is the Gell-Mann Low equation for the ‘charge’ u and β(u) is the Beta-
function providing the flow of u. For d = 4, we have β(u) = −(9/2π 2 )u2
(see Eq. (16.56) with = 0)
9 2
∂ln R u = − u , (16.71)
2π 2
and a simple integration provides us with the solution
u0
u(R) = . (16.72)
1 + (9u0 /2π 2 ) ln(R/a)
Hence, when integrating out degrees of freedom, the initial charge u0 de-
creases and renormalizes to zero as R/a → ∞—as a result, the theory
becomes ‘free at large distances’.
In a similar spirit, we can integrate out the degrees of freedom at scales
below R and obtain a renormalized mass r or a renormalized temperature
τ : in the absence of rescaling, there is no trivial term 2r in (16.56) and we
drop the u term in (16.56) by absorbing it in a shifted Tc , i.e., we replace
Σ1 → Σ1 (r) − Σ1 (r = 0) in the derivation of the RG equation. These
manipulations leave us with the equation
3
∂ln R r = − ur (16.73)
2π 2
276 CHAPTER 16. K-SPACE RG
that is solved by
r0
r(R) = . (16.74)
[1 + (9u0 /2π 2 ) ln(R/a)]1/3
1 [− ln |τ |]1/3
χ∼ ∼ , (16.76)
τren |τ |
1 [− ln |τ |]1/6
ξ∼ 1/2
∼ p . (16.77)
τren |τ |
We can use the above ‘old-style’ analysis also away from d = 4. We start
from Eq. (16.52) without rescaling cwcouplings and fields,
r0 = r + 12 uI1 , (16.78)
0 2
u = u − 36 u I2 ,
i.e., we only account for the integration over short-wave degrees of freedom.
For the integral I2 , we have (in d = 4)
Z Λ
S4 1
I2 ≈ dkk −1 ≈ ln(RΛ). (16.79)
(2π)4 1/R 8π 2
Redefining
a
u = ũ , (16.82)
R
this equation assumes the form of (16.56),
h 9 i
∂R ũ = ũ − 2 ũ (16.83)
2π
with a fixed point ũ∗ ∝ . We find several new results: First, the original
u ∝ /R flows to zero algebraically in R rather than logarithmically u ∝
1/ ln R, see Eq. (16.72). Second, we note that the u in (16.56) has been field-
rescaled with a factor b , what corresponds exactly to the step in (16.82),
hence our ũ above corresponds to the u in (16.56). So while it is the fixed
point plus exponent defined nearby that is of central interest in the Wilson
RG, it is the flow of the coupling that is in the focus of the present approach.
Repeating the above analysis for the coupling r, we note that the reduced
integral I1 (r) − I1 (0) that is relevant after absorption of the Tc shift has a
similar form as I2 ,
Z Λ
1 ru R
I1 (r) − I1 (0) = −ru (d4− k) 2 ≈ − , (16.84)
1/R rk + k 4 8π 2
that leads to the equation
3
∂R r = − rũ. (16.85)
2π 2
Its solution near the fixed point ũ∗ assumes the form
a /3
r(R) = r0 . (16.86)
R
Comparing massive and elastic terms in the Hamiltonian, we find that
r(ξ)/c ∼ 1/ξ 2 , from which we obtain the scaling for ξ in the form,
1/(2−/3)
ξ ∼ 1/τ0 , (16.87)
implying a new exponent ν = 1/(2 − /3) ≈ 1/2 + /12, in agreement with
the result of the Wilson approach. Similarly, we obtain for the susceptibility
χ(R) ∼ 1/τ (R) ∼ R/3 /τ0 at the scale R and stopping the flow at the scale
1+ν/3
ξ, we find that χ ∼ 1/τ0 , resulting in an exponent γ = 1 + /6.
As already explained in Sec. 12.3, we can find a perturbative series for the
correlator G(~k ). Expressing the Green’s function with the help of the self-
energy Σ,
G−1 (~k ) = G−1 ~ ~
0 (k ) − Σ(k ), (16.88)
278 CHAPTER 16. K-SPACE RG
The symmetry factor 96 is obtained by choosing the first, second, and third
leg in each vertex for connection with the second vertex, and dividing by 6
to avoid double-counting (we can order the sequence of 3 lines in 6 ways),
resulting in (4 · 3 · 2)2 /6 = 96. The first term Σ1 contributes a constant
and shifts Tc downwards due to large-k fluctuations, see Sec. 12.3. Here,
we are interested in the renormalization of the elastic energy k 2 and hence
determine the dispersive (i.e., ~k-dependent shift δΣ(k) = Σ2 (k) − Σ2 (0).
Since this term is already ∝ u2 ∼ 2 , we can work in d = 4. At Tc , we
have G0 (q) = 1/q 2 and transforming to real space, we obtain (using spatial
isotropy and ~r = (r, 0, 0))
d4 q ei~q·~r dp d3 k eipr
Z Z
G(r) = = (16.91)
(2π)4 q 2 (2π)4 p2 + k 2
π
Z −|k|r Z
3 e 1 1
= d k = dk k e−kr = 2 2 .
(2π)4 |k| 4π 2 4π r
For the dispersive shift, we obtain
96u2
Z Z
~
δΣ(k) = −96u2 d4 r G3 (r) [eik·~r − 1] ≈ d4 r (~k · ~r )2 G3 (r),
2
96u2
Z Z
3 3 dΩ ~
= dr r G (r) (k · ~r )2 (16.92)
2 2π 2
96u2 k 2
Z
= dr r5 G3 (r), (16.93)
8
where we have expanded exp(i~k · ~r ) − 1 ≈ i~k · ~r + (~k · ~r )2 /2 and performed
an angular average (~k · ~r )2 = k 2 r2 /4. The remaining integral is logarithmic;
the divergence at small x is cut on 1/Λ ∼ a, the one at large x is cut on 1/k
~
(due to the factor [eik·~r − 1]). We then find the result (with S 4 = 2π 2 the
surface area of a 4D sphere)
2π 2 Λ
δΣ(k) = 12u2 k 2 6
ln . (16.94)
64π k
The correlator thus picks up a logarithmic correction that we can transform
to a small power k η (since 1 − η ln k ≈ exp(−η ln k) = k −η ),
G−1 (k) = k 2 [1 − (3u2 /8π 4 ) ln(k/Λ)] = k 2−η Λη (16.95)
16.5. TWO-LOOP EXPANSION FOR G AND FINITE η 279
with
3u2 2
η= = , (16.96)
8π 4 54
where we have used the fixed point u∗WF = (2π 2 /9) in the last equation, see
Eq. (16.57). Note the numerical smallness of η. The anomalous scaling of
G(k) ∼ aη /k 2−η is in agreement with the discussion in Chap. 13, see Eq.
(13.17); as expected, the microscopic scale a enters the result through cut-
ting off an integral in the perturbation series. Translating the k-dependence
to a real-space dependence, we find that G(r) ∼ 1/rd−2+η , what produces
the field scaling zr = b(d−2+η)/2 , see Eq. (14.7).
Philosophically, the above renormalization of G(k) should be viewed as
a renormalization of the elasticity. Let us assume that τ is small but finite,
i.e., we are slightly away from the critical point. Then, the renormalization
produces a dispersive elasticity c(k) ∝ k −η on scales k > 1/ξ (‘short’ dis-
tances), what leads to a non-zero anomalous exponent η. On the other hand,
on scales k < 1/ξ (‘large’ distances), we should stop the flow at k ∼ 1/ξ and
the elasticity is renormalized to a new constant c → c[1 + (2 /54) ln(Λξ)].
Finally, the anomalous exponent η is defined at τ = 0 or ξ = ∞ and c is
dispersive on all scales.
The systematic calculation of the RG equations for the couplings r and
u (at ~k = 0) to second order (2-loop RG) is pretty involved. In particular,
care has to be taken not to double-count logarithmic terms that already
have been accounted for in the one-loop calculation.
280 CHAPTER 16. K-SPACE RG
Chapter 17
Berezinskii-Kosterlitz-
Thouless Transition
d2 R ~
Z
K
KXY ≈ [∇ϑsw ]2 (17.1)
2 a2
Z Z
− πK d2 R d2 R0 nv (R ~ 0 ) ln[|R
~ )nv (R ~ 0 |/a] + Ec
~ −R ~ ),
d2 R nv (R
z = exp(−Ec ). (17.2)
281
282 CHAPTER 17. BKT TRANSITION
Fig. 17.1: The interaction between two charges (constituting a large dipole)
is screened by smaller dipoles. The latter are polarized in the field of the
charges of the large dipole with angles θ preferentially in the forward direc-
tion; temperature fluctuations reduce this polarization.
with θ the angle between the dipole- and the electric-field vector and the
average has to be taken with the Boltzmann weight exp(−E qR cos θ/kB T ),
q2 2 q 2 R2 π
α(R) = R hcos2 θi = = KR2 , (17.4)
kB T 2kB T 2
where we have used that πK = q 2 /kB T (as follows from comparing the
prefactors of the log-interaction). Such pairs of extension R appear with a
density
~ d2 R = z 2 −2πK ln(R/a) d2 R
n(R) e . (17.5)
a2 a2
283
In (17.10), we integrate out pairs only in a thin shell a < R0 < ab, a typical
decimation step where we integrate over ‘fast’ degrees of freedom, here, small
vortex–anti-vortex pairs. Next, we perform the usual rescaling step in the
cutoff ab → a in order to recover the original relation (17.7)
∞ 3−2πKb
dR0 R0
Z
1 1
= + 4π 3 zb2 (17.11)
K∞ Kb a a a
but with renormalized couplings Kb and zb ,
zb = b2−πK z. (17.12)
~ d2 R = z 2 −2π R ln R d ln R0 K(R0 ) d2 R
n(R) e ln a . (17.15)
a2 a2
The perturbative RG equation is transformed into a self-consistent equation
for KL (we rewrite R → L = ln(R/a)),
Z L Z L0
1 1 3 2 0 0 00
= + 4π z dL exp 4L − 2π dL KL00 . (17.16)
KL K 0 0
z separatrix
critical
line of fixed point
initial points
separatrix
0
0 π/2 −1
K
T 0 T
8
∂` x = y 2 , (17.21)
∂` y = xy.
Taking their ratio, we find that xdx = ydy or d(x2 − y 2 ) = 0, which defines
a hyperbola
y 2 = x2 + C (17.22)
C = b2 (T − TBKT ). (17.23)
287
Next, we solve the flow for both cases above and below TBKT and derive
some physical results. For T > TBKT , we deal with a positive parameter C
and the flow takes off to infinity, after approaching the critical point. We
solve the differential equationR for x = x` , ∂` x = y 2 = x2 + C, by simple
2
R
integration, dx/(x + C) = dl and obtain (with a trajectory starting at
x0 , ` = 0 and ending at x, `)
√ √ √
C` = arctan(x/ C) − arctan(x0 / C). (17.24)
√ √
√ x − −C x0 + −C
2 −C` = ln √ √ . (17.25)
x + −C x0 − −C
The first solution above
√ TBKT provides us with an expression
√ for ξ(T ). On
approaching TBKT , C |x0 | and we have arctan(x0 / C) ≈ −π/2 since
x0 < 0. As ` increases, the trajectory approaches the critical point but then
deviates away from it on a scale ξ. We √ choose a symmetric trajectory
√ with
xξ = x(`ξ ) = −x0 such that arctan(xξ / C) ≈ π/2 and `ξ = π/ C. Using
the relation between distances R and the curve parameter `, ` = ln(R/a),
we find that the coherence length
√
T −TBKT
ξBKT (T ) ∼ a eπ/b (17.26)
J J(TBKT ) bp 2 bp
= + TBKT − T = + TBKT − T . (17.27)
kB T kB TBKT π π π
The second relation makes use of the universal jump in the stiffness at TBKT ,
J(TBKT ) 2 J∞
= ≡ , (17.28)
kB TBKT π kB TBKT
i.e., the stiffness collapses at TBKT from a finite value J∞ = (2/π)kB TBKT to
zero (Nelson-Kosterlitz universal jump, in the superfluid density).
A nice application is the establishment of true order and a finite Tc
in a layered system, e.g., a layered superconductor as is the case in the
288 CHAPTER 17. BKT TRANSITION
289
290 INDEX
Bosegase configurational, 74
BE-Kondensation, 125 integrating factor, 6
gauge invariance, 106 interaction
Gauss integral, 206 statistical, 85
matrix, 206 inversion, 7
Gauss model, 183, 204 involutory, 25
Gaussian approximation, 212 irrelevant variable, 239
Gaussian fixed point, 271 dangerous, 272
generating function, 97 Ising model, 136
Gibbs surface, 32 RG, 249
Gibbs’ paradoxon, 50 isochore
Gibbs’ phase rule, 31 critical, 173
Gibbs-Duhem, 28 isotherm
Ginzburg criterion, 183, 221, 226 van der Waals, 166
Ginzburg temperature, 220
jumps, 34
glass, 47
Goldstone bosons, 157 Kadanoff
Goldstone modes, 133 block spins, 235
grand potential, 60 kink, 25
Grassmann variables, 94
group Lambda transition, 187
semi-, 240 Landau
Fermi Liquid Theory, 79
harmonic approximation, 96 Landau diamagnetism, 112
Hartree approximation, 218 Landau expansion, 181
heat, 11 Landau levels, 110
latent, 33 Landau-Ginzburg model, 263
specific, 18, 54 Langevin dynamics, 48
at constant p, 12 Larmor radius, 110
at constant V , 12 latent heat, 33
Heisenberg model, 134 lattice gas model, 136, 177
Helium layered system, 287
nucleation, 42 Lee-Yang
Helium mixture, 198 circle theorem, 179
Hilbert space, 67 Legendre transformation, 25
Hohenberg-Mermin-Wagner, 133, 151, Lifshitz transition, 201
191 linear response, 4, 48
homogeneous function, 233 linear response theory, 211
hysteresis, 193 Liouville theorem, 69
Liouville’s theorem, 46
imaginary time, 90 liquid, 30
information, 25 overheated, 37
integrability, 6, 21, 26, 28 superheated, 168
integral liquid–crystal transition, 194
INDEX 293
Gibbs’, 50 pressure, 4
paramagnetismc, 107 primitve, 5
parameter probability
equilibrium, 52 a priori, 47, 69
particle number, 4 process
partition function adiabatic, 8
canonical, 57 adjabatic, 13
grand canonical, 60 cyclic, 11
path integral, 207 isobaric, 8
free particle, 90 isochoric, 8
particle in a potential, 92 isothermal, 8
Pauli reversible, 7
repulsion, 30 stochastic, 48
Pauli exclusion, 81 prozess
Pauli spin susceptibility, 115 irreversible, 8
Peierls, 152 pseudopotential, 86
percolation, 221
perturbation theory in u, 226 quantum gas
phase coexistence, 32 bosons, 81
phase diagram, 30, 32 fermions, 81
phase line, 35 partition sum, 81
phase rule quantum Hall effect, 110
Gibbs’, 31 quantum mechanics, 48
phase space, 45, 49 quantum simulation, 79
Γ-space, 49
µ-space, 74 random phase approximation, 218
average, 46 random phases, 68
probability, 46 Rayleigh-Jeans radiation law, 120
phase transition, 30 recursion
gas-liquid, 66 RG, 240
order, 179 relation
phenomenological theory, 3 Maxwell, 26, 28
photons, 117 thermal–caloric, 17
Planck relevant variable, 239
radiation law, 121 renormalization
Planck-Nernst, 17 in ~k-space, 260
polarizability, 282 renormalization group, 183
potential rescaling, 260
chemical, 20, 65 reservoir, 48, 67
grand, 60 response, 48
potentials linear, 4
convex, 168 resummation methods, 226
thermodynamic, 21 RG equations
power counting, 265, 266 differential, 245
INDEX 295
Wenzel-Kramers-Brillouin, 76
Wick rotation, 90
Widom scaling, 231
Wien
displacement law, 120
Wilson-Fisher fixed point, 271
winding number, 159
WKB, 76
work, 11
XY model, 136