5Quantizing the Dirac field
III Quantum Field Theory
5.1 Fermion quantization
Quantization
We now start to quantize our spinor field! As it will turn out, this is complicated.
Recall that for a complex scalar field, classically a solution of momentum
p
can be written as
b
p
e
−ip·x
+ c
∗
p
e
ip·x
for some constants
b
p
, c
p
. Here the first term is the positive frequency solution,
and the second term is the negative frequency solution. To quantize this field,
we then promoted the b
p
, c
p
to operators, and we had
φ(x) =
Z
d
3
p
(2π)
3
1
p
2E
p
(b
p
e
ip·x
+ c
†
p
e
−ip·x
).
Similarly, classically, every positive frequency solution to Dirac’s equation of
momentum p can be written as
(b
1
p
u
1
p
+ b
2
p
u
2
p
)e
−ip·x
for some b
s
p
, and similarly a negative frequency solution can be written as
(c
1
p
v
1
p
+ c
2
p
v
2
p
)e
ip·x
for some
c
s
p
. So when we quantize our field, we could promote the
b
p
and
c
p
to
operators and obtain
ψ(x) =
2
X
s=1
Z
d
3
p
(2π)
3
1
p
2E
p
b
s
p
u
s
p
e
ip·x
+ c
s†
p
v
s
p
e
−ip·x
ψ
†
(x) =
2
X
s=1
Z
d
3
p
(2π)
3
1
p
2E
p
b
s†
p
u
s†
p
e
−ip·x
+ c
s
p
v
s†
p
e
ip·x
In these expressions,
b
s
p
and
c
s
p
are operators, and
u
s
p
and
v
s
p
are the elements of
the spinor space we have previously found, which are concrete classical objects.
We can compute the conjugate momentum to be
π =
∂L
∂
˙
ψ
= i
¯
ψγ
0
= iψ
†
.
The Hamiltonian density is then given by
H = π
˙
ψ − L
= iψ
†
˙
ψ − i
¯
ψγ
0
˙
ψ − i
¯
ψγ
i
∂
i
ψ + m
¯
ψψ
=
¯
ψ(−iγ
i
∂
i
+ m)ψ.
What are the appropriate commutation relations to impose on the
ψ
, or, equiva-
lently, b
s
p
and c
s†
p
? The naive approach would be to impose
[ψ
α
(x), ψ
β
(y)] = [ψ
†
α
(x), ψ
†
β
(y)] = 0
and
[ψ
α
(x), ψ
†
β
(y)] = δ
αβ
δ
3
(x − y).
These are in turn equivalent to requiring
[b
r
p
, b
s†
q
] = (2π)
3
δ
rs
δ
3
(p − q),
[c
r
p
, c
s†
q
] = −(2π)
3
δ
rs
δ
3
(p − q),
and all others vanishing.
If we do this, we can do the computations and find that (after normal
ordering) we obtain a Hamiltonian of
H =
Z
d
3
p
(2π)
3
E
p
(b
s†
p
b
s
p
− c
s†
p
c
s
p
).
Note the minus sign before
c
s†
p
c
s
p
. This is a disaster! While we can define a
vacuum
|0i
such that
b
s
p
|0i
=
c
s
p
|0i
= 0, this isn’t really a vacuum, because we
can keep applying c
†
p
to get lower and lower energies. This is totally bad.
What went wrong? The answer comes from actually looking at particles in
the Real World
TM
. In the case of scalar fields, the commutation relation
[a
†
p
, a
†
q
] = 0.
tells us that
a
†
p
a
†
q
|0i = |p, qi = |q, pi.
This means the particles satisfy Bose statistics, namely swapping two particles
gives the same state.
However, we know that fermions actually satisfy Fermi statistics. Swapping
two particles gives us a negative sign. So what we really need is
b
r
p
b
s
q
= −b
s
q
b
r
p
.
In other words, instead of setting the commutator zero, we want to set the
anticommutator to zero:
{b
r
p
, b
s
q
} = 0.
In general, we are going to replace all commutators with anti-commutators. It
turns out this is what fixes our theory. We require that
Axiom. The spinor field operators satisfy
{ψ
α
(x), ψ
β
(y)} = {ψ
†
α
(x), ψ
†
β
(y)} = 0,
and
{ψ
α
(x), ψ
†
β
(y)} = δ
αβ
δ
3
(x − y).
Proposition. The anti-commutation relations above are equivalent to
{c
r
p
, c
s†
q
} = {b
r
p
, b
s†
q
} = (2π)
3
δ
rs
δ
3
(p − q),
and all other anti-commutators vanishing.
Note that from a computational point of view, these anti-commutation
relations are horrible. When we manipulate our operators, we will keep on
introducing negative signs to our expressions, and as we all know, keeping track
of signs is the hardest thing in mathematics. Indeed, these negative signs will
come and haunt us all the time, and we have to insert funny negative signs
everywhere.
If we assume these anti-commutation relations with the same Hamiltonian
as before, then we would find
Proposition.
H =
Z
d
3
p
(2π)
3
E
p
b
s†
p
b
s
p
+ c
s†
p
c
s
p
.
We now define the vacuum as in the bosonic case, where it is annihilated by
the b and c operators:
b
s
p
|0i = c
s
p
|0i = 0.
Although the
b
’s and
c
’s satisfy the anti-commutation relations, the Hamiltonian
satisfies commutator relations
[H, b
r†
p
] = Hb
r†
p
− b
r†
p
H
=
Z
d
3
q
(2π)
3
E
q
(b
s†
q
b
s
q
+ c
s†
q
c
s
q
)b
r
†
p
− b
r†
p
Z
d
3
q
(2π)
3
E
q
(b
s†
q
b
s
q
+ c
s†
q
c
s
q
)
= E
p
b
r†
p
.
Similarly, we have
[H, b
r
p
] = −E
p
b
r
p
and the corresponding relations for c
r
p
and c
r†
p
.
Heisenberg picture
As before, we would like to put our Dirac field in the Heisenberg picture. We
have a spinor operator at each point x
µ
ψ(x) = ψ(x, t).
The time evolution of the field is given by
∂ψ
∂t
= i[H, ψ],
which is solved by
ψ(x) =
2
X
s=1
Z
d
3
p
(2π)
3
1
p
2E
p
b
s
p
u
s
p
e
−ip·x
+ c
s†
p
v
s
p
e
ip·x
ψ
†
(x) =
2
X
s=1
Z
d
3
p
(2π)
3
1
p
2E
p
b
s†
p
u
s†
p
e
ip·x
+ c
s
p
v
s†
p
e
−ip·x
.
We now look at the anti-commutators of the fields in the Heisenberg picture.
We define
Definition.
iS
αβ
(x − y) = {ψ
α
(x),
¯
ψ
β
(y)}.
Dropping the indices, we write
iS(x − y) = {ψ(x),
¯
ψ(y)}.
If we drop the indices, then we have to remember that
S
is a 4
×
4 matrix with
indices α, β.
Substituting ψ and
¯
ψ in using the expansion, we obtain
iS(x − y) =
X
s,r
Z
d
3
p d
3
q
(2π)
6
1
p
4E
p
E
q
{b
s
p
, b
r†
q
}u
s
p
¯u
r
q
e
−i(p·x−q·y)
+ {c
s†
p
, c
r
q
}v
s
p
¯v
r
q
e
i(p·x−q·y)
=
X
s
Z
d
3
p
(2π)
3
1
2E
p
u
s
p
¯u
s
p
e
−p·(x−y)
+ v
s
p
¯v
s
p
e
ip·(x−y)
=
Z
d
3
p
(2π)
3
1
2E
p
h
(
/
p + m)e
−ip·(x−y)
+ (
/
p − m)e
ip·(x−y)
i
Recall that we had a scalar propagator
D(x −y) =
Z
d
3
p
(2π)
3
1
2E
p
e
−ip·(x−y)
.
So we can write the above quantity in terms of the D(x − y) by
iS(x − y) = (i
/
∂
x
+ m)(D(x −y) − D(y − x)).
Some comments about this: firstly, for spacelike separations where
(x − y)
2
< 0,
we had D(x − y) −D(y − x) = 0.
For bosonic fields, we made a big deal of this, since it ensured that
[
φ
(
x
)
, φ
(
y
)] = 0. We then interpreted this as saying the theory was causal.
For spinor fields, we have anti-commutation relations rather than commutation
relations. What does this say about causality? The best we can say is that all
our observables are bilinear (or rather quadratic) in fermions, e.g.
H =
Z
d
3
x ψ
†
(−iγ
i
∂
i
+ m)ψ.
So these observables will commute for spacelike separations.
Propagators
The next thing we need to figure out is the Feynman propagator. By a similar
computation to that above, we can determine the vacuum expectation value
h0|ψ(x)
¯
ψ(y) |0i =
Z
d
3
p
(2π)
3
1
2E
p
(
/
p + m)e
−ip·(x−y)
h0|
¯
ψ(y)ψ(x) |0i =
Z
d
3
p
(2π)
3
1
2E
p
(
/
p − m)e
ip·(x−y)
.
We define the Feynman propagator by
Definition
(Feynman propagator)
.
The Feynman propagator of a spinor field
is the time-ordered product
S
F
(x − y) = h0|T ψ
α
(x)
¯
ψ
β
(y) |0i =
(
h0|ψ
α
(x)
¯
ψ
β
(y) |0i x
0
> y
0
−h0|
¯
ψ
β
(y)ψ
α
(x) |0i y
0
> x
0
Note that we have a funny negative sign! This is necessary for Lorentz
invariance — when (
x − y
)
2
<
0, then there is no invariant way to determine
whether one time is bigger than the other. So we need the expressions for the two
cases to agree. In the case of a boson, we had
h0|φ
(
x
)
φ
(
y
)
|0i
=
h0|φ
(
y
)
φ
(
x
)
|0i
,
but here we have anti-commutation relations, so we have
ψ(x)
¯
ψ(y) = −
¯
ψ(y)ψ(x).
So we need to insert the negative sign to ensure time-ordered product as defined
is Lorentz invariant.
For normal ordered products, we have the same behaviour. Fermionic opera-
tors anti-commute, so
:ψ
1
ψ
2
: = −:ψ
2
ψ
1
: .
As before, the Feynman propagator appears in Wick’s theorem as the contraction:
Proposition.
ψ(x)
¯
ψ(y) = T (ψ(x)
¯
ψ(y)) −:ψ(x)
¯
ψ(y): = S
F
(x − y).
When we apply Wick’s theorem, we need to remember the minus sign. For
example,
:ψ
1
ψ
2
ψ
3
ψ
4
: = −:ψ
1
ψ
3
ψ
2
ψ
4
: = −ψ
1
ψ
3
:ψ
2
ψ
4
: .
Again, S
F
can be expressed as a 4-momentum integral
S
F
(x − y) = i
Z
d
4
p
(2π)
4
e
−ip·(x−y)
/
p + m
p
2
− m
2
+ iε
.
As in the case of a real scalar field, it is a Green’s function of Dirac’s equation:
(i
/
∂
x
− m)S
F
(x − y) = iδ
4
(x − y).