Path Integrals in Quantum Mechanics
Path Integrals in Quantum Mechanics
2010-2011
2
Chapter 1
1.1.1 Functionals
What is a functional? Let us first consider a real function f : [a, b] → R, mapping all elements x ∈ [a, b] in the
domain [a, b] to real numbers f (x). A functional is similar to a function in that it maps all elements in a certain
domain to real numbers, however, the nature of its domain is very different. Instead of acting on all points of an
interval or some other subset of the real numbers, the domain of functionals consists of (suitably chosen) classes
of functions itself. In other words, given some class {f } of functions, a functional F : {f } → R assigns to each
such f a real number F [f ].
Differentiation and integration are the two main operations on which the calculus of (real or complex) functions
is based. We would like to understand how to perform analogous operations on functionals, that is, we would like
to make sense of the expressions
δF [f ]
Z
, Df F [f ] (1.1)
δf
for functional differentiation and integration. Here, we denote a functional derivative by δ/δf , replacing the
ordinary d/dx or ∂/∂x for functions, and the measure of a functional integral by Df , replacing dx in the case of
integrals over functions. A common application of ordinary differentiation is to find the extrema of a function f
by solving the equation f ′ (x) = 0. Of course, one can also ask about the extrema of a functional F , that is the
functions f which minimise (or maximise) the value F [f ]. One important application of functional differentiation
which we will discuss is how these extrema can be obtained as solutions to the equation
δF [f ]
=0. (1.2)
δf
To make our discussion more concrete we will now introduce two (well-known) examples of functionals.
What is the minimum of this functional? Of course, we know that the correct answer is ”a line”, but how do we
prove this formally? Let us approach this problem in a somewhat pedestrian way first before re-formulating it in
3
4 CHAPTER 1. PATH INTEGRALS IN QUANTUM MECHANICS
terms of functional derivatives along the lines of Eq. (1.2). Consider a small but otherwises abritrary perturbation
ǫ of a path f which vanishes at the endpoints, that is, which satisfies ǫ(x1 ) = ǫ(x2 ) = 0. We can then define the
function Z x2 p
l(λ) ≡ F [f + λǫ] = dx 1 + (f ′ (x) + λǫ′ (x))2 , (1.4)
x1
where λ is a real parameter. A necessary condition for f to be a minimum of the functional F is then
dl
(0) = 0 (1.5)
dλ
Hence, this simple trick has reduced our functional minimisation problem to an ordinary one for the function l.
The derivative of l at λ = 0 can be easily worked out by differentiating ”under the integral” in Eq. (1.4). One finds
Z x2
dl f ′ (x)ǫ′ (x)
(0) = dx p . (1.6)
dλ x1 1 + f ′ (x)2
The derivative on ǫ in this expression can be removed by partial integration, keeping in mind that the boundary
term vanishes due to ǫ being zero at x1 and x2 . This results in
Z x2 " #
dl d f ′ (x)
(0) = − dx ǫ(x) p . (1.7)
dλ x1 dx 1 + f ′ (x)2
From Eq. (1.5) this integral needs to vanish and given that ǫ is an arbitrary function this is only the case if the
integrand is zero pointwise in x. This leads to the differential equation
" #
d f ′ (x)
p =0, (1.8)
dx 1 + f ′ (x)2
for the function f . The desired extrema of the length functional F must be among the solutions to this differential
equation which are given by the lines
f ′ (x) = const . (1.9)
Note that the differential equation (1.8) is second order and consequently its general solution has two integration
constant. They are precisely ”used up” by implementing the boundary conditions f (x1 ) = y1 and f (x2 ) = y2 , so
that we remain with a unique solution, the line between our two chosen points.
It depends on classical paths q(t) between times t1 and t2 satisfying the boundary conditions q(t1 ) = q1 and
q(t2 ) = q2 . Hamilton’s principle states that the classical path of the system minimises the action functional S. To
work out the implications of Hamilton’s principle we follow the same steps as in the previous example. We define
the function Z t2
l(λ) ≡ S[q + λǫ] = dt L(q + λǫ, q̇ + λǫ̇) . (1.11)
t1
where we have introduced small but arbitrary variations ǫ = (ǫ1 , . . . , ǫn ) of the coordinates, satisfying ǫ(t1 ) =
ǫ(t2 ) = 0. As before, we work out the derivative of l at λ = 0
Z t2
dl ∂L ∂L
(0) = dt (q, q̇) ǫi + (q, q̇) ǫ̇i (1.12)
dλ t1 ∂qi ∂ q̇i
Z t2
∂L d ∂L
= dt ǫi (q, q̇) − (q, q̇) , (1.13)
t1 ∂qi dt ∂ q̇i
1.1. MATHEMATICAL BACKGROUND 5
where the last step follows from partial integration (and the boundary term vanishes due to ǫ(t1 ) = ǫ(t2 ) = 0 as
previously). The ǫi vary independently and, hence, for the above integral to vanish the bracketed expression in
(1.13) must be zero for each index value i. This gives rise to the Euler-Lagrange equations
d ∂L ∂L
(q, q̇) − (q, q̇) = 0 . (1.14)
dt ∂ q̇i ∂qi
The solutions to these equations represent the classical paths of the system which minimise the action functional.
They contain 2n integration constants which are fixed by the boundary conditions q(t1 ) = q1 and q(t2 ) = q2 .
In the bracket we recognise the desired left-hand-side of the differential equation (1.8). Our last step consisted
of removing the integral due to the presence of an arbitrary variation ǫ in the integrand. Here, we can formally
reproduce this step by demanding the relation
δf (x̃)
= δ(x̃ − x) , (1.17)
δf (x)
which can be viewed as a continuous version of the equation ∂qi /∂qj = δij for a set of n coordinates qi . Using
the relation (1.17) in Eq. (1.16) we finally obtain the desired result
" #
δF [f ] d f ′ (x)
=− p . (1.18)
δf (x) dx 1 + f ′ (x)2
To summarise, functional derivation δ/δf (x) can be understood as a linear operation which satisfies the product
and chain rules of ordinary differentiation, commutes with ordinary integrals and derivatives and is subject to the
”normalisation” (1.17).
Armed with this ”definition”, let us quickly re-consider our second example, Hamilton’s principle of classical
mechanics. We find
Z t2
δS[q] δ
= dt̃ L(q(t̃), q̇(t̃)) (1.19)
δqi (t) δqi (t) t1
Z t2
∂L δqj (t̃) ∂L δ q̇j (t̃)
= dt̃ (q, q̇) + (q, q̇) (1.20)
t1 ∂qj δqi (t) ∂ q̇j δqi (t)
Z t2
∂L d ∂L δqj (t̃)
= dt̃ (q, q̇) − (q, q̇) (1.21)
t1 ∂qj d t̃ ∂ q̇j δqi (t)
Z t2
∂L d ∂L ∂L d ∂L
= dt̃ (q, q̇) − (q, q̇) δij δ(t̃ − t) = (q, q̇) − (q, q̇) (1.22)
t1 ∂qj d t̃ ∂ q̇j ∂qi dt ∂ q̇i
6 CHAPTER 1. PATH INTEGRALS IN QUANTUM MECHANICS
where, in the second last step, we have used the obvious generalisation
δqj (t̃)
= δij δ(t̃ − t) (1.23)
δqi (t)
of Eq. (1.17).
We also need a discrete version F̃ of the functional such that F̃ (y) ∼ F [f ] for y ∼ f . Then, we can formally
introduce functional integrals by setting
Z Z
Df F [f ] ∼ lim dn y F̃ (y) . (1.24)
∆x→0
In practise, this relation means we can often represent functional integrals by multiple ordinary integrals or, at
least, deduce properties of functional integrals from those of multiple integrals.
In general, functional integrals are prohibitively difficult to evaluate and often a lattice calculation (that is, the
numerical evaluation of the RHS of Eq. (1.24) on a computer) is the only option. However, there is one class of in-
tegrals, namely Gaussian integrals and integrals closely related to Gaussian ones, which is accessible to analytical
computations. It turns out that such Gaussian functional integrals are of central importance in physics and we will
focus on them in the remainder of this section. In keeping with the spirit of Eq. (1.24), the discussion will center
on ordinary multiple Gaussian integrals from which we deduce properties of Gaussian functional integrals.
As a reminder, we start with a simple one-dimensional Gaussian integral over a single variable y. It is given by
∞
r
1 2π
Z
I(a) ≡ dy exp(− ay 2 ) = , (1.25)
−∞ 2 a
I(a)2 as a two-dimensional
where a is real and positive. The standard proof of this relation involves writing p
integral over y1 and y2 and then introducing two-dimensional polar coordinates r = y12 + y22 and ϕ. Explicitly,
∞ Z ∞ Z ∞ Z ∞
1 1 1
Z
I(a)2 = dy1 exp(− ay12 ) dy2 exp(− ay22 ) = dy1 dy2 exp(− a(y12 + y22 ))(1.26)
−∞ 2 −∞ 2 −∞ −∞ 2
Z 2π Z ∞
1 2π
= dϕ dr r exp(− ar2 ) = . (1.27)
0 0 2 a
1.1. MATHEMATICAL BACKGROUND 7
over variables y = (y1 , . . . , yn ), where A is a symmetric, positive definite matrix. This integral can be easily
reduced to a product of one-dimensional Gaussian integrals by diagonalising the matrix A. Consider an orthogonal
rotation O such that A = ODOT with a diagonal matrix D = diag(a1 , . . . , an ). The Eigenvalues ai are strictly
positive since we have assumed that A is positive definite. Introducing new coordinates ỹ = Oy we can write
n
X
yT Ay = ỹT Dỹ = ai ỹi2 , (1.29)
i=1
where the property OT O = 1 of orthogonal matrices has been used. Note further that the Jacobian of the coordi-
nate change y → ỹ is one, since |det(O)| = 1. Hence, using Eqs. (1.25) and (1.29) we find for the integral (1.28)
n Z
Y 1
W0 (A) = dỹi exp(− ai ỹi2 ) = (2π)n/2 (a1 a2 . . . an )−1/2 = (2π)n/2 (detA)−1/2 . (1.30)
i=1
2
To summarise, we have found for the multidimensional Gaussian integral (1.28) that
Here J = (J1 , . . . , Jn ) is an n-dimensional vector which can be viewed as an external source for y. With the
inverse matrix ∆ = A−1 and a change of variables y → ỹ, where
y = ∆J + ỹ (1.33)
The exponential in Gaussian integrals can be considered as a probability distribution and, hence, we can intro-
duce expectation values for functions O(y) by
1
Z
hOi0 ≡ N dn y O(y) exp − yT Ay . (1.36)
2
The sub-script 0 refers to the fact that the expectation value is computed with respect to a Gaussian distribution.
As we will see later, such expectation values play a central role when it comes to extracting physical information
from a quantum field theory. Of particular importance are the moments of the Gaussian probability distribution (or
8 CHAPTER 1. PATH INTEGRALS IN QUANTUM MECHANICS
l-point functions) which are given by the expectation values of monomials, hyi1 yi2 . . . yil i0 . From Eqs. (1.32) and
(1.35) they can be written as
∂ ∂ ∂ ∂ 1 T
hyi1 yi2 . . . yil i0 = N ... W0 (A, J) = ... exp J ∆J . (1.38)
∂Ji1 ∂Jil J=0 ∂Ji1 ∂Jil 2 J=0
The first equality above suggests an interpretation of the Gaussian integral W0 (A, J) with source J as a generating
function for the l-point functions. The second equality provides us with a convenient way of calculating l-point
functions explicitly. It is easy to see, for example, that the two-point and four-point functions are given by
hyi yj i0 = ∆ij , hyi yj yk yl i0 = ∆ij ∆kl + ∆ik ∆jl + ∆il ∆jk . (1.39)
Every differentiation with respect to Ji in Eq. (1.38) leads to a factor of ∆ij Jj in front of the exponential. Unless
Jj is removed by another differentiation the contribution vanishes after setting J to zero. In particular, this means
that all odd l-point functions vanish. For the even l-point functions one deduces that the ∆ij terms which appear
are in one-to-one correspondence with all ways of pairing up the indices i1 , . . . , il . Hence, for even l we have
X
hyi1 . . . yil i0 = ∆ip1 ip2 . . . ∆ipl−1 ipl (1.40)
pairings p
where the sum runs over all pairings p = {(p1 , p2 ), . . . , (pl−1 , pl )} of the numbers {1, . . . , l}. This statement is
also known as Wick’s theorem. For l = 2 and l = 4 this formula reproduces the results (1.39) for the two- and
four-point functions.
We would now like to use some of the above formulae for multidimensional Gaussian integrals to infer analo-
gous formulae for Gaussian functional integrals, following the corrrespondence (1.24). We start with the analogue
of W0 (A, J) in Eq. (1.32), the generating functional
1
Z Z Z Z
W0 [A, J] = Df exp − dx dx̃ f (x̃)A(x̃, x)f (x) + dx J(x)f (x) . (1.41)
2
Note, that the summations in the expondent of Eq. (1.32) have been replaced by integrations. The ”kernel” A(x̃, x)
is the generalisation of the matrix A and, in our applications, will be usually a differential operator. It is straight-
forward to write down the analogue of the result (1.35) for the generating functional
Z
1 1
Z
W0 [A, J] = exp − tr ln A exp dx̃ dxJ(x̃)∆(x̃, x)J(x) . (1.42)
2 2
where ∆ = A−1 is the inverse of the operator A and we have used the well-known matrix identity detA =
exp tr ln A to re-write the determinant of A in terms of its trace. (We have absorbed the powers of 2π in front of
Eq. (1.35) into the definition of the path integral measure Df .) It is not clear, at this stage, how to compute the
inverse and the trace of the operator A, but we will see that, in many cases of interest, this can be accomplished by
looking at the Fourier transform of A. For now, we proceed formally and define the l-point functions
1
Z Z Z
hf (x1 ) . . . f (xl )i0 = N Df f (x1 ) . . . f (xl ) exp − dx dx̃ f (x̃)A(x̃, x)f (x) , (1.43)
2
where the normalisation N is given by N = exp( 12 tr ln A). In analogy with Eq. (1.38) we have
Z
δ δ 1
Z
hf (x1 ) . . . f (xl )i0 = ... exp dx̃ dx J(x̃)∆(x̃, x)J(x) . (1.44)
δJ(x1 ) δJ(xl ) 2 J=0
This expression can be used to compute the l-point functions in terms of ∆, where Eq. (1.17) should be taken into
account. Since functional derivatives have the same general properties than ordinary derivatives Wick’s theorem
applies and the l-point functions take the form
X
hf (x1 ) . . . f (xl )i0 = ∆(xp1 , xp2 ) . . . ∆(xpl−1 , xpl ) (1.45)
pairings p
1.1. MATHEMATICAL BACKGROUND 9
with a real function s = s(y) and a real positive parameter ξ. (In our physical examples the role of ξ will be played
by ~.) We would like to study the value of this integral for small ξ. For such values, the main contributions to the
integral come from regions in y close to the minima of s, that is points y (c) satisfying s′ (y (c) ) = 0. Around such a
minimum we can expand
1
s(y) = s(y (c) ) + ξaǫ2 + . . . , (1.47)
2
√
where ǫ = (y − y (c) )/ ξ, a = s′′ (y (c) ) and the dots stand for terms of order higher than ξ. Then, we can
approximate
1
p Z
I(ξ) ≃ ξ exp(−s(y (c) )/ξ) dǫ exp − aǫ2 . (1.48)
2
The remaining integral is Gaussian and has already been computed in Eq. (1.25). Inserting this result we finally
find r
2πξ
I(ξ) = exp(−s(y (c) )/ξ) . (1.49)
a
To leading order the integral is simply given by the exponential evaluated at the minimum y (c) and the square root
pre-factor can be interpreted as the first order correction.
Let us generalise this to the multi-dimensional integrals
Z
W (ξ) = dn y exp(−S(y)/ξ) , (1.50)
∂S
where y = (y1 , . . . , yn ). Around a solution y(c) of ∂yi (y) = 0 we expand
1
S(y) = S(y(c) ) + ξǫT Aǫ + . . . , (1.51)
2
√ ∂2S
where ǫ = (y − y(c) )/ ξ and Aij = ∂yi ∂yj (y
(c)
). Using this expansion and the result (1.31) for the multi-
dimensional Gaussian integral we find
−1/2
W (ξ) = (2πξ)n/2 (detA) exp(−S(y(c) )/ξ) . (1.52)
We would now like to apply the steepest descent method to the generating functional for l-point functions. In
analogy with our procedure for Gaussian integrals, we can introduce this generating function by
1
Z
n T
W (J, ξ) = d y exp − (S(y) − J y) (1.53)
ξ
so that the l-point functions are given by
∂ ∂
Z
hyi1 . . . yil i ≡ N dn y yi1 . . . yil exp(−S(y)/ξ) = ξ l N ... W (J, ξ) . (1.54)
∂Ji1 ∂Jil J=0
Minima y(c) of the exponent in Eq. (1.53) are obtained from the equations
∂S
Ji = . (1.55)
∂yi
Applying the result (1.52) to the generating function (1.53) one finds
−1/2 1 (c) T (c)
W (J, ξ) ∼ (detA) exp − (S(y ) − J y ) (1.56)
ξ
10 CHAPTER 1. PATH INTEGRALS IN QUANTUM MECHANICS
From Eq. (1.55), we should think of y(c) as a function of the source J which has to be inserted into the above result
for W (J, ξ). The l-point functions in the steepest descent approximation can then be computed by inserting the
so-obtained W (J, ξ) into Eq. (1.54).
Finally, we apply the steepest descent method to functional integrals of the form
Z
Wξ = Df exp(−S[f ]/ξ) , (1.57)
where S[f ] is a functional. In our examples S will be the action of a physical system. The steepest descent
approximation now amounts to an expansion around the ”classical solution” f (c) characterised by
δS (c)
[f ] = 0 . (1.58)
δf
With
δ2 S
A(x̃, x) = [f (c) ] (1.59)
δf (x̃)δf (x)
we find in analogy with the multi-dimensional steepest descent approximation (1.52)
1
Wξ ∼ exp − tr ln A exp −S[f (c)]/ξ . (1.60)
2
Note that the leading contribution to the functional integral in this approximation comes from the classical path
f (c) . We can now introduce the generating functional
1
Z Z
Wξ [J] = Df exp − S[f ] − dx J(x)f (x) , (1.61)
ξ
which generates the l-point functions
δ δ
Z
l
hf (xi1 ) . . . f (xil )i ≡ N Df f (xi1 ) . . . f (xil ) exp(−S[f ]/ξ) = ξ N ... Wξ [J] .
δJ(xi1 ) δJ(xi1 ) J=0
(1.62)
The steepest descent approximation can then be applied around the solution f (c) = f (c) (J) of
δS
J(x) = (1.63)
δf (x)
and leads to the result
1 1
Z
(c) (c)
Wξ [J] ∼ exp − tr ln A exp − S[f ] − dx J(x)f (x) (1.64)
2 ξ
(We can set ξ to one for the purpose of this section.) Here, we take the ”interaction term” V (y) to be a polynomial
in y and λ is the ”coupling constant”. Expanding the part of the exponential which involves the ”perturbation”
λV , we can write
∞
X (−λ)k
Wλ (A) = W0 (A) hV (y)k i0 , (1.66)
k!
k=0
where the expectation value h. . .i0 is defined with respect to the Gaussian measure, as in Eq. (1.36). Given that
V (y) is polynomial each of these expectation values can be calculated using Wick’s theorem (1.40). Of course,
the resulting expressions will be progressively more complicated for increasing powers in λ but, provided λ is
1.1. MATHEMATICAL BACKGROUND 11
sufficiently small, we can hope that the integral is well-approximated by cutting the infinite series in (1.66) off at
some finite value of k. Let us see how this works in practise by focusing on the example of a quartic interaction
term
n
1 X 4
V (y) = y . (1.67)
4! i=1 i
Neglecting terms of order λ3 and higher we have
Wλ (A) λX 4 λ2 X 4 4
= 1− hyi i0 + hy y i0 + O(λ3 ) (1.68)
W0 (A) 4! i 2(4!)2 i,j i j
λX 2 λ2 X 2 X 2 λ2 X
= 1− ∆ii + ∆ii ∆jj + (3∆ii ∆jj ∆2ij + ∆4ij ) + O(λ3 ) (1.69)
8 i 128 i j
48 i,j
where Wick’s theorem (1.40) has been used in the second step. There is a very useful representation of expressions
of the form ∆i1 i2 . . . ∆il−1 il in terms of diagrams. First, recall that the indices i, j, . . . label points xi , xj , . . . of
our discrete space-time. This suggests we should represent each such index by a point (or dot). Each ∆ij is then
represented by a line connecting dot i and dot j. One may think of the quantity ∆ij as ”moving the system” from
xi to xj and it is, therefore also called propagator. Given these rules, the diagrams for the terms at order λ and
order λ2 in Eq. (1.69) are shown in Figs. 1.1 and 1.2, respectively. Obviously, these diagrams are constructed by
Figure 1.1: Feynman diagram contributing to Wλ (A) at first order in perturbation theory for a quartic interaction.
i
i j i j
j
Figure 1.2: Feynman diagrams contributing to Wλ (A) at second order in perturbation theory for a quartic interac-
tion.
joining up four-leg vertices, one for each power of the coupling λ. The reason we are dealing with four-leg vertices
is directly related to the quartic nature of the interaction term, as an inspection of Wick’s theorem (1.40) shows.
The diagrams associated to the order λk terms in the perturbation expansion can then be obtained by connecting
the legs of k four-vertices in all possible ways.
We are also interested in computing the l-point functions
1
Z
(l) (l)
hyi1 . . . yil iλ = Wλ (A)−1 Gi1 ...il , Gi1 ...il ≡ dn y yi1 . . . yil exp − yT Ay − λV (y) (1.70)
2
(l)
in perturbation theory. Here, we have defined the Green functions Gi1 ...il . Since we have already discussed Wλ (A)
the Green functions are what we still need to compute, in order to fully determine the l-point functions. One way
to proceed is to simply expand the Green function in powers of λ and write the individual terms as Gaussian
expectation values (1.36). This leads to
∞
(l)
X (−λ)k
Gi1 ...il = W0 (A) hyi1 . . . yil V (y)k i0 . (1.71)
k!
k=0
Each Gaussian expectation value in this expansion can be calculated using Wick’s theorem. Alternatively, in
analogy with the Gaussian case, we may also introduce the generating function
1 T
Z
n T
Wλ (A, J) = d y exp − y Ay − λV (y) + J y (1.72)
2
12 CHAPTER 1. PATH INTEGRALS IN QUANTUM MECHANICS
around the Gaussian generating function W0 . The partial differential operator V (∂/∂J) is obtained by formally
replacing yi with ∂/∂Ji in the function V = V (y). Combining the last two equations we find
∞
" k #
(l)
X (−λ)k ∂ ∂ ∂
Gi1 ...il = ... V W0 (A, J) (1.75)
k! ∂Ji1 ∂Jil ∂J
k=0 J=0
∞
" k #
(−λ)k
X ∂ ∂ ∂ 1 T
= W0 (A) ... V exp J ∆J , (1.76)
k! ∂Ji1 ∂Jil ∂J 2
k=0 J=0
where we have used Eq. (1.35) in the second step. From our result (1.38) for Gaussian l-point functions this
expression for the Green function is of course equivalent to the earlier one (1.71). Either way, we end up having
to compute each term in the perturbation series applying Wick’s theorem. Let us see how this works explicitly by
computing the two-point function for the quartic perturbation (1.67) up to first order in λ. We find for the Green
function
(2)
Gij /W0 (A) = hyi yj i0 − λhyi yj V (ym )i0 + O(λ2 ) (1.77)
λ X
4 λX
= ∆ij − ∆ij < ym >0 − ∆mi ∆mm ∆mj + O(λ2 ) . (1.78)
24 m
2 m
The three terms in this expression can be represented by the Feynman diagrams shown in Fig. 1.3. At zeroth order
i j
i j m
m
i j
Figure 1.3: Feynman diagram contributing to two-leg Green function at zeroth and first order in perturbation theory
for a quartic interaction.
in λ the two-point function is simply given by the propagator and at first order we have two diagrams. The first
of these, however, is disconnected and consists of a propagator and a ”vacuum bubble” which we have earlier (see
Fig. 1.1) recognised as one of the contributions to Wλ (A). If we compute the two-point function from Eq. (1.70)
by inserting the Green function (1.78) and our earlier result (1.69) for Wλ (A) (dropping terms of order λ2 ) we
find that this disconnected term cancels, so that
λX
hyi yj iλ = ∆ij − ∆mi ∆mm ∆mj + O(λ2 ) . (1.79)
2 m
Hence, the two-point function up to first order in λ consists of the two connected diagrams in Fig. 1.3 only.
The cancellation of disconnected diagrams in the two-point function is an attractive feature since it reduces the
number of diagrams and removes ”trivial” contributations. However, in general, l-point functions still do contain
disconnected diagrams. For example, the four-point function at order λ2 has a disconnected contribution which
consists of a propagator plus the third diagram in Fig 1.3. It is, therefore, useful to introduce a new generating
function Zλ (A, J) whose associated Green functions
(l) ∂ ∂
Gi1 ...il = ... Zλ (A, J) , (1.80)
∂Ji1 ∂Jil J=0
1.2. QUANTUM MECHANICS 13
the so-called connected Green functions, only represent connected diagrams. It turns out that the correct definition
of this new generating function is
Zλ (A, J) = ln Wλ (A, J) . (1.81)
One indication that this is the right choice comes from the result (1.69) for the vacuum contribution Wλ (A) which
can also be written as
2 X
Wλ (A) λ X λ
= exp − ∆2ii + (3∆ii ∆jj ∆2ij + ∆4ij ) + O(λ3 ) . (1.82)
W0 (A) 8 i 48 i,j
Note that the third term in Eq. (1.69) which corresponds to the disconnected diagram in Fig. 1.2 has indeed been
”obsorbed” by the exponential. It turns out that this is a general feature. As a result, connected Green functions
can be evaluated in terms of Feynman diagrams in the same way ordinary Green function can, except that only
connected Feynman diagrams are taken into account.
As before, the above results for multiple integrals can be formally re-written for functional integrals. We
(l) (l)
can simply carry out the replacements
P R yip → f (xp ), Gi1 ...il → G (x1 , . . . , xl ), ∆ip iq → ∆(xp , xq ),
∂/∂Jip → δ/δJ(xp ) and ip → dxp . This leads to the formalism for perturbative quantum (field) the-
ory based on functional integrals. Perturbation theory is arguably the most important method to extract physical
information from quantum field theories and we will come back to it more explicitly later in the course.
∂ψ(t)
i~ = Hψ(t) .
∂t
14 CHAPTER 1. PATH INTEGRALS IN QUANTUM MECHANICS
ψ(t) = exp(−iHt/~)ψ(0)
Here, e−iHt/~ is the time-evolution operator. It is a function (the exponential) of another operator (the Hamilto-
nian), and we need to remember that functions of operators can be defined via the Taylor series for the function
involved.
Next, we’d like to think of the particle’s path as passing through points xn at the successive time-steps nt/N . We
can do this by introducing position eigenstates |xi and inserting the identity in the form
Z ∞
1= dx|xihx| (1.84)
−∞
between each factor on the right of Eq. (1.83). In this way we get
Z Z
hx0 |e−iHt/~ |xN i = dx1 . . . dxN −1 hx0 |e−iHt/N ~ |x1 ihx1 |e−iHt/N ~ |x2 i . . . hxN −1 |e−iHt/N ~ |xN i .
The way to do this is by separating the contributions from T and V , using the fact that, for N → ∞
N N
e−iHt/N ~ = e−iT t/N ~ e−iV t/N ~ .
Now we use two versions of the resolution of the identity: the one in Eq. (1.84) based on position eigenstates,
which diagonalise V , and a second one based on momentum eigenstates |pi, which diagonalise T . We have
Z
−iT t/N ~ −iV t/N ~
hxn |e e |xn+1 i = dpn hxn |e−iT t/N ~ |pn ihpn |e−iV t/N ~ |xn+1 i .
The action of the operator V on the state |xn+1 i is simple: V |xn+1 i = V (xn+1 )|xn+1 i, where V (x) denotes the
c-number function. Similarly for p on |pn i: p|pn i = pn |pn i, where p is the operator and pn is a number. Then
Z Z
2
dpn hxn |e−iT t/N ~ |pn ihpn |e−iV t/N ~ |xn+1 i = dpn e−ipn t/2mN ~ hxn |pn ihpn |xn+1 ie−iV (xn+1 )t/N ~ .
Next, the matrix element between states from our two bases is hx|pi = (2π~)−1/2 eipx/~ (check that you can justify
this), and so
1/2
mN
Z Z
2 2 2
dpn e−ipn t/2mN ~ hxn |pn ihpn |xn+1 i = (2π~)−1 dpn e−ipn t/2mN ~−ipn (xn+1 −xn )/~ = eimN (xn+1−xn ) /2~t .
2πi~t
Putting everything together, and writing the interval between time-slices as ǫ = t/N , we get
N −1
" 2 #!
m N/2 Z iǫ X m xn+1 − xn
Z
−iHt/~
hx0 |e |xN i = lim dx1 . . . dxN −1 exp − V (xn+1 ) .
N →∞ 2πiǫ~ ~ n=0 2 ǫ
(1.85)
1.2. QUANTUM MECHANICS 15
Finally, we interpret the sum in the argument of the exponential as a time integral, writing
N −1
" 2 #
iǫ X m xn+1 − xn i t i
Z
− V (xn+1 ) = L dt′ ≡ S[x(t′ )]
~ n=0 2 ǫ ~ 0 ~
The notation here is as follows: D[x(t′ )] indicates a functional integral over paths x(t′ ), in this case starting from
R
xi at t′ = 0 and ending at xf at t′ = t, and N is a normalisation factor, which is independent of what the start and
end points are.
The significance of Eq (1.86) is that it expresses the quantum amplitude — for a particle to propagate from a
starting point xi to a finishing point xf in time t — in terms of an integral over possible paths, with the contribution
from each path weighted by a phase factor eiS/~ that is expressed in terms of a classical quantity, the action S for
the path. In this sense, we have succeeded in extending the treatment of interference in the context of Young’s slits
to a general quantum system.
Correspondence principle
One of the simplest and most attractive features of the path integral formulation of quantum mechanics is that it
provides a framework in which we can see classical mechanics arising as a limiting case. We expect classical
mechanics to emerge in the limit ~ → 0 (or more precisely, when ~ is small compared to the characteristic scale in
S). Within the path integral, this limit is just the one in which we can apply the steepest descents approximation,
since then the argument of the exponential in Eq. (1.86) is multiplied by the large quantity ~−1 . Using that
approximation, the dominant contribution to the path integral will come from the vicinity of the path for which the
action S is stationary. But we know from Hamilton’s principle that this path is precisely the classical one, and so
we have reached the desired conclusion.
with " #
β~ 2
m dx
Z
S[x(τ )] = dτ + V (x) . (1.88)
0 2 dτ
The crucial point to note about Eq. (1.88) is that because of the substitution τ = it′ , there is a change in the relative
sign of the terms in the action coming from kinetic and potential energies. With this sign change, S[x(τ )] is often
referred to as the Euclidean action, by analogy with the change from a Lorentz metric to a Euclidean one, produced
by a similar factor of i. Our result, Eq. (1.87), gives matrix elements of the Boltzmann factor for a quantum system
as a functional integral over classical paths, this time with paths weighted by a real factor, e−S[x(τ )]/~.
Quantities we want to calculate in the context of statistical mechanics are the partition function Z and thermal
averages, for example the average hf (x)i of a function f (x) of the coordinate x of the particle. These are related
in a simple way to the matrix elements of the Boltzmann factor that we have discussed.
P −βEConsider first the partition
function, and recall its definition for a system with energy levels Ei : Z = i e i
. Clearly, this is simply
16 CHAPTER 1. PATH INTEGRALS IN QUANTUM MECHANICS
Z = Tre−βH , written in the eigenbasis of H. Since the trace of a matrix is invariant under change of basis, we can
use this second form, which leads us to
Z Z
Z = dx hx|e−βH |xi ≡ N D[x(τ )]e−S[x(τ )]/~ ,
where the path integral is now over periodic paths, satisfying x(β~) = x(0). In a similar fashion, we can evaluate
thermal averages as Z
hf (x)i = Z −1 D[x(τ )]f (x(0))e−S[x(τ )]/~ .
[Prove this by completing the square in the argument of the exponential.] We use this result to integrate out xN −1 ,
then xN −2 and so on. After completing k of the integrals we obtain
N −1
! 1/2
m N/2 Z Z
m X m (N −k−1)/2 m
2
dx1 . . . dxN −1 exp − (xn+1 − xn ) = · ×
2πiǫ~ 2iǫ~ n=0 2πiǫ~ 2πi[k + 1]ǫ~
N −k−2
!
m X m
Z Z
2 2
× dx1 . . . dxN −k−1 exp − (xn+1 − xn ) − (xN − xN −k−1 )
2iǫ~ n=0 2i[k + 1]ǫ~
ip2 t
Z
−1
= (2π~) dp exp − + ip(x0 − xN )/~
2m~
m 1/2 m
= exp − (xN − x0 )2
2πi~t 2i~t
— confirming our earlier result, Eq. (1.90).
In terms of the matrix element of the Boltzmann factor for which we have a path integral representation, this
thermal average is
dx x2 hx|e−βH |xi
R
2
hx i = R . (1.91)
dx hx|e−βH |xi
For the harmonic oscillator S[(τ )] is quadratic in x(τ ), and so the path integral is Gaussian. To evaluate it, we
should diagonalise the Euclidean action, which is
Z β~ " 2 #
m dx(τ ) κ
S[x(τ )] = dτ + x(τ )2 .
0 2 dτ 2
Since x(τ ) is real, φ−n = (φn )∗ , and so we can take as independent variables the real and imaginary parts of φn
for n > 0, together with φ0 which is real. The Jacobian for this change of variables is unity, because with our
choice of normalisation the Fourier transform is an expansion of x(τ ) in an orthonormal basis. In this way we have
" 2 #
2πn 1
Z Z YZ Z X
D[x(τ )] ≡ dφ0 dℜφn dℑφn and S[x(τ )] = m + κ |φn |2 + κφ20 .
n>0 n>0
~β 2
p
Averaging the amplitudes φn with the weight e−S[x(τ )]/~, and introducing the oscillator frequency ω = κ/m,
we have
~/κ
hφm φ−n i = δm,n .
1 + (2πn/β~ω)2
We obtain finally (using a contour integral to evaluate a sum on n)
* !2 +
1 X 1 X 1 ~ω/κ ~ω 1 1
hx2 (0)i = φn = = = + .
β~ n
βκ n 1 + (2πn/β~ω)2 2 tanh(β~ω/2) κ eβ~ω − 1 2
(1.92)
The rightmost form of the result is recognisable as the Planck formula for the average energy in the oscillator
multiplied by κ−1 , which is to be expected since this energy has equal kinetic and potential contributions, and the
potential energy is κhx2 i/2.
To summarise, we’ve succeeded in writing the Boltzmann factor for the harmonic oscillator as a path integral,
and in evaluating this path integral, and the result we’ve obtained matches what is familiar from more elementary
treatments. Of course, the approach seems rather laborious for an elementary problem like this one, but it brings
real advantages for more advanced ones.
Be warned that, except for the first one, these references were not written with undergraduate readers in mind.
18 CHAPTER 1. PATH INTEGRALS IN QUANTUM MECHANICS
Chapter 2
2.1 Introduction
In this chapter we turn to a class of problems from classical statistical physics known as stochastic processes. We
will discuss stochastic processes from several points of view, including the path integral formulation.
Recalling our discussion of path integrals in quantum mechanics, it is worth stressing a distinction between the
path integral for the time evolution operator and the one for the Boltzmann factor. For the time evolution operator,
the weight for a path x(t) is the phase factor exp(iS[x(t)]/~), while for the Boltzmann factor, the weight for a path
x(τ ) is the real, positive quantity exp(−S[x(τ )]/~). In this second case, we have the possibility of viewing the
weights as probabilities attached to the paths. For stochastic processes, one literally has probabilities associated
with various paths that a system may follow as a function of time. A path integral treatment of stochastic processes
therefore gives us a context in which functional integrals take on a very clear and concrete meaning.
An example of a stochastic process that we will discuss in some detail is Brownian motion: a small particle,
suspended in a fluid and observed under a microscope, moves in a random way because it is buffeted by collisions
with water molecules. The theory of Brownian motion, developed by Einstein in 1905, gave important confirmation
of ideas of statistical mechanics and a way to determine from observations the value of Avogadro’s number.
2.2.1 Moments of X
Z
hX n i = xn PX (x) dx
mean hXi
variance hX 2 i − hXi2 ≡ σ 2
standard deviation σ
19
20 CHAPTER 2. STOCHASTIC PROCESSES AND PATH INTEGRALS
1 dn
hX n i = φX (k) .
in dk n k=0
2.2.3 Cumulants
Write as a definition of the cumulants Cn (X)
!
X (ik)n
φX (k) = exp Cn (X) .
n
n!
Then
1 dn
Cn (X) = ln(φX (k)) .
in dk n k=0
We have
C1 (X) = hXi
C2 (X) = hX 2 i − hXi2 etc
In a rough sense, the size of the nth cumulant tells you how much the nth moment deviates from the value you
would have expected for it, given only knowledge of the lower order moments.
YN = X1 + X2 + . . . + XN .
We can think of N as labelling discrete time, and of YN as a random walk in one dimension. The mean of YN is
hYN i = N hXi
and so
Cm (SN ) = Cm (X) · N 1−m .
This shows that for large N , only C1 (SN ) and C2 (SN ) are non-vanishing. Such a feature of the cumulants is
unique to a Gaussian distribution, and so SN must be Gaussian distributed.
To demonstrate explicitly that this is the case, we can consider the calculation of the probability distribution for
SN , which we write as PSN (s). It can be obtained as the inverse Fourier transform of the characteristic function
for SN , and so we have Z ∞
1
PSN (s) = e−iks φSN (k)dk . (2.1)
2π −∞
Now, taking account only of the first two cumulants, we have
C2 (X)k 2
|φSN (k)| ≈ exp − .
2N
√
This makes it clear that only the range |k| . N gives a significant contribution to the integral. But in this range,
the corrections involving higher cumulants are small: they involve
m
k
N· . N 1−m/2 ,
N
which goes to zero as N → ∞ for m ≥ 3. Retaining only C1 (SN ) and C2 (SN ) in φSN (k), we can evaluate the
Fourier transform to obtain
1/2
N (s − hXi)2
N
PSN (s) ≈ exp − .
2πσ 2 2σ 2
Note that what we have done is a steepest descents evaluation of Eq. (2.1).
since Y must take some value at t2 . Conditional and unconditional probabilities are related by
P2 (y1 , t1 ; y2 , t2 ) = P1|1 (y2 , t2 |y1 , t1 ) · P1 (y1 , t1 ) .
The idea of conditional probabilities generalises in the obvious way: the quantity
P1|n−1 (yn , tn |y1 , t1 ; y2 , t2 . . . yn−1 , tn−1 )
is the probability to find Y = yn at tn , given that Y = y1 at t1 , Y = y2 at t2 etc.
and the evolution equation, which is the condition that, for any t2 > t1
Z
P1 (y2 , t2 ) = dy1 P1|1 (y2 , t2 |y1 , t1 )P1 (y1 , t1 ) .
2.4. MARKOV CHAINS 23
2.4.1 An example
Consider a two-state process, in which Y = 1 or Y = 2. Evolution from time t to time t + 1 can be specified as
follows:
a(t + 1) = T · a(t)
with
1−q r
T= .
q 1−r
Note that, in general, T is a square matrix with non-negative entries, but that it is not necessarily symmetric.
Hence hLβ |Rα i = 0 for λβ 6= λα , while by choice of normalisation (and of basis in the case of degeneracy) we
can set hLβ |Rα i = δαβ , so including the case λβ = λα .
A useful idea is that of spectral decomposition, meaning that we can write M in terms of its eigenvectors and
eigenvalues, as X
M= |Rα iλα hLα | .
α
To avoid divergences, one requires |λα | ≤ 1 for all α. in addition, for a large class of problems λ1 is the unique
eigenvalue with largest magnitude, so that |λα | < 1 for α 6= 1. Then in the limit n → ∞ we have Tn → |R1 ihL1 |.
From this we can conclude that the limiting distribution is given by |R1 i.
where the term involving the constant γ represents viscous drag exerted by the fluid on the particle, which tends to
reduce its velocity, and the function η(t) represents a random force on the particle, due to collisions with molecules
of the fluid. Since η(t) is random, it can by characterised by its correlation functions. In view of the central limit
theorem, it is often sufficient to know just the first two moments. We take this force to be zero on average, and take
its values at times separated by more than a microscopic interval to be uncorrelated, so that
In this way, the parameter Γ characterises the strength of the noise. Note that averages over realisations of the
noise are given by a Gaussian functional integral:
R
Z
1 2
h. . .i ≡ D[η(t)] . . . e− 2Γ dt η (t) .
Now, we can solve Eq. (2.2) explicitly for any noise history η(t) in the standard way for a first-order ordinary
differential equation, using an integrating factor. We find
Z t
′
−γt
v(t) = v(0)e + dt′ e−γ(t−t ) η(t′ ) (2.3)
0
which can be checked by substituting back. From this solution, we can calculate averages. We find
hv(t)i = v(0)e−γt
and
Z tZ t
h[v(t)]2 i = [v(0)]2 e−2γt + e−2γt dt1 dt2 eγ(t1 +t2 ) hη(t1 )η(t2 )i
0 0
Γ
= [v(0)]2 e−2γt + [1 − e−2γt ] , (2.4)
2γ
and also
Z t−τ /2 Z t+τ /2
lim hv(t − τ /2)v(t + τ /2)i = lim e−2γt dt1 dt2 eγ(t1 +t2 ) hη(t1 )η(t2 )i
t→∞ t→∞ 0 0
Γ −γ|τ |
= e . (2.5)
2γ
particular hx2 (t)i. One obvious approach to this calculation is to obtain x(t) in terms of η(t) by integrating the
expression for v(t) given in Eq. (2.3), and to average moments of x(t) over realisations of η(t). An alternative is
to derive and solve equations for the time evolution of the moments, as follows.
Consider
d 2 d
hx (t)i = 2hx(t) x(t)i .
dt dt
Differentiating again, we have
d d d
hx(t) x(t)i = hx(t) v(t)i + h[v(t)]2 i .
dt dt dt
Now, we can simplify the first term on the right side of this equation by using the Langevin equation to substitute
for v(t), and we know the second term from equipartition. In this way we obtain
d d d kB T
hx(t) x(t)i = −γhx(t) x(t)i + + hx(t)η(t)i . (2.6)
dt dt dt m
Since hx(t)η(t)i = 0, we can integrate Eq. (2.6), obtaining
d kB T
hx(t) x(t)i = Ce−γt + ,
dt γm
d kB T
hx(t) x(t)i = [1 − e−γt ] .
dt γm
To appreciate this result, it is useful to think about limiting cases. At short times (γt ≪ 1)
2kB T 1 (γt)2 kB T 2
hx2 (t)i ≈ [t − (γt − . . .)] = t ,
γm γ 2 m
which, reasonably enough, is ballistic motion at the mean square speed expected from equipartition. At long times
(γt ≫ 1)
2kB T
hx2 (t)i = t.
γm
In this case, the mean square displacement grows only linearly with time, which we recognise as diffusive motion.
The diffusion coefficient, 2kB T /γm involves system-specific constants γ and m: as indicated above, γ can be
calculated in terms of the size of particles and the viscosity of the fluid, and m can likewise be determined inde-
pendently. A measurement of this diffusion constant therefore constitutes a measurement of Boltzmann’s constant
kB , which is related to the gas constant R via Avogadro’s number. Since R is known from elementary measure-
ments on nearly ideal gases, we have a determination of Avogadro’s number. The relevant theory was published
by Einstein in 1905, and early experimental results were provided by Perrin in 1910.
so that the velocity change ∆v ≡ v(t + ∆t) − v(t) is Gaussian ditributed if η(t′ ) is. Now, following the general
approach to Markov processes, as introduced above, the probability distribution P1 (v, t) for the velocity v of a
Brownian particle at time t satisfies the (integral) evolution equation
Z
P1 (v, t + ∆t) = duP1|1 (v, t + ∆t|u, t)P1 (u, t) . (2.8)
We would like to get from this to a differential equation, by taking the limit of small ∆t, in which we expect that
|v − u| will also be small. Some care is necessary in the derivation, because we should think of P1|1 (v, t + ∆t|u, t)
as giving the probability distribution of v, for a fixed u, while in Eq. (2.8) the quantity v is fixed and the integral is
on u. To deal with this, we change variables from u to ∆v = v − u, and then Taylor expand in ∆v, obtaining
Z
P1 (v, t + ∆t) = d(∆v) P1|1 (v, t + ∆t|v − ∆v, t)P1 (v − ∆v, t)
Z
= d(∆v) P1|1 (v − ∆v + ∆v, t + ∆t|v − ∆v, t)P1 (v − ∆v, t)
∞
(−∆v)n ∂ n
Z X
= d(∆v) P1|1 (v + ∆v, t + ∆t|v, t)P1 (v, t) .
n=0
n! ∂v n
Note in the middle line of this equation the substitution v = v − ∆ + ∆v, used so that the combination v − ∆v
appears uniformly in the integrand, facilitating the Taylor expansion. The moments of ∆v which appear here,
Z
d(∆v)(∆v)n P1|1 (v + ∆v, t + ∆t|v, t) ≡ h(∆v)n i ,
Hence
∂ Γ ∆t ∂ 2
P1 (v, t + ∆t) = P1 (v, t) + γ∆t vP1 (v, t) + P1 (v, t) + O(∆t2 )
∂v 2 ∂v 2
and hence
Γ 2
∂t P1 (v, t) = γ∂v vP1 (v, t) + ∂ P1 (v, t) . (2.9)
2 v
which is the Fokker-Planck equation. Its solution is simplest to find in the limit t → ∞, because then P1 (v, t) is
independent of initial conditions and of t. It is simple to show that
γ 1/2
exp −γv 2 /Γ ,
P1 (v, ∞) =
πΓ
which is the Maxwell distribution familiar from kinetic theory (recall from the fluctuation-dissipation relation that
γ/Γ = m/(2kB T )).
dx(t)
= v(t) . (2.10)
dt
From Eq. (2.5) the moments of the velocity, for t ≫ γ −1 , are
Γ −γ|t−t′ | Γ
hv(t)i = 0 and hv(t)v(t′ )i = e ∼ 2 δ(t − t′ )
2γ γ
28 CHAPTER 2. STOCHASTIC PROCESSES AND PATH INTEGRALS
′
where the replacement of (γ/2)e−γ|t−t | by δ(t − t′ ) is appropriate if other quantities that depend on t − t′ vary
much more slowly. Clearly, Eq. (2.10) is like the Langevin equation, Eq. (2.2), but with x(t) replacing v(t) and
with γ = 0. We can therefore adapt many of the results we have already derived. In particular, by comparison with
the Fokker-Planck equation, Eq. (2.9), the evolution equation for P1 (x, t), the probability distribution for position,
is
Γ
∂t P1 (x, t) = 2 ∂x2 P1 (x, t) , (2.11)
2γ
which we recognise as the diffusion equation, with diffusion coefficient D = Γ/2γ 2. Its solution, with an initial
condition that the particle is at position x0 at time t0 , is
[x − x0 ]2
−1/2
P1 (x, t) ≡ P2 (x, t; x0 , t0 ) = (2πD|t − t0 |) exp − . (2.12)
4D|t − t0 |
γ2
Z
exp − dt[v(t)]2 .
2Γ
This is equivalent to the Euclidean path integral for a free particle, on making the identification
γ2 1 m
≡ = .
Γ 2D ~
Alternatively, we can consider evolution through a sequence of discrete intervals, dividing a total time t into N
equal timeslices, each of duration ǫ = t/N . Extending Eq. (2.12), we have
N −1 2 !
−N/2 ǫ X xn+1 − xn
Pn (x0 , 0; x1 , ǫ; . . . : xn , nǫ; . . . xN , t) = (4πDǫ) exp − , (2.13)
4D n=0 ǫ
which diverges as ǫ → 0. In other words, paths are continuous (since h(xn+1 − xn )2 i → 0 as ǫ → 0), but not
differentiable. They are in fact fractal objects, with dimension two, in contrast to one-dimensional smooth curves.
• F. Reif, Fundamentals of Statistical and Thermal Physics (McGraw-Hill). Chapter 15 covers Master equa-
tions, Langevin equations and Fokker-Planck equations.
• N. G. van Kampen, Stochastic Processes in Physics and Chemistry (North-Holland). A reasonably complete
and quite advanced reference book.
• L. S. Schulman, Techniques and Applications of Path Integration (Wiley). Chapter 9 gives a careful account
of the link between Brownian motion and the so-called Wiener integral.
Chapter 3
In this chapter we treat some problems from statistical mechanics that are selected partly because the mathematical
methods used to study them have close links with the ideas we have met in the context of the path integral formu-
lation of quantum mechanics and in the theory of stochastic processes. One common feature of the systems we
discuss, which makes their behaviour interesting, is that they are built from particles (or other microscopic degrees
of freedom) which have interactions between them. In this crucial sense, they differ from the simplest problems
of statistical mechanics, which involve only non-interacting particles: the ideal classical gas of kinetic theory, and
the ideal Fermi and Bose gases of elementary quantum statistical mechanics. Interparticle interactions can lead
to behaviour quite different from that of non-interacting systems, with the possibility of spontaneous order and
symmetry breaking at low temperature, and phase transitions as temperature is varied. In general, it is impossible
or very difficult to treat the statistical mechanics of interacting systems exactly. It can be done, however, for a
variety of one-dimensional models, and the necessary techniques are introduced in this chapter: the criterion of
tractability is the second reason for the selection of problems we make here. As we will see, it is a general feature
of one-dimensional systems that they do not show phase transitions or symmetry-breaking at non-zero temperature.
Nevertheless, they serve to show how interactions can have a controlling influence. We will return to the topic of
phase transitions in interacting systems in higher dimensions, using approximate methods, later in the course. The
one-dimensional models we define in this chapter all have obvious extensions to higher dimensions.
It is sometimes useful to specify periodic boundary conditions for a system of N sites, by setting SN +1 ≡ S1 and
taking the sums on i in Eq. (3.1) to run from 1 to N .
29
30 CHAPTER 3. STATISTICAL MECHANICS IN ONE DIMENSION
In a ground state all spins have the same orientation, so that the first term of H is minimised. With h = 0 there
are two such states: (i) Si = +1 for all i, and (ii) Si = −1 for all i. Hence, at zero temperature, when the system
is in a ground state, there is long-range order. Also, at h = 0, the model has symmetry under global spin inversion
(the energy is unchanged if the directions of all spins are reversed via Si → −Si ): this symmetry is broken in each
of the ground states. With h 6= 0 there is a unique ground state: Si = sign(h).
Both the Ising model and the lattice gas have microscopic degrees of freedom with two possible states, and one
can be mapped onto the other: for a lattice in which each site has z neighbours (z=2 in one dimension) the
replacements Si = 2ni − 1, 4J = V , and 2h − 2zJ = µ give HIsing = Hlattice gas + constant.
~i+1 − ~h ·
~i · S ~i .
X X
H = −J S S (3.3)
i i
~i , writing
For the XY model we can alternatively use an angle θi to specify the orientation of spin S
X X
H = −J cos(θi+1 − θi ) − h cos(θi ) . (3.4)
i
At ~h = ~0 both these models are symmetric under global rotations of the spins; they both also have continuous sets
of ground states, in which all spins are aligned but the common direction is arbitrary.
where we use the notation δ(σi , σi+1 ) for a Kronecker delta. With ferromagnetic interactions (J > 0), all σi are
the same in a ground state and there are q such states. For q = 2 the Potts model is equivalent to the Ising model.
where the potential V (ϕ(x)) is chosen to have minima at ϕ(x) = ±ϕ0 , and in its simplest form is
t 2 u
V (ϕ(x)) = ϕ (x) + ϕ4 (x) . (3.7)
2 4
We require u > 0 so that the energy has a lower bound; if in addition t > 0, then V (ϕ(x)) has a single minimum
at ϕ(x) = 0, whereas for t < 0 it has two minima at ϕ(x) = ±(−t/u)1/2 .
Here, the sum on states indicates literally a sum on discrete states for the Ising and Potts models, and multiple
integrals on spin orientations at each site for the lattice XY and Heisenberg models, while for continuum models it
denotes a functional integral over field configurations. We will be concerned with thermal averages of observables:
averages over configurations, weighted with the Boltzmann factor. We use the notation
X
h. . .i = Z −1 . . . e−βH (3.10)
states
where . . . stands for the observable. One example is the internal energy E, the average energy of the system, which
can be calculated from the partition function via
∂
E ≡ hHi = − ln(Z) . (3.11)
∂β
Other thermodynamic quantities can also be obtained from the partition function. In particular, the (Helmholtz)
free energy is
F = −kB T ln(Z) (3.12)
and the entropy S is
E−F
S= . (3.13)
T
32 CHAPTER 3. STATISTICAL MECHANICS IN ONE DIMENSION
When we come to develop intuition about the behaviour of models in statistical mechanics, it is useful to
remember the expression for the entropy of a system in the microcannonical ensemble, with W accessible states:
S = kB ln(W ). In addition, it is helpful to recall that free energy F is minimised for a system in equilibrium.
To characterise the behaviour of interacting systems, we will be particularly concerned with thermal averages
of products of microscopic variables. For example, for the Ising model we are interested in the magnetisation hSi i
and the two-point correlation function hSi Sj i. The partition function, viewed as a function of external field h, is
the generating function for these correlation functions, provided we allow the external field to take independent
values hi at each site. In particular, for the Ising model we have
1 ∂ 1 ∂2
hSi i = ln(Z) and hSi Sj i − hSi ihSj i = ln(Z) . (3.14)
β ∂hi β 2 ∂hi ∂hj
Moreover, the magnetisation m and the magnetic susceptibility χ can be obtained from derivatives with respect to
the strength of an external field that is uniform: for a system of N sites
N N
1 X 1 ∂ ∂m β X
m≡ hSi i = ln(Z) and χ≡ = [hSi Sj i − hSi ihSj i] . (3.15)
N i=1 βN ∂h ∂h N i,j=1
In a similar way, for the continuum theories we have introduced, we allow the external field to be a function
of position, h(x), so that the partition function is a functional, and then functional derivatives give correlation
functions, as in Chapter 1.
The dependence of the two-point correlation function on separation provides a measure of the influence of
interactions on the behaviour of the system. In the high temperature limit (β → 0) all configurations are weighted
equally in a thermal average, spins fluctuate independently at each site, and the correlation function is short range,
meaning that it falls to zero at large separation (in fact, in this limit hSi Sj i=0 unless i=j). In the opposite limit
of zero temperature (β → ∞), the system is in its ground state: as we have seen, all spins then adopt the same
configuration and the correlation function is infinite-ranged. We would like to understand in detail how behaviour
interpolates between these two limiting cases as temperature varies.
The Boltzmann factor, being the exponential of this sum, is a product of terms:
Y
e−βH = T (Si , Si+1 ) with T (Si , Si+1 ) = e−βH(Si ,Si+1 ) . (3.16)
i
Because interactions couple only adjacent slices, only two terms in this product depend on a given Si : T (Si−1 , Si )
and T (Si , Si+1 ). Moverover, summation on Si acts just like matrix multiplication. This leads us to define the
transfer matrix T: for the case in which Si takes M values (M =2 for the one-dimensional Ising model with only
nearest-neighbour interactions) it is an M × M matrix with rows and columns labelled by the possible values of
Si and Si+1 , and matrix elements T (Si , Si+1 ) as defined in Eq. (3.16). For a system of N slices and periodic
boundary conditions (so that i runs from 1 to N ), the partition function is simply the matrix trace
Z = Tr TN .
Alternatively, with fixed configurations for S1 and SN , and without periodic boundary conditions, the partition
function is the matrix element
Z = TN −1 (S1 , SN ) .
3.5. TRANSFER MATRIX SOLUTION OF THE ISING MODEL IN ONE DIMENSION 33
The advantage of the transfer matrix approach is that in this way calculations are reduced to study of an M ×M
matrix, independent of system size N . We will assume for simplicity that H(Si , Si+1 )=H(Si+1 , Si ) so that T is
symmetric and its eigenvectors |αi can be chosen to form a complete, orthonormal set. We order the associated
eigenvalues by magnitude (the first is in fact necessarily positive and non-degenerate): λ0 > |λ1 | ≥ . . . |λM−1 |,
with X
T|αi = λα |αi , hα|βi = δαβ and T= |αiλα hα| .
α
Using this notation, it is simple to write down a power of the transfer matrix: we have
X
TN = |αiλN α hα| .
α
The free energy per slice with periodic boundary conditions is therefore
M−1
!
F kB T kB T X
=− ln(Z) = −kB T ln(λ0 ) − ln 1 + [λα /λ0 ]N .
N N N α=1
In the thermodynamic limit (N → ∞), [λα /λ0 ]N → 0 for α ≥ 1, and so the free energy density f = limN →∞ F/N
is simply
f = −kB T ln(λ0 ) . (3.17)
Clearly, we can obtain other thermodynamic quantities, including the energy density E/N , the entropy density
S/N , the magnetisation m and magnetic susceptibility χ by this route.
To determine correlation functions, one might imagine we should first evaluate a generating function dependent
on field values hi at each site. Within the transfer matrix apporach, however, this is not convenient quantity to
consider, because if field values vary at different sites, the transfer matrices are different for each slice and transfer
matrix products no longer have simple expressions in terms of powers of the eigenvalues. Instead, we extend the
transfer matrix formalism by defining diagonal matrices C with diagonal elements C(Si , Si ) that are functions of
the degrees of freedom within a slice, chosen to reproduce the required correlation function. For example, for the
one-dimensional Ising model with N sites and periodic boundary conditions, taking C(Si , Si ) = Si we have
As happened for the free energy density, these expressions simplify greatly in the thermodynamic limit when
written in terms of the eigenvalues and eigenvectors of the transfer matrix: they reduce to
Tr Ti CTN −i
lim = h0|C|0i (3.18)
N →∞ Tr [TN ]
and x
Tr Ti CTx CTN −i−x
X λα
lim = h0|C|αihα|C|0i . (3.19)
N →∞ Tr [TN ] α
λ0
In summary, diagonalisation of the transfer matrix provides a path to calculating all quantities of physical interest.
For simplicity, we will set h=0. Then the eigenvalues are λ0 = 2 cosh(βJ) and λ1 = 2 sinh(βJ), and the
eigenvectors are |0i = 2−1/2 (1, 1)T and |1i = 2−1/2 (1, −1)T . The matrix representing the spin operator is
1 0
C= ,
0 −1
and its matrix elements in the basis of eigenstates are h0|C|0i = h1|C|1i = 0 and h0|C|1i = h1|C|0i = 1.
We are now in a position to write down some results. From Eq. (3.17), the free energy density is
As a check, we should examine its behaviour in the high and low temperature limits. At high temperature (β → 0),
f ∼ −kB T ln(2). This is as expected from the relation between free energy, energy and entropy, Eq. (3.13): in the
high-temperature limit neighbouring spins are equally likely to be parallel or antiparallel, and so the entropy per
spin (from the general formula S = kB ln(W ), with W the number of accessible states) is S = kB ln(2), while the
average energy is hHi=0. Conversely, at low temperature (β → ∞), f = −J, which arises because in this limit
the system is in a ground state, with neighboring spins parallel, so that hHi = −N J and S = 0. Beyond these
checks, the most interesting and important feature of our result for the free energy density is that it is analytic in
temperature for all T > 0. As we will discuss later in the course, phase transitions are associated with singularities
in the free energy density as a function of temperature, and analyticity of f in the one-dimensional Ising model
reflects the absence of a finite-temperature phase transition. The model in fact has a critical point at T = 0, and f
is non-analytic there (compare the limits T → 0+ and T → 0− ).
What happens in the model at low temperatures is most clearly revealed by the form of the correlation functions,
although for this we have to go to the two-point function. The one-point function, or magnetisation is trivial:
hSi i = h0|C|0i = 0 ,
which is a consequence of symmetry at h = 0 under global spin reversal. The two-point correlation function
between spins separated by a distance x is
|x|
λ1 1
hSi Si+x i = = exp(−|x|/ξ) with = ln[coth(βJ)] .
λ0 ξ
We see that correlations decay exponentially with separation, with a lengthscale ξ. This lengthscale, termed the
correlation length, increases monotonically with decreasing temperature and diverges as T → 0. Its asymptotic
form at low temperature (βJ ≫ 1) is
1 1 1
ξ= ∼ ∼ e2βJ . (3.21)
ln[coth(βJ)] ln[1 + 2e−2βJ ] 2
Useful physical insight into this result comes from a simple picture of typical states at low temperature. It is
clear that they consist of long sequences of parallel spins, with occasional reversals in orientation, as in Fig. 3.1.
In these circumstances it is natural to focus on the reversals, called domain walls or kinks, as the elementary
excitations in the system. Their average spacing sets the correlation length. Moreover, each domain wall costs
an energy 2J, the difference between the two possible values of JSi Si+1 . It also has an entropy associated with
it, since it can be placed between any one of O(1/ξ) neighbouring pairs of spins, without disrupting the overall
arrangement of irregularly spaced kinks. An estimate of the free energy for a chain of L sites is therefore
L L
F ∼ 2J − kB T ln(ξ) . (3.22)
ξ ξ
The actual value of ξ at a given temperature can be estimated as the one that minimises Eq. (3.22), yielding
ξ ∼ exp(2βJ) for βJ ≫ 1, in reasonable agreement with our earlier, detailed calcuation.
3.6. STATISTICAL MECHANICS IN ONE DIMENSION AND QUANTUM MECHANICS 35
over functions ϕ(x) satisfying the boundary conditions. Referring back to Eqns. (1.87) and (1.88), we see that βH
in the classical statistical mechanics problem plays the same role as the Euclidean action S/~ in a path integral
expression for the Boltzmann factor arising in a quantum problem: we use HQ to denote the Hamiltonian of the
quantum problem. The translation dictionary is as follows.
The action is that for a particle of mass J moving in a potential V (x): reversing the steps of section 1.2.3, we can
read off the quantum Hamiltonian as
1 d2
HQ = − 2 + V (ϕ) . (3.23)
2β J dϕ2
Knowledge of the eigenvalues ǫα and eigenfunctions |αi of HQ , which satisfy HQ |αi = ǫα |αi , gives access to
thermodynamic quantities and correlation functions for the classical system. In particular, for a classical system
extending over −∞ < x < ∞, the arguments leading to Eq. (3.19) also give
X
hϕ(x1 )ϕ(x2 )i = h0|ϕ|αihα|ϕ|0i e−β(ǫα −ǫ0 )|x1 −x2 |
α
Now, although we cannot find the eigenfunctions of HQ exactly, we know quite a lot about them for the case of
interest, in which V (ϕ) is a quartic double-well potential. In particular, since V (−ϕ) = V (ϕ), all eigenfunctions
have definite parity. The ground state wavefunction, hϕ|0i has even partiy and is nodeless, while the first excited
state wavefunction, hϕ|1i, is odd, having a single node at ϕ = 0. From this it follows that h0|ϕ|0i = 0 and
h0|ϕ|1i = 6 0. The inverse correlation length, governing the correlation function at large separation |x1 − x2 | is
therefore ξ −1 = β(ǫ1 − ǫ0 ). At low temperatures, the form of the lowest and first excited eigenstates is as sketched
in Fig. 3.2. In the language of quantum mechanics, the splitting between eigenvalues ǫ1 − ǫ0 arises because of
tunneling through the barrier between the two minima of V (ϕ) and is exponentially small for small T , which leads
to a correlation length exponentially large in β, as for the Ising model, Eq. (3.21).
1.5 1.5
0.25
0.2 1 1
0.15
0.05
0 0 0
−0.05
−0.15
−0.2 −1 −1
−0.25
−1.5 −1.5
−2 −1.5 −1 −0.5 0 0.5 1 1.5 2 −2 −1.5 −1 −0.5 0 0.5 1 1.5 2 −2 −1.5 −1 −0.5 0 0.5 1 1.5 2
Figure 3.2: Left: the potential V (ϕ) appearing in HQ , Eq. (3.23). Centre and right: form of the lowest two
eigenfunctions of HQ , hϕ|0i and hϕ|1i, for small T
The correlation functions we use to characterise behaviour of the system should be constructed in a way that
respects the periodic nature of the coordinate θ. An obvious candidate is hei[θ(x1 )−θ(x2 )] i. As in our previous
examples, we can express this in terms of the eigenfunctions and eigenvalues of HQ . We find
X
hei[θ(x1 )−θ(x2 )] i = h0|e−iθ |nihn|eiθ |0ie−β(ǫn−ǫ0 )|x1 −x2 | .
n
The correlation length for this model is therefore ξ = [β(ǫ1 − ǫ0 )]−1 = 2J/(kB T ). As for the Ising model,
ξ diverges in the limit T → 0, which is expected since for T = 0 the system adopts a ground state with θ(x)
indepndent of x. The divergence of ξ, howover, is much less rapid in the XY model than in the Ising model. The
reason for this is that, whereas excitations in the Ising model cost a minimum energy, the kink energy 2J, long
wavelength twists of θ(x) in the XY model can have arbitrarily low energy. As a consequence, thermal excitations
are more effective in disordering the XY model at low temperature than for the Ising model, leading to a shorter
correlation length in the XY model.
3.7. FURTHER READING 37
• J. J. Binney, N. J. Dowrick, A. J. Fisher, and M. E. J. Newman, The Theory of Critical Phenomena (OUP).
Chapter 3 gives a self-contained introduction to models and to transfer matrix methods.
38 CHAPTER 3. STATISTICAL MECHANICS IN ONE DIMENSION
Chapter 4
In this chapter, we will develop the Lagrangian approach to the classical theory of fields, focusing on field the-
ories with scalar and vector fields. Our discussion will culminate in the discussion of scalar electrodymanics, a
theory which couples Maxwell’s theory of electrodymanics to scalar fields. As we will see, symmetries play an
important role in constructing all these theories. Traditionally, symmetries of physical theories have often been
identified only after a proper mathematical formulation of the theory. For example, the fact that Maxwell’s the-
ory of electrodynamics is covariant under Lorentz transformations was only discovered significantly after its the
first formulation. In modern theoretical physics, this traditional relation between theories and their symmetries
is frequently reversed. One starts by identifying the symmetries of the given physical situation and then writes
down the (most general) theory compatible with these symmetries. This approach has been immensely fruitful
and has played a major role, for example, in constructing the standard model of particle physics. It is, therefore,
crucial to understand the relevant symmetries (or groups in Mathematical language) and the objects they act on
(representations, in Mathematical language) first.
4.1 Symmetries
4.1.1 Definition of groups and some examples
The word ”symmetry” in Physics usually (although not always) refers to the Mathematical structure of a group, so
we begin with the following definition.
Definition A group G is a set with a map · : G × G → G (”multiplication”) satisfying the three conditions
1) g1 · (g2 · g3 ) = (g1 · g2 ) · g3 for all g1 , g2 , g3 ∈ G (associativity)
2) There exists an e ∈ G such that g · e = g for all g ∈ G (neutral element)
3) For each g ∈ G there exist a g −1 ∈ G such that g · g −1 = e (inverse element)
It is easy to prove from the above axioms that the neutral element e is unique, that it is also the neutral element
when acting from the left, that is e · g = g for all g ∈ G, that the right-inverse g −1 is uniquely defined for each
g ∈ G and that it is also the left-inverse, that is g −1 · g = e. If, in addition to the three axions above, g1 · g2 = g2 · g1
is satisfied for all g1 , g2 ∈ G the group is called Abelian.
Well-known groups are the integers with respect to addition and the real and complex numbers with respect to
addition and multiplication. All these groups are Abelian. Here are some more interesting groups which will play
a role in our field theory constructions. Consider first the group ZN = {0, 1, · · · , N − 1} with ”multiplication”
defined by n1 · n2 = (n1 + n2 ) mod N . This group is obviously finite (that is, it has a finite number of elements)
and Abelian. Another Abelian example is given by the complex numbers of unit length, U (1) = {z ∈ C| |z| = 1},
with group multiplication the ordinary multiplication of complex numbers. Not only is this group infinite but, as
it corresponds to the unit circle in the complex plane, it is also ”continuous” and one-dimensional. Examples for
non-Abelian groups are provided by the unitary groups SU(n) which consist of all complex n × n matrices U
satisfying U † U = 1 and det(U ) = 1, with ordinary matrix multiplication as the group multiplication and the unit
matrix as the neutral element. Matrix multiplication does in general not ”commute” which causes the non-Abelian
character of the unitary groups. The simplest non-trivial example of a unitary group on which we will focus later
is SU(2). Solving the unitary conditions U † U = 1 and det(U ) = 1 by inserting an arbitrary 2 × 2 matrix with
39
40 CHAPTER 4. CLASSICAL FIELD THEORY
This shows that we can think of SU(2) as the unit sphere in four-dimensional Euklidean space and, hence, that
it is a three-dimensional continuous group. Such continuous groups are also called Lie groups in Mathematical
parlance and we will discuss somepof their elementary properties in due course. We can solve the constraint on α
and β in Eq. (4.1) by setting α = 1 − |β|2 eiσ and β = −β2 + iβ1 which leads to the explicit parameterization
p
1 − |β|2 eiσ p−β2 + iβ1
U= , (4.2)
β2 + iβ1 1 − |β|2 e−iσ
by real two-dimensional rotation matrices over the vector space V = R2 , where n = 0, . . . , N − 1. We can
restrict the value of q to the range 0, . . . , N − 1 (as two values of q which differ by N lead to the same set
of matrices) and this provides, in fact, a complete list of representations for ZN . Equivalently, we can write
down the same representations over a one-dimensional complex vector space V = C where they take the form
Rq (n) = exp(2πiqn/N ). If, for a given representation Rq , we denote elements of the vector space V = C by Φ
then the group acts on them as Φ → Rq (n)Φ = exp(2πiqn/N )Φ. In this case, Φ is said to have ”charge” −q in
physics language.
Representations Rq for U (1) (on V = C) are just as easily obtained by writing
where α ∈ [0, 2π]. For Rq to be continuous when going around the circle, q must be an integer, however, unlike in
the ZN case it is not otherwise restricted. The above representations Rq for q an arbitrary integer, in fact, provide
all (continuous) representations of U (1). As before, a (complex) field transforming as Φ → Rq (eiα )Φ is said
to have charge −q. Also note that charge q = 0 corresponds to the trivial representation. The groups SU(n)
are already given by matrices, so we can think about them as representing themselves. This representation is n-
complex dimensional and is also called the fundamental representation. On a complex vector Φ = (φ1 , . . . , φn ) is
acts as Φ → U Φ, where U ∈ SU(n). However, this is by no means the only representation of SU(n), in fact, there
is an infinite number of them, as we will see.
4.1. SYMMETRIES 41
There are a number of general ways of constructing new representations from old ones which should be men-
tioned. For a representation R : G → Gl(n) of a group G there is a complex conjugate representation R⋆ defined
by R⋆ (g) = R(g)⋆ , that is, each group element is now represented by the complex conjugate of the original
representation matrix. Applying this to SU(n) leads to the complex conjugate of the fundamental representation
U → U ⋆ . For two representations R1 and R2 of a group G with dimensions n1 and n2 one can consider the direct
sum representation R1 ⊕ R2 with dimension n1 + n2 defined by the block-diagonal matrices
R1 (g) 0
(R1 ⊕ R2 )(g) = . (4.5)
0 R2 (g)
A representation such as this is called reducible and, conversely, a representation which cannot be split into smaller
blocks as in (4.5) is called irreducible. For example, the direct sum representation
of U (1) consisting of a charge +1 and −1 representation realises an explicit embedding of U (1) into SU(2).
Another, less trivial way of combining the two representations R1 and R2 to a new one is the tensor representation
R1 ⊗ R2 with dimension n1 n2 defined by (R1 ⊗ R2 )(g) = R1 (g) × R2 (g) 1 . In general, a tensor representation
R1 × R2 is not irreducible and decomposes into a sum of irreducible representations R(i) , so one can write
M
R1 ⊗ R2 = R(i) (4.7)
i
The matrices Ti are called the generators of the Lie group and the vector space L(G) = {ti Ti } spanned by
these matrices is referred to as Lie algebra. In the case of SU(2), the generators are given by (i times) the Pauli
matrices, as differentiating Eq. (4.2) shows. In general, there is a theorem which states that the group (or, rather, a
neighbourhood of the identity of the group) can be reconstructed from the Lie algebra by the exponential map
Now consider two matrices M (t) and M (s) and the product
X
M (t)−1 M (s)−1 M (t)M (s) = 1 + ti sj [Ti , Tj ] + · · · , (4.10)
i,j
where [·, ·] is the ordinary matrix commutator. Since the product on the LHS of Eq. (4.10) is an element of the
group, we conclude that the commutators [Ti , Tj ] must be elements of the Lie algebra and can, hence, be written
as
[Ti , Tj ] = fij k Tk . (4.11)
The coefficients fij k are called the structure constants of the Lie algebra L(G). More accurately, the Lie-algebra
L(G) is then the vector space L(G) = {ti Ti } together with the commutator bracket [·, ·]. The concept of a
representation can now also be defined at the level of the Lie algebra.
Definition A representation r of a Lie algebra L is a linear map which assigns to elements T ∈ L matrices r(T )
such that [r(T ), r(S)] = r([T, S]) for all T, S ∈ L.
1 For two matrices M and N the product M ×N can be thought of as the matrix obtained by replacing each entry in M by a block consisting
of that entry times N . A useful property of this product is (M1 × N1 )(M2 × N2 ) = (M1 M2 ) × (N1 N2 ).
42 CHAPTER 4. CLASSICAL FIELD THEORY
Note this is equivalent to saying that the representation matrices r(Ti ) commute in the same way as the generators
Ti , so [r(Ti ), r(Tj )] = fij k r(Tk ), with the same structure constants fij k as in Eq. (4.11). It is usually easier to
find representations of Lie algebras than representations of groups. However, once a Lie-algebra representation has
been found the associated group representation can be re-constructed using the exponential map (4.9). Concretely,
i i
for a Lie-algebra representation Ti → r(Ti ) the corresponding group representation is et Ti → et r(Ti ) . Recall
that the dimension of the representation r is defined to be the dimension of the vector space on which the repre-
sentation matrices r(T ) (or the associated group elements obtained after exponentiating) act, that is, it is given by
the size of the matrices r(T ). This dimension of the representation r is not to be confused with the dimension of
the Lie-algebra itself, the latter being the dimension of the Lie algebra L(G) as a vector space of matrices.
Example SU(2)
Let us see how all this works for our prime example SU(2). Consider an SU(2) matrix U close to the identity
matrix and write 2 U = 1 + iT + . . . , where T is a Lie algebra element. Then, evaluating the conditions U † U = 1
and det(U ) = 1 at linear level in T , one finds the constraints T † = T and tr(T ) = 0. In other words, the
Lie algebra L(SU(2)) of SU(2) consists of all traceless, hermitian 2 × 2 matrices. Note that this space is three-
dimensional. A convenient basis of generators τi for this Lie algebra is obtained from the Pauli matrices σi . Recall
that they satisfy the useful identities
1
τi = σi with [τi , τj ] = iǫijk τk , (4.13)
2
and the structure constants are simply given by the Levi-Civita tensor. While the dimension of the Lie algebra
L(SU(2)) is 3 (as it is spanned by three Pauli matrices), the dimension of the SU(2) representation defined by the
Pauli matrices is 2 (since they are 2 × 2 matrices).
Finite SU(2) matrices are then obtained by exponentiating
Note that the generator τ3 corresponds to the U (1) subgroup (4.6) of SU(2). The commutation relations (4.13) of
the SU(2) Lie algebra are identical to the commutation relations of the angular momentum operators in quantum
mechanics. Hence, we already know that the finite-dimensional representations of this algebra can be labelled
by a ”spin” j, that is an integer or half-integer number j = 0, 1/2, 1, 3/2, . . .. For a given j the dimension of the
representation is 2j +1 and the representation space is spanned by states |jmi, where m = −j, −j +1, . . . , j −1, j.
The two-dimensional representation for j = 1/2 of course corresponds to the explicit representation of the algebra
in terms of Pauli matrices which we have written down above. The complex conjugate of the fundamental is also
a two-dimensional representation and, on purely dimensional grounds, must also be identified with the j = 1/2
representation.
We also know from quantum mechanics that the tensor product of two representations characterised by j1 and
j2 decomposes into the irreducible representations with j in the range |j1 − j2 |, |j1 − j2 | + 1, . . . , j1 + j2 . This
is an explicit example of a Clebsch-Gordon decomposition. It is customary to refer to representations by their
dimensions, that is, write for example, the j = 1/2 representation as 2 (or 2̄ for the conjugate) and the j = 1
representation as 3. With this notation, examples of SU(2) Clebsch-Gordon decompositions are
Actions should be invariant under a symmetry group and, hence, it is of particular importance to understand the
singlets which occur in a Clebsch-Gordan decomposition. They will tell us about the invariant terms which are
allowed in an action. For example, if we have a field Φ which transforms as a doublet under SU(2), the first of
Eqs. (4.15) tells us that we should be able to write a quadratic term in Φ, corresponding to the direction of the
singlet on the right-hand side.
Example SO(3)
Another important Lie group is SO(3), the group of three-dimensional rotations, consisting of real 3 × 3 matrices
2 In the physics literature it is conventional to include a factor of i in front of the generators T .
4.1. SYMMETRIES 43
Table 4.1: The four disconnected components of the Lorentz group. The union L+ = L↑+ ∪ L↓+ is also called the
proper Lorentz group and L↑ = L↑+ ∪L↑− is called the orthochronos Lorentz group (as it consists of transformations
preserving the direction of time). L↑+ is called the proper orthochronos Lorentz group.
O satisfying OT O = 1 and det(O) = 1. Writing O = 1 + iT with (purely imaginary) generators T , the relation
OT O = 1 implies T = T † and, hence, that the Lie-algebra of SO(3) consists of 3 × 3 anti-symmetric matrices
(multiplied by i). A basis for this Lie algebra is provided by the three matrices Ti defined by
ΛT ηΛ = η . (4.18)
Special Lorentz transformations are the identity 14 , parity P = diag(1, −1, −1, −1), time inversion T = diag(−1, 1, 1, 1)
and the product P T = −14 . We note that the four matrices {14 , P, T, P T } form a finite sub-group of the Lorentz
group. By taking the determinant of the defining relation (4.18) we immediately learn that det(Λ) = ±1 for all
Lorentz transformations. Further, if we write out Eq. (4.18) with indices
where Ti are the generators (4.16) of the rotation group. Given the embedding (4.20) of the rotation group into
the Lorentz group the appearance of the Ti should not come as a surprise. It is straightforward to work out the
commutation relations
The above matrices can also be written in a four-dimensional covariant form by introducing six 4 × 4 matrices σµν ,
labelled by two anti-symmetric four-indices and defined by
By explicit computation one finds that Ji = 21 ǫijk σjk and Ki = σ0i . Introducing six independent parameters ǫµν ,
labelled by an anti-symmetric pair of indices, a Lorentz transformation close to the identity can be written as
i
Λρ σ ≃ δ ρ σ − ǫµν (σµν )ρ σ = δ ρ σ + ǫρ σ ; . (4.24)
2
The commutation relations (4.22) for the Lorentz group are very close to the ones for SU(2) in Eq. (4.13). This
analogy can be made even more explicit by introducing a new basis of generators
1
Ji± = (Ji ± iKi ) . (4.25)
2
In terms of these generators, the algebra (4.22) takes the form
that is, precisely the form of two copies (a direct sum) of two SU(2) Lie-algebras. Irreducible representations of the
Lorentz group can therefore be labelled by a pair (j+ , j− ) of two spins and the dimension of these representations
is (2j+ +1)(2j− +1). A list of a few low-dimensional Lorentz-group representations is provided in Table 4.2. Field
theories in Minkowski space usually require Lorentz invariance and, hence, the Lorentz group is of fundamental
importance for such theories. Since it is related to the symmetries of space-time it is often also referred as external
symmetry of the theory. The classification of Lorentz group representations in Table 4.2 provides us with objects
which transform in a definite way under Lorentz transformations and, hence, are the main building blocks of
such field theories. In these lectures, we will not consider spinors in any more detail but focus on scalar fields
φ, transforming as singlets, φ → φ under the Lorentz group, and vector fields Aµ , transforming as vectors,
Aµ → Λµ ν Aν .
classical mechanics is a function of the generalised coordinates and their first (time) derivatives. Analogously, we
start with a field theory Lagrangian density L = L(φa , ∂µ φa ) which is a function of the fields φa and their first
space-time derivatives ∂µ φa . The field theory action can then be written as
Z
S = d4 x L(φa (x), ∂µ φa (x)) , (4.27)
where the integration ranges over all of Minkowski space. Our first task is to derive the Euler-Lagrange equations
for such a general field theory by applying the variational principle to the above action. One finds
δS δ ∂L δφb (x̃) ∂L δ∂µ φb (x̃)
Z Z
0 = = d4 x̃ L(φb (x̃), ∂µ φb (x̃)) = d4 x̃ + (4.28)
δφa (x) δφa (x) ∂φb δφa (x) ∂(∂µ φb ) δφa (x)
∂L ∂L δφb (x̃) ∂L ∂L
Z
= d4 x̃ − ∂µ = − ∂µ (x) (4.29)
∂φb ∂(∂µ φb ) δφa (x) ∂φa ∂(∂µ φa )
where Φaα (x) are functions which encode the particular type of symmetry action. We assume that these transfor-
mations leave the action (4.27) invariant and, hence, change the Lagrangian density by total derivatives only. This
means the transformation of the Lagrangian density is of the form
with certain functions Λµα (x) which can be computed for each type of symmetry action. Let us now compare this
variation of L with the one induced by transformation (4.34) of the fields. One finds
∂L ∂L
L → L+ Φaα ǫα + ∂µ Φaα ǫα (4.36)
∂φa ∂(∂µ φa )
∂L ∂L ∂L
= L + ǫα ∂µ Φaα − ǫα ∂µ − Φaα (4.37)
∂(∂µ φa ) ∂(∂µ φa ) ∂φa
The last term vanishes thanks to the Euler-Lagrange equations (4.31) and equating the remaining variation of L
with the one in Eq. (4.35) it follows that
∂L
∂µ jαµ = 0 where jαµ = Φaα − Λµα . (4.38)
∂(∂µ φa )
46 CHAPTER 4. CLASSICAL FIELD THEORY
Hence, for each symmetry generator ǫα we obtain a conserved current jαµ , that is, a current with vanishing diver-
gence. Each such current can be used to define a conserved charge Qα by
Z
Qα = d3 x jα0 (4.39)
Using Eq. (4.38) and assuming that fields fall off sufficiently rapidly at infinity it follows that
Z Z
Q̇α = d x ∂0 jα = d3 x∂i jαi = 0 ,
3 0
(4.40)
The role of the symmetry parameters ǫα is here played by the infinitesimal translations aν . Therefore, the index α
which appears in the general equations above becomes a space-time index ν. Under a translation the Lagrangian
density changes as
L → L + aµ ∂µ L = L + aν ∂µ (δνµ L) . (4.42)
Comparing the last two equations with the general formulae (4.34) and (4.35) we learn that Φaν = ∂ν φa and
Λµν = δνµ L. Inserting this into the general result (4.38) leads to four currents T µ ν = jνµ given by
∂L
T µν = ∂ν φa − δ µ ν L . (4.43)
∂(∂µ φa )
For a translation-invariant theory they are conserved, that is, they satisfy ∂µ T µ ν = 0. The tensor T µ ν is called the
energy-momentum tensor and its associated charges
Z
Pν = d3 x T 0 ν (4.44)
represent the conserved energy and momentum of the system. In particular, the conserved energy P0 is explicitly
given by Z
∂L
Z
P0 = d3 x ∂0 φa − L = d3 x H = H , (4.45)
∂(∂0 φa )
that is, by the Hamiltonian (4.33).
• Choose a group which corresponds to the symmetries of the theory. Normally, the symmetries include the
external symmetry, that is Lorentz symmetry. In addition, there may be internal symmetries which do not
act on space-time indices but internal indices. (We will see explicit examples of such internal symmetries
shortly.)
• Choose a set of representations of the symmetry group. This fixes the field content of the theory and the
transformation properties of the fields.
• Write down the most general action invariant under the chosen symmetry (with at most two derivatives in
each term) for the fields selected in the previous step. Normally, only polynomial terms in the fields are
considered and an upper bound on the degree of the polynomials is imposed (for example by requiring that
the theory does not contain (coupling) constants with negative energy dimension).
4.3. SCALAR FIELD THEORY 47
1 2 2 1 1
V = m φ + λ3 φ3 + λφ4 . (4.47)
2 3! 4!
(A possible linear term in φ can be removed through a re-definition of φ by a shift.) Note that m and λ3 have
dimension 1 and λ is dimensionless. The quadratic term in V is called a mass term with mass m and the other
terms represent couplings. Applying the Euler-Lagrange equations (4.31) to the above Lagrangian leads to the
equation of motion
2φ + V ′ (φ) = 0 (4.48)
for φ, where 2 = ∂µ ∂ µ and the prime denotes the derivative with respect to φ. For non-vanishing λ3 or λ solutions
to this equation are not easy to find.
Since the Fourier transform can be inverted we conclude that (k 2 − m2 )φ̃(k) = 0 and, hence, that φ̃ can be
written in the form φ̃(k) = δ(k 2 − m2 )ϕ̃(k) for some function ϕ̃. Inserting this result for φ̃ into the Fourier
transform (4.50) and using 3
1 p
δ(k 2 − m2 ) = (δ(k0 − wk ) + δ(k0 + wk )) , wk = k2 + m2 (4.52)
2wk
3 This
P 1
follows from the well-known delta-function identity δ(f (x)) = x0 :f (x0 )=0 |f ′ (x0 )| δ(x − x0 ).
48 CHAPTER 4. CLASSICAL FIELD THEORY
2wk δ 3 (k − q) . (4.55)
is Lorentz invariant as well and can be viewed as a covariant version of the three-dimensional delta function. Up
to this point we have, effectively, solved the Klein-Gordon equation for a complex scalar field. However, imposing
a reality condition on the solutions is easy and leads to a+ (k) = a− (k) ≡ a(k). The final result for the general
solution of the Klein-Gordon equation for a real scalar field then reads
Z
φ(x) = d3 k̃ a(k)e−ikx + a⋆ (k)eikx .
(4.56)
Our theory is also invariant under Lorentz transformations and from Noether’s theorem we expect associated con-
served currents which we will now derive. First, recall from Eq. (4.24) that an infinitesimal Lorentz transformation
on xµ can be written as xµ → xµ + ǫµ ν xν where ǫµν is anti-symmetric. On the field φ this transformation acts as
Comparing with Eqs. (4.34) and (4.35) we learn that the symmetry parameters ǫα are here given by ǫµν , so we
have to replace the index α in our general equations with an anti-symmetric pair of space-time indices. Further,
ρ
we have Φµν = 2x[µ ∂ν] φ and Λρµν = 2x[µ δν] ρ
L. Inserting this into Eq. (4.38) we find the conserved currents Mµν
given by
ρ
Mµν = xµ Tνρ − xν Tµρ , (4.63)
with the energy momentum tensor for a scalar field theory defined in Eq. (4.59). From Noether’s theorem we know
ρ
that these currents are divergence-free, ∂ρ Mµν = 0, and imply the existence of conserved charges
Z Z
Mµν = d3 x Mµν
0
= d3 x (xµ Tν0 − xν Tµ0 ) . (4.64)
V ′ (v) = 0 , (4.66)
so we are instructed to look at extrema of the scalar potential V . In fact, to minimise the energy (4.60) we should
be focusing on minima of the scalar potential V . We will also refer to such a solution of the classical theory as
a vacuum. If the quartic coupling λ is negative the scalar potential is unbounded from below and the energy of a
constant field configuration tends to minus infinity for large field values. To avoid such an unphysical situation we
assume that λ > 0 in the following. Then we should distinguish two cases which are illustrated in Fig. 4.1.
• m2 ≥ 0 : In this case there is a single minimum at φ = v = 0. This solution is mapped into itself under the
action φ → −φ of the Z2 symmetry and we say that the symmetry is unbroken in this vacuum.
• m2 < 0 : In this case, φ = 0 is a maximum of the potential and there are two minima at
r
−6m2
φ=v=± . (4.67)
λ
Neither minimum is left invariant under the Z2 action φ → −φ (in fact the two minima are mapped into each
other under Z2 ) and we say that the symmetry is spontaneously broken. In general, spontaneous breaking
of a symmetry refers to a situation where a symmetry of a theory is partially or fully broken by a vacuum
solution of the theory. The potential value at the minima is given by
1 1
V (v) = V0 + m2 v 2 = V0 − λv 4 . (4.68)
4 24
Just as the constant V0 which we have included earlier, the potential value at the minima does not affect
any of the physics discussed so far. However, if we couple our theory to gravity, it turns out that V (v) acts
like a cosmological constant Λ in the Einstein equations. Cosmological constraints tell us that Λ cannot
be much bigger than O(meV 4 ). On the other hand, there is no obvious constraint on V (v). Unless there
is a cancellation of the two terms in Eq. (4.68), one would expect V (v) to be of the order of the symmetry
breaking scale v to the fourth power. Electroweak symmetry is broken spontaneously by a mechanism similar
to the above (and we will study a model related to this in Section (4.3.4)) at a scale of v ∼ TeV. Hence,
the ”natural” cosmological constant which arises at electroweak symmetry breaking is about 60 orders of
magnitude larger than the observational limit. So, we have to assume that the two terms in Eq. (4.68) cancel
each other to a precision of 60 digits. This enormous ”fine tuning” is one of the manifestations of what is
referred to as the cosmological constant problem. The question of why the cosmological constant is as small
as it is is one of the most important unresolved problems in modern physics.
Φ
-v v
Figure 4.1: Shape of scalar potential (4.65) for m2 ≥ 0 (solid line) and m2 < 0 (dashed line). In the latter case the
position v of the minima is given by Eq. (4.67).
Such a symmetry, which does not act on space-time indices but on internal indices is also called an internal symme-
try. For now, we will study the case of global internal symmetries, that is symmetries for which the transformation
is the same everywhere in space-time. This means the group parameter α ∈ [0, 2π] is independent of the space-
time coordinates xµ . To discuss the two-scalar field theory with this global SO(2) symmetry and scalar charge q
explicitly it proves convenient to arrange the two real scalars into a single complex one
1
φ = √ (φ1 + iφ2 ) . (4.70)
2
On this complex scalar, the SO(2) symmetry acts via the charge q representation (4.4) of U(1), that is
φ → exp(−iqα)φ . (4.71)
The complex conjugate φ⋆ transforms as φ⋆ → exp(iqα)φ⋆ and, hence, corresponds to a representation with
charge −q. Allowed terms in the Lagrangian density have to be U(1) invariant which is equivalent to saying that
their total charge needs to be zero. For example, the term φ2 has total charge 2q and cannot appear while the term
φ⋆ φ has charge zero and is allowed. In general, we can only allow terms with the same number of φ and φ⋆ , so
that the general U(1) invariant Lagrangian density reads
λ ⋆ 2
L = ∂µ φ⋆ ∂ µ φ − V (φ, φ⋆ ) , V = V0 + m2 φ⋆ φ + (φ φ) . (4.72)
4
Note that it is essential for the invariance of the kinetic term that the group parameter α is space-time indepen-
dent, that is, that the symmetry is global. For the equation of motion for φ we find from the Euler-Lagrange
equation (4.31)
∂V λ
2φ + ⋆
= 2φ + m2 φ + (φ⋆ φ)φ = 0 . (4.73)
∂φ 2
For λ = 0 this is the Klein-Gordan equation for a complex scalar field whose general solution has already been
obtained in Eq. (4.53).
∂L ∂L
π= = ∂0 φ⋆ , π⋆ = = ∂0 φ . (4.74)
∂(∂0 φ) ∂(∂0 φ⋆ )
Being translation and Lorentz invariant the above theory has conserved energy-momentum and angular momentum
tensors which can be obtained in complete analogy with the single scalar field case in sub-section 4.3.1. In addition,
4.3. SCALAR FIELD THEORY 51
the presence of the internal U(1) symmetry leads to a new type of conserved current which we will now derive.
From Eq. (4.71), infinitesimal U(1) transformations are given by
φ → φ − iqαφ , φ⋆ → φ⋆ + iqαφ⋆ . (4.76)
Comparing with the general transformation (4.34) we conclude that α plays the role of the (single) symmetry
parameter and Φ = −iqφ, Φ⋆ = iqφ⋆ . Since the Lagrangian density is invariant under U(1) the total derivative
terms in Eq. (4.35) vanish and we can set Λµ to zero. Inserting this into the general formula (4.38) for the conserved
current we find
jµ = iq(φ⋆ ∂µ φ − φ∂µ φ⋆ ) . (4.77)
Spontaneous symmetry breaking
As we did before, we would now like to discuss the vacua of the theory, that is solutions to the equation of
motion (4.73) with φ = v = const. For m2 ≥ 0 there is a single minimum at φ = 0. This solution is left invariant
by the transformations (4.71) and, hence, the U(1) symmetry is unbroken in this case. For m2 < 0 the shape of
the potential is shown in Fig. (4.2). In this case, there is a whole circle of minima
Φ2
Φ1
∂2V
Mab = (v) . (4.83)
∂φa ∂φb
The eigenvalues of the mass matrix M are the mass squares of the fields around the vacuum v. Now let us assume
that our scalar field theory is invariant under a continuous symmetry group G and that the scalar fields transform
as φ → R(g)φ under the representation R of G. In particular, this means the scalar potential is invariant, that is
for all g ∈ G. The vacuum v will in general not respect the full symmetry group G but will spontaneously break
it to a sub-group H ⊂ G, so that R(g)v = v for g ∈ H and R(g)v 6= v for g ∈ / H. Now introduce infinitesimal
transformations R(g) ≃ 1 + itI TI with generator TI (in the representation R) and small parameters tI . We can
split these generators into two sets, {TI } = {Hi , Sα }, where Hi are the generators of the unbroken sub-group
H and Sα are the remaining generators corresponding to the broken part of the group. Hence, these two types of
generators can be characterised by
Hi v = 0 , Sα v 6= 0 . (4.85)
Now, write down the infinitesimal version of Eq. (4.84)
T
I I ∂V
V (φ) = V (φ − it TI φ) = V (φ) − it (φ) TI φ , (4.86)
∂φ
∂V
differentiate one more time with respect to φ and evaluate the result at φ = v using that ∂φa (v) = 0. This leads
to
M TI v = 0 , (4.87)
where M is the mass matrix defined above. Every broken generator Sα satisfies Sα v 6= 0 and, hence, leads to an
eigenvector of the mass matrix with eigenvalue zero. In other words, every broken generator leads to one massless
scalar which is precisely the statement of Goldstone’s theorem.
Note that the invariance of this Lagrangian density under SU(2) is due to the appearance of a singlet in the Clebsch-
Gordan decomposition 2̄ ⊗ 2 = 1 + 3. Provided that m2 < 0, the scalar potential is minimised for
m2
φ† φ = v02 = − (4.91)
2λ
4.4. VECTOR FIELDS, GAUGE SYMMETRY AND SCALAR ELECTRODYNAMICS 53
τ 1 v 6= 0 , τ 2 v 6= 0 , (τ 3 − Y )v 6= 0 , (τ 3 + Y )v = 0 . (4.93)
Hence, three of the four generators of SU(2) × UY (1) are broken, while the generator τ 3 + Y remains unbroken.
This last generator corresponds to a combination of the U(1) ⊂ SU(2) and the additional UY (1) and defines the
unbroken U(1) subgroup. So the induced breaking pattern can be summarised as
This is precisely the symmetry breaking pattern which arises in the electro-weak sector of the standard model of
particle physics. There, SU(2) × UY (1) is the electro-weak (gauge) symmetry and the unbroken U(1) corresponds
to electromagnetism. In the present case we are working with a global symmetry and Goldstone’s theorem tells
us that we should have three massless scalars from the three broken generators. In the case of the electro-weak
theory, the SU(2) × UY (1) symmetry is actually promoted to a local (or gauge) symmetry where the symmetry
parameters are allowed to depend on space-time. In this case, it turns out that the Goldstone bosons are absorbed
by three vector bosons which receive masses from symmetry breaking. This phenomenon is also called the Higgs
effect and to investigate this in more detail we need to introduce vector fields and gauge symmetries.
Aµ → Aµ + ∂µ Λ (4.97)
on Aµ , parameterized by an arbitrary function 5 Λ = Λ(x), leaves the field strength tensor Fµν unchanged (as
can be easily seen by transforming the RHS of Eq. (4.95)). Let us now write down the most general Lagrangian
density for Aµ (up to second order in derivatives) which is Lorentz invariant and invariant under gauge transfor-
mations (4.97). Gauge invariance implies that the Lagrangian should depend on Aµ only through the field strength
F and, since F contains one derivative the most we should consider is quadratic terms in F . In addition, Lorentz
invariance means all indices should be contracted in L. Basically, this leaves only one allowed term 6 , namely
1
L = − Fµν F µν , (4.98)
4
where we think of Fµν as being given by Eq. (4.95). One finds
∂L 1 ∂(∂ρ Aσ − ∂σ Aρ ) 1 ∂L
= − F ρσ = − F ρσ δρµ δσν − δσµ δρν = −F µν ,
=0. (4.99)
∂(∂µ Aν ) 2 ∂(∂µ Aν ) 2 ∂Aµ
5 Notto be confused with a Lorentz transformation!
6 The term ǫµνρσ Fµν Fρσ is also consistent with all stated requirements. However, it can be written as the total derivative
4∂µ (ǫµνρσ Aν ∂ρ Aσ ) and, hence, does not effect the equations of motion.
54 CHAPTER 4. CLASSICAL FIELD THEORY
where the second equation is a trivial consequence of the definition (4.95). These are the free Maxwell’s equations
in covariant form. Splitting indices up into space and time components and inserting Eqs. (4.96) they can be easily
shown to be equivalent to the better-known version in terms of the electric and magnetic fields E and B. This
example illustrates the power of the Lagrangian formulation of field theories. Starting with a simple set of assump-
tions about the symmetries (Lorentz symmetry and gauge invariance in the present case) and the field content (a
single vector field Aµ ) one is led to the correct theory by writing down the most general Lagrangian consistent
with the symmetries.
2Aµ − ∂µ ∂ν Aν = 0 , (4.101)
which follows from the first Eq. (4.100) after inserting the definition (4.95) (Note that the second equation (4.100)
is automatically satisfied once Fµν is written in terms of Aµ ). This equation can be further simplified by exploiting
the gauge symmetry (4.97). Gauge invariance allows one, via transformation with an appropriate gauge parameter
Λ, to impose a gauge condition on Aµ . There are several possibilities for such gauge conditions and here we
consider the Lorentz gauge defined by
∂µ Aµ = 0 . (4.102)
This condition has the obvious benefit of being covariant (unlike, for example, the so-called temporal gauge which
requires A0 = 0) and it simplifies the equation of motion for Aµ to
2Aµ = 0 . (4.103)
Note, however, that the Lorentz gauge does not fix the gauge symmetry completely but leaves a residual gauge
freedom with gauge parameters Λ satisfying
2Λ = 0 . (4.104)
Eq. (4.103) is a massless Klein-Gordon equation for a vector field and, hence, can be easily solved using our earlier
result (4.56) with an additional µ index attached. This leads to
Z
Aµ (x) = d3 k̃ aµ (k)e−ikx + a⋆µ (k)eikx .
(4.105)
where wk = |k| and (kµ ) = (wk , k). In addition, the Lorentz gauge condition demands that the coefficents aµ
satisfy
k µ aµ (k) = 0 . (4.106)
(α)
To exploit this constraint in detail it is useful to introduce a set of polarisation vectors ǫµ (k), where α = 0, 1, 2, 3,
with the following properties. The vectors ǫ(1) (k) and ǫ(2) (k) are orthogonal to both k and a vector n with n2 = 1
and n0 > 0 in the time direction and they satisfy
′ ′
ǫ(α) (k) · ǫ(α ) (k) = −δ αα for α, α′ = 1, 2 . (4.107)
Further, ǫ(3) (k) is chosen to be in the (n, k) plane, orthogonal to n and normalised, that is, n · ǫ(3) (k) = 0 and
(ǫ(3) (k))2 = −1. Finally, we set ǫ(0) = n. With these conventions we have an orthogonal set of vectors satisfying
′ ′
ǫ(α) · ǫ(α ) = η αα (4.108)
where a(α) (k) are four expansion coefficients. The idea of introducing this basis of polarisation vectors is to
separate the directions transversal to k, corresponding to ǫ(1) (k) and ǫ(2) (k), from the other two directions ǫ(0) (k)
4.4. VECTOR FIELDS, GAUGE SYMMETRY AND SCALAR ELECTRODYNAMICS 55
and ǫ(3) (k). As an example, if we choose a spatial momentum k pointing in the z-direction, the above vectors are
explicitly given by
1 0 0 0
(0)
0 (1)
1 (2)
0 (3)
0
ǫ = 0 , ǫ = 0 , ǫ = 1 , ǫ = 0 .
(4.110)
0 0 0 1
Returning to the general case, it is easy to see that the unique choice for ǫ(3) (k) is
1 1
ǫ(3)
µ (k) = kµ − nµ so that ǫ(0) (3)
µ + ǫµ = kµ . (4.111)
k0 k0
Let us now return to the gauge condition (4.106) and insert the expansion (4.109). The two transverse directions
ǫ(1) (k) and ǫ(2) (k) drop out of the equation, so that
0 = k µ aµ (k) = (k · ǫ(0) (k))a(0) (k) + (k · ǫ(3) (k))a(3) (k) = k0 (a(0) (k) − a(3) (k)) . (4.112)
Hence, we are left with the two transversal polarisations and a longitudinal one along the direction of k. Recall
that we still have a residual gauge freedom from gauge parameters Λ satisfying the Klein-Gordon equation (4.104).
The most general such parameters can be written as
Z
Λ(x) = d3 k̃ λ(k)e−ikx + λ⋆ (k)eikx .
(4.115)
A gauge transformation (4.97) with Λ of this residual form changes the expansion coefficients aµ (k) for the vector
field as
aµ (k) → a′µ (k) = aµ (k) − ikµ λ(k) . (4.116)
Comparing with Eq. (4.114), it is clear that we can use this residual gauge freedom to remove the longitudinal
degree of freedom in aµ (k). We are then left with the two transversal polarisations only and we conclude that
the number of physical degrees of freedom for a vector field is two. This reduction from four apparent degrees of
freedom to two is directly related to the gauge invariance of the theory.
∂µ F µν + m2 Aν = 0 . (4.117)
Applying ∂ν to this equation we conclude that ∂ν Aν = 0 and, hence, that Eq. (4.117) can equivalently be written
as
(2 + m2 )Aµ = 0 , ∂µ Aµ = 0 . (4.118)
The first of these equations is a massive Klein-Gordon equation with general solution given by Eq. (4.105) but
now with k 2 = m2 instead of k 2 = 0. To also satisfy the second equation above we need k µ aµ (k) = 0 which
reduces the number of degrees from four to three. Since gauge invariance is not available in this case, no further
reduction occurs and we conclude that the number of physical degrees of freedom of a massive vector field is three.
Including a source
Let us return to the massless theory and ask how can we include a source Jµ into this Lagrangian formulation of
56 CHAPTER 4. CLASSICAL FIELD THEORY
Maxwells’ theory. This current should appear linearly in the equations of motion (and therefore in the Lagrangian
density), so the obvious generalisation of Eq. (4.98) is
1
L = − Fµν F µν − Aµ J µ . (4.119)
4
The additional term depends on Aµ explicitly which leads to a non-trivial gauge variation
Z Z
S → S − d4 x ∂µ ΛJ µ = S + d4 x Λ∂µ J µ , (4.120)
of the action. This apparent breaking of the gauge invariance of the theory can be avoided if we require that
∂µ J µ = 0 , (4.121)
∂L
that is, if the current Jµ is conserved. With ∂A µ
= −J µ and the first Eq. (4.99) we find from the Euler-Lagrange
equations
∂µ F µν = J ν , ∂[µ Fνρ] = 0 , (4.122)
which are indeed Maxwell’s equations in the presence of a source Jµ . In a fundamental theory, the current Jµ
should arise from fields, rather than being put in ”by hand” as an external source. We will now study an example
for such a theory in which vector fields are coupled to scalars.
The additional term proportional to ∂µ Λ in this transformation is reminiscent of the transformation law (4.97) for
the gauge field and Aµ can indeed be used to cancel this term. Define the covariant derivative
Dµ = ∂µ + iqAµ . (4.125)
Dµ φ → e−iqΛ(x) Dµ φ , (4.126)
and the modified kinetic term (Dµ φ)⋆ Dµ φ is gauge invariant. With this modification, we can now combine
Eqs. (4.72) and (4.98) to obtain the gauge invariant Lagrangian density
1 λ ⋆ 2
L = (Dµ φ)⋆ Dµ φ − V (φ, φ⋆ ) − Fµν F µν , V = V0 + m2 φ⋆ φ + (φ φ) (4.127)
4 4
This is the Lagrangian for scalar electrodynamics. It shows the gauge field Aµ in a new role, facilitating the
invariance of the scalar field theory under local symmetry transformations. In fact, had we started with just the
globally symmetric scalar field theory (4.72) with the task of finding a locally U(1) invariant version we would
4.5. FURTHER READING 57
have been led to introducing a gauge field Aµ . We also note that the covariant derivative in (4.127) has introduced
a non-trivial coupling between the scalar field and the vector field.
We will now use scalar electrodynamics as a toy model to study spontaneous breaking of a local U(1) symme-
try. The scalar potential is unchanged from the globally symmetric model and, hence, we can apply the results of
Sec. 4.3.2 . For m2 ≥ 0 the minimum is at φ = 0 and the symmetry is unbroken. Therefore, we focus on the case
m2 < 0 where the potential has a ”Mexican hat” shape as in Fig. 4.2. In this case, there is a circle of minima given
by Eq. (4.78). Instead of choosing the parameterisation (4.80) for φ around the minimum on the φ1 axis it proves
more useful in the present context 7 to use
1
φ = √ (v0 + H)eiχ (4.128)
2
where χ = χ(x) is the Goldstone mode and H = H(x) is the massive mode. Next we perform a gauge transfor-
mation with parameter Λ = χ/q, that is,
1 1
φ → φ′ = e−iχ φ = √ (v0 + H) , Aµ → A′µ = Aµ + ∂µ χ . (4.129)
2 q
We can now write the Lagrangian density (4.127) in terms of the transformed fields (this leaves the action un-
changed) and then insert the explicit expression for φ′ from the previous equation. One finds for the covariant
derivative
1
Dµ φ′ = √ (∂µ H + iq(v0 + H)A′µ ) . (4.130)
2
and inserting into Eq. (4.127) then leads to
1 1 ′ ′µν 1 2 ′ ′µ 2
L = ∂µ H∂ µ H − V (H) − Fµν F + q Aµ A (v0 + 2v0 H + H 2 ) (4.131)
2 4 2
1 2 2 λ
V (H) = V0 + m v0 − m H + (4v0 H + H 4 ) .
2 2 3
(4.132)
4 16
The most striking feature about this result is that the Goldstone mode χ has completely disappeared from the
Lagrangian and we are left with just the real, massive scalar H and the vector field A′µ . However, the vector field
is now massive with mass
m(A′ ) = qv0 , (4.133)
end, hence, has three degrees of freedom as opposed to just two for a massless vector field. This explains the
disappearance of the scalar χ: It has been ”absorbed” by A′µ to provide the additional (longitudinal) degree of
freedom necessary for a massive vector field. This can also be seen from the transformation of Aµ in Eq. (4.129).
Hence, we have learned that a spontaneously broken local symmetry leads to a mass for the associated vector
boson and the conversion of the Goldstone field into the longitudinal mode of the vector. This is also called the
Higgs effect. The same phenomenon but in its generalisation to non-Abelian gauge groups occurs in the breaking
of the electro-weak SU(2) × UY (1) → U(1) gauge group to U(1) of electromagnetism. The three W ± and Z
vector bosons of the broken part of the electro-weak symmetry receive a mass (4.133) proportional to the symmetry
breaking scale v0 and absorb the three Goldstone bosons (see Section 4.3.4) which arise. A detailed discussion of
electro-weak symmetry breaking requires introducing non-Abelian gauge symmetries which is beyond the scope
of this lecture.
As usual, be warned that these references were not written with undergraduate readers in mind.
7 One can think of Eq. (4.80) as a linearized version of the more accurate parameterization (4.128).
58 CHAPTER 4. CLASSICAL FIELD THEORY
Chapter 5
Canonical Quantization
[q̂i (t), q̂j (t)] = [p̂i (t), p̂j (t)] = 0 , [p̂i (t), q̂j (t)] = −iδij . (5.1)
The dymanics of this system is governed by the operator version of Hamilton’s equations
[φ̂ia (t), φ̂bj (t)] = [π̂ai (t), π̂bj (t)] = 0 , [π̂ai (t), φ̂bj (t)] = −iδab δij . (5.3)
The continuum version of these equations is obtained by replacing i → x, j → y, δij → δ 3 (x − y), φ̂ai (t) →
φ̂(t, x) and similarly for the conjugate momenta. This results in
[φ̂a (t, x), φ̂b (t, y)] = [π̂ a (t, x), π̂ b (t, y)] = 0 , [π̂ a (t, x), φ̂b (t, y)] = −iδba δ 3 (x − y) (5.4)
Note that these commutators are taken at equal time but at generally different points in space. The canonical
commutation relations (5.4) together with the continuum version
∂0 φ̂a (t, x) = i[Ĥ, φ̂a (t, x)] , ∂0 π̂ a (t, x) = i[Ĥ, π̂ a (t, x)] , (5.5)
of Hamilton’s equations (5.2) provide the starting point for the canonical quantization of field theories. We will
now investigate the consequences of this quantization procedure for the simplest case, a single real, free scalar
field. To avoid cluttering the notation, we will drop hats on operators from now on. It will usually be clear from
the context whether we refer to a classical object or its operator version.
1 We are using operators in the Heisenberg picture.
2 More formally, the transition between classical and quantum mechanics can also be understood as a replacement of canonical coordinates
with operators and Poisson brackets with commutator brackets.
59
60 CHAPTER 5. CANONICAL QUANTIZATION
2φ + m2 φ = 0 . (5.7)
It should now be thought of as an operator equation for the field operator φ. The canonical commutation rela-
tions (5.4) in this case read
[φ(t, x), φ(t, y)] = [π(t, x), π(t, y)] = 0 , [π(t, x), φ(t, y)] = −iδ 3 (x − y) . (5.8)
What do the canonical commutation relations for φ and π imply for the commutators of a(k) and a† (k)? To answer
this question, we Rwould like to express a(k) in terms of φ(x) and π(x) by inverting the above relations. We start
by applying 2wq d3 x e−iq·x to the equation for φ(x). After carrying out the integration over k 3 one finds
Z
2wq d3 x e−iq·x φ(x) = a(q)e−iwq t + a† (−q)eiwq t . (5.10)
Z
2i d3 x e−iq·x π(x) = a(q)e−iwq t − a† (−q)eiwq t (5.11)
Adding the last two equations we find the desired expression for a(q) and its conjugate
Z Z
a(q) = d x e (wq φ(x) + iπ(x)) , a (q) = d3 x e−iqx (wq φ(x) − iπ(x)) .
3 iqx †
(5.12)
Combining these results and the canonical commutation relations (5.8) one finds for the commutators of a(q) and
a† (q)
[a(k), a(q)] = [a† (k), a† (q)] = 0 , [a(k), a† (q)] = (2π)3 2wk δ 3 (k − q) . (5.13)
For each k = q these equations are reminiscent of the commutation relations for the creation and annihila-
tion operators of a harmonic oscillator. This should not come as a surprise given a single plain wave Ansatz
φ(t, x) = φk (t)eik·x turns the Hamiltonian (5.6) (or, more precisely, as the Ansatz is complex, its complex
2 2 2
√ (4.75)) into the Hamiltonian Hk = πk + wk φk for a single harmonic oscillator with frequency
counterpart
wk = m2 + k2 . We should, therefore, think of the free scalar field as an infinite collection of (decoupled)
harmonic oscillators labelled by three-momentum k and with frequency wk .
Single particle states |ki are obtained by acting on this vacuum state with a single creation operator, so that
Using the commutation relations (5.13) and Eq. (5.14), one finds for their normalization
We note that the RHS of this relation is the ”covariant delta function” (4.55). A basis set of states is provided by
all n particle states |k1 , . . . , kn i, obtained by acting on the the vacuum with n creation operators
We can repeat this procedure and commute N with all creation operators, picking up at each step the term
|k1 , . . . , kn i. In the last step we use N |0i = 0 and find
form Hk = wk (a† (k)a(k) + 21 ) using the commutator (5.13) one encounters an infinite zero point energy, propor-
tional to δ 3 (0), which can be interpreted as the energy of the vacuum state |0i. This is one of the many infinities
in quantum field theory and it should not come as a surprise: After all, we are summing over an infinite number
of harmonic oscillators each with finite zero-point energy wk /2. To deal with this infinity, we define the concept
of normal ordering of operators. The normal ordered version : O : of an operator O is obtained by writing all
creation operators to the left of all annihilation operators. So, for example, : a(k)a† (q) : = a† (q)a(k). An operator
and its normal-ordered counterpart differ by a (usually infinite) number, often referred to as a c-number. With this
definition we have
Z Z Z
: H : = d k̃ wk a (k)a(k) , : P : = d k̃ k a (k)a(k) or : Pµ : = d3 k̃ kµ a† (k)a(k) .
3 † 3 †
(5.29)
and, hence, these states are eigenstates of : Pµ : with energy P wk and spatial momentum k. Likewise, multi-particle
n
states |k1 , . . . , kn i are eigenstates of : Pµ : with eigenvalue i=1 kiµ . To summarise, we have constructed a basis
of Fock space for the free real scalar field and, by looking at the action of conserved operators on these states, we
have found an interpretation of this basis as n particle states with definite four-momentum.
d4 k
Z
′
= {δ(k0 − wk ) + δ(k0 + wk )} ǫ(k0 )e−ik(x−x )
(2π)3 2wk
d4 k
Z
′
= δ(k 2 − m2 )ǫ(k0 )e−ik(x−x ) ≡ ∆(x − x′ ) , (5.32)
(2π)3
where ǫ(k0 ) = k0 /|k0 |. All elements in the final integral are Lorentz invariant except for the function ǫ(k0 ).
However, if we restrict to Lorentz transformations Λ which preserve the sign of the time-component of four vectors
(these are the orthochronos Lorentz transformations, see Table 4.1) then ǫ(k0 ) remains unchanged and we have
∆(x) = ∆(Λx). We already know that ∆(0, x) = 0 since this corresponds to an equal-time commutator. Now
consider a space-like vector x, that is x2 < 0. With a suitable orthochronos Lorentz transformation Λ this vector
can be written in the form x = Λ(0, y) for some vector y and, hence, ∆(x) = ∆(Λ(0, y)) = ∆(0, y) = 0. We
conclude that ∆(x) vanishes for all space-like vectors x. This means that the commutator [φ(x), φ(x′ )] vanishes
whenever x and x′ are space-like separated. This fact is also referred to as micro causality. Two field operators at
points with space-like separation should not be able to causally effect one another and one would, therefore, expect
their commutator vanishes. We have just shown that this is indeed the case.
This concludes the canonical quantization of the free real scalar field and we move on to the next, more
complicated case, the free complex scalar field.
5.3. CANONICAL QUANTIZATION OF THE FREE COMPLEX SCALAR FIELD 63
2φ + m2 φ = 0 . (5.34)
Canonical quantisation of this system can be performed by simply imposing the general commutation relations (5.4)
on the two real scalars φ1 and φ2 and their associated conjugate momenta π1 = ∂0 φ1 and π2 = ∂0 φ2 . The only
non-zero commutation relations are then
[π1 (t, x), φ1 (t, y)] = −iδ 3 (x − y) , [π2 (t, x), φ2 (t, y)] = −iδ 3 (x − y) . (5.35)
√ √
Using φ = (φ1 + iφ2 )/ 2 and π = (π1 + iπ2 )/ 2 they can be easily translated to the complex field and its
conjugate with the only non-zero commutation relations
[π(t, x), φ(t, y)] = [π † (t, x), φ† (t, y)] = −iδ(x − y) . (5.36)
As for the real scalar field below, we can invert these relations and compute the commutators of a± (k) and a†± (k)
from the canonical commutation relations (5.36). The only non-zero commutators one finds in this way are
This shows that we have two sets, {a†+ (k), a+ (k)} and {a†− (k), a− (k)}, of creation and annihilation operators
and, hence, two different types of one-particle states
As usual, the vacuum |0i is defined by a± (k)|0i = 0 and h0|0i = 1. Multi-particle states |(k1 , ǫ1 ), . . . , (kn , ǫn )i =
a†ǫ1 (k1 ) . . . a†ǫn (kn )|0i are now labelled by n momenta ki and, in addition, n signs ǫi ∈ {+1, −1} to distinguish
the two types of quanta.
Hence N+ (N− ) acting on a multi-particle state counts the number of quanta of + type (− type). The (normal-
ordered) conserved four-momentum can be computed as for the real scalar field and one finds
Z
: Pµ : = d3 k̃ kµ a†+ (k)a+ (k) + a†− (k)a− (k) , (5.42)
64 CHAPTER 5. CANONICAL QUANTIZATION
In particular, this shows that the states |(k, +)i have charge +q and the states |(k, −)i have charge −q. Therefore,
it is sensible to identify + states with particles and − states with anti-particles. With this terminology, we see from
Eq. (5.37) that φ(x) annihilates a particle or creates an anti-particle while φ† (x) creates a particle or annihilates an
anti-particle.
T (φ(t1 , x1 ) . . . φ(tn , xn )) = φ(tin , xin ) . . . φ(ti1 , xi1 ) where ti1 ≤ ti2 ≤ · · · ≤ tin . (5.46)
and their vacuum expectation values. The above equation means time ordering T rearranges fields operators so
that time increases from the right to the left.
where the Heaviside function θ is defined by θ(x0 ) = 1 for x0 ≥ 0 and θ(x0 ) = 0 for x0 < 0. To compute the
Feynman propagator we first evaluate the above expression for the case x0 > y0 . Inserting the field expansion (5.9)
we have
∆F (x − y) = h0|φ(x)φ(y)|0i (5.48)
Z
h0| d3 k̃ d3 q̃ a(k)e−ikx + a† (k)eikx a(q)e−iqy + a† (q)eiqy |0i
= (5.49)
Z Z
= h0| d3 k̃ d3 q̃ [a(k), a† (q)]eiqy−ikx |0i = d3 k̃ e−ik(x−y) . (5.50)
d3 k
Z n o
−iwk (x0 −y0 ) iwk (x0 −y0 )
∆F (x − y) = θ(x0 − y 0 )e + θ(y 0 − x0 )e eik·(x−y) (5.51)
(2π)3 2wk
d4 k i
Z
= 4
e−ik(x−y) (5.52)
(2π) (k0 − wk + iǫ̃)(k0 + wk − iǫ̃)
d4 k i
Z
= e−ik(x−y) . (5.53)
(2π) k − m2 + iǫ
4 2
The small quantity ǫ̃ > 0 in the second integral is to indicate that the pole at k0 = wk − iǫ̃ is slightly below the
real k0 axis and the pole at k0 = −wk + iǫ̃ is slightly above (The quantity ǫ > 0 in the final integral serves the
same purpose). With this understanding about the position of the poles the equality of the first and second line
above can be shown by a contour integration as indicated in Fig. 5.1. The above result provides us with a simple
x0 < y0
-wk + i ε~
wk - i ε~ k0
x0 > y0
Figure 5.1: Location of poles for the Feynman propagator and contours to prove the equality between Eqs. (5.51)
and (5.52). For x0 > y0 the integration along the real k0 axis can be closed for Im(k0 ) < 0 since the real part
of the exponent in Eq. (5.52) is negative in this case. Only the pole at k0 = wk − iǫ̃ constributes and leads to the
first term in Eq. (5.51). Analogously, for x0 < y0 the contour can be closed for Im(k0 ) > 0. Only the pole at
k0 = −wk + iǫ̃ contributes and leads to the second term in Eq. (5.51)
˜ F is given by
representation of the Feynman propagator and tells us that its Fourier transform ∆
˜ F (k) = i
∆ . (5.54)
k 2 − m2 + iǫ
The Feynman propagator, particularly in the momentum space form (5.54), is central in the formulation of Feyn-
man rules for the perturbation theory of interacting fields as we will see later.
with a delta-function source. Functions with this property are also called Green functions. There are other Green
functions of the Klein-Gordon equation which are given by an integral such as in Eq. (5.53) but with the poles in
a different position relative to the real k0 axis. For example, the case where both poles are below the real k0 axis
leads to the so-called retarded Green function. It vanishes for x0 < y0 since the upper contour in Fig. 5.1 contains
no poles in this case.
66 CHAPTER 5. CANONICAL QUANTIZATION
∂2 ∂
∆F (x − y) = {δ(x0 − y0 )h0|[φ(x), φ(y)]|0i + h0|T (π(x), φ(y))|0i}
∂x20 ∂x0
= δ(x0 − y0 )h0|[π(x), φ(y)]|0i + h0|T (∂02 φ(x), φ(y))|0i
= −iδ 4 (x − y) + h0|T (∂02 φ(x), φ(y))|0i ,
∂
where we used that ∂x 0
θ(x0 ) = δ(x0 ) and δ(x0 − y0 )h0|[φ(x), φ(y)]|0i = δ(x0 − y0 )∆(0, x − y) = 0 (see the
discussion at the end of Section 5.2). It follows that
2
∂2
2 ∂ 2
(2x + m )∆F (x − y) = − + m ∆F (x − y) (5.56)
∂x20 ∂x2
= −iδ 4 (x − y) + h0|T (2x + m2 )φ(x)φ(y) |0i .
(5.57)
The field φ(x) satisfies the (free) Klein-Gordon equation and, hence, the second term vanishes. This completes the
proof.
Since φ+ (x) only contains annihilation operators it is clear that φ+ (x)|0i = 0. Likewise, h0|φ− (x) = 0. We begin
with the time ordered product of two fields φ(x) and φ(y) for the case x0 > y0 and write
The point about introducing the commutator in the second line is that the first four terms have the creation operators
to the left of the annihilation operators and, hence, correspond to the normal ordering of the field product. A similar
calculation for the case x0 < y0 leads to T (φ(x)φ(y)) =: φ(x)φ(y) : +h0|φ(y)φ(x)|0i. Combining these two
results we have
T (φ(x)φ(y)) = : φ(x)φ(y) + ∆F (x − y) : . (5.59)
Hence, we see that time and normal ordering are related via a Feynman propagator. When two field operators in a
time-ordered product are combined into a Feynman propagator as in the above equation they are said to have been
contracted. With this terminology we can formulate Wick’s theorem as
By ”all possible contraction” we refer to all possible ways of pairing up the field operators φ(x1 ), . . . , φ(xn )
into Feynman propagators, including partial pairings. When taking the vacuum expectation value of the above
equation the first term and all partially contracted terms vanish due to normal ordering. This shows that the
vacuum expectation of an odd number of time-ordered fields vanishes. For an even number of fields we are left
with X
h0|T (φ(x1 ) . . . φ(xn )) |0i = ∆F (xi1 − xi2 ) . . . ∆F (xin−1 − xin ) , (5.61)
pairings p
where the sum runs over all pairings p = {(i1 , i2 ), . . . , (in−1 , in )} of the numbers 1, . . . , n. This is precisely the
same structure as the one we encountered in the context of n-point functions for Gaussian multiple integrals and
Gaussian functional integrals. This is of course not an accident and the precise relation will become clear when we
discuss the path integral quantisation of field theories.
done this in Eq. (5.59). The general case can be proven by induction in n. Rather than presenting this general
argument it is probably more instructive to consider the example for n = 3. After a suitable relabelling we can
arrange that t1 ≥ t2 ≥ t3 . With the short-hand notation φi = φ(xi ) and ∆ij = ∆F (xi − xj ) we have
T (φ1 φ2 φ3 ) = φ1 φ2 φ3 = φ1 : φ2 φ3 + ∆23 :
= (φ1+ + φ1− ) : φ2+ φ3+ + φ2− φ3+ + φ3− φ2+ + φ2− φ3− + ∆23 :
= : φ1− φ2 φ3 + φ1+ φ2+ φ3+ + φ2− φ1+ φ3+ + φ3− φ1+ φ2+ + φ2− φ3− φ1+ + ∆23 φ1 :
+[φ1+ , φ2− ]φ3+ + [φ1+ , φ3− ]φ2+ + [φ1+ , φ2− φ3− ]
= : φ1 φ2 φ3 + ∆12 φ3+ + ∆13 φ2+ + ∆12 φ3− + ∆13 φ2− + ∆23 φ1 :
= : φ1 φ2 φ3 + ∆12 φ3 + ∆13 φ2 + ∆23 φ1 :
In the first line we have simply used Wick’s theorem (5.59) for a product of two fields, applied to φ2 and φ3 . In
the second line, we need to move φ1+ and φ1− into the normal ordering. This is easy for φ1− since it consists of
creation operators and, hence, has to be on the left of a product anyway. The annihilation part φ1+ of the field, on
the other hand, has to be commuted to the right of any component φi− . The key is that the commutators which
arise in this way precisely lead to the necessary contractions of φ1 with φ2 and φ3 . The proof for general n is
analogous to the above calculation. It applies the induction assumption (that is, the validity of Wick’s theorem
for n − 1 fields) to φ2 , . . . , φn . Then, moving φ1 into the normal ordering and commuting φ1+ with all negative
frequency parts φi− generates all the missing contractions of φ1 with the other fields.
h0|T φ(x)φ† (y) |0i = ∆F (x − y) , h0|T (φ(x)φ(y)) |0i = h0|T φ† (x)φ† (y) |0i = 0 .
(5.62)
For a product of operators φ and φ† Wick’s theorem can be applied straightforwardly, but the above equations tell
us that only contractions of φ with φ† need to be taken into account.
Aµ → Aµ + ∂µ Λ . (5.64)
However, Eq. (5.65) implies that the canonical momentum π 0 vanishes and this is inconsistent with the µ = ν = 0
part of Eq. (5.66). Clearly, viewing Aµ as fourP scalar fields is too simple. In fact, while typical kinetic terms
for four scalar fields Aµ would be of the form ν ∂µ Aν ∂ µ Aν and, hence, depend on the symmetric and anti-
symmetric parts of ∂µ Aν , the Lagrangian (5.63) only depends on the anti-symmetric part, that is, on the field
68 CHAPTER 5. CANONICAL QUANTIZATION
strength Fµν . This special form of the Maxwell Lagrangian is of course responsible for the existence of the gauge
symmetry (5.64) as well as for the vanishing of the conjugate momentum π 0 and it links these two features. In
essence, gauge symmetry is the crucial difference to the scalar field theory. There are various viable methods to
quantise a gauge theory, but here we will follow the most obvious approach of fixing a gauge before quantisation.
Gauge fixing
Since we would like to manifestly preserve covariance we use the Lorentz gauge condition (4.102), which we im-
pose on the theory by means of a Lagrange multiplier λ. This means that, instead of the Maxwell Lagrangian (5.63),
we start with
1 λ
L = − Fµν F µν − (∂µ Aµ )2 , (5.67)
4 2
For the conjugate momenta we now find
∂L
πµ = = F µ0 − λη µ0 ∂ν Aν , (5.68)
∂(∂0 Aµ )
and, hence, in particular π 0 = −λ∂ν Aν is no longer zero. The equation of motion (4.101) for Aµ is now modified
to
2Aµ − (1 − λ)∂µ ∂ν Aν = 0 , (5.69)
and, in addition, we have to impose the gauge condition ∂µ Aµ = 0 (which formally arises from the action (5.67)
as the λ equation of motion). However, clearly we should not impose this condition as an operator equation since
this would lead us back to a situation where π 0 = 0. Instead, it will later be imposed as a condition on physical
states. At any rate, the obstruction to imposing canonical quantisation conditions has been removed and we require
the canonical commutation relations (5.66) for Aµ and the conjugate momenta (5.68).
Oscillator expansion
We should now solve the equation of motion (5.69) to work out the properties of the creation and annihilation
operators. While this can be done for arbitrary λ, we adopt the so-called Feynman gauge 4 , λ = 1, which simplifies
the equation of motion (5.67) to 2Aµ = 0. The general solution to this equation has already been written down in
Section 4.4.1 and is given by
3 Z
X
Aµ (x) = d3 k̃ ǫ(α)
µ (k) a (α)
(k)e −ikx
+ a (α)†
(k)e ikx
. (5.70)
α=0
(α)
Here, we have used the polarisation vectors ǫµ (k) defined before Eq. (4.109) and the four-momentum is given
(α)
by (kµ ) = (wk , k) with wk = |k|. We recall from our discussion in Section 4.4.1 that classically ǫµ (k) for
α = 1, 2 are the two transversal, physical polarisations while the other two polarisations can be removed through
gauge transformations.
Inserting Eq. (5.70) into Eq. (5.68) gives the expansion for the conjugate momenta. Similar to the scalar field
case (see below Eq. (5.9)), we can now invert these expansions to express the oscillators a(α) (k) and their hermitian
conjugates in terms of the field Aµ and the momenta π µ . As for scalar fields, these results can be inserted into
the quantisation condition (5.66) to determine the commutation relations of the oscillators. After a straightforward
(but slightly lengthy) calculation one finds that
′ ′ ′ ′
[a(α) (k), a(α ) (q)] = [a(α)† (k), a(α )† (q)] = 0 , [a(α) (k), a(α )† (q)] = −η αα (2π)3 2wk δ 3 (k − q) (5.71)
remark that independence of the quantisation procedure on the choice of λ can be explicitly demonstrated.
5.5. CANONICAL QUANTIZATION OF A VECTOR FIELD 69
and, hence, has negative norm. Clearly this is physically unacceptable and indicates that the space F0 still contains
unphysical states and cannot be the proper Fock space of the theory. A related problem emerges when looking at
the conserved four-momentum. Calculating the four-momentum for the theory (5.67) from the general formalism
in Section 4.2.3 and inserting the field expansion (5.70) one finds after some calculation
Z " 3 #
X
3 (α)† (α) (0)† (0)
: Pµ : = d k̃ kµ a (k)a (k) − a (k)a (k) . (5.73)
α=1
The negative sign in front of the last term means that the state |(k, 0)i has ”negative energy”.
The existence of unphysical states is not at all surprising given that we have not yet imposed the Lorentz gauge
condition. We have seen above that requiring the operator equation ∂µ Aµ = 0 is too strong and leads to problems
with quantisation. Instead we define a space F1 ⊂ F0 of physical states on which the gauge condition is satisfied,
that is we require
hΦ̃|∂µ Aµ |Φi = 0 . (5.74)
between two physical states |Φi, |Φ̃i ∈ F1 . To guarantee this condition is satsified it is sufficient that (∂µ Aµ )(+) |Φi =
0 where the annihilation part (∂µ Aµ )(+) of ∂µ Aµ is proportional to
Z
d3 k̃ k0 a(0) (k) − a(3) (k) e−ikx . (5.75)
(α)
To obtain this last result, we have taken the ∂µ derivative of the first term in Eq. (5.70) and used that k µ ǫµ (k) = 0
(0) (3)
for α = 1, 2 and k µ ǫµ (k) = −k µ ǫµ (k) = k0 (see Section 4.4.1). This means physical states |Φi ∈ F1 are
defined by the condition
b− (k)|Φi = 0 . (5.76)
(3) (0)
√
where we have defined a new basis b± (k) = (a (k) ± a (k))/ 2 for the non-transversal operators. Clearly
all transversal states, that is, states created from the vacuum by acting only with transversal creation operators
a(α)† (k), where α = 1, 2 along with their linear combinations, satisfy the condition (5.76) and, hence, are elements
of F1 . It would be nice if F1 consisted of such states only but things are not quite so simple. To analyse the
condition (5.76) for non-transversal states we note that
while [b− (k), b†+ (q)] 6= 0. This means, in addition to a(α)† (k), where α = 1, 2, physical states can also contain
b†− (k) but not b†+ (k). This means F1 is spanned by states created from the vacuum acting with any number of
operators a(α)† (k), where α = 1, 2 and b†− (k). Note that the commutators of these operators with their associated
annihilation operators are all non-negative and, hence, no states with negative norm are left in F1 . However, if a
state contains at least one operator b†− (k) its norm vanishes, as is clear from the vanishing of the commutator (5.77).
Physically, we should discard such zero norm states and the formal way of doing this is to identify each two states
in F1 if their difference has zero norm. In this way, we obtain the proper Fock space F2 , whose elements are the
classes of states obtained from this identification. In particular, in each class there is a ”representative” with only
transverse oscillators. In conclusion, we see that the proper Fock space F2 can be thought of as spanned by states
of the form |(k1 , α1 ), . . . , (kn , αn )i = a(α1 )† (k1 ) . . . a(αn )† (kn )|0i, where αi = 1, 2 are transverse oscillators
only.
Taken between a physical state |Φi ∈ F1 , that is hΦ| : Pµ : |Φi, the last two terms in Eq. (5.78) do not contribute
(using that all a type operators commutate with all b type operators as well as Eqs. (5.76), (5.77)). This means
Z 2
X
hΦ| : Pµ : |Φi = hΦ| d3 k̃ kµ a(α)† (k)a(α) (k)|Φi . (5.79)
α=1
70 CHAPTER 5. CANONICAL QUANTIZATION
This shows that the four-momentum only depends on the transverse modes and, since every class in F2 has exactly
one representative with transverse modes only, the desired independence on the choice of representative.
Feynman propagator
Finally, we should look at the Feynman propagator for a vector field, defined, as usual, as the vacuum expectation
value h0|T (Aµ (x)Aν (y))|0i of a time-ordered product of two fields. Inserting the field expansion (5.70) and using
the commutation relations (5.71) we can perform a calculation completely analogous to the one in Section 5.4 to
obtain the vector field propagator in Feynman gauge, λ = 1. This leads to
In this chapter, we will develop the formalism for perturbatively interacting quantum field theory. The main goal
is to derive the Feynman rules for bosonic field theories from first principles and show how to calculate matrix
elements and cross sections. To do this we will first introduce the S-matrix and then show how cross sections
and decay rates are expressed in terms of S-matrix elements. We then need to understand how to compute S-
matrix elements from field theory. The first step is to derive the reduction formula which shows how to compute
S-matrix elements from vacuum expectation values of time-ordered products of interacting fields. These products
of interacting fields are then re-written in terms of free fields using the evolution operator. From there, Wick’s
theorem will lead to the Feynman rules. Whenever the formalism needs to be developed with reference to a
particular theory we will first focus on the real scalar field theory (4.47) with λφ4 interaction for simplicity and
subsequently generalise to more complicated theories.
|k1 , . . . , kn iin = a†in (k1 ) . . . a†in (kn )|0i , |k1 , . . . , kn iout = a†out (k1 ) . . . a†out (kn )|0i (6.1)
Given an initial state |iiin = |k1 , . . . , kn iin with n particles and momenta k1 , . . . , kn and a final state |f iout =
|q1 , . . . , qm i with m particles and momenta q1 , . . . , qm , we are interested in computing the amplitude
which provides the probability for a transition from |iiin to |f iout . With the operator S defined by
and, hence, the S matrix is unitary. We also require that the vacuum state is invariant and unique so that
71
72 CHAPTER 6. INTERACTING QUANTUM FIELDS
where ki and qf are the total four-momenta of the initial and final state and M is the invariant matrix element. We
need to understand how to compute M from the underlying field theory and, to make contact with experimentally
measurable quantities, we need to express cross sections and decay rates in terms of M. We start with the latter.
where p0 = wp in the second integral, so that fa (x) are solutions of the Klein-Gordon equation. By writing
fa (x) = e−ika x Fa (x) with slowly varying functions Fa (x) we ensure that these solutions are ”close” to momen-
tum eigenstates with momenta ka . The probability currents jaµ for these solutions are given by the usual formula
in relativistic quantum mechanics
jaµ ≡ i (fa⋆ (x)∂µ fa (x) − fa (x)∂µ fa⋆ (x)) ≃ 2kaµ |fa (x)|2 . (6.10)
They provide us with expressions for the number of particles per volume and the flux per volume which will be
needed in the calculation of the cross section.
In a first instance, we are interested in the transition probability
from the initial state |ı̃iin as defined above to some final state |f iout . Inserting Eqs. (6.7), (6.8) and (6.9) into this
expression 2 and keeping only the non-trivial scattering part, related to the T matrix, one finds
n
Z Y ! n ! !
d p̃a f˜ (pa ) d p̃ f˜b (p ) δ
Y X X
8 3 ⋆ 3 ′ ′ 4 ′
W (ı̃ → f ) = (2π) a bpa − b p (6.12) b
a=1 b=1 a b
!
X
4
× δ qf − pa M(p1 , . . . , pn → f )⋆ M(p′1 , . . . , p′n → f ) . (6.13)
a
The wave functions f˜a are peaked around momenta ka and we can, hence, approximate M(p1 , . . . , pn → f ) ≃
M(p′1 , . . . , p′n → f ) ≃ M(k1 , . . . , kn → f ) in the above integral. Together with the integral representation
! Z !!
X X X X
4 4 ′ 4 ′
(2π) δ pa − pa = d x exp −ix pa − pa (6.14)
a a a a
2 When doing this, we have to keep in mind that the ”in” states |ii in the definition (6.8) of the matrix element are the exact momentum
in
eigenstates |k1 , . . . , k1 iin , whereas the ”in” states |ı̃iin which enter the transition probability (6.11) are the wave packets (6.9).
6.2. CROSS SECTIONS AND DECAY RATES 73
We have obtained our first important result, the transition probability from state |ı̃iin to state |f iout per unit time
and volume
n
! !
dW (ı̃ → f ) Y
2 4 4
X
= |fa (x)| (2π) δ qf − ka |M(k1 , . . . , kn → f )|2 . (6.16)
dV dt a=1 a
Decay rate
Let us now be more specific and first consider a single particle with mass M , momentum k = (M, 0) and wave
function f (x) in the initial state, that is, we study the decay of a particle at rest. Let us define the decay rate
Γ(k → f ) as the transition probability per unit time and volume (6.16) divided by the number density of initial
particles and integrated over the momenta q1 , . . . , qm of the final state |f iout with energy resolution ∆. From
Eq. (6.10), the number density of initial particles is given by 2M |f (x)|2 and, hence, from Eq. (6.16) we find for
the decay rate into m particles
m
!
1
Z Y
Γ(k → f ) = d q̃b (2π)4 δ 4 (k − qf ) |M(k → f )|2 .
3
(6.17)
2M ∆
b=1
Here, we have re-written the original kinematical pre-factor 1/(4m2 |k1 |) in a covariant way using the identity
q q
m2 |k1 | = m2 k10 2 − m2 = (k1 · k2 )2 − m21 m22 . (6.19)
1
In this way, the above expression for the cross section is manifestly covariant.
where E = k10 + k20 = q10 + q20 is the total center of mass energy. We write the integral in Eq. (6.18) as
d3 q1 d3 q2
Z Z
d3 q̃1 d3 q̃2 (2π)4 δ 4 (k1 + k2 − q1 − q2 ) = 3 3
(2π)4 δ(E − q10 − q20 )δ 3 (q1 + q2 )
∆ ∆ (2π) 2q10 (2π) 2q20
1 dΩ |q1 |2 d|q1 |
Z
= 2
δ(E − q10 − q20 ) (6.21)
16π ∆ q10 q20
p
where we have carried out the q2 integral in the last step, so that q10 = q20 = m2 + q21 . Also dΩ = sin θ dθ dφ
is the solid-angle differential for q1 . Since
d |q1 | |q1 | |q1 |E
(q10 + q20 ) = + = (6.22)
d|q1 | q10 q20 q10 q20
74 CHAPTER 6. INTERACTING QUANTUM FIELDS
θ
k 1 = (k 10 , k) k 2 = (k 20 , -k)
q 1 = (q , q)
10
Figure 6.1: Kinematics of 2 → 2 scattering in the center of mass frame. The total center of mass energy is given
by E = k10 + k20 = q10 + q20 .
we can use the delta function to replace the integral over |q1 | in Eq. (6.21) by the inverse of this factor. Inserting
Eq. (6.21)
R in this form together with Eq. (6.20) into Eq. (6.18) and assuming good angular resolution, so that the
integral ∆ can be dropped we have
dσ(k1 , k2 → q1 , q2 ) 1
= |M(k1 , k2 → q1 , q2 )|2 . (6.23)
dΩ 64π 2 E 2
To summarise, this expression provides the differential cross section for a 2 → 2 scattering where all four particles
have the same mass. Note that E is the total center of mass energy. Since we have integrated out the delta-function
in Eq. (6.18) the matrix element should be evaluated for conserved four-momentum.
† †
out hf |i, kiin = out hf |ain (k)|ii = out hf |aout (k)|iiin + out hf |(a†in (k) − a†out (k))|iiin (6.24)
←→
Z
= out hf − k|iiin − out hf |i d3 x e−ikx ∂0 (φin (x) − φout (x)|iiin (6.25)
The creation operators have been expressed in terms of ”in” and ”out” fields by means of Eq. (5.12) which can also
be written in the form
←→
Z
a†in/out (k) = −i d3 x e−ikx ∂0 φin/out (x) , (6.26)
←→
where k0 = wk , as usual. Here, we have used the short-hand notation f (t) ∂t g(t) = f (t)∂t g(t) − ∂t f (t)g(t).
The state |f − kiout refers to the ”out” state |f iout with a particle of momentum k removed (or zero if there
is no such particle in |f iout ). Hence, the first term in Eq. (6.25) either vanishes or corresponds to a forward
scattering contribution where one of the momenta has not changed. We will usually drop this term in the subsequent
calculation since we are interested in the non-trivial scattering part of the amplitude. Taking time limits x0 → ±∞
6.3. THE REDUCTION FORMULA 75
we can convert the ”in” and ”out” fields above into a full interacting field φ by writing
←→
Z
out hf |i, kiin = out hf | lim − lim i d3 x e−ikx ∂0 φ(x)|iiin (6.27)
x0 →∞ x0 →−∞
←→
Z
= out hf |i d4 x ∂0 e−ikx ∂0 φ(x) |iiin (6.28)
Z
= out hf |i d4 x e−ikx ∂02 φ(x) − ∂02 e−ikx φ(x) |iiin
(6.29)
The wave function e−ikx (with k0 = wk ) satisfies the Klein-Gordon equation, so that ∂02 e−ikx = (∇2 −m2 )e−ikx .
Using this in the second term in Eq. (6.29) and subsequently integrating by parts in the spatial directions one finds
Z
out hf |i, kiin = i d4 x e−ikx (2x + m2 ) out hf |φ(x)|iiin . (6.30)
Hence, we have succeeded in removing one particle from the ”in” state of our S-matrix element and replacing it
by a field operator φ(x). The idea is to now repeat this process until all particles, both from the ”in” and the ”out”
state are removed and one is left with a vacuum expectation value of field operators.
This is precisely the same structure as in the first reduction (6.30) apart from the opposite sign in the exponent
eiqy . This sign difference arises because we have reduced a particle in the ”out” state as opposed to a particle in
the ”in” state. In addition, we learn that products of field operators which arise in this way arrange themselves into
time-ordered products.
The vacuum expectation values of time-ordered products which arise in the above formula are also referred to as
N -point Green functions
G (N ) (z1 , . . . , zN ) = h0|T (φ(z1 ) . . . φ(zN ))|0i . (6.35)
Eq. (6.34) shows that the Green functions determine the S-matrix elements and are the crucial objects we need to
calculate. It also displays a remarkable property, the so-called crossing symmetry: All that changes in Eq. (6.34) if
a particle with momentum k is shifted from the ”in” to the ”out” state is the sign of the momentum k in the wave
function e−kx . This means, for example, that a 2 → 2 scattering and a 1 → 3 decay related by transferring one
76 CHAPTER 6. INTERACTING QUANTUM FIELDS
”in” particle to an ”out” particle are described by the same four-point Green function G (4) .
The reduction formula for a complex scalar can now be derived exactly as above. This leads to Eq. (6.34) with one
modification: Positively charged ”in” particle, generated by a†+,in (k), lead to an operator φ† in the time-ordered
product, while negatively charged ”in” states, generated by a†−,in (k), lead to φ. For ”out” states the situation is
reversed, with operators φ for positive particles and operators φ† for negative particles.
(α)
while for a photon with momentum q and polarisation ǫµ (k) in the ”out” state we have in analogy with Eq. (6.33)
Z
out hf, (q, α)|iiin = i d4 x ǫ(α)
µ (q)e
iqx
2x out hf |Aµ (x)|iiin . (6.38)
Repeated reduction leads to time-ordered operator product as for scalar fields. It should also be clear that for
theories with different types of fields, for example for theories with scalar and vector fields, the various reduction
formulae above can be combined and applied successively until all particles are reduced.
where H0 is the free Hamiltonian (quadratic in the fields and conjugate momenta) and H1 = d3 x H1 contains
R
the interactions. For example, for our simple toy model, a real scalar field φ, we have
1 λ 4 λ
Z Z
µ 2 2 3
L= ∂µ φ∂ φ − m φ − φ , H1 = d x H1 = d3 x φ(x)4 . (6.40)
2 4! 4!
So far, we have worked in the Heisenberg picture where the field operator φ(x) is time-dependent and Fock space
states |ai are time-independent. In the Schrödinger picture, on the other hand, the field operator φS (x) is time-
independent while the corresponding Fock space states |a, tiS are time-dependent. The two pictures are related by
the time-evolution operator eiHt in the usual way
For the purpose of perturbation theory it is useful to introduce a third picture, the interaction picture, with fields
and states denoted by φI (x) and |a, tiI . It is intermediate between the two previous pictures and both operators
and states depend on time. In terms of the Schrödinger picture, it is defined through time evolution with the free
Hamiltonion H0 , that is
φI (t, x) = eiH0 t φS (x)e−iH0 t , |a, tiI = eiH0 t |a, tiS . (6.42)
φI (t, x) = U (t, 0)φ(t, x)U −1 (t, 0) , |a, tiI = U (t, 0)|ai , (6.43)
where
U (t, 0) = eiH0 t e−iHt (6.44)
is called the evolution operator. It encodes the difference in time evolution between interaction and Heisenberg
picture due to the interaction Hamiltonian H1 . More generally, we can define the evolution operator U (t, t0 )
between two times t0 and t by |a, tiI = U (t, t0 )|a, t0 iI . This means from Eq. (6.43) that
so, in particular, U −1 (t, t0 ) = U (t0 , t) and U (t0 , t0 ) = 1. For the composition of two evolution operators we have
the rule
U (t, t′ )U (t′ , t0 ) = U (t, 0)U −1 (t′ , 0)U (t′ , 0)U −1 (t0 , 0) = U (t, 0)U −1 (t0 , 0) = U (t, t0 ) . (6.46)
From Eqs. (6.45) and (6.44), the time derivative of the evolution operator is given by
∂ ∂
U (t, 0) U −1 (t0 , 0) = eiH0 t He−iH0 t − H0 U (t, t0 ) .
i U (t, t0 ) = i (6.47)
∂t ∂t
∂
i U (t, t0 ) = H1,I (t)U (t, t0 ) where H1,I = eiH0 t H1 e−iH0 t . (6.48)
∂t
Note that H1,I , the interaction Hamiltonian written in the interaction picture, has the same form as H1 but with φ
replaced by φI .
Perturbative solution for evolution operator
Let us pause for a moment and see what we have achieved. The field φI in the interaction picture evolves with
the free Hamiltonion, as is clear from Eq. (6.42), and we should, hence, think of it as a free field. The evolution
operator relates this free field to the full, interacting field φ in the Heisenberg picture via Eq. (6.43). Therefore, if
we can find an explicit expression for the evolution operator we have succeeded in writing the interacting field in
terms of a free field. To do this, we need to solve the differential equation (6.48) subject to the initial condition
U (t0 , t0 ) = 1. It turns out that the solution is given by
Z t ∞
(−i)p t 4
X Z Z t
U (t, t0 ) = T exp −i dt1 H1,I (t1 ) = 1 + d x1 . . . d4 xp T (H1,I (x1 ) . . . H1,I (xp )) .
t0 p=1
p! t0 t0
(6.49)
The exponential form of this solution is as expected for a differential equation of the form (6.48) and the initial
condition is obviously satisfied. However, the appearance of the time-ordering operator is perhaps somewhat
surprising. To verify the precise form of this solution, we focus on the terms up to second order in H1,I . This
will be sufficient to illustrate the main idea and justify the appearance of the time-ordering. While the term linear
Rt
in H1,I in Eq. (6.49) contains a single integral t0 dt1 and is easy to differentiate with respect to t the situation
Rt Rt
is more complicated for the quadratic term which contains two integrals t0 dt1 t0 dt2 . We re-write this double-
Rt Rt
integral by splitting the integration region into two parts, for t1 ≥ t2 and t1 < t2 , so that t0 dt1 t0 dt2 =
78 CHAPTER 6. INTERACTING QUANTUM FIELDS
Rt R t1 Rt R t2
t0 dt1 t0 dt2 + t0 dt2 t0 dt1 . Then we have
Z t
1 t
Z t
∂ ∂
Z
i U (t, t0 ) = i 1−i dt1 H1,I (t1 ) − dt1 dt2 T (H1,I (t1 )H1,I (t2 )) + O(3)
∂t ∂t t0 2 t0 t0
Z t Z t1 Z t Z t2
i ∂
= H1,I (t) − dt1 dt2 H1,I (t1 )H1,I (t2 ) + dt2 dt1 H1,I (t2 )H1,I (t1 ) + O(3)
2 ∂t t0 t0 t0 t0
Z t Z t
i
= H1,I (t) − H1,I (t) dt2 H1,I (t2 ) + dt1 H1,I (t1 ) + O(2)
2 t0 t0
Z t
= H1,I (t) 1 − i dt1 H1,I (t1 ) + O(3) = H1,I (t)U (t, t0 ) + O(3)
t0
The main point is that the order of the interaction Hamiltonians in the two integrals in the second line is reversed as
a consequence of time-ordering. This allows us, after differentiating in the third line, to factor H1,I (t) to the left for
both terms in the bracket. Without time-ordering H1,I (t) would have been to the left of the first term in the bracket
and to the right of the second. Since interaction Hamiltonians at different times do not necessarily commute this
would have been a problem. The general proof is a straightforward generalisation of the above calculation to all
orders.
G (N ) (z1 , . . . , zn ) = h0|T (φ(z1 ) . . . φ(zN )) |0i = h0|T U −1 (t, 0)φI (z1 ) . . . φI (zN )U (t, 0) |0i
= lim h0|U −1 (t, 0)T (φI (x1 ) . . . φI (xn )U (t, −t)) U (−t, 0)|0i . (6.50)
t→∞
The vacuum should be invariant under the action of the evolution operator up to a constant so we have
We can now insert the explicit solution (6.49) for the evolution operator in the Taylor expanded form into this
result. This leads to a formula for the Green functions which only depends on vacuum expectation values of time-
ordered free-field operator product, each of which can be evaluated using Wick’s theorem. As in Section 1.1.4 each
term which arises in this way can be represented by a Feynman diagram. Feynman diagrams which arise form the
numerator of Eq. (6.53) have N external legs, Feynman diagrams from the denominator are vacuum bubbles. It
turns out that the denominator in Eq. (6.53) precisely cancels all Feynman diagrams from the numerator which
contain disconnected vacuum bubbles. In fact, we have seen an explicit example of this in Eq. (1.79). With this
additional information we can write
∞
(i)p
X Z
G (N ) (z1 , . . . , zN ) = d4 y1 . . . d4 yp h0|T (φI (z1 ) . . . φI (zN )Lint (y1 ) . . . Lint (yp )) |0ino bubbles
p=0
p!
(6.54)
with the interaction Lagrangian Lint = −H1,I . The interaction Lagrangian is proportional to a coupling constant
(λ in the case of our scalar field theory example), so the above expression for the Green functions can be seen as
an expansion in this coupling constant. If its value is sufficiently small it should be a good approximation to only
6.5. FEYNMAN RULES AND EXAMPLES 79
compute a finite number of terms up to a certain power in the coupling. Only in this case is Eq. (6.54) of direct
practical use. In the above derivation we have used notation appropriate for scalar field theory but it is clear that
the basic structure of Eq. (6.54) remains valid for other types of fields. To summarise, combining the reduction
formula (6.34) (and its generalisations to other types of fields) with the perturbation expansion (6.54) and Wick’s
theorem provides us with a practical way of calculating S-matrix elements and, via Eqs. (6.17) and (6.18) decay
rates and cross sections, in terms of Feynman diagrams.
where we denote ”in” momenta ka and ”out” momenta qb collectively by (pA ) = (ka , −qb ). To see where this
leads we should rewrite the LSZ reduction formula (see Eq. (6.34))
Z Y
N
out hq1 , . . . , qm |k1 , . . . , kn iin = i d4 zA e−pA zA (2zA + m2 )G (N ) (z1 , . . . , zN ) (6.56)
A
of Eq. (6.55). Inserting this into the LSZ formula (6.56) we can explicitly carry out the (2zA + m2 ) operations
which now only act on the exponentials in Eq. (6.57) and produce factors −i/∆ ˜ F (QB ) of inverse Feynman prop-
agators, one for each external leg of the Green function. It is, therefore, useful to introduce amputated Green
functions
(N ) G̃ (N ) (p1 , . . . , pN )
G̃amp (p1 , . . . , pN ) = QN , (6.58)
∆˜ F (pA )
A=1
which are related to the ordinary Green functions by removing the propagators for the external legs. The remaining
integerations over zA and QB can then trivially be carried out and we remain with
!
X X
4 4 (N )
out hq1 , . . . , qm |k1 , . . . , kn iin = (2π) δ ka − qb G̃amp (k1 , . . . , kn , −q1 , . . . , −qm ) . (6.59)
a b
Comparison with Eq. (6.8) shows that for real scalar fields the amputated Green function equals the matrix element
M. For other types of fields there are slight modifications. For complex scalar fields, we have to consider Green
functions defined with the appropriate number of fields φ and their conjugates φ† , depending on the number of
particles and anti-particles in the ”in” and ”out” state. Apart from this the above formulae remain valid. For each
vector particle with momentum k and polarisation ǫµ (k) in the ”in” or ”out” state, the Green function carries a
Lorentz index µ and from Eqs. (6.37) and (6.38) this should be contracted into the corresponding polarisation
vector ǫµ (k) to obtain the matrix element M. With this small modification, the above formulae also apply to
vector fields.
p2
p3
p -i λ
1
p4
Figure 6.2: Feynman diagram for 4-point function in real scalar field theory with λφ4 /4! interaction to order λ.
Lint = −λφ4 /4!. Let us first calculate some examples before formulating the general Feynman rules. We start
with the 4-point function to order λ. Dropping disconnected parts, we have from Eq. (6.54)
iλ
Z
G (4) (z1 , z2 , z3 , z4 ) = − d4 y h0|T φ(z1 )φ(z2 )φ(z3 )φ(z4 )φ(y)4 |0i
4!
Z
= −iλ d4 y ∆F (z1 − y)∆F (z2 − y)∆F (z3 − y)∆F (z4 − y)
Z 4
d k1 d4 k4
Z
= −iλ ... ˜ F (k1 ) . . . ∆
d4 y eiy(k1 +···+k4 ) e−i(k1 z1 +···+k4 z4 ) ∆ ˜ F (k4 )
(2π)4 (2π)4
Z 4
4 d k1 d4 k4 4
= −iλ(2π) . . . δ (k1 + · · · + k4 )e−i(k1 z1 +···+k4 z4 )
(2π)4 (2π)4
×∆ ˜ F (k1 ) . . . ∆
˜ F (k4 ) , (6.60)
From the first to the second line we have used Wick’s theorem which only leads to one type of term (although
with multiplicity 24) which can be represented by the Feynman diagram in Fig. 6.2. Then we have used the
representation (5.53) of the Feynman propagator to re-write the expression in momentum space. From this form it
is easy to see that the amputated momentum space Green function is
(4)
G̃amp (p1 , . . . , p4 ) = −iλ . (6.61)
Inserting this into Eq. (6.23) we find the cross section for 2 → 2 scattering
dσ λ2 λ2
= , σ= , (6.62)
dΩ 64π 2 E 2 16πE 2
where E is the total center of mass energy.
Two-point function
Next, we discuss the 2-point function to order λ. From Eq. (6.54) we have
iλ
Z
G (2) (z1 , z2 ) = h0|T (φ(z1 )φ(z2 )) |0i − d4 y h0|T φ(z1 )φ(z2 )φ(y)4 |0i + O(λ2 )
4!
iλ
Z
= ∆F (z1 − z2 ) − ∆F (0) d4 y ∆F (z1 − y)∆F (z2 − y) + O(λ2 )
2
The two corresponding Feynman diagrams are depicted in Fig. 6.3. It is straightforward to Fourier transform this
expression (introducing integration variables z = z1 − z2 and z̃ = (z1 + z2 )/2 in the first term and z̃1 = z1 − y
and z̃2 = z2 − y in the second term). Dividing by the external propagators one then finds for the amputated Green
function
(2)
G̃amp (p, −p) = ˜ F (p)−1 − iλ ∆F (0) + O(λ2 )
∆
2
λ i
Z
2 2 4
= −i p − m + d k 2 + O(λ2 ) (6.63)
2 k − m2 + iǫ
6.5. FEYNMAN RULES AND EXAMPLES 81
+ + ...
p p -i λ p
Figure 6.3: Feynman diagrams for 2-point function in real scalar field theory with λφ4 /4! interaction to order λ.
The diagram at order λ is a loop diagram and, as a result, we have to carry out an integration over the internal loop
momentum k in Eq. (6.63).
Feynman rules
Before we discuss this result, we would like to summarise our experience so far and formulate the Feynman rules
for the amputated Green functions in real scalar field theory.
• For the (connected part of the) N point function at order λp draw all possible (connected) Feynman graphs
with N external legs and p four-vertices. Assign a directed momentum to each line in those Feyman graphs
such that momentum is conserved at each vertex. Then, to each graph associate an expression obtained by:
˜ F (k) for each internal line with momentum k
• writing down a Feynman propagator ∆
• writing down −iλ for each vertex
• writing down an integral d4 k /(2π)4 for each momentum k which is not fixed by momentum conservation
R
at the vertices
• dividing the expression by a symmetry factor which equals the number of permutations of internal lines one
can make for fixed vertices.
The first two of these rules are obvious from Eq. (6.54) and Wick’s theorem. The fact that momentum needs to be
conserved at each vertex can be seen from Eq. (6.60): The calculation which led to this result would be exactly
the same for a single vertex within a larger Feynman diagram. The appearance of the delta function in Eq. (6.60)
then signals momentum conservation at this vertex. Our result (6.61) for the amputated 4-point function shows
that −iλ is the correct factor to include for a vertex. Finally, the precise symmetry factor for each graph follows
from the number of pairings which arise in Wick’s theorem.
scalar field theory with λφ4 interaction should be renormalizable and this is indeed the case. Let us briefly discuss
how one might go about dealing with the above singularity. The first step is to regularise the amplitude. This refers
to some sort of prescription which allows one to assign a finite value to the diverging integral in question. There
are many different ways of doing this but perhaps the simplest and most intuitive one is to introduce a cut-off, that
is, to modify the upper limit of the integral (6.64) to a finite value Λ. If we do this we find
π2 λ 2
G̃(2) 2 2 2 2 2
amp (p, −p) = −i p − m + Λ − m ln(1 + Λ /m ) . (6.65)
2
π2 λ 2
m2R = m2 − Λ − m2 ln(1 + Λ2 /m2 ) .
(6.66)
2
The main idea is that expressed in terms of this mass (and the other remormalized parameters of the theory)
physical quantities are finite provided the theory is renormalizable.
Eq. (6.66) reveals another interesting feature which is characteristic for renormalizable scalar field theories:
the appearance of divergencies quadratic in the cut-off Λ. It turns out that in renormalizable theories without
scalar fields divergencies are at most logarithmic. Quadratic divergencies in scalar field theory lead to a serious
”naturalness” problem called the hierarchy problem. Physically, we can think about the cut-off Λ as the energy
scale above which the theory in question ceases to be valid and new physics becomes relevant. Assume that we
have a theory with a very high cut-off scale Λ but a small physical scalar mass mR . Then Eq. (6.66) requires
a very precise cancelation between the two terms on the right-hand side. For example, in the standard model of
particle physics, a theory shown to be renormalizable, the mass of the Higgs particle should be of the order of
the electroweak symmetry breaking scale, so a few hundred GeV. If the standard model was valid all the way up
to the Planck scale, ∼ 1018 GeV, this would require a 32 digit cancellation in Eq. (6.66). The hierarchy problem
may, therefore, be seen as a reason for why we should expect new physics not far above the electroweak scale.
Supersymmetry is one of the main candidates for such physics and it is indeed capable of resolving the hierarchy
problem.
p q
Figure 6.4: Vertices for scalar electrodynamics. From left to right the corresponding vertex Feynman rules are
−ie(pµ + qµ ), 2iq 2 ηµν and −iλ.
1 1 † λ
L = − Fµν F µν − (∂µ Aµ )2 + [(∂µ + ieAµ )φ] [(∂ µ + ieAµ )φ] − m2 φ† φ − (φ† φ)2 , (6.67)
4 2 4
so we have for the interaction Lagrangian
λ
Lint = ieAµ φ∂µ φ† − φ† ∂µ φ + e2 Aµ Aµ φ† φ − (φ† φ)2 .
(6.68)
4
So, in addition to a quartic scalar field vertex similar to what we have encountered in real scalar field theory, we
expect a triple vertex coupling a vector field to two scalars and another quartic vertex coupling two vectors to two
6.6. FURTHER READING 83
scalars. To find the Feynman rules for these vertices let us compute the appropriate Green functions. Using the
first, cubic interaction in Eq. (6.68) we find for the 3-point function
Z
Gµ(3) (z1 , z2 , z3 ) = −e d4 y h0|T φ(z1 )φ† (z2 )Aµ (z3 )Aν (y)(φ(y)∂ν φ† (y) − φ† (y)∂ν φ(y)) |0i
d4 p d4 q d4 k 4
Z
= −ie(2π)4 δ (p − q + k)e−i(pz1 −qz2 +kz3 )
(2π)4 (2π)4 (2π)4
×∆˜ F (p)∆
˜ F (q)∆
˜ F,0 (k)(pµ + qµ ) , (6.69)
˜ F,0 (k) = i/(k 2 + iǫ) is the Feynman propagator for zero mass. For the amputated Green function this
where ∆
means
(3)
G̃amp,µ (p, −q, k) = −ie(pµ + qµ ) (6.70)
and this is precisely the expression for the triple vertex. From a very similar calculation using two external vector
fields, two external scalars and the second interaction term in (6.68) one finds
(4)
G̃amp,µν (p1 , p2 , p3 , p4 ) = 2ie2 ηµν . (6.71)
Finally, the quartic scalar field interaction from the last term in (6.68) comes with a vertex factor −iλ, just as in
the case of a real scalar field. These Feynman rules for the interactions in scalar electrodynamics are summarised
in Fig. 6.4. Finally, the propagators for internal scalar and vector field lines are
˜ F (k) = i −i
∆ , D̃F (k)µν = ηµν . (6.72)
k 2 − m2 + iǫ k 2 + iǫ
Space-time indices on internal vector field propagators and vertices have to be contracted in the obvious way.
This completes the Feynman rules for the amputated Green functions in scalar electrodynamics. To obtain the
matrix element the µ index for each external photon with momentum k and polarisation α has to be contracted into
(α)
ǫµ (k).
In the previous chapters we have followed the ”traditional”, canonical approach to quantise field theories. We
would now like to show how field theories can be quantised, arguably in a more elegant way, by using path
integrals. We begin by recalling some features of path integrals in quantum mechanics.
for the matrix element between two position eigenstates. Here, |x, ti = eiHt |xiS are Heisenberg picture states
which coincide with Schrödinger picture states |x(t)iS at a given, fixed time t. The central objects to compute
in quantum field theory are vacuum expectation values of time-ordered field operator product. The analogous
quantities in quantum mechanics, for simplicity written for just two operators, are 1 hxf , tf |T (x̂(t1 )x̂(t2 ))|xi , ti i.
Focusing on the case t1 > t2 and inserting complete sets of states we have
hxf , tf |T (x̂(t1 )x̂(t2 ))|xi , ti i = hxf |e−iH(ti −t1 ) x̂S e−iH(t1 −t2 ) x̂S e−iH(t2 −tf ) |xi i
Z
= dx1 dx2 hxf |e−iH(ti −t1 ) |x1 ihx1 |x̂S e−iH(t1 −t2 ) |x2 ihx2 |x̂S e−iH(t2 −tf ) |xi i
Using x̂|xi = x|xi, replacing every expectation value by Eq. (7.1) and combining the three path integrals, together
with the integrations over x1 and x2 into a single path integral this can be written as
Z Z tf
hxf , tf |T (x̂(t1 )x̂(t2 ))|xi , ti i ∼ Dx x(t1 )x(t2 ) exp i dτ L(x, ẋ) . (7.2)
ti
For t2 > t1 the result is actually the same: time-ordering is automatic in path integrals. It is also clear that the
above argument can be repeated for a product of an arbitrary number of operators x̂(t). Further, it can be shown
that limti →−∞,tf →∞ hxf , tf |T (x̂(t1 ) . . . x̂(tN ))|xi , ti i ∼ h0|T (x̂(t1 ) . . . x̂(tN ))|0i, so we have the final result
Z
h0|T (x̂(t1 ) . . . x̂(tN ))|0i ∼ Dx x(t1 ) . . . x(tN )eiS[x] . (7.3)
85
86 CHAPTER 7. PATH INTEGRALS IN QUANTUM FIELD THEORY
In analogy with Eq. (7.3) we can now express these Green functions in terms of a path integral as
Z
G (x1 , . . . , xN ) = N Dφ φ(x1 ) . . . φ(xN )eiS[φ] ,
(N )
(7.5)
where N is a normalization to be fixed shortly. Analogous to what we did in Chapter 1, it is useful to introduce a
generating functional Z Z
W [J] = N Dφ exp i d4 x [L(∂µ φ, φ) + J(x)φ(x)] (7.6)
To fix the normalization N we require that W [J]|J=0 = h0|0i = 1. For the purpose of explicit calculations
it is useful to introduce a Eulkidean or Wick rotated version of the generating functional. To do this we define
Euklidean four-vectors by x̄ = (x̄0 , x̄) = (ix0 , x), associated derivatives ∂¯µ = ∂ x̄∂ µ and a Euklidean version of
the Lagrangian density LE = LE (∂¯µ φ, φ). We can then re-write the generating functional W [J] and obtain its
Euklidean counterpart
Z Z
4 ¯
WE [J] = N Dφ exp d x̄ LE (∂µ φ, φ) + J(x̄)φ(x̄) (7.8)
In analogy with Chapter 1 (see Eq. (1.81) we also define the generating functional Z[J] by
which correspond to connected Feynman diagrams. We will refer to G(n) as connected Green functions and to
Z[J] as the generating functional for connected Green functions.
The full information about the quantum field theory is now encoded in the generating functional, so this is the
primary object to compute. We begin by doing this in the simplest case, the free theory.
with the inverse ∆E (ȳ − x̄) = A−1 (ȳ, x̄) and a suitable normalization Ñ . How do we compute the inverse of the
R d4 k̄ ik̄(x̄−ȳ)
operator A? With the representation δ 4 (x̄ − ȳ) = (2π) 4e of the delta function, we write A as a Fourier
transform
d4 k̄ 2
Z
A(ȳ, x̄) = (k̄ + m2 )eik̄(x̄−ȳ) . (7.16)
(2π)4
Then, we invert A by taking the inverse inside the Fourier transform, that is
d4 k̄ ik̄(x̄−ȳ) 1
Z
∆E (ȳ − x̄) = e . (7.17)
(2π)4 k̄ 2 + m2
Now, we would like to revert to Minkowski space with coordinates x and y. To obtain a Minkowski product in the
exponent in Eq. (7.17) we also need to introduce the Minkowski momentum k = (k0 , k) = (ik̄0 , k̄). By inserting
all this into Eq. (7.15) we find for the generating functional in Minkowski space
1
Z
4 4
W [J] = exp − d x d y J(x)∆F (x − y)J(y) , (7.18)
2
where ∆F is the Feynman propagator precisely as introduces earlier (see Eq. (5.53). We have also chosen Ñ = 1
so that W [J]|J=0 = 1. Eq. (7.18) is the general result for the generating functional of the free scalar field theory.
All the free Green functions can now be calculated from Eq. (7.7) and we already know from Chapter 1 that the
result can be obtained by applying Wick’s theorem. In particular, we have for the 2-point function
δ 2 W [J]
G (2) (x, y) = − = ∆F (x − y) . (7.19)
δJ(x)δJ(y) J=0
Wick’s theorem appears both in the context of canonical quantisation (see Eq. (5.61)) and path integral quantisation.
Applying it in either case shows that both types of Green functions are indeed identical, at least for the case of free
fields. This confirms that canonical and path integral quanitsation are equivalent.
For later reference, we note that from Eqs. (7.18) and (7.10) the generating functional Z for free fields is given
by
i
Z
Z[J] = d4 x d4 y J(x)∆F (x − y)J(y) . (7.20)
2
i δW [J] h0|φ̂(x)|0iJ
φc (x) = − = , (7.22)
W [J] δJ(x) h0|0iJ
where we have suggestively defined the vacuum expectation values h0|0iJ = W [J] and h0|φ̂(x)|0iJ = −iδW [J]/δJ(x)
in the presence of the source J. Eq. (7.22) shows that φc is the suitably normalised vacuum expectation value of
the field operator φ̂, so its interpretation as a classical field is sensible. Note that φc is a function of the source
J. From Eq. (7.10), the generating functional Z is ”on the same footing” as the exponential in the path integral
and can, hence, be seen as some sort of effective action. However, it still contains the effect of the source term
J(x)φ(x) in Eq. (7.6). To remove this source term we define the effective action Γ by a Legendre transform
Z
Γ[φc ] = Z[J] − d4 x J(x)φc (x) . (7.23)
Differentiating the left-hand side of this equation by δ/δJ(y) and using the definition (7.21) of the classical field
it follows that Γ is independent of the source J, as the notation suggests.
88 CHAPTER 7. PATH INTEGRALS IN QUANTUM FIELD THEORY
To see that these definitions conform with our intuition, let us first discuss free fields. In Eq. (7.20) we have
calculated the generating functional Z for free fields and inserting this into Eq. (7.21) we find for the classical field
Z
φc (x) = i d4 y ∆F (x − y)J(y) . (7.24)
Since the Feynman propagator satisfies the equation (2 + m2 )∆F (x) = −iδ 4 (x) (see Eq. (5.55)) we find for the
above classical field that
(2 + m2 )φc (x) = J(x) , (7.25)
so it is a solution to the Klein-Gordon equation with source J(x), as one would expect for a classical field coupled
to a source J. Inserting (7.20) and (7.24) into the effective action (7.23) it follows
1 1 1
Z Z Z
d4 x φc (x)J(x) = − d4 x φc (x)(2 + m2 )φc (x) = d4 x ∂µ φc ∂ µ φc − m2 φ2c
Γ[φc ] = − (7.26)
2 2 2
and, hence, the effective action coincides with the classical action as one would expect for a free theory.
For interacting theories, the generating functional can of course not be calculated exactly. However, we can
proceed to evaluate the path integral
Z Z
4
W [J] = N Dφ exp iS[φ] + i d x J(x)φ(x) (7.27)
in the saddle point approximation as discussed in Chapter 1. The solution φ0 to the classical equations of motion
is determined from the classical equations of motion
δS
[φ0 ] = −J(x) . (7.28)
δφ(x)
With this result for Z and Eqs. (7.21), (7.23) and (7.28) we immediately conclude that φc = φ0 and Γ[φc ] = S[φc ].
Hence, in the lowest order saddle point approximation the effective action Γ is simply the classical action. Beyond
this leading order, the effective action of course receives corrections due to quantum effects and differs from the
classical action. This leads to a systematic approach to calculate these quantum corrections to the effective action
the details of which are beyond the scope of the lecture.
The above formalism also sheds light on another point which we have glossed over so far. Our discussion of
spontaneous symmetry breaking in Chapter 4 has been carried in the context of classical fields and it has not been
obvious what its status should be in the quantum theory. Spontaneous symmetry breaking of a quantum theory
should be analysed using the above effective action (or, more precisely, the effective potential, which is the scalar
potential of the effective action). Hence, we see that the results of Chapter 4 make sense in quantum theory, but
have to be viewed as a leading order approximation.
where Z
δ
N −1 = exp i d4 x Lint −i W0 [J] (7.31)
δJ(x) J=0
to ensure that W [J]|J=0 = 1. This is a perturbative series for the full generating functional in terms of the free
one, W0 [J]. We recall from Eq. (7.18) that the free generating functional is given by
1
Z
W0 [J] = exp − d4 x d4 y J(x)∆F (x − y)J(y) , (7.32)
2
so all the functional differentiations in Eq. (7.30) can be carried out explicitly and lead to Feynman propagators.
The Green functions (7.7) can then written as
δ δ
G (N ) (z1 , . . . , zN ) = N ...
δJ(z1 ) δJ(z1 )
∞ p Z
" #
X i 4 4 δ δ
× 1+ d y1 . . . d yp Lint −i . . . Lint −i W0 [J]
p=1
p! δJ(y1 ) δJ(yp )
J=0
We know from Chapter 1 (see Eq. (1.76) that this expression can be worked out using Wick’s theorem. The
result is a sum over products of Feynman propagators, suitably integrated, and each term can be associated to a
Feynman diagram. This is precisely the same structure as for the perturbative Green function (6.54) obtained from
canonical quantization (after diagrams with disconnected vacuum bubbles are cancelled due to the normalization
factor N ). We have, therefore, explicitly verified that canonical and path integral approach lead to the same
perturbative Green functions. From hereon, working out the Green functions explicitly and calculating decay rates
and cross sections works exactly like in the canonical formalism: We calculate the space-time Green functions
from the above formula and then derive the Fourier transformed and amputated Green functions from Eqs. (6.55)
and (6.58), respectively. Is is clear that this also leads to the same set of Feynman rules, so there is no need to
repeat their derivation. Compared to the hard work in the canonical approach it is remarkable how relatively easily
the path integral formalism delivers the same results.
There are two routes which take us to problems in many-particle quantum physics. One is to start from classical
field theory and quantise, as in Chapter 5 of these notes. The other is to start from one-body or few-body quantum
mechanics, and consider the special aspects which become important when many particles are involved. Clearly,
the first route is the natural one to take if, for example, we want to begin with Maxwell’s equations and arrive at
a description of photons. Equally, the second route is the appropriate one if, for instance, we want to begin with
a model for liquid 4 He and arrive at an understanding of superfluidity. In this chapter we will set out the second
approach and illustrate it using applications from condensed matter physics. Although the problems we cover
can all be formulated using functional integrals, we will use Hamiltonians, operators and operator transformations
instead. This choice is made partly for simplicity, and partly in order to introduce a useful set of techniques.
N N
~2 X 2 X X
H=− ∇i + V (ri ) + U (ri − rj ) .
2m i=1 i=1 i<j
Here there are three contributions to the energy: the kinetic energy of each particle (∇2i operates on the coordinates
ri ); the one-body potential energy V (r); and the two-particle interaction potential U (ri −rj ). To discuss symmetry
under exchange of particles, we define the exchange operator Pij via its action on wavefunctions:
Since [H, Pij ] = 0, we can find states that are simultaneous eigenstates of H and Pij . Moreover, a system that
is initially in an eigenstate of Pij will remain in one under time evolution with H. For these reasons we examine
the eigenvalues of Pij . Since (Pij )2 = 1, these are +1 and −1. Now, it is an observational fact (explained in
relativistic quantum field theory by the spin-statistics theorem) that particles come in two kinds and that particles
of a given kind are always associated with the same eigenvalue of the exchange operator: +1 for bosons and −1
for fermions.
91
92 CHAPTER 8. MANY-PARTICLE QUANTUM SYSTEMS
Next we would like to write down a wavefunction representing an N -particle system with one particle in state l1 ,
one in state l2 and so on. The choice
φl1 (r)φl2 (r) . . . φlN (r)
is clearly unsatisfactory because for general l1 , l2 . . . it has no particular exchange symmetry. Instead we take
X
ψ(r1 . . . rN ) = N (±1)P φk1 (r1 ) . . . φkN (rN ) . (8.1)
distinct perms.
Several aspects of the notation in Eq. (8.1) require comment. The sign inside the brackets in (±1)P is +1 for
bosons and −1 for fermions. The set of labels {k1 . . . kN } is a permutation of the set {l1 . . . lN }. The permutation
is called even if it can be produced by an even number of exchanges of adjacent pairs of labels, and is odd otherwise;
the integer P is even or odd accordingly. The sum is over all distinct permutations of the labels. This means that
if two or more of the labels ln are the same, then permutations amongst equal labels do not appear as multiple
contributions to the sum. Finally, N is a normalisation, which we determine next.
To normalise the wavefunction, we must evaluate
Z Z
d r1 . . . dd rN ψ ∗ (r1 . . . rN )ψ (r1 . . . rN ) .
d
Substituting from Eq. (8.1), we obtain a double sum (over permutations k1 . . . kN and h1 . . . hN ) of terms of the
form Z Z
dd r1 φ∗k1 (r)φh1 (r1 ) . . . dd rN φ∗kN (r)φhN (r1 ) .
These terms are zero unless k1 = h1 , k2 = h2 , and . . . kN = hN , in which case they are unity. Therefore only the
diagonal terms in the double sum contribute, and we have
N!
Z Z X
. . . |ψ|2 = |N |2 (±1)2P = |N |2
n1 ! n2 ! . . .
dist. perms.
where the n1 , n2 . . . are the numbers of times that each distinct orbital appears in the set {l1 . . . lN }, and the ratio
of factorials is simply the number of distinct permutations. Hence we normalise the wavefunction to unity by
taking
1/2
n1 ! n2 ! . . .
N = .
N!
Note that this determinant is zero either if two orbitals are the same (li = lj ) or if two coordinates coincide
(ri = rj ), so the Pauli exclusion principle is correctly built in. Note also that, since the sign of the determinant is
changed if we exchange two adjacent rows, it is necessary to keep in mind a definite ordering convention for the
single particle orbitals φl (r) to fix the phase of the wavefunction.
For bosons, we should use an object similar to a determinant, but having all terms combined with a positive
sign: this is known as a permanent.
A detail to note is that the sum in Eq. (8.4) is over all permutations, while that in Eq. (8.1) included only distinct
permutations. The difference (which is significant only for bosons, since it is only for bosons that we can have
nl > 1), is the reason for the factor (n1 ! n2 ! . . .)1/2 appearing on the right of Eq. (8.3). This choice anticipates
what is necessary in order for boson creation and annihilation operators to have convenient commutation relations.
Annihilation operators appear when we take the Hermitian conjugate of Eq. (8.3), obtaining h0| clN . . . cl2 cl1 .
Let’s examine the effect of creation and annihilation operators when they act on various states. Since c†l |0i is the
state with coordinate wavefunction φl (r), we know that h0|cl c†l |0i = 1, but for any choice of the state |φi other
than the vacuum, c†l |φi contains more than one particle and hence h0|cl c†l |φi = 0. From this we can conclude that
demonstrating that the effect of cl is to remove a particle from the state |nl =1i ≡ c†l |0i. We also have for any |φi
the inner products h0|c†l |φi = hφ|cl |0i = 0, and so we can conclude that
cl |0i = h0|c†l = 0 .
where the upper sign is for bosons and the lower one for fermions. From this we conclude that boson creation
operators commute, and fermion creation operators anticommute: that is, for bosons
[c†l , c†m ] = 0
where we use the standard notation for an anticommutator of two operators A and B: {A, B} = AB +BA. Taking
Hermitian conjugates of these two equations, we have for bosons
[cl , cm ] = 0
Note for fermions we can conclude that (cl )2 =(c†l )2 =0, which illustrates again how the Pauli exclusion principle
is built into our approach.
Finally, one can check that to reproduce the values of inner products of states appearing in Eq. (8.3), we require
for bosons
[cl , c†m ] = δlm
To illustrate the correctness of these relations, consider for a single boson orbital the value of |[(c† )n |0i]|2 . From
Eq. (8.3) we have |[(c† )n |0i]|2 = n!. Let’s recover the same result by manipulating commutators: we have
Of course, manipulations like these are familiar from the theory of raising and lowering operators for the harmonic
oscillator.
or zero for fermions if nl =1. Similarly, by considering the Hermitian conjugate of a similar equation, we have
√
cl |n1 . . . nl . . .i = (±1)n1 +...+nl−1 nl |n1 . . . nl − 1 . . .i ,
or zero for both bosons and fermions if nl =0. In this way we have
c†l cl | . . . nl . . .i = nl | . . . nl . . .i
where the possible values of nl are nl =0, 1, 2 . . . for bosons and nl =0, 1 for fermions. Thus the combination c†l cl ,
which we will also write as n̂l , is the number operator and counts particles in the orbital φl .
Now let c†l create a particle in orbital φl (r), and let d†α create a particle in orbital ρα (r). We can read off from
Eq. (8.5) an expression for d†α in terms of c†l :
X †
d†α = cl Ulα .
l
{dα , d†β } =
X
∗
Ulα Umβ {cl , c†m } = (U† U)αβ = δαβ .
lm
Similarly, for boson operators commutation relations are preserved under unitary transformations.
We want to represent  using creation and annihilation operators. As a first step, we can characterise A(r) by
its matrix elements, writing Z
Alm = φ∗l (r)A(r)φm (r)dd r .
Then X
A(r)φm (r) = φl (r)Alm . (8.6)
l
The second-quantised representation is X
 = Apq c†p cq . (8.7)
pq
To justify this, we should verify that reproduces the correct matrix elements between all states from the Fock space.
We will simply check the action of  on single particles states. We have
X
Â|φm i = Apq c†p cq c†m |0i .
pq
c†p cq c†m |0i = c†p ([cq , c†m ] + c†m cq )|0i = c†p δqm |0i
96 CHAPTER 8. MANY-PARTICLE QUANTUM SYSTEMS
so
X
Â|φm i = |φp iApm ,
p
Again, to justify this one should check matrix elements of the second-quantised form between all states in Fock
space. We will content ourselves with matrix elements for two-particle states, evaluating
by two routes. In a first-quantised calculation with ± signs for bosons and fermions, we have
1
Z Z
hAi = [φ∗x (r1 )φ∗y (r2 ) ± φ∗x (r2 )φ∗y (r1 )] · [A(r1 , r2 ) + A(r2 , r1 )] · [φa (r1 )φb (r2 ) ± φa (r2 )φb (r1 )]dd r1 dd r2
2
1
= [Axyba ± Axyab + Ayxab ± Ayxba + Axyba ± Axyab + Ayxab ± Ayxba ]
2
= (Axyba + Ayxab ) ± (Axyab + Ayxba ) . (8.9)
We can simplify the vacuum expectation value of products of creation and annihilation operators such as the one
appearing here by using the appropriate commutation or anticommutation relation to move annihilation operators
to the right, or creation operators to the left, whereupon acting on the vacuum they give zero. In particular
and
h0|cy cx c†l c†m = h0|(δym δxl ± δyl δxm ) .
Note that in order for the operator H to be Hermitian, we require the matrix H to be Hermitian. Since the matrix
H is Hermitian, it can be diagonalised by unitary transformation. Denote this unitary matrix by U and let the
eigenvalues of H be εn . The same transformation applied to the creation and annihilation operators will diagonalise
H. The details of this procedure are as follows. Let
α†l =
X †
ai Uil .
i
Thus the eigenstates of H are the occupation number eigenstates in the basis generated by the creation operators
α†n .
Fermions
Consider for fermion operators the Hamiltonian
which arises in the BCS theory of superconductivity. Note that λ must be real for H to be Hermitian (more
generally, with complex λ the second term of H would read λc†1 c†2 + λ∗ c2 c1 ). Note as well the opposite ordering
of labels in the terms c†1 c†2 and c2 c1 , which is also a requirement of Hermiticity.
The fermionic Bogoliubov transformation is
where u and v are c-numbers, which we can in fact take to be real, because we have restricted ourselves to real
λ. The transformation is useful only if fermionic anticommutation relations apply to both sets of operators. Let us
suppose they apply to the operators d and d† , and check the properties of the operators c and c† . The coefficients
of the transformation have been chosen to ensure that {c†1 , c†2 } = 0, while
The remaining step is to substitute in H for c† and c in terms of d† and d, and pick θ so that terms in d†1 d†2 +d2 d1
have vanishing coefficient. The calculation is perhaps clearest when it is set out using matrix notation. First, we
can write H as
c1
ǫ λ 0 0
1 λ −ǫ 0 †
0 c
H= c†1 c2 c†2 c1 2 +ǫ
2 0 0 ǫ −λ c2
0 0 −λ −ǫ c†1
We pick θ so that
cos θ − sin θ ǫ λ cos θ sin θ ǫ̃ 0
= ,
sin θ cos θ λ −ǫ − sin θ cos θ 0 −ǫ̃
√
where ǫ̃ = ǫ2 + λ2 . Including the other 2 × 2 block of H, we conclude that
H = ǫ̃(d†1 d1 + d†2 d2 ) + ǫ − ǫ̃ .
Bosons
The Bogoliubov transformation for a bosonic system is similar in principle to what we have just set out, but
different in detail. We are concerned with a Hamiltonian of the same form, but now written using boson creation
and annihilation operators:
H = ǫ(c†1 c1 + c†2 c2 ) + λ(c†1 c†2 + c2 c1 ) .
We use a transformation of the form
Note that one sign has been chosen differently from its counterpart in Eq. (8.10) in order to ensure that bosonic
commutation relations for the operators d and d† imply the result [c†1 , c†2 ] = 0. We also require
and hence u2 − v 2 = 1. The bosonic Bogoliubov transformation may therefore be parameterised as u = cosh θ,
v = sinh θ.
We can introduce matrix notation much as before (but note some crucial sign differences), with
c1
ǫ λ 0 0
1 λ ǫ 0 0 c†2
H= c†1 c2 c†2 c1 0 0 ǫ λ c
−ǫ,
2
2
0 0 λ ǫ c†1
where for bosons we have used the commutator to write c† c = c c† − 1. Again, we focus on one 2 × 2 block
ǫ λ c1
c†1
c2
λ ǫ c†2
8.3. DENSITY CORRELATIONS IN IDEAL QUANTUM GASES 99
Substituting for c and c† in terms of d and d† , this block of the Hamiltonian becomes
u v ǫ λ u v d1
d†1 d2 .
v u λ ǫ v u d†2
In the fermionic case the matrix transformation was simply an orthogonal rotation. Here it is not, and so we should
examine it in more detail. We have
ǫ[u2 + v 2 ] + 2λuv 2ǫuv + λ[u2 + v 2 ]
u v ǫ λ u v
= .
v u λ ǫ v u 2ǫuv + λ[u2 + v 2 ] ǫ[u2 + v 2 ] + 2λuv
It is useful to recall the double angle formulae u2 + v 2 = cosh 2θ and 2uv = sinh 2θ. Then, setting tanh 2θ =
−λ/ǫ we arrive at
H = ǫ̃(d†1 d1 + d†2 d2 ) − ǫ + ǫ̃ .
with p
ǫ̃ = ǫ2 − λ2 . (8.11)
Note that in the bosonic case the transformation requires ǫ > λ: if this is not the case, H is not a Hamiltonian
for normal mode oscillations about a stable equilibrium, but instead represents a system at an unstable equilibrium
point.
and
N N (N − 1)
Z Z Z
hρ(r)ρ(0)i = d3 ri δ(r − ri )δ(ri ) + d3 ri d3 rj δ(r − ri )δ(rj )
L3 L6
N 1
= hρ(r)i2 + 3 δ(r) − 3 .
L L
In the large volume limit we have simply
Now we move to the quantum calculation, taking the average h. . .i to mean an expectation value in number
eigenstates for the orbitals φk (r), weighted by Boltzmann factors at finite temperature. To evaluate hρ(r)i we need
The final term is negligible in the limit L → ∞ with N/L3 fixed, and we discard it. In this limit we can also make
the replacement
1 X 1
Z
→ d3 p
L3 p (2π)3
so that we obtain finally
2
1
Z
2
hρ(r)ρ(0)i = hρ(r)i + hρ(r)iδ(r) ± d3 phnp ieip·r . (8.13)
(2π)3
The final term on the right-hand side of Eq. (8.13) is the correction to our earlier classical result, Eq. (8.12), and
represents the consequences of quantum statistics. Its detailed form as a function of r depends on the momentum
distribution of particles, but the most important features are quite general: there is an enhancement in density
correlations for bosons and a suppression for fermions, on a lengthscale of order ~/∆p, where ∆p is a characteristic
momentum for the gas. In summary: bosons bunch and fermions exclude.
We reproduce them in order to emphasise two points: first, they are more complicated than those for creation
and annihilation operators, since the commutator is itself another operator and not a number; and second, spin
operators acting at different sites commute. We will also make use of spin raising and lowering operators, defined
in the usual way as S + = S x + iS y and S − = S x − iS y .
The Heisenberg Hamiltonian with nearest neighbour ferromagnetic exchange interactions of strength J is
X X 1 + −
Srz Srz′ + Sr Sr′ + Sr− Sr+′
H = −J Sr · Sr′ ≡ −J . (8.14)
2
hrr i
′ ′ hrr i
P
Here hrr′ i denotes a sum over neighbouring pairs of sites on the lattice, with each pair counted once. Thinking
of the spins as classical, three component vectors with length S, the lowest energy states are ones in which all
spins are parallel. The model is simple enough that we can write down the exact quantum ground states as well:
for the case in which the spins are aligned along the positive z-axis, the ground state |0i is defined by the property
Srz |0i = S|0i for all r. Other ground states can be obtained by a global rotation of spin direction. Individual
ground states in both the classical and quantum descriptions break the rotational symmetry of the Hamiltonian,
and so we expect excitations which are Goldstone modes and therefore gapless. There are known as spin waves or
magnons. Eigenstates with a single magnon excitation can in fact be found exactly, but in general to go further we
need to make approximations. As a next step, we set out one approximation scheme.
S z = S − b† b .
From this we might guess S + ∝ b and S − ∝ b† . In an attempt to identify the proportionality constants we can
compare the commutator [S + , S − ] = 2S z with [b, b† ] = 1. Since the commutator is an operator in the first case
and a number in the second, our guessed proportionality cannot be exact, but within states for which hS z i ≈ S
(meaning hS z i − S ≪ S, which can be satisfied only if S ≫ 1) we can take
In an exact treatment, corrections to these expressions form a series in powers of b† b/S. Using this transformation
and omitting the higher order terms, the Hamiltonian may be rewritten approximately as
Xh i
b†r br′ + b†r′ br − b†r br − b†r′ br′ .
X
H = −J S 2 − JS (8.16)
hrr′ i hrr′ i
We use these transformations, and introduce the notation d for vectors from a site to its nearest neighbours, and
z for the coordination number of the lattice (the number of neighbours to a site: six for the simple cubic lattice), to
102 CHAPTER 8. MANY-PARTICLE QUANTUM SYSTEMS
obtain
XX 1
H = −JS 2 N z − JS eir·(k−q) [eid·q − 1]b†q bk
N
rd kq
X
2
= −JS N z − JS ǫq b†q bq ,
q
where
ǫq = 2JS(3 − cos qx − cos qy − cos qz ) .
In this way we have approximated the original Heisenberg Hamiltonian, involving spin operators, by one that is
quadratic in boson creation and annihilation operators. By diagonalising this we obtain an approximate description
of the low-lying excitations of the system as independent bosons. The most important feature of the result is the
form of the dispersion a small wavevectors. For q ≪ 1 we have ǫq = JSq 2 + O(q 4 ), illustrating that excitations
are indeed gapless. The fact that dispersion is quadratic, and not linear as it is, for example for phonons, reflects
broken time-reversal symmetry in the ground state of the ferromagnet.
M v = M v − M ∆v + ~k
Short range repulsive interactions of strength parameterised by u are represented in first-quantised form by
uX
Hint = δ(ri − rj ) .
2
i6=j
uN02 uN0 X h † † † †
i
Hint = + 2c k c k + 2c −k c −k + c k c −k + c k c −k + O([N0 ]0 ) .
2L3 2L3
k6=0
At this stage N0 is unknown, but we can write an operator expression for it, as
c†k ck .
X
N0 = N −
k6=0
It is also useful to introduce notation for the average number density ρ = N/L3 . Substituting for N0 we obtain
uρ uρ X h † i
Hint = N+ ck ck + c†−k c−k + c†k c†−k + ck c−k + O([N0 ]0 )
2 2
k6=0
and hence
uρ 1 Xh i
H= N+ E(k) c†k ck + c†−k c−k + uρ c†k c†−k + ck c−k + . . . (8.17)
2 2
k6=0
with
~2 k 2
E(k) = + uρ .
2m
At this order we have a quadratic Hamiltonian, which we can diagonalise using the Bogoliubov transformation for
bosons set out in Section 8.2.2. From Eq. (8.11), we find that the dispersion relation for excitations in the Bose gas
is
" 2 #1/2
~2 k 2
ǫ(k) = + uρ − (uρ)2 .
2m
At large k (~2 k 2 /2m ≫ uρ), this reduces to the dispersion relation for free particles, but in the opposite limit it
has the form
r
uρ
ǫ(k) ≃ ~vk with v = .
m
In this way we obtain a critical velocity for superfluid flow, which is proportional to the interaction strength u,
illustrating how interactions can lead to behaviour quite different from that in a non-interacting system.
104 CHAPTER 8. MANY-PARTICLE QUANTUM SYSTEMS
• J-P Blaizot and G. Ripka Quantum Theory of Finite Systems (MIT Press) is a useful, clear and complete
advanced reference book.
Chapter 9
Phase Transitions
In this chapter we will examine how the statistical mechanics of phase transitions can be formulated using the
language of field theory. As in relativistic applications, symmetry will be an important guide, but since phenomena
in a condensed matter setting are not constrained by Lorentz invariance, we encounter a variety of new possibilities.
9.1 Introduction
To provide a context, we start by summarising some basic facts about phase transitions. Take first a substance such
as water, that can exist as a solid, liquid or vapour. Consider its phase diagram in the plane of temperature and
pressure, as sketched in Fig. 9.1. The solid is separated from both liquid and vapour by phase boundaries, and one
liquid
solid
p
c
vapour
T
Tc
Figure 9.1: Schematic phase diagram in temperature T and pressure p for a material with solid, liquid and vapour
phases. Line denote phase boundaries; the liquid-vapour critical point has coordinates (Tc , pc ).
cannot get from the solid to one of these other phases without crossing the phase boundary. By contrast, the liquid-
vapour phase boundary ends at a critical point, and while some paths between the liquid and vapour cross this phase
boundary, others do not. On crossing a phase boundary in this phase diagram, the system undergoes a discontinuous
phase transition: properties such as density change discontinuously and there is a latent heat. On the other hand,
if one follows a path between liquid and vapour that avoids the phase boundary by going around the critical point,
properties vary smoothly along the path. As an intermediate case, we can consider a path between liquid and
vapour that goes through the critical point. This turns out to involve a continuous (but sharp) phase transition,
and it is partly behaviour at such critical points that will be focus of this chapter. To emphasise the distinction
between discontinuous, or first order, transitions and continuous ones, it is useful to compare the behaviour of the
heat capacity as a function of temperature in each case, as sketched in Fig. 9.2. At a first order transition there is a
105
106 CHAPTER 9. PHASE TRANSITIONS
latent heat, represented as a δ-function spike in the heat capacity, while at a continuous transition there is no latent
heat, but the heat capacity shows either a cusp or a divergence, depending on the specific example.
C C
T T
Figure 9.2: Schematic dependence of heat capacity on temperature: at a pressure p < pc (left); and at p = pc
(right).
An alternative way of viewing the liquid-vapour transition is to examine behaviour as a function of density
rather than pressure. A phase diagram of this kind is shown in Fig. 9.3. From this viewpoint the consequence of
the first-order transition is that at temperatures T < Tc there is intermediate range of densities which the system
can attain only as a mixture of two distinct, coexisting phases. One of these phases is the liquid, which has a
higher density than that of the system on average, and the other is the vapour, with lower density than average.
As the critical point is approached, the density difference between liquid and vapour reduces. Near the critical
point, as the two phases become more similar, microscopic density fluctuations grow in size: roughly speaking,
such fluctuations involve microscopic regions of vapour appearing within the liquid, or vice-versa. The correlation
length represents the maximum size of such fluctuations, and diverges at the critical point. Sufficiently close to the
critical point it is larger than the wavelength of light, and in these circumstance density fluctuations scatter light
strongly, leading to a phenomenon known as critical opalescence: a cloudiness that appears in fluids close to their
critical point.
ρ liquid
vapour T
Tc
Figure 9.3: Phase diagram in the density-temperature plane. The region of two-phase coexistence is shaded;
elsewhere the system exists in a single phase.
behaviour as a function of magnetic field (in place of pressure) and temperature. The properties of the system,
including most importantly its magnetisation, vary smoothly with applied field above the critical temperature.
Below the critical temperature, and neglecting hysteresis effects, there is a discontinuous change on sweeping
applied field through zero, represented by a first-order phase boundary in the figure. In this case there is a symmetry
under reversal of field and magnetisation that was not evident for the liquid vapour transition. Also shown in
Fig. 9.4 are the two-phase coexistence region (centre) and the behaviour of the magnetic susceptibility (right). The
latter diverges at the critical point: the divergence again reflects the presence of large fluctuations in the critical
region, which are readily polarised by an applied field.
M χ
H
T
Tc T
T
Tc
Figure 9.4: Left: phase diagram for a ferromagnet in the plane of field H vs temperature T , with Tc the critical
point. Centre: the region of two-phase coexistence, in the plane of magnetisation M vs temperature. Right:
behaviour of magnetic susceptibility χ.
defined as
T − Tc
t= .
Tc
The cusp or divergence in the heat capacity C is represented by the exponent α, via C ∼ |t|−α . The way in
which the order parameter m decreases as the critical point is approached is described by the exponent β, with
m ∼ |t|β for t < 0, while the divergences in the susceptibility χ and the correlation length ξ are represented as
χ ∼ |t|−γ and ξ ∼ |t|−ν . Finally, at the critical temperature, the variation in the order parameter with field is
written as m ∼ |h|1/δ .
and obtain the free energy for the sample from this partition function in the usual way, as −kB T ln Z. To take a
specific example, consider again the paramagnet-ferromagnet transition, with a real, scalar order parameter field
ϕ(r) that represents the local magnetisation (which we take to have orientations only parallel or antiparallel to
a preferred crystal axis). Time-reversal takes ϕ(r) to −ϕ(r), and we expect this and also spatial inversion to be
symmetries. As a result, an expansion of the free energy density should contain only even powers of ϕ(r) (apart
from a linear coupling to an applied magnetic field h) and only even order derivatives of ϕ(r). Hence we have
a 2 b 4 1
F (ϕ) = ϕ + ϕ + |∇ϕ|2 + . . . − hϕ . (9.1)
2 4 2
What can we say about the coefficients a, b, . . . that appear in this expansion? First, if it makes sense to truncate
the expansion at order ϕ4 , we must have b > 0, so that F (ϕ) is bounded below. (Alternatively, if b < 0, we would
need to include a term in ϕ6 .) Secondly, we note that the minimum of F (ϕ) is at ϕ = 0 for a > 0, and at non-zero
ϕ for a < 0. This suggests that, to describe a phase transition, a should vary with temperature and change sign
in the vicinity of the critical point. We therefore postulate that a = At, with t the reduced temperature and A an
unimportant coefficient.
Unless boundary conditions impose a spatially-varying solution, we expect ϕ to be independent of r. Then with
h = 0 Eq. (9.2) has the solutions r
−a
ϕ=0 and ϕ=± .
b
The first is the minimum of F (ϕ) for a > 0 and a maximum for a < 0; the other solutions are real only for a < 0
when they are minima. From this we can conclude that the order parameter varies as
Thus the critical exponent β takes the value β = 1/2. At the critical point ϕ is non-zero only for non-zero h, with
the dependence
ϕ ∝ |h|1/3 .
From this we recognise the critical exponent value δ = 3. To evaluate the susceptibility χ ≡ ∂ϕ/∂h we should
consider Eq. (9.2) for non-zero h. Rather than attempting to find the solution for ϕ explicitly as a function of h, it
is more convenient to differentiate the equation directly, giving
∂ϕ
(a + 3bϕ2 ) =1
∂h
and hence
1
At t>0
χ=
1
− 2At t<0
We do so by diagonalising the quadratic form in the argument of the exponentials: since our system is translation-
invariant, this is done by Fourier transform. To be explicit about the Fourier transforms, we consider a system
of linear size L in each direction with periodic boundary conditions. Then we can introduce wavevectors k =
2π
L (l1 , . . . ld ) with li ’s integer, and write
1 1
X Z
ϕ(r) = ϕk eik·r and ϕk = dd re−ik·r .
Ld/2 Ld/2
k
Also, we note from the definition of ϕk that ϕ∗k = ϕ−k , so that, as independent quantities we can take the real
and imaginary parts of ϕk for one half of all wavevectors – say those with component k1 > 0. This means we can
write the functional integration as multiple integrals over just these components, with
Z Y Z Z
Dϕ = dReϕk dImϕk .
k1 >0
For general d one finds that this correlation function falls off on a scale set by the correlation length ξ given by
ξ −2 = a. In d = 3 the integral gives
1 −r/ξ
hϕ(0)ϕ(r)i = e .
2πr
Since we have made the connection ξ ∝ t−1/2 , we have obtained the value for another critical exponent: ν = 1/2.
110 CHAPTER 9. PHASE TRANSITIONS
1.4 1.4
1.2 1.2
1
1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2 0
0 −0.2
−0.5 0 0.5 1 1.5 2 2.5 −0.5 0 0.5 1 1.5 2 2.5
Figure 9.5: Sketch of Eq. (9.3). Left: for a large and positive. Right: for a small and positive. The minimum of
F (ϕ) jumps discontinuously from ϕ = 0 to a finite value of ϕ as a is reduced, representing a first-order transition.
a 2 b 4 1
F (ψ) = |ψ| + |ψ| + |∇ψ|2 + . . . .
2 4 2
a 2 b 4 1 ∗
F (ψ) = |ψ| + |ψ| + ψ (−i~∇ − qA)2 ψ + . . . .
2 4 2m
+ + +
− −
+ +
+ +
− −
+ +
Figure 9.6: Typical low-temperature configuration of a two-dimensional Ising model. The system has positive
magnetisation, but domains of spins with the opposite orientation appear as thermal excitations.
Our conclusions about order in the Ising model in fact also apply to other systems in which the broken symmetry is
discrete. Systems with a continuous broken symmetry behave differently, however: their low energy excitations are
long-wavelength Goldstone modes, which are more effective at disrupting long-range order than are sharp domain
walls in the Ising model. To understand the effect of Goldstone modes on long-range order, we should calculate
the correlation function for fluctuations ϕ(r) in the order parameter. A calculation similar to the one set out above
in Section 9.2.2 leads to the result
kB T eik·r
Z
hϕ(0ϕ(r)i ∼ dd k . (9.4)
(2π)d J k2
If the notion that the system has long range order is to be self-consistent, fluctuations should be small. In particular
a divergence in hϕ2 (0i would signal an instability. The integral on the right-hand side of Eq. (9.4) is divergent at
small k in one and two dimensions, and this indicates the absence of spontaneously broken symmetry at non-zero
temperature in one and two-dimensional systems. (The integral is divergent at large k is two and more dimensions,
but this divergence does not have the same physical significance because in a condensed matter system there is
always an upper limit to wavevectors, set by the inverse atomic spacing.)
Summarising, we have found that the value of the lower critical dimension is dl = 1 for systems with discrete
symmetry, and dl = 2 for systems with continuous symmetry.
= ξ −d χ . (9.5)
9.4. FURTHER READING 113
The last equality in Eq. (9.5) follows from the definition of susceptibility, which gives
R
Dϕ ϕ(0) e−F +h ϕ(r)
R
∂ ∂
χ= hϕ(0)i = R R
∂h h=0 ∂h Dϕ e−F +h ϕ(r)
Z h=0
= dd r [hϕ(0)ϕ(r)i − hϕ(0)ihϕ(r)i] .
We are now in a position to express both the mean square fluctuations and the order parameter in terms of the
reduced temperature t and the critical exponents. We have ξ −d χ ∼ |t|dν−γ and hϕi2 ∼ |t|2β . If fluctuations close
to the critical point are to be small compared to the mean, we require
for |t| small. This is the case if νd − γ > 2β, which is to say
2β + γ
d> =4,
ν
where the ratio of exponents has been evaluated using the results of mean field theory presented earlier in this
chapter. Our conclusion, then, is that while mean field theory provides a qualitatively correct treatment of phase
transitions, for systems in three dimensions it is not quantitatively accurate to neglect fluctuations. In fact, accurate
calculations of exponent values for systems in fewer than four dimensions require a more serious treatment of
interactions in field theory, using renormalisation group methods.