Quantum Field Theory: Roberto Casalbuoni
Quantum Field Theory: Roberto Casalbuoni
Quantum Field Theory: Roberto Casalbuoni
Roberto Casalbuoni
Dipartimento di Fisica
Universit`a di Firenze
1 Lectures
Contents
Index
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1 Introduction
1.1 Major steps in quantum eld theory . . . . . . . . . . . . . . . . . .
1.2 Many degrees of freedom . . . . . . . . . . . . . . . . . . . . . . . . .
1.3 Linear atomic string . . . . . . . . . . . . . . . . . . . . . . . . . . .
3
3
4
5
Klein-Gordon field
Relativistic quantum mechanics and its problems
Quantization of the Klein-Gordon eld . . . . . .
The Noethers theorem for relativistic elds . . .
Energy and momentum of the Klein-Gordon eld
Locality and causality in eld theory . . . . . . .
The charged scalar eld . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
25
25
28
34
38
41
46
4 The
4.1
4.2
4.3
4.4
4.5
4.6
4.7
4.8
Dirac field
The Dirac equation . . . . . . . . . . . . . . . . . . . . . . . .
Covariance properties of the Dirac equation . . . . . . . . . .
Free particle solutions of the Dirac equation . . . . . . . . . .
Wave packets and negative energy solutions . . . . . . . . . .
Electromagnetic interaction of a relativistic point-like particle
Non relativistic limit of the Dirac equation . . . . . . . . . . .
Charge conjugation, time reversal and PCT transformation . .
Dirac eld quantization . . . . . . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
51
51
54
59
65
67
73
76
81
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
101
101
107
111
113
115
119
. . . . .
. . . . .
. . . . .
in QED
. . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
136
136
138
146
149
157
9 Applications
9.1 The cross-section . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9.2 The scattering e+ e + . . . . . . . . . . . . . . . . . . . . .
9.3 Coulomb scattering . . . . . . . . . . . . . . . . . . . . . . . . . . .
164
. 164
. 166
. 171
10 One-loop renormalization
10.1 Divergences of the Feynman integrals . . .
10.2 Dimensional regularization of the Feynman
10.3 Integration in arbitrary dimensions . . . .
10.4 One loop regularization of QED . . . . . .
10.5 One loop renormalization . . . . . . . . . .
10.6 Lamb shift and g 2 . . . . . . . . . . . .
.
.
.
.
.
.
. . . . .
integrals
. . . . .
. . . . .
. . . . .
. . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
175
175
183
184
187
193
198
Chapter 1
Introduction
1.1
1924 Bose and Einstein introduce a new statistics for light-quanta (photons).
1925
January - Pauli formulates the exclusion principle.
July - Heisenbergs rst paper on quantum mechanics (matrix mechanics).
September - Born and Jordan extend Heisenbergs formulation of quantum mechanics to electrodynamics.
1926
January - Schrodinger writes down the wave equation.
February - Fermi introduces a new statistics (Fermi-Dirac).
August - Dirac relates statistics and symmetry properties of the wave
function, and shows that the quantized electromagnetic eld is equivalent
to a set of harmonic oscillators satisfying the Bose-Einstein statistics.
1927
March - Davisson and Germer detect the electron diraction by a crystal.
October - Jordan and Klein show that quantum elds satisfy commutation rules.
1928
January - The Dirac equation.
January - Jordan and Wigner introduce anticommuting elds for describing particles satisfying Fermi-Dirac statistics.
3
January - Pauli and Heisenberg develop the analog for elds of the Lagrangian and Hamiltonian methods of mechanics.
Klein and Nishina complete the theory of the scattering Compton based
on the Dirac equation.
1929
March - Weyl formulates gauge invariance and its relation to charge conservation.
December - Dirac introduces the notion of hole theory, identifying a hole
with a proton.
1931 Dirac proposes the positron to interpret the energy negative solutions of his
equation and Heisenberg introduces the idea of antiparticles.
1932 Anderson detects the positron.
1934 Dirac and Heisenberg evaluate the vacuum polarization of the photon. First
battle with innities in quantum eld theory.
1936 Serber introduces the concept of renormalized charge.
1947 Bethe evaluates the Lamb-shift.
1948 Schwinger ends the calculation of the Lamb-shift and the renormalization
program starts.
1.2
Aim of this course is to extend ordinary quantum mechanics, which describes non
relativistic particles in interaction with given forces, to the relativistic case where
forces are described by elds, as for the electromagnetic case. The most relevant
dierences between the two cases are that the forces become dynamical degrees of
freedom, and that one needs a relativistic treatment of the problem. In order to get
a consistent description we will need to quantize the eld degrees of freedom.
The concept of eld is a very general one. A eld represents a physical quantity
depending on the space-time point. Examples are the distribution of temperatures
in a room, the distribution of the pressure in the atmosphere, the particle velocities
inside a uid, the electric and magnetic elds in a given region of space. The common
physical feature of these systems is the existence of a fundamental state, for example:
pressure or temperature = state with T = constant, or P = constant
particle velocities in a uid = state at rest
electromagnetic eld = state of vacuum.
4
1.3
qn-1
a
qn
qn+1
Fig. 1.1 - In the upper line the atoms are in their equilibrium position, whereas in
the lower line they are displaced by the quantities qn .
and the inter-atomic distance is a. Therefore L = N a. The equations of motion are
the following
qn = 2 [(qn+1 qn ) + (qn1 qn )] = 2 [qn+1 + qn1 2qn ]
(1.1)
as it follows immediately from the expression of the potential energy of the system
N
1 2
(qn qn+1 )2
U=
2 n=1
(1.2)
In order to dene the problem one has to specify the boundary conditions, although
in the N limit we do not expect that they play any role. Usually one considers
two possible boundary conditions
Periodic boundary conditions, that is qN +1 = q1 .
Fixed boundary conditions, that is qN +1 = q1 = 0.
To quantize the problem is convenient to go to the hamiltonian formulation. The
hamiltonian is given by (pn = qn )
H =T +U =
N (
)
1
p2n + 2 (qn qn+1 )2
2 n=1
(1.3)
The equations of motion can be diagonalized by looking for the eigenmodes. Let us
put
qn(j) = Aj eikj an eij t
(1.4)
where the index j enumerates the possible eigenvalues. Notice that in this equation
the dependence on the original equilibrium position has been made explicit through
qn(j) = q (j) (xn ) = q (j) (na) eikj xn
(1.5)
where xn = na is the equilibrium position of the nth atom. By substituting eq. (1.4)
into the equations of motion we get
(
j2 = 4 2 sin2
kj a
2
kj a
2
(1.6)
(1.7)
m = 1, 2, ...
(1.8)
kj = kj + 2m ,
a
correspond to the same j . This allows us to restrict kj to be in the so called rst
Brillouin zone, that is |kj | /a. Let us now take into account the boundary
6
2
4 j
3
2
1
-3
-2
-1
ak j
(1.9)
from which
kj aN = 2j
(j = integer)
(1.10)
2
2
j=
j
aN
L
j = 0, 1, 2, . . .
N
2
(1.11)
where we have taken N even. The restriction on j follows from considering the rst
Brillouin zone (|kj | /a). Notice that the possible values of kj are 2(N/2) + 1 =
N + 1, and that j = 0 corresponds to a uniform translation of the string (with zero
frequency). Since we are interested only in the oscillatory motions, we will omit this
solution in the following. It follows that we have N independent solutions
qn(j) = Aj eij t eiakj n
(1.12)
pn =
2
jn Qj
e N
N
2
jn P
N j
N
(1.13)
(1.14)
Pj = Pj
(1.15)
(1.16)
n=1
2
2
(j
j)n
i
(j
j)
N
N
= N
1
n=1 e
n=0 e
2
(j j)(N + 1)
1e N
=
1=0
2
i (j j)
1e N
i
(1.17)
whereas for j = j the sum gives N . By using this equation we can invert the
previous expansions
N
n=1
qn e
2
2
j n i (j j)n Qj
N
= N Qj
=
e N
N
n
j
(1.18)
obtaining
2
N
i jn
1
Qj =
qn e N
N n=1
(1.19)
2
N
i jn
1
Pj =
pn e N
N n=1
(1.20)
N/2 (
|Pj |2 + j2 |Qj |2
(1.21)
j=1
This is nothing but the hamiltonian of N decoupled harmonic oscillators each having
a frequency j , as it can be seen by putting
Pj = Xj + iYj
Qj = Zj + iTj
(1.22)
The result we have obtained so far shows that the string of N atoms is equivalent to
N decoupled harmonic oscillators. The oscillator modes are obtained through the
8
L = Na
(1.23)
(1.24)
(1.25)
Let us now consider the continuum limit of this system. Physically this is equivalent
to say that we are looking at the system at a scale much bigger than the inter-atomic
distance. We will dene the limit by taking a 0, by keeping xed the length of
the string, that is to say
a 0,
N ,
aN = L sso
(1.26)
The quantity u(xn , t) goes to a function of the variable x dened in the interval
(0, L). Furthermore
u(xn , t) u(xn1 , t)
u (x, t)
(1.27)
a
and
(u(xn+1 , t) u(xn , t)) (u(xn , t) u(xn1 , t)) a(u (xn+1 , t) u(xn , t))
a2 u (x, t)
(1.28)
The equation of motion becomes
u(x, t) = a2 2 u (x, t)
(1.29)
Let us recall that the quantity appearing in the equation of motion for the eld
is the elastic constant divided by the mass of the atom. In order to give a sense to
the equation of motion in the previous limit we need that diverges in the limit.
One could say that in order the string has a nite mass in the continuum, the mass
of each atom must go to zero. That is, we will require
lim a = v
a0
nite
(1.30)
where v has the dimensions of a velocity. We see that in the limit we get the equation
for the propagation of waves with velocity given by v
u(x, t) = v 2 u (x, t)
9
(1.31)
N
+1
n=1
1 L
dx
a 0
(1.32)
from which
[
]
1 L
2
dx (u(x,
t))2 + v 2 (u (x, t))
2a 0
To get nite energy we need also a redenition of the eld variable
u(x, t) = a(x, t)
H=
getting nally
(
)
1 L
t) 2 + v 2 ( (x, t))2
dx (x,
2 0
The normal modes decomposition becomes
u(x, t)
qn (t) ikj an Qj
e
a
a
aN
j
2
j,
L
(1.34)
H=
(x, t) =
(1.33)
(1.35)
(1.36)
< j < +
(1.37)
2
i
jx
1 +
(x, t) =
e L Qj (t)
L j=
(1.38)
kj =
giving rise to
j2 = 4 2 sin2
j
a
j 4 2
L
L
)2
a2 = (akj )2 v 2 kj2
(1.39)
In the continuum limit the frequency is a linear function of the wave vector. The
relation between the normal modes Qj (t) and the eld (x, t) can be inverted by
using the following relation
L
0
dx eix(k k ) = Lk,k
(1.40)
which holds for k and k of the form (1.37). The hamiltonian is easily obtained as
H=
|Q j |2 + v 2 kj2 |Qj |2
(1.41)
j=1
The main result here is that in the continuum limit the hamiltonian of the system
describes an innite set of decoupled harmonic oscillators. In the following we will
show that the quantization of eld theories of the type described in this Section
gives rise, naturally, to a description in terms of particles.
10
Chapter 2
Lagrangian formalism for
continuum systems and
quantization
2.1
String quantization
N/2 (
|Pj |2 + j2 |Qj |2
(2.1)
j=1
where
Qj = Qj ,
kj a
,
2
whereas in the continuum case
j2 = 4 2 sin2
H=
Pj = Pj
kj =
(
2
j,
L
|j| = 1, 2, . . . ,
|Pj |2 + j2 |Qj |2
(2.2)
N
2
(2.3)
(2.4)
j=1
with
j = v|kj |
(2.5)
and kj given by eq. (2.3). In both cases the quantization is trivially done by
introducing creation and annihilation operators
aj =
j
1
Qj + i
Pj ,
2
2j
aj
11
j
1
Qj i
Pj
2
2j
(2.6)
(2.7)
[aj , ak ] = [aj , ak ] = 0
(2.8)
and
aj =
j
1
Qj + i
P =
2
2j j
implying
aj
(2.9)
j
1
Qj + i
Pj
2
2j
(2.10)
j
1
Qj i
Pj = aj
2
2j
(2.11)
We see that aj e aj are 2N (in the discrete case) independent operators as Qj and
Pj . The previous relations can be inverted to give
1
Qj =
(aj + aj ),
2j
Pj = i
j
(aj aj )
2
(2.12)
H=
j [aj aj + aj aj + 1] =
j=1
N/2
j aj aj +
j=N/2
1
2
(2.13)
j
j
(2.14)
(2.15)
In the continuum limit the energy of the fundamental state is innite (we will come
back later on this point). The generic energy eigenstate is obtained by applying
to the fundamental state the creation operators (the space generated in this way is
called the Fock space)
|nN/2 , , nN/2 =
1
(a
)nN/2 (aN/2 )nN/2 |0
(nN/2 ! nN/2 !)1/2 N/2
12
(2.16)
The state given above can be thought of being formed by nN/2 quanta of type
N/2 of energy N/2 , up to nN/2 quanta of type N/2 of energy N/2 . In this kind
of interpretation the nj quanta (or particles) of energy j are indistinguishable one
from each other. Furthermore, in a given state we can put as many particles we
want. We see that we are describing a set of particles satisfying the Bose-Einstein
statistics. Formally this follows from the commutation relation
[ai , aj ] = 0
(2.17)
from which the symmetry of the wave-function follows. For instance a two-particle
state is given by
|i, j = ai aj = |j, i
(2.18)
As we have already noticed the energy of the fundamental state becomes innite
in the continuum limit. This is perhaps the most simple of the innities that we will
encounter in our study of eld quantization. We will learn much later in this course
how it is possible to keep them under control. For the moment being let us notice
that in the usual cases only relative energies are important, and then the value of
E0 (see eq. (2.15)) is not physically relevant. However there are situations, as in the
Casimir eect (see later) where it is indeed relevant. Forgetting momentarily these
special situations we can dene a new hamiltonian by subtracting E0 . This can be
done in a rather formal way by dening the concept of normal ordering. Given an
operator which is a monomial in the creation and annihilation operators, we dene
its normal ordered form by taking all the annihilation operators to the right of the
creation operators. We then extend the denition to polynomials by linearity. For
instance, in the case of the hamiltonian (2.4) we have
: H : N (H) =
j
j
N (aj aj + aj aj ) =
j aj aj
(2.19)
Coming back to the discrete case, recalling eqs. (1.13) and (1.14)
qn =
2
jn Qj
e N ,
N
i
j
pn =
2
jn P
N j
N
(2.20)
2
2
(jn km) 1
i i j(n m)
e N
[Qj , Pk ] =
= inm
e N
N
N j
i
jk
(2.21)
2
i
jx
1 +
(x, t) =
e L Pj
L j=
13
(2.22)
from which
2
i
(jx ky)
1
t)] =
ijk
[(x, t), (y,
e L
L jk
2
i
j(x y)
i +
=
= i(x y)
e L
L j=
(2.23)
This relation could have been obtained from the continuum limit by recalling that
(x, t)
u(x, t)
qn
a
a
(2.24)
implying
m , t)] = i
[(xn , y), (x
nm
a
(2.25)
In the limit
nm
= (x y)
a
if for a 0, xn x, and xm y. In fact
lim
(2.26)
a0
1=
nm
a
a
nm
dx lim
a0 a
(2.27)
L
qn
14
(2.29)
pn qn L
(2.30)
At the classical level the time evolution of the observables is obtained through the
equation
A = {A, H}
(2.31)
where the Poisson brackets can be dened starting form the brackets between the
canonical variables {qn , pm } = nm . The theory is then quantized through the rule
[., .] i{., .}
(2.32)
In the next Section we will learn how to extend the lagrangian and hamiltonian
formalism to the continuum case.
2.2
We will now show how to construct the lagrangian starting from the equations of
motion. For the string this can be simply done by starting from the kinetic energy
and the potential energy. Let us start recalling the procedure in the discrete case.
In this case the kinetic energy is given by
T =
=
N
N
1
1
p2n =
u 2 (xn , t)
2 n=1
2 n=1
N
1
1 L 2
a 2 (xn , t)
(x, t)dx
2 n=1
2 0
(2.33)
N
1
2 (qn qn+1 )2
2 n=1
N
1
2 (u(xn , t) u(xn+1 , t))2
2 n=1
N
1
2 a((xn , t) (xn+1 , t))2
2 n=1
(2.34)
N
1
(xn , t) (xn+1 , t)
U=
av 2
2 n=1
a
15
)2
v2 L 2
(x, t)dx
2 0
(2.35)
and
(2.36)
[
]
1 L
2
dx 2 (x, t) v 2 (x, t)
(2.37)
L=T U =
2 0
The important result is that in the continuum limit, the lagrangian can be written
as a spatial integral of a function of the eld and its rst derivatives, which will
be called lagrangian density, and having the expression
)
1 ( 2
2
v 2
2
L=
(2.38)
L=
0
Ldx
(2.39)
Of course, this is not the most general situation one can envisage, but we will
consider only the case in which the lagrangian density is a local function of the eld
and its derivatives
, x, t)dx
L = L(, ,
(2.40)
Furthermore, we will consider only theories in which the lagrangian contains at
most the rst derivatives of the elds. The reason is that otherwise one can run into
problems with the conservation of probability.
Given the lagrangian, the next step is to build up the action functional. The
extrema of the action give rise to the equations of motion. The action is given by
t2
S=
t1
t2
Ldt =
dt
, x, t)
dxL(, ,
(2.41)
t1
We require that S is stationary with respect to those variations that are consistent
with the boundary conditions satised by the elds. If is the spatial surface
delimiting the region of spatial integration (for the string reduces to the end
points), we will ask that
(x, t) = 0 on
(2.42)
Furthermore we will require that the variations at the times t1 and t2 are zero at
any space point x
(x, t1 ) = (x, t2 ) = 0, at any x
(2.43)
In the discrete case we have only boundary conditions of the second type, but here
the rst ones are necessary in order to be consistent with the boundary conditions
16
for the eld. Let us now require the stationarity of S with respect to variations
satisfying the previous boundary conditions (2.42) and (2.43)
t2
0 = S =
dt
t2
dxL =
L
L L
dx
+
+
dt
t1
t1
(2.44)
Integrating by parts
t2
[ L
]L
L
L
0 =
dt dx
+
t
t
t1
(
) (
)
]
L
L
+
x
x
[
]t2
L
dx
t2
dt
t1
t2
+
t1
t1
L
dt
L
L
L
dx
t
x
(2.45)
The boundary terms are zero due to eqs. (2.42) and (2.43). Then from the arbitrariness of within the region of integration, we get the Euler-Lagrange equations
L
L
L
=0
t x
(2.46)
In fact can be chosen to be zero everywhere except for a small region around any
given point x (see Fig. (2.1).
Fig. 2.1 - Here the arbitrary variation (x) is chosen to be zero all along the
string, except for a small region around the point x.
This discussion can be easily extended to the case of N elds i , i = 1, . . . , N
(think, as an example, to the electromagnetic eld), and to the case of n spatial
dimensions with points labelled by x , = 1, . . . , n. In this case the structure of
the action will be
)
(
t2
i
(2.47)
S=
dt
dn xL i , i ,
x
t1
V
17
Here V is the spatial volume of integration. We will require again the stationarity of
the action with respect to variations of the elds satisfying the boundary conditions
i (x , t) = 0,
on , for any t, t1 t t2
(2.48)
for any x V
(2.49)
The rst boundary conditions are required because, in the general case, one requires
the elds to go to zero at the boundary of the spatial region (usually the innite).
The Euler-Lagrange equations one gets in this case are
L
L
L
= 0,
i t i x i
i = 1, . . . , N,
= 1, . . . , n
(2.50)
To go to the hamiltonian description one introduces the momentum densities conjugated to the elds i :
L
i =
(2.51)
i
and the hamiltonian density
H=
i i L
(2.52)
= ,
L
= v 2 ,
L
=0
(2.53)
From which one recovers the equations of motion for the eld . Furthermore
=
implying
(2.54)
)
)
1 2 1 2 2
1( 2
2
H = L = 2
v
+ v 2
(2.55)
=
2
2
2
which coincides with the energy density given in eq. (2.36).
A big merit of the lagrangian formalism is the possibility to formulate in a simple
way the symmetry properties of the theory. We shall see later on that this is due
to the rst theorem of Emmy Noether which allows to put in a direct relation
the symmetry properties of the lagrangian and the conservation laws. Due to this
correspondence it is also possible to make use of the theorem in a constructive way,
that is to restrict the possible forms of the lagrangian from the requirement of a given
set of symmetries. We will discuss later on the theorem. For the moment being we
will show how the equations of the vibrating string give rise to conservation laws.
18
+
=
dx + v = v
dx
dt
a
a
b
[
]
[
]
= v 2
b
= v2
dx
a
x
a
(2.57)
where we have made use of the equations of motion in the second step. Dening the
local quantity
P (x, t) = v 2
(2.58)
which is the analogous of the Poyntings vector in electrodynamics, we get
dE(a, b)
= [P (b, t) P (a, t)]
dt
(2.59)
This is the classical energy conservation law, expressing the fact that if the energy
decreases in the segment [a, b], then there must be a ux of energy at the end
points a and b. The total energy is conserved due to the boundary conditions,
P (0, t) = P (L, t). But the previous law says something more, because it gives us a
local conservation law, as it follows by taking the limit b a. In fact, in this limit
E(a, b) (b a)H
(2.60)
H P
+
=0
(2.61)
t
x
This conservation law can be checked by using the explicit expressions of H and P ,
and the equations of motion.
2.3
As we have seen in Section 2.2, in a eld theory one denes the density of conjugated
momenta as
L
(2.62)
i =
i
it is then natural to assume the following commutation relations
[i (x , t), j (y , t)] = iij n (x y )
= 1, . . . , n,
19
i, j = 1, . . . , N
(2.63)
and
[i (x , t), j (y , t)] = 0,
[i (x , t), j (y , t)] = 0
(2.64)
In the string case we have = and we reproduce eq. (2.23). Starting from the
previous commutation relations and expanding the eld in terms of normal modes
one gets back the commutation relations for the creation and annihilation operators.
Therefore we reconstruct the particle interpretation. Using the Heisenberg representation (but omitting from now on the corresponding index for the operators),
the expansion of the string eld in terms of creation and annihilation operators is
obtained through the eqs. (2.22) and (2.12).
2
1 i jx
(x, t) =
e L Qj
L j
2
)
1 1 i jx (
e L
=
aj (t) + aj (t)
2j
L j
1
2j
2
2
i
jx
i
jx
L
aj (t) + e L aj (t)
e
(2.65)
(2.66)
(2.67)
(2.68)
from which
aj (t) = aj (0)eij t aj eij t ,
and
(x, t) =
1
2j
i
e
(2.69)
2
2
jx j t
i
jx j t
L
L
aj + e
aj
(2.70)
From this equation one gets immediately the commutation rules for the creation
and annihilation operators, but before doing that let us notice the structure of the
previous expansion. This can be written in the following way
(x, t) =
20
(2.71)
with
1
fj (x, t) =
2j L
or
2
i
jx ij t
e L
(2.72)
1
fj (x, t) =
ei(kj x ij t)
2j L
(2.73)
2
j
L
(2.74)
The functions fj (x, t) and their complex conjugated satisfy the wave equation
2
2 fj (x, t)
2 fj (x, t)
v
=0
t2
x2
(2.75)
(2.76)
(2.77)
L
0
(2.78)
where
()
At B = A(t B) (t A)B
(2.79)
1
eikj (x y) = (x y)
j
(2.80)
Evaluating this expression with two fj (x, t)s or two fj (x, t)s one gets zero. As far
as the second relation is concerned we get
L
0
()
dxfj (x, t)it fl (x, t)
=
0
1 l + j L
=
dxeix(kl kj ) ei(j l )t
L 2 j l 0
1 l + j i(j l )t
=
e
Ljl = jl
(2.81)
L 2 j l
21
Also in this case, by taking two fj (x, t)s or two fj (x, t)s, the result is zero due to
the factor l j .
We repeat that the set fj (x, t) is a complete set of orthonormal solutions of the
wave equation with periodic boundary conditions. A legitimate question is why
()
the operator t appears in these relations. The reason is that the scalar product
should be time independent (otherwise two orthonormal solutions at a given time
could loose these feature at a later time). For instance, in the case of the Schrodinger
equation, we dene the scalar product as
d3 x (x, t)(x, t)
(2.82)
[
]
d 3
3
d x (x, t)(x, t) =
d x + =
d3 x [i(H) i (H)] = 0
dt
(2.83)
In the present case we can dene a time independent scalar product, by considering two solutions f and f of the wave equation, and evaluating the following two
expressions
]
[
L
2
2
f
2
v
dxf
=0
(2.84)
t2
x2
0
2
2 f
2 f
v
dx
f =0
t2
x2
(2.85)
f
f
f
f
dx
f
f v2
f
f
t
t
t
x
x x
)]
(2.86)
If both f and f satisfy periodic boundary conditions, the second term is zero, and
it follows that the quantity
L
()
dxft f
(2.87)
0
is a constant of motion. Using eq. (2.78), we can invert the relation between eld
and creation and annihilation operators. We get
L
0
()
and therefore
aj =
22
(2.89)
and
aj =
L
0
()
(2.90)
k )(y,t) ]
dxdy[(ifj ifj )(x,t) , ifk if
dxfj it fk = jk
=
=
0
()
(2.91)
[aj , ak ] = [aj , ak ] = 0
(2.92)
We have seen that the total energy of the string is a constant of motion. There
is another constant which corresponds to the total momentum of the string, dened
by
L
L
P =
dx P =
dx
(2.93)
0
We will show in the following that this expression is just the total momentum of the
string, by showing that its conservation derives from the invariance of the theory
under spatial translations. For the moment being let us check that this is in fact a
conserved quantity:
L
L
dP
1 2
2
+ ) =
=
dx(
dx
( + v 2 ) = 0
dt
x
2
0
0
(2.94)
where we have used the equations of motion of the string and the boundary conditions. By using the eld expansion
P =
j,l
dx j kl [fj aj fj aj ][fl al fl al ]
j kl
dx
ei(kj x j t) aj ei(kj x j t) aj
2L j l
j,l
ei(kl x l t) al ei(kl x l t) al
1
l
(2.95)
The rst two terms in the last step give zero contribution because they are antisymmetric in the index of summation(kl l). Therefore
P =
1
kj [aj aj + aj aj ] =
kj aj aj
2 j
j
23
(2.96)
kj = 0
(2.97)
for the antisymmetry j. We see that P has an expression similar to that of H (see
eq. (2.19). We deduce that the states
(aN/2 )nN/2 (aj )nj |0 = |
(2.98)
H| = (nN/2 N/2 + + nj j + )|
(2.99)
P | = (nN/2 kN/2 + + nj kj + )|
(2.100)
have energy
and a momentum
as it follows from
[H, aj ] = j aj ,
[P, aj ] = kj aj
(2.101)
j aj aj
(2.102)
[A, aj ] = j aj
(2.103)
A=
we have
24
(2.104)
Chapter 3
The Klein-Gordon field
3.1
The extension of quantum mechanics to the relativistic case gives rise to numerous
problems. The diculties originate from the relativistic dispersion relation
E 2 = |p|2 + m2
(3.1)
E = |p|2 + m2
(3.2)
It is not dicult to convince himself that the solutions with negative energy have
unphysical behaviour. For instance, increasing the momentum, the energy decreases!
But their presence is not a real problem at a classical level. In fact, we see from
eq. (3.2) that there is a gap of at least 2m between the energies of the two types of
solutions. At the classical level, the way in which the energy is transferred is always
a continuous one. So there is no way to start with an energy positive particle and
nish with a negative energy one. On the contrary, in quantum mechanics one can,
through the emission of a quantum of energy E > 2m, go from positive energy to
negative energy states. Since a system behaves in such a way to lower its energy, all
the positive energy states would migrate to negative energy ones, causing a collapse
of the usual matter. In fact we shall see that it is not possible to ignore this kind
of solutions, but they will be reinterpreted in terms of antiparticles. This will allow
us to get rid of the problems connected with the negative energy solutions, but it
will cause another problem. In fact, one of the properties of antiparticles is that
they may be annihilated or pair created. Let us suppose now to try to localize
a particle on a distance of the order of its Compton wave-length, that is of order
1/m. By doing that we will allow an uncertainty on the momentum of about m,
due to the uncertainty principle. This means that the momentum (and the energy)
25
of the particle could reach values of order 2m, enough to create a pair particleantiparticle. This will be possible only violating the conservation of energy and
momentum. Again, this is the case if the violation of energy conservation is on a
time-scale of order t x 1/m. But this is the scale of the Compton wavelength, therefore the attempt of localization will be nullied by the fact that at the
same scale we start pair creating particles and antiparticles, meaning that we will
be unable to dene the concept of a localized single particle. At the Compton scale
there is no such a thing as a particle, but the picture we get from the previous
considerations is the one of a cloud of particles and antiparticles surrounding our
initial particle, and there is no way to distinguish our particle from the many around
it.
These considerations imply that the relativistic theories cannot be seen as theories at a xed number of particles, which is the usual way of describing things in
ordinary quantum mechanics. In this sense a eld theory, as far as we have seen till
now, looks as the most natural way to describe such systems. In fact, it embeds,
in a natural way, the possibility of describing situations with variable number of
particles.
One can look also at dierent ways leading to the necessity of using eld theories.
For instance, by looking at the quantization of the electromagnetic eld, physicists
realized that this gives a natural explanation of the particle-wave duality, and that
in the particle description one has to do with a variable number of photons. On the
contrary, physical entities as the electrons, were always described in particle terms
till 1927, when Davisson and Germer showed experimentally their wave-like behavior. This suggested that the particle-wave duality would be a feature valid for any
type of waves or particles. Therefore, based on the analogy with the electromagnetic
eld, it is natural to introduce a eld for any kind of particle.
Historically, the attempt of making quantum mechanics a relativistic theory
was pursued by looking for relativistic generalizations of the Schrodinger equation.
Later it was realized that these equations should be rather used as equations for
the elds describing the corresponding particles. As we shall see, these equations
describe correctly the energy dispersion relation and the spin of the various particles.
Therefore they can be used as a basis for the expansion of the eld in terms of
creation and annihilation operators. In order to illustrate this procedure, let us
start considering the Schrodinger equation for a free particle
i
= H
t
(3.3)
|p|2
1 2
=
||
2m
2m
(3.4)
(3.5)
H=
|p|2 + m2
(3.7)
2
2
i
=
|| + m
t
(3.8)
2 + m2
||
v
u
u
)
2
(
||
t
2 k
= m 1 2 = m
ck ||
(3.9)
)
2 ( 2
2
=
|
|
+
m
t2
(3.10)
This equation is both local and invariant under Lorentz transformations, in fact we
can write it in the following form
(
2 + m2 = 0
(3.11)
where
2
2
||
(3.12)
t2
is the DAlembert operator in (3 + 1) dimensions. Notice that in order to solve
the diculties we have listed above
we have been obliged to consider both
types of
solutions: positive energy, E = |p|2 + m2 , and negative energy E = |p|2 + m2 .
The equation we have obtained in this way is known as the Klein-Gordon equation.
As relativistic extension of the Schrodinger theory it was initially discarded because
it gives rise to a non denite positive probability. In fact, if and are two
solutions of such an equation, we can write the following identity
2 =
0 = 2 + m2 ) 2 + m2 ) = [ ( )]
27
(3.13)
(3.14)
d x J0 =
d3 x( )
(3.15)
p = (E, p),
= 0, 1, 2, 3
(3.16)
The metric tensor g is diagonal with components (+1, 1, 1, 1). The fourmomentum operator in coordinate space is given by
(
p i
= i , i
x
t
(3.17)
= 2
x x
(3.18)
x p = Et p x
3.2
(3.19)
In this Section we will discuss the quantization of the Klein-Gordon eld, that is a
eld satisfying the equation (3.11). The quantization will be performed by following
the steps we have previously outlined, that is
construction of the lagrangian density and determination of the canonical momentum density (x);
quantization through the requirement of canonical commutation relations
[(x, t), (y, t)] = i 3 (x y),
We start by the construction of the lagrangian, requiring that the related EulerLagrangian equation gives rise to the Klein-Gordon equation. To this end let us
recall how one proceeds in the discrete case. Suppose to have a system of N degrees
of freedom satisfying the following equations of motion
mi qi =
V
qi
(3.21)
t2
dt
t1
mi qi qi =
i=1
t2
dt
t1
i=1
qi
V
qi
(3.23)
Integrating by parts
{
N
1
dt
mi qi2 V
2 i=1
t2
t1
]}
]t2
mi qi qi
i=1
=0
(3.24)
t1
Using the boundary conditions we see that if the equations of motion are satised,
than the lagrangian, as dened by
t2
S=
t1
N
1
mi qi2 V dt
2 i=1
(3.25)
is stationary. Conversely from the requirement that the action is stationary under
variations satisfying eq. (3.22), the equations of motion follow. Analogously, in the
Klein-Gordon case, we multiply the equation by arbitrary local variations of the
lim (x, t) = 0
(3.26)
t2
0=
dt
t1
)
( ) (
dx
+
+ m2
t
3
(3.27)
t2
0=
dt
t1
1
1
1 m2 2
d x 2
2
2
2
3
(3.28)
d3 xL
L=
29
(3.29)
with
]
1[
m2 2
2
In fact, we have just shown that the quantity (the action)
L=
(3.30)
t2
S=
dtL
(3.31)
t1
is stationary at the point in which the equations of motion are satised. We can
now write down the canonical momentum density
=
L
=
(3.32)
t), (y,
t)] = 0
[(x, t), (y, t)] = [(x,
(3.33)
Let us now construct a complete set of solutions of the Klein-Gordon equation. First
of all we need a scalar product. But we have already one, because we have shown in
the previous Section that the Klein-Gordon equation admits a conserved quantity
(see eq. (3.15)), therefore, if f and g are two solutions, the scalar product is
f |g = i
d3 xf t g
()
(3.34)
f = A(k)eikx = A(k)ei(k0 x0 k x)
(3.35)
(3.36)
k 2 = m2 = k02 = |k|2 + m2
(3.37)
from which
To x the normalization, we proceed as in the one-dimensional case by taking a nite
volume and requiring periodic boundary conditions (normalization in the box). By
taking a cube of side L we require
(x + L, y, z, t) = (x, y + L, z, t) = (x, y, z + L, t) = (x, y, z, t)
it follows
(3.38)
k = 2 n
L
(3.39)
(3.40)
where
30
d3 xfk t fk = k,k
()
(3.41)
ni ,ni
(3.42)
i=1
with n e n are two vectors with integer components, related to k and k , by the
relation (3.39). It follows
d3 xAk Ak ei(k0 k0 )x0 i(k k ) x (k0 + k0 ) = k,k
(3.43)
2
(n1 n1 )x
dxe L
= Ln1 ,n1
(3.44)
d3 xei(k k ) x = L3 k,k
(3.45)
(3.46)
Using
we get
from which
()
i
V
where
2 2 2
|n| + m2
(3.47)
L
By considering, for the moment being, the positive solution of this equation, we
obtain
)
(
1
1
2 2 2
2
Ak = 3/2
, k =
|n| + m = |k|2 + m2
(3.48)
L
L
2k
k02 =
1 ikx
1
e
L3/2 2k
(3.49)
Often we will make use also of the so called normalization in the continuum.The
space integration is then extended to all of R3 and we require
fk |fk = i
()
d3 xfk t fk = 3 (k k )
(3.50)
In this case the spatial momentum can assume all the possible values in R3 . It
follows
(3.51)
d3 xAk Ak eikx ik x (k0 + k0 ) = (2)3 3 (k k )|Ak |2 2k0
31
Ak =
(2)3 2k
where
k =
(3.52)
|k|2 + m2
(3.53)
We see that one goes from the normalization in the box to the normalization in the
continuum through the formal substitution
1
1
V
(2)3
(3.54)
1
(2)3
1 ikx
e
2k
(3.55)
(3.56)
k0 = |k|2 + m2 = k
(3.57)
As a consequence we get two kind of wave functions having positive and negative
energy and behaving as eik x0 and eik x0 , k > 0, respectively. The second kind
of solutions has negative norm in the scalar product we have dened. This would
be a big problem if this equation had the same interpretation as the Schrodinger
equation. In the eld theory, no such a problem exists. In fact, the physical Hilbert
space is the Fock space, where the scalar product is between the states build up
in terms of creation and annihilation operators. Having two types of solutions the
most general expansion for the eld operator (in the Heisenberg representation) is
1
(x) =
(2)3
d k
a(k)eik x0 + ik x + a
(k)eik x0 + ik x
2k
(3.58)
1
1 3
(k)]
d3 k[fk a(k) + fk a
d k
a(k)eikx + a
(k)eikx
(x) =
3
2k
(2)
(3.59)
32
Notice that the energy positive and negative solutions are orthogonal (remember
the one-dimensional case discussed in Section 2.3). We can then invert the previous
expansion with the result
a(k) = i
a
(k) = i
()
d3 x(x)t fk (x)
(3.60)
(3.61)
(x) =
(3.62)
From these equations one can evaluate the commutators among the operators a(k)
e a (k), obtaining
[a(k), a (k )] = 3 (k k )
(3.63)
[a(k), a(k )] = [a (k), a (k )] = 0
(3.64)
These commutation relations depend on the normalization dened for the fk s. For
instance, if we change this normalization by a factor Nk
fk |fk = i
()
d3 xfk t fk = Nk 3 (k k )
(3.65)
=
we get
i 3 ()
a(k) =
d xfk t ,
Nk
(3.66)
i 3
()
a (k) =
d xt fk
Nk
(3.67)
i 3 ()
1 3
d xfk t fk =
(k k )
Nk Nk
Nk
(3.68)
and therefore
[a(k), a (k )] =
For instance, a normalization which is used very often is the covariant one
1 3 1
(x) =
dk
[A(k)eikx + A (k)eikx ]
3
(2)
2k
(3.69)
The name comes from the fact that the factor 1/2k makes the integration over the
three-momentum Lorentz invariant. In fact one has
1 4
1 3 1
dk
d k(2)(k 2 m2 )(k0 )
=
(2)3
2k
(2)4
33
(3.70)
|k|2 + m2 )
(3.71)
i
and therefore
3.3
1
1 ikx
e
(2)3 2k
d3 xfk t fk =
()
1
1 3
(k k )
3
(2) 2k
[A(k), A (k )] = (2)3 2k 3 (k k )
(3.72)
(3.73)
(3.74)
We will now review the Noethers theorem. This allows to relate symmetries of the
action with conserved quantities. More precisely, given a transformation involving
both the elds and the coordinates, if it happens that the action is invariant under
this transformation, then a conservation law follows. When the transformations
are limited to the elds one speaks about internal transformations. When both
types of transformations are involved, it is convenient to evaluate, in general, the
variation of a local quantity F (x) (that is a function of the space-time point)
F (x) = F (x ) F (x) = F (x + x) F (x)
F (x)
= F (x) F (x) + x
x
(3.75)
The total variation keeps into account both the variation of the reference frame
and the form variation of F . It is then convenient to dene a local variation F ,
depending only on the form variation
F (x) = F (x) F (x)
(3.76)
Then we get
F (x)
x
Let us now start form a generic four-dimensional action
F (x) = F (x) + x
S=
V
d4 x L(i , x),
i = 1, . . . , N
(3.77)
(3.78)
and let us consider a generic variation of the elds and of the coordinates, x =
x + x
(3.79)
i (x) = i (x ) i (x) i (x) + x
x
34
(3.80)
SV
=
V
V
V
=
V
i , x )
d4 x L(
d4 xL(i , x)
V
(x )
4 i
d xL( , x + x)
d4 xL(i , x)
(x)
V
i , x + x)(1 + x )
d4 xL(
d4 xL(i , x)
i , x + x) L(i , x)] +
d4 x[L(
[
d4 xL(i , x) x
]
L i
L i
L
dx
+ i , + x +
d4 xL x
i
,
x
V
4
L
L
L
i
dx
+
d4 x Lx + i i
i
i
,
,
V
4
=
V
(3.81)
The rst term in the last line is zero due to the Euler-Lagrange equations of motion
L
L
i = 0
i
(3.82)
d x
L
L
Lx + i i x i i, = 0
,
,
(3.83)
This is the general result expressing the local conservation of the quantity in parenthesis. According to the choice one does for the variations x and i , and of
the corresponding symmetries of the action, one gets dierent kind of conserved
quantities.
Let us start with an action invariant under space and time translations. In the
case we take x = a with a independent on x e i = 0. From the general result
in eq. (3.83) we get the following local conservation law
T =
L i
Lg ,
i, ,
T = 0
(3.84)
T is called the energy-momentum tensor of the system. From its local conservation
we get four constant of motion
d3 xT0
P =
35
(3.85)
L
L i
i =
,
i
,
i,
J = 0
(3.86)
Q=
d3 xJ 0
(3.87)
In general, if the system has more that one internal symmetry, we may have more
that one conserved charge Q, that is we have a conserved charge for any .
The last case we will consider is the invariance with respect to Lorentz transformations. Let us recall that they are dened as the transformations leaving invariant
the norm of a four-vector
2
x2 = x
(3.88)
For an innitesimal transformation
x = x + x
(3.89)
x2 x2 + 2x x = x x = 0
(3.90)
it follows
Since Lorentz transformations are linear
x = x x + x
(3.91)
x x = 0 = x x = 0
(3.92)
we get
The most general solution for the parameters of the transformation is that the
form an antisymmetric second order tensor
=
(3.93)
We see that the number of independent parameters characterizing a Lorentz transformation is six. As well known, three of them correspond to spatial rotations,
whereas the remaining three correspond to Lorentz boosts. In general, the relativistic elds are chosen to belong to a representation of the Lorentz group ( for instance
the Klein-Gordon eld belongs to the scalar representation). This means that under
a Lorentz transformation the components of the eld mix together, as, for instance,
a vector eld does under rotations. Therefore, the transformation law of the elds
i under an innitesimal Lorentz transformation can be written as
1
j
i = ij
2
36
(3.94)
where we have required that the transformation of the elds is of rst order in the
Lorentz parameters . The coecients (antisymmetric in the indices (, ))
dene a matrix in the indices (i, j) which can be shown to be the representative of
the innitesimal generators of the Lorentz group in the eld representation. Using
this equation and the expression for x we get the local conservation law
[(
0 =
L i
1 L ij j
, Lg x +
i
,
2 i,
[
)
(
L
1
j
=
T x T x + i ij
2
,
(3.95)
L ij j
i,
(3.96)
it follows the existence of six locally conserved currents (one for each Lorentz transformation)
M = 0
(3.97)
and consequently six constants of motion (notice that the lower indices are antisymmetric)
M =
d3 xM0
(3.98)
Three of these constants ( the ones with and assuming spatial values) are nothing
but the components of the angular momentum of the eld.
In the case of Klein-Gordon
)
1(
T =
m2 2 g
(3.99)
2
from which
)
1
1( 2
T00 = 2 +
|| + m2 2
(3.100)
2
2
This current corresponds to the invariance under time translations, and it must be
identied with the energy density of the eld (compare with the equation (2.55 for
the one-dimensional case). In analogous way
Ti0 = i
(3.101)
x
is the momentum density of the eld. Using = , the energy and momentum of
the Klein-Gordon eld can be written in the form
)
1 3 ( 2
2 + m2 2
d x + ||
(3.102)
P0 = H =
d3 xT 00 =
2
P =
d xT
0i
d x i ,
x
3
37
(P =
d3 x)
(3.103)
3.4
It is very easy to verify that the energy density found previously coincides with the
hamiltonian density evaluated in the canonical way through the Legendre transformation of the lagrangian density
H = L
(3.104)
We will verify now, that the momentum P is the generator, as it should be, of
the space-time translations. Which amounts to say that it satises the following
commutation relation with the eld
[(x), P ] = i
(3.105)
In fact
[(y , t), H] =
1 3
y , t)
d x[(y , t), 2 (x, t)] = i(y , t) = i(
2
(3.106)
Analogously
[(y , t), P i ] =
(x, t)
(y , t)
(y , t)
]
=
i
=
i
xi
y i
yi
U = eia P
(3.107)
(3.108)
(x)
(x + a) (3.109)
x
With a calculation completely analogue to the one done in Section 2.3 we can
evaluate the hamiltonian and the momentum in terms of the creation and annihilation operators
1 3
d kk [a (k)a(k) + a(k)a (k)]
(3.110)
H=
2
P =
d3 kka (k)a(k)
(3.111)
(3.112)
the corresponding states. In the case of the box normalization, for any k = (2/L)n
(that is for any choice of the three integer components of the vector n, (n1 , n2 , n3 )),
one can build up a state |nk such that
1 ( )nk
|nk =
a (k)
|0
nk !
(3.113)
1
n
|nki =
(a (k1 )) k1 (a (k ))nk |0
nk1 ! nk !
(3.114)
The fundamental state is the one with zero particles in any cell of the momentum
space (vacuum state)
|0 =
|0i
(3.115)
where |0i is the fundamental state for the momentum in the cell i. That is
aki |0i = 0
(3.116)
1
k [a (k)a(k) + a(k)a (k)]
2
(3.117)
and therefore
1
k |0
2
H|0 =
(3.118)
This sum is innite. Recalling that k = (2/L)n, it follows that the cell in the
k-space has a volume
(2)3
Vk =
(3.119)
L3
from which
1
1 Vk k
1 L3 3 2
d k |k| + m2
k =
=
2
2 Vk
2 (2)3
k
(3.120)
which is divergent.
Let us recall that this problem can be formally avoided through the use of the
normal product. In other words by subtracting the innite energy of the vacuum
from the hamiltonian. In the box normalization we have
: H :=
k a (k)a(k)
39
(3.121)
: H :=
d3 kk a (k)a(k)
(3.122)
As we see, the energy of the vacuum depends on the quantization volume. This
implies that it depends on the boundary conditions of the problem. In the real
vacuum this is not a diculty, but it must be considered when one quantize elds
which are inside a nite given volume. In this case this dependence produces measurable eects, as it was pointed out theoretically by Casimir in 1948, and then
proved experimentally by Sparnay in 1958.
L
Fig. 3.1 - The Casimir eect
We will discuss very briey the Casimir eect arising when we have an electromagnetic eld conned between two large perfectly conducting plates. We idealize
the two plates as two large parallel squares of side L at a distance R L. The
theory shows that there is an attractive force per unit surface between the two plates
given by
2 /hc
0.013
p=
=
dyn/cm2
(3.123)
4
240 R
(Rm )4
We can understand the origin of this force in a very qualitative way by quantizing
the electromagnetic eld (that we will take here as a Klein-Gordon eld with zero
mass,m = 0) in a box of side L. The vacuum energy will be
E0 L3
kmax
k d3 k
(3.124)
1/L
with the integration between a lower momentum of order 1/L and an arbitrary upper
momentum which is necessary in order to make nite the integral. If we insert two
plates of side L, as shown in Fig. 3.1, at a distance R, the energy of the eld in this
region, before the introduction of the plates is
E L2 R
kmax
1/L
40
k d3 k
(3.125)
When we insert the plates we get an analogous result, but the lower momentum will
be of order 1/R. therefore the variation of the energy results to be
E L R
2
1/L
1/L
kd k=L R
1/R
1/R
L2 R
k dk =
4
3
[(
1
L
)4
)4 ]
(3.126)
L2
(3.127)
R3
The energy per unit surface behaves as 1/R3 , and the pressure is given by
E
3.5
E/L2
1
4
R
R
(3.128)
For a free particle there are generally three conserved quantum numbers, as the spatial momentum, or energy, angular momentum and its third component. All these
quantities can be expressed as spatial integrals of local functions of the elds. The
locality property is a crucial one and is connected with the causality. To understand
this point let us consider the following example. For a Klein-Gordon free eld there
is a further constant of motion, the number of particles
d3 ka (k)a(k)
N=
We will show
quantity, and
particles. We
theory admits
(3.129)
(3.130)
However this expression vanishes for a hermitian eld. But it turns out that the
operator N can be expressed in terms of the positive energy
(+) (x) =
d3 k
1
2k (2)3
ei(k t k x) a(k)
(3.131)
(3.132)
N=
()
d3 x(+) it ()
41
(3.133)
This is a constant of motion, because both (+) and () are solutions of the equation
of motion, and therefore
j = () ( (+) ) ( () )(+)
(3.134)
is a conserved current. However this current is not a local expression in the eld
. This is because (+) and () are not local functions of . In fact, in order to
project out these components from the eld we need a time integration. In fact, by
dening
(x) = d4 k(k)eikx
(3.135)
with
(k) =
(
)
2k
2
2
(k
m
)
a(k)(k
)
+
a
(k)(k
)
0
0
(2)3
one has
(+)
(x) =
d4 k(+) (k)eikx
(3.136)
(3.137)
with
(+) (k) = (k0 )(k)
(3.138)
(+) (x) =
x )(x )
d4 x (x
(3.139)
But
1 4 ik(x x )
d ke
(k0 )
(2)4
x ) =
(x
Therefore
(+)
(x, x0 ) =
0 x )(x, x )
dx0 (x
0
0
(3.140)
(3.141)
To show the implications of having to do with a non local current, let us dene a
particle density operator
()
N (x) = i(+) t ()
(3.142)
This operator does not commute with itself at equal times and dierent space points
[N (x, t), N (y , t)] = 0,
x = y
(3.143)
However, for local operators, O(), this commutator is automatically zero, due to
the canonical commutation relations
[O((x, t)), O((y , t))] = 0,
42
x = y
(3.144)
We want to argue that the vanishing of this commutator is just the necessary condition in order that O represents an observable quantity. In fact, if the commutator
of a local operator with itself is not zero at space-like distances, then the measure of
the observable at some point, x, would inuence the measures done at points with
space-like separation from x, because we cannot measure the operator simultaneously at two such points. But this would imply the propagation of a signal at a
velocity greater than the light velocity, in contrast with the causality principle. We
see that the vanishing of the commutator of a local observable with itself at spacelike distances is a necessary condition in order to satisfy the causality principle. We
show now that this is automatically satised if the operator under consideration is
a local function of the elds. We will start showing that the commutator of the
eld with itself is a Lorentz invariant function. Therefore, from the vanishing of the
commutator for separations between points of the type x = (t, x), and y = (t, y ),
it follows the vanishing for arbitrary space-like separations. Let us evaluate the
commutator
[(x), (y)] =
[
d3 k1 d3 k2
d3 k
eik(x y) eik(x y)
=
(2)3 2k
d3 k
= 2i
sin(k (x0 y0 ))eik(x y )
3
(2) 2k
(3.145)
d4 k
(k0 )(k 2 m2 ) eik(x y) eik(x y)
(2)3
d4 k
=
(k0 )(k 2 m2 )eik(x y)
(3.146)
3
(2)
[(x), (y)] =
Since the sign of the fourth component of a time-like fourvector is invariant under
proper Lorentz transformations, we see that by putting
[(x), (y)] = i(x y)
(3.147)
the function
(x y) = i
d4 k
(k0 )(k 2 m2 )eik(x y)
(2)3
(3.148)
Therefore the canonical commutation relations make sure the observability for the
Klein-Gordon eld. For the negative and positive energy components we get
(+)
()
(x),
d3 k
eik(x y)
(2)3 2k
d4 k
=
(k0 )(k 2 m2 )eik(x y)
(2)3
(+) (x y)
(y)] =
(3.149)
Also in this case we have a Lorentz invariant function, and therefore it is enough to
study its equal times behaviour:
(+)
d3 k
eik x
3
(2) 2k
2
k dkd(cos )d ikr cos
=
e
(2)3 2k
[
]
1 kdk ikr
ikr
e
e
= i 2
4 r 0 2k
1 d +
eikr
= 2
dk
8 r dr
|k|2 + m2
(0, x) =
(3.150)
(+)
Since
1 d +
(0, x) = 2
deimr sinh
8 r dr
(3.151)
(3.152)
(3.153)
m (1)
H (imr)
8r 1
(3.154)
(1)
The asymptotic behaviour of the Hankels function H1 (imr) for large and small
values of r is given by
(1)
lim H (imr)
r 1
from which
m
lim (+) (0, x)
r
8r
2 mr
e
,
mr
2 mr
e
,
mr
44
(1)
2
mr
(3.155)
1
4r2
(3.156)
lim H1 (imr)
r0
r0
We see that for space-like separations this commutator does not vanish. But for
space separations larger than the Compton wave length 1/m, (+) is practically
zero. Remember that for an electron the Compton wave length is about 3.9 1011
cm. Clearly, an analogous result is obtained for the commutator of the particle
density operator. From this we can derive the impossibility of localize a KleinGordon particle (but the result can be extended to any relativistic particle) over
distances of the order of 1/m. We start dening the following operators
d3 xN (x)
N (V ) =
(3.157)
Vi
ri
V0
x0
R
3.6
We have shown that a hermitian Klein-Gordon eld describes a set of identical scalar
particles. If we want to describe dierent kind of particles we need to introduce
dierent kind of elds. Let us begin with two dierent hermitian scalar elds. The
free lagrangian is a simple sum
L=
2 [
]
1
( i )( i ) m2i 2i
2 i=1
(3.158)
(3.159)
(3.160)
All the considerations done up to now can be easily extended to the case of two elds.
However we notice that there are two kind of creation and annihilation operators,
a i (k), i = 1, 2, and as a consequence
a 1 (k1 )a 2 (k2 )|0 = a 2 (k1 )a 1 (k2 )|0
(3.161)
The two particle state is not any more symmetric, since it is built up in terms of two
dierent types of creation operators. That means that the two elds correspond to
distinguishable particles.
Something really new comes out when the two elds have the same mass term
in the lagrangian
L=
]
[
1
1
[( 1 )( 1 ) + ( 2 )( 2 )] m2 21 + 22
2
2
(3.162)
Then the theory acquires a symmetry under rotations in the plane of the two elds
1 e 2
1 = 1 cos + 2 sin
2 = 1 sin + 2 cos
46
(3.163)
= ( 1 , 2 )
(3.164)
and the norm is invariant under rotations. For innitesimal transformations we have
2 = 1
1 = 2 ,
(3.165)
(3.166)
(3.167)
From the Noethers theorem, we have a conserved current, associated to this symmetry, given by (see eq. (3.86))
J =
L
i = i, ij j
i,
(3.168)
It is convenient to factorize out the angle of the innitesimal rotation and dene a
new current
1
j = J = i, ij j = 1, 2 2, 1
(3.169)
The conservation of the current follows from the equality of the masses of the two
elds, as one can also verify directly
j = ( 2 1 )2 ( 2 2 )1 = (m21 m22 )1 2
(3.170)
Q=
d3 x j 0 =
d3 x( 1 2 2 1 )
(3.171)
[Q, 2 ] = i1
(3.172)
(3.173)
i2 2
iQ
iQ
e
1 e
= 1 + i[Q, 1 ] + [Q, [Q, 1 ]] + . . .
2!
1
2
= 1 + 2 1 + . . .
2
= 1 cos + 2 sin
(3.174)
In fact
47
In analogous way one can show the transformation properties of 2 . The invariance
of L under rotations in the plane (1 , 2 ) is referred to as the invariance under the
group O(2). The real basis for the elds used so far is not the most convenient one.
In fact, the charge Q mixes the two elds. One can understand better the properties
of the charge operator in a basis in which the elds are not mixed. This basis is a
complex one and it is given by the combinations
1
= (1 + i2 ),
2
1
= (1 i2 )
2
(3.175)
(the factor 1/ 2 has been inserted for a correct normalization of the elds). It
follows
1
1
[Q, ] = [Q, 1 + i2 ] = (i2 1 ) =
(3.176)
2
2
and analogously
[Q, ] =
(3.177)
Therefore the eld lowers the charge of an eigenstate of Q by one unit, whereas
increases the charge by the same amount. In fact, if Q|q = q|q
Q(|q) = ([Q, ] + Q)|q = (1 + q)|q
(3.178)
|q |q 1
(3.179)
|q |q + 1
(3.180)
and
In analogous way
i
2 = ( )
2
(3.181)
from which
[
)2
(
)2
1 (
+
4
[
) ]
)
(
1 2 (
2
2
+
m
4
= m2
L =
(3.182)
and
j = 1, 2 2, 1
i
i
= ( + )( ) + ( )( + )
2
[2
]
= i ( ) ( )
48
(3.183)
Q=i
d3 x t
()
(3.184)
(3.186)
Let us notice that these commutation relations could have also been obtained directly
from the lagrangian (3.182), since
L
= ,
L
=
(3.187)
The original O(2) symmetry becomes now an invariance of the lagrangian (3.182)
under a phase transformation of the elds. This follows from (3.182) but it is seen
also from the change of variables
1
(1 + i2 )
2
1
= (1 cos + 2 sin + i(1 sin + 2 cos ))
2
(
)
1
i
i
1 e
+ i2 e
= ei
=
2
and
(3.188)
ei
(3.189)
In this basis we speak of invariance under the group U (1) (the group of unitary
transformations on the complex vectors of dimensions d = 1).
Using the expansion for the real elds
i (x) =
(3.190)
we get
(x) =
1
1
d3 k fk (x) (a1 (k) + ia2 (k)) + fk (x) (a 1 (k) + ia 2 (k))
2
2
(3.191)
1
b(k) = (a1 (k) ia2 (k))
2
49
(3.192)
it follows
(x) =
(x) =
(3.193)
from which we can evaluate the commutation relations for the creation and annihilation operators in the complex basis
[a(k), a (k )] = [b(k), b (k )] = 3 (k k )
(3.194)
[a(k), b(k )] = [a(k), b (k )] = 0
(3.195)
We get also
: P :=
d3 kk
a i (k)ai (k) =
d3 kk a (k)a(k) + b (k)b(k)
(3.196)
i=1
Therefore the operators a (k) e b (k) both create particles states with momentum
k, as the original operators a i . The charge Q is given by
Q = i
()
d3 x
= i
d3 x t
]
()
()
(3.197)
where we have used the orthogonality relations (3.50). For the normal ordered
charge operator we get
: Q :=
d3 k a (k)a(k) b (k)b(k)
(3.198)
(x y)2 < 0
(3.199)
Q(Vi ) =
Vi
d3 x j0 (x)
(3.200)
a(k), and creation, b (k), operators, meaning that a local theory deals in a symmetric
way with the annihilation of a particle and the creation of an antiparticle. For
instance, to annihilate a charge +1 is equivalent to the creation of a charge 1.
50
Chapter 4
The Dirac field
4.1
In 1928 Dirac tried to to solve the problem of a non positive probability density,
present in the Klein-Gordon case, formulating a new wave equation. Dirac thought,
correctly, that in order to get a positive quantity it was necessary to have a wave
equation of the rst order in the time derivative (as it happens for the Schrodinger
equation). Therefore Dirac looked for a way to reduce the Klein-Gordon equation
(of the second order in the time derivative) to a rst-order dierential equation.
The Pauli formulation of the electron spin put Dirac on the right track. In fact,
Pauli showed that in order to describe the spin, it was necessary to generalize the
Schrodinger wave function (a complex number) to a two components object
=
1
2
(4.1)
i
=
H
t
=1
(4.2)
where the hamiltonian H is, in general, a 2 2 matrix. The electron spin is then
described by a special set of 2 2 matrices, the Pauli matrices,
= 1
S
2
(4.3)
Dirac realized that it was possible to write the squared norm of a spatial vector as
|k|2 = ( k)2
(4.4)
[i , j ]+ = 2ij
(4.5)
as it follows from
51
+ m H
= i
t
(4.6)
where
and are matrices. The requirements that this equation should satisfy are
the wave function , solution of the Dirac equation, should satisfy also the
Klein-Gordon equation in order to get the correct dispersion relation between
energy and momentum;
the equation should admit a conserved current with the fourth component
being positive denite;
the equation should be covariant with respect to Lorentz transformations (see
later)
In order to satisfy the rst requirement, we iterate the Dirac equation and ask
that the resulting second order dierential equation coincides with the Klein-Gordon
equation
2
+ m)2
= (i
2
t
(
)
2
i j
2 2
=
+ m i(
+
)
xi xj
(
)
1[ i j]
2
2 2
= ,
+ m i(
+
)
+ xi xj
2
(4.7)
i , j
]
+
= 2ij ,
i ,
]
+
= 0,
2 = 1
(4.8)
(4.9)
(4.10)
where we have made use of the cyclic property of the trace. The consequence is that
the matrices i and can be realized only in a space of even dimensions. This is
perhaps the biggest diculty that Dirac had to cope with. In fact, the i s enjoy
52
the same properties of the Pauli matrices, but in a 2 2 matrix space, a further
anticommuting matrix does not exist. It required some time to Dirac before he
realized that the previous relations could have been satised by 4 4 matrices.
An explicit realization of the Dirac matrices is the following
[
0
=
i
i
i
,
0
1 0
=
0 1
(4.11)
as it can be checked
[
]
+
[i , j ]+
=
0
[
]
+
0
=
i
0
= 2ij
[i , j ]+
[
i
0
+
0
i
(4.12)
i
=0
0
(4.13)
Let us now show that also the second of our requirements is satised. We multiply
the Dirac equation by at the left, and then we consider the equation for
i
)
= i(
+ m
t
(4.14)
+ m) (i
(
+i
= (i
+ m) = i
) (4.15)
t
t
that is
i
We see that the current
( ) + i j ( j ) = 0
t
x
(4.16)
j = ( , i )
(4.17)
is a conserved one.
j
=0
(4.18)
x
Furthermore its fourth component j 0 = is positive denite. Of course we have
still to prove that j is a four-vector, implying that
d3 x
53
(4.19)
4.2
+ m
= i
t
(4.20)
1 0
==
,
0 1
0
= =
i
i
0
(4.21)
i
+ i i i m = 0
0
x
x
0
(4.22)
(4.23)
where
= =
x
i, j
= i j + j i = i , j
0, i
]
+
= , i
]
+
]
+
= 2ij
= i + i = 0
(4.24)
(4.25)
(4.26)
or
Notice that
[ , ]+ = 2g
(4.27)
( i ) = (i ) = i = i
(4.28)
( i )2 = 1
(4.29)
and
The covariance of the Dirac equation means that the following two conditions are
satised
given the Dirac wave function (x) in the Lorentz frame, S, an observer in
a dierent frame, S should be able to evaluate, in terms of (x), the wave
function (x ) describing the same physical state as (x) in S;
54
m (x ) = 0
x
(4.30)
The matrices should satisfy the same algebra as the matrices , because in
both cases the wave functions should satisfy the Klein-Gordon equation (which is
invariant in form). Therefore, neglecting a possible unitary transformations, the two
sets of matrices can be identied. As a consequence, the Dirac equation in S will
be
(
)
i
(4.31)
m (x ) = 0
x
Since both the Dirac equation and the Lorentz transformations are linear, we will
require that the wave functions in two dierent Lorentz frames are linearly correlated
(x ) = (x) = S()(x)
(4.32)
where S() is a 4 4 matrix operating on the complex vector (x) and is the
Lorentz transformation. On physical ground, the matrix S() should be invertible
(x) = S 1 () (x )
(4.33)
but using the relativity principle, since one goes from the frame S to the frame S
through the transformation 1 , we must have
from which
(x) = S(1 ) (x )
(4.34)
S 1 () = S(1 )
(4.35)
i
m (x) = 0
x
we can write
(4.36)
i m S 1 () (x ) = 0
x
(4.37)
=
= ,
x
x x
x
it follows
x = x
(4.38)
iS() S 1 () m (x ) = 0
x
55
(4.39)
or
S() S 1 () =
(4.40)
S 1 () S() =
(4.41)
(4.42)
from which
i
i
1 + 1 = (g +
4
4
(4.44)
1
i
[ , ] = = (g g )
4
2
(4.45)
[ , ] = 2i(g g )
(4.46)
and nally
It is not dicult to show that the solution of this equation is given by
i
= [ , ]
2
(4.47)
In fact
i
i
[ , ] = [ + ]
2
2
i
=
[(2g ) (2g ) + ]
2
= i[ ] = i[ (2g ) (2g ) ]
= 2i[g g ]
(4.48)
[ , ] =
i
= [ , ]
2
56
(4.49)
(4.50)
We can now verify that the current j , dened in eq. (4.17) transforms as a four
vector. To this end we introduce the following notation
(x)
= (x) = (x)0
It follows
0 ,
j 0 = =
(4.51)
i =
i
j i = i =
(4.52)
or
j =
(4.53)
it follows
and therefore
from which
i
i
= [ , ] = [ , ]
2
2
(4.55)
0 0 =
(4.56)
0 S ()0 = S 1 ()
(4.57)
(x ) = (x)S
()
(4.58)
We get
(x) = j (x)
j (x ) = (x ) (x ) = (x)S
() S()(x) = (x)
(4.59)
We see that j has the desired transformation properties. The representation for
the Lorentz generators, in the same basis used previously for the matrices, is
i
i
0
0i = [0 , i ] = ( 2 i i ) = ii = i
i
2
2
i
0
i
i
i
ij = [i , j ] = (i j j i ) = [i , j ]
2
2
2[
]
i [i , j ]
0
=
0
[i , j ]
2
[
]
0
= ijk k
0 k
(4.60)
(4.61)
We see that the generators of the spatial rotations are nothing but the Pauli matrices,
as one should expect for spin 1/2 particles.
57
The behaviour of the Dirac wave function under parity x x can be obtained
in analogous way. In this case
P =
and therefore
1
1
(4.62)
S 1 (P ) S(P ) =
(4.63)
(4.64)
x = (x0 , x)
(4.65)
We are now in the position to classify the bilinear expressions in the Dirac wave
function under Lorentz transformations. Let us consider expressions of the type
A,
where A is an arbitrary 4 4 matrix. As a basis for the 4 4 matrices we can
take the following set of 16 linearly independent matrices
S = 1
V =
A
= 5
T =
P = 5
(4.66)
(4.67)
[5 , ]+ = 0
(4.68)
(4.69)
0 1
(4.70)
1 0
One can easily verify that the bilinear expressions have the following behaviour
under Lorentz transformations
scalar
four vector
5 pseudoscalar
(4.71)
5 =
58
(x)
5 (x) (x )5 (x ) = P P (x)0 5 0 (x) = (x)5 (x)
4.3
(4.72)
In this Section we will study the wave plane solutions of the Dirac equation. In the
rest frame of the particle we look for solutions of the type
(t) = ueimt
(4.73)
where u is a four components complex vector (usually called a spinor). This solution
has positive energy. Substituting inside the Dirac equation we get
(i m)(t) = (m0 m)ueimt = 0
(4.74)
(0 1)u = 0
(4.75)
that is
Since 0 has eigenvalues 1, we see that the Dirac equation admits also solutions
of the type eimt , corresponding to a negative energy state. More generally we can
look for solutions of the form
(+) (x) = eikx u(k),
() (x) = eikx v(k),
positive energy
negative energy
(4.76)
(4.77)
(0 1)u(m, 0) = 0
(0 + 1)v(m, 0) = 0
(4.78)
There are two independent spinors of type u and two of type v satisfying these
equations. In the basis where 0 is a diagonal matrix we can choose the following
solutions
1
0
u(1) (m, 0) = ,
0
0
0
0
v (1) (m, 0) =
,
1
0
1
u(2) (m, 0) =
0
0
0
0
v (2) (m, 0) =
0
1
59
(4.79)
In a general Lorentz frame the solutions could be obtained by boosting the solutions
in the rest frame. Or, we can notice that the following expression
(k m)(k + m) = k 2 m2
(4.80)
vanishes for k 2 = m2 . Therefore we can solve our problem (except for a normalization constant), by putting
u() (k) = c (k + m)u() (m, 0)
v () (k) = d (k + m)v () (m, 0)
(4.81)
(4.82)
(4.83)
()
By taking into account that u() (m, 0) and u (m, 0) are eigenstates of 0 with
eigenvalue +1, we get
()
u() (m, 0) u() (m, 0) = u (m, 0) u() (m, 0)
()
()
= u (m, 0)0 0 u() (m, 0) = u (m, 0) u() (m, 0)
(4.84)
from which
()
u() (m, 0) u() (m, 0) = u (m, 0) u() (m, 0)
()
= g 0 u (m, 0)0 u() (m, 0) = g 0
(4.85)
(4.86)
that is
Than we choose
c =
2m(m + E)
E=
|k|2 + m2
(4.87)
(4.88)
()
and using the fact that v () (m, 0) and v (m, 0) are eigenstates of 0 with eigenvalue
1, we get
()
v() (m, 0) v () (m, 0) = v (m, 0) v () (m, 0)
()
()
= v (m, 0)0 0 v () (m, 0) = v (m, 0) v () (m, 0)
(4.89)
that is
()
v() (m, 0) v () (m, 0) = v (m, 0) v () (m, 0)
()
= g 0 v (m, 0) 0 v () (m, 0) = g 0
(4.90)
(4.91)
1
d = c =
2m(m + E)
(4.92)
and nally
k + m
v () (k) =
v () (m, 0)
2m(m + E)
(4.93)
Notice that positive and negative energy spinors are orthogonal. In the following it
will be useful to express our solutions in terms of two component spinors, () (m, 0)
and () (m, 0)
u
()
[ ()
]
(m, 0)
(m, 0) =
,
0
()
(m, 0) =
0
()
(m, 0)
(4.94)
E
k =
0
]
k
0
0
E
E
k
0
k
k
E
(4.95)
from which
k
m + E ()
0)
() (m, 0)
(m,
2m
2m(m + E)
()
,
v
(k)
=
u() (k) =
()
m
+
E
(m, 0)
()
(m,
0)
2m(m + E)
2m
(4.96)
61
In the following we will need the explicit expression for the projectors of the
positive and negative energy solutions. To this end, let us observe that
1
0
2
0
1
1 + 0
u() (m, 0)
u() (m, 0) = [ 1 0 0 0 ] + [ 0 1 0 0 ] =
(4.97)
0
0
2
=1
0
0
and analogously
1 0
v () (m, 0)
v () (m, 0) =
2
=1
2
(4.98)
(4.99)
u() (k)
u() (k) =
2m(m + E) 2
2m(m + E)
=1
(k + m)2 + 2E(k + m)
1
k 2 + m2 + 2mk + 2E(k + m)
=
=
2m(m + E)
2
4m(m + E)
(2E + 2m)(k + m)
k + m
=
=
(4.100)
4m(m + E)
2m
+ (k) =
Analogously
k m
1 0
k m
v () (k)
v () (k) =
2m(m + E) 2
2m(m + E)
=1
(k m)2 2E(k m)
1
k 2 + m2 2mk 2E(k m)
=
=
2m(m + E)
2
4m(m + E)
(2E + 2m)(k + m)
k + m
=
=
(4.101)
4m(m + E)
2m
(k) =
It is easy to verify that the matrices (k) verify all the properties of a complete
set of projection operators
2 = ,
+ = 0,
+ + = 1
(4.102)
In this normalization the density has the correct Lorentz transformation properties
(+)
(+)
1
u() (m, 0)(k + m)0 (k + m)u() (m, 0)
2m(m + E)
1
=
2E u() (m, 0)(k + m)u() (m, 0)
2m(m + E)
2E(m + E)
E
=
=
(4.103)
2m(m + E)
m
Therefore the density for positive energy solutions transforms as the fourth component of a four vector. The same is true for the negative energy solutions
()
()
() (x)()
We nd also
u
()
(4.105)
In fact
u
()
(4.106)
It follows that solutions with opposite energy and same three momentum are orthogonal
0
() = e+i(Ex + k x) v(k),
k = (E, k)
, k = (E, k)
(4.107)
12 =
3
0
0
3
(4.108)
Clearly u(1) (m, 0) and v (1) (m, 0) are eigenstates of this operator (and therefore of
the third component of the spin operator) with eigenvalues +1, whereas u(2) (m, 0)
63
and v (2) (m, 0) belong to the eigenvalue 1. The projector for the eigenstates with
eigenvalues +1 is given by
1 + 12 n3R
1 + 12
(4.109)
2
2
where nR = (0, 0, 0, 1) is a unit space-like four-vector. Also we have
i
12 = [1 , 2 ] = i 1 2 = 0 5 3 = 5 3 0
2
(4.110)
12 n3R = 5 n
R 0
(4.111)
and
The presence of 0 forbids a simple extension of this expression to a generic Lorentz
frame. We can avoid this, by changing the denition of the projection in the rest
frame system. Let us put
[
(nR ) =
1 12 n3R 0
1 1 3
=
0
2
2
0
1 3
(4.112)
In this case (nR ) and (nR ) project out u(1) (m, 0), v (2) (m, 0) and u(2) (m, 0),
v (1) (m, 0), respectively. That is, (nR ) projects out the positive energy solutions
with spin 1/2 and the negative energy solutions with spin 1/2. Then, we have
(nR ) =
1 5 n
R
2
(4.113)
1 5 n
,
2
n2 = 1,
nk = 0
(4.114)
The projector (n) projects out energy positive states that in the frame system
n = 1/2, and the negative energy states with
have a polarization given by S
polarization S n = 1/2.
In the following we will use the following notation
u(kR , nR )
u(kR , nR )
v(kR , nR )
v(kR , nR )
=
=
=
=
u(1) (m, 0)
u(2) (m, 0)
v (1) (m, 0)
v (2) (m, 0)
(4.115)
(4.116)
and
(nR )u(kR , nR ) = (nR )v(kR , nR ) = 0
64
(4.117)
(n)v(k, n) = v(k, n)
(n)u(k, n) = (n)v(k, n) = 0
(4.118)
(4.119)
(n)2 = (n),
(n)(n) = 0
(4.120)
1 + 5 n
)2
1 + (5 n
)2 + 25 n
2 + 25 n
=
=
= (n)
4
4
(4.121)
1 + 5 n
1 5 n
1 (5 n
)2
=
=0
2
2
4
4.4
(4.122)
As we have shown the Dirac equation leads to a positive probability density. This
solves the problem one had with the Klein-Gordon equation. On the other hand the
Dirac equation does not solve the problem of the negative energy solutions (and it
should not, as we have seen their importance in the Klein-Gordon case). In fact,
the completeness of the spinors involve all the solutions
2 [
u() (k)
u() (k) v () (k)
v () (k) = + (k) + (k) = 1
(4.123)
=1
In the case of a non interacting theory there are no possibilities of transitions among
positive and negative energy states but, when an interaction is turned on, such a
possibility cannot be excluded. In fact, if we try to localize a Dirac particle within
distances of order 1/m the negative energy solutions cannot be ignored. To clarify
this point let us consider the time evolution of a gaussian wave packet, assigned at
time t = 0,
|x|2
2
1
e 2d w
(4.124)
(x, 0) =
(d2 )3/4
where w is a xed spinor, w = (, 0), with w w = 1. As one can check, the wave
packet is normalized to one
|x|2
1
d3 x e d2
d3 x =
(d2 )3/2
65
1
(d2 )3/2
i=1
3
1
(d2 )3/2
x2
2i
dxi e d
(d2 )1/2 = 1
(4.125)
i=1
The solution of the Dirac equation with this boundary condition is obtained by
expanding over all the wave plane solutions
d k
(2)3
(x, t) =
2
m
b(k, )u() (k)eikx + d (k, )v () (k)eikx
E =1
(4.126)
and evaluating the expansion coecients b(k, ) and d (k, ), by requiring that the
solution coincides with eq. (4.124) at time t = 0. We get
d k
(2)3
(x, 0) =
2 [
]
m
)v () (k)
eik x
b(k, )u() (k) + d (k,
E =1
|x|2
1
=
e 2d2 w
2
3/4
(d )
(4.127)
(2)3
2 [
]
m
)v () (k)
|x|2
d3 x e 2d2 eik x
w
(d2 )3/4
|k|2 d2
w
2
(2d2 )3/2 e
(d2 )3/4
(4.128)
From which
2 [
m
)v () (k)
=
b(k, )u() (k) + d (k,
=1
)3/4
|k|2 d2
2 w
e
(4.129)
b(k, ) =
m
E
) =
d (k,
m
E
d2
)3/4
|k|2 d2
()
2 u (k)w
e
(4.130)
|k|2 d2
()
2 v (k)w
e
(4.131)
)3/4
66
b(k, ) =
) =
d (k,
m
E
m
E
)3/4
d2
)3/4
|k|2 d2
m + E ()
2
(m, 0)
e
2m
|k|2 d2
1
2
() (m, 0)k
e
2m(m + E)
(4.132)
(4.133)
from which we can evaluate the ratio of the negative energy amplitudes to the
positive energy ones
)
d (k,
|k|
(4.134)
b(k, )
m+E
The amplitudes (for both signs of the energy) contribute only if |k| 1/d (due to the
gaussian exponential). Suppose that we want to localize the particle over distances
larger than 1/m, that is we require d 1/m. Since the negative energy state
amplitudes are important only for |k| > m 1/d, their contribution is depressed
by the gaussian exponential. On the other hand, if we try to localize the particle
over distances d 1/m, the negative energy states contribution becomes important
for values of |k| of order m, or of order 1/d, that is in the momentum region in
which the corresponding amplitudes are not negligible. We see that the physical
interpretation is essentially the same following from the uncertainty principle.
4.5
Before continuing our discussion about the properties of the Dirac equation, let us
describe the interaction of a point-like particle with the electromagnetic eld in the
relativistic formalism.
Let us recall that the classical expression for the electromagnetic four current is
given by
j = (, v )
(4.135)
where is the charge density, and v the velocity eld. In the case of a pointlike particle which follows the world line describe in a parametric form by the four
functions x ( ), with an arbitrary line parameter, the charge density at the time
t is localized at the position x( ), evaluated at the parameter value such that
t = x0 ( ) (see Fig. 4.1). Therefore
(y , t) = e 3 (y x( ))|t=x0 ( )
67
(4.136)
t
t = x0( )
x ( )
Fig. 4.1 - The space-time trajectory of a point-like particle
It follows
dx 3
(y x( ))|y0 =x0 ( )
(4.137)
dx0
This expression can be put in a covariant form, after taking into account the following
relation
)1
(
dx0
0
0
f ( ) 0
d f ( )(y x ( )) =
(4.138)
x ( )=y 0
d
j (y) = e
From this
j (y) = e
dx 4
d
(y x( ))
d
(4.139)
+
dx 4
d
(y x( )) = e
d
d 4 (y x( )) (4.140)
d x
d
The expression vanishes at any space-time point y, except at the end points x().
We recall also that the equations of motion for a free relativistic scalar particle
can be derived by the following action
S = m
x 2 ,
x =
dx
d
(4.141)
= e(E
dt 1 |v |2
(4.142)
We will show that the lagrangian depends on the four-vector potential A and not
and B.
In fact we will verify that the following lagrangian reproduces
on the elds E
68
if
d
i
x 2
x 2 e
d4 yA (y)j (y)
f
i
d A (x( ))x ( )
(4.143)
Using
L
A
= e x
x
x
L
x
= m eA
x
x 2
(4.144)
=0
x
d x
(4.145)
we get
e
A
A
d x
x + m + e x = 0
x
d x 2
x
(4.146)
d x
= e( A A )x
d x 2
(4.147)
Therefore
m
(4.148)
d
U = eF U
ds
(4.149)
(4.150)
Since the denition of the elds in terms of the vector and scalar potential is given
by
= A
0 A ,
=
A
E
B
(4.151)
t
we get
= (F 10 , F 20 , F 30 ), B
= (F 23 , F 31 , F 12 )
E
(4.152)
69
d
v k
dxi
= eFk0 + eFki
= eE k kij B j v i
dt 1 |v |2
dt
(4.153)
(4.156)
0 =
V
=
V
d4 xA ( 2 A A j )
[
d4 x ( A ) A + (A A )
+ ( A ) A (A A ) (A j )
[
1
1
=
d4 x ( A ) A ( A ) A
2
2
V
]
1
1
+
( A ) A + ( A ) A (A j ) + surface terms
2
2
1
4
d x(F )F
d4 x(A j ) + surface terms
=
2
V
V
[
]
1
4
=
d x ( F F A j ) + surface terms
(4.157)
4
V
70
We see that the action for an electromagnetic eld interacting with an external
current j is given by (here F must be though as a function of A )
1
4
S=
d x F F
d4 x j A
4 V
V
(4.158)
Notice that the interacting term has the same structure we found for the point-like
particle.
We stress again that the A s are the canonical variables of the electrodynamics.
In principle, one could reintroduce the elds by inverting the relations between elds
and potentials. However, in this way, one would end up with a non-local action.
From these considerations one can argue that the potentials play an important
role in quantum mechanics, much more than in the classical case, where they are
essentially a convenient trick. Recall also that the canonical variables satisfy local
commutation relations (the commutator vanishes at space-like distances), implying
that local observables should be local functions of the potentials. This is going to
create us some problem, because the theory is invariant under gauge transformations,
whereas the potentials are not
A (x) A (x) + (x)
(4.159)
(4.160)
Adding a four divergence to the lagrangian density does not change the equations
of motion
t2
t2
dt d3 x =
dt
d3 x0
(4.161)
t
t1
t1
Therefore the invariance of the lagrangian under gauge transformations (neglecting
a four divergence) is guaranteed, if the potentials are coupled to a gauge invariant
and conserved current
j = 0
(4.162)
We have shown that this condition is indeed satised for the point-like particle.
In order to derive the general prescription to couple the electromagnetic potentials to a charged particle, let us go back to the action for the point-like particle.
This prescription is known as the minimal substitution. By choosing x0 = in
(4.143), we get
(4.163)
L = m 1 |v |2 e(A0 v A)
71
from which
p =
L
mv
=
+ eA
v
1 |v |2
(4.164)
= p v L = m
|v |2
1 |v |2
+ m 1 |v |2 + e(A0 v A)
+ ev A
m
=
+ eA0
1 |v |2
(4.165)
= m v
,
p eA
1 |v |2
(4.166)
where the quantities in the left hand sides are the same as in the free case. It follows
that we can go from the free case to the interacting one, by the simple substitution
(minimal substitution)
p p eA
(4.167)
In the free case, inverting the relations between momenta and velocities
|v |2 =
|p|2
,
m2 + |p|2
1 |v |2 =
m2
m2 + |p|2
(4.168)
Hfree =
m2 + |p|2
(4.169)
from which
H = eA0 +
2
m2 + (p eA)
(4.170)
2
m2 + (p eA)
(4.171)
which is nothing but eq. (4.165), after using eq. (4.164). From the point of view
of canonical quantization, the minimal substitution corresponds to the following
substitution in the derivatives
+ ieA
72
(4.172)
4.6
In order to understand better the role of the spin in the Dirac equation we will study
now the non relativistic limit in presence of an electromagnetic eld.
(i m)(x) = 0 = (i eA m)(x) = 0
(4.173)
Notice that the Dirac equation is invariant under the transformation (4.159)
A (x) A (x) + (x)
(4.174)
if we perform also the following local phase transformation on the wave function
(x) eie(x) (x)
(4.175)
Also, eq. (4.173) is invariant under Lorentz transformations, if in going from the
frame S to the frame S (x x = x), the eld A is transformed as
A (x) A (x ) = (1 ) A (x)
(4.176)
=
= (1 )
x
x
x x
(4.177)
Eq. (4.176) says simply that A transforms as a four vector under Lorentz transformations.
In order to study the non relativistic limit is better to write (x) in the following
form
[
]
(x)
(x) =
(4.178)
(x)
where (x)
and (x)
= p eA
(4.179)
and using the representation in 22 blocks of the Dirac matrices given in eq. (4.11),
we get
i
[
]
[
][
] [
][
]
[
]
(x)
0
(x)
m
0
(x)
(x)
=
+
+ eA0
(4.180)
0
(x)
0 m (x)
(x)
t (x)
73
(4.181)
In the non relativistic limit, and for weak elds, the mass term is the dominant
one, and the energy positive solution will behave roughly as eimt . With this
consideration in mind we put
[
]
(x)
(x)
= eimt
(x)
(x)
(4.182)
= + eA0
t
i
= (2m eA0 )
t
i
(4.183)
from which
2m
(4.184)
]
( )2
i
=
+ eA0
t
2m
(4.185)
Aj
Ai
ie
xi
xj
(4.186)
f (x)
xi
(4.187)
[ i , j ] = ieijk B k
(4.188)
and
)
1
1
[i , j ] + [i , j ]+ i j
2
2
1
1
= | |2 + [i , j ] i j = | |2 + [i , j ][ i , j ]
2
4
i
= | |2 + ijk k (ie)ijl B l
2
( )2 = i j i j =
(4.189)
Finally
( )2 = | |2 e B
74
(4.190)
2
(p eA)
e
+ eA0
i
=
B
t
2m
2m
(4.191)
This is nothing but the Pauli equation for an electron interacting with an electromagnetic eld. In particular, the term proportional to the magnetic eld represents
the interaction with a magnetic dipole given by
e
e
=
= S
(4.192)
2m
m
= /2. This shows that the Dirac
where we have introduced the spin matrices S
equation predicts a gyromagnetic ratio equal to two. We may see this also in a
slightly dierent way, by considering the interaction with an uniform and weak
magnetic eld. In this case the vector potential is given by
= 1B
x
A
2
(4.193)
(4.194)
A
=0
[pi , Ai ] = i
(4.195)
it follows
2 |p|2 2ep A
= |p|2 ep (B
x)
(p eA)
= |p|2 epi ijk B j xk = |p|2 ekij xk pi B j
= |p|2 eL
B
= |p|2 e(x p) B
and nally
(4.196)
|p|2
e
B
+ eA0
i
=
(L + 2S)
t
2m 2m
(4.197)
which shows explicitly the value of the gyromagnetic ratio. Experimentally this is
very close to two, and we shall see, in the following, that the dierence is explained by
the quantum electrodynamics (QED). This is in fact, one of the biggest successes of
this theory. However, let us notice that, from the point of view of the Dirac equation,
to nd a value of the gyromagnetic ratio so close to the experimental value is not a
real prediction. In fact, one could think to add to the theory a further interaction
. This term is both Lorentz and gauge invariant. It is
term of the kind F
also possible to show that such a term gives a contribution to the magnetic moment
of the electron, and therefore it changes the gyromagnetic ratio. We shall see that
the requirement that QED is a nite theory forbids, in fact, the appearance of such
a term.
75
4.7
Dirac equation had a great success in explaining the ne structure of the hydrogen atom, but the problem of negative energy solutions, that, in principle, make
the theory unstable, was still there. Dirac looked for a solution to this problem
by taking advantage of the exclusion Pauli principle, which applies to half-integer
spin particles. Dirac made the hypothesis that all the negative energy states were
occupied by electrons. In such a situation, the Pauli principle forbids to any electron in a positive energy state to make a transition to a negative energy state. This
solves the stability problem, but at the same time new phenomena may happen.
For, instance, an electron in a negative energy state could get enough energy (bigger
than 2m which is the minimal energy gap between the negative and positive energy
states) to make a transition to a state of positive energy. If we imagine that in the
state of energy E are present N electrons (we are simplifying things, because due
to the momentum degeneracy there actually an innite number of electrons), and
that one of these electrons undergoes the transition, the energy of the state changes
as follows
E N E E (N 1)E = E N E + E
(4.198)
where E is the energy of all the other electrons (with energy dierent from E) in
the fundamental state. Notice that in the Dirac theory the fundamental state is the
one with all the negative energy states occupied and zero electrons in the energy
positive states. In a sense this is the physical explanation of the innite energy of
the vacuum that we found in the case of the Klein-Gordon eld, and one nds also
for the Dirac case (see later). In a complete analogous way, also the charge of the
vacuum is innite and its variation in the previous transition is given by
Q + N e Q + (N 1)e = Q + N e e
(4.199)
where e is the charge of the electron (e < 0). We see that the vacuum energy and
the charge increase respectively by E and e in the transition. We can interpret
this by saying that the hole left in the vacuum by the electron has charge e and
energy E. That is we can think to the hole being a particle of positive energy and
positive charge. This is the way in which the idea of antiparticles came around.
That the hole is thought as the antiparticle of the electron. The transition of an
electron of negative energy to a state of positive energy is then seen as the creation
of particle antiparticle (the hole) pair. In the same way, once we have a hole in the
vacuum, it may happen that a positive energy electron makes a transition to the
hole state. In this case both the electron and the hole disappears. This is the pair
annihilation phenomenon. Of course this happens with some energy released, that
usually is under electromagnetic form.
The hole theory is nowadays reinterpreted in terms of antiparticles, but this way
of thinking has been extremely fruitful in many elds, as the study of electrons in
metals, in nuclear physics and so on.
76
Assuming seriously the hole theory means that the Dirac equation should admit,
beyond the positive energy solutions corresponding to an electron, other positive
energy solutions with the same mass of the electron, but with opposite charge. To
see this point in a formal way, we look for a transformation of the electron wave
function, (x), to the antielectron (positron) wave function C (x), such that, if
satises
(i eA m)(x) = 0
(4.200)
then C satises
(i + eA m) C (x) = 0
(4.201)
We will require that the transformation is a local one, and that the transformed of the
antiparticle wave function gives back, except for a possible phase factor, the electron
wave function. To build up C we will start by taking the complex conjugate of .
This is clearly the only possibility to change a negative energy solution, described
by eiEt , in a positive energy solution, described by eiEt . By taking the hermitian
conjugate, multiplying by 0 and transposing, we get
eA m = 0
(i eA m)(x) = 0 i
[ T (i eA ) m]T = 0
(4.202)
T = 0 T
(4.203)
where
If there exists a matrix, C, such that
C T C 1 =
(4.204)
(4.205)
This describes a particle with charge e, therefore, apart a phase factor C , we can
identify C with C T :
C = C C T
(4.206)
In the representation where 0 is diagonal we have
0 T = 0 ,
1 T = 1 ,
3 T = 3
2 T = 2 ,
(4.207)
C = i =
In this way, C satises
0
i2
i2
0
C = C 1 = C T = C
77
(4.208)
(4.209)
To understand how the transformation works let us consider, in the rest frame, a
negative energy solution with spin down
()
0
= eimt
0
(4.210)
1
Then
()C = C C T = C C0 () = C i 2 ()
0
0 0
0 1
0
0
0
1
0
= C eimt
0 1
0 00
1
1 0
0 0
and
(4.211)
()C
down
(+)
= C eimt
= C up
0
(4.212)
0
That is, given an energy negative wave function describing an electron with spin
down, its charge conjugated is a positive energy wave function describing a positron
with spin up. For an arbitrary solution with dened energy and spin, by using the
projectors of Section 4.3, we write
=
p + m 1 + 5 n
2m
2
(4.213)
where p0 > 0. = 1 selects the energy sign. Noticing that C commutes with 5 ,
that 5 = 5 , and that
0 0 = T
(4.214)
following from
0 0 =
(4.215)
we obtain
pT + m 1 5 n
p + m 1 + 5 n
T 0
= C C
2m
2
2m
2
p + m 1 + 5 n
C
=
(4.216)
2m
2
C = C C0
We see that C is described by the same four vectors p and n appearing in , but
with opposite sign of the energy. Then
u(p, n) = C v C (p, n),
v(p, n) = C uC (p, n)
78
(4.217)
Since the spin projector selects the states of spin 1/2 along n according to the
sign of the energy, it follows that the charge conjugation inverts the spin projection
of the particle. Notice also that, being C a solution of the Dirac equation with
e e, it follows that the following transformation
C ,
A A
(4.218)
(4.219)
eA)
+ 0m
H = eA0 + 0 (i
(4.220)
i
with
t = t
(4.221)
1
K (t ) = HK 1 (t )
t
(4.222)
K(i)K 1 (t ) = KHK 1 (t )
t
(4.223)
K(i)K 1 = i;
79
KHK 1 = H
(4.224)
KHK 1 = H
(4.225)
The second possibility can be excluded immediately, since under time reversal we
have
,
A,
A
A0 A0
(4.226)
As it follows recalling that the vector potential is generated by a distribution of currents (changing sign under time reversal), whereas the scalar potential is generated
by a distribution of charges. Let us put
K = T (4 4 matrix) (complexconjugation)
(4.227)
(t ) = T H T 1 (t )
(4.228)
T H T 1 = H
(4.229)
with
(4.230)
T 0 T 1 = 0
(4.231)
T 1 T 1 = 1 ,
T 2 T 1 = 2 ,
T 3 T 1 = 3 (4.233)
(4.234)
T2 = 1
(4.235)
To understand the correspondence with the classical results, where momentum and
angular momentum change sign under time reversal, let us study how a positive
energy solution transforms:
[
p + m 1 + 5 n
(t)
K
2m
2
p + m 1 + 5 n
= T
(t)
2m
2
[
]
[
]
p + m 1
1 + 5 n
= T
T T
T 1 T (t)
2m
2
p + m 1 + 5 n
=
(t )
(4.236)
2m
2
80
n
= (n0 , n)
(4.237)
The three discrete symmetry operations described so far, parity, P , charge conjugation, C, and time reversal, (T ), can be combined together into a symmetry
transformation called P CT . Omitting all the phases, this transformation is rather
simple
T
p + m 1 + 5 n
(x)
2m
2
p + m 1 5 n
p + m 1 5 n
=
(i5 (x)) =
P CT (x) (4.239)
2m
2
2m
2
P CT (x) = i5
Comparison with eq. (4.216), giving the charge conjugated of an energy negative
state
p + m 1 + 5 n
C
C =
(4.240)
2m
2
we see that the two expressions dier only in the spin direction. Similar conclusion
can be reached by starting from the Dirac equation multiplied by i5 . We get
(x = x, and A (x ) = A(x))
i5 (ix eA(x)
m)(x) = (ix + eA(x)
m)P CT (x )
) m)P CT (x )
= (ix + eA(x
(4.241)
showing that a positron moving backward in time satises the same equation as
an electron moving forward. Eq. (4.241) tells us that the P CT transformation
on , combined with the P CT transformation on the four-vector potential, that is
A (x) A (x), is a symmetry of the theory.
The interpretation of the positrons as negative energy electrons moving backward
in time is the basis of the positron theory formulated by St
uckelberg and Feynman.
In this approach it is possible to formulate the scattering theory without using eld
theory. In fact, the pair creation and pair annihilation processes can be reinterpreted
in terms of scattering processes among electrons moving forward and backward in
time.
4.8
In this Section we abandon the study of the Dirac wave equation thought as a generalization of the Scrhodinger equation, due to its diculties to cope with many
81
particle states. We will adopt here the point of view of the quantum eld theory. That is the relativistic equation is the equation describing the eld operator.
However, we have shown in the Klein-Gordon case, that after quantization we get
a many particle system satisfying Bose-Einstein statistics. On the other hand we
have seen that the Dirac equation describes spin 1/2 particles, which should satisfy
the Fermi-Dirac statistics. It is quite clear that we will run into troubles insisting in
quantizing the Dirac eld as we did for the Klein-Gordon case. However, in order to
understand the problems and the way to deal with them, we will follow the canonical
way of quantization, showing that this leads to problems with the positivity of the
energy. Looking for a solution of this problem we will nd also the solution to the
problem of the wrong statistics.
We will begin our study by looking for the action giving rise to the Dirac equation.
We will take the quantities and as independent ones. Following the usual
procedure, we multiply the Dirac equation by (in such a way to form a Lorentz
scalar) and integrate over the space-time volume V
0=
V
m) =
d4 x (i
m)
d4 x (i
(4.242)
S=
V
m)
d4 x (i
(4.243)
It is simply veried that this action gives rise to the correct equation of motion for
In fact,
.
L
L
= m,
= i
(4.244)
,
from which
=0
m i
(4.245)
The canonical momenta result to be
L
=
= i ,
L
=0
(4.246)
The canonical momenta do not depend on the velocities. In principle, this creates a
problem in going to the hamiltonian formalism. In fact a rigorous treatment requires
an extension of the classical hamiltonian treatment, which was performed by Dirac
himself. In this particular case, the result one gets is the same as proceeding in
a naive way. For this reason we will avoid to describe this extension, and we will
proceed as in the standard case. Then the hamiltonian density turns out to be
0 0 + i k k m) = (i
+ m)
H = L = i (i
(4.247)
If one makes use of the Dirac equation, it is possible to write the hamiltonian density
as
(4.248)
H = i
t
82
Contrarily to the Klein-Gordon case (see eq. (3.100), the hamiltonian density is
not positive denite. Let us recall the general expression for the energy momentum
tensor (see eq. 3.84)
L i
T =
g L
(4.249)
i, ,
In our case we get
, g ((i
m))
T = i
(4.250)
(4.251)
We verify immediately that this expression has vanishing four divergence. Also
Tk0 = i k
(4.252)
Pk =
d3 x T 0k = P = i
d3 x
(4.253)
In analogous way, by using the general expression for the angular momentum density
(see eq. (3.96))
L
j
M = x T x T i ij
(4.254)
,
The matrices ij
are dened in terms of the transformation properties of the
eld(see eq. (3.94))
1
j
i = ij
(4.255)
2
In our case from eq. (4.32), and from eq. (4.43) in the case of an innitesimal
Lorentz transformation, we get
i
(x) = (x ) (x) = [S() 1](x) (x)
4
from which
(4.256)
i
1
= = [ , ]
2
4
(4.257)
x x + 1 [ , ] (4.258)
x x i = i
= i
2
4
+ 1 12
d3 x ix
2
83
(4.259)
where 12 is the identity matrix in 2 dimensions, and using eq. (4.61) we have dened
12 =
0
0
(4.260)
d3 p
(2)3
(x) =
d3 p
(2)3
]
m [
b(p, n)u(p, n)eipx + d (p, n)v(p, n)eipx
Ep
(4.261)
]
m [
d(p, n)
v (p, n)eipx + b (p, n)
u(p, n)eipx 0 (4.262)
Ep
where Ep = |p|2 + m2 . We will collect here the various properties of the previous
spinors, which are nothing but a trivial extension of the case in which, in the rest
frame, the spin is quantized along the z axis:
Dirac equation
(
p m)u(p, n) = u(p, n)(
p m) = 0
(
p + m)v(p, n) = v(p, n)(
p + m) = 0
(4.263)
u(p, n)u(p, n ) =
v (p, n)v(p, n ) = nn
Ep
u (p, n)u(p, n ) = v (p, n)v(p, n ) =
nn
m
v(p, n)u(p, n ) = v (p, n)u(
p, n ) = 0
(4.264)
Orthogonality
p + m
2m
n
p m
v(p, n)
v (p, n) =
2m
n
u(p, n)
u(p, n) =
84
(4.265)
We are now in the position to express the hamiltonian in terms of the operators
b(p, n) e d(p, n). Using eq. (4.248) we nd
H =
dx
d3 p
m d 3 p
Ep (2)3
m
Ep
Ep
(2)3
]
d(p, n)v (p, n)eipx + b (p, n)u (p, n)eipx
n,n
n,n
Ep m
d3 p d 3 p
Ep Ep
d(p, n)b(
p, n )ei(Ep + Ep )t v (p, n)u(
p, n ) 3 (p + p )
+b (p, n)b(p, n )e+i(Ep Ep )t u (p, n)u(p, n ) 3 (p p )
d(p, n)d (p, n )ei(Ep Ep )t v (p, n)v(p, n ) 3 (p p )
b (p, n)d(
p, n )e+i(Ep + Ep )t u (p, n)v(
p, n ) 3 (p + p )
(4.266)
Performing one of the integrations and using the orthogonality relations, we get
H=
(4.267)
(4.268)
(4.269)
That is, in the Dirac vacuum, the operator d(p, n) behaves as an annihilation operator (thing that we have anticipated in writing). Since the Dirac vacuum is obtained
by the original vacuum acting with d(p, n), it follows that these operators should
satisfy the following algebraic identity
(d(p, n))2 = 0
85
(4.270)
We can satisfy this relation in a uniform algebraic way by requiring that the operators d(p, n) anticommute among themselves
[d(p, n), d(p , n )]+ = 0
(4.271)
This leads to the Jordan and Wigner idea of quantizing the Dirac eld in terms of
anticommutators
[
b(p, n), b (p , n )
= d(p, n), d (p , n )
= nn 3 (p p )
(4.272)
The problem of positivity is then solved automatically, since the four momentum
operator can be written as
P =
] ]
+
(4.273)
Due to the anticommutation relations, the last term turns out to be an innite
negative constant, which, physically, can be associated to the energy of the innite
electrons lling up the Dirac vacuum (called also the Dirac sea). If we ignore this
constant (as we did in the Klein-Gordon case, and with the same warnings), the
energy operator is positive denite. The use of the anticommutators solves also the
problem of the wrong statistics. In fact, the wave functions are now antisymmetric
in the exchange of two Dirac particles (from now on we will put |0Dirac = |0):
b (p1 , n1 )b (p2 , n2 )|0 = b (p2 , n2 )b (p1 , n1 )|0
(4.274)
Therefore, the quanta of the Dirac eld satisfy the Fermi-Dirac statistics. Once we
have realized all that, we can safely forget about the hole theory and related stu.
In fact, looking at the four momentum operator, we can simply say that d (p, n)
creates and d(p, n) annihilates a positron state. then we think to the vacuum as a
state with no electrons and/or positrons (that is without electrons and /or holes).
In the Klein-Gordon case we interpreted the conserved current as the electromag , has
netic current. We shall show now that in the Dirac case the expression
the same interpretation. Let us start evaluating the spatial integral of the density
d3 x
(4.275)
d3 x =
(4.276)
S=
V
eA m)
d4 x(i
86
(4.277)
(4.278)
This forces us to say that the integral of the fourth component of the current should
be the charge operator, and as such it should not be positive denite. In fact, we
nd
Q = e
=
d3 x
[
] ]
The subtraction of the innite charge associated to the Dirac sea leaves us with
an operator which is not anymore positive denite. We see also that the operators
b create particles of charge e (electrons) whereas d create particles of charge e
(positrons). Notice that the interpretation of Q as the charge operator would not
have worked by using commutation relations.
A further potential problem is connected with the causality. We have seen in
the Klein-Gordon case that the causality properties is guaranteed, for the local
observable, by the canonical commutation relations for the elds. But this is just
the property we have given up in the Dirac case. In order to discuss this point, let
us start evaluating the equal time anticommutator for the Dirac eld
[
(x, t), (y , t)
n,n
d3 p
(2)3
m d3 p
Ep
(2)3
m
nn 3 (p p )
Ep
u(p, n)
u(p, n)0 eip (x y ) + v(p, n)
v (p, n)0 eip (x y )
[(
]
]
p m ip (x y )
p + m ip (x y )
e
+
e
0
=
2m
2m
) (
)]
[(
d3 p m
p + m
p m
=
+
0 eip (x y )
(2)3 Ep
2m
2m
d3 p m 2Ep ip (x y )
=
e
= 3 (x y )
(2)3 Ep 2m
d3 p m
(2)3 Ep
(4.280)
]
+
=0
(4.281)
(4.282)
This shows that also in the Dirac case one can use the canonical formalism, but with
anticommutators in place of the commutators. For arbitrary space-time separations
87
we get
[
[(
(x), (y)
p + m ip(x y)
p m ip(x y)
=
e
+
e
0
2m
2m
[
]
[(
)
]
d3 p 1
ip(x y) eip(x y)
=
i + m 0
e
x
(2)3 2Ep
[(
)
]
=
i + m 0 i(x y)
(4.283)
d3 p m
(2)3 Ep
where (x) is the invariant function dened in eq. (3.148), in the evaluation of
commutator for he Klein-Gordon eld. From the properties of the (x) function, it
follows that the anticommutator of the Dirac elds vanishes at space-like distances.
It follows that the Dirac eld cannot be an observable quantity. In fact one could do
more, by evaluating the commutator of the Dirac eld. It is easy to show that the
result does not vanish at space-like distances. This by itself would put in a serious
trouble the idea of quantizing Dirac eld via commutation relations. But how do
we solve the causality problem. The crucial observation is in the following identity
[AB, C] = A[B, C] + [A, C]B = A[B, C]+ [A, C]+ B
(4.284)
which holds for arbitrary operators. The identity shows that AB commutes with
C if A and B separately commute or anticommute with C. An immediate consequence is that a local quantity containing an even number of Dirac elds commutes
with itself at space-like distances. So, in order to reconcile the causality with the
quantization of the Dirac eld we have to give up with its property of being an
observable. However, all the important physical quantities, as energy-momentum
tensor and electromagnetic current are bilinear in the Fermi elds, and therefore
they are observable quantities.
What we have shown here is that, in order to give a sense to the quantization
of the Dirac eld, we have been forced to use anticommutation relations, which, in
turn, imply that the corresponding quanta obey the Fermi-Dirac statistics. This
is nothing but an example of the celebrated spin statistics theorem that was
proved by Pauli in 1940. This theorem says that in a Lorentz invariant local eld
theory, integer and half-integer particles must satisfy respectively Bose-Einstein and
Fermi-Dirac statistics.
88
Chapter 5
The electromagnetic field
5.1
In Section 4.5 we have shown that the action for the electromagnetic eld must be
expressed in terms of the four-vector potential A . We recall that the lagrangian
density for the free case is given by (see. eq. (4.158)
1
L = F F
4
(5.1)
F = A A
(5.2)
where
The resulting equations of motion are
2 A ( A ) = 0
(5.3)
(5.4)
In fact, A and A satisfy the same equations of motion (and give rise to the same
electromagnetic eld). In other words the action of electrodynamics is gauge invariant. Up to now, we have considered only symmetries depending on a nite number
of parameters. For instance, in the case of the O(2) symmetry for the charged scalar
eld, the transformation symmetry depends on a single parameter, the rotation angle. In the case of the gauge symmetry one deals with a continuous number of
parameters, given by the function (x). In fact, in each space-time point, we can
change the denition of A by adding the gradient of evaluated at that point.
The main consequence of this type of invariance is to reduce the eective degrees of
freedom of the theory from 4 to 2. Let us start from the classical theory. We recall
that it is possible to use the gauge invariance to require some particular condition
89
on the eld A . For instance, we can perform a gauge transformation in such a way
that the transformed eld satises
A = 0
(5.5)
(5.7)
(5.8)
In the Lorentz gauge the equations of motion simplify and reduce to the wave equation, or to the Klein-Gordon equation with m = 0. This, together the fact, that
in the Lorentz gauge the covariance of the theory is explicitly preserved (the gauge
condition is Lorentz invariant), makes this gauge a very popular one. However, the
counting of the eective degrees of freedom is not so evident. From this point of
view a better choice is the Coulomb gauge, which is dened as the gauge where
the scalar potential and the spatial divergence of the vector potential vanish. To see
that such a gauge exists, let us perform the following gauge transformation
A (x)
Clearly
= A (x)
A0 (x, t )dt
A0 = 0
(5.9)
(5.10)
(5.11)
(5.12)
|x|
90
(5.13)
Then
from which
1
d3 x
(x, t) =
A (x , t)
4
|x x |
(5.14)
1
(x, t)
d3 x
=
A (x , t)
t
4
|x x |
(5.15)
(5.16)
= 0
E
=
A
(5.17)
and, in terms of A
Therefore
(x, t)
=0
(5.18)
t
from which A0 = A0 = 0 showing that the second gauge transformation does not
destroy the vanishing of the scalar potential. In conclusion, we have shown that it
is possible to choose a gauge such that
A
=0
A0 =
(5.19)
It follows that the independent degrees of freedom are only two. Unfortunately in
this gauge the explicit Lorentz covariance is lost. Another way of showing that A
has only two degrees of freedom is through the equations of motion. Let us consider
the four dimensional Fourier transform of A (x)
A (x) =
d4 k eikx A (k)
(5.20)
(5.21)
Let us now decompose A (k) in terms of four independent four vectors, which can
be chosen as k = (E, k), k = (E, k), and two further four vectors e (k), = 1, 2,
orthogonal to k
= 1, 2
(5.22)
k e = 0,
The decomposition of A (k) reads
A (k) = a (k)e + b(k)k + c(k)k
(5.23)
(5.24)
The term in b(k) cancels, therefore it is left undetermined by the equations of motion.
For the other quantities we have
k 2 a (k) = c(k) = 0
(5.25)
where
(x) =
d4 k eikx (k)
(5.27)
(5.28)
(5.29)
showing that k A (k) = 0. Therefore the choice b(k) = 0 is equivalent to the choice
of the Lorentz gauge.
Let us consider now the quantization of this theory. If we would like to require the
explicit covariance of the theory we would require non trivial commutation relations
for all the component of the eld. That is
[A (x, t), (y , t)] = ig 3 (x y )
(5.30)
(5.31)
with
=
L
A
(5.32)
To evaluate the conjugated momenta is better to write the lagrangian density (see
eq. (5.1) in the following form
1
1
1
L = [A, A, ][A, A, ] = A, A, + A, A,
4
2
2
Therefore
L
= A, + A, = F
A,
92
(5.33)
(5.34)
implying
=
L
= F 0
A
(5.35)
It follows
L
=0
A 0
We see that it is impossible to satisfy the condition
0 =
(5.36)
(5.37)
(5.38)
(5.39)
(just think to the Klein-Gordon case). We will see that the minus sign is necessary
to recover a positive hamiltonian density. We now express this lagrangian density
in terms of the gauge invariant one, given in eq. (5.1). To this end we observe that
the dierence between the two lagrangian densities is nothing but the second term
of eq. (5.33)
[
1
1
1
A, A, = A, A ( A, )A
2
2
2
[
]
[
]
1
1
=
A, A
( A )A + ( A )2
2
2
2
(5.40)
Then, up to a four divergence, we can write the new lagrangian density in the form
1
1
L = F F ( A )2
4
2
(5.41)
One can check that this form gives the correct equations of motion. In fact from
L
= A, + A, g ( A ),
A,
L
=0
A
(5.42)
we get
0 = 2 A + ( A ) ( A ) = 2 A
93
(5.43)
the term
1
( A )2
(5.44)
2
which is not gauge invariant, is called the gauge fixing term. More generally, we
could add to the original lagrangian density a term of the form
( A )2
2
(5.45)
(5.46)
These equations are the same as the Maxwell equations in the Lorentz gauge. Therefore, in the following we will use = 1. From eq. (5.42) we see that
0 =
L
= A
A 0
(5.47)
(5.48)
The price to pay to quantize the theory in a covariant way is to work in a Hilbert
space much bigger than the physical one. The physical states span a subspace which
is dened by the previous relation. A further bonus is that in this way one has to
do with local commutation relations. On the contrary, in the Coulomb gauge, one
needs to introduce non local commutation relations for the canonical variables. We
will come back later to the condition (5.48).
Since we dont have to worry any more about the operator condition 0 = 0, we
can proceed with our program of canonical quantization. The canonical momentum
densities are
L
= F 0 g 0 ( A )
(5.49)
=
A
or, explicitly
A
0 = A = A 0
i = i A0 0 Ai = A i + i A0
(5.50)
Since the spatial gradient of the eld commutes with the eld itself at equal time,
the canonical commutator (5.30) gives rise to
[A (x, t), A (y , t)] = ig 3 (x y )
94
(5.51)
To get the quanta of the eld we look for plane wave solutions of the wave equation.
We need four independent four vectors in order to expand the solutions in the
momentum space. In a given frame, let us consider the unit four vector which
denes the time axis. This must be a time-like vector, n2 = 1, and we will choose
n0 > 0. For instance, n = (1, 0, 0, 0). Then we take two four vectors ()
, = 1, 2,in
2
the plane orthogonal to n and k . Notice that now k = 0,since we are considering
solutions of the wave equation. Therefore
()
k ()
= n = 0,
= 1, 2
(5.52)
( )
()
=
(5.53)
Next, we dene a unit space-like four vector, orthogonal to n and lying in the plane
(k, n)
.
n (3) = 0
(5.54)
with
(3)
(3)
= 1
(5.55)
()
By construction (3)
is orthogonal to . This four vector is completely xed by
the previous conditions, and we get
k (n k)n
(n k)
(5.56)
(0)
= n
(5.57)
(3)
=
A last unit four vector we choose n
( )
()
= g
(5.58)
( )
()
g = g
(5.59)
(3) = (0, 0, 0, 1)
(5.60)
A (x) =
d3 k
2k (2)3
()
(k)
=0
95
a (k)eikx + a (k)eikx
(5.61)
where we have included the hermiticity condition for A (x). For any xed , this
expansion is the same as the one that we wrote for the Klein-Gordon eld, with the
substitution ()
a (k) a(k). Then, from eq. (3.60)
()
(k)a (k) = i
d3 x fk (x)t A (x)
()
(5.62)
with the functions fk (x) dened as in Section 3.2. Using the orthogonality of the
( ) s we nd
()
a (k) = ig
d3 x ( ) (k)fk (x)t A (x)
(5.63)
and analogously
a (k) = ig
()
(5.64)
d3 x d3 y fk (x)fk (y) ig g (
(k)(
(k )g 3 (x y )
)]
fk (x)fk (y) ig g ( ) (k)( ) (k )g 3 (x y )
d3 x fk (x)it fk (x)g
()
(5.65)
Analogously
[a (k), a (k )] = g 3 (k k )
(5.66)
[a (k), a (k )] = [a (k), a (k )] = 0
(5.67)
(5.68)
with the invariant function (x) dened as in eq. (3.148), but with m = 0. The
commutation rules we have derived for the operators a (k) create some problem.
Let us consider a one-particle state
|1, =
d3 k f (k)a (k)|0
(5.69)
= g
d3 k |f (k)|2
96
(5.70)
Therefore the states with = 0 have negative norm. This problem does not come
out completely unexpected. In fact, our expectation is that only the transverse
states ( = 1, 2), are physical states. For the moment being we have ignored the
gauge xing condition phys| A |phys = 0, but its meaning is that only part of
the total Hilbert space is physical. Therefore the relevant thing is to show that the
states satisfying the Lorentz condition have positive norm. To discuss the gauge
xing condition, let us notice that formulated in the way we did, being bilinear in
the states, it could destroy the linearity of the Hilbert space. So we will try to
modify the condition in a linear one
A |phys. = 0
(5.71)
But this would be a too strong requirement. Not even the vacuum state satises it.
However, if we consider the positive and negative frequency parts of the eld
A(+)
(x)
d3 k
2k (2)3
ikx ,
()
(k)a (k)e
(+)
A()
(x)) (5.72)
(x) = (A
=0
(5.73)
(5.74)
To make this condition more explicit let us evaluate the four divergence of A(+)
A(+)
(x)
d3 k
2k (2)3
eikx
k ()
(k)a (k)
(5.75)
=0,3
k (0)
= (n k)
(5.76)
from which
[a0 (k) a3 (k)]|phys. = 0
(5.77)
Notice that
[a0 (k) a3 (k), a0 (k ) a3 (k )] = 3 (k k ) + 3 (k k ) = 0
(5.78)
Let us denote by k (n0 , n3 ) the state with n0 scalar photons (that is with polarization
(3)
(0)
(k)), and with n3 longitudinal photons (that is with polarization (k)). Then
the following states satisfy the condition (5.73)
(m)
1
(a0 (k) a3 (k))m k (0, 0)
m!
97
(5.79)
||k ||2 = 0
(5.80)
More generally we can make the following observation. Let us consider the number
operator for scalar and longitudinal photons
N=
(5.81)
Notice the minus sign that is a consequence of the commutation relations, and it
ensures that N has positive eigenvalues. For instance
N a0 (k)|0 =
(5.82)
Let us consider a physical state with a total number n of scalar and longitudinal
photons. Then
n |N |n = 0
(5.83)
since a0 and a3 act in the same way on a physical state (see eq. (5.77)). It follows
nn |n = 0
(5.84)
Therefore all the physical states with a total denite number of scalar and longitudinal photons have zero norm, except for the vacuum state (n = 0). Then
n |n = n,0
(5.85)
A generic physical state with zero transverse photons is a linear superposition of the
previous states
ci |i
(5.86)
| = c0 |0 +
i=0
(5.87)
The proof that a physical state has a positive norm can be extended to the case in
which also transverse photons are present. Of course, the coecients ci , appearing
in the expression of a physical state, are completely arbitrary, but this is not going
to modify the values of the observables. For instance, consider the hamiltonian, we
have
H =
d3 x : [ A L] :
1
1
d3 x : F 0 A ( A )A 0 + F F + ( A )2 :
4
2
98
(5.88)
One can easily show the hamiltonian is given by the sum of all the degrees of freedom
appearing in A (see the Klein-Gordon case, eq. (3.102))
[
3 (
)
1 3
i )2 A 2 A
2 :
H =
dx :
A 2i + (A
0
0
2
i=1
d3 k k :
a (k)a (k)
a0 (k)a0 (k)
(5.89)
=1
Since on the physical states a0 and a3 act in the same way, we get
phys.|H|phys. = phys.|
d3 k k
a (k)a (k)|phys.
(5.90)
=1
cn cm n |A (x)|m
(5.91)
n,m
Since A change the occupation number by one unit and all the states |n have
zero norm (except for the state with n = 0), the only non vanishing contributions
come from n = 0, m = 1 and n = 1, m = 0
|A (x)| = c0 c1 0|
d3 k
2k
(2)3
(0)
eikx [(3)
(k)a3 (k) + (k)a0 (k)]|1 + c.c.
(5.92)
In order to satisfy the gauge condition the state |1 is of the form
|1 =
(5.93)
(0)
(5.94)
and therefore
|A (x)| =
d3 k
2k
(2)3
99
k
(k n)
(5.95)
from which
|A (x)| = (x)
with
(x) =
d3 k
1
(ic0 c1 eikx f (k) + c.c.)
3
n
k
2k (2)
(5.96)
(5.97)
100
Chapter 6
Symmetries in field theories
6.1
In this Section and in the following we will study, from a classical point of view,
some eld theory with particular symmetry properties. We will start examining the
linear -model. This is a model for N scalar elds, with a symmetry O(N ). The
lagrangian is given by
(
N
N
N
1
1
i i
L=
i i 2
i i
2 i=1
2 i=1
4 i=1
)2
(6.1)
This lagrangian is invariant under linear transformations acting upon the vector
= (1 , , N ) and leaving invariant its norm
2=
||
i i
(6.2)
i=1
(6.3)
(6.4)
=0
(6.5)
i ij j = 0
(6.6)
ij = ji
(6.7)
from which
or, in components,
This is satised by
101
i = Sij j
(6.9)
SS T = 1
(6.10)
implying
In fact, by exponentiating the innitesimal transformation one gets
S = e
(6.11)
T =
(6.12)
with
implying that S is an orthogonal transformation. The matrices S form the rotation
group in N dimensions, O(N ).
The Noethers theorem implies a conserved current for any symmetry of the
theory. In this case we will get N (N 1)/2 conserved quantities. It is useful for
further generalizations to write the innitesimal transformation in the form
i
i = ij j = AB TijAB j ,
2
i, j = 1, , N,
A, B = 1, , N
(6.13)
which is similar to what we did in Section 3.3 when we discussed the Lorentz transformations. By comparison we see that the matrices T AB are given by
TijAB = i(iA jB jA iB )
(6.14)
(6.15)
This is nothing but the Lie algebra of the group O(N ), and the T AB are the innitesimal generators of the group. Applying now the Noethers theorem we nd
the conserved current
j =
i
L
i = i, AB TijAB j
( i )
2
(6.16)
(6.18)
in agreement with (3.169). One can easily check that the charges associated to the
conserved currents close the same Lie algebra as the generators T AB . More generally,
if we have conserved currents given by
jA = ii, TijA j
(6.19)
[T A , T B ] = if ABC T C
(6.20)
with
then, using the canonical commutation relations, we get
[QA , QB ] = if ABC QC
with
QA =
d3 x j0A (x) = i
d3 x i TijA j
(6.21)
(6.22)
(6.23)
(6.24)
where are the Pauli matrices. Noticing that 2 is pure imaginary, 1 and 3 real,
and that 2 anticommutes with 1 and 3 , we get
M = 2 M 2
(6.25)
2 = 2 + | |2 = 1 T r(M M )
||
2
(6.26)
Using this it is easy to write the lagrangian for the -model in the form
)2
1
1
1 (
L = T r( M M ) 2 T r(M M ) T r(M M )
4
4
16
(6.27)
(6.28)
where L and R are two special (that is with determinant equal to 1) unitary matrices,
that is L, R SU (2). The reason to restrict these matrices to be special is that
only in this way the transformed matrix satisfy the condition (6.25). In fact, if A is
a 2 2 matrix with detA = 1, then
2 AT 2 = A1
103
(6.29)
2 M 2 = 2 L M RT 2 = 2 L 2 (2 M 2 )2 RT 2
(6.30)
2 L 2 = 2 L 2 = L
2 RT 2 = R1 = R
=L
(6.31)
since the L and R are independent transformations, the invariance group in this
basis is SU (2)L SU (2)R . In fact this group and O(4) are related by the following
observation: the transformation M LM R is a linear transformation on the
matrix elements of M , but from the relation (6.26) we see that M LM R leaves
= ( , ) invariant and therefore the same must be true
the norm of the vector
for the linear transformation acting upon the matrix elements of M , that is on
and . Therefore this transformation must belongs to O(4). This shows that the
two groups SU (2) SU (2) and O(4) are homomorphic (actually there is a 2 to 1
relationship, since L and R dene the same S as L and R).
We can evaluate the eect of an innitesimal transformation. To this end we
will consider separately left and right transformations. We parameterize the transformations as follows
i
L 1 L ,
2
i
R 1 R
2
(6.32)
then we get
i
i
L M = ( L )M = ( L )( + i ) =
2
2
1
i
L + (L L ) (6.33)
2
2
(6.34)
L M = L + iL
(6.35)
Since
we get
1
L = L ,
2
Analogously we obtain
1
R = R ,
2
1
L = (L L )
2
1
R = (R + R )
2
(6.36)
(6.37)
1
= [(L + R ) (L R )]
2
104
(6.38)
(6.39)
(6.40)
These transformations span a subgroup SU (2) of SU (2)L SU (2)R called the diagonal subgroup. In this case we have
=
= 0,
(6.41)
We see that the transformations corresponding to the diagonal SU (2) are the rotations in the 3-dimensional space spanned by . These rotations dene a subgroup
O(3) of the original symmetry group O(4). From the Noethers theorem we get the
conserved currents
1
1
jL = , L + , (L L )
2
2
(6.42)
JL = , , ,
(6.43)
JR = , + , ,
(6.44)
and dividing by L /2
and analogously
Using the canonical commutation relations one can verify that the corresponding
charges satisfy the Lie algebra of SU (2)L SU (2)R
[QLi , QLj ] = iijk QLk ,
R
R
[QR
i , Qj ] = iijk Qk ,
[QLi , QR
j ] = 0
(6.45)
1
JA = (JL JR )
2
(6.46)
one has
JV = ,
(6.47)
JA = , ,
(6.48)
and
The corresponding algebra of charges is
[QVi , QVj ] = iijk QVk ,
A
[QVi , QA
j ] = iijk Qk ,
A
V
[QA
i , Qj ] = iijk Qk
(6.49)
These equations show that QVi are the innitesimal generators of a subgroup SU (2)
of SU (2)L SU (2)R which is the diagonal subgroup, as it follows from
[QVi , j ] = iijk k ,
105
[QVi , ] = 0
(6.50)
H=
d xH =
3
]
[ N
]
L
1 2
3
2
2
dx
i L = d x
(i + |i | ) + V (|| )
3
i=1
2 i=1
(6.51)
Since in the last member of this relation the rst two terms are positive denite,
it follows that the absolute minimum is obtained for constant eld congurations,
such that
2)
V (||
=0
(6.52)
i
Let us call by vi the generic solution to this equation (in general it could happen that
the absolute minimum is degenerate). Then the condition for getting a minimum is
that the eigenvalues of the matrix of the second derivatives of the potential at the
stationary point are denite positive. In this case we dene new elds by shifting
the original elds by
i i = i vi
(6.53)
The lagrangian density becomes
1
+ v |2 )
L = i i V (|
2
(6.54)
1 2 V
i j +
v
2 i j =
(6.55)
This equation shows that the particle masses are given by the eigenvalues of the
second derivative of the potential at the minimum. In the case of the linear -model
we have
1 2 2 2
V = 2 ||
+ (|| )
(6.56)
2
4
Therefore
V
2
= 2 i + i ||
(6.57)
i
In order to have a solution to the stationary condition we must have i = 0, or
2=
||
(6.58)
This equation has real solutions only if 2 / < 0. However, in order to have a
potential bounded from below one has to require > 0, therefore we may have non
106
zero solutions to the minimum condition only if 2 < 0. But, notice that in this
case, 2 cannot be identied with a physical mass, these are given by the eigenvalues
of the matrix of the second derivatives of the potential at the minimum and they
are positive denite by denition. We will study this case in the following Sections.
In the case of 2 > 0 the minimum is given by i = 0 and one can study the
2 )2 as a small perturbation (that is requiring that
theory by taking the term (||
both and the values of i , the uctuations, are small). The free theory is given
by the quadratic terms in the lagrangian density, and they describe N particles of
common mass m. Furthermore, both the free and the interacting theories are O(N )
symmetric.
6.2
In this Section we will see that the linear -model with 2 < 0, is just but an example
of a general phenomenon which goes under the name of spontaneous symmetry
breaking of symmetry. This phenomenon lies at the basis of the modern description
of phase transitions and it has acquired a capital relevance in the last years in all
eld of physics. The idea is very simple and consists in the observation that a
theory with hamiltonian invariant under a symmetry group may not show explicitly
the symmetry at the level of the solutions. As we shall see this may happen when
the following conditions are realized:
The theory is invariant under a symmetry group G.
The fundamental state of the theory is degenerate and transforms in a non
trivial way under the symmetry group.
Just as an example consider a scalar eld described by a lagrangian invariant under
parity
P :
(6.59)
The lagrangian density will be of the type
1
L = V (2 )
2
(6.60)
If the vacuum state is non degenerate, barring a phase factor, we must have
P |0 = |0
(6.61)
(6.62)
0||0 = 0
(6.63)
from which
107
Things change if the fundamental state is degenerate. This would be the case in the
example (6.60), if
2
(6.64)
V (2 ) = 2 + 4
2
4
with 2 < 0. In fact, this potential has two minima located at
= v,
v=
(6.65)
By denoting with |R e |L the two states corresponding to the classical congurations = v, we have
P |R = |L = |R
(6.66)
Therefore
(6.67)
which now does not imply that the expectation value of the eld vanishes. In the
following we will be rather interested in the case of continuous symmetries. So let
us consider two scalar elds, and a lagrangian density with symmetry O(2)
1 1 2 2
L =
( )
2
2
4
(6.68)
= 2 + 2
1
2
(6.69)
where
= 0,
For 2 > 0 there is a unique fundamental state (minimum of the potential)
whereas for 2 < 0 there are innite degenerate states given by
2 = 2 + 2 = v 2
||
1
2
(6.70)
with v dened as in (6.65). By denoting with R() the operator rotating the elds
in the plane (1 , 2 ), in the non-degenerate case we have
and
R()|0 = |0
(6.71)
(6.72)
(6.73)
where | is one of the innitely many degenerate fundamental states lying on the
2 = v 2 . Then
circle ||
0|i |0 = 0|R1 ()R()i R1 ()R()|0 = |i |
108
(6.74)
with
i = R()i R1 () = i
(6.75)
Again, the expectation value of the eld (contrarily to the non-degenerate state)
does not need to vanish. The situation can be described qualitatively saying that the
existence of a degenerate fundamental state forces the system to choose one of these
equivalent states, and consequently to break the symmetry. But the breaking is only
at the level of the solutions, the lagrangian and the equations of motion preserve
the symmetry. One can easily construct classical systems exhibiting spontaneous
symmetry breaking. For instance, a classical particle in a double-well potential.
This system has parity invariance x x, where x is the particle position. The
equilibrium positions are around the minima positions, x0 . If we put the particle
close to x0 , it will perform oscillations around that point and the original symmetry
is lost. A further example is given by a ferromagnet which has an hamiltonian
invariant under rotations, but below the Curie temperature exhibits spontaneous
magnetization, breaking in this way the symmetry. These situations are typical for
the so called second order phase transitions. One can describe them through the
Landau free-energy, which depends on two dierent kind of parameters:
Control parameters, as 2 for the scalar eld, and the temperature for the
ferromagnet.
Order parameters, as the expectation value of the scalar eld or as the
magnetization.
The system goes from one phase to another varying the control parameters, and the
phase transition is characterized by the order parameters which assume dierent
values in dierent phases. In the previous examples, the order parameters were zero
in the symmetric phase and dierent from zero in the broken phase.
The situation is slightly more involved at the quantum level, since spontaneous
symmetry breaking cannot happen in nite systems. This follows from the existence
of the tunnel eect. Let us consider again a particle in a double-well potential, and
recall that we have dened the fundamental states through the correspondence with
the classical minima
x = x0 |R
x = x0 |L
(6.76)
But the tunnel eect gives rise to a transition between these two states and as a
consequence it removes the degeneracy. In fact, due to the transition the hamiltonian
acquires a non zero matrix element between the states |R and |L. By denoting
with H the matrix of the hamiltonian between these two states, we get
[
H= 0
1
109
1
0
(6.77)
(6.78)
(6.79)
(6.80)
with eigenvalue EA = 0 1 . One can show that 1 < 0 and therefore the fundamental state is the symmetric one, |S. This situation gives rise to the well known
eect of quantum oscillations. We can express the states |R and |L in terms of
the energy eigenstates
|R =
|L =
1 (|S + |A)
2
1 (|S |A)
2
(6.81)
Let us now prepare a state, at t = 0, by putting the particle in the right minimum.
This is not an energy eigenstate and its time evolution is given by
(
1
1
|R, t = eiES t |S + eiEA t |A = eiES t |S + eitE |A
2
2
(6.82)
and the energy dierence between the maximum at = 0, and the minimum at
= v, becomes innite in the limit of innite volume
H( = 0) H( = v) =
d3 x
V
4
4
2 2 4
v + v =
d3 x =
V
2
4
4 V
4
110
(6.85)
6.3
From our point of view, the most interesting consequence of spontaneous symmetry
breaking is the Goldstone theorem. This theorem says that for any continuous
symmetry spontaneously broken, there exists a massless particle (the Goldstone
boson). The theorem holds rigorously in a local eld theory, under the following
hypotheses
The spontaneous broken symmetry must be a continuous one.
The theory must be manifestly covariant.
The Hilbert space of the theory must have a denite positive norm.
We will limit ourselves to analyze the theorem in the case of a classical scalar eld
theory. Let us start considering the lagrangian for the linear -model with invariance
O(N )
1 2 2
L =
(6.86)
( )
2
2
4
The conditions that V must satisfy in order to have a minimum are
V
2=0
= 2 l + l ||
l
with solutions
(6.87)
2
(6.88)
The character of the stationary points can be studied by evaluating the second
derivatives
2V
2 ) + 2l m
= lm (2 + ||
(6.89)
l m
l = 0,
2 = v2,
||
v=
111
(6.91)
V
V ()
minimum
1 2 V
(l vlN )(m vmN )
2 l m minimum
(6.92)
If we are going to make a perturbative expansion, the right elds to be used are
l vlN , and their mass is just given by the coecient of the quadratic term
2
=
Mlm
2 V
= 22 lN mN
minimo
l m
0 0
0
0 0
0
0 0 22
(6.93)
(6.95)
N
1
m2
2a + 2 +
2
4
a=1
(N 1
)2
2a
(6.96)
a=1
In this form the original symmetry O(N ) is broken. However a residual symmetry
N 1 2
a , and it is
O(N 1) is left. In fact, V depends only on the combination a=1
invariant under rotations around the axis we have chosen as representative for the
fundamental state, (0, , v). It must be stressed that this is not the most general
potential invariant under O(N 1). In fact the most general potential (up to the
fourth order in the elds) describing N scalar elds with a symmetry O(N 1)
would depend on 7 coupling constants, whereas the one we got depends only on the
two parameters m and . Therefore spontaneous symmetry breaking puts heavy
constraints on the dynamics of the system. We have also seen that we have N 1
massless scalars. Clearly the rotations along the rst N 1 directions leave the
potential invariant, whereas the N 1 rotations on the planes a N move away
from the surface of the minima. This can be seen also in terms of generators. Since
the eld we have chosen as representative of the ground state is i |min = viN , we
have
(6.97)
Tijab j |min = i(ia jb ib ja )vjN = 0
since a, b = 1, , N 1, and
TijaN j |min = i(ia jN iN ja )vjN = ivia = 0
(6.98)
(here O(N 1)), and in the so called broken generators, each of them corresponding
to a massless Goldstone boson. In general, if the original symmetry group G of the
theory is spontaneously broken down to a subgroup H (which is the symmetry of
the vacuum), the Goldstone bosons correspond to the generators of G which are left
after subtracting the generators of H. Intuitively one can understand the origin of
the massless particles noticing that the broken generators allow transitions from a
possible vacuum to another. Since these states are degenerate the operation does
not cost any energy. From the relativistic dispersion relation this implies that we
must have massless particles. One can say that Goldstone bosons correspond to at
directions in the potential.
6.4
Many eld theories possess global symmetries. These are transformations leaving
invariant the action of the system and are characterized by a certain number of
parameters which are independent on the space-time point. As a prototype we can
consider the free Dirac lagrangian
/ m](x)
L0 = (x)[i
(6.99)
(6.100)
If one has more than one eld, Q is a diagonal matrix having as eigenvalues the
charges of the dierent elds measured in unit e. For instance, a term as 2 1 ,
with a scalar eld, is invariant by choosing Q(1 ) = Q() = 1, and Q(2 ) = 2.
This is a so called abelian symmetry since
eiQ eiQ = ei( + )Q = eiQ eiQ
(6.101)
(6.102)
(6.103)
will be invariant as the mass term under the local phase transformation. To construct the covariant derivative, we need to enlarge the eld content of the theory,
by introducing a vector eld, the gauge field A , in the following way
D = + ieQA
(6.104)
1
= eiQ(x) + ieQ(A )
e
from which
1
A = A +
e
(6.105)
(6.106)
(6.107)
is then invariant under gauge transformations, or under the local group U (1). We see
also that by requiring local invariance we reproduce the electromagnetic interaction
as obtained through the minimal substitution we discussed before.
In order to determine the kinetic term for the vector eld A we notice that
eq. (6.102) implies that under a gauge transformation, the covariant derivative
undergoes a unitary transformation
D D = eiQ(x) D eiQ(x)
(6.108)
(6.109)
F = A A
(6.110)
F eiQ(x) F eiQ(x) = F
(6.111)
with
transforms in the same way
The last equality follows from the commutativity of F with the phase factor. The
complete lagrangian density is then
( + ieQA ) m] 1 F F
L = L + LA = [i
4
114
(6.112)
The gauge principle has automatically generated an interaction between the gauge
eld and the charged eld. We notice also that gauge invariance prevents any
mass term, 12 M 2 A A . Therefore, the photon eld is massless. Also, since the
local invariance implies the global ones, by using the Noethers theorem we nd the
conserved current as
L
(Q)
j =
=
(6.113)
,
from which, eliminating the innitesimal parameter ,
Q
J =
6.5
(6.114)
The approach of the previous section can be easily extended to local non-abelian
symmetries. We will consider the case of N Dirac elds. The free lagrangian
L0 =
/ m)a
a (i
(6.115)
a=1
(6.116)
U G
(6.117)
where T A denote the generators of the Lie algebra associated to G, Lie(G), (that is
the vector space spanned by the innitesimal generators of the group) in the fermion
representation. The generators T A satisfy the algebra
[T A , T B ] = ifCAB T C
(6.118)
where fCAB are the structure constants of Lie(G). For instance, if G = SU (2), and
we take the fermions in the fundamental representation,
[
1
=
2
we have
TA =
A
,
2
A = 1, 2, 3
115
(6.119)
(6.120)
where A are the Pauli matrices. In the general case the T A s are N N hermitian
matrices that we will choose normalized in such a way that
1
T r(T A T B ) = AB
2
(6.121)
(6.123)
Let us now proceed to the case of the local symmetry by dening again the concept
of covariant derivative
D (x) [D (x)] = U (x)[D (x)]
(6.124)
D = + igB
(6.125)
Dab
= ab + ig(B )ab
(6.126)
(6.127)
(6.128)
and
i
B (x) = U (x)B (x)U 1 (x) + ( U (x))U 1 (x)
g
For an innitesimal transformation
(6.129)
U (x) 1 iA (x)T A
(6.130)
1
B (x) = iA (x)[T A , B (x)] + ( A (x))T A
g
(6.131)
we get
(6.132)
(6.133)
The dierence with respect to the abelian case is that the eld undergoes also a
homogeneous transformation.
The kinetic term for the gauge elds is constructed as in the abelian case. In
fact the quantity
[D , D ] igF
(6.134)
in virtue of the eq. (6.124), transforms as under gauge transformations, that is
= igF
U (x)
([D , D ]) = igF
= U (x)([D , D ]) = U (x) (igF )
(6.135)
This time the tensor F is not invariant but transforms homogeneously, since it
does not commute with the gauge transformation as in the abelian case
F
= U (x)F U 1 (x)
(6.136)
(6.137)
(6.138)
F = ( B B ) + ig[B , B ]
(6.139)
F = FC T C
(6.140)
FC = AC AC gfCAB AA AB
(6.141)
or
in components
with
The main feature of the non-abelian gauge theories is the bilinear term in the previous expression. Such a term comes because fCAB = 0, expressing the fact that G is
not abelian. The kinetic term for the gauge eld, expressed in components, is given
by
1
LA =
FA FA
(6.142)
4 A
Therefore, whereas in the abelian case LA is a free lagrangian (it contains only
quadratic terms), now it contains interaction terms cubic and quartic in the elds.
117
The physical motivation lies in the fact that the gauge elds couple to everything
which transforms in a non trivial way under the gauge group. Therefore they couple
also to themselves (remember the homogeneous piece of transformation).
To derive the equations of motion for the gauge elds, let us consider the total
action
[
]
/ m) g
B + SA
d4 x (i
(6.143)
V
where
SA =
1 4
d x T r(F F )
2 V
(6.144)
d4 xT r(F F )
(6.145)
(6.146)
from which
SA = 2
(6.147)
= 2
=
V
d4 xT r [( F + ig[B , F ]) B ]
d4 x ( F + ig[B , F ])A AA
(6.148)
where we used the cyclic property of the trace. By taking into account also the free
term for the Dirac elds and the interaction we nd the equations of motion
T A
F A + ig[B , F ]A = g
/ m) = g B
(i
(6.149)
(6.150)
The reason is that under a global transformation of the symmetry group, the gauge
elds are not invariant, said in dierent words they are charged elds with respect
118
to the gauge elds. In fact we can verify immediately that the previous currents are
precisely the Noethers currents. Under a global variation we have
AC = fCAB A AB ,
= iA T A
(6.151)
and we get
j =
L
L
+
AC
,
A, C
(6.152)
from which
A T A FC f AB A AB
j =
C
(6.153)
In the case of simple compact Lie groups one can dene f ABC = fCAB with the
property f ABC = f BCA . It follows
FC f ABC AB T A = i[B , F ] i[B , F ]A T A
(6.154)
A T A i[B , F ]A A
j =
(6.155)
Therefore
After division by A we get the Noethers currents (6.150). The contribution of the
gauge elds to the currents is also crucial in order they are conserved quantities. In
fact, the divergence of the fermionic contribution is given by
T A ) = ig
T A B + ig
B T A = ig
[T A , B ]
(
(6.156)
which vanishes for abelian gauge elds, whereas it is compensated by the gauge
elds contribution in the non abelian case.
6.6
We have seen that the spontaneous symmetry breaking mechanism, in the case of
continuous symmetry leads to massless scalar particles, the Goldstone bosons. Also
gauge theories lead to massless vector bosons, in fact, as in the electromagnetic case,
gauge invariance forbids the presence in the lagrangian of terms quadratic in the
elds. Unfortunately in nature the only massless particles we know are the photon
and perhaps the neutrinos, which however are fermions. But once one couples
spontaneous symmetry breaking to a gauge symmetry, things change. In fact, if we
look back at the hypotheses underlying a gauge theory, it turns out that Goldstone
theorem does not hold in this context. The reason is that it is impossible to quantize
a gauge theory in a way which is at the same time manifestly covariant and has a
Hilbert space with positive denite metric. This is well known already for the
electromagnetic eld, where one has to choose the gauge before quantization. What
happens is that, if one chooses a physical gauge, as the Coulomb gauge, in order
to have a Hilbert space spanned by only the physical states, than the theory looses
the manifest covariance. If one goes to a covariant gauge, as the Lorentz one, the
119
theory is covariant but one has to work with a big Hilbert space, with non-denite
positive metric, and where the physical states are extracted through a supplementary
condition. The way in which the Goldstone theorem is evaded is that the Goldstone
bosons disappear, and, at the same time, the gauge bosons corresponding to the
broken symmetries acquire mass. This is the famous Higgs mechanism.
Let us start with a scalar theory invariant under O(2)
1 2 2
L =
( )
2
2
4
(6.157)
and let us analyze the spontaneous symmetry breaking mechanism. If 2 < 0 the
symmetry is broken and we can choose the vacuum as the state
= (v, 0),
v=
(6.158)
After the translation 1 = +v, with 0||0 = 0, we get the potential (m2 = 22 )
1
m4
+ m 2 2 +
V =
16 2
m2
(22 + 2 ) (22 + 2 )2
2
4
(6.159)
In this case the group O(2) is completely broken (except for the discrete symmetry
2 2 ). The Goldstone eld is 2 . This has a peculiar way of transforming
under O(2). In fact, the original elds transform as
1 = 2 ,
2 = 1
(6.160)
2 = + v
(6.161)
from which
= 2 ,
We see that the Goldstone eld undergoes a rotation plus a translation, v. This
is the main reason for the Goldstone particle to be massless. In fact one can have
invariance under translations of the eld, only if the potential is at in the corresponding direction. This is what happens when one moves in a way which is tangent
to the surface of the degenerate vacuums (in this case a circle). How do things change
if our theory is gauge invariant? In that case we should have invariance under a
transformation of the Goldstone eld given by
2 (x) = (x)(x) + (x)v
(6.162)
Since (x) is an arbitrary function of the space-time point, it follows that we can
choose it in such a way to make 2 (x) vanish. In other words it must be possible
to eliminate the Goldstone eld from the theory. This is better seen by using polar
coordinates for the elds, that is
21 + 22 ,
2
sin =
2
1 + 22
120
(6.163)
(6.164)
It should be also noticed that the two coordinate systems coincide when we are close
to the vacuum, as when we are doing perturbation theory. In fact, in that case we
can perform the following expansion
22 + 2 + 2v + v 2 v + ,
2
2
v+
v
(6.165)
Again, if we make the theory invariant under a local transformation, we will have
invariance under
(x) (x) + (x)
(6.166)
By choosing (x) = (x) we can eliminate this last eld from the theory. The only
remaining degree of freedom in the scalar sector is (x).
Let us study the gauging of this model. It is convenient to introduce complex
variables
1
1
= (1 + i2 ),
= (1 i2 )
(6.167)
2
2
The O(2) transformations become phase transformations on
ei
(6.168)
(6.169)
We know that it is possible to promote a global symmetry to a local one by introducing the covariant derivative
( + igA )
(6.170)
from which
1
L = ( igA ) ( + igA ) 2 ( )2 F F
4
(6.171)
1
= ei
2
(6.172)
(6.173)
(6.174)
the lagrangian will depend only on the elds and A (we will put again A = A )
1
2
1
L = ( igA )( + igA ) 2 4 F F
2
2
4
4
(6.175)
In this way the Goldstone boson disappears. We have now to translate the eld
0||0 = 0
= + v,
(6.176)
and we see that this generates a bilinear term in A , coming from the covariant
derivative, given by
1 2 2
g v A A
(6.177)
2
Therefore the gauge eld acquires a mass
m2A = g 2 v 2
(6.178)
It is instructive to count the degrees of freedom before and after the gauge transformation. Before we had 4 degrees of freedom, two from the scalar elds and two
from the gauge eld. After the gauge transformation we have only one degree of
freedom from the scalar sector, but three degrees of freedom from the gauge vector,
because now it is a massive vector eld. The result looks a little bit strange, but
the reason why we may read clearly the number of degrees of freedom only after the
gauge transformation is that before the lagrangian contains a mixing term
A
(6.179)
between the Goldstone eld and the gauge vector which makes complicate to read
the mass of the states. The previous gauge transformation realizes the purpose of
making that term vanish. The gauge in which such a thing happens is called the
unitary gauge.
We will consider now the further example of a symmetry O(N ). The lagrangian
invariant under local transformations is
2
1
L = (D )ij j (D )ik k i i (i i )2
2
2
4
(6.180)
where
g
(6.181)
(D )ij = ij + i (T AB )ij WAB
2
A B
B
l ). In the case of broken symmetry (2 < 0), we
m
where (T AB )lm = i(lA m
choose again the vacuum along the direction N , with v dened as in (6.158)
i = viN
122
(6.182)
Recalling that
Tijab j
min
= 0,
TijaN j
min
= ivia ,
a, b = 1, , N 1
(6.183)
(6.184)
Therefore, the elds WaN associated to the broken directions T aN acquire a mass
g 2 v 2 /2, whereas Wab , associated to the unbroken symmetry O(N 1), remain massless.
In general, if G is the global symmetry group of the lagrangian, H the subgroup
of G leaving invariant the vacuum, and GW the group of local (gauge) symmetries,
GW G, one can divide up the broken generators in two categories. In the rst
category fall the broken generators lying in GW ; they have associated massive vector
bosons. In the second category fall the other broken generators; they have associated
massless Goldstone bosons. Finally the gauge elds associated to generators of GW
lying in H remain massless. From the previous derivation this follows noticing that
the generators of H annihilate the minimum of the elds, leaving the corresponding
gauge bosons massless, whereas the non zero action of the broken generators generate
a mass term for the other gauge elds.
The situation is represented in Fig. 6.1.
We can now show how to eliminate the Goldstone bosons. In fact we can dene
new elds a and as
(
)
aN
iT
a
= e
( + v)
(6.185)
i
iN
aN
eiT a
)
iN
iN i(T aN )iN a = iN + ia a
(6.186)
i (va , + v)
(6.187)
from which
showing that the a s are really the Goldstone elds. The unitary gauge is dened
through the transformation
(
aN
i eiT a
j = iN ( + v)
ij
123
(6.188)
G_W
Fig. 6.1 -This figure shows the various groups, G, the global symmetry of the
lagrangian, H G, the symmetry of the vacuum, and GW , the group of local symmetries. The broken generators in GW correspond to massive vector bosons. The
broken generators do not belonging to GW correspond to massless Goldstone bosons.
Th unbroken generators in GW correspond to massless vector bosons.
(
aN
aN
aN
aN
i
W eiT a W eiT a
eiT a eiT a
g
(6.189)
This transformation eliminates the Goldstone degrees of freedom and the resulting
lagrangian depends on the eld , on the massive vector elds WaN and on the
massless eld Wab . Notice again the counting of the degrees of freedom N +2N (N
1)/2 = N 2 in a generic gauge, and 1 + 3(N 1) + 2(N 1)(N 2)/2 = N 2 in the
unitary gauge.
124
Chapter 7
Time ordered products
7.1
One of the most relevant quantities in perturbative eld theory is the propagator,
that is the vacuum expectation value of a time ordered product of two elds. To
introduce the propagator from a physical point of view we will consider a charged
Klein-Gordon. As we know from Section 3.6, the eld destroys a particle of charge
+1 and creates a particle of charge 1. In any case the net variation of the charge
is 1. In analogous way the eld gives rise to a net variation of the charge equal
to +1. Let us now construct a state with charge +1 applying to the vacuum
|(y , t) = (y)|0 =
d3 k
2k
(2)3
eiky |k, m, 1
(7.1)
where |k, m, 1 is the single particle state with charge +1, momentum k and mass m.
We want to evaluate the probability amplitude for the state, |(y , t), to propagate
to the same state, |(x, t ), at a later time t > t. This is given by the matrix
element
(t t)(x, t )|(y , t) = (t t)0|(x, t ) (y , t)|0
(7.2)
It turns out to be convenient to think to this matrix element as the one corresponding
to the creation of a charge +1 at the point y and time t, and to its annihilation at
the point x and time t . This interpretation is a correct one, since the state |(y , t)
is eigenstate of the charge density operator
(z, t)|(y , t) = [(z, t), (y , t)]|0 = + 3 (z y ) (y , t)|0 = 3 (z y )|(y , t)
(7.3)
where we have used eq. (3.197)
=: i t :
()
(7.4)
.
n
t'
t
+1
-1
t'
Fig. 7.1 - The two probability amplitudes contributing to the process np np.
However for t < t we could reach the same result by creating a particle of charge
-1 at (x, t ), and annihilating it at (y , t). The corresponding amplitude is
(t t )0| (y , t)(x, t )|0
(7.5)
The situation is represented in Fig. 7.1, where we have considered the case of a
charged particle exchanged between a proton (charge +1) and a neutron(charge 0).
The total amplitude is obtained by adding the two contributions together. We dene
the vacuum expectation value of the time ordered product (T product) of two elds
as
0|T ((x) (y)|0 = 0|T ( (y)(x))|0
= (x0 y0 )0|(x) (y)|0 + (y0 x0 )0|(y) (x)|0
iF (x y)
(7.6)
The function F (x y) is called the Feynman propagator, and we will show
immediately that it depends indeed on the dierence of the two coordinates x and
y. Using the expressions (3.193) for the elds, we nd
iF (x y) =
d3 k 1 [
(x0 y0 )eik(x y)
=
3
(2) 2k
]
+ (y0 x0 )eik(x y)
(7.7)
This representation can be veried immediately by noticing that for t < 0 the integral is convergent in the upper complex half-plane of . Since there no singularities
in this region (the integral has a pole at = i)), we see that the integral vanishes.
In the case t > 0 the integral is convergent in the lower half-plane. Then we pick
up the contribution of the pole (in a clockwise direction) and we nd
iF (x y) = i
d3 k
1 [ ei(x0 y0 ) ik(x y)
d
e
+ i
(2)4
2k
]
i(x0 y0 )
+ e + i eik(x y)
(7.9)
d4 k 1 eik(x y)
eik(x y)
+
(2)4 2k k0 k + i k0 k + i
[
1
1
d4 k 1 ik(x y)
e
= i
(2)4 2k
k0 k + i k0 + k i
d4 k eik(x y)
= i
(7.10)
(2)4 k 2 m2 + i
where we have dened = 2k . Notice that is a positive quantity. Then
F (x y) =
d4 k eik(x y)
(2)4 k 2 m2 + i
(7.11)
From this representation it follows that F (x) is a Green function for the KleinGordon operator
( 2 + m2 )F (x) = 4 (x)
(7.12)
That the T -product is a Green function for the Klein-Gordon operator is a simple
consequence of its very denition, by using the canonical commutators
( 2 + m2 )x 0|T ((x) (y))|0 = 02 0|T ((x) (y))|0
2 + m2 )(x) (y))|0
+ 0|T ((
x
2
+ m )(x) (y))|0
+ 0|T ((
x
= 0|(x0 y0 )[(x),
(y)]|0 + 0|T (( 2 + m2 )(x) (y))|0
x
= i 4 (x y)
(7.13)
It is easily seen that an analogous result holds for the hermitian Klein-Gordon eld,
that is
F (x y) = i0|T ((x)(y))|0
(7.14)
127
.
Im()
Re()
- k
k
C+
Re()
- k
C_
C (x) =
d4 k eikx
(2)4 k 2 m2
(7.15)
(x) = i
128
d4 k eikx
(2)4 k 2 m2
(7.16)
Im()
Re()
- k
d4 k
eikx
d3 k eikx
(x) = i
=
= [(+) (x), () (y)]
4
3
(2) 2k
C+ (2) (k0 k )(k0 + k )
(7.17)
where we have used eq. (3.149). Also
(+)
()
(x) =
d3 k eikx
(7.18)
(7.19)
It follows
and using eq. (3.147), i(x y) = [(x), (y)]. Therefore the commutator can be
represented as
d4 k eikx
(x) =
(7.20)
C (2)4 k 2 m2
with C given in Fig. 7.4.
In the case of F , the poles position is the one in Fig. 7.5. Then it can be also
dened by taking the poles on the real axis and choosing an integration path CF ,
as specied in Fig. 7.6. That is
F (x) =
CF
d4 k eikx
(2)4 k 2 m2
129
(7.21)
.
Im()
Re()
- k
CF
Re()
k
- k
(7.22)
Let us also notice that all the C invariant functions dened on a close path C,
satisfy the homogeneous Klein-Gordon equation. In fact, the action of the KleinGordon operator removes the singularities from the integrand, leaving the integral
of an analytic function, which vanishes due to the Cauchy theorem.
For a free Dirac eld, the Feynman propagator has a similar denition
SF (x y) = i0|T ( (x) (y))|0
(7.23)
(7.24)
Notice that
T ( (x) (y)) = T ( (y) (x))
130
(7.25)
(an analogous property holds for the Klein-Gordon case, but with a plus sign). The
minus sign introduced in the denition of the T -product for the Dirac eld, is needed
because only in this way it may represent the Green function for the Dirac operator.
In fact
(ix m) T ( (x) (y))
[
]
= i(0 ) (x0 y0 ) (x), (y) = i(x0 y0 )(0 ) (0 ) 3 (x y )
+
= i (x y)
4
(7.26)
that is
(i m)SF (x) = 4 (x)
(7.27)
It is clear that the choice of sign is related to the way in which we perform the
canonical quantization. From the property of SF of being the Green function of the
Dirac operator, we can see that
SF (x) = (i + m)F (x)
(7.28)
and
SF (x) =
k + m
d4 k ikx k + m
d4 k ikx
e
=
e
(2)4
k 2 m2 + i
k 2 m2
CF (2)4
(7.29)
Finally we consider the photon propagator. The only dierence with the KleinGordon case is that the polarization vector give an extra factor g , and therefore
0|T (A (x)A (y)|0 = ig
Dening
D(x) =
d4 k eikx
(2)4 k 2 + i
d4 k eikx
(2)4 k 2 + i
(7.30)
(7.31)
we get
0|T (A (x)A (y)|0 = +ig D(x y)
7.2
(7.32)
The choice of the integration paths in eq. (7.15) allows us to dene various types of
Greens functions according to the boundary conditions we require. Suppose that
we want to solve the Klein-Gordon in a given external source
( 2 + m2 )(x) = j(x)
(7.33)
(7.34)
.
Im()
C ret
Re()
k
- k
C adv
Fig. 7.7 - The integration paths for Gret and Gadv .
In fact,
d4 y G(x y)j(y)
(0)
(x) = (x) +
(7.35)
where (0) (x) satisfy the homogeneous Klein-Gordon equation and it is chosen in
such a way that (x) satises the boundary conditions of the problem. For instance,
if we give the the function (x) at t = we may require
lim (x) = lim (0) (x)
(7.36)
Then, to satisfy the boundary conditions it is enough to choose for G(x) the retarded
solution dened by
Gret (x, x0 < 0) = 0
(7.37)
Such a solution can be found easily by applying the method we have illustrated
in the previous Section. By choosing the integration path as in Fig. 7.7, that is,
leaving both poles below the path we get
Gret (x) =
Cret
d4 k eikx
(2)4 k 2 m2
(7.38)
Clearly Gret (x) vanishes for x0 < 0. In fact, in this case, we can close the path
on the half-plane Im > 0 without hitting any singularity, therefore the function
vanishes. In analogous way we can dene a function Gadv (x) vanishing for x0 > 0, by
choosing a path below the poles (see Fig. 7.7). By integrating explicitly over one
sees easily that the retarded solution propagates forward in time both the positive
and negative energy solutions, whereas the advanced one propagates both solutions
backward in time. By using the expression (7.7) for the Feynman propagator we see
that it propagates the positive energy solutions forward in time and the negative
energy ones backward in time. In fact, Feynman and Stueckelberg showed that the
backward propagation of the negative energy solutions is equivalent to the forward
propagation of the anti-particles. The Feynman propagator acquires its full meaning
only in the quantum theory, where. as we shall prove, it represents the central
element of the perturbation theory.
132
In order to stress the relevance of the Feynman propagator we will consider now
a simple application of what we have learned so far. Let us consider two static
point-like electric charges placed at x1 and x2 . They can be described by the charge
density
j0 (x) =
em 3 (x xm )
(7.39)
m=1,2
where em are the values of the electric charges. The current density vanishes because
we have supposed the charges to be static. Therefore the interaction hamiltonian is
Hint = Lint =
d3 x j (x)A (x) =
em A0 (xm , 0)
(7.40)
m=1,2
In this equation we have taken the electromagnetic eld operator A0 at the time
t = 0, because in this case it will be more convenient to use the Schrodinger representation. We have worked so far with the Heisenberg representation because it
turns out to be more convenient from the point of view of the relativistic covariance
of the formalism, but the problem we will be interested here is the evaluation of the
interaction energy between the two charges. We recall that the Heisenberg and the
Schrodinger representation, as far as the operators are concerned, are related by
AH (t) = eiHt AS eiHt
(7.41)
where AH (t) is the operator in the Heisenberg representation and AS is the operator
in the Schrodinger one. Therefore the two operators are the same at t = 0. To
evaluate the interaction energy we can do a perturbative calculation by evaluating
the energy shift induced by the interaction hamiltonian. Since the electric charges
are classical we are quantizing only the photon eld A0 (x, 0), and the state that we
are perturbing is the vacuum state (the state without photons). Since
0|Hint |0 = 0
(7.42)
we must evaluate the energy shift at the second order in the perturbative theory.
We will have
E =
0|Hint |nn|Hint |0
E0 En
d3 k
0|Hint |kk|Hint |0
k
(7.43)
The only state which contributes in the sum is the state with a single photon of
energy k . By using the expression for Hint and the following representation for
1/k
1
= lim+ i
dt ei(k i)t dt
(7.44)
k 0
0
we can write E in the form (notice that we have suppressed the limit because the
resulting expression is regular at for 0+ ).
E = i
m,n
em en
d3 k
0
em en
d3 k
m,n
= i
(7.45)
m,n
where H0 is the free hamiltonian for the electromagnetic eld. To get this expression
we have used the completeness (recalling again that A0 couples only states that dier
by a photon), and now the operators appearing in the last line can be interpreted as
operators in the Heisenberg representation. By writing explicitly the various terms
in the sum we get
E = i
dt(t) e21 0|A0 (x1 , t)A0 (x1 , 0)|0 + e1 e2 0|A0 (x1 , t)A0 (x2 , 0)|0
]
+ e1 e2 0|A0 (x2 , t)A0 (x1 , 0)|0 + e22 0|A0 (x2 , t)A0 (x2 , 0)|0
(7.46)
The intermediate states description of these four contributions is given in Fig. 7.8.
Notice that the rst and the fourth diagram describe a correction to the intrinsic
properties of the charges, and since we are interested in the evaluation of the interaction energy we can omit them from our calculation. Then, the energy interaction
E12 is given by
E12 = ie1 e2
(7.47)
(7.48)
and using H0 |0 = 0,
0|A0 (x2 , t)A0 (x1 , 0)|0 = 0|eiH0 t A0 (x2 , 0)eiH0 t A0 (x1 , 0)eiH0 t |0
= 0|A0 (x2 , 0)A0 (x1 , t)|0
(7.49)
134
(7.50)
We see that the energy interaction is expressed in terms of the Feynman propagator.
Recalling the eqs.(7.31) and (7.32) we get
d4 k eik(x1 x2 )
(2)4
k 2 + i
d3 k eik(x1 x2 )
e1 e2
1
= e1 e2
=
3
2
(2)
4 |x1 x2 |
k
E12 = e1 e2
dt
(7.51)
and
from which
d3 k eik x
= 3 (x)
3
2
(2) |k|
(7.52)
2 1 = 4 3 (x)
|x|
(7.53)
d3 k eik x
1 1
=
(2)3 |k|2
4 |x|
(7.54)
135
Chapter 8
Perturbation theory
8.1
(8.2)
Lint. = ie ( ) A + e2 A2
(8.3)
j = ie ( )
(8.4)
but another interacting term appears. This term is a straight consequence of the
gauge invariance. In fact the current j which was conserve in absence of the interaction is now neither conserved, neither gauge invariant. Consider the innitesimal
gauge transformation
(x) = ie(x)(x),
then
A (x) = (x)
(8.5)
Lfree = ie, ie , = j
(8.6)
136
we nd
Lint = j A + e2 A2
(8.7)
Lint = j , (j )A + 2e2 , A
(8.8)
The rst term cancels with the variation of Lfree , whereas the other two terms cancel
among themselves
[
(8.9)
This shows that the A2 term is necessary to compensate the fact that j = 0 since
the current is not gauge invariant. In fact the conserved and gauge invariant current
comes by using the Noethers theorem
[
(8.10)
The situation is far more simple in the case of the Dirac equation where
m) 1 F F (i
eA m) 1 F F
Lfree = (i
4
4
(8.11)
(8.12)
Here the gauge eld is coupled to a conserved and gauge invariant current. As a
consequence j A is the only interaction term. In fact,
, = j ,
Lfree = e
(8.13)
Lint = j A = j ,
(8.14)
and
The canonical quantization for an interacting system follows the same procedure
as in the non interacting case. We require canonical commutation and/or anticommutation relations at equal times for the various elds. For dierent elds we
require equal time vanishing commutation (anticommutation) relations for spin integer (half-integer) elds, whereas we require zero commutation relations among elds
of integer spin and elds of half-integer spin. Usually the canonical commutation
relations among the elds are not changed by the interactions with respect to the
free case. However this is not the case if the interaction term involves derivatives of
the elds. This follows from the denition of the canonical momentum densities
i =
Lfree Lint
L
=
+
i
i
i
137
(8.15)
L
= ie A0 ,
= = + ieA0
(8.16)
Since the canonical momenta contain the time component of the gauge eld, one can
verify that the canonical commutators among the scalar elds and their derivatives
are changed by the interaction. Also the propagators are modied. However we
will not insist on this point, because in practice it has no consequences on the
perturbation theory (see later). When derivative interactions are not present the
canonical momentum densities coincide with the free ones, and we get
H = L = Lfree Lint = Hfree Lint
(8.17)
Hint = Lint
(8.18)
and therefore
This is what happens for the interaction between a Dirac and the electromagnetic
elds. The corresponding theory is called QED (Quantum Electro Dynamics). We
recall also that in general the hamiltonian and the electromagnetic current are normal ordered in such a way that the vacuum is an eigenstate of these operators with
vanishing eigenvalues. Therefore the interaction term is written as
: A
Lint = e :
(8.19)
(8.20)
] ]
=: Q :
(8.21)
8.2
The scattering processes are a central element in the study of the elementary particles, since they are the only experimental technique available. In the typical scattering process the incoming particles are prepared in a state of denite momentum,
138
after that the scattering process os some target has taken place, one looks at the nal
states. In ordinary quantum mechanics this situation is well described by using free
wave functions for the initial and nal states. This description is certainly correct
if one has to do with short-range potentials. In eld theory this representation is
not really correct, since also in absence of reciprocal interactions the particles have
self interactions as we have already noticed. For instance, a real electron can be
thought of as if it would have a surrounding cloud of photons which can be emitted and absorbed also when very far from other electrons. A rigorous treatment of
these problems is highly non trivial and it is outside of the scopes of this course.
Therefore we will conne ourselves to a rather intuitive treatment of the problem.
On the other side the imitations of the method will be rather obvious so it may well
constitute the basis for a more rened approach. To simplify the matter we will
make use of the adiabatic hypothesis. This consists in looking at a scattering
process in the following way. At time t = we will suppose that our system
can be described in terms of free particles, that is with the interaction turned o.
Between t = and a time t = T , much before the scattering process takes
place, we let the coupling describing the interaction grow from zero to its actual
value. In the interval T < t < +T , the coupling stays at this value, and then from
t = +T and t = + the coupling goes again to zero (see Fig. 8.1)
.
g(t)
-T
(8.22)
performing all the calculations and taking the limit 0+ at the end. The consistency of this procedure has been shown by various authors and a detailed discussion
139
can be found, for instance, in the book by Jauch and Rohrlich, Theory of Photons
and Electrons.
By using the adiabatic hypothesis we can now discuss the perturbative calculation of the scattering amplitudes. The perturbative expansion will be possible only
if the interaction term is small. For in stance in QED one gets a series of powers
in the ne structure constant e2 /4 1/137, therefore, if the coecients of the
expansion do not grow too much, the expansion is justied. Let us start with the
equation of motion for the states in the Schrodinger representation
i
|S (t)
= HS |S (t)
t
(8.23)
Suppose also that we have two interacting elds A and B. Then we can write
HS = HS0 + HSI
(8.24)
(8.25)
(8.26)
with
and
where HS0 (A) and HS0 (B) are the free hamiltonians for the elds A and B, and
HSI is the interaction hamiltonian. It turns out convenient to introduce a new
representation for the vectors of state, the interaction representation. This is dened
by the following unitary transformation upon the states and on the operators in the
Schrodinger representation
0
0
O(t) = eiHS t OS eiHS t
(8.27)
Of course the matrix elements of any operator in the interaction representation are
the same as in the Schrodinger representation
(t)|O(t)|(t) = S (t)|OS |S (t)
(8.28)
We have also HS0 = H 0 , where H 0 is the free hamiltonian in the interaction representation. Notice also that the interaction representation coincides with the Heisenberg
representation when we switch o the interaction. In the interaction representation
the time evolution of the states is dictated by the interaction hamiltonian
(
|(t)
= i
eiHS t |S (t)
i
t
t
0
0
= HS0 eiHS t |S (t) + eiHS t (HS0 + HSI )|S (t)
0
0
= eiHS t H I eiHS t |(t)
S
from which
i
|(t)
= H I |(t)
t
140
(8.29)
(8.30)
(8.32)
We will dene the S matrix as the operator that give us |(+) once we know
|()
|(+) = S|()
(8.33)
The amplitude Sf i is then
Sf i = f |S|i
(8.34)
ThereforeSf i is the S matrix element between free states. To evaluate the S matrix
we rst transform the Schrodinger equation in the interaction representation in an
integral equation
|(t) = |() i
(8.35)
One can verify that this indeed a solution, and furthermore it satises explicitly the
boundary condition at t = . The perturbative expansion consists in evaluating
|(t) by iterating this integral equation
|(t) = |() i
t1
1i
]
+ |()
dt1
t1
141
(i)
dt1
n=1
t1
dt2
tn1
(8.38)
(i)n
n!
n=1
dt1
(8.39)
The T -product of n terms means that the factors have to be written from left to
right with decreasing times. For instance, if t1 t2 tn , then
(
(8.40)
The equality of the two expressions (8.38) and (8.39) holds term by term. As an
example, consider n = 2. The term in eq. (8.39) can be written as
t2
t2
A =
t1
t2
dt ds T H I (t)H I (s)
t1
(
I
dt H (t)
t1
ds H (s) +
t1
t2
t2
dt
t1
)
I
By looking at Figs. 8.2 and 8.3 one sees easily that exchanging the integrations on
s and t one gets
t2
t2
t2
s
dt
ds =
ds
dt
(8.42)
t1
t1
t1
s
t2
t2
t
t1
t1
t2
142
t2
t1
dt
t2
t
ds.
s
t2
t1
t1
t1
t2
t2
t1
ds
s
t1
dt.
Therefore
t2
A=
dt
t1
ds H I (t)H I (s) +
t1
t2
ds
t1
dt H I (s)H I (t)
(8.43)
t1
t2
A=2
dt
t1
ds H I (t)H I (s)
(8.44)
t1
The result for the nth term in the series can be obtained in a completely analogous
way.
Since the S matrix connects the set of free states at t = with a set of free
states at t = +, it should represent simply a change of basis and as such it should
be unitary. From this point of view the unitarity property of the S matrix is a
very fundamental one because it has to do with the very fundamental properties of
quantum mechanics. So it is important to check that at least formally (this refers
to the fact that we dont really know if the series which we have found for the S
matrix is a convergent one) the expression (8.39) represents a unitary operator. In
order to do that we start rewriting S in the form
( i
S=T e
dt H I (t) )
(8.45)
This expression is a symbolic one and it is really dened by its series expansion
S=
(i)n
n=0
n!
)n
I
dt H (t)
(i)n
n=0
n!
(8.46)
143
The motivation for introducing the T -ordered exponential is that it satises the
following factorization property
T e
t3
O(t)dt )
t3
=T e
t1
t2
O(t)dt ) (
T e
t2
O(t)dt )
(8.47)
t1
)n
t3
O(t)dt
t1
t1
n!
k=0 (n k)!k!
t2
t1
n
t3
t1
t2
t1
t2
t3
t2 )
t2
t1
t3
t1
t2
t2 ) ( t3
t3
t3
t2
t2 )
t3
+
t2
t1
ds1 dsnk
t3
t3
n!
k)!k!
t
t
2
2
k=0
t2
t1
t2
t1
(8.48)
In the last term we have used the fact that all the times zi are smaller than the
times si . What we have proved is the relation
(
t3
)n
dt O(t)
t1
n!
=
T
k=0 (n k)!k!
)nk
t3
dt O(t)
)k
t2
dt O(t)
(8.49)
t1
t2
The factorization property (8.47) follows immediately if we remember that the analogous property for the ordinary exponential
ea + b = ea eb
(8.50)
1
1
n!
a
+
b
n
e
=
(a + b) =
ank bk
n=0 n!
n=0 n! k=0 (n k)!k!
using
n=0 k=0
(8.51)
(8.52)
k=0 n=k
and putting h = n k. Since eq. (8.49) generalizes the binomial formula to T products of powers of time integrals of operators, by the same token we get the
formula (8.47). With this property we can now prove the unitarity of any operator
of the form
tf
( i
dt O(t) )
ti
U =T e
(8.53)
144
with O(t)a hermitian operator. To this end let us divide the time interval (ti , tf ) in
N innitesimal intervals t with
ti t1 t2 tN = tf
(8.54)
(8.55)
(8.56)
from which
S=T e
dt H I (t) )
( +i
=T e
d4 x Lint )
(8.57)
It follows that if the theory is Lorentz invariant, also the S matrix enjoys the same
property. One could think that for theories with derivative interactions the Lorentz
invariance is lost. However it is possible to show that also in these theories the
S matrix is given by the same equation. To see that the S matrix is Lorentz
invariant, notice that the statement would be trivial but for the presence of the T product. However this is invariant under proper Lorentz transformations for timelike separations (remember that a proper Lorentz transformation cannot change the
sign of the time component of a four-vector). In the case of space-like separations
consider, for instance, the second order term in the series for S
d4 x1 d4 x2 (x01 x02 )Lint (x1 )Lint (x2 ) + (x02 x01 )Lint (x2 )Lint (x1 )
(8.58)
Since we integrate over x1 and x2 , the separation between the two points can be either space-like either time-like. However, if the lagrangian density is a local function
of the elds, it follows that
[Lint (x1 ), Lint (x2 )] = 0
per
(x1 x2 )2 < 0
(8.59)
(8.60)
This shows that the T -product of local invariant Lorentz operators is Lorentz invariant. We had an example of this property when we evaluated the propagator for
a scalar eld.
145
8.3
In the previous Section we have shown that the S matrix can be evaluated in terms
of matrix elements of T -products. As we shall see in the applications, the matrix
elements of the S matrix between free particle states can in turn be expressed
as vacuum expectation values (VEVs) of T -products. These V EV s satisfy an
important theorem due to Wick that states that the T -products of an arbitrary
number of free elds (the ones we have to do in the interaction representation) can
be expressed as combinations of T -products among two elds, that is in terms of
Feynman propagators. In order to prove the theorem we will use the technique of
generating functionals. That is we we will start by proving the following identity
( i
d4 x j(x)(x) )
T e
=: e
d4 x j(x)(x)
1 4 4
d x d y j(x)j(y)0|T ((x)(y))|0
:e 2
(8.61)
where (x) is a free real scalar eld and j(x) an ordinary real function. The previous
formula can be easily extended to charged scalar, fermionic and photon elds. The
Wicks theorem is then obtained by expanding both sides of this equation in powers
of j(x) and taking the VEV of both sides. Let us start by expanding the left hand
side of the formula, by using the factorization property in eq. (8.47). Let us also
dene
O(t) =
d3 x j(x)(x)
(8.62)
and notice that for free elds [O(t), O(t )] is just an ordinary number (called a cnumber, to be contrasted with operators which are called q-numbers). Dividing
again the interval (ti , tf ) in N pieces of amplitude t as in the previous Section, we
get
( i
T e
tf
dt O(t) )
ti
(8.63)
and using
1
+ [A, B]
A
B
A
+
B
e e =e
e 2
valid if [A, B] commutes con A e B, we get
( i
T e
= lim
tf
dt O(t) )
ti
N
1 2
t [O(t2 ), O(t1 )]
i=3
146
(8.64)
= lim
i=4
1
1
t2 [O(t3 ), O(t2 ) + O(t1 )] t2 [O(t2 ), O(t1 )]
2
e 2
N
1 2
[O(tj ), O(ti )]
it
O(ti ) t
2
1ijN
i=1
e
= lim e
N
i
=e
from which
( i
T e
i
=e
tf
ti
1 tf
dt O(t)
dt1 dt2 (t1 t2 )[O(t1 ), O(t2 )]
e 2 ti
(8.65)
d4 x j(x)(x) )
1 4 4
d x d y j(x)j(y)(x0 y0 )[(x), (y)]
d4 x j(x)(x)
e 2
(8.66)
The next step is to expand the rst exponential on the right hand side of this
equation in normal products. We have
e
d x j(x)(x)
= e
= e
i
i
()
d x j(x)
(x) i
e
d4 x j(x)(+) (x)
1 4 4
+
d x d y j(x)j(y)[() (x), (+) (y)]
e 2
(8.67)
d x j(x)(x)
= :e
d4 x j(x)(x)
:
1 4 4
+
d x d y j(x)j(y)[((x)() , (+) (y)]
e 2
(8.68)
T e
1
+
e 2
d x j(x)(x) )
4
i
=: e
d4 x j(x)(x)
(8.70)
T (1 2 3 ) = : 1 2 3 : +
: i : 0|T (j k )|0
(8.73)
(8.74)
i=j=k=1
T (1 2 3 4 ) =
: 1 2 3 4 : +
i=j=k=l=1
: i j : 0|T (k l )|0
]
(8.75)
and so on. By taking the VEV of these expression, and recalling that the VEV of
a normal product is zero, we get the Wicks theorem. The T -product of two eld
operators is sometimes called the contraction of the two operators. Therefore to
evaluate the VEV of a T -product of an arbitrary number of free elds, it is enough
to consider all the possible contractions of the elds appearing in the T -product.
For instance, from the last of the previous relations we get
4
0|T (1 2 3 4 )|0 =
(8.76)
i=j=k=l=1
An analogous theorem holds for the photon eld. For the fermions one has to
remember that the T -product is dened in a slightly dierent way. This gives a
minus sign any time we have a permutation of the fermion elds which is odd with
respect to the original ordering. As an illustration the previous formula becomes
0|T (1 2 3 4 )|0 =
(8.77)
i=j=k=l=1
where P = 1 is the sign of the permutation (i, j, k, l) with respect to the fundamental one (1, 2, 3, 4) appearing on the right hand side. More explicitly
0|T (1 2 3 4 )|0 = 0|T (1 2 )|00|T (3 4 )|0
0|T (1 3 )|00|T (2 4 )|0
+ 0|T (1 4 )|00|T (2 3 )|0
148
(8.78)
8.4
(+i)n
n=1
n!
(8.79)
(8.80)
We have now to understand how to use the Wicks theorem in the actual situation.
Consider, for simplicity, two scalar elds. From eq. (8.73) we get (at equal times
the T -product and the usual product coincides)
a (x)b (x) =: a (x)b (x) : +0|T (a (x)b (x))|0
(8.81)
(8.82)
from which
Therefore
T (: a (x)b (x) : 1 (x1 ) n (xn )) = T (a (x)b (x)1 (x1 ) n (xn ))
0|T (a (x)b (x))|0T (1 (x1 ) n (xn ))
(8.83)
Since the second term subtracts the contraction between the two operators taken at
the same point, we can generalize the Wick expansion by saying that when normal
product are contained inside a T -product the Wicks expansion applies with the
further rule that the contractions of operators at the same point, inside the normal
product, are vanishing. With this convention we can write
T (: A(x1 )B(x1 ) : : A(xn )B(xn ) :) = T (A(x1 )B(x1 ) A(xn )B(xn ) )
(8.84)
In the case of QED one can get convinced more easily by recalling that
:= 1 [,
]
:
(8.85)
2
and noticing that inside a T -product the elds can be freely commuted except for
taking into account of their statistics. For instance
:) = 1 T ([,
]) = T (
)
T (:
(8.86)
2
For the following analysis it is useful to remember how the various eld operators
act on the kets
+
+
A+
annihilates
annihilates
annihilates
e
e+
149
creates
creates
creates
e+
e
(8.87)
LI = e : (+ + )(A+
+ A ) ( + ) :
(8.88)
+
: + A
:= A
(8.89)
(8.90)
This process corresponds to a positron emitting a photon. This and the other seven
processes described by the S matrix at the rst order
S (1) = ie
A(x)(x)
d4 x : (x)
:
(8.91)
e
e
Fig. 8.4 -Diagrams for the processes described by the S matrix at the first order.
Fig. 8.5 -Diagrams for the processes described by the S matrix at the first order.
However none of these contributions corresponds to a physically possible process
since the four momentum is not conserved. We will show later that the conservation
of the four momentum is a consequence of the theory. For the moment we will
assume it and we will show that for real particles (that is for particles on the mass
shell p2 = m2 ) these processes cannot happen. For instance consider
e (p) e (p ) + (k)
(8.92)
(8.93)
from which
m2 = m2 2p k
(8.94)
where we have used k 2 = 0 for the photon. In the rest frame of the electron we get
mk0 = 0. Therefore the process is possible only for a photon with vanishing four
momentum. Let us now consider the 20 order contribution
S
(2)
(
)
(ie)2 4
1 )A(x
2 )A(x
1 )(x1 )(x
2 )(x2 )
=
d x1 d4 x2 T (x
2!
(8.95)
We can expand S (2) with the Wicks theorem, and classify the various contributions
according to the number of contractions. If we associate to a contraction of two
elds at the points x1 and x2 a line, we see that the terms originating from S (2) can
be obtained by connecting among them the diagrams depicted in Figs. 8.4 and 8.5
in all the possible ways. In this case the only non vanishing contractions are the
and between A andA . Recalling from Section 7.1
ones between and ,
0|T ( (x) (y))|0 = iSF (x y)
0|T (A (x)A (y))|0 = ig D(x y)
(8.96)
we get
S
(2)
(2)
Si
(8.97)
i=1
where
(2)
S1 =
(2)
S2
(2)
S3
(ie)2 4
1 )A(x
2 )A(x
1 )(x1 )(x
2 )(x2 ) :
d x1 d4 x2 : (x
2!
(8.98)
(ie)2 4
1 )A(x
1 )iSF (x1 x2 )A(x
2 )(x2 ) :
d x1 d4 x2 : (x
2!
(ie)2
2 )A(x
2 )iSF (x2 x1 )A(x
1 )(x1 ) :
+
d4 x1 d4 x2 : (x
2!
1 )A(x
1 )iSF (x1 x2 )A(x
2 )(x2 ) : (8.99)
= (ie)2
d4 x1 d4 x2 : (x
(ie)2 4
1 ) (x1 )ig D(x1 x2 )(x
2 ) (x2 ) : (8.100)
=
d x1 d4 x2 : (x
2!
(ie)2 4
1 ) iSF (x1 x2 )ig D(x1 x2 ) (x2 ) :
d x1 d4 x2 : (x
2!
(ie)2 4
2 ) iSF (x2 x1 )ig D(x2 x1 ) (x1 ) :
+
d x1 d4 x2 : (x
2!
1 ) iSF (x1 x2 )ig D(x1 x2 ) (x2 ) : (8.101)
d4 x1 d4 x2 : (x
= (ie)2
(2)
S4 =
151
(ie)2 4
2 )iSF (x2 x1 )A(x
1 )] : (8.102)
d x1 d4 x2 (1) : T r[iSF (x1 x2 )A(x
2!
(ie)2 4
(2)
S6 =
d x1 d4 x2 (1) : T r[iSF (x1 x2 ) iSF (x2 x1 ) ig D(x2 x1 )] :
2!
(8.103)
(2)
(1)
The term S1 is nothing but the product of two processes of type S and it does not
(2)
give rise to real processes. The term S2 is obtained by contracting a two fermionic
elds, and this means to connect with a fermionic line two of the vertices of Figs.
8.4 and 8.5. Th possible external particles are two , two e , two e+ , or a pair e+ e .
Selecting the external states we can get dierent physical processes. One of these
processes is the Compton scattering + e + e . In this case we must select in
(2)
S2 + (x2 ) to destroy the initial electron and to create the nal electron. As for
the photons are concerned, since A is a real eld, we can destroy the initial photon
both in x2 and x1 and create the nal photon in the other point. Therefore we get
two contributions
(2)
S2 (e e ) = Sa + Sb
(8.104)
(2)
S5 =
with
Sa = (ie)2
and
Sb = (ie)
+
+
d4 x1 d4 x2 (x1 ) iSF (x1 x2 ) A
(x1 )A (x2 ) (x2 )
(8.105)
+
+
d4 x1 d4 x2 (x1 ) iSF (x1 x2 ) A
(x2 )A (x1 ) (x2 )
(8.106)
x2
e
x2
x1
Sa
x1
Sb
d4 x1 d4 x2
+
(x1 ) iSF (x1 x2 ) A+
(x1 )A (x2 ) (x2 )
152
(8.107)
(8.108)
x2
e
x1
x2
e
x1
The term
(8.109)
(8.110)
gives rise to two contributions, and other two come from the nal state. The corresponding diagrams are given in Fig. 8.8.
The terms a) and d) dier only for the exchange x1 x2 and therefore they
are equal after having exchanged the integration variables. The same is true for the
terms b) and c). In this way we get a factor 2 which cancels the 2! in the denom(2)
inator. This the same phenomenon as for the term S2 . That is there diagrams
giving rise to the same contribution to the scattering amplitude. At the order n on
has n! equivalent diagrams which cancel the factor n! coming from the expansion of
the S matrix. This means that it is enough to draw all the inequivalent diagrams.
For the electron scattering we have two such diagrams diering for a minus sign
due to the exchange of the fermionic lines. This reects the fact that eld theory
takes automatically into account the statistics of the particles. In the present case
by giving rise to a properly antisymmetrized amplitude.
For the process e+ e e+ e we get the diagrams the diagrams of Fig. 8.9.
The inequivalent diagrams are those of Fig. 8.10 and correspond to the following
contributions
S3 (e+ e e+ e ) = Sa (e+ e e+ e ) + Sb (e+ e e+ e )
(2)
153
(8.111)
a)
b)
d)
c)
d4 x1 d4 x2
(8.112)
d4 x1 d4 x2
(8.113)
(2)
S4
= (ie)
In analogous way the term S2 contributes to the self-energy of the photon. However,
as we shall show in the following, it cannot change the mass of the photon which
154
a)
b)
e
e
d)
c)
S5 = (ie)2
d4 x1 d4 x2
+
(1)T r[iSF (x1 x2 ) (x2 )iSF (x2 x1 ) (x1 )]A
(x1 )A (x2 ) (8.115)
Notice the minus sign which is due to the fact that in a fermionic loop we have to
invert two fermionic elds inside the T -product. The last term is the one where
all the elds are contracted. There are no external particles and the corresponding
diagram of Fig. 8.13 is called a vacuum diagram. These diagrams can be generally
ignored because they contribute to a simple phase factor for the vacuum state.These
diagrams can be generally ignored because they contribute to a phase factor for the
vacuum state.
155
x1
x2
x1 e
x2
Sa
Sb
e
x2
x1
156
x2
x1
x2
x1
x2
x1
8.5
1
D(x) =
(2)4
where
SF (p) =
d4 p eipx D(p)
1
p m + i
(8.116)
(8.117)
(8.118)
and
1
(8.119)
p2 + i
From the elds expansion in terms of creation and annihilation operators, one can
evaluate the action of the positive frequency part of , e A on the one particle
states. We get
D(p) =
(x)|e (p, r) =
dk
m
eikx u(k, s)b(k, s)b (p, r)|0
3
(2) Ek
157
(8.120)
and using
it follows
(x)|e (p, r) =
]
+
= rs 3 (p k)
(8.121)
m
u(p, r)eipx |0
3
(2) Ek
(8.122)
For later applications it is more convenient to use the normalization in a box than
the continuous one. This amounts to the substitution
1
V
(2)3
(8.123)
m
u(p, r)eipx |0
V Ep
(8.124)
m
v(p, r)eipx |0
V Ep
(8.125)
1 ()
(k)eikx |0
2V Ek
(8.126)
(x)|e (p, r) =
+ (x)|e+ (p, r) =
and for the photon
A+
(x)|(k, )
(8.127)
(1)
|i = ie
= ie
(
d4 x e (p ), (k)| (x)A
(x) (x)|e (p)
dx
m
u(p )eip x
V Ep
1
(k)eikx
2V Ek
m
= ie
3 E E 2E
p p
k
) (
m
u(p)eipx
V Ep
(2)4 4 (p + k p)
u(p )(k)u(p)
(8.128)
f |S
(1)
m
|i = (2) (p + k p)
V Ep
4 4
158
m
V Ep
1
M
2V Ek
(8.129)
where
M = ie
u(p )(k)u(p)
(8.130)
is called the Feynman amplitude for the process. Notice that M is a Lorentz invariant quantity. The term (2)4 4 (p + k p) gives the conservation of the four
momentum in the process, whereas the other factors are associated to the various
external particles (incoming and outgoing). As we said previously this process is
not physically possible since it does not respect the four momentum conservation.
This structure is quite general and for any process we will have
the delta-function
expressing the four momentum conservation, and factors as m/V Ep for each ex
ternal fermion and 1/2V Ek for each boson. Now we have to investigate the rules
for evaluating M. To this end let us consider the Compton scattering
|i = |e (p), (k) |f = |e (p ), (k )
(8.131)
The S matrix element for the Compton scattering is given in eqs. (8.104)-(8.105)
S2 (e e ) = Sa + Sb
(2)
(8.132)
and
m
m
1
1
f |Sa |i = (ie)
d x1 d x2
V Ep V Ep 2V Ek 2V Ek
i 4 iq(x1 x2 )
ip x1
ik x1
d qe
SF (q)
u(p )e
(k )e
(2)4
(k)eikx2 u(p)eipx2
2
d4 q
(2)4
= (ie)
d x 1 d x2
m
V Ep
m
V Ep
1
2V Ek
1
2V Ek
i
u(p)
ei(p + k q)x1 ei(p + k q)x2 (k ) (k)
u(p )
q m + i
m
m
1
1
4 4
= (2) (p + k p k)
V Ep V Ep 2V Ek 2V Ek
i
(k ) (k)
u(p )(ie )
(ie )u(p)
(8.133)
p + k m + i
2
m
f |Sa |i == (2)4 4 (p + k p k)
V Ep
m
V Ep
1
2V Ek
1
Ma
2V Ek
(8.134)
with
Ma = (k ) (k)
u(p )(ie )
i
(ie )u(p),
q m + i
159
q =p+k
(8.135)
q=p+k
i
(ie )u(p) (k) (k ),
q m + i
q = p k
(8.136)
If we write down the Compton amplitude for the positrons we see that we must
associate v(p), v(p) to the initial and nal states respectively. This is seen from eq.
(8.125) which shows that the annihilation operator for the positrons is associated
to v(p). For this reason, when drawing the diagrams in momentum space, is often
convenient to invert the direction of the positron lines, in such a way that the barred
spinors are always written to the left and the unbarred to the right (see Fig. 8.15).
160
q=pk
m
M
V Ei
(8.137)
(8.138)
Mb =
u(p2 )(ie )u(p1 )ig D(p2 p1 )
u(p1 )(ie )u(p2 )
(8.139)
(2)
f |S3 |i
4 4
= (2)
(p1
p2
p1 p2 )
i=1
where M = Ma + Mb , with
and
The relative minus sign comes fro the exchange of the two electrons in the initial
state, that is from the Fermi statistics.
From this example we get the further rule for the Feynman diagrams
- for each internal photon line there is a factor ig D(p) (propagator).
As a last example let us consider the electron self-energy
|i = |e (p) |f = |e (p )
the S matrix element is
(2)
f |S4 |i
= e
d x1 d x2
m
V Ep
m
u(p )eip x1
V Ep
161
(8.140)
e
v(q)
u(p)
final state
initial state
v(q)
e
u(p)
p1
p2p2
p1
p2p1
p2
p 2
p 2
p2
a)
b)
d4 q1 iq1 (x1 x2 )
e
iSF (q1 )
(2)4
d4 q2 iq2 (x1 x2 )
e
iD(q2 )u(p)eipx2
(2)4
m
m d 4 q1 d 4 q2
=
(2)4 4 (p q1 q2 )(2)4 4 (p q1 q2 )
V Ep V Ep (2)4 (2)4
m
m
4 4
= (2) (p p)
M
V Ep V Ep
162
(8.141)
d4 q2
u(p )(ie )iSF (p q2 )(ie )u(p)ig D(q2 )
(2)4
(8.142)
q2
p =p
q1=p-q2
In this way the factor (2)4 4 ( pext ) is automatically produced. We can verify
the previous rule in the self-energy case
= (2)4 4 (p p)
d4 q1 d4 q2
(2)4 4 (q1 + q2 p)(2)4 4 (p q1 q2 )
(2)4 (2)4
d4 q2
(2)4
(8.143)
163
Chapter 9
Applications
9.1
The cross-section
Let us consider a scattering process with a set of initial particles with four momenta
pi = (Ei , pi ) which collide and produce a set of nal particles with four momenta
pf = (Ef , pf ). From the rules of the previous Chapter we know that each external
photon line contributes with a factor (1/2V E)1/2 , whereas each external fermionic
line contributes with (m/V E)1/2 . Furthermore the conservation of the total four
Sf i = f i + (2)4 4
pi
pf
fermioni
m
VE
)1/2 (
bosoni
1
2V E
)1/2
(9.1)
where M is the Feynman amplitude which can be evaluated by drawing the corresponding Feynman diagrams and using the rules of the previous Chapter.
Let us consider the typical case of a two particle collision giving rise to an N
particles nal state. Since we are interested to a situation with the nal state
dierent from the initial one, the probability for the transition will be the modulus
square of the second term in eq. (9.1). In doing this operation we encounter the
square of the Dirac delta which is not a denite quantity. However we should recall
that we are quantizing the theory in a box, and really considering the system in a
nite, although large, time interval that we parameterize as (T /2, T /2). Therefore
we have not really to do with the delta function but rather with (Pi,f = i,f pi,f )
(2) (Pf Pi )
4 4
dx
V
T /2
T /2
dtei(Pf Pi )x
(9.2)
Consider one of the factors appearing in this equation, for instance the time integral.
By performing the integration we get (E = Ef Ei ))
T /2
T /2
2 sin(T E/2)
dteiEt =
E
164
(9.3)
4 sin2 (T E/2)
(E)2
(9.4)
On the right hand side we have a function of E, whose integral holds 2T , and
has a peak at E = 0. Therefore in the T limit we have a delta-convergent
sequence
4 sin2 (T E/2)
lim
= 2T (Ef Ei )
(9.5)
T
(E)2
By doing the same operations also for the space integrals we get
2
(2)4 4 (Pf Pi ) (2)4 L3 T 4 (Pf Pi )
(9.6)
where L is the side of the volume V = L3 . Therefore, the probability per unit time
of the transition is
w = V (2)4 4 (Pf Pi )
1 1
2V Ei f 2V Ef
(2m)|M|2
(9.7)
fermioni
(9.8)
w gives the probability per unit time of a transition toward a state with well dened
quantum numbers, but we are rather interested to the nal states having momenta
between pf and pf + dpf . Since in the volume V the momentum is given by p =
2n/L, the number of nal states is given by
(
L
2
)3
d3 p
(9.9)
The cross-section is dened as the probability per unit time divided by the ux of
the ingoing particles, and has the dimensions of a length to the square
[t1 2 t] = [2 ]
(9.10)
from which
m
u(p)eipx
EV
d3 x
d x(x) =
V
165
m
u (p)u(p) = 1
VE
(9.11)
(9.12)
Then the cross-section for getting the nal states with momenta between pf e pf +dpf
is given by
V
V V d3 p f
d = w
dNF = w
(9.13)
vrel
vrel f (2)3
We obtain
d =
V V d 3 pf
1
1
4 4
V
(2)
(P
P
)
f
i
3
2
vrel f (2)
4V E1 E2 f 2V Ef
= (2)4 4 (Pf Pi )
1
4E1 E2 vrel
(2m)
fermioni
(2m)|M|2
fermioni
d pf
|M|2
(2)3 2Ef
(9.14)
(9.15)
E1 E2
but in the frame where the particle 2 is at rest (laboratory frame) we have p2 =
(m2 , 0) and vrel = v1 , from which
|p1 |
E1 E2 |vrel | = E1 m2
= m2 |p1 | = m2 E12 m21
E1
m22 E12
m21 m22
(9.16)
9.2
LI = e e e + A
(9.17)
Fig. 9.1 describes the Feynman diagram for this process at the second order (in this
diagram the arrows are oriented according to the direction of the momenta). Notice
that in contrast to the process e+ e e+ e the diagram in the crossed channel of
Fig. 9.2 is now missing.
The Feynman amplitude is
ig
v(p1 , r1 )(ie )u(p2 , r2 )
(p1 + p2 )2
1
v(p1 , r1 ) u(p2 , r2 )
= ie2 u(p4 , r4 ) v(p3 , r3 )
(p1 + p2 )2
166
(9.18)
p2
p4
p1
e
p3
1
|M|2
4 ri
(9.19)
where
1
u(p2 , r2 ) v(p1 , r1 )
2
(p1 + p2 )
This expression can be written in the following form
M = ie2 v(p3 , r3 ) u(p4 , r4 )
ie2
M=
Amuons
(r3 , r4 )Aelectrons (r1 , r2 )
2
(p1 + p2 )
with
= u(p4 , r4 ) v(p3 , r3 )
Amuons
167
(9.20)
(9.21)
Aelectrons
= v(p1 , r1 ) u(p2 , r2 )
(9.22)
We get
X=
(
1
e4
Aelectrons A
4 (p1 + p2 )4 r1 ,r2
electrons
(Amuons Amuons
)
(9.23)
r3 ,r4
A
Aelectrons
electrons
v(p1 , r1 ) u(p2 , r2 )
u(p2 , r2 ) v(p1 , r1 )
r1 ,r2
r1 ,r2
(9.24)
and using the eqs. (4.100) and (4.101) for the positive and negative energy projectors, we obtain
[
]
p
m
p
+
m
1
e
2
e
(9.25)
Aelectrons
= Tr
2me
2me
and the analogous quantity for the muons
[
Amuons
p4 + m p3 m
= Tr
2m
2m
(9.26)
To evaluate the trace of Dirac matrices we may use several theorems. Let us
start showing that the trace of an odd number of gamma matrices is zero. In fact
for odd n
T r [
a1 a
n ] = T r [
a1 a
n 5 5 ] = T r [5 a
1 a
n 5 ]
n
= (1) T r [
a1 a
n ]
(9.27)
where we have used the cyclic property of the trace and the anticommutativity of
5 and . Obviously we have
T r[1] = 4
(9.28)
Furthermore
1
1
T r[
ab] = T r[
ab + b
a] = a b T r([ , ]+ ) = 4a b g = 4a b
2
2
(9.29)
Then, using a
b = b
a + 2a b we can evaluate
T r[
a1 a
2 a
3 a
4 ] = T r[(
a2 a
1 + 2a1 a2 )
a3 a
4 ]
= T r[
a2 (
a3 a
1 + 2a1 a3 )
a4 ] + 8(a1 a2 )(a3 a4 )
= T r[
a2 a
3 (
a4 a
1 + 2a1 a4 )] 8(a1 a3 )(a2 a4 ) + 8(a1 a2 )(a3 a4 )
= T r[
a2 a
3 a
4 a
1 ] + 8(a1 a4 )(a2 a3 )
8(a1 a3 )(a2 a4 ) + 8(a1 a2 )(a3 a4 )
(9.30)
that is
T r[
a1 a
2 a
3 a
4 ] = 4[(a1 a2 )(a3 a4 ) (a1 a3 )(a2 a4 ) + (a1 a4 )(a2 a3 )]
168
(9.31)
This relation can be easily extended by induction. Other useful properties are
= 4
(9.32)
a
= (
a + 2a ) = 2
a
(9.33)
a
b = 4a b
(9.34)
and
Let us go back to our process. Evaluating the trace we get
1
T r[(
p1 p2 ) m2e ]
4m2e
1
=
[p1 p2 g (p1 p2 ) + p1 p2 m2e g ]
2
me
1
=
[p p + p1 p2 g (p1 p2 + m2e )]
m2e 1 2
A
electrons =
(9.35)
Analogously we get
A
muons =
1
[p p + p3 p4 g (p3 p4 + m2 )]
m2 3 4
(9.36)
e4
1
1
[2(p1 p3 )(p2 p4 ) + 2(p1 p4 )(p2 p3 ) 2(p1 p2 )(p3 p4 + m2 )
4
2
4 (p1 + p2 ) me m2
e4
[(p1 p3 )(p2 p4 ) + (p1 p4 )(p2 p3 )
2m2e m2 (p1 + p2 )4
(9.37)
This process is studied in the circular colliders where two beams, one of electrons
and the other of positron with same energy are made to collide, and looking for a
nal pair + . Therefore it is convenient to use the frame of the center of mass
for the pair e+ e . We will choose the momentum variables as in Fig. 9.3 with
p1 = (E, p),
p2 = (E, p),
p3 = (E, p ),
p4 = (E, p )
(9.38)
p1 p4 = p2 p3 = E 2 + pp cos
p 3 p4 = E 2 + p ,
2
169
(p1 + p2 )2 = 4E 2
(9.39)
(9.40)
p3
p1
p2
p4
2
2
2
2
2
2
(E
pp
cos
)
+
(E
+
pp
cos
)
+
m
(E
+
p
)
2m2e m2 16E 4
X =
e4
2
[2E 4 + 2E 2 p cos2 + 2E 2 m2 ]
2
2
4
32me m E
1 2
e4
2
[E + p cos2 + m2 ]
=
2
2
2
16me m E
(9.41)
from which
d =
e4
1
1 2
2
2
2
4 4
[E
+
p
cos
+
m
](2)
(p
+
p
p
)
1
2
3
4
16m2e m2 E 2
4(E 2 + p2 )
d3 p3 d3 p4
(2)6 4E 2
e4 1 2
2
= 4 (p1 + p2 p3 p4 )
[E + p cos2 + m2 ]d3 p3 d3 p4
128 2 E 6
(2me )2 (2m )2
(9.42)
where we have used E 2 = p2 , since we are neglecting the electron mass. We can integrate this expression over four variables using the conservation of the four momentum
given by the delta function. We will integrate over p4 e |p3 |. Using d3 p3 = p 2 dp d
we get
e4 1
d
2
2
=
p dp (E1 + E2 E3 E4 )
(E 2 + p cos2 + m2 )
2
6
d
128 E
170
(E3 + E4 )
p
]1
e 4 p 2 2
2
(E + p cos2 + m2 )
2
6
128 E
(9.43)
(9.44)
+
m
=
(E
+
p
(E + p cos2 + m2 )
d
8 E 6 2p
16E 4 E
(9.45)
(9.46)
(9.47)
2 p
2 p
2 2
2
2
2
2
d(E
+
p
cos
+
m
)
=
2 2E 2 + p + 2m2
4
4
16E E
16E E
3
)
2 (
p
1 2
=
E 2 + p + m2
(9.48)
4
4E E
3
In the high energy limit we can easily estimate the total cross-section. Recalling
that
1 GeV 2 = 0.389 mbarn
(9.49)
we get
(mbarn)
5.6 105
2.17 105
0.389
(E(GeV ))2
(E(GeV ))2
(9.50)
9.3
20 (nbarn)
(E(GeV ))2
(9.51)
Coulomb scattering
as the Coulomb eld of a heavy nucleus. The full electromagnetic will be of course,
the sum of the classical part and of the quantized part. The expansion of the S
matrix is still given by
S =1+
(i)n
n=1
n!
(9.52)
with
LI (x) = e : (x)
(x)[A (x) + Aext
(x)] :
(9.53)
.
p p
Ze
Fig. 9.4 - The Feynman diagram for the Coulomb scattering of an electron.
We will consider here the scattering of an electron from nucleus, that will be
thought as innitely heavy. Therefore it will give rise to a static Coulomb potential.
Let us introduce the Fourier transform of this eld
Aext
x) =
(
d3 q iq x ext
e
A (q)
(2)3
(9.54)
x) :
d4 x : (x)
(x)Aext
(
(9.55)
(1)
|i = ie
(9.56)
)1/2 (
)1/2
)1/2 (
)1/2
m
m
= ie
d4 x eip x u(p , s)
Ep V
Ep V
3
d q iq x ext
e
A (q) u(p, s)eipx
(2)3
m
m
= ie
(2)(E E)
Ep V
Ep V
d3 q
q )u(p, s)
(2)3 3 (p + q p )
u(p , s) Aext
(
(2)3
(
)
m
= (2)(E E)
M
EV
with
(9.57)
(9.58)
Notice that in this case we have only the conservation of the energy, whereas the
spatial momentum is not conserve. In fact the external eld violates the translational
invariance of the theory, and, as a consequence, the nucleus absorbs the momentum
p p from the electron. From the [previous expression we see also that when there
are external elds the Feynman rules are modied, and we have to substitute the
wave function of a photon
(
)1/2
1
()
(9.59)
(q)
2Eq V
with the Fourier transform of the external eld
Aext
q)
(
(9.60)
2
1
m
w = f |S (1) |i = 2(E E)
T
EV
)2
|M|2
(9.61)
(9.62)
V wdNf
m
=
d =
v
2
)2
dE (E E)|M|2 d
173
(9.64)
The dierential cross-section is obtained by integrating over the nal energy of the
electron
(
)
(
)
2
d
m 2
me 2
2
=
|M| =
u
(p , s) Aext
(
q
)u(p,
r)
(9.65)
d
2
2
where q = p p. Averaging over the initial polarizations and summing over the
nal ones we obtain
(
d
me
=
d
2
)2
1 ext
p + m p + m
A (q)Aext
(
q
)T
r
2
2m
2m
(9.66)
2
p
p
g
(p
p)
+
p
p
+
g
m
m2
Assuming now that the external eld is of Coulomb type, we have
T r[...] =
Aext
x)
(
Ze
=
,0
4|x|
(9.67)
(9.68)
and
Ze
|q|2
We see that we need only the terms with = = 0 from the trace
Aext
q) =
0 (
T r[...] =
(9.69)
1
[E 2 + m2 + p p ]
m2
(9.70)
1
[E 2 + m2 + p2 cos ]
2
m
(9.71)
(9.72)
Using
|q|2 = |p p|2 = 4p2 sin2
(9.73)
(9.75)
which is the classical Rutherford formula for the Coulomb scattering with T =
mv 2 /2.
174
Chapter 10
One-loop renormalization
10.1
p p
p p
p
p
pk
p pk
p
p p
p p
p
a)
p
b)
p
pk
k
c)
k+q
p
k
p k
p p
p p
d)
Fig. 10.1 - The Feynman diagram for the Coulomb scattering at the third order in
the electric charge and at the first order in the external field.
175
The Coulomb scattering can be used, in principle, to dene the physical electric
charge of the electron. This is done assuming that the amplitude is linear in ephys ,
from which we get an expansion of the type
ephys = e + a2 e3 + = e(1 + a2 e2 + )
(10.1)
in terms of the parameter e which appears in the original lagrangian. The rst
problem we encounter is that we would like to have the results of our calculation in
terms of measured quantities as ephys . This could be done by inverting the previous
expansion, but, and here comes the second problem, the coecient of the expansion are divergent quantities. To show this, consider, for instance the self-energy
contribution to one of the external photons (as the one in Fig. 10.1a). We have
Ma = u(p )(ie )Aext
p p)
(
[
]
i
ie2 (p) u(p)
p m + i
(10.2)
where
d4 k
ig
i
(ie
)
(ie )
4
2
(2)
k + i p k m + i
1
d4 k
1
= e2
2
4
(2)
k + i p k m + i
ie (p) =
or
(p) = i
(10.3)
d4 k
p k + m
1
(2)4
(p k)2 m2 + i k 2 + i
(10.4)
For large momentum, k, the integrand behaves as 1/k 3 and the integral diverges
linearly. Analogously one can check that all the other third order contributions
diverge. Let us write explicitly the amplitudes for the other diagrams
Mb = u(p )ie2 (p )
Mc = u(p )(ie )
ig 2
ie (q)Aext
p p)u(p),
(
q 2 + i
where
ie (q) = (1)
i
(ie )Aext
p p)u(p)
(
p m + i
(10.5)
q = p p
(10.6)
]
i
i
d4 k
Tr
(ie )
(ie )
4
(2)
k + q m + i
k m + i
(10.7)
(the minus sign originates from the fermion loop) and therefore
(q) = i
1
d4 k
1
Tr
4
(2)
k + q m k m + i
(10.8)
(10.9)
where
e2 (p , p) =
d4 k
i
i
ig
(ie
)
(ie
) 2
(2)4
k + i
p k m + i p k m + i
(10.10)
or
(p , p) = i
d4 k
1
1
1
4
(2)
p k m + i p k m + i k + i
(10.11)
The problem of the divergences is a serious one and in order to give some sense to
the theory we have to dene a way to dene our integrals. This is what is called the
regularization procedure of the Feynman integrals. That is we give a prescription in
order to make the integrals nite. This can be done in various ways, as introducing
an ultraviolet cut-o, or, as we shall see later by the more convenient means of
dimensional regularization. However, we want that the theory does not depend on
the way in which we dene the integrals, otherwise we would have to look for some
physical meaning of the regularization procedure we choose. This bring us to the
other problem, the inversion of eq. (10.1). Since now the coecients are nite we
can indeed perform the inversion and obtaining e as a function of ephys and obtain all
the observables in terms of the physical electric charge (that is the one measured in
the Coulomb scattering). By doing so, a priori we will introduce in the observables
a dependence on the renormalization procedure. We will say that the theory is
renormalizable when this dependence cancels out. Thinking to the regularization in
terms of a cut-o this means that considering the observable quantities in terms of
ephys , and removing the cut-o (that is by taking the limit for the cut-o going to
the innity), the result should be nite. Of course, this cancellation is not obvious
at all, and in fact in most of the theories this does not happen. However there is a
restrict class of renormalizable theories, as for instance QED. We will not discuss the
renormalization at all order and neither we will prove which criteria a theory should
satisfy in order to be renormalizable. We will give these criteria without a proof but
we will try only to justify them in a physical basis. As far QED is concerned we will
study in detail the renormalization at one-loop.
The previous way of dening a renormalizable theory amounts to say that the
original parameters in the lagrangian, as e, should be innite and that their divergences should compensate the divergences of the Feynman diagrams. Then one can
try to separate the innite from the nite part of the parameters (this separation
is ambiguous, see later). The innite contributions are called counterterms, and
by denition they have the same operator structure of the original terms in the
lagrangian. On the other hand, the procedure of regularization can be performed
by adding to the original lagrangian counterterms cooked in such a way that their
contribution kills the divergent part of the Feynman integrals. This means that
the coecients of these counterterms have to be innite. However they can also be
regularized in the same way as the other integrals. We see that the theory will be
177
renormalizable if the counterterms we add to make the theory nite have the same
structure of the original terms in the lagrangian, in fact, if this is the case, they can
be absorbed in the original parameters, which however are arbitrary, because they
have to be xed by the experiments (renormalization conditions).
In the case of QED all the divergences can be brought back to the three functions
(p), (q) e (p , p). This does not mean that an arbitrary diagram is not
divergence, but it can be made nite if the previous functions are such. In such
a case one has only to show that eliminating these three divergences (primitive
divergences) the theory is automatically nite. In particular we will show that
the divergent part of (p) can be absorbed into the denition of the mass of the
electron and a redenition of the electron eld (wave function renormalization).
The divergence in , the photon self-energy, can be absorbed in the wave function
renormalization of the photon (the mass of the photon is not renormalized due to
the gauge invariance). And nally the divergence of (p , p) goes into the denition
of the parameter e. To realize this program we divide up the lagrangian density in
two parts, one written in terms of the physical parameters, the other will contain the
counterterms. We will call also the original parameters and elds of the theory the
bare parameters and the bare elds and we will use an index B in order to distinguish
them from the physical quantities. Therefore the two pieces of the lagrangian should
look like as follows: the piece in terms of the physical parameters
Lp
m) e
A 1 F F 1 ( A )2
Lp = (i
(10.12)
4
2
and the counter terms piece Lc.t.
A
C F F E ( A )2 eDA
Lc.t. = iB
(10.13)
4
2
and we have to require that la sum of these two contributions should coincide with
the original lagrangian written in terms of the bare quantities. Adding together Lp
and Lc.t. we get
(m + A)
e(1 + D)
A 1 + C F F + gauge xing
L = (1 + B)i
4
(10.14)
where, for sake of simplicity, we have omitted the gauge xing term. Dening the
renormalization constant of the elds
Z1 = (1 + D),
Z2 = (1 + B),
Z3 = (1 + C)
(10.15)
1/2
AB = Z3 A
B = Z2 ,
(10.16)
we obtain
B m + A B B eZ1 B B A 1 FB, F + gauge xing
L = iB
B
B
1/2
Z2
4
Z 2 Z3
(10.17)
178
and putting
mB =
m+A
,
Z2
eB =
eZ1
1/2
Z2 Z3
(10.18)
we get
B mB B B eB B B A 1 FB, F + gauge xing (10.19)
L = iB
B
B
4
So we have succeeded in the wanted separation.Notice that the division of the parameters in physical and counter term part is well dened, because the nite piece
is xed to be an observable quantity. This requirement gives the renormalization
conditions. The counter terms A, B, ... are determined recursively at each perturbative order in such a way to eliminate the divergent parts and to respect the
renormalization conditions. We will see later how this works in practice at one-loop
level. Another observation is that Z1 and Z2 have to do with the self-energy of the
electron, and as such they depend on the electron mass. Therefore if we consider
the theory for a dierent particle, as the muon, which has the same interactions as
the electron and diers only for the value of the mass (m 200me ), one would
get a dierent bare electric charge for the two particles. Or, phrased in a dierent
way, one would have to tune the bare electric charge at dierent values in order to
get the same physical charge. This looks very unnatural, but the gauge invariance
of the theory implies that at all the perturbative orders Z1 = Z2 . As a consequence
1/2
eB = e/Z3 , and since Z3 comes from the photon self-energy, the relation between
the bare and the physical electric charge is universal (that is it does not depend on
the kind of charged particle under consideration).
Summarizing, one starts dividing the bare lagrangian in two pieces. Then we
regularize the theory giving some prescription to get nite Feynman integrals. The
part containing the counter terms is determined, order by order, by requiring that
the divergences of the Feynman integrals, which come about when removing the
regularization, are cancelled out by the counter term contributions. Since the separation of an innite quantity into an innite plus a nite term is not well dened, we
use the renormalization conditions, to x the nite part. After evaluating a physical
quantity we remove the regularization. Notice that although the counter terms are
divergent quantity when we remove the cut o, we will order them according to the
power of the coupling in which we are doing the perturbative calculation. That is we
have a double limit, one in the coupling and the other in some parameter (regulator)
which denes the regularization. The order of the limit is rst to work at some order
in the coupling, at xed regulator, and then remove the regularization.
Before going into the calculations for QED we want to illustrate some general
results about the renormalization. If one considers only theory involving scalar,
fermion and massless spin 1 (as the photon) elds, it is not dicult to construct
an algorithm which allows to evaluate the ultraviolet (that is for large momenta)
divergence of any Feynman diagram. In the case of the electron self-energy (see
Figs. 10.1a and 10.1b) one has an integration over the four momentum p and a
179
behaviour of the integrand, coming from the propagators, as 1/p3 , giving a linear
divergence (it turns out that the divergence is only logarithmic). From this counting
one can see that only the lagrangian densities containing monomials in the elds
with mass dimension smaller or equal to the number of space-time dimensions have
a nite number of divergent diagrams. It turns out also that these are renormalizable
theories ( a part some small technicalities). The mass dimensions of the elds can be
easily evaluated from the observation that the action is dimensionless in our units
/ = 1). Therefore, in n space-time dimensions, the lagrangian density, dened as
(h
dn x L
(10.20)
has a mass dimension n. Looking at the kinetic terms of the bosonic elds (two
derivatives) and of the fermionic elds (one derivative), we see that
dim[] = dim[A ] =
n
1,
2
dim[] =
n1
2
(10.21)
In particular, in 4 dimensions the bosonic elds have dimension 1 and the fermionic
3/2. Then, we see that QED is renormalizable, since all the terms in the lagrangian
density have dimensions smaller or equal to 4
= 3,
dim[]
A ] = 4,
dim[
dim[( A )2 ] = 4
(10.22)
The condition on the dimensions of the operators appearing in the lagrangian can
be translated into a condition over the coupling constants. In fact each monomial
Oi will appear multiplied by a coupling gi
L=
gi Oi
(10.23)
therefore
dim[gi ] = 4 dim[Oi ]
(10.24)
dim[Oi ] 4
(10.25)
dim[gi ] 0
(10.26)
(10.27)
Here has dimension 1 and is dimensionless. We see that the linear -models are
renormalizable theories.
180
Giving these facts let us try to understand what makes renormalizable and non
renormalizable theories dierent. In the renormalizable case, if we have written
the most general lagrangian, the only divergent diagrams which appear are the
ones corresponding to the processes described by the operators appearing in the
lagrangian. Therefore adding to L the counter term
Lc.t. =
gi Oi
(10.28)
we can choose the gi in such a way to cancel, order by order, the divergences. The
theory depends on a nite number of arbitrary parameters equal to the number of
parameters gi . Therefore the theory is a predictive one. In the non renormalizable
case, the number of divergent diagrams increase with the perturbative order. At
each order we have to introduce new counter terms having an operator structure
dierent from the original one. At the end the theory will depend on an innite
number of arbitrary parameters. As an example consider a fermionic theory with an
2 . Since this term has dimension 6, the relative coupling
interaction of the type ()
has dimension -2
2
Lint = g2 ()
(10.29)
.
a)
b)
c)
2.
Fig. 10.2 - Divergent diagrams coming from the interaction ()
At one loop the theory gives rise to the divergent diagrams of Fig. 10.2. The
divergence of the rst diagram can be absorbed into a counter term of the original
type
2
g2 ()
(10.30)
The other two need counter terms of the type
3 + g4 ()
4
g3 ()
181
(10.31)
a)
b)
2 and ()
4.
Fig. 10.3 - Divergent diagrams coming from the interactions ()
These counter terms originate new one-loop divergent diagrams, as for instance the
ones in Fig. 10.3. The rst diagram modies the already introduced counter term
4 , but the second one needs a new counter term
()
5
g5 ()
(10.32)
c2
E2
(10.34)
(10.35)
E
182
c
g2
(10.36)
It follows that we can give a meaning also to non renormalizable theories, but only
for a limited range of values of the energy. This range is xed by the value of the
non renormalizable coupling. It is not diculty to realize that non renormalizability
and bad behaviour of the amplitudes at high energies are strictly connected.
10.2
d4 p F (p, k)
I4 (k) =
(10.37)
d2 p F (p, k)
I(, k) =
(10.38)
(z) =
dt et tz1
(10.39)
(10.40)
(z) =
dt et tz1 +
dt et tz1
(10.41)
Expanding the exponential in the rst integral and integrating term by term we get
(1)n
dt et tz1
n!
0
n=0
(1)n n+z
=
+
dt et tz1
n! n + z
n=0
(z) =
n+z1
dt t
(10.42)
The second integral converges for any z since > 0. This expression coincides with
the representation for the function for Re z > 0, but it is dened also for Re z < 0
where it has simple poles located at z = n. Therefore it is a meaningful expression
on all the complex plane z. Notice that in order to isolate the divergences we need
to introduce an arbitrary parameter . However the result does not depend on the
particular value of this parameter. This the Weierstrass representation of the Euler
(z). From this example we see that we need the following three steps
Find a domain where I(, k) is convergent. Typically this will be for Re < 2.
Construct an analytic function identical to I(, k) in the domain of convergence, but dened on a larger domain including the point = 2.
At the end of the calculation take the limit 2.
10.3
dN pF (p2 )
IN =
(10.43)
(10.44)
IN = SN
pN 1 F (p2 )dp
(10.45)
The value of the sphere surface can be evaluated by the following trick. Consider
I=
ex dx =
2
184
(10.46)
IN =
(10.47)
N/2
= SN
2
N 1 e d
(10.48)
By putting x = 2 we have
(
N/2
1
1
N
xN/21 ex dx = SN
= SN
2
2
2
0
(10.49)
where we have used the representation of the Euler function given in the previous
Section. Therefore
2 N/2
(10.50)
SN = ( N )
2
and
N/2
IN =
( )
N
2
xN/21 F (x)dx
(10.51)
with x = p2 .
The integrals we will be interested in are of the type
(M )
IN
dN p
(p2 a2 + i)A
(10.52)
with p a vector in a N dimensional Minkowski space. We can perform an anticlockwise rotation of 900 (Wicks rotation) in the complex plane of p0 without hitting
any singularity . Then we do a change of variables p0 ip0 obtaining
(M )
IN
dN p
= i(1)A IN
(p2 a2 )A
=i
(10.53)
N/2
( )
N
2
xN/21
dx
(x + a2 )A
(10.55)
By puttingx = a2 y we get
2 N/2A
IN = (a )
N/2
( )
N
2
185
y N/21 (1 + y)A dx
(10.56)
and recalling the integral representation for the Euler B(x, y) function (valid for
Rex, y > 0)
(x)(y) x1
B(x, y) =
=
t (1 + t)(x+y) dt
(10.57)
(x + y)
0
it follows
IN = N/2
(A N/2)
1
2
AN/2
(A)
(a )
(10.58)
We have obtained this representation for N/2 > 0 and Re(A N/2) > 0. But we
know how to extend the Euler gamma-function to the entire complex plane, and
therefore we can extend this formula to complex dimensions N = 2
I2 =
(A )
1
2
(A) (a )A
(10.59)
(10.60)
d2 p
(A )
1
= i (1)A
2
2
A
2
(p a )
(A) (a )A
(10.61)
For the following it will be useful to derive another formula. Let us put in the
previous equation p = p + k and b2 = a2 + k 2 , then
1
(A )
d2 p
= i (1)A
2
2
2
A
2
(p + 2p k + k a )
(A) (a )A
(10.62)
d2 p
1
A (A )
=
i
(1)
(p2 + 2p k + b2 )A
(A) (k 2 b2 )A
(10.63)
from which
d2 p
p
k
A (A )
=
i
(1)
(p2 + 2p k + b2 )A
(A) (k 2 b2 )A
(10.64)
and
p p
(1)A
=
i
(p2 + 2p k + b2 )A
(A)(k 2 b2 )A
]
[
1
2
2
(A )k k g (k b )(A 1)
2
d2 p
186
(10.65)
Since at the end of our calculation we will have to take the limit 2, it will be
useful to recall the expansion of the Gamma function around its poles
() =
1
+ O()
(10.66)
where
= 0.5772...
(10.67)
(n + ) =
where
(n + 1) = 1 +
10.4
(1)n 1
+ (n + 1) + O()
n!
1
1
+ +
2
n
(10.68)
(10.69)
(10.71)
= (2 d)
(10.72)
and
Other relations can be obtained by starting from the algebraic properties of the
-matrices. Let us start with the electron self-energy which we will require to have
dimension 1 as in d = 4. From eq. (10.4) we have
(p) = i
42
p k + m
1
d2 k
2
2
2
(2)
(p k) m + i k + i
(10.73)
187
(10.74)
1
1
1 1
1
=
=
ab
ba a b
ba
dx
x2
(10.75)
(10.76)
We get
(p) = i42
dz
0
p k + m
d2 k
(10.77)
(10.78)
and the term p k can be eliminated through the following change of variables
k = k + pz. We nd
[...] = (p2 m2 )z + (k + pz)2 2p (k + pz)z = k m2 z + p2 z(1 z) (10.79)
2
(p) = i42
dz
0
d2 k
p(1 z) k + m
2
2
(2)
[k m2 z + p2 z(1 z)]2
(10.80)
42
(p) = i
dz
0
d2 k
p(1 z) + m
2
2
2
2
2
(2)
[k m z + p z(1 z)]
(10.81)
(p) = i
42
dz
0
1
(2 )
p(1 z) + m
i
2
2
(2)
(2)
[m z p2 z(1 z)]2
(10.82)
By dening = 4 2 we get
(p) =
dz
0
1
p(1 z) + m
(4)/2 (/2) 2
4
(2)
[m z p2 z(1 z)]/2
(10.83)
1
( 2)
p(1 z) + (4 )m
1
(/2)
dz(42 )/2
2
2
16
[m z p2 z(1 z)]/2
0
188
(10.84)
we obtain
(p) =
1
0
1
(/2)
16 2
m2 z p2 z(1 z)
dz[2
p(1 z) 4m (
p(1 z) m)]
42
]/2
(10.85)
Dening
A = 2
p(1 z) 4m,
B =
p(1 z) + m,
m2 z p2 z(1 z)
C=
42
(10.86)
(p) =
=
=
][
1
1
2A
dz
+ 2B A 1 log C
2
16 0
2
[
]
1 1
2A
dz
A
log
C
+
2B
A
16 2 0
1
1
(
p
4m)
[
p 2m + (
p 4m)]
8 2
16 2
m2 z p2 z(1 z)
1 1
dz[
p
(1
z)
2m]
log
8 2 0
42
1
(
p 4m) + nite terms
8 2
(10.87)
1
d2 k
1
(q) = i42
Tr
2
(2)
k + q m k m
d2 k T r[ (k + m) (k + q + m)]
42
= i
(2)2 (k 2 m2 )((k + q)2 m2 )
(10.88)
(q) = i42
dz
0
d2 k
T r[ (k + m) (k + q + m)]
(2)2 [(k 2 m2 )(1 z) + ((k + q)2 m2 )z]2
(10.89)
(10.90)
(10.91)
42
dz
(q) = i
d2 k T r[ (k qz + m) (k + q(1 z) + m)]
(10.92)
(2)2
[k 2 + q 2 z(1 z) m2 ]2
189
Since the integral of the odd terms in k is zero, it is enough to evaluate the contribution of the even term to the trace
T r[...]even = T r[ (k qz + m) (k + q(1 z)]pari + m2 T r[ ]
k]
T r[ q q]z(1 z) + m2 T r[ ]
= T r[ k
(10.93)
If we dene the as matrices of dimension 2 2 we can repeat the calculation of
Section 9.2 obtaining a factor 2 instead of 4. Therefore
T r[...]even = 2 [2k k g k 2 (2q q g q 2 )z(1 z) + m2 g ]
= 2 [2k k 2z(1 z)(q q g q 2 ) g (k 2 m2 + q 2 z(1 z)](10.94)
and we nd
d2 k [
2k k
2
2
2
(2) [k + q z(1 z) m2 ]2
0
]
2z(1 z)(q q g q 2 )
g
[k 2 + q 2 z(1 z) m2 ]2 [k 2 + q 2 z(1 z) m2 ]
(q) = i42 2
dz
(10.95)
d2 p
whereas
2p p
(1 )
= ig 2 1
2
2
a ]
(a )
(10.96)
(1 )
1
= i 2 1
2
a ]
(a )
(10.97)
[p2
d2 p
[p2
Therefore the rst and the third contribution to the vacuum polarization cancel out
and we are left with
d2 k
1
2
2
2
(2) [k + q z(1 z) m2 ]2
0
(10.98)
Notice that the original integral was quadratically divergent, but due to the previous
cancellation the divergence is only logarithmic. The reason is again gauge invariance. In fact it is possible to show that this implies q (q) = 0. Performing the
momentum integration we have
(q) = i42 2 (q q g q 2 )
dz 2z(1 z)
(2 )
i
2
2
(2) [m q 2 z(1 z)]2
0
(10.99)
By putting again = 4 2 and expanding the previous expression
(q) = i2
2 (q q g q )
42
(q) = 2 22/2 (q q g q 2 )
dz z(1 z)
dz z(1 z)
190
2/2
(/2)
4
2
(2) [m q 2 z(1 z)]/2
]/2
1
2 2/2
m2 q 2 z(1 z)
2
=
2
(q
q
g
q
)
dz
z(1
z)(/2)
16 2
42
0
2
=
(4 2 log 2)(q q g q 2 )
2
16
[
]
1
2
2
0
][
8
1
2
(q
q
g
q
)
dz
z(1
z)
4 4 log 2
8 2
1
=
(q q g q 2 )
2
[2
[
]]
m2 q 2 z(1 z)
1
dz z(1 z) log
3 6
22
(q) =
log C
2
(10.101)
1
1
(q q g q 2 ) + nite terms
2
6
(10.102)
(p , p) = i
d2 k p k + m p k + m
1
2
2
2
2
2
2
(2)
(p k) m
(p k) m
k
42
(10.103)
1
i=1
ai
= (n 1)!
(1 n i )
di n i=1 n
[ i=1 i ai ]
i=1
n
1
0
(10.104)
1
i=1
ai
1=
0
di e
i ai
(10.105)
i=1
i=1
d(
i )
(10.106)
i=1
1
i=1
ai
=
0
i=1
di d(
i )e
i=1
i ai
=
i=1
i=1
di n
d
e
i ai
i=1
(10.107)
191
The integration over i can be restricted to the interval [0, 1] due to the delta
function, and furthermore
(n 1)!
dn1 e =
n
(10.108)
d2 k 1 1x
dx
dy
(2)2 0
0
(
p k + m) (
p k + m
(p , p) = 2i42
(10.109)
The denominator is
[...] = k 2 m2 (x + y) + p2 x + p y 2k (px + p y)
2
(10.110)
Changing variable, k = k + px + py
[...] = (k + px + p y)2 m2 (x + y) + p2 x + p y 2(k + px + p y) (px + p y)
2
= k 2 m2 (x + y) + p2 x(1 x) + p y(1 y) 2p p xy
(10.111)
2
Letting again k k
d2 k
(2)2
0
0
(
p (1 y) px k + m) (
p(1 x) p y k + m)
2
(10.112)
[k m2 (x + y) + p2 x(1 x) + p 2 y(1 y) 2p p xy]3
(p , p) = 2i42
1x
dx
dy
The odd term in k is zero after integration, the term in k 2 is logarithmically divergent, whereas the remaining part is convergent. Separating the divergent piece,
(3)
(1)
, from the convergent one, ,
(2)
= (1)
+
(10.113)
(1)
(p , p)
= 2i
42
dx
0
1x
dy
0
d2 k
(2)2
kk
[k 2 m2 (x + y) + p2 x(1 x) + p 2 y(1 y) 2p p xy]3
(
)
1
1x
i (1)3
1
42
= 2i
dx
dy
(2 )
(3)
2
0
0
[m2 (x + y) p2 x(1 x) p 2 y(1 y) + 2p p xy]2
(
)
1
1x
1 42 1
=
(2 )
dx
dy
2
4
0
0
(10.114)
[m2 (x + y) p2 x(1 x) p 2 y(1 y) + 2p p xy]2
192
(10.115)
we obtain ( = 4 2)
(
1 1
(1)
(p , p) =
2
4
)2/2
(/2)( 2)2
1x
dx
0
dy
0
1
[m2 (x + y) p2 x(1 x) p 2 y(1 y) + 2p p xy]/2
[
]
1x
1
1
2
=
dy
dx
[4
4]
32 2
0
0
]/2
m2 (x + y) p2 x(1 x) p 2 y(1 y) + 2p p xy
42
1
1x
1
1
1
dx
dy
= 2
( + 2) 2
8 [
16 2
8
0
0
]
m2 (x + y) p2 x(1 x) p 2 y(1 y) + 2p p xy
log
42
(10.116)
and nally
1
+ nite terms
8 2
In the convergent part we can put directly = 2
(1)
(p , p) =
(2)
(p , p)
10.5
(10.117)
i 1 1x i 2 (1)3
= 4
dx
dy
8 0
(3)
0
(
p (1 y) px + m) (
p(1 x) p y + m)
m2 (x + y) p2 x(1 x) p 2 y(1 y) + 2p p xy
1 1 1x
dx
dy
=
16 2 0
0
(
p (1 y) px + m) (
p(1 x) p y + m)
(10.118)
m2 (x + y) p2 x(1 x) p 2 y(1 y) + 2p p xy
1
(
p 4m) + f (p).
2
8
]
1
+ f (q) (q q g q 2 )(q 2 ).
6 2
(p , p) =
(10.119)
1
+ f (p , p),
2
8
193
(10.120)
(10.121)
where the functions with the superscript f represent the nite contributions. Let us
start discussing the electron self-energy. As shown in eq. (10.2), the eect of (p)
is to correct the electron propagator. In fact we have (see Fig. 10.4):
SF (p) =
i
i
i
+
ie2 (p)
+ ,
p m p m
p m
(10.122)
e2 (p)
i
1+
SF (p) =
p m
p m
)1
i
.
p m + e2 (p)
(10.123)
Therefore the eect of the divergent terms is to modify the coecients of p and m:
(
iSF1 (p)
e2
e2
= p m + e (p) = p 1 + 2 m 1 + 2 + nite terms. (10.124)
8
2
2
(L1 )ct = iB
(10.125)
(L1 )ct modies the Feynman rules adding two operators parameterized by A and B.
These coecients can be evaluated noticing that the expression of the propagator,
taking into account (L1 )ct , is
(
i
i
B p A
1
(1 + B)
p (m + A)
p m
p m
i
i
i
+
(iB p iA)
,
p m p m
p m
(10.126)
where, consistently with our expansion we have taken only the rst order terms in
A and B. We can associate to these two terms the diagrams of Fig. 10.5, with
contributions iA to the mass term, and iB p to p.
The propagator at the second order in the coupling constant is then obtained by
adding the diagrams of Fig. 10.6 to the free part. We get
SF (p) =
)
i ( 2
i
i
+
ie (p) + iB p iA
p m p m
p m
194
(10.127)
.
-iA
iBp/
e (p) + B p A =
2
e2
e2
+
B
p
m + A + nite terms.
8 2
2 2
(10.128)
1
m
+ F2
))
m
1
+ Fm
(10.129)
))
(10.130)
with F2 and Fm nite for 0. Notice that these two functions are dimensionless
and for the moment being completely arbitrary. However they can be determined by
the renormalization conditions, that is by xing the arbitrary constants appearing
in the Lagrangian. In fact, given
iSF1 (p) (2) (p) = p m + B p A + e2 (p)
)
(
me2
e2
=
1 2 F2 p m +
Fm + e2 f ,
8
2 2
195
(10.131)
we can require that at the physical pole, p = m, the propagator coincides with the
free propagator
i
SF (p)
, for p = m.
(10.132)
p m
From here we get two conditions. The rst one is
(2) (
p = m) = 0
from which
e2 f (
p = m)
The second condition is
me2
me2
F
+
Fm = 0.
2
8 2
2 2
(2) (p)
= ,
p p=m
giving
e2
f (p)
e2
F2 = 0.
p p=m 8 2
(10.133)
(10.134)
(10.135)
(10.136)
D
(q) =
ig ig 2 ig
+
ie (q) 2 + ,
q2
q2
q
(10.137)
from which
] ig
ig ig [ 2
2
2
+
(ie
)(q
q
g
q
)(q
)
+
2
2
2
q
q
q
]
ig [
q q
2
2
=
1
e
(q
)
i 4 e2 (q 2 ) + .
(10.138)
2
q
q
D
(q) =
We see that the one loop propagator has a divergent part in g , as well in the term
proportional to the momenta. Therefore the propagator does not correspond any
more to the Lorenz gauge and we need to add to the following terms in Lp
1
1
L2 = F F ( A )2 ,
4
2
(10.139)
2
= C
F F + ( A )
( A )2 .
4
2
2
(10.140)
As for the electron propagator, we can look at these two contributions as perturbations to the free Lagrangian, and evaluate the corresponding Feynman rules, or
evaluate the eect on the propagator. The modication in the equation dening the
propagator due to these two terms is
[(1 + C)q 2 g (C E)q q ]D (q) = ig .
(10.141)
(10.142)
q 2 (1
i
,
+ C)
i
C E
.
4
q (1 + C)(1 + E)
(10.143)
The free propagator, including the corrections at the rst order in C and E is
D (q) =
q q
ig
(1 C) i 4 (C E)
2
q
q
(10.144)
D
(q) =
q q
ig
[1 e2 (q 2 ) C] i 4 [e2 (q 2 ) + C E].
2
q
q
(10.145)
m
e2
.
C = 2 + F3
6
197
(10.146)
In fact, we are free to choose the nite term in the gauge xing, since this choice
does not change the physics. This is because the terms proportional to q q , as it
follows from gauge invariance and the consequent conservation of the electromagnetic
current, do not contribute to the physical amplitude. For instance, if we have a
vertex with a virtual photon (that is a vertex connected to an internal photon line)
and two external electrons, the term proportional to q q is saturated with
u(p ) u(p),
q = p p.
(10.147)
The result is zero, by taking into account the Dirac equation. Let us now consider
the mass of the photon. We have
]
]
ig [
q q [
2 f
2 f
1
i
F
+
e
3
3
q2
q4
ig
i q q ,
[1 ]
2
q
q4
D
(q) =
(10.148)
where
= e2 f + F3 .
(10.149)
At the second order in the electric charge we can write the propagator in the form
[
D
(q)
i
q q
= 2
g + 2 (q)
q
q (1 + (q))
(10.150)
2 ) is nite for q 2 0.
and we see that the propagator has a pole at q 2 = 0, since (q
Therefore the photon remains massless after renormalization. This is part of a rather
general aspect of renormalization which says that, if the regularization procedure
respects the symmetries of the original Lagrangian, the symmetries are preserved at
any perturbative order. However, there are cases where it is not possible to devise
a regularization procedure such to preserve a given symmetry. This the case of
the anomalous symmetries which are symmetries only at the classical level but
broken by quantum corrections.
Also for the photon we will require that at the physical pole, q 2 0, the propagator has the free form, that is
(0)
= 0,
(10.151)
from which
F3 = e2 f (0).
10.6
(10.152)
In this Section we will use the technology developed in this Chapter to evaluate two
important physical quantities. The rst one is a contribution to the Lamb shift. The
levels 2S1/2 and 2P1/2 of the hydrogen atom are degenerate except for corrections
198
coming from QED. Among these corrections there is one coming from the radiative
correction to the photon propagator and this is the one that we will evaluate here.
The second quantity that we will evaluate comes from the vertex corrections and it
is a contribution to the gyromagnetic factor of the electron which is equal to 2 in
the Dirac theory. Let us start with the vacuum polarization.
We will consider the expression (10.149) for small momenta:
d (q)
(q)
e2 q 2
,
dq 2 q2 =0
(10.153)
where use has been made of F3 being constant. Using the expression for f from
the previous Section, we get
1
(q) 2
2
f
1
m2
+ log
6 6
22
))
1 1
q2
+ 2
dz z 2 (1 z)2 2 +
2 0
m
and
(q)
It follows
e2 q 2
.
60 2 m2
(10.154)
(10.155)
g
e2 q 2
= i 2 1
+ gauge terms.
q
60 2 m2
(10.156)
The rst term, 1/q 2 , gives rise to the Coulomb potential, e2 /4r. The second gives
a correction by a term proportional to a delta function in the real space (we are
using the same notations of Section 7.2)
12 = e
dt
d4 q iq(x1 x2 ) 1
e2
e
(2)4
q 2 60 2 m2
e2
e4
3 (r).
4r 60 2 m2
(10.157)
This modication of the Coulomb potential changes the energy levels of the hydrogen
atom, and it is one of the contributions to the Lamb shift, which produces a splitting
of the levels 2S1/2 and 2P1/2 . The total Lamb shift is the sum of all the self-energy
and vertex corrections, and turns out to be about 1057.9 MHz. The contribution
we have just evaluated is only 27.1 MHz, but it is important since the agreement
between experiment and theory is at the level of 0.1 MHz.
We have now to discuss the vertex corrections. We have seen that the divergent
contribution is (1)
(p , p), and this is proportional to . The counter term to add
to the interacting part of Lp ,
Lint
p = e A ,
(10.158)
A .
(L3 )ct = eD
(10.159)
is
199
(counter-term)
(10.160)
= + T .
(10.162)
(p)
= D + e2 (p, p),
p
(10.163)
where and are the self-energy and vertex corrections evaluated from the Feynman diagrams associated to Lp . One can check that the Ward-Takahashi identity
holds for the the regularized one-loop expressions for and in equations (10.73)
and (10.103), and therefore it must be satised order by order in the expansion.
Using this result one gets
B = D,
(10.164)
at this order in perturbation theory. This equality is certainly satised by the
divergent parts of these counterterms, since they are devised in order to eliminate
the terms in 1/ in and which are equal (see eqs. (10.119) and 10.121)).
However, in order to satisfy (10.164) also the nite parts of B and D should be
equal. As a consequence we x the counter term D by
[
e2 1
+ F2 ,
D= 2
8
(10.165)
with F2 the same as in eq. (10.136). Notice that in this way the wave function
renormalization terms Z1 and Z2 are made equal, and the charge renormalization
depends only on Z3 , that is from the vacuum polarization.
200
= + e
(f
F2
2 ) .
8
(10.166)
By using again the Ward-Takahashi identity for the regularized expressions, valid at
any order in 1/ and, therefore, also for the nite parts, and the condition (10.136)
p=m
f
F2
=
,
p p=m 8 2
(10.167)
(10.168)
(10.169)
at the classical level, besides being invariant under a phase transformation. Both
symmetries are preserved by making the Dirac eld to interact with the electromagnetic eld. However, it turns out that it is impossible to dene a renormalization
procedure (choice of the counterterms) in order to satisfy both the Ward-Takahashi
identities following from the conservations of the two currents associated to the two
symmetries. Usually one makes the choice of satisfying the Ward-Takahashi identity
associated to the phase symmetry and the chiral symmetry is broken at the quantum
level (axial anomaly).
We will now evaluate the radiative corrections to the g 2 of the electron. Here g
is the gyromagnetic ratio, which is predicted to be equal to 2 by the Dirac equation
(see Section 4.6). To this end we rst need to prove the Gordon identity for the
current of a Dirac particle
[
i
p + p
u(p ) u(p) = u(p )
+
q u(p),
2m
2m
(10.170)
(10.171)
pu(p) = m u(p).
(10.172)
and
201
u(p) =
p
i
p u(p).
m m
(10.173)
An analogous operation on the barred spinor leads to the result. We observe also
that the Gordon identity shows immediately that the value of the gyromagnetic
ratio is 2, because it implies that the coupling with the electromagnetic eld is just
e
F (q).
2m
(10.174)
To evaluate the correction to this term from the one loop diagrams, it is enough to
evaluate the matrix element
e2 u(p )(2)
(p , p)u(p)
(10.175)
(10.176)
(10.177)
(10.178)
Making use of
= 2 ,
(10.179)
= 4g ,
(10.180)
p p = 2
p p,
(10.181)
we get
[
(10.182)
V = 4m
u(p ) p (y xy x2 ) + p (x xy y 2 ) u(p).
(10.183)
dx
dy
u
(p
)
p
16 2 0
m
(x + y)2
0
2 ]
(x xy y )
+p
u(p).
(10.184)
(x + y)2
e2 u(p )(2)
(p , p)u(p)
dy
=
0
0
1
1x
dx
1
y xy x2 1
1x
x
dz
dx
=
2
2
(x + y)
z
z
0
x
[(1 x) log x + x 1]
x2
x2 x2
x log x + x +
log x
+
x
2
4
2
]1
0
1
= ,
4
(10.185)
from which
e2 u(p )(2)
(p , p)u(p)
e2
u(p )[p + p ]u(p).
16 2 m
(10.186)
Using the Gordon identity in this expression, and eliminating the further contribution in , we obtain the correction to the magnetic moment
e2 u(p )(2)
(p , p)u(p)|magn.
mom.
ie2
u(p ) q u(p).
16 2 m
(10.187)
Finally we have to add this correction to the vertex part taken at p = p, which
coincides with the free vertex
ie2
q
]u(p)
p p
16 2 m
(
)
]
[
2
p
+
p
i
e
+
1 + 2 q u(p).
u(p )
2m
2m
8
u(p )[ +
(10.188)
1+
.
2m
2m
2
(10.189)
B
and e/2m, we get
Recalling that g is the ratio between S
g
=1+
+ O(2 ).
2
2
203
(10.190)
This correction was evaluated by Schwinger in 1948. Actually we know the rst
three terms of the expansion
( )2
ath
1
1
=
(g 2) =
0.32848
2
2
+ 1.49( )3 +
(10.191)
204
(10.192)