The Dirac Equation

November 25, 2020

Contents
1 Introduction 2

2 The Bjorken-Drell convention 4

3 The Probability Density and Current 5

4 Examples 7
4.1 Particle at rest . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
4.2 Particle in Motion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
4.3 Charge Conjugation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
4.4 Helicity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
4.5 Chirality . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

5 What you should know 14

6 Futhur reading 14

1
1 Introduction
In non-relativistic quantum mechanics, wave functions are descibed by the time-dependent Schrodinger
equation :
1 2 ∂ψ
− ∇ ψ+Vψ =i (1)
2m ∂t
p 2
This is really just energy conservation ( kinetic energy ( 2m ) plus potential energy (V) equals total
energy (E)) with the momentum and energy terms replaced by their operator equivalents
∂
p → −i∇; E → i (2)
∂t
In relativistic quantum theory, the energy-momentum conservation equation is E 2 −p2 = m2 (note
that we are working in the standard particle physics units where h̄ = c = 1). Proceeding with the
same replacements, we can derive the Klein-Gordon equation :
∂2ψ
E 2 − p2 − m2 → − 2
+ ∇2 ψ − m2 ψ = 0 (3)
∂t
In covariant notation this is
−∂ µ ∂µ ψ − m2 ψ = 0 (4)
Suppose ψ is solution to the Klein-Gordon equation. Multiplying it by −iψ ∗ we get
∂2ψ
iψ ∗ − iψ ∗ ∇2 ψ + iψ ∗ ψm2 = 0 (5)
∂t2
Taking the complex conjugate of the Klein-Gordon equation and multiplying by −iψ we get
∂2 ψ∗
iψ − iψ∇2 ψ ∗ + iψψ ∗ m2 = 0 (6)
∂t2
If we subtract the second from the first we obtain
∂ψ ∗
  
∂ ∗ ∂ψ
i ψ −ψ + ∇ · [−i (ψ ∗ ∇ψ − ψ∇ψ ∗ )] = 0 (7)
∂t ∂t ∂t
This has the form of an equation of continuity
∂ρ
+∇·j=0 (8)
∂t
with a probability density defined by
∂ψ ∗
 
∗ ∂ψ
ρ=i ψ −ψ (9)
∂t ∂t
and a probability density current defined by

j = i (ψ ∗ ∇ψ − ψ∇ψ ∗ ) (10)

Now, suppose a solution to the Klein-Gordon equation is a free particle with energy E and mo-
mentum p
ψ = Ne−ipµ x
µ
(11)

2
Substitution of this solution into the equation for the probability density yields

ρ = 2E|N|2 (12)

The probability density is proportional to the energy of the particle. Now, why is this a problem? If
you substitute the free particle solution into the Klein-Gordon equation you get, unsurpisingly, the
relation
E 2 − p2 = m2 (13)
so the energy of the particle could be
p
E = ± p2 + m2 (14)

The fact that you have negative energy solutions is not that much of a problem. What is a problem
is that the probability density is proportional to E. This implies a possibility for negative probability
densities...which don’t exist.
This (and some others) problem drove Dirac to think about another equation of motion. His
starting point was to try to factorise the energy momentum relation. In covariant formalism

E 2 − p2 − m2 → pµ pµ − m2 (15)

where pµ is the 4-momentum : (E, px , py , pz ). Dirac tried to write

pµ pµ − m2 = (β κ pκ + m)(γ λ pλ − m) (16)

where κ and λ range from 0 to 3. This notation looks a bit confusing. It may help if we write both
sides out long-hand. On the left-hand side we have

pµ pµ − m2 = p20 − p · p = p20 − p21 − p22 − p23 − m2 (17)

On the right-hand side we have

(β κ pκ + m)(γ λ pλ − m) = (β 0 p0 − β 1 p1 − β 2 p2 − β 3 p3 + m)(γ 0 p0 − γ 1 p1 − γ 2 p2 − γ 3 p3 − m) (18)

Expanding the right-hand side

(β κ pκ + m)(γ λ pλ − m) = β κ γ λ pκ pλ − m2 + mγ λ pλ − mβ κ pκ (19)

This should be equal to pµ pµ − m2 , so we need to get rid of the terms that are linear in p. We can do
this just be choosing β κ = γ κ . Then

pµ pµ − m2 = γ κ γ λ pκ pλ − m2 (20)

Now, since κ, λ runs from 0 to 3, the right hand side can be fully expressed by

γ κ γ λ pκ pλ − m2 = (γ 0 )2 p20 + (γ 1 )2 p21 + (γ 2 )2 p22 + (γ 3 )2 p23 (21)

+ (γ 0 γ 1 + γ 1 γ 0 )p0 p1 + (γ 0 γ 2 + γ 2 γ 0 )p0 p2 (22)
+ (γ 0 γ 3 + γ 3 γ 0 )p0 p3 + (γ 1 γ 2 + γ 2 γ 1 )p1 p2 (23)
+ (γ 1 γ 3 + γ 3 γ 1 )p1 p3 + (γ 2 γ 3 + γ 3 γ 2 )p2 p3 − m2 (24)

3
This must be equal to
pµ pµ − m2 = p20 − p21 − p22 − p23 − m2 (25)
One could make the squared terms equal by choosing (γ 0 )2 = 1 and (γ 1 )2 = (γ 2 )2 = (γ 3 )2 =-1, but
this would not eliminate the cross-terms. Dirac realised that what you needed was something which
anticommuted : i.e. γ µ γ ν + γ ν γ µ = 0 if µ 6= ν. He realised that the γ factors must be 4x4 matrices,
not just numbers, which satisfied the anticommutation relation

{γ µ , γ ν } = 2g µν (26)
where  
1 0 0 0
0 −1 0 0
g µν =
0 0 −1 0 
 (27)
0 0 0 −1
We also define, for later use, another γ matrix called γ 5 = iγ 0 γ 1 γ 2 γ 3 . This has the properties (you
should check for yourself) that

(γ 5 )2 = 1, {γ 5 , γ µ } = 0 (28)
We will want to take the Hermitian conjugate of these matrices at various times. The Hermitian
conjugate of each matrix is

γ 0† = γ 0 γ 5† = γ 5 γ µ† = γ 0 γ µ γ 0 = −γ for µ 6= 0 (29)
If the matrices satisfy this algebra, then you can factor the energy momentum conservation equa-
tion
pµ pµ − m2 = (γ κ pκ + m)(γ λ pλ − m) = 0 (30)
The Dirac equation is one of the two factors, and is conventionally taken to be

γ λ pλ − m = 0 (31)
Making the standard substitution, pµ → i∂µ we then have the usual covariant form of the Dirac
equation

(iγ µ ∂µ − m)ψ = 0 (32)

∂ ∂ ∂ ∂
where ∂µ = ( ∂t , ∂x , ∂y , ∂z ), m is the particle mass and the γ matrices are a set of 4-dimensional
matrices. Since these are matrices, ψ is a 4-element column matrix called a “bi-spinor”. This bi-spinor
is not a 4-vector and doesn’t transform like one.

2 The Bjorken-Drell convention

Any set of four 4x4 matrices that obey the algebra above will work. The one we normally use includes
the Pauli spin matrices. Recall that each component of the spin operator S for spin 1/2 particles is
defined by its own Pauli spin matrix :
     
1 1 0 1 1 1 0 −i 1 1 1 0
Sx = σ1 = Sy = σ2 = Sz = σ3 = (33)
2 2 1 0 2 2 i 0 2 2 0 −1

4
 In terms of the Pauli spin matrices, the γ matrices have the form
     
0 1 0 µ 0 σµ 5 0 1
γ = γ = γ = (34)
0 −1 −σ µ 0 1 0
Each element of the matrices in Equations 34 are 2x2 matrices. 1 denotes the 2 x 2 unit matrix,
and 0 denotes the 2 x 2 null matrix.

3 The Probability Density and Current

In order to understand the probability density and probability flow we will want to derive an equation
of continuity for the probability. The first step is to write the Dirac equation out longhand :
∂ψ ∂ψ ∂ψ ∂ψ
iγ 0
+ iγ 1 + iγ 2 + iγ 3 − mψ = 0 (35)
∂t ∂x ∂y ∂z
We want to take the Hermitian conjugate of this :
∂ψ ∂ψ ∂ψ ∂ψ
+ iγ 1
[iγ 0 + iγ 2 + iγ 3 − mψ]† (36)
∂t ∂x ∂y ∂z
Now, we must remember that the γ are matrices and that ψ is a 4-component column vector.
Thus

   ∂ψ1 
1 0 0 0 ∂t
∂ψ † 0 1 0 0  ∂ψ2  †
[γ 0 ] = [   ∂t ] (37)
∂t 0 0 −1 0   ∂ψ∂t
3

0 0 0 −1 ∂ψ 4
∂t
 ∂ψ1 †  †
∂t 1 0 0 0
 ∂ψ2  0 1 0 0
=  ∂t  
 ∂ψ3  0
 (38)
∂t
0 −1 0 
∂ψ4 0 0 0 −1
∂t
 
1 0 0 0
∂ψ4† 0 1 0 0
 
∂ψ1† ∂ψ2† ∂ψ3†
=  (39)
∂t ∂t ∂t ∂t  0 0 −1 0 
0 0 0 −1
∂ψ † 0
= γ (40)
∂t
Using the Hermitian conjugation properties of the γ matrices defined in the previous section we
can write Equation 36 as

∂ψ † 0† ∂ψ † 1† ∂ψ † 2† ∂ψ† 3†
−i γ −i γ −i γ −i γ − mψ † (41)
∂t ∂x ∂y ∂z
which, as γ µ† = −γ µ for µ 6= 0 means we can write this as

∂ψ † 0 ∂ψ † ∂ψ † ∂ψ†
−i γ −i (−γ 1 ) − i (−γ 2 ) − i (−γ 3 ) − mψ † (42)
∂t ∂x ∂y ∂z

5
 We have a problem now - this is no longer covariant. That is, I would like to write this in the form
← ←
ψ † (−i ∂ µ γ µ − m) where the derivative ∂ µ now operates to the left, as indicated by the arrow above
the partial derivative symbol1 I can’t because the minus signs on the spatial components coming for
the −γ µ spoils the scalar product.
To deal with this we have to multiply the equation on the right by γ 0 , since we can flip the sign
of the γ matrices using the relationship −γ µ γ 0 = γ 0 γ µ . Thus

∂ψ † 0 0 ∂ψ † 1 0 ∂ψ † ∂ψ†
−i γ γ −i (−γ γ ) − i (−γ 2 γ 0 ) − i (−γ 3 γ 0 ) − mψ † γ 0 (43)
∂t ∂x ∂y ∂z
or

∂ψ † 0 0 ∂ψ † 0 1 ∂ψ † 0 2 ∂ψ† 0 3
−i γ γ −i (γ γ ) − i (γ γ ) − i (γ γ ) − mψ † γ 0 (44)
∂t ∂x ∂y ∂z
Defining the adjoint spinor ψ = ψ † γ 0 , we can rewrite this equation as

∂ψ 0 ∂ψ ∂ψ ∂ψ
−i γ − i (γ 1 ) − i (γ 2 ) − i (γ 3 ) − mψ (45)
∂t ∂x ∂y ∂z
and finally we get the adjoint Dirac equation
←
ψ(i ∂ µ γ µ + m) = 0 (46)

Now, we multiply the Dirac equation on the left by ψ :

ψ(iγ µ ∂µ − m)ψ = 0 (47)

and the adjoint Dirac equation on the right by ψ :

←
ψ(i ∂ µ γ µ + m)ψ = 0 (48)

and we add them together, which eliminates the two mass terms
←
ψ(γ µ ∂µ ψ) + (ψ ∂ µ γ µ )ψ = 0 (49)

or, more succinctly,

∂µ (ψγ µ ψ) = 0 (50)
.
If we now identify the bit in the brackets as the 4-current :

j µ = (ρ, j) (51)

where ρ is the probability density and j is the probability current, we can write the conservation
equation as
∂µ j µ = 0 with j µ = ψγ µ ψ (52)
which is the covariant form for an equation of continuity. The probability density is then just ρ =
ψγ 0 ψ = ψ † γ 0 γ 0 ψ = ψ † ψ and the probability 3-current is j = ψ † γ 0 γ µ ψ.

6
 Figure 1: Feynman diagram for the electromagnetic scattering of one fermion off another.

This current is the same one which appears in the Feynman diagrams. It is called a Vector current,
and is the current responsible for the electromagnetic interaction.
For the interaction in Figure 1, with two electromagnetic interactions, the matrix element is then

e2 1
M= [ψf γ µ ψi ] 2 [φf γ ν φi ] (53)
4π q

4 Examples
Now let’s look at some solutions to the Dirac Equation. The first one we will look at is for a particle
at rest.

4.1 Particle at rest

Suppose the fermion wavefunction is a plane wave. We can write this as

ψ(p) = e−ixµ p u(p)

µ
(54)

where pµ = (E, −px , −py , −pz ) and xµ = (t, x, y, z) and so −ixµ pµ = −i(Et − x · p). This is just the
phase for an oscillating plane wave.
For a particle at rest, the momentum term disappears : −ixµ pµ = −i(Et). Furthur, since the
∂ψ
momentum is zero, the spatial derivatives must be zero : ∂x,y,z = 0. The Dirac equation therefore
reads
1
This notation is fairly standard and is there to reflect the symmetry of ψ and ψ in the problem. arrow above the
operator to show which direction it acts. The notation comes from quantum mechanics. For example we define the
operator Ô operating on a ket |v > to be some vector Ô|v >. The operator acting to the left on a bra < a| : < a|Ô is
defined to be that function which, when acting on a ket |v > returns the number < a|Ô|v >

7
 ∂ψ
iγ 0 − mψ = 0 (55)
∂t
or
iγ 0 (−iE)u(p) − mu(p) = 0 (56)
which gives us
Eγ 0 u(p) = mu(p) (57)
Now, u(p) is a 4-component bispinor, so properly expanding the gamma matrix in Equation 57,
we get    
1 0 0 0 u1
0 1 0 0 u2 
  
E0 0 −1 0  = m u3  (58)
0 0 0 −1 u4
This is an eigenvalue equation. There are four independent solutions. Two with energy E = m
and two with E = −m. The solutions are just
       
1 0 0 0
0 1 0 0
0 u2 = 0 u3 = 1 u4 = 0
u1 =         (59)
0 0 0 1

with eigenvalues m, m, -m and -m respectively.

The first two solutions can be interpreted as positive energy particle solutions with spin up and
spin down. One can see this be checking if the function is an eigenfunction of the spin matrix : Sz .
  
1 0 0 0 1
1 0 −1 0 0  0
  
 = 1 u1
Sz u 1 =  (60)
2 0 0 1 0   0 2
0 0 0 −1 0
and similarly for u2 .
Note that we have yet to normalise these solutions. We won’t either for the purposes of this
discussion. What about these negative energy solutions? We are required to keep them (we can’t just
say that they are unphysical and throw them away) since quantum mechanics requires a complete set
of states.
What happens if we allow the particle to have momentum?

4.2 Particle in Motion

Again, let’s look for a plane wave solution

ψ(p) = e−ixµ p u(p)

µ
(61)

Plugging this into the Dirac equation means that we can replace the ∂µ by pµ and remove the
oscillatory term to obtain the momentum-space version of the Dirac Equation

(γ µ pµ − m)u(p) = 0 (62)

8
 Notice that this is now purely algebraic and can be easily(!) solved

γ µ pµ − m = Eγ 0 − px γ 1 − py γ 2 − pz γ 3 − m (63)
     
1 0 0 σ 1 0
= E− ·p−m (64)
0 −1 −σ 0 0 1
 
(E − m) −σ · p
= (65)
σ·p −(E + m)

where each element in this 2x2 matrix is actually a 2x2 matrix itself. In this spirit, let’s write the
4-component bispinor solution as 2-component vector
 
uA
u= (66)
uB

then the Dirac Equation implies that

    
µ (E − m) −σ · p uA 0
(γ pµ − m)u(p) = 0 =⇒ = (67)
σ·p −(E + m) uB 0

leading two coupled simultaneous equations

(σ · p)uB = (E − m)uA (68)

(σ · p)uA = (E + m)uB (69)

Now, let’s expand that (σ · p) :

     
0 1 0 −1 1 0
(σ · p) = p + py + p (70)
1 0 x i 0 0 −1 z
 
pz px − ipy
= (71)
px + ipy −pz

and right back to the Dirac Equation

(σ · p)uB = (E − m)uA (72)

(σ · p)uA = (E + m)uB (73)

gives

(σ · p)
uB = uA (74)
E+m 
1 pz px − ipy
= uA (75)
E + m px + ipy −pz

Now, we just need to make choices for the form of uA . Let’s make the obvious choice and remember
that uA is a 2-component spinor so we need to specify two solutions:
   
1 0
uA = or uA = (76)
0 1

9
 These give
1 0
   
 0   1 
u1 = N1  pz 
 E+m  and u2 = N2 
 px −ipy 
 (77)
E+m
px +ipy −pz
E+m E+m

where N1 and N2 are normalisation factors we won’t go into. Note that, if p = 0 these correspond to
the E > 0 solutions we found
 for
 a particle 
at rest,
 which is good.
1 0
Repeating this for uB = and uB = which gives solutions u3 and u4
0 1
 pz   px −ipy 
E−m E−m
 px +ipy   −pz 
u3 = N3  E−m 
 1  and u4 = N4  E−m 
 0  (78)
0 1
and collecting them together
pz   px −ipy 
1 0
    
E−m E−m
 0   1   px +ipy   −pz 
u1 = N1  p z
 u2 = N2  px −ip  u3 = N3 
 E−m  u4 = N4  E−m  (79)
1   0 
 E+m   y
E+m
px +ipy −pz
E+m E+m 0 1

All of which, if put back into the Dirac Equation, yields : E 2 = p2 + m2 as you might expect.
Comparing
p with the p = 0 case we can identify u1 and u2 as positive energy p solutions with energy
E = p2 + m2 , and u3 , u4 as negative energy solutions with energy E = − p2 + m2 .
How do we interpret these negative energy solutions? The conventional interpretation is called
the “Feynman-Stuckelberg” interpretation : A negative energy solution represents a negative energy
particle travelling backwards in time, or equivalently, a positive energy antiparticle going forwards in
time.
With this interpretation we tend to rewrite the negative energy solutions to represent positive
antiparticles. Starting from the E < 0 solutions
 pz   px −ipy 
E−m E−m
 px +ipy   −pz 
u3 = N3  E−m 
 1  and u4 = N4  E−m 
 0  (80)
0 1

Here we implicitly assume that E is negative. We define antiparticle states by just flipping the sign
of E and p following the Feynman-Stuckelberg convention.

v1 (E, p)e−i(Et−x·p) = u4 (−E, −p)ei(Et−x·p) (81)

v2 (E, p)e−i(Et−x·p) = u3 (−E, −p)ei(Et−x·p) (82)
p
where E = |p|2 + m2 . Note that v1 is associated with u4 and v2 is associated with u3 . We do
this because the spin matrix, Santiparticle , for anti-particles is equal to −Sparticle i.e. the spin flips for
antiparticles as well and the spin eigenvalue for v1 = u4 is spin-up and v2 = u3 is spin-down.
In general u1, u2 , v1 and v2 are not eigenstates of the spin operator (Check for yourself). In
fact we should expect this since solutions to the Dirac Equation are, by definition, eigenstates of the

10
Hamiltonian operator, Ĥ, and the Sz does not commute with the Hamiltonian : [Ĥ, Ŝz ] 6= 0, and
hence we can’t find solutions which are simultaneously solutions to Sz and Ĥ. However if the z-axis
is aligned with particle direction : px = py = 0, pz = ±|p| then we have the following Dirac states
       ±|p| 
1 0 0 E+m
 0   1   ∓|p|   0 
u1 = N  ±|p|  u2 = N  0  v1 = N  0  v2 = N  1 
    E+m    (83)
E+m
∓|p|
0 E+m 1 0

These are eigenstates of Sz

1 1
Sz u 1 = + u 1 Szanti v1 = −Sz v1 = + v1 (84)
2 2
1 1
Sz u 2 = − u 2 Szanti v2 = −Sz v2 = − v2 (85)
2 2
So we have found solutions of the Dirac Equation which are also spin eigenstates....but only if the
particle is travelling along the z-axis.

4.3 Charge Conjugation

Classical electrodynamics is invariant under a change in the sign of the electric charge. In particle
physics, this concept is represented by the “charge conjugation operator” that flips the signs of all
the charges. It changes a particle into an anti-particle, and vice versa:

Ĉ|p >= |p > (86)

. Application of Ĉ twice brings us back to the the orignal state : Ĉ 2 = 1 and so eigenstates of Ĉ are
±1. In general most particles are not eigenstates of Ĉ. If |p > were an eigenstate of Ĉ then

Ĉ|p >= ±|p >= |p > (87)

implying that only those particles which are their own antiparticles are eigenstates of Ĉ. This leaves
us with photons and the neutral mesons like π 0 , η and ρ0 .
Ĉ changes a particle spinor into an anti-particle spinor. The relevant operation is

ψ → ψC = Ĉγ 0 ψ ∗ (88)

with Ĉ = iγ 2 γ 0 .
We can apply this to, say, the u1 solution to the Dirac Equation. if ψ = u1 e−i(Et−x·p) , the
ψC = iγ 2 ψ ∗ = iγ 2 u∗1 ei(Et−x·p)
   ∗  px −ipy 
0 0 0 −i 1 E+m
2 ∗
0 0 i 0   0   −pz 
iγ u1 = i   pz  = N1  E+m  = v1
0 −i 0 0  N1  E+m

  0  (89)
px +ipy
i 0 0 0 E+m 1
or, in full, Ĉ(u1 e−i(Et−x·p) ) = v1 ei(Et−x·p) .
Hence charge conjugation changes the particle eigenspinors into their respective anti-particle
spinors.

11
4.4 Helicity
The fact that we can find spin eigenvalues for states in which the particles are travelling along the
spin-direction indicates that the quantity we need is not spin but helicity. The helicity is defined as
the projection of the spin along the direction of motion:
 
σ 0
ĥ = Σ · p̂ = 2S · p̂ = · p̂ (90)
0 σ
and has eigenvalues equal to +1 (called right-handed where the spin vector is aligned in the same
direction as the momentum vector) or -1 (called left-handed where the spin vector is aligned in the
opposite direction as the momentum vector), corresponding to the diagrams in Figure2.

Figure 2: The two helicity states. On the left the spin vector is aligned in the direction of motions of
the particle. This is a right-handed helicity state and has helicity +1. On the right the spin vector
is antiparallel to the particle momentum. This is a left-handed helicity state with helicity -1.

It can be shown that the helicity does commute with the Hamiltonian and so one can find eigen-
states that are simultaneously states of helicity and the Hamiltonian.
The problem, and it is a big problem, is that helicity is not Lorentz invariant in the case of a
massive particle. If the particle is massive it is possible to find an inertial reference frame in which
the particle is going in the opposite direction. This does not change the direction of the spin vector,
so the helicity can change sign.
The helicity is Lorentz invariant only in the case of massless particles.

4.5 Chirality
We’d rather have operators which are Lorentz invariant, than commute with the Hamiltonian. In
general wave functions in the Standard Model are eigenstates of a Lorentz invariant quantity called
the chirality. The chirality operator is γ 5 and it does not commute with the Hamiltonian. Due to this,
it is somewhat difficult to visualise. The best picture I can get comes from the definition : something
is chiral if it’s mirror image does not superimpose on itself. Think of your left hand. It’s mirror image
(from the point of view of someone in that mirror looking back at you) is actually a right hand. It
and it’s mirror image cannot be superimposed andit is therefore an intrinsically chiral object.
In the limit that E >> m, or that the particle is massless, the chirality is identical to the helicity.
For a massive particle this is no longer true.
In general the eigenstates of the chirality operator are
γ 5 uR = +uR γ 5 uL = −uL γ 5 vR = −vR γ 5 vL = +vL (91)
, where we define uR and uL are right- and left-handed chiral states. We can define the projection
operators
1 1
PL = (1 − γ 5 ) PR = (1 + γ 5 ) (92)
2 2
12
such that PL projects outs the left-handed chiral particle states and right-handed chiral anti-particle
states. PR projects out the right-handed chiral particle states and left-handed chiral anti-particle
states. The projection operators have the following properties :
2
PL,R = PL,R ; PR PL = PL PR = 0; PR + PL = 1 (93)

Any spinor can be written in terms of it’s left- and right-handed chiral states:

ψ = (PR + PL )ψ = PR ψ + PL ψ = ψR + ψL (94)

Chirality holds an important place in the standard model. Let’s take a standard model probability
current
ψγ µ φ (95)
We can decompose this into it’s chiral states

ψγ µ φ = (ψL + ψR )γ µ (φL + φR )
= ψL γ µ φL + ψL γ µ φR + ψR γ µ φL + ψR γ µ φR

Now, ψL = ψL† γ 0 = ψ † 12 (1 − γ 5 )γ 0 = ψ † γ 0 12 (1 + γ 5 ) = ψ 21 (1 + γ 5 ) = ψPR where I have used the

fact that γ 5† = γ 5 and γ 5 γ 0 = −γ 0 γ 5 . Similarly ψR = ψPL .
Using this, it is easy to show that the terms ψL γ µ φR and ψR γ µ φL are equal to 0:

ψL γ µ φR = ψPR γ µ PR φ (96)
1 1
= ψ (1 + γ 5 )γ µ (1 + γ 5 )φ (97)
2 2
1
= ψ (1 + γ )(γ + γ µ γ 5 )φ
5 µ
(98)
4
1
= ψ (1 + γ 5 )(γ µ − γ 5 γ µ )φ (99)
4
1
= ψ (1 + γ 5 )(1 − γ 5 )γ µ φ (100)
4
1
= ψ (1 + γ 5 − γ 5 − (γ 5 )2 )γ µ φ (101)
4
= 0 (102)

since (γ 5 )2 = 1 and γ µ γ 5 = −γ 5 γ µ . Similarly for the other cross term, ψR γ µ φL .

Hence,

ψγ µ φ = ψL γ µ φL + ψR γ µ φR (103)

So left-handed chiral particles couple only to left-handed chiral fields, and right-handed chiral fields
couple to right-handed chiral fields.
One must be very careful with how one interprets this statement. What it does not say is that
there are left-handed chiral electrons and right-handed chiral electrons which are distinct particles.
The “particles” are those states which propagate with fixed mass under the Dirac equation. What this
means is that the Standard Model does not mix chiralities or, to put it another way, no fundamental
interaction is capable of turning a left-chiral particle into a right-chiral particle. At the same time, it

13
also shows that the electromagnetic interaction (embodied in the vector current) treats both chiralities
the same. In this case, this part of the Standard Model exhibits chiral symmetry. We will find that
this is not the case for the weak interaction.
Useful though it is when describing the interaction of fields in the Standard Model, chirality is
not conserved in the propagation of a free particle. In fact the chiral states φL and φR do not even
satisfy the Dirac equation. Since chirality is not a good quantum number it can evolve with time. A
massive particle starting off as a completely left handed chiral state can evolve a right-handed chiral
component. By contrast, helicity is a conserved quantity during free particle propagation. Only in
the case of massless particles, for which helicity and chirality are identical and are conserved in free-
particle propagation, can left- and right-handed particles be considered distinct. For neutrinos this
mostly holds.

5 What you should know

• Understand the covariant form of the Dirac equation, and know how to manipulate the γ
matrices.

• Know where the definition of a current comes from.

• Understand the difference between helicity and chirality.

• Be able to manipulate gamma matrices.

6 Futhur reading
Griffiths Chapter 7, Sections 7.1, 7.2 and 7.3 Griffiths Chapter 10, Section 10.7

14