LectureNotes SM and Higgs
LectureNotes SM and Higgs
Preliminaries i
Introduction 1
i.1 Quantum fields . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
i.2 The Yukawa interaction . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
i.3 Feynman diagrams . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
i.4 The Standard Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
i.5 Units in particle physics . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
i.6 Four-vector notation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
i
ii CONTENTS
6 Spin-1/2 Electrodynamics 93
6.1 Feynman rules for fermion scattering . . . . . . . . . . . . . . . . . . . . 93
6.2 Electron-muon scattering . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
6.3 Crossing: the process e− e+ → µ− µ+ . . . . . . . . . . . . . . . . . . . . . 101
6.4 Summary of QED Feynman rules . . . . . . . . . . . . . . . . . . . . . . 104
Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
These are the lecture notes for the Particle Physics 1 (PP1) master course that is
taught at Nikhef in the autumn semester of 2014. These notes contain 14 chapters, each
corresponding to one lecture session. The topics discussed in this course are:
• Lecture 1 - 4: Electrodynamics of spinless particles
• Lecture 5 - 6: Electrodynamics of spin 1/2 particles
• Lecture 7: The weak interaction
• Lecture 8 - 10: Gauge symmetries and the electroweak theory
• Lecture 11-14: Electroweak symmetry breaking
Each lecture of 2 × 45 minutes is followed by a 1.5 hour problem solving session. The
exercises are included in these notes, at the end of each chapter.
The notes mainly follow the material as discussed in the books of Halzen and Martin.
The first ten chapters have been compiled by Marcel Merk in the period 2000-2011, and
updated by Wouter Hulsbergen for the PP1 courses of 2012 and 2013. The last four
chapters, written by Ivo van Vulpen, were added in 2014.
Literature
The following is a non-exhaustive list of course books on particle physics. (The comments
reflect a personel opinion of your lecturers!)
Thomson: “Modern Particle Physics”:
This is a new book (2013) that covers practically all the material in these lectures. If
you do not have another particle physics book yet, then we recommend that you acquire
this book.
Halzen & Martin: “Quarks & Leptons: an Introductory Course in Modern Particle
Physics ”:
This is the book that your lecturers used when they did their university studies. Though
most of the theory is timeless, it is a bit outdated when it comes to experimental results.
i
ii PRELIMINARIES
The book builds on earlier work of Aitchison (see below). Most of the course follows
this book, but it is no longer in print.
Griffiths: “Introduction to Elementary Particle Physics”, second, revised ed.
The text is somewhat easier to read than H & M and is more up-to-date (2008) (e.g.
neutrino oscillations) but on the other hand has a somewhat less robust treatment in
deriving the equations. The introduction chapter of this book gives a very readable
popular history of particle physics.
Aitchison & Hey: “Gauge Theories in Particle Physics”
Meanwhile in its 4th edition(2012), this 2-volume book provides a thorough theoreti-
cal introduction to particle physics, including field theory. It is excellent (notably its
’comments’ and appendices), but a bit more formal than needed for this course.
Perkins: “Introduction to High Energy Physics”, (1987) 3-rd ed., (2000) 4-th ed.
The first three editions were a standard text for all experimental particle physics. It is
dated, but gives an excellent description of, in particular, the experiments. The fourth
edition is updated with more modern results, while some older material is omitted.
Aitchison: “Relativistic Quantum Mechanics”
(1972) A classical, very good, but old book, often referred to by H & M.
Burcham & Jobes: “Nuclear & Particle Physics”
(1995) An extensive text on nuclear physics and particle physics. It contains more
(modern) material than H & M. Formula’s are explained rather than derived and more
text is spent to explain concepts.
Das & Ferbel: “Introduction to Nuclear and Particle Physics”
(2006) A book that is half on experimental techniques and half on theory. It is more
suitable for a bachelor level course and does not contain a treatment of scattering theory
for particles with spin.
Martin and Shaw: “Particle Physics ”, 2-nd ed.
(1997) A textbook that is somewhere inbetween Perkins and Das & Ferbel. In my
opinion it has the level inbetween bachelor and master.
Particle Data Group: “Review of Particle Physics”
This book appears every two years in two versions: the book and the booklet. Both of
them list all aspects of the known particles and forces. The book also contains concise,
but excellent short reviews of theories, experiments, accellerators, analysis techniques,
statistics etc. There is also a version on the web: http://pdg.lbl.gov
The Internet:
In particular Wikipedia contains a lot of information. However, one should note
that Wikipedia does not contain original articles and they are certainly not re-
viewed! This means that they cannot be used for formal citations.
In addition, have a look at google books, where (parts of) books are online avail-
able.
iii
About Nikhef
Nikhef is the Dutch institute for subatomic physics. (The name was originally an
acronym for ”Nationaal Instituut voor Kern en Hoge Energie Fysica”.) The name
Nikhef is used to indicate simultaneously two overlapping organisations:
• Nikhef is a national research lab (“institute”) funded by the foundation FOM; the
dutch foundation for fundamental research of matter.
• Nikhef is also a collaboration between the Nikhef institute and the particle physics
departements of the UvA (A’dam), the VU (A’dam), the UU (Utrecht) and the
RU (Nijmegen) contribute. In this collaboration all dutch activities in particle
physics are coordinated.
In addition there is a collaboration between Nikhef and the Rijksuniversiteit Groningen
(the former FOM nuclear physics institute KVI) and there are contacts with the Uni-
versities of Twente, Leiden and Eindhoven. For more information see the Nikhef web
page: http://www.nikhef.nl.
The research at Nikhef includes both accelerator based particle physics and astro-particle
physics. A strategic plan, describing the research programmes at Nikhef can be found
on the web, from: www.nikhef.nl/fileadmin/Doc/Docs & pdf/StrategicPlan.pdf .
The accelerator physics research of Nikhef is currently focusing on the LHC experiments:
Alice (“Quark gluon plasma”), Atlas (“Higgs”) and LHCb (“CP violation”). Each of
these experiments search answers for open issues in particle physics (the state of matter
at high temperature, the origin of mass, the mechanism behind missing antimatter) and
hope to discover new phenomena (eg supersymmetry, extra dimensions).
The LHC has started taking data in 2009 and the first LHC physics run has officially
ended in winter 2013. In a joint meeting in summer 2012 the ATLAS and CMS experi-
ments presented their biggest discovery so far, a new particle, called the Higgs particle,
with a mass of approximately 125 GeV. The existence of this particle was a prediction
of the Standard Model. Further research concentrates on signs of physics beyond the
Standard Model, sometimes called “New Physics”.
In preparation of these LHC experiments Nikhef has also been active at other labs:
STAR (Brookhaven), D0 (Fermilab) and Babar (SLAC). Previous experiments that
ended their activities are: L3 and Delphi at LEP, and Zeus, Hermes and HERA-B at
Desy.
A more recent development is the research field of astroparticle physics. It includes
Antares & KM3NeT (“cosmic neutrino sources”), Pierre Auger (“high energy cosmic
rays”), Virgo & ET (“gravitational waves”) and Xenon (”dark matter”).
Nikhef houses a theory departement with research on quantum field theory and gravity,
string theory, QCD (perturbative and lattice) and B-physics.
iv PRELIMINARIES
Driven by the massive computing challenge of the LHC, Nikhef also has a scientific
computing departement: the Physics Data Processing group. They are active in the
development of a worldwide computing network to analyze the large datastreams from
the (LHC-) experiments (“The Grid”).
The book of Griffiths starts with a nice historical overview of particle physics in the
previous century. This is a summary of key events:
Atomic Models
1897 Thomson: Discovery of Electron. The atom contains electrons as “plums in
a pudding”.
1911 Rutherford: The atom mainly consists of empty space with a hard and heavy,
positively charged nucleus.
1913 Bohr: First quantum model of the atom in which electrons circled in stable
orbits, quatized as: L = ~ · n
1932 Chadwick: Discovery of the neutron. The atomic nucleus contains both
protons and neutrons. The role of the neutrons is associated with the binding
force between the positively charged protons.
The Photon
1900 Planck: Description blackbody spectrum with quantized radiation. No inter-
pretation.
1905 Einstein: Realization that electromagnetic radiation itself is fundamentally
quantized, explaining the photoelectric effect. His theory received scepticism.
1916 Millikan: Measurement of the photo electric effect agrees with Einstein’s
theory.
1923 Compton: Scattering of photons on particles confirmed corpuscular character
of light: the Compton wavelength.
Mesons
1934 Yukawa: Nuclear binding potential described with the exchange of a quan-
tized field: the pi-meson or pion.
1937 Anderson & Neddermeyer: Search for the pion in cosmic rays but he finds a
weakly interacting particle: the muon. (Rabi: “Who ordered that?”)
1947 Powell: Finds both the pion and the muon in an analysis of cosmic radiation
with photo emulsions.
Anti matter
1927 Dirac interprets negative energy solutions of Klein Gordon equation as energy
levels of holes in an infinite electron sea: “positron”.
1931 Anderson observes the positron.
vi PRELIMINARIES
1940-1950 Feynman and Stückelberg interpret negative energy solutions as the positive
energy of the anti-particle: QED.
Neutrino’s
1930 Pauli and Fermi propose neutrino’s to be produced in β-decay (mν = 0).
1958 Cowan and Reines observe inverse beta decay.
1962 Lederman and Schwarz showed that νe 6= νµ . Conservation of lepton number.
Strangeness
1947 Rochester and Butler observe V 0 events: K 0 meson.
1950 Anderson observes V 0 events: Λ baryon.
The Eightfold Way
1961 Gell-Mann makes particle multiplets and predicts the Ω − .
1964 Ω − particle found.
The Quark Model
1964 Gell-Mann and Zweig postulate the existence of quarks
1968 Discovery of quarks in electron-proton collisions (SLAC).
1974 Discovery charm quark (J/ψ) in SLAC & Brookhaven.
1977 Discovery bottom quarks (Υ ) in Fermilab.
1979 Discovery of the gluon in 3-jet events (Desy).
1995 Discovery of top quark (Fermilab).
Broken Symmetry
1956 Lee and Yang postulate parity violation in weak interaction.
1957 Wu et. al. observe parity violation in beta decay.
1964 Christenson, Cronin, Fitch & Turlay observe CP violation in neutral K meson
decays.
The Standard Model
1978 Glashow, Weinberg, Salam formulate Standard Model for electroweak inter-
actions
1983 W-boson has been found at CERN.
1984 Z-boson has been found at CERN.
1989-2000 LEP collider has verified Standard Model to high precision.
Introduction
1
2 PRELIMINARIES
~
rBohr = . (i.3)
αme c
(A proper treatment in QM tells you that the expectation value for the radius is not
exactly the Bohr radius, but it comes close.) Hence, the velocity of the electron is
which indeed makes the electron in the hydrogen atom notably non-relativistic.
The second distance scale is the Compton wavelength of the electron. Suppose that
you study electrons by shooting photons at zero-velocity electrons. The smaller the
wavelength of the photon, the more precise you look. However, at some point the
energy of the photons becomes large enough that you can create a new electron. (In
our real theory, you can only create pairs, but that factor 2 is not important now.) The
energy at which this happens is when ~ω = me c2 , or at a wavelength
2π~
λe = . (i.5)
me c
Usually, we divide both sides by 2π and speak of the reduced Compton wavelength λ̄e ,
just like ~ is usually called the reduced Planck’s constant. Note that λ̄e = αrBohr . In
electromagnetic collisions at this energy, classical quantum mechanics no longer suffices:
as soon as collisions involve the creation of new particles, one needs QFT.
Finally, consider the collisions of two electrons at even higher energy. If the electrons
get close enough, the Coulomb energy is sufficient to create a new electron. (Again,
ignore the factor two required for pair production.) Expressing the Coulomb potential
as V (r) = α~c/r, and setting this equal to me c2 , one obtains for the distance
α~
re = . (i.6)
me c
Note that, taking into account the definition of α, this expression does not explicitly
depend on ~: you do not need quantization to compute this distance, which is why it
is usually called the classical radius of the electron. At energies this high lowest order
perturbation theory may not be sufficient to compute a cross-section. The effect of
‘screening’ becomes important, amplitudes described by Feynman diagrams with loops
contribute and QED needs renormalization to provide meaningful answers.
Fortunately for most of us, we will not discuss renormalization in this course. In fact,
we will hardly discuss quantum field theory at all! Do not be disappointed: there are
two pragmatic reasons for this. First, a proper treatment requires a proper QFT course
with some non-trivial math. This would leave insufficient time for other things that
we do need to address. Second, if you accept a little handwaving here and there, then
I.1. QUANTUM FIELDS 3
Table i.1: Values for the Bohr radius, the reduced Compton wavelength of the electron and
the classical radius of the electron, and the corresponding energy.
we do not actually need QFT: starting from quantum mechanics and special relativity
we can derive the ’Born level’ — that is, ‘leading order’ — cross-sections, following a
route that allows us to introduce new concepts in a somewhat historical, and hopefully
enlightening, order.
However, before continuing and setting aside the field theory completely until chapter 8,
it is worthwhile to briefly discuss some relevant features of QFT, in particular those that
distinguish it from ordinary quantum mechanics. In QM particles are represented by
waves, or wave packets. Quantization happens through the ‘fundamental postulate’ of
quantum mechanics that says that the operators for space coordinates and momentum
coordinates do not commute,
[x̂, p̂] = i~ (i.7)
The dynamics of the waves is described by the Schrödinger equation. Scattering cross-
sections are derived by solving, in perturbation theory, a Schrödinger equation with a
Hamiltonian operator that includes terms for kinetic and potential energy. Usually we
expand the solution around the solution for a ‘free’ particle and write the solution as a
sum of plane waves. This is exactly what you have learned in your QM course and we
will come back to this in Lecture 2.
In QFT particles are represented as ‘excitations’ (or ‘quanta’) of a field q(x), a function
of space-time coordinates x. There are only a finite number of fields, one for each type
of particle, and one for each force carier. This solves one imminent problem, namely
why all electrons are exactly identical. In its simplest form QED has only two fields: one
for a spin- 21 electron and one for the photon. The dynamics of these field are encoded
in a Lagrangian density L. Equations of motions are obtained with the principle of
least action. Those for the free fields (in a Lagrangian without interaction terms) lead
to wave equations, reminiscent of the Schrödinger equation, but now Lorentz covariant.
Again, solutions are written as superpositions of plane waves. The fields are quantized
by interpreting the fields as operators and imposing a quantization rule similar to that
in ordinary quantum mechanics, namely
[q, p] = i~ (i.8)
where the momentum p = ∂L/∂ q̇ is the so-called adjoint coordinate to q. (You may
remember that you used similar notation to arrive at Hamilton’s principle in your classi-
cal or quantum mechanics course.) The Fourier components of the quantized fields can
be identified as operators that create or destruct field excitations, exactly what we need
4 PRELIMINARIES
for a theory in which the number of particles is not conserved. The relation to classical
QM can be made by identifying the result of a ‘creation’ operator acting on the vacuum
as the QM wave in the Schrödinger equation.
That was a mouth full and you can forget most of it. One last thing, though: one very
important aspect of quantum field theory is the role of symmetries in the Lagrangian. In
fact, as we shall see in Lecture 8 and 9, the concept of phase invariance allows to define
the standard model Lagrangian by specifying only the matter fields and the symmetries:
once the symmetries are defined, the dynamics (the force carriers) come for free.
That said, we leave the formal theory of quantum fields alone. In the remainder of
this chapter we briefly discuss some concepts and summarize our ‘Standard Model’ of
elementary particles. In Lecture 1 we formulate a relativistic wave equation for a spin-0
particle. In Lecture 2 we discuss classical QM perturbation theory and Fermi’s Golden
rule, which allows us to formalize the computation of a cross-section. In Lecture 3, we
show how the Maxwell equations take a very simple form when expressed in terms of a
new spin-1 field, which we identify as the photon. In lecture 4 we apply the developed
tools to compute the scattering of spin-0 particles. In Lectures 5 and 6 we turn to
spin- 12 fields, which are considerably more realistic given that all SM matter fields are
indeed fermions. In Lectures 7 through 10, we introduce the weak interaction, gauge
theory and electroweak unification. Finally, in Lectures 11-14 we look in more detail at
electroweak symmetry breaking.
Figure i.1: The discovery of the positron as reported by Anderson in 1932. Knowing the
direction of the B field Anderson deduced that the trace was originating from an anti electron.
Question: how?
Furthermore, though the constituents of atoms were fairly well established, there was
something puzzling about atoms: What was keeping the nucleus together? It clearly
had to be a new force, something beyond electromagnetism. Rutherford’s scattering
experiments had given an estimate of the size of the nucleus, of about 1 fm. With
protons packed this close, the new force had to be very strong to overcome the repulsive
coulomb interaction of the protons. (Being imaginative, physicists simply called it the
strong nuclear force.) Yet, to explain scattering experiments, the range of the force had
to be small, bound just to the nucleus itself.
In an attempt to solve this problem Japanese physicist Yukawa published in 1935 a fun-
damentally new view of interactions. His idea was that forces, like the electromagnetic
force and the nuclear force, could be described by the exchange of virtual particles, as
illustrated in Fig. i.3. These particles (or rather, their field ) would follow a relativistic
wave-equation, just like the electromagnetic field.
In this picture, the massless photon was the carrier of the electromagnetic field. As we
will see in exercise 1.4.3 the relativistic wave equation for a massless particle leads to
an electrostatic potential of the form (in natural units, ~ = c = 1)
1
V (r) = −α . (i.9)
r
Because of its 1/r dependence, the force is said to be of ‘infinite range’.
6 PRELIMINARIES
1.0
0.6 0.00008
Mass = 0
0.00004
0.4 Mass = 30 eV
0
18.45 18.50 18.55 18.60
0.2
0
2 6 10 14 18
Energy (keV)
Figure i.2: The beta spectrum as observed in tritium decay to helium. The endpoint of the
spectrum can be used to set a limit of the neutrino mass. Question: how?
e−r/R
U (r) = −g 2 (i.10)
r
which is called the one-pion-exchange-potential. Since it falls of exponentially, it has a
finite range. The range R is inversely proportional to the mass of the force carrier and
for a massless carrier the expression reduces to that for the electrostatic potential.
In the exercise you will derive the relation between mass and range properly, but it
can also be obtained with a heuristic argument, following the Heisenberg uncertainty
principle. (As you can read in Griffiths, whenever a physicists refers to the uncertainty
principle to explain something, take all results with a grain of salt.) In some interpreta-
tion, the principle states that we can borrow the energy ∆E = mc2 to create a virtual
particle from the vacuum, as long as we give it back within a time ∆t ≈ ~/∆E. With
the particle traveling at the speed of light, this leads to a range R = c∆t = c~/mc2 .
From the size of the nucleus, Yukawa estimated the mass of the force carrier to be ap-
proximately 100 MeV/c2 . He called the particle a meson, since its mass was somewhere
in between the mass of the electron and the nucleon.
In 1937 Anderson and Neddermeyer, as well as Street and Stevenson, found that cosmic
rays indeed consist of such a middle weight particle. However, in the years after, it
became clear that this particle could not be Yukawa’s meson, since it did not interact
I.3. FEYNMAN DIAGRAMS 7
Figure i.3: Illustration of the interaction between protons and neutrons by charged pion
exchange. (From Aichison and Hey.)
strongly, which was very strange for a carrier of the strong force. In fact this particle
turned out to be the muon, the heavier brother of the electron.
Only in 1947 Powell and Perkins found Yukawa’s pion in cosmic rays. They took photo-
graphic emulsion plates to mountain tops to study the contents of cosmic rays (see Fig.
i.4). In a cosmic ray event a cosmic proton scatters with high energy on an atmospheric
nucleon and produces many secondary particles. Pions produced in the atmosphere
decay long before they reach sea level, which is why they had not been observed before.
As a carrier of the strong force Yukawa’s meson did not stand the test of time. We now
know that the pion is a composite particle and that the true carrier for the strong force is
the massless gluon. The range of the strong force is small, not because the force carrier
is massive, but because gluons carry a strong interaction charge themselves. However,
even if Yukawa’s original meson model did not survive, his interpretation of forces as
the exchange of virtual particles is still central to the description of particle interactions
in quantum field theory.
Figure i.3 is an example of a Feynman diagram. You have probably seen Feynman
diagrams before and already know that they are not just pictures that help us to ‘visu-
alize’ a scattering process: they can actually be translated efficiently into mathematical
expressions for the computation of quantum mechanical transition amplitudes.
In this course, we will always draw Feynman diagrams such that time runs from left to
right. This is just a convention: the diagrams in Fig. i.3 are equally valid if time runs
from right to left, or from top to bottom, etc.
With this convention, the two diagrams in Fig. i.3 represent two different ways of scat-
8 PRELIMINARIES
Figure i.4: A pion entering from the left decays into a muon and an invisible neutrino.
I.4. THE STANDARD MODEL 9
tering a proton and a neutron via pion exchange: In case (a) a negative virtual pion is
first emitted by the neutron and then absorbed by the proton, while in the case (b) a
positive virtual pion is first emitted by the proton and then absorbed by the neutron.
As usual in quantum mechanics these two contributions are represented by two complex
amplitudes, which need to be added in order to obtain the total amplitude. It turns
out that only if both amplitudes are taken into account, Lorentz-covariant results can
be obtained in a quantum theory.
However, now that we know that both amplitudes must be taken into account, we no
longer need to draw both of them! In fact, in the remainder of this course, we will
always draw only one diagram, with the line that represents the pion exchange drawn
vertically. By convention, Feynman diagrams always present all possible time orderings
for the ‘internal’ lines, the virtual particles.
charge Quarks
2 u (up) c (charm) t (top)
3 1.5–4 MeV 1.15–1.35 GeV (174.3 ± 5.1) GeV
d (down) s (strange) b (bottom)
− 13
4–8 MeV 80–130 MeV 4.1–4.4 GeV
charge Leptons
0 νe (e neutrino) νµ (µ neutrino) ντ (τ neutrino)
< 3 eV < 0.19 MeV < 18.2 MeV
−1 e (electron) µ (muon) τ (tau)
0.511 MeV 106 MeV 1.78 GeV
Table i.2: Matter particles in the Standard Model, with their approximate mass.
Apart from the charged weak interaction (mediated by the W ), all interactions conserve
quark and lepton flavour: They do not change one type of fundamental fermion into
another type. On the other hand, the charged weak interaction allows for transitions
between an up-type quark and a down-type quark, and between charged leptons and
neutrinos. Some of the fundamental diagrams are shown in figure i.5.
The fundamental particles carry a charge (essentially a quantum number) to characterize
10 PRELIMINARIES
Table i.3: Standard Model forces, the mediating bosons, and the associated strength of the
coupling at an energy of about 1 GeV. (The latter are taken from Thomson, 2013.)
e+ µ− νe µ− q q
γ W g
a: b: c:
e− µ+ e− νµ q q
Figure i.5: Feynman diagrams of fundamental lowest order perturbation theory processes in
a: electromagnetic, b: weak and c: strong interaction.
how they couple to the force carriers. For the electromagnetic interaction this charge
is the electric charge, for the weak interaction it is called weak isospin and for the
strong interaction it is called colour. The strength of an interaction is determined by
the charges of the particles involved, by a coupling strength that is a characteristic of
the interaction (see table i.3), and by the mass of the vector boson.
There is an important difference between the electromagnetic force on one hand, and
the weak and strong force on the other hand. The photon does not carry electric charge
and, therefore, does not interact with itself. The gluons, however, carry colour and do
interact amongst each other. Also, the weak vector bosons carry weak isospin charge
and undergo this so-called self-coupling.
The coupling strengths are not actually constant, but vary as a function of energy, which
is called the running of the coupling constants. At a momentum transfer of 1015 GeV
the couplings of electromagnetic, weak and strong interaction all obtain approximately
the same value. (See figure i.6.) Grand unification refers to the hypothesis that at high
energy there is actually only a single force, originating from a single gauge symmetry
with a single coupling constant.
Due to the self-coupling of the force carriers the running of the coupling constants of the
weak and strong interaction are opposite to that of electromagnetism. Electromagnetism
becomes weaker at low momentum (i.e. at large distance), the weak and the strong force
become stronger at low momentum or large distance. The strong interaction coupling
becomes so large at momenta less than a few 100 MeV that perturbation theory is
no longer applicable. (The coupling constant is larger than 1.) Although this is not
rigorously proven, it is assumed that the self-coupling of the gluons is also responsible
for confinement: the existence of free coloured objects (i.e. objects with net strong
charge) is forbidden.
Confinement means that free quarks do not exist, at least, not at time-scales longer
I.4. THE STANDARD MODEL 11
Figure i.6: Running of the coupling constants and possible unification point. On the left:
Standard Model. On the right: Supersymmetric Standard Model.
than that corresponding to the range of the strong interaction. Quarks always appear
in bound states, either as combinations of three quarks (baryons) or as combinations
of a quark an an anti-quark (mesons). Together these are called hadrons. In the quark
model the various species of hadrons are organized by exploiting quark flavour symmetry,
the fact that equally charged quarks of different families are indistinguisable except for
their mass. Due to lack of time, we will not discuss the quark model in this course.
For reference table i.4 gives a list of common hadrons, some of which we encounter in
examples in the lectures.
Table i.4: Name, quark content and approximate mass of common baryons and mesons. The
complete list of all known hadrons, together with a lot of experimental data, can be found in
the particle data book, http://pdglive.lbl.gov.
Finally, the Standard Model includes a scalar boson field, the Higgs field, which provides
mass to the vector bosons and fermions in the Brout-Englert-Higgs mechanism. The
motivation for the Higgs particle and corresponding precision tests of the SM are the
subject of the last four lectures of this course.
12 PRELIMINARIES
Despite the success of the standard model in describing all physics at ’low’ energy scale,
there are still many open questions, such as:
• why are the masses of the particles what they are?
• why are there 3 generations of fermions?
• are quarks and leptons truly fundamental?
• is there really only one Higgs particle?
• why is the charge of the electron exactly opposite to that of the proton? Or phrased
differently: why is the total charge of leptons and quarks in one generation exactly
zero?
• is a neutrino its own anti-particle?
• can all forces be described by a single gauge symmetry (unification)?
• why is there no anti-matter in the universe?
• what is the source of dark matter?
• what is the source of dark energy?
Particle physicists try to address these questions with scattering experiments in the
laboratory and by studying high energy phenomena in the cosmos.
such that all factors ~ and c can be omitted. As a consequence (see textbooks), there is
only one basic unit for length (L), time (T), mass (M), energy and momentum. In high
energy physics this basic unit is often chosen to be the energy in MeV or GeV, where
1 eV is the kinetic energy an electron obtains when it is accelerated over a voltage of
1 V. Momentum and mass then get units of energy, while length and time get units of
inverse energy.
To confront the result of calculation with experiments the factors ~ and c usually need to
be reintroduced. There are two ways to do this. First one can take the final expressions
in natural units and then use the table i.5 to convert the quantities for space, time, mass,
I.5. UNITS IN PARTICLE PHYSICS 13
energy and momentum back to their original counterparts. (For the positron charge,
see below.) Alternatively, one can express all results in GeV, then use the table i.6 with
conversion factors to translate it into SI units.
Table i.5: Conversion of basic quantities between natural and ordinary units.
e2 1
F = , (i.13)
4π0 r2
where 0 is the vacuum permittivity. The dimension of the factor e2 /0 is fixed — it is
[L3 M/T2 ] — but this still leaves a choice of what to put in the charges and what in the
vacuum permittivity.
In the SI system the unit of charge is the Coulomb. (It is currently defined via the
Ampére, which in turn is defined as the current leading to a particular force between
two current-carrying wires. In the near future, this definition will probably be replaced
by the charge corresponding to a fixed number of particles with the positron charge.)
The positron charge expressed in Coulombs is about
As we shall see in Lecture 3 the Maxwell equations look much more neat if, in addition
to c = 1, we choose 0 = 1. This is called the Heaviside-Lorentz system. Obviously, this
choice affects the numerical value of e. However, note that coupling constant α, defined
14 PRELIMINARIES
in equation i.2, is dimensionless and hence independent of the system of units. In this
course we will often write e2 , when in fact we mean α.
Finally, it is customary to express scattering cross-sections in barn: one barn is equal
to 10−24 cm2 .
xµ = (x0 , x1 , x2 , x3 ) (i.16)
where the first component x0 = ct is the time coordinate and the latter three components
are the spatial coordinates (x1 , x2 , x3 ) = x. Under a Lorentz transformation along the
x1 axis with velocity β = v/c, xµ transforms as
0
x0 = γ(x0 − βx1 )
0
x1 = γ(x1 − βx0 )
0 (i.17)
x2 = x2
0
x3 = x3
p
where γ = 1/ 1 − β 2 .
A ‘contravariant four-vector’ is defined to be any set of four quantities Aµ = (A0 , A1 , A2 , A3 ) =
(A0 , A) that transforms under Lorentz transformations exactly as the corresponding
components of the coordinate four-vector xµ . Note that it is the transformation prop-
erty that defines what a ‘contravariant’ vector is.
Lorentz transformations leave the quantity
2
|A|2 = A0 − |A|2 (i.18)
invariant. This expression may be regarded as the scalar product of Aµ with a related
‘covariant vector’ Aµ = (A0 , −A), such that
X
A · A ≡ |A|2 = Aµ Aµ . (i.19)
µ
From now on we omit the summation sign and implicitly sum over any index that
appears twice. Defining the metric tensor
1 0 0 0
0 −1 0 0
gµν = g µν =
0 0 −1 0
(i.20)
0 0 0 −1
I.6. FOUR-VECTOR NOTATION 15
Note that the position of the minus sign is opposite to that of the coordinate four-vector
itself.
16 PRELIMINARIES
Lecture 1
1 ∂2
2
− ∇ φ(x, t) = 0 . (1.1)
c2 ∂t2
The solution to this equation is given by plane waves of the form
where the wave-vector k and the angular frequency ω are related by the dispersion
relation
ω = c|k| . (1.3)
Since the equation above is real, we can restrict ourselves to real solutions. In fact,
the photon field is described by just the real component of the wave. However, it is
often more convenient to work with complex waves. The real solutions can always be
identified with one component of the complex wave.
Maxwell identified propagating electromagnetic fields with light, and thereby firmly
established what everybody already knew: light behaves as a wave. However, to explain
the photo-electric effect Einstein hypothesized in 1904 that light is also a particle with
zero mass. For a given frequency, lights comes in packets (‘quanta’) with a fixed energy.
The energy of a quantum is related to the frequency by
E = hν = ~ω , (1.4)
p = ~k . (1.5)
17
18 LECTURE 1. WAVE EQUATIONS AND ANTI-PARTICLES
In terms of energy and momentum the dispersion relation takes the familiar form E = pc.
The idea of light as a particle was received with much skepticism and only generally
accepted after Compton showed in 1923 that photons scattering of electrons behave as
one would expect from colliding particles.
So, by 1923 light was a wave and a particle: it satisfied a wave equation, yet it only came
about in packets of discrete energy. That lead De Broglie in 1924 to make another bold
preposition: if light is both a wave and a particle, then why wouldn’t matter particles
be waves as well? It took another few years before physicists established the wave-like
character of electrons in diffraction experiments, but well before that people took De
Broglie hypothesis seriously and started looking for a suitable wave-equation for massive
particles.
The crucial element is to establish the dispersion relation for the wave. Schrödinger
started with the relativistic equation for the total energy
E 2 = m2 c4 + p2 c2 , (1.6)
but abandoned the idea, for reasons we will discuss later. He then continued with the
equation for the kinetic energy in the non-relativistic limit
p2
E = , (1.7)
2m
which, as we shall see now, led to his famous equation.
Since total probability is conserved, the density must satisfy a so-called continuity equa-
tion
∂ρ
+∇·j =0 (1.12)
∂t
where j is the density current or flux. When considering charged particles you can
think of ρ as the charge per volume and j as the charge times velocity per volume. The
continuity equation can then be stated in words as “The change of charge in a given
volume equals the current through the surrounding surface”.
What is the current corresponding to a quantum mechanical wave ψ? It is straightfor-
ward to obtain this current from the continuity equation by writing ∂ρ/∂t = ψ∂ψ ∗ /∂t +
ψ ∗ ∂ψ/∂t and inserting the Schrödinger equation. However, because this is useful later
on, we follow a slightly different approach. First, rewrite the Schrödinger equation as
∂ i~ 2
ψ= ∇ ψ.
∂t 2m
Now multiply both sides on the left by ψ ∗ and add the expression to its complex conju-
gate
∗ ∂ψ ∗ i~
ψ = ψ ∇2 ψ
∂t 2m
∂ψ ∗
−i~
ψ = ψ ∇2 ψ ∗
∂t 2m
+
∂ i~
(ψ ∗ ψ) = −∇ · (ψ ∇ψ ∗ − ψ ∗ ∇ψ) (1.13)
∂t | {z } 2m
ρ | {z }
j
where in the last step we have used that ∇ · (ψ ∗ ∇ψ − ψ∇ψ ∗ ) = ψ ∗ ∇2 ψ − ψ∇2 ψ ∗ . In the
result we can recognize the continuity equation if we interpret the density and current
as indicated.
Plane waves of the form
ψ = N ei(p·x−Et)/~ (1.14)
with E = p2 /2m are solutions to the free Schrödinger equation. (In fact, starting
from the idea of particle-wave duality, Schrödinger took the plane wave form above and
‘derived’ his equation as the equation that described its time evolution.) We will leave
the definition of the normalization constant N for the next lecture: as the plane wave
is not localized in space (it has precise momentum, and infinitely imprecise position!),
it can only be normalized on a finite volume.
To get rid of the inconvenient factor ~ in the exponent, we can express energy and
momentum in terms of the wave vector and angular frequency defined above or work
in natural units. Since it is sometimes useful to verify expressions with a dimensional
analysis, we shall keep the factors ~ for now.
20 LECTURE 1. WAVE EQUATIONS AND ANTI-PARTICLES
with E = p2 /2m. For t = 0 this is just the usual Fourier transform. For the exercises,
remember that in one dimension the Fourier transform and its inverse are given by
(Plancherel’s theorem),
Z +∞ Z +∞
1 1
f (x) = √ ikx
F (k)e dk ⇐⇒ F (k) = √ f (x)e−ikx dx (1.18)
2π −∞ 2π −∞
1 ∂2 2 m 2 c2
− φ = −∇ φ + φ (1.20)
c2 ∂t2 ~2
1.3. THE KLEIN-GORDON EQUATION 21
This equation is called the Klein-Gordon equation. Having once seen it with the factors
~ and c included, we will from now on omit them. The Klein-Gordon equation can then
be efficiently written in four-vector notation as
2 + m2 φ(x) = 0 ,
(1.21)
where
1 ∂2
2 ≡ ∂µ ∂ µ ≡ 2 2
− ∇2 (1.22)
c ∂t
is the so-called d’Alembert operator.
Unlike the Schrödinger equation, the KG equation does not contain factors i. Conse-
quently, it can have both real and complex solutions. These have different applications.
In chapter 3 we shall see an example of the KG equation for a real field. In this section
we assume that the waves are complex.
Planes waves of the form
µ
φ(x) = N ei(px−Et) = e−ipµ x (1.23)
µ
with p = (E, p) are solutions of the KG equation provided that they satisfy the dis-
persion relation E 2 = p2 + m2 . Nothing restricts solution to have positive energy: we
discuss the interpretation of negative energy solutions later in this lecture.
Any solution to the KG equation can be written as a superposition of plane waves, like
for the Schrödinger equation. However, in contrast to the classical case, the complex
conjugate of the plane wave above
µ
φ∗ (x) = N ei(−px+Et) = eipµ x (1.24)
is also a solution to the KG equation and need to be accounted for in the decomposition.
It is not independent though, since φ∗ (p, E) = φ(−p, −E). Consequently, we can write
the generic decomposition restricting ourselves to positive energy solutions, if we write
Z
µ µ
d3 p A(p) e−ipµ x + B(p) eipµ x
φ(x) = (1.25)
p
with E = + p2 + m2 . By popular convention, motivated later, we identify the first
exponent as an incoming particle wave, or an outgoing anti-particle wave, and vice-versa
for the second exponent.
In analogy to Eq. 1.13 we now derive a continuity equation. We multiply the Klein Gor-
den equation for φ from the left by −iφ∗ , then add to the complex conjugate equation:
2
∗ ∂ φ
= −iφ∗ −∇2 φ + m2 φ
−iφ − 2
∂t
2 ∗
∂ φ
= iφ −∇2 φ∗ + m2 φ∗
iφ − 2
∂t
+
∗
∂ ∂φ ∂φ
i φ∗ −φ = ∇ · [i (φ∗ ∇φ − φ ∇φ∗ )] (1.26)
∂t ∂t ∂t | {z }
| {z } −j
ρ
22 LECTURE 1. WAVE EQUATIONS AND ANTI-PARTICLES
where we can recognize again the continuity equation. In four-vector notation the con-
served current becomes
∂µ j µ = 0 (1.28)
You may wonder why we introduced the factor i in the current: this is in order to make
the density real.
Substituting the plane wave solution gives
ρ = 2 |N |2 E
(1.29)
j = 2 |N |2 p
or in four-vector notation
j µ = 2 |N |2 pµ . (1.30)
Like for the the classical Schrödinger equation, the ratio of the current to the density
is still a velocity since v = p/E. However, in contrast to the non-relativistic case,
the density of the Klein-Gordon wave is proportional to the energy. This is a direct
consequence of the Klein-Gordon equation being second order in the time derivative.
We write the conserved current as a four-vector assuming that it transforms under
Lorentz transformation the way four-vector are supposed to do. It is not so hard to show
this by looking at how a volume and velocity change under Lorentz transformations (see
e.g. the discussion in Feynman’s Lectures, Vol. 2, sec. 13.7.) The short argument
is that since φ is a Lorentz-scalar, and ∂ µ a Lorentz vector, their product must be a
Lorentz vector.
You may remember that conservation rules in physics are related to symmetries. That
makes you wonder which symmetry leads to the conserved currents for the Schrödinger
and Klein-Gordon equations. In Lecture 8 we discuss Noether’s theorem and show that
it is the phase invariance of the Lagrangian, a so-called U (1) symmetry. The phase of
the wave functions is not a physical observable. For QM wave functions the conserved
current implies that probability is conserved. For the QED Lagrangian it implies that
charge is conserved.
E E
+m +m
−m −m
In a feat that is illustrative for his ingenuity Dirac turned to Pauli’s exclusion principle.
The exclusion principle states that identical fermions cannot occupy the same quantum
state. Dirac’s picture of the vacuum and of a particle are schematically represented in
Fig. 1.1.
The plot shows all the available energy levels of an electron. Its lowest absolute energy
level is given by |E| = m. Dirac imagined the vacuum to contain an infinite number
of states with negative energy which are all occupied. Since an electron is a spin-1/2
particle each state can only contain one spin “up” electron and one spin-”down” electron.
All the negative energy levels are filled. Such a vacuum (“sea”) is not detectable since
the electrons in it cannot interact, i.e. go to another state.
If energy is added to the system, an electron can be kicked out of the sea. It now gets
a positive energy with E > m. This means this electron becomes visible as it can now
interact. At the same time a “hole” in the sea has appeared. This hole can be interpreted
as a positive charge at that position: an anti-electron! Dirac’s original hope was that he
could describe the proton in such a way, but it is essential that the anti-particle mass
is identical to that of the electron. Thus, Dirac predicted the positron, a particle that
can be created by ’pair production’. The positron was discovered in 1931 by Anderson.
24 LECTURE 1. WAVE EQUATIONS AND ANTI-PARTICLES
There is one problem with the Dirac interpretation: it only works for fermions!
Pauli and Weiskopf proposed in 1934 that the density should be regarded as a charge
density. For an electron the charge density is written as
j µ = −ie(φ∗ ∂ µ φ − φ∂ µ φ∗ ). (1.31)
Stückelberg and later Feynman took this approach one step further. Consider the current
for a plane wave describing an electron with momentum p and energy E. Since the
electron has charge −e, this current is
Now consider the current for a positron with momentum p. Its current is
Consequently, the current for the positron is identical to the current for the electron
but with negative energy and traveling in the opposite direction. Or, in terms of the
plane waves, to go from the positron current to the electron current, we just need to
µ
change the sign in the exponent of eixµ p . By our earlier convention, this is equivalent
to saying that the incoming plane wave of a positron is identical to the outgoing wave
of an electron.
Now consider what happens to the electron wave function if we change the direction of
time: We will have ct → −ct and p → −p. You immediately notice that this has exactly
the same effect on the plane wave exponent as the transformation (E, p) → (−E, −p).
In other words, we can interprete the negative energy current of the electron as an
electron moving backward in time. This current is identical to that of a positron moving
forward in time.
This interpretation, illustrated in Fig. 1.2, is very convenient when computing scattering
amplitudes: in our calculations with Feynman diagrams we can now express everything
in terms of particle waves, replacing every anti-particle with momentum pµ by a particle
with momentum −pµ , as if it were traveling backward in time. For example, the process
of an absorption of a positron with energy E is the same as the emission of an electron
1.4. INTERPRETATION OF NEGATIVE ENERGY SOLUTIONS 25
e+ e−
E>0 E<0
Figure 1.2: A positron travelling forward in time is an electron travelling backwards in time.
with energy −E (see Fig.1.3). Likewise, the process of an incoming positron scattering
off a potential will be calculated as that of a scattering electron travelling back in time
(see Fig. 1.4).
−e
(+E,p)
absorption
time
emission
+e
(−E,−p)
Figure 1.3: There is no difference between the process of an absorption of a positron with
pµ = (−E, −p) and the emission of an electron with pµ = (e, p).
e+ −
e
time x
µ µ
Figure 1.4: In terms of the charge current density j+(E,p) (+e) ≡ j−(E,p) (−e)
The advantage of this approach becomes more apparent when one considers higher order
corrections to the amplitudes. Consider the scattering of an electron on a localized
potential, illustrated in Fig. 1.5. To first order the interaction of the electron with the
perturbation is described by the exchange of a single photon. When the calculation is
extended to second order the electron interacts twice with the field. It is important to
note that this second order contribution can occur in two time orderings as indicated in
26 LECTURE 1. WAVE EQUATIONS AND ANTI-PARTICLES
the figure. These two contributions are different and both of them must be included in
a relativistically covariant computation.
The time-ordering on the right can be viewed in two ways:
• The electron scatters at time t2 runs back in time and scatters at t1 .
• First at time t1 “spontaneously” an e− e+ pair is created from the vacuum. Later-
on, at time t2 , the produced positron annihilates with the incoming electron, while
the produced electron emerges from the scattering process.
The second interpretation would allow the process to be computed in terms of particles
and anti-particles that travel forward in time. However, the first interpretation is just
more economic. We realize that the vacuum has become a complex environment since
particle pairs can spontaneously emerge from it and dissolve into it!
e−
time
e− e−
t2 x t2 x
t1 t1 x
x
e − e−
Exercises
Exercise 1.1 (Conversion factors)
Derive the conversion factors for mass, length and time in table i.6.
∂2
2V =0 ; 2 ≡ ∂ µ∂ µ ≡ − ∇2
∂t2
which in the static case can be written in the form of Laplace equation:
∇2 V = 0
Now consider a point charge in vacuum. Exploiting spherical symmetry, show that
this equation leads to a ‘potential’ V (r) ∝ 1/r.
Hint: look up the expression for the Laplace operator in spherical coordinates.
(b) The wave equation for a massive field is the Klein Gordon equation:
2 U + m2 U = 0
∇2 U − m2 U = 0
(c) Estimate the mass of the π-meson assuming that the range of the nucleon force is
1.5 × 10−15 m = 1.5 fm.
π − + p → K 0 + anything
Exercise 1.6 (From A&H, chapter 3. See also Griffiths, exercise 7.1)
In this exercise we derive expression Eq. (i.22) of the Introductory chapter.
(a) Start with the expressions for a Lorentz transformation along the x1 axis in
Eq. (i.17). Write down the inverse transformation ( i.e. express (x0 , x1 ) in
0 0
(x0 , x1 ))
0 0
(b) Use the chain rule to express the derivatives ∂/∂x0 and ∂/∂x1 in the derivatives
∂/∂x0 and ∂/∂x1 .
(c) Use the result to show that (∂/∂x0 , −∂/∂x1 ) transforms in the same way as
(x0 , x1 ).
with a real and positive. Compute the normalization constant A such that
Z +∞
|ψ(x, 0)|2 dx = 1.
−∞
Hint: ∞ √
Z
2
e−y dy = π
∞
(b) Take the Fourier transform to derive the wave function in momentum space at
t = 0,
1/4
1 2
ψ(k) = e−(k−k0 ) /4a
2aπ
and then move the integration boundaries by −b/2a. (Don’t mind that b is com-
plex.)
(c) Use this result and Eq. (1.17) to show that the solution to the Schrödinger equation
(with E(p) = p2 /2m or ω(k) = ~k 2 /2m) is given by
1/4
−ax2 + ik0 x − iνtk02 /4a
2a −1/2
ψ(x, t) = (1 + iνt) exp
π 1 + iνt
with ν ≡ 2~a/m.
(d) Compute |ψ(x, t)|2 . Qualitatively, what happens to ψ 2 as time goes on?
(e) Now compute the same for a solution to the massless Klein-Gordon equation (ω =
ck). Note that the wave packet maintains its size as a function of time.
(a) Quarks are fermions with spin 1/2. Show that the spin of a meson (2 quarks) can
be either a triplet of spin 1 or a singlet of spin 0.
Hint: Remember the Clebsch Gordon coefficients in adding quantum numbers.
In group theory this is often represented as the product of two doublets leads to
the sum of a triplet and a singlet: 2 ⊗ 2 = 3 ⊕ 1 or, in terms of quantum numbers:
1/2 ⊗ 1/2 = 1 ⊕ 0.
(b) Show that for baryon spin states we can write: 1/2 ⊗ 1/2 ⊗ 1/2 = 3/2 ⊕ 1/2 ⊕ 1/2
or equivalently 2 ⊗ 2 ⊗ 2 = 4 ⊕ 2 ⊕ 2
30 LECTURE 1. WAVE EQUATIONS AND ANTI-PARTICLES
(c) Let us restrict ourselves to two quark flavours: u and d. We introduce a new
quantum number, called isospin in complete analogy with spin, and we refer to
the u quark as the isospin +1/2 component and the d quark to the isospin -1/2
component (or u= isospin “up” and d=isospin “down”). What are the possible
isospin values for the resulting baryon?
(d) The ∆++ particle is in the lowest angular momentum state (L = 0) and has
spin J3 = 3/2 and isospin I3 = 3/2. The overall wavefunction (L⇒space-part,
S⇒spin-part, I⇒isospin-part) must be anti-symmetric under exchange of any of
the quarks. The symmetry of the space, spin and isospin part has a consequence
for the required symmetry of the Colour part of the wave function. Write down
the colour part of the wave-function taking into account that the particle is colour
neutral.
(e) In the case that we include the s quark the flavour part of the wave function
becomes: 3 ⊗ 3 ⊗ 3 = 10 ⊕ 8 ⊕ 8 ⊕ 1. In the case that we include all 6 quarks it
becomes: 6 ⊗ 6 ⊗ 6. However, this is not a good symmetry. Why not?
Lecture 2
In this chapter we discuss Fermi’s golden rule, which allows us to compute cross-sections
and decay rates. A very readable account of this is given in Griffiths chapter 6 and
Thomson section 2.3.6 and chapter 3.
Most species of particles do not live long. This holds for all baryons except the proton
(even the neutron decays, when it is not inside a nucleus), but also for the muon and
the tau. As particles do not age, the probability to decay is independent of time. Given
a large number of particles N0 , the number of surviving particles is hence given by the
exponential law
N (t) = N0 e−t/τ , (2.1)
where τ is the mean lifetime. For particles that decay via the weak interaction, the mean
lifetime is typically 10−12 − 10−9 seconds. A notable exception is the neutron which lives
for about 15 minutes.
The mean lifetime is inversely proportional to the decay width
~
Γ = , (2.2)
τ
which has units of energy. If the particle can decay through different decay channels
(e.g. a charged pion can decay to µ− ν̄µ and to e− ν̄µ ), the decay width can be written
as the sum of the decay widths to the individual channels
X
Γ = Γi . (2.3)
i
31
32 LECTURE 2. PERTURBATION THEORY AND FERMI’S GOLDEN RULE
The ratio Γi /Γ is called the branching fraction. The particle data book is full of branching
fractions of species in the particle zoo. The decay rate to a particular final state i, i.e.
the average number of decays per unit of time into i, is just Γi /~.
If the decay is to more than two particles, the distribution of angles and energies of
particles in the final state becomes an observable as well. That is why we often consider
partial or differential decay widths,
dΓ
, (2.4)
dp1 · · · dpN
where p1 , . . . , pN are the momenta of the N particles in the final state.
Besides decay widths we also measure scattering cross-sections. (In fact, in our com-
putations, decays and scattering are quite similar, so we deal with both at once.) In
scattering experiments we collide beams of particles and study the collision rate. Con-
sider an experiment in which we scatter a beam of particles A on a target of particles
B. If nA is the particle number density in the beam, and vA is the particle velocity, the
number of collisions per second per unit volume of B is
dN
= vA nA nB σtot . (2.5)
dt
The quantity σtot is called the total scattering cross-section. It has units of ‘surface’.
In most cases we do not study the total collision rate, but rather the rate of particular
final states. The total cross-section is a sum of cross-sections for all possible final states,
such that X
σtot = σi . (2.6)
i
Since the energy and direction of final state particles can be measured as well, we usually
consider differential scattering cross-sections,
dσ(A + B → f1 + · · · + fN )
. (2.7)
dp1 · · · dpN
The expression for the calculation of a (differential) cross-section can be written schemat-
ically as
Wfi
dσ = dΦ (2.8)
flux
The ingredients to this expression are:
1. the transition rate Wfi . You can think of this as the probability per unit time and
unit volume to go from an initial state i to a final state f ;
2. a flux factor that accounts for the ‘density’ of the incoming states;
3. the Lorentz invariant phase space factor dΦ, sometimes referred to as ‘dLIPS’.
It accounts for the density of the outgoing states. (It takes care of the fact that
experiments cannot observe individual states but integrate over a number of states
with nearly equal momenta.)
2.2. NON-RELATIVISTIC SCATTERING 33
The ‘physics’ (the dynamics of the interaction) is contained in the transition rate Wfi .
The flux and the phase space factors are just ‘bookkeeping’, required to compare the
result with the measurements.
The rigorous computation of the transition rate requires quantum field theory, which
is outside the scope of this course. However, to illustrate the concepts we discuss non-
relativistic scattering of a single particle in a time-dependent potential and formulate
the result in a Lorentz covariant way. In the next chapter we will derive the lowest order
amplitude for the scattering of A + B → A + B, which can still be done without field
theory. We can link that result to the ‘Feynman rules’ derived in field theory.
H0
t=−T/2 t=T/2
t=0 ψ
f
H0
ψi V(x,t)
Figure 2.1: Scattering of a single particle in a potential.
Consider the scattering of a particle in a potential as depicted in Fig. 2.1 Assume that
both long before and long after the interaction takes place, the system is described by
the free Schrödinger equation,
∂ψ
i~ = H0 ψ (2.9)
∂t
where H0 is the unperturbed, time-independent Hamiltonian for a free particle. Let
φm (x) be a normalized eigenstate of H0 with eigenvalue Em ,
is a solution to the Schrödinger equation. Since these states form a complete set, any
other wave function can be written as a superposition of the wave functions ψm .
Now consider a Hamiltonian that includes a time-dependent perturbation,
∂ψ
i~ = (H0 + V (x, t)) ψ . (2.13)
∂t
Any solution ψ can be written as
∞
X
ψ= an (t) φn (x) e−iEn t . (2.14)
n=0
where we have used that the ψm are solutions of the free Schrödinger equation. Multiply
the resulting equation from the left with ψf∗ = φ∗f (x) eiEf t and integrate over x to obtain
∞ Z
dan (t)
d3 x φ∗f (x) φn (x) e−i(En −Ef )t~ =
X
i~
n=0
dt
| {z }
δf n
∞ Z
d3 x φ∗f (x) V (x, t) φn (x) e−i(En −Ef )t/~ (2.16)
X
an (t)
n=0
Using the orthonormality relation for φm we then arrive at the following coupled linear
differential equation for ak (t),
∞
dak (t) X
i~ = an (t) Vkn eiωkn t , (2.17)
dt n=0
In some cases the set of equations (2.17) can be solved explicitly. A general solution is
obtained in perturbation theory, by expanding in Vkn . The approximation of order p + 1
can be obtained by inserting the p-th order result on the right hand side of Eq. (2.17),
(p+1)
dak (t) X
i~ ≈ a(p)
n (t)Vkn (t)e
iωkn t
(2.20)
dt n
Without loss of generality we now assume that the incoming wave is prepared in eigen-
state i of the free Hamiltonian, i.e. ak (−∞) = δki . The zeroeth order approximation
(0)
then is ak (t) = δki (no interaction occurs) and the first order result becomes
(1)
da (t)
i~ k = Vki (t)eiωki t (2.21)
dt
(1)
Using that af (−∞) = 0 and integrating this equation we obtain for the coefficient ak (t)
at time t,
Z t Z t
(1) daf (t0 ) 0 1 0
ak (t) = dt = Vki (t0 )eiωki t dt0 for k 6= i (2.22)
−∞ dt i~ −∞
Higher order approximations can be obtained by inserting the lowest order solution in
the right side of Eq. (2.20). (See textbooks.) A graphical illustration of the first and
second order perturbation is given in Fig. 2.2. The lowest order approximation makes
one ‘quantum step’ from the initial state i to the final state f , while the second order
approximation includes all amplitudes i → n → f .
1−st order 2−nd order
f f
time
Vfn
space Vfi Vni
i i
Figure 2.2: First and second order approximation in scattering.
In the following we only consider the first order approximation (Born approximation).
We define the transition amplitude Tfi as the amplitude to go from a state i to a final
state f at large times,
Z ∞ Z
1
Tfi ≡ af (t → ∞) = dt d3 x ψf∗ (x, t) V (x, t) ψi (x, t) (2.23)
i~ −∞
where we substituted the definitions of Vkn and ωkn . We can write the result more
compactly as
Z
1
Tfi = d4 x ψf∗ (x) V (x) ψi (x) (2.24)
i~
36 LECTURE 2. PERTURBATION THEORY AND FERMI’S GOLDEN RULE
Somewhat deceptively, the expression for Tfi seems to have a Lorentz covariant form.
However, as we have seen in the previous lecture, the ‘classical’ free particle waves cor-
responds to a density that does not correctly transform under Lorentz transformations.
Therefore, Tfi is actually not yet a proper Lorentz scalar.
We now make a simplification and consider a potential that is time-independent. The
expression for the transition amplitude then becomes
Vfi ∞ iωfi t
Z
Tfi = e dt = −2π i Vfi δ (Ef − Ei ) (2.25)
i~ −∞
where we have used that the integral is an important representation of the Dirac δ
function Z +∞
1
δ(x) = eikx dk (2.26)
2π −∞
and substituted our definition of ωfi . The δ function expresses conservation of energy.
Note that Tfi is dimensionless.
Can we interprete |Tfi |2 as a probability? Well, there is one conceptual problem and one
pragmatic problem. The conceptual problem is that if the potential is time-independent,
then this probability will just grow with time. The pragmatic problem is that there is
the δ function. These issues can be solved by considering a potential that is turned on
for a ‘finite time’ T . We define the mean transition rate in the limit for large T as
|Tfi |2
Wfi ≡ lim . (2.27)
T →∞ T
For an interaction that is turned on at time −T /2 and turned off at time T /2, the
equation above can be integrated to give for the transition amplitude at T /2,
Z T /2
Vfi 0 2Vfi sin(ωfi T /2)
af (T /2) = eiωfi t dt0 = . (2.28)
i~ −T /2 i~ ωfi
The function on the right is strongly peaked near ωfi = (Ef − Ei )/~ = 0, again enforcing
energy conservation. In fact, for T → ∞ it is yet another representation of the Dirac δ
function,
1 sin2 αx
δ(x) = lim . (2.30)
α→∞ π αx2
2π 2π
Wfi = 2
|Vfi |2 δ(ωfi ) = |Vfi |2 δ(Ef − Ei ) . (2.31)
~ ~
2.2. NON-RELATIVISTIC SCATTERING 37
You can verify that Wfi is indeed a rate: Vfi is an energy, one of the factors of energy is
canceled by the δ function and the other one is divided by ~ to turn it into reciprocal
time.
As indicated before we can never actually probe final states with definite energy in a
measurement with finite duration. In general, there will be a number of states with
energy close to Ei that can be reached. Assuming that these states can be numbered by
a continuos variable n, the total transition rate can be written as an integral over these
final states
Z
W fi ≡ Wfi dn
Z (2.32)
2π 2
= |Vfi | δ(En − Ei ) dn .
~
If ρ(Ef ) is the density of states per unit energy near Ef , the number of final states with
energy between Ef and Ef + dEf is given by
Inserting this in the expression above, we obtain Fermi’s (Second) Golden Rule,
Z
W fi ≡ Wfi ρ (Ef ) dEf
(2.34)
2π
= |Vfi |2 ρ (Ei ) .
~
Note that in this expression ρ(Ei ) is really the density of final states at the energy Ei .
Some textbooks therefore write this as ρ(Ef )|Ef =Ei .
Above, we encountered a δ function in the transition amplitude. To deal with the square
of that δ function we considered a finite time interval and went back to the expression
(0)
for ak (T /2) for finite times T , taking the limit T → ∞ only after taking the square. To
make the final step you need to recognize the special representation of the δ function. For
future applications it is useful to know that one can also solve this problem differently,
namely by taking the limit T → ∞ one integral at a time:
! !∗
1 Vfi T /2 iωfi t Vfi T /2 iωfi t0 0
Z Z
|Wfi | = lim e dt e dt
T →∞ T i~ −T /2 i~ −T /2
Z T /2
|Vfi |2 1
= 2πδ(ωfi ) lim dt0
~2 T →∞ T −T /2
| {z }
T
The final result is of course identical. We will encounter this trick at various places
when going from a transition amplitude to a transition rate.
38 LECTURE 2. PERTURBATION THEORY AND FERMI’S GOLDEN RULE
You may wonder why we need to consider a finite time interval T . The reason is that
when we assume that the initial state is an eigenstate of the free Hamiltonian with fixed
momentum (or energy), we have lost track of where a particle is in both space and
time. A moving wave packet would see the static potential during a finite time, but the
plane waves do not. Just like we will need to normalize the wave functions on a finite
volume, we will need to normalize the potential to a finite time. A proper treatment is
rather lengthy and relies on the use of wave packets. (See e.g. the book by K.Gottfried,
“Quantum Mechanics” (1966), Volume 1, sections 12, 56.) In the end, we can write
transition probabilities in terms of plane waves, provided that we normalize to T and
V . We discuss the normalization in more detail below.
C
e−f
Aµ
A
B
e−i
e−i
D e−f
Figure 2.3: Scattering of two electrons in an electromagnetic potential.
Such scattering processes can be described by the exchange of virtual particles, Yukawa’s
force carriers. Even without understanding the details of the interaction, we can readily
identify one place where it should differ from the discussion above: the result must
somehow encode four-momentum conservation and not just energy conservation.
Our master formula for the differential cross-section, Eq. (2.8), is essentially a gen-
eralization to problems with more than one particle in the initial or final state. We
cannot derive the expressions for a scattering cross-section at high energies without
going through the machinery of quantum field theory. (This is not entirely true: see
Thomson, chapter 3 and section 5.1.) Instead, we will sketch the main results, then work
through the electrodynamics of spin-less particles as an example in the next lectures.
2.3. RELATIVISTIC SCATTERING 39
In quantum electrodynamics with scalar particles the transition amplitude Tfi for the
process A + B → C + D still takes the form in Eq. (2.24). Performing the integral using
incoming and outgoing plane waves φ = N e−ipx the result can be written as
Tfi = −i NA NB NC ND (2π)4 δ 4 (pA + pB − pC − pD ) M . (2.35)
where Ni are the plane wave normalization factors, which we will discuss shortly. The
δ-function takes care of energy and momentum conservation in the process. (Note that
the momentum vectors are four-vectors).
The quantity M is called the (Lorentz) invariant amplitude. It is computed using Feyn-
man diagrams. For topologies with n particles (counting both incident and final state),
the dimension of M is p4−n . Using the convention for the wave function normalization
described below, the invariant amplitude does not depend on arbitrary time intervals T
or normalization volumes V .
To find the transition probability we square the expression for Tfi ,
Z Z
2 2 −i(pA +pB −pC −pD )x 0
|Tfi | = |NA NB NC ND | |M| 2 4
d xe × d4 x0 e−i(pA +pB −pC −pD )x
(2.36)
Z
2 2 4 4
= |NA NB NC ND | |M| (2π) δ (pA + pB − pC − pD ) × lim d4 x (2.37)
T,V →∞ TV
2 2 4 4
= |NA NB NC ND | |M| (2π) δ (pA + pB − pC − pD ) × lim T V (2.38)
T,V →∞
Since we now have a δ-function over 4 dimensions (the four-momentum rather than just
the energy), the integral becomes proportional to both T and V . To get rid of them we
consider a transition probability per unit time and per unit volume:
|Tfi |2
Wfi ≡ lim
T,V →∞ T V
To use this result in our master formula, we now need to discuss a few remaining
ingredients, namely the normalization of the wave functions, the flux factor and the
phase space factor.
Above we defined the eigenstates of the free Hamiltonian to have unit normalization.
As we have seen in lecture 2 the eigenstates for free particles (for both the Schrödinger
equation and the Klein-Gordon equation) are plane waves
ψ(x, t) = N e−i(Et−x·p) . (2.40)
In contrast to wave packets the plane waves cannot be normalized over full space x
(which further on leads to problems when computing the square of δ-functions as above).
40 LECTURE 2. PERTURBATION THEORY AND FERMI’S GOLDEN RULE
The solution is to apply so-called box normalization: we choose a finite volume V and
normalize all wave functions such that
Z
ψ ∗ (x, t)ψ(x, t)d3 x = 1 . (2.41)
V
√
For the plane waves this gives N = 1/ V . Like the time interval T , the volume V is
arbitrary and must drop out once we compute an observable cross-section or decay rate.
For the classical wave function the density ρ = |ψ|2 so that the normalization gives one
particle per volume V . This normalization is not Lorentz invariant: under a Lorentz
transformation the volume element d3 x shrinks with a factor γ = E/m.
2
√ Klein-Gordon equation, we had ρ = 2|N | E. It is
For the plane wave solutions of the
now customary to still use N = 1/ V , but define the box normalization such that there
are 2E particles in the volume V . In other words,
ρ = 2E/V. (2.42)
We refer to this as 2E particles per unit volume. The factor 2E exactly cancels the
contraction of the volume, such that the number of particles in a given volume is now
Lorentz invariant.
Above we have introduced the Lorentz invariant amplitude without an explicit definition,
which is how we have found it in text books that do not derive the formalism with
field theory. Thomson takes an alternative approach: the plane wave functions √ in the
classical and relativistic case only differ by the normalization √ constant 2E. If we
label classical by ψ and relativistic by ψ 0 , then we have ψ 0 = 2Eψ. For a process
A + B + · · · → 1 + 2 + · · · , we now define the Lorentz-invariant matrix element in terms
of the wave functions with relativistic normalization,
where V is the perturbation to the free Hamiltonian (and not the volume!). As the name
suggests, with this construction M is Lorentz invariant. The non-relativistic transition
element that appears in Fermi’s golden rule is then related to M by
p
Mfi = 2E1 · 2E2 · · · 2EA · 2EB Vfi . (2.44)
In the final step to Fermi’s golden rule we introduced the density of final states ρ(E).
In the more general expression for the cross-section, it is the phase space factor that ac-
counts for the density of final states. It depends on the volume V and on the momentum
p of each final state particle.
2.3. RELATIVISTIC SCATTERING 41
where the product runs over all final state particles. To compute an actual number
for our experiment, we now convolute with experimental resolutions and integrate over
eventual particles or momentum components that we do not measure. (For example,
we often just measure the number of particles in a solid angle element dΩ.) For the
differential cross-section the question of the number of accessible states should then be
rephrased as “how many states fit in the momentum-space volume V d3 p”.
Assume that our volume V is rectangular with sides Lx , Ly , Lz . Using periodic boundary
conditions to ensure no net particle flow out of the volume we need to require that
Lx px = 2π~nx with nx integer. Hence, the total number of states in the range px to
px + dpx is dnx = Lx dpx /2π~. Since the total number of available states is n = nx ny nz ,
we find that the number of states with momentum between p and p + dp (i.e. between
(px , py , pz ) and (px + dpx , py + dpy , pz + dpz ) ) is:
V d3 p
dn = . (2.46)
(2π~)3
L
n
n= x
λx
Lz 2
1
Ly
Lx
Figure 2.4: Schematic calculation of the number of states in a box of volume V .
As explained above, in the relativistic case the wave functions are normalized such that
the volume V contains 2E particles. Therefore, the number of states per particle is:
V d3 p
# states/particle = (2.47)
(2π~)3 2E
If there are more particles in the final state, the density of states in Fermi’s rule must
account for each of those. Consequently, the phase space factor for a process with N
final state particles becomes
N
Y V d3 pf
dΦ = dLIPS = . (2.48)
f =1
(2π~)3 2Ef
42 LECTURE 2. PERTURBATION THEORY AND FERMI’S GOLDEN RULE
In exercise 2.3 you will show that (ignoring V ) the phase space factor is indeed Lorentz
invariant. We will omit the factors ~ in what follows.
The flux factor corresponds to the number of particles that pass each other per unit
area and per unit time. It can be most easily computed in a frame in which one of the
particles is not moving. Consider the case that a beam of particles (A) is shot on a
target (B), see Fig. 2.5.
target
beam
A B
The number of beam particles that pass through unit area per unit time is given by
|vA | nA . The number of target particles per unit volume is nB . For relativistic plane
waves the density of particles n is ρ = 2E
V
such that
2|pA | 2mB
flux = |vA | na nb = (2.49)
V V
(Remember that in relativity v = p/E, modulo a factor c. For the KG waves we had
indeed that the current density was j = ρ p/E.) In exercise 2.2 you will show that the
kinematic factor |pA |mB is actually Lorentz invariant and that this expression can be
rewritten as
q
flux = 4 (pA,µ pµB )2 − m2A m2B / V 2 (2.50)
The volume factor is not Lorentz invariant, but it will drop out later, as explained above.
The incident flux as defined here is not actually a certain number of particles per unit
surface per unit time per unit volume: we need to account for the fact that it is pro-
portional to the square of an energy. The factors of energy will be accounted for by the
other ingredients to the cross-section formula.
2.3. RELATIVISTIC SCATTERING 43
Putting this all together, we arrive at the formula to calculate a cross-section for the
process Ai + Bi → Cf + Df + ...:
1
dσfi = Wfi dΦ
fluxZ
1
Tfi = d4 x ψf∗ (x) V (x) ψi (x)
i~
|Tfi |2
Wfi = lim (2.51)
V,T →∞ T V
N
Y V d3 pf
dΦ =
f =1
(2π)3 2Ef
q
flux = 4 (pA · pB )2 − m2A m2B / V 2
In exercise 2.5 you will show that the cross-section is indeed independent on the volume
V.
Inserting the expression for the transition rate per unit time and volume, Eq. (2.39), we
find for the differential cross-section of the process A + B → C + D
(2π)4 δ 4 (pA + pB − pC − pD ) d 3 pC d3 pD
dσ = · |M|2 · (2.52)
(2π)3 2EC (2π)3 2ED
q
4 (pA · pB )2 − m2A m2B
The integrals of the flux factors are only over the spatial part of the outgoing four-
momentum vectors. The energy component has been integrated out, using the fact
that the outgoing particles
q are on the mass shell. Therefore, Ef is not an independent
variable, but equal to |pf |2 + m2f . This is important when performing integrals over
phase space.
In exercise 2.4 we calculate the integrals and flux factors in the centre-of-momentum
system, where pA + pB = pC + pD = 0. The result is
dσ 1 |pf |
= 2
|M|2 (2.53)
dΩ cm 64π s |pi |
(2π)4 δ 4 (pA − pC − pD ) d3 pC d3 pD
dΓ = · |M|2 · (2.54)
2EA (2π)3 2EC (2π)3 2ED
44 LECTURE 2. PERTURBATION THEORY AND FERMI’S GOLDEN RULE
√
which after integration of one of the momenta gives (4pi s → 2EA = 2mA )
dΓ 1
= 2 2
|pf | |M|2 (2.55)
dΩ cm 32π mA
Exercises
Exercise 2.1 (The Dirac δ-Function)
infinite
Consider a function defined by the following prescription
1/∆ for |x| < ∆/2 surface = 1
δ(x) = lim
∆→0 0 otherwise
0
The integral of this function is normalized
Z ∞
δ(x) dx = 1 (2.56)
−∞
These last two properties define the Dirac δ-function. The prescription above gives an
approximation of the δ-function. We shall encounter more of those prescriptions which
all have in common that they are the limit of a sequence of functions whose properties
converge to those given here.
(a) Starting from the defining properties of the δ-function, prove that
1
δ(kx) = δ(x) . (2.58)
|k|
Hint: Since the left hand side is Lorentz invariant, you can compute it in any frame.
Note that pA · pB is an inner product of the four-vectors, not the three-vectors.
Exercise 2.4 (AB → CD cross-section in the c.m.s. See also H&M, Ex. 4.2)
In this exercise we derive a simplified expression for the A + B → C + D cross-section
in the center-of-momentum frame.
(a) Start with the expression:
d3 pC d3 pD
Z
dΦ = (2π)4 δ 4 (pA + pB − pC − pD ) (2.63)
(2π)3 2EC (2π)3 2ED
Do the integral over d3 pD using the δ function and show that we can write:
1 p2f dpf dΩ
Z
dΦ = δ (EA + EB − EC − ED ) (2.64)
(2π)2 4EC ED
where we have made use of spherical coordinates (i.e. d3 pC = |pC |2 d|pC | dΩ) and
defined pf ≡ |pC |.
(b) In the C.M. frame
√ we have |pA | = |pB | = pi and |pC | = |pD | = pf . Furthermore,
in this frame s ≡ |pA + pB | = EA + EB ≡ W . Show that the expression becomes
(hint: calculate dW/dpf ):
Z
1 pf 1
dΦ = dW dΩ δ (W − EC − ED ) (2.65)
(2π)2 4 EC + ED
So that we finally get:
1 pf
dΦ = √ dΩ (2.66)
4π 2 4 s
(c) Show that the flux factor in the C.M. frame is:
√
flux = 4pi s (2.67)
46 LECTURE 2. PERTURBATION THEORY AND FERMI’S GOLDEN RULE
and hence that the differential cross-section for a 2 → 2 process in the center-of-
momentum frame is given by
dσ 1 pf
= 2
|M|2 (2.68)
dΩ cm 64π s pi
(a) The delta-function can have many forms. One of them is:
1 sin2 αx
δ(x) = lim (2.69)
α→∞ π αx2
Make this plausible by sketching the function sin2 (αx)/(παx2 ) for two relevant
values of α.
(b) Remember the Fourier transform,
Z +∞
1
f (x) = g(k) eikx dk
2π −∞
Z +∞ (2.70)
−ikx
g(k) = f (x) e dx
−∞
Use this to show that another (important!) representation of the Dirac delta
function is given by Z +∞
1
δ(x) = eikx dk (2.71)
2π −∞
Lecture 3
47
48 LECTURE 3. THE ELECTROMAGNETIC FIELD
waves travel, see the Feynman lectures, Vol.2, section 18. From now on we choose units
of charge such that we can set 0 = 1 and velocities such that c = 1. (We use so-called
’Heaviside-Lorentz rationalised units’: see section i.5.)
For what follows it is convenient to write the Maxwell equations in a covariant way (i.e.
in a manifestly Lorentz invariant way). As shown below we can formulate them in terms
of a single 4 component vector field, which we denote by Aµ = (V /c, A). As suggested
by the notation, the components of this field transform as a Lorentz vector.
You may prove for yourself that for any vector field A and scalar field V
From your electrostatics course you may remember that, because the rotation of E is
zero (which is the same as saying that E is a conservative vector field), all physics can be
derived by considering a scalar potential field V . The electric field becomes the gradient
of the potential, E = −∇V . The potential V is not unique: we can add an arbitrary
constant and the physics will not change. Likewise, because the divergence of the B
field is zero, we can always find a vector field A such that B is the rotation of A.
So, let’s choose a vector field A such that
B = ∇×A (3.8)
∂µ ∂ µ Aν − ∂ ν ∂µ Aµ = j ν . (3.10)
The current for electric charge j µ is a conserved current and transforms as a Lorentz
vector. (See also Feynman, Vol.2, section 13.6.) The derivative ∂ µ also transforms as
a Lorentz vector. Therefore, if the equation above is Lorentz covariant, then Aν must
transform as a Lorentz vector as well. Showing that the electromagnetic field indeed
transform this way is outside the scope of these lecture, but you may know that the
transformation properties of the electric and magnetic fields were an important clue
when Einstein formulated his theory of special relativity.
The expressions can be made even more compact by introducing the tensor
F µν ≡ ∂ µ Aν − ∂ ν Aµ . (3.11)
3.2. GAUGE TRANSFORMATIONS 49
such that
∂µ F µν = j ν . (3.12)
Just as the potential V in electrostatics was not unique, neither is the field Aµ . Imposing
additional constrains on Aµ is called choosing a gauge. In the next section we shall
discuss this in more detail. Written out in terms of the components E and B the
(4 × 4) matrix for the electromagnetic field tensor F µν is given by
0 −Ex −Ey −Ez
Ex 0 −Bz By
F µν = . (3.13)
Ey Bz 0 −Bx
Ez −By Bx 0
The field tensor is uniquely specified in terms of E and B. In other words, it does not
depend on the choice of the gauge.
∂λ
V0 = V +
∂t (3.14)
0
A = A − ∇λ.
or in terms of four-vectors
Aµ → A0µ = Aµ + ∂ µ λ (3.15)
do not change E and B.
If the laws of electrodynamics only involve the electric and magnetic fields, then, when
expressed in terms of the field A, the laws must be gauge ‘invariant’: physical observables
should not depend on λ. Sometimes we choose a particular gauge in order to make the
expressions in calculations simpler. In other cases, we exploit gauge invariance to impose
constraints on a solution, as with the photon below.
A common gauge choice is the so-called Lorentz gauge 1 . In exercise 3.3 you will show
that it is always possible to choose the gauge field λ such that Aµ satisfies the condition
With this choice Aµ becomes a conserved current. In the Lorentz gauge the Maxwell
equations simplify further:
Maxwell equations in the Lorentz gauge: ∂µ ∂ µ Aν = j ν (3.17)
However, as you will see in the exercise, Aµ still has some freedom since the Lorentz
condition fixes only ∂µ (∂ µ λ) and not ∂ µ λ itself. In other words a gauge transformation
of the form
Aµ → A0µ = Aµ + ∂ µ λ with 2λ = ∂µ ∂ µ λ = 0 (3.18)
is still allowed within the Lorentz gauge ∂µ Aµ = 0. Consequently, we can in addition
impose the Coulomb condition:
Coulomb condition: A0 = 0 (3.19)
(In combination with the Lorentz condition, also ∇ · A = 0 with this choice of gauge.)
This choice of gauge is not Lorentz invariant. This is allowed since the choice of the
gauge is irrelevant for the physics observables, but it is sometimes considered less elegant.
(Note that the second term is the complex conjugate of the first.) The four-vector aµ (p)
depends only on the momentum vector. It has four components but due to the gauge
transformation not all of those are physically meaningful. The Lorentz condition gives
which leads to
pµ aµ = 0. (3.24)
The Lorentz condition therefore reduces the number of independent complex components
to three. However, as explained above, we have not yet exhausted all the gauge freedom:
we are still free to make an additional shift Aµ → Aµ +∂ µ λ, provided that λ itself satisfies
the Klein-Gordon equation. If we choose it to be
With a bit of algebra we see that the result of the gauge transformation corresponds to
Note that aµ0 still satisfies the Lorentz condition only because p2 = 0 for a massless
photon.
As we have already seen, this additional freedom allows us to apply the Coulomb con-
dition and choose A0 = 0, or equivalently a0 (p) = 0. In combination with the Lorentz
condition this leads to
a·p=0 (3.28)
or p · A = 0.
At this point it is customary to write a(p) as a product of two terms
where is a vector of unit length and N (p) is real. The normalization N (p) depends
only on the magnitude of the momentum and corresponds to the energy density of the
wave. The vector depends only on the direction of p and is called the polarization
vector. Choosing the z axis along the direction of the momentum vector and imposing
the gauge conditions, the latter can be parameterized as
where ci and φi are all real and c21 + c22 = 1. We can remove one phase by moving the
origin. (Just look at how a shift of the origin affects the factors e±ipx .) Therefore, only
52 LECTURE 3. THE ELECTROMAGNETIC FIELD
two parameters of the polarization vector are physically meaningful: these are the two
polarization degrees of freedom of the photon.
Any polarization vector can be written as a (complex) linear combination of the two
transverse polarization vectors
You will show in exercise 3.4 that the circular polarization vectors + and − transform
under a rotation with angle θ around the z-axis (the momentum direction) as
+ → 0+ = e−iθ +
(3.33)
− → 0− = eiθ −
We now show that this means that these polarization states correspond to the two
helicity eigenstates of the photon.
You may remember from your QM course that the z component of the angular momen-
tum operator Jz is the generator of rotations around the z-axis. That means that for a
wavefunction ψ(x) the effect of an infinitesimal rotation around the z axis is given by
Comparing this to the effect of rotations on the polarization states above we now identify
+ with an m = +1 state and − with an m = −1 state.
Apparently, the polarization states belong to a representation of the rotation group:
they are spin states. Since we find ±1 for the Jz quantum number the photon must be
a spin-1 representation: it could not be spin zero, because than you would have only
3.4. ELECTRODYNAMICS IN QUANTUM MECHANICS 53
have a state with m = 0. And it could not have higher spin state, because there are no
degrees of freedom in the photon field that could be identified with larger values of m.
Since the photon is spin-1, one could have expected to find 3 spin states, namely for
mz = −1, 0, +1. You may wonder what happened to the mz = 0 component. This
component was removed when we applied the Coulomb gauge condition, exploiting
p2 = 0, leading to A · p = 0. For massive vector fields (or virtual photons!), there is no
corresponding gauge freedom and a component parallel to the momentum (a longitudinal
polarization) remains. Massive vector fields have one spin degree of freedom more.
Another way to look at this is to say that to define spin properly one needs to boost to
the rest frame of the particle. For the massless photon this is not possible. Therefore,
we can talk only about helicity (spin projection on the momentum) and not about spin.
The equivalent of the mz = 0 state does not exist for the photon.
Finally, we compute the electric and magnetic fields. Substituting the generic expression
for Aµ in the definitions of E and B and exploiting the coulomb condition A0 ≡ V = 0,
we find
E = i a p0 e−ipx + c.c.
(3.37)
B = −i (p × a) e−ipx + c.c.
Indeed, for the electromagnetic waves, the E and B fields are perpendicular to each
other and to the momentum, while the ratio of their amplitudes is 1 (or rather, c).
pµ → pµ − qAµ . (3.38)
in the equations of motion of the free particle. Written out in terms of the potential V
and vector potential A, the free Hamiltonian is then replaced by
1
H= (p − qA)2 + qV (3.39)
2m
It can be shown (see e.g. Jackson §12.1, page 575) that this indeeds leads to the Lorentz
force law, Eq. (3.1).
Performing the operator substitution, the Schrödinger equation for the Hamiltonian
above becomes
1 2 ∂
(−i∇ − qA) + qV ψ(x, t) = i ψ(x, t) (3.40)
2m ∂t
54 LECTURE 3. THE ELECTROMAGNETIC FIELD
Comparing this to the Schrödinger equation for the free particle, we note that we have
essentially made the substitution
−∇ → −∇ + iqA
(3.41)
∂/∂t → ∂/∂t + iqV
∂ µ → Dµ ≡ ∂ µ + iqAµ (3.42)
The gauge transformation leads to a change of the phase of the wave function. If λ
is not constant, then the change in phase is different at different points in space-time.
That is why we also call the gauge transformation a local phase transformation.
This result is at the heart of the application of gauge symmetries in quantum field
theory. Because, as we will see in more detail in Lecture 8, one can turn this argument
around: Since the phase of the wave function is not an observable, the equations that
describe the dynamics (a Schrödinger equation, or a Lagrangian) must be invariant
to such arbitrary phase transformations. If we impose this requirement, then we are
forced to introduce an Aµ field in the Hamiltonian via the substitution above and with
transformation properties defined above. In other words, the requirement of local phase
invariance imposes the form of the interaction!
path integral picture, the particle follows all possible trajectories to get from point x1
to point x2 , accumulating a phase eiS/~ , where S is the action along the path. Different
paths have different phases. It is only around the classical trajectory (obtained by
requiring the action to be minimal) that these phases interfere constructively. The size
of deviations along the classical trajectory is determined by ~.
As we have seen above the vector potential appears in the Schrödinger equation and
affects the wave function. In the presence of a magnetic field, the phase of the wave
function is changed along a trajectory according to
q x2
Z
∆α(A) = A(r, t) · dr (3.44)
~ x1
where the integral runs along the trajectory. (We do not prove this here. See also
Feynman Vol 2, section 15-5.) Although we do not need it here, for completeness we
also mention that the change in phase due an electric field is given by the integral of
the potential over the time:
q t2
Z
∆α(V ) = − V (r, t)dt (3.45)
~ t1
This last equation can easily be derived from the SE for a constant electric field. When
combined these two equations lead to a Lorentz covariant formulation if the integral is
performed over space and time.
Let us now consider Feynman’s famous two-slit experiment demonstrating the interfer-
ence between two electron trajectories. In the absence of external fields, the intensity at
a detection plate positioned behind the two slits shows an interference pattern. This is
most easily understood by considering the two ‘classical’ trajectories, depicted by ψ1 and
ψ2 in Fig. 3.1. The relative length of these trajectories differs as a function of the posi-
tion along the detection screen. The resulting phase difference leads to the interference
pattern. For a great description see chapter 1 of the “Feynman Lectures on Physics”
volume 3 (“2-slit experiment”) and pages 15-8 to 15-14 in volume 2 (“Bohm-Aharanov”).
Now consider the presence of a magnetic field in the form of vector field A. (We choose
the electric field zero, so A0 = 0.) Due to the A field, the phases of the two contributions
to the wave functions change,
ψ = ψ1 eiα1 (r,t) + ψ2 eiα2 (r,t) = ψ1 ei(α1 −α2 ) + ψ2 eiα2 .
(3.46)
The extra contribution to the relative phase is given by
Z Z I
q 0 0 q
α1 − α2 = dr1 A1 − dr2 A2 = dr 0 · A(r 0 , t)
~ ~
Z r1 r2
Z
q q q
= ∇ × A(r 0 , t) · dS = B · dS = Φ (3.47)
~ S ~ S ~
where we have used Stokes’ theorem to relate the integral around a closed loop to the
magnetic flux Φ through the surface. The magnetic field shifts the interference pattern
56 LECTURE 3. THE ELECTROMAGNETIC FIELD
detector
slits
Intensity
ψ2
source ψ
1 coil
Figure 3.1: The schematical setup of an experiment that investigates the effect of the presence
of an A field on the phase factor of the electron wave functions.
on the screen. In exercise 3.5 you will show that for a homogenous magnetic field this
leads to the same deflection as the classical force law.
Let us now consider the case that a very long and thin solenoid is positioned in the setup
of the two-slit experiment. Inside the solenoid the B-field is homogeneous and outside
it is zero (or sufficiently small). However, the A field is not zero outside the coil, as
illustrated in Fig. 3.2. The classical trajectories do not pass through the B field, but
they do pass through the A field, leading to a shift in the relative phase. Experimentally
it has been verified (in a technically difficult experiment) that the interference pattern
indeed shifts.
formulation. It is the only correct way to implement the Maxwell equations in quantum
mechanics. The gauge freedom may seem an undesirable feature now, but will turn out
to be a fundamental concept in our description of interactions.
Exercises
Exercise 3.1 (Maxwell equations)
Using the vector identity
∇ × (∇ × A) = −∇2 A + ∇ (∇ · A) (3.48)
(which one can prove using εijk εklm = δil δjm − δim δjl ) show that (with c = 1 and 0 = 1)
Maxwell’s equations can be written as:
∂µ ∂ µ Aν − ∂ ν ∂µ Aµ = j ν (3.49)
Exercise 3.5 (Deflection in magnetic field. Feynman, Vol II, sec. 15-5.)
We have stated above that with the minimal substitution recipe the Schrödinger equation
leads to the Lorentz force law. We have also stated (not proven) how the change of the
phase of a wavefunction due to a vector field A can be obtained by integrating the vector
58 LECTURE 3. THE ELECTROMAGNETIC FIELD
field along the trajectory. Let’s take these things for given and see if we can reproduce
the deflection of a particle in a magnetic field. Feynman beautifully illustrates that by
looking at the famous two-split experiment.
Consider the setup in Fig. 3.3. Particles with charge q, mass m and momentum p travel
from a source, via two slits, to a photographic plate. The interference of the two paths
leads to a diffraction pattern. The distance between the slits and the plate is L. Directly
behind the slits is a thin strip of magnetic field. The thickness of the strip is w and
w L. The B field is homogenous, coming out of the plane of the figure. We label the
coordinate along the photographic plate by x.
(a) For very small deflections, compute the deflection of the particles of the particles
using the Lorentz force law. Translate this into the displacement ∆x at the pho-
tographic plate. Hint: Assume that the plate is thin enough that direction of the
force is along the x-axis. The force lasts for a time w/v.
(b) Consider two classical (shortest distance) trajectories through the two slits (indi-
cated by 1 and 2 ). For small deflections, compute the phase shift between the
two trajectories as a function of x, in the absence of a magnetic field. Compute
the distance between two maxima in the diffraction pattern. Hint: The reduced
wavelength of the particles is λ/2π = ~/p.
(c) Assuming again small deflections use equation (3.47) to compute the increase in
phase shift between the two trajectories as a result of the B field. Translate the
phase shift in a shift ∆x of the diffraction pattern.
Lecture 4
Electromagnetic Scattering of
Spinless Particles
59
60LECTURE 4. ELECTROMAGNETIC SCATTERING OF SPINLESS PARTICLES
The sign of V is chosen such that compared to the kinetic energy it gets the same sign as
in the Schrödinger equation, Eq. (3.39). Since e2 is small (α = e2 /4π = 1/137) and we
only consider the Born level cross-section, we neglect the second order term, e2 A2 ≈ 0.
From Lecture 2 we take the general expression for the transition amplitude in the Born
approximation and insert the expression for V (x),
Z
Tf i ≡ −i d4 x φ∗f (x) V (x) φi (x)
Z
= −i d4 x φ∗f (x) (−ie) (Aµ ∂ µ + ∂µ Aµ ) φi (x). (4.7)
The second ∂µ operator on the right hand side acts on both Aµ and φ. However, we can
use integration by parts to write
Z Z
∗
∗ µ ∞
4 µ
∂µ φ∗f Aµ φi d4 x
d x φf ∂µ (A φi ) = φf A φi −∞ − (4.8)
Requiring the A field to be zero at t = ±∞, the first term on the left vanishes, such
that the transition amplitude becomes
Z
Tf i = −i (−ie) φ∗f (x) (∂µ φi (x)) − ∂µ φ∗f (x) φi (x) Aµ d4 x .
(4.9)
In this expression the derivatives no longer act on the field Aµ . Remember the definition
of the charge current density for the Klein-Gordon field of the electron, Eq. (1.31),
This is the expression for the transition amplitude for going from free particle solution i
to free particle solution f in the presence of a perturbation caused by an electromagnetic
field. Restricting ourselves to plane wave solutions of the unperturbed Klein-Gordon
equation,
φi = Ni e−ipi x and φ∗f = Nf∗ eipf x , (4.12)
4.2. COULOMB SCATTERING 61
Ze2
Tf i = i Ni Nf∗ (Ei + Ef ) 2π δ (Ef − Ei ) . (4.19)
|pf − pi |2
As we have seen before for a time-independent potential, we consider a time-averaged
transition rate,
|Tf i |2
Wf i = lim (4.20)
T →∞ T
where the time-averaging effectively takes care of one of the δ functions when taking the
square of the amplitude. The result is
!2
2
Ze (Ei + E f )
Wf i = |Ni Nf |2 2π δ (Ef − Ei ) 2 (4.21)
|pf − pi |
62LECTURE 4. ELECTROMAGNETIC SCATTERING OF SPINLESS PARTICLES
Working with normalization of the plane waves over a box with volume V , we have
|Ni Nf |2 = V 2 . The flux factor for a single particle is given by
2Ei 2|pi |
flux = |vi | = , (4.22)
V V
while the phase space factor is
V d3 pf
dLips = . (4.23)
(2π)3 2Ef
Inserting these expressions in our master formula for the cross-section, Eq. (2.8), we find
!2
2π δ (Ef − Ei ) Ze2 (Ei + Ef ) d 3 pf
dσ = (4.24)
2|pi | |pf − pi |2 (2π)3 2Ef
where pf,i now refers to the size of the three-momentum. Since Ef2 = m2 + p2f , we have
pf dpf = Ef dEf , and therefore,
Ef
δ (Ef − Ei ) dpf = δ (Ef − Ei ) dEf (4.26)
pf
dσ Z 2 Ei2 e4 Z 2 |Ei |2 α2
= = (4.28)
dΩ 64π 2 |pi |4 sin4 θ/2 4 |pi |4 sin4 θ/2
dσ Z 2 α2
= , (4.29)
dΩ 2
16Ekin sin4 θ/2
and an incident wave with momentum along the z-axis, the wave function for large r
takes the form
r→∞ eikr
φ(r) −−−→ eikz + f (θ, φ) (4.30)
r
where the first term is the incident ‘unscattered’ wave and the second term the ‘scattered’
wave. Expressed in terms of the f (θ, φ), the differential cross-section can be written as
dσ
= |f (θ, φ)|2 (4.31)
dΩ
The interference between the scattered and the unscattered wave leads to a ‘shadow’
behind the scattering potential. The flux that is missing in the shadow is exactly the
total scattered flux. This is expressed in the optical theorem, which states that the total
cross-section is proportional to value of f in the forward direction,
k
Imf (0) = σ. (4.32)
4π
See also appendix H of Aichison and Hey, and references therein.
B : K− D : K−
Aµ γ
A : π− C : π−
Figure 4.1: Leading order diagram for electromagnetic scattering of a charged kaon and a
charged pion.
We have seen above how a particle scatters in an external field. In this case the field is
not external as the particles scatter in each others field. How do we deal with this?
64LECTURE 4. ELECTROMAGNETIC SCATTERING OF SPINLESS PARTICLES
First consider a pion scattering in the vector field Aµ generated by the current of the
kaon. The transition current of the kaon is given by (see Eq. (4.13))
µ
jBD = −eNB ND∗ (pµB + pµD ) ei(pD −pB )x (4.33)
We now assume that the field generated by the kaon can be computed by inserting this
current in the Maxwell equations for the vector potential, i.e.
µ
∂ν ∂ ν Aµ = jBD (4.34)
where we have adopted the Lorentz gauge. (A proof that this indeed works requires the
full theory.) Since ∂ν ∂ ν eiqx = −q 2 eiqx , we can easily verify that the solution is given by
1 µ
Aµ = − j , (4.35)
q 2 BD
µ −1 BD 4 µ −gµν ν
Z Z Z
µ 4
Tf i = −i jAC Aµ d x = −i jAC 2 jµ d x = −i jAC j d4 x . (4.36)
q q 2 BD
−1
Z
Tf i = −ie 2
(NA NC∗ ) (pµA + pµC ) ei(pC −pA )x · 2 ·(NB ND∗ ) (pµB + pµD ) ei(pD −pB )x d4 x (4.37)
q
−1
Tf i = −ie2 (NA NC∗ ) (pµA + pµC ) ∗ B D
(2π)4 δ 4 (pA + pB − pC − pD )
2
(N B ND ) p µ + p µ
q
(4.38)
where the δ-function that takes care of four-momentum conservation appears. Usually
this is written in terms of the invariant amplitude M (sometimes called ‘matrix element’)
as
Tf i = −i NA NB NC∗ ND∗ (2π)4 δ 4 (pA + pB − pC − pD ) · M (4.39)
4.3. SPINLESS π − K SCATTERING 65
The signs and factors i are assigned such that the expressions for vertex factors and
propagator are also appropriate for higher orders. These are, in fact, our first set of
Feynman rules!
B D
ie(pB + pD )µ
−igµν
q2
ie(pA + pC )µ
A C
Figure 4.2: Feynman rules for the t-Channel contribution to electromagnetic scattering of
spinless particles.
not really need field theory anymore to compute cross-sections. We discuss the role of the
Lagrangian in more detail in Lecture 8. In appendix B we sketch how the propagators
and vertex factors are obtained from the Lagrangian.
We can now insert the invariant amplitude into the expression for the A + B → C + D
cross-section that we derived in lecture 2,
(2π)4 δ 4 (pA + pB − pC − pD ) d3 pC d 3 pD
dσ = q |M|2 3 3 . (4.41)
2 2
4 (pA · pB ) − mA mB 2 (2π) 2EC (2π) 2E D
In the centre-of-momentum frame (pA = −pB ) this expression became (see Eq. 2.53)
dσ 1 1 pf
= 2
|M|2 . (4.42)
dΩ 64π s pi
where s = (pA + pB )2 . As we consider elastic scattering, |pf |/|pi | = 1.
Finally, consider the limit of massless particles. Define p ≡ pA and p0 ≡ pC . In the
centre-of-momentum frame the four-vectors are given by
pµA = (|p|, p)
pµB = (|p|, −p)
pµC = (|p0|, p0 )
pµD = (|p0|, −p0 )
Define p ≡ |p| which, by four-vector conservation is also equal to |p0|. Define θ as the
angle between pA and pC , which means that cos θ = pA · pC /|pA ||pC | = pA · pC /p2 . We
then have
(pA + pC )µ gµν (pB + pD )ν = (pA )µ (pB )µ + (pA )µ (pD )µ + (pC )µ (pB )µ + (pC )µ (pD )µ
= 2p2 + p2 (1 + cos θ) + p2 (1 + cos θ) + 2p2
= p2 (6 + 2 cos θ)
Likewise, we get for q 2 ,
q 2 = (pA − pC )2
= p2A + p2C − 2(pA )µ (pC )µ
= −2p2 (1 − cos θ)
Consequently, we obtain for the invariant amplitude defined above
2
2 p (6 + 2 cos θ) 2 3 + cos θ
M = e =e . (4.43)
2p2 (1 − cos θ) 1 − cos θ
Inserting this in Eq. (4.42) gives (with α = e2 /4π),
2 2
α2 3 + cos θ
dσ 1 1 2 2 3 + cos θ
= e = (4.44)
dΩ 64π 2 s 1 − cos θ 4s 1 − cos θ
This is the leading order QED cross-section for the scattering of massless spin-0 particles
in the centre-of-momentum frame. In exercise 4.1 you will derive the formula for particles
with non-zero mass.
4.4. FORM FACTORS 67
In the previous sections we studied ‘point charges’, objects with their charge located in
an infinitely small region. If the charge distribution has a finite size, the differential cross-
section is different from that of a point source. Consequently, the measured differential
cross-section can tell us important information over the substructure of particles. For
example, most information about the structure of the proton has been obtained in
electron-proton scattering experiments, most notably at the Hera collider in Hamburg.
By following the same procedure as above for the static source, one can show that the
differential cross-section can be written as
dσ dσ
= |F (pi − pf )|2 (4.46)
dΩ dΩ point
where F (q) is called the form factor. It is given by the Fourier transform of the charge
distribution
Z
F (q) = ρ(x)eiq·x d3 x (4.47)
In real electron-proton scattering we also need to account for the spin and the magnetic
moment of the proton. The form factor will then become more complicated. You will
learn more about this in the Particle Physics II course.
We have seen that the negative energy state of a particle can be interpreted as the
positive energy state of its anti-particle. How does this effect energy conservation that
we encounter in the δ-functions? We have seen that the invariant amplitude takes the
form
Z
M∝ φ∗f (x) V (x) φi (x) dx (4.48)
Z
∗
M ∝ e−ipf x e−ikx e−ipi x dx
pi Z
= e−i(pi +k−pf )x dx
pf
= (2π)4 δ (Ei + ω − Ef ) δ 3 (pi + k − pf )
k
⇒ Energy and momentum conservation are
enforced by the δ-function.
Z
-p+ ∗
M ∝ e−ip− x e−i(−p+ +k)x dx
k
Z
= e−i(k−p+ −p− )x dx
p−
= (2π)4 δ (k − p− − p+ )
Z
p− ∗
M ∝ e−i(k−p+ )x e−i(p− )x dx
k
Z
= e−i(p− +p+ −k)x
-p+
= (2π)4 δ (p− + p+ − k)
4.5. PARTICLES AND ANTI-PARTICLES 69
Exercises
0 0 µ
(a) Express (ppA + pA )µ (pB + pB ) inpmA , mB , p and cos θ. (You may of course also
use EA = mA + p and EB = m2B + p2 .)
2 2
Hint: one method is to first just write down all fourvectors in these symbols. Take
the x, y-coordinates together, because we do not care about the azimuthal angle
(φ).Now you can either first add and them take the inner product, or vice versa.
(b) Do the same for q 2 = (pA − p0A )2 and for s = (pA + pB )2 .
(c) Write down the differential cross-section dσ/dΩ using Eq. (4.42). Note that this
result is more general than our ‘massless particle’ result in Eq. (4.44).
(d) Take the limit mB p and mB mA . Compare to the result for scattering of a
static source, Eq. (4.28).
Hint: Look at figure 4.2. Can you imagine additional ‘leading order’ diagrams? If so,
draw them!
M ∝ fπ α.
where α is the dimensionless coupling constant, and additional factors are dimen-
sionless as well. If you do not know anything else about the π 0 decay constant but
its dimension, what value would you use?
(e) Assuming that the π 0 is a uū + dd¯ state
(i) give the expression for the decay width (by adding up the amplitudes);
(ii) calculate the decay width expressed in GeV;
(iii) convert the rate into a mean lifetime in seconds.
(iv) How does the value compare to the Particle Data Group (PDG) value?
Remark: Do not be disappointed if your prediction is completely wrong! It turns
out that the π 0 lifetime is quite hard to compute.
Lecture 5
In Lecture 1 we have seen how the Klein-Gordon equation leads to solutions with neg-
ative energy and negative ‘probability density’. This is a consequence of the fact that
the wave equation is quadratic in ∂/∂t. In 1928 in an attempt to avoid this problem
Dirac developed a relativistic wave equation that is linear in ∂/∂t. Lorentz invariance
requires that such a wave equation is also linear in ∂/∂x. What Dirac found was an
equation that describes particles with spin- 21 , just what was needed for electrons. We
now think that all fundamental fermions are described by this wave equation. Dirac
also predicted the existence of anti-particles, an idea that was not taken seriously until
1932, when Anderson discovered the positron.
71
72 LECTURE 5. THE DIRAC EQUATION
So far, this is just classical electrodynamics. The classical spin S is nothing but the total
angular momentum of all bits and pieces that the particle is made up from. However,
as you remember from your quantum mechanics course, elementary particles also carry
intrinsic spin. Though we sometimes imagine it as a result of a charged particle spinning
around an axis, this interpretation actually falls short. In particular, the prediction of
the gyromagnetic ratio that would come out of this picture is wrong.
On the other hand, elementary particles do feel a torque in a magnetic field, as demon-
strated in the Stern-Gerlach experiment in 1922. So, in 1927 Pauli tried to address the
question of how to describe their magnetic moment in quantum mechanics.
Pauli considered a spin- 21 system. As you know, such a system has two values for the
eigenvalue of spin, namely ± 12 ~. An arbitrary spin wave function is a superposition of
the two eigenstates. Pauli represented it as a two-component vector, called a spinor,
a
χ = = aχ(1) + bχ(2) (5.4)
b
where the basis vectors
(1) 1 (2) 0
χ ≡ and χ ≡ (5.5)
0 1
represent the spin-up and spin-down state respectively. The hermitian operator that
measures spin, the spin operator S, satisfies the same algebra as for orbital angular
momentum in quantum mechanics, namely [Si , Sj ] = i~ijk Sk . In the basis above, S is
represented by 2x2 matrices. Choosing the z-axis as the quantization axis, S is given
by
~
S = σ (5.6)
2
where the Pauli spin matrices are
0 1 0 −i 1 0
σ1 = σ2 = σ3 = . (5.7)
1 0 i 0 0 −1
The matrices σi all have zero trace, are hermitian, and satisfy
σi σj = δij + iijk σk (5.8)
which implies as well that they anti-commute ({σi , σj } = 2δij ).
Pauli then proposed to replace the momentum operator in the Schrödinger equation by
a new operator in spinor space,
p −→ σ · p (5.9)
Take a careful look at what is written on the right-hand-side: it is an inner product of a
vector of 2x2 matrices with the momentum operator. The result is again a 2x2 matrix,
which is more apparent when written out,
pz px − ipy
σ·p= . (5.10)
px + ipy −pz
5.2. THE DIRAC EQUATION 73
You will show in an exercise that (σ · p)2 = |p|2 112 , where 112 is the 2x2 identity matrix.
Therefore, the Schrödinger equation for free spinors is just the ordinary Schrödinger
equation,
p2
d a 1 2 a a
i~ = (σ · p) = (5.11)
dt b 2m b 2m b
and the two spin states are degenerate in energy.
This is no longer the case if we introduce the vector field. Using again minimal substi-
tution, the Hamiltonian (a matrix in spinor space) for a particle in a vector field (V, A)
becomes
1
H= [σ · (p − qA)]2 + qV (5.12)
2m
It is a not entirely trivial exercise to show1 that this equation can be rewritten as
1 ~q
H= (p − qA)2 + qV − σ·B (5.13)
2m 2m
The Schrödinger equation with this Hamiltonian is called the Pauli-Schrödinger equa-
tion, or simply Pauli equation. Comparing this Hamiltonian to the Hamiltonian of the
classical spin, we find that the gyromagnetic ratio for a spin- 12 particle is
q
γspin-1/2 = (5.14)
m
exactly a factor 2 larger than for the classical picture of a spinning charge distribution.
The ratio of the magnetic moment relative to that of the classical case is called the
g-factor. For spin- 12 particles the Schrödinger-Pauli equation predicts g = 2. In QED
the magnetic moment is modified by higher order corrections. The predictions and
measurements of the magnetic moment of the electron and muon are so precise that
they make QED the most precisely tested theory in physics.
Pauli introduced the spin matrices in the Hamiltonian on purely phenomenological
grounds. As we shall see in the rest of this Lecture and the next, Dirac found a the-
oretical motivation: His construction of a wave equation that is linear in space and
time derivatives, leads (in its simplest form) to the description of spin- 12 particles and
anti-particles. As you will prove in exercise 5.8, the Pauli-Schrödinger equation can be
obtained as the non-relativistic limit of the equation of motion of Dirac particles in a
vector field Aµ .
closer to what Dirac actually did, see Griffiths, §7.1.) Consider the usual form of the
Schrödinger equation,
∂
i ψ = Hψ . (5.15)
∂t
The classical Hamiltonian is quadratic in the momentum. Dirac searched for a Hamil-
tonian that is linear in the momentum. We start from the following ansatz2 :
H = (α · p + βm) (5.16)
H 2 ψ = p2 + m2 ψ
(5.17)
where p2 + m2 is the eigenvalue. What should H look like such that these eigenvectors
exist? Squaring Dirac’s ansatz for the Hamiltonian gives
! !
X X
H2 = αi pi + βm αj pj + βm
i j
!
X X X
= αi αj pi pj + αi βpi m + βαi pi m + β 2 m2 (5.18)
i,j i i
!
X X X
= αi2 p2i + (αi αj + αj αi )pi pj + (αi β + βαi )pi m + β 2 m2
i i>j i
where we on purpose did not impose that the coefficients (αi , β) commute. In fact,
comparing to equation (5.17) we find that the coefficients must satisfy the following
requirements:
• α12 = α22 = α32 = β 2 = 1
• α1 , α2 , α3 , β anti-commute with each other.
With the following notation of the anti-commutator
Clearly, the αi and β cannot be ordinary numbers. At this point Dirac had a brilliant
idea, possibly motivated by Pauli’s picture of fermion wave functions as spinors: what
if the αi and β are matrices that act on a wave function that is a column vector? As
2
Note that α · p = αx px + αy py + αz pz .
5.3. COVARIANT FORM OF THE DIRAC EQUATION 75
we require the Hamiltonian to be hermitian (such that its eigenvalues are real), the
matrices αi and β must be hermitian as well,
Furthermore, we can show using just the anti-commutation relations and normalization
above that they all have eigenvalues ±1 and zero trace. It then follows that they must
have even dimension.
It can be shown that the lowest dimensional matrices that have the desired behaviour
are 4x4 matrices. (See exercise 5.6 and also Aitchison (1972) §8.1). The choice of the
matrices αi and β is not unique. Here we choose the Dirac-Pauli representation,
0 σ 11 0
α= and β= . (5.22)
σ 0 0 −11
Of course, we may expect that the final expressions for the amplitudes are independent
of the representation: all the physics is in the anti-commutation relations themselves.
Another frequently used choice is the Weyl representation,
−σ 0 0 11
α= and β= . (5.23)
0 σ 11 0
(iγ µ ∂µ − m) ψ = 0 (5.27)
This equation is called the Dirac equation. Note that ψ is a four-element vector. We
call it a bi-spinor or Dirac spinor. We shall see later that the solutions of the Dirac
76 LECTURE 5. THE DIRAC EQUATION
equation have four degrees of freedom, corresponding to spin-up and spin-down for a
particle and its anti-particle.
The Dirac equation is actually a set of 4 coupled differential equations,
4
" 3 #
for each X X
: i (γ µ )jk ∂µ − mδjk (ψk ) = 0
j=1,2,3,4
k=1 µ=0
. . . . 1 0 0 0 ψ1 0
. . . . 0 1 0 0 ψ2 0
or : i · ∂µ − · m
ψ3 = 0
. . . . 0 0 1 0
. . . . 0 0 0 1 ψ4 0
| {z }
γµ
Take note of the use of the Dirac (or spinor) indices (j, k = 1, 2, 3, 4) simultaneously
with the Lorentz indices (µ = 0, 1, 2, 3). As far as it concerns us, it is a coincidence that
both types of indices assume four different values.
To simplify notation even further we define the ‘slash’ operator of a four-vector aµ as
6 a = γµ aµ . (5.28)
The wave equation for spin- 12 particles can then be written very concisely as
(i 6 ∂ − m) ψ = 0 . (5.29)
{γ µ , γ ν } ≡ γ µ γ ν + γ ν γ µ = 2 g µν 114 (5.30)
5.5. ADJOINT SPINORS AND CURRENT DENSITY 77
where the identity matrix on the right-hand side is the 4 × 4 identity in bi-spinor space.
Text books usually leave such identity matrices away. However, it is important to realize
that the equation above is a matrix equation for every value of µ and ν. In particular,
g µν is not a matrix in spinor space. (In the equation, it is just a number!)
Using this result we find
2 2 2 2
γ 0 = 114 γ1 = γ2 = γ3 = −114 . (5.31)
γ µ† = γ 0 γ µ γ 0 (5.34)
γ 5 = iγ 0 γ 1 γ 2 γ 3 (5.35)
As we now work with matrices, we use hermitian conjugates rather than complex con-
jugates and find for the conjugate equation
∂ψ † 0 X ∂ψ †
k
− mψ † = 0
−i γ −i k
−γ (5.39)
∂t k=1,2,3
∂x
78 LECTURE 5. THE DIRAC EQUATION
However, we now see a potential problem: the additional minus sign in −γ k dis-
turbs the Lorentz invariant form of the equation. We can restore Lorentz covariance by
multiplying the equation from the right by γ 0 . Therefore, we define the adjoint Dirac
spinor
ψ = ψ†γ 0 . (5.40)
Note that the adjoint spinor is a row-vector:
ψ1
ψ2
Dirac spinor : ψ3
Adjoint Dirac spinor: ψ1, ψ2, ψ3, ψ4
ψ4
∂ψ 0 X ∂ψ
−i γ −i k
γ k − mψ = 0 (5.41)
∂t k=1,2,3
∂x
i∂µ ψγ µ + mψ = 0. (5.42)
Now we multiply the Dirac equation from the left by ψ and we multiply the adjoint
Dirac equation from the right by ψ:
i∂µ ψγ µ + mψ ψ = 0
ψ (i∂µ γ µ ψ − mψ)= 0
+
µ µ
ψ (∂µ γ ψ) + ∂µ ψγ ψ = 0
j µ = ψγ µ ψ (5.43)
then this current satisfies a continuity equation, ∂µ j µ = 0. The first component of this
current is simply
X4
0 0 †
j = ψγ ψ = ψ ψ = |ψi |2 (5.44)
i=1
which is always positive. This property was the original motivation of Dirac’s work.
The form Eq. (5.43) suggests that the Dirac probability current density transforms as
a contravariant four-vector. In contrast to the Klein-Gordon case, this is not so easy to
show since the Dirac spinors transform non-trivially. We will leave the details to the
textbooks.
5.6. BILINEAR COVARIANTS 79
The names ‘axial’ and ‘pseudo’ refer to the behaviour of these objects under the parity
transformation, x → −x. The scalar is invariant under parity, while the pseudo scalar
changes sign. The space components of the vector change sign under parity, while
those of the axial vector do not. We shall discuss the bilinear covariants and their
transformation properties in more detail in Lecture 7.
We now consider explicit expressions for the solutions of the Dirac equation, Eq. 5.27.
In exercise 5.1 you will show that each of the components of the Dirac wave satisfies
the Klein-Gordon equation. Therefore, we try to construct the solutions as plane wave
solutions
ψ(x) = u(p) e−ipx (5.47)
where u(p) is a 4-component column-vector that does not depend on x. After substitu-
tion in the Dirac equation we find the Dirac equation in the momentum representation,
(γ µ pµ − m) u(p) = 0 (5.48)
where uA and uB are two-component spinors. We can rewrite this as a set of two
equations
(σ · p) uB = (E − m) uA
(5.51)
(σ · p) uA = (E + m) uB ,
where σ ≡ (σ1 , σ2 , σ3 ).
Now consider a particle with non-zero mass in its restframe, p = 0. In this case, the two
equations decouple,
E uA = m uA
(5.52)
E uB = −m uB .
In order to extend the solution to particles with non-zero momentum, consider two Dirac
spinors for which the two upper coordinates uA (p) of u(p) are given by
(1) (2)
uA = χ(1) and uA = χ(2) . (5.54)
with the basis spinors χ1,2 defined in Eq. 5.5. Substituting this into the second equation
of (5.51) gives for the lower two components
(1,2) σ · p (1,2) σ · p (1,2)
uB = uA = χ . (5.55)
E+m E+m
To prove that these are indeed solutions of the equations, one can use the identity
(1,2) (1,2)
such that uA and uB also satisfy the first equation in (5.51). (See also exercise 5.2.)
In the Pauli-Dirac representation the Hamiltonian is given by
0 σ·p m112 0
H = α · p + βm = + (5.57)
σ·p 0 0 −m112
With a bit of algebra we obtain for our solution
h i !
p2 (1)
m + uA
Hu(1) = E+m
(1)
, (5.58)
E uB
which illustrates two things: In order that u(1) be a solution we need indeed that
E 2 = p2 + m2 . Furthermore, in the limit that p → 0, the energy eigenvalue is +m,
such that this is a positive energy solution. The calculation for u(2) is identical. Hence,
two orthogonal positive-energy solutions are
χ(1) χ(2)
(1) (2)
u =N σ·p and u =N σ·p (5.59)
E+m
χ(1) E+m
χ(2)
To gain slightly more insight, let’s write them out in momentum components. Using
the definition of the Pauli matrices we have
0 1 0 −i 1 0
σ·p= px + py + pz (5.64)
1 0 i 0 0 −1
82 LECTURE 5. THE DIRAC EQUATION
to find
(1) pz px − ipy 1 pz
(σ · p) uA = = (5.65)
px + ipy −pz 0 px + ipy
(2) (3) (4)
and similar for uA , uB , uB . The solutions can then be written as
1 0
0 1
E > 0 spinors u(1) (p) = N
pz
, u(2) (p) = N
px −ipy
E+m E+m
px +ipy −pz
E+m E+m
−pz −px +ipy
−E+m −E+m
−px −ipy pz
E < 0 spinors u(3) (p) = N −E+m , u(4) (p) = N −E+m
1 0
0 1
You can verify that these solutions are indeed orthogonal, i.e. that
The spinors u(p) of matter waves are solutions of the Dirac equation in momentum
space, Eq. (5.48). Replacing p with −p in the Dirac equation we find that our positive
energy anti-particle spinors satisfy another Dirac equation,
(6 p + m) v(p) = 0 (5.68)
5.9. NORMALIZATION OF THE WAVE FUNCTION 83
As for the Klein-Gordon case we choose a normalization such that there are 2E particles
per unit volume. The density is given by the time-component of the bi-spinor current
(see Eq. 5.44)
ρ(x) = ψ † (x) ψ(x) . (5.69)
Substituting the plane wave solution ψ = u(p) e−ipx , and integrating over a volume V
we find Z Z
3
ρd x = u† (p) eipx u(p) e−ipx d3 x = u† (p) u(p) V (5.70)
V V
Consequently, to find 2E particles per unit volume we must normalize such that
Explicit calculation for the positive energy solutions (s ∈ {1, 2}) gives
!
†
† T T (σ · p) (σ · p) (s)
u(s) u(s) = N 2 χ(s) χ(s) + χ(s) χ
(E + m)2
p2
2 2E
= N 1+ 2
= N2
(E + m) E+m
The computation for the positive energy anti-particle waves v(p) leads to the same
normalization. We can now write the orthogonality relations as (with r, s ∈ {1, 2})
†
u(r) u(s) = 2E δrs
†
(5.73)
v (r) v (s) = 2E δrs
We now consider the Dirac equation for the adjoint spinor u, v. Taking the hermitian
conjugate of Eq. (5.48) and multiplying on the right by γ 0 we have
u† γ µ† γ 0 pµ − u† γ 0 m = 0 (5.74)
Using that γ µ† γ 0 = γ 0 γ µ we then find for the Dirac equation of the adjoint spinor
u = u† γ 0 ,
u (6 p − m) = 0 (5.75)
84 LECTURE 5. THE DIRAC EQUATION
v (6 p + m) = 0 (5.76)
Using these results you will derive in exercise 5.4 the so-called completeness relations
X
u(s) (p) u(s) (p) = (6 p + m)
s=1,2
X (5.77)
v (s) (p) v (s) (p) = (6 p − m)
s=1,2
These relations will be used later on in the calculation of amplitudes with Feynman
diagrams. Note that the left-hand side is not an inner product. Rather, on both sides
we have a (4 × 4) matrix, or schematically
.
.
· (....) =
. γµ · pµ + 11 ·m
5.11 Helicity
The Dirac spinors for a given momentum p have a two-fold degeneracy. This implies
that there must be an additional observable with an operator that commutes with H
and p, and the eigenvalues of which distinguish between the degenerate states. It is
tempting to identify this extra quantum number with the intrinsic spin, but the relation
is subtle.
The spin operator in bi-spinor space is derines as S = 12 ~Σ, with
σ 0
Σ= . (5.78)
0 σ
In exercise 5.3 you will show that Σ does not commute with the Hamiltonian in
Eq. (5.57). We can also realize this by looking directly at our Dirac spinor solutions: If
spin is a good quantum number then those solutions should be eigenstates of the spin
operator,
Σ u(i) = s u(i) ?
5.11. HELICITY 85
where s is the spin eigenvalue. Now insert one of the solutions, for example u(1) ,
1 1
σ 0 0 ? s
0
0 σ pz / (E + m) = pz / (E + m)
(px + ipy ) / (E + m) (px + ipy ) / (E + m)
and you realize that this could never be true for arbitrary px , py , pz .
The orbital angular momentum operator is defined as usual as
L=r×p (5.79)
You will also show in exercise 5.3 that the total angular momentum
J = L + 21 Σ (5.80)
does commute with the Hamiltonian. Now, as we can choose an arbitrary axis to get the
spin quantum numbers, we can choose an axis such that the orbital angular momentum
vanishes, namely along the direction of the momentum. Consequently, we define the
helicity operator λ as
1 1 σ · p̂ 0
λ = Σ · p̂ ≡ (5.81)
2 2 0 σ · p̂
where p̂ ≡ p/|p|. We can interpret the helicity as the “spin component in the direction
of movement”. One can verify that indeed λ commutes with the Hamiltonian in (5.57).
As λ and H commute, they have a common set of eigenvectors. However, that does
not necessarily mean that our solutions u(i) are indeed also eigenvectors of λ. In fact,
with our choice above, they are only eigenvectors of λ if we choose the momentum along
the z-axis. The reason is that the two-component spinors χ(s) are eigenvectors of σ3
only. For other directions of the momentum, we would need to choose a different linear
combination of the u(i) to form a set of states that are eigenvectors for both H and λ.
Now, consider a momentum vector p = (0, 0, p). Applying the helicity operator on u(i)
gives
1 1 1 ±
(σ · p̂) u±
A = σ3 u±
A = ± uA
2 2 2
1 1 1 ±
(σ · p̂) u±
B = σ3 u±
B = ± uB
2 2 2
where the plus sign holds for u(1,3) and the minus sign for u(2,4) . So you see that indeed
u is an eigenvector of λ with eigenvalues ±1/2. Positive helicity states have spin and
momentum parallel, while negative helicity states have them anti-parallel.
It is not so difficult to derive the spinors that are eigenvectors of both λ and the Dirac
Hamilitonian for arbitrary momentum p. (See for instance §4.8.1 in Thomson.) We
86 LECTURE 5. THE DIRAC EQUATION
save you the algebra and just give the result. To simplify notation we switch to polar
coordinates,
p = (p sin θ cos φ, p sin θ sin φ, p cos θ) . (5.82)
The particle spinors for helicity +1/2 and helicity −1/2 become, respectively,
cos 2θ − sin 2θ
eiφ sin θ eiφ cos 2θ
u↑ = N 2 u = N , (5.83)
p θ ↓ p θ
E+m
cos 2
E+m
sin 2
p θ p θ
iφ iφ
E+m
e sin 2 − E+m e cos 2
We have seen above that although the probability density of the Dirac fields is positive,
the negative energy solutions just remain. Following the Feynman-Stückelberg interpre-
tation the solution with negative energy is again seen as the anti-particle solution with
positive energy. However, when it comes to the Feynman rules, there is an additional
subtlety for fermions.
In the case of Klein-Gordon waves the current of an anti-particle (j µ = 2|N |2 pµ ) gets a
minus sign with respect to the current of the particle, due to reversal of 4-momentum.
This cancels the change in the sign of the charge and that is how we came to the nice
property of ‘crossing’: simply replace any anti-particle by a particle with opposite mo-
mentum. For fermions this miracle does not happen: the current does not automatically
change sign when we go to anti-particles. As a result, if we want to keep the convention
that allows us to replace anti-particles by particles, we need an additional ‘ad-hoc’ minus
sign in the Feynman rule for the current of the spin- 12 anti-particle.
5.13. THE CHARGE CONJUGATION OPERATION 87
This additional minus sign between particles and anti-particles is only required for
fermionic currents and not for bosonic currents. It is related to the spin-statistics con-
nection: bosonic wavefunctions are symmetric, and fermionic wavefunctions are anti-
symmetric. In field theory3 the extra minus sign is a result of the fact that bosonic
field operators follow commutation relations, while fermionic field operators follow anti-
commutation relations. This was realized first by Pauli in 1940.
γ µ∗ = M −1 γ µ M . (5.90)
where we have used the definition of the adjoint spinor (see Lecture 6) and defined the
charge conjugation matrix C = M γ 0 . It can be shown (see Halzen and Martin exercise
5.6) that in the Pauli-Dirac representation a possible choice of M is
1
−1
M = Cγ 0 = iγ 2 = . (5.93)
−1
1
3
See Aitchison & Hey, 3rd edition §7.2
88 LECTURE 5. THE DIRAC EQUATION
Interpreting the probability current as a charge current, we define the electron current
as
µ
j−e = −eψγ µ ψ (5.94)
The current of the charge conjugate wave function is then
µ C
j−e = −eψC γ µ ψC = . . . = eψγ µ ψ (5.95)
Exercises
Exercise 5.1 (From Dirac to Klein-Gordon)
Each of the four components of the Dirac equation satisfies the Klein Gordon equation,
(∂µ ∂ µ + m2 ) ψi = 0. Show this explicitly by operating on the Dirac equation from the
left with (iγ ν ∂ν + m).
Hint: For any aµ and bν we can write γ ν γ µ aν bµ = 12 (γ ν γ µ aν bµ + γ µ γ ν aµ bν ) by just
’renaming’ indices’. Now take the special case that b = a and the aµ commute. Then
we can write γ ν γ µ aν aµ = 21 (γ ν γ µ + γ µ γ ν )aµ aν . Now use the anti-commutation relation
of the γ-matrices.
Exercise 5.3 (See also exercise 5.4 of H& M and exercise 7.8 of Griffiths)
The purpose of this problem is to demonstrate that particles described by the Dirac
equation carry “intrinsic” angular momentum (S) in addition to their orbital angular
momentum (L). We will see that L and S = Σ/2 are not conserved individually but
that their sum is.
(a) Consider the Hamiltonian that leads to the Dirac equation,
H = α · p + βm
[H, L] = −i α × p (5.96)
5.13. THE CHARGE CONJUGATION OPERATION 89
where L = x × p.
Hint: To do
P this efficiently use the Levi-Civita tensor to write out the cross product
as Li = j,k ijk xj , pk . Now evaluate the commutator [H, Li ].
where the operator Σ (see also Eq. (5.78)) and α in the Pauli-Dirac representation
were
σ 0 0 σ
Σ= and α =
0 σ σ 0
Hint: Use the commutation relation
P for the Pauli spin matrices [σi , σj ] = 2iijk σk
(which follows from σi σj = i k ijk σk ).
(c) Use the result in (b) to show that
[H, Σ] = 2 i α × p (5.97)
We see from (a) and (c) that the Hamiltonian commutes with J = L + 12 Σ.
Hint: For s = 1, 2 take the solution u(s) from Eq. (5.59) and write out the row-
vector for us using the explicit form of γ 0 P
in the Dirac-Pauli representation. Then
write out the matrix u(s) u(s) and use that s=1,2 χ(s) χ(s)† = 112 . Finally, note that
E 112 −σ · p
6p = (5.98)
σ · p −E 112
Tr γ 5 6 a 6 b 6 c 6 d = −4 i εαβγδ aα bβ cγ dδ
where εαβγδ = +1 (−1) for an even (odd) permutation of 0,1,2,3; and 0 if two
indices are the same.
Hint:
P Write the innerPproducts as sums over i and j. Use that σi σj = δij +
i k i,jk σk . Use that i,j ijk ai bj = (a × b)k .
(b) Prove the identity in Eq. (5.56), i.e.
(σ · p) (σ · p) = |p|2 112
In the following we concentrate on the uA component because in exercise 5.2 you have
shown that in the non-relativistic limit the other component is small.
(a) Starting from the Dirac equation in momentum space Eq. (5.51),write down the
equations for uA and uB after minimal substitution.
(b) In coordinate space p and E are operators. Therefore, they do not commute with
A and V . That means that unlike in the case for free particles (exercise 5.2)
you cannot just eliminate uB ! However, in the non-relativistic limit and assuming
|qV | m we can use the approximation E − qV + m ≈ 2m. Use this to eliminate
uB and obtain the Dirac equation for uA
where Ekin = E − m.
(c) Use the Pauli vector identity, Eq. (5.99), to show that
σ · (p × A + A × p) uA = −i~σ · (∇ × A) uA (5.102)
where on the right hand side the derivative ∇ works only on A and not on uA .
Therefore, we can replace it with B = ∇ × A.
92 LECTURE 5. THE DIRAC EQUATION
(e) Using these results, show that the Dirac equation for an electron with charge q = −e
in the non-relativistic limit in an electromagnetic field Aµ = (A0 , A) reduces to
the Schrödinger-Pauli equation
d 1 2 e 0
i ψA = (p + eA) + σ · B − eA ψA , (5.103)
dt 2m 2m
Spin-1/2 Electrodynamics
(H0 + V ) ψ = E ψ . (6.2)
One can either start from the Dirac equation in terms of α and β, or, work towards
that form by multiplying the Dirac equation on the left by γ 0 . The result is
Eψ = γ 0 γ k pk + γ 0 m ψ − eγ 0 γµ Aµ ψ .
(6.3)
| {z } | {z }
H0 =α·p+βm V
Compare this to the transition amplitude for spin-0 scattering, Eq. 4.7: The wave
function now has four components and the perturbation potential V (x) is a (4 × 4)
93
94 LECTURE 6. SPIN-1/2 ELECTRODYNAMICS
matrix. We take a hermitian conjugate of the wave ψ, rather than its complex conjugate.
The transition amplitude is still just a scalar.
Substituting the expression for V (x) we obtain
Z
ψf† (x) −eγ 0 γµ Aµ (x) ψi (x)d4 x
Tf i = −i
Z (6.6)
= −i ψ f (x) (−e) γµ ψi (x)Aµ (x) d4 x
After inserting the plane wave decomposition ψ(x) = u(p)e−ipx , the transition current
becomes
jfµi = −euf γ µ ui ei(pf −pi )x . (6.9)
jfµi = ( uf ) γµ ui (6.10)
µ
jAC
uA uC
Just as we did for the scattering of spinless particles, we obtain the vector potential Aµ
by using the Maxwell equation with the transition current of one of the two particles
(say ‘particle AC’) as a source. That is, we take
µ
2Aµ = jAC .
6.1. FEYNMAN RULES FOR FERMION SCATTERING 95
BD −1 µ −gµν ν 4
Z Z
4 µ
Tf i = −i jµ jAC d x = −i jBD j dx (6.11)
q 2 q 2 AC
which is symmetric in terms of particle BD and AC. Inserting the expressions of the
plane wave currents using Eq. (6.9) we obtain
−gµν
Z
Tf i = −i −e uC γ µ uA ei(pC −pA )x · 2
· −e uD γ ν uB ei(pD −pB )x d4 x (6.12)
q
Performing the integral (and realizing that nothing depends on x except the exponen-
tials) leads us to the expression
From the matrix element we can now read off the Feynman rules. Again, as for the
spinless case, the various factors are defined such that the rules can also be applied to
higher order diagrams.
1 1 ui uf
ie (pf + pi )µ ieγ µ
Figure 6.1: Diagrams for a spin-0 (left) and spin- 12 (right) particle with charge e interacting
with the electromagnetic field.
The rules for the vertex factors for spin-0 and spin- 12 particles are shown side-by-side
in Fig. 6.1. A spinless electron can interact with Aµ only via its charge. The coupling
is proportional to (pf + pi )µ . However, an electron with spin can also interact with
the magnetic field, via its magnetic moment. As you will prove in exercise 6.1, we can
rewrite the Dirac current as
1
uf γ µ ui = uf (pf + pi )µ + iσ µν (pf − pi )ν ui
(6.15)
2m
96 LECTURE 6. SPIN-1/2 ELECTRODYNAMICS
where the tensor σ µν was defined in Eq. 5.46. This formulation of the current is called the
‘Gordon decomposition’. We observe that in addition to the contribution that appears
for the spinless wave, there is a new contribution that involves the factor iσ µν (pf − pi ).
In the non-relativistic limit this leads indeed to a term proportional to the magnetic
field component of Aµ , just as you would expect from a magnetic moment.
−iM = −igµν
q2
ieuC γ µ uA
e− : uA e− : u C
Figure 6.2: Lowest order Feynman diagram for e− µ− scattering.
For a given value of µ and ν the currents are just complex numbers. (The γ-matrices
are sandwiched between the bi-spinors.) Therefore, we can reorder them and write the
amplitude as
e4 X
|M|2 = [(uC γ µ uA ) (uC γ ν uA )∗ ] [(uD γµ uB ) (uD γν uB )∗ ] (6.18)
q 4 µν
We have factorized the right hand side into two tensors, each of which only depends on
one of the leptons. We call these the polarized lepton tensors.
Up to now we have ignored the fact that the particle spinors come in two flavours,
namely one for positive and one for negative helicity. Assuming that we do not measure
the helicity (or spin) of the incoming and outgoing particles, the cross-section that we
need to compute is a so-called ‘unpolarized cross-section’:
6.2. ELECTRON-MUON SCATTERING 97
The spin summation is unfortunately rather tedious. The rest of the lecture is basically
just the calculation to do this!
First, take a look at the complex conjugate of the transition current that appears in the
tensor. Since it is just a (four-vector of) numbers, complex conjugation is the same as
hermitian conjugation. Consequently, we have
[uC γ ν uA ]∗ = [uC γ ν uA ]†
h i† h i
† 0 ν † ν† 0
= uC γ γ uA = uA γ γ uC (6.22)
= uA γ 0 γ ν † γ 0 uC = [uA γ ν uC ]
In other words, by reversing the order of the spinors, we can get rid of the complex
conjugation and find X
Lµν
e = (uC γ µ uA ) (uA γ ν uC ) (6.23)
e spin
98 LECTURE 6. SPIN-1/2 ELECTRODYNAMICS
Next, we apply what is called Casimir’s trick. Write out the matrix multiplications in
the tensors explicitly in terms of the components of the matrices and the incoming spins
s and outgoing spins s0 ,
X X X (s0 ) µ (s) (s) (s0 )
Lµν
e = ν
uC,k γkl uA,l uA,m γmn uC,n (6.24)
s0 s klmn
All of the factors on the right are just complex numbers, so we can manipulate their
order and write this as
X X (s0 ) (s0 ) µ X (s) (s)
Lµν
e = uC,n uC,k γkl ν
uA,l uA,m γmn (6.25)
klmn s0 s
Now remember the completeness relation, Eq. (5.77), that we derived in the previous
lecture1 , X
u(s) u(s) =6 p + m (6.26)
s
Lµν µ ν
e = Tr [(6 pC + m) γ (6 pA + m) γ ] (6.28)
You now realize why we made you compute the traces of products of γ-matrices in
exercise 5.5. We briefly repeat here the properties that we need:
• In general, for matrices A, B and C and any complex number z
– Tr(zA) = z Tr(A)
– Tr (A + B) = Tr(A) + Tr(B)
– Tr (ABC) = Tr (CAB) = Tr (BCA)
• For γ-matrices (from the anti-commutator γ µ γ ν + γ ν γ µ = 2g µν ):
– Tr(odd number of γ-matrices = 0)
– Tr (γ µ γ ν ) = 4 g µν
1 µν
for anti-fermions this gives an overall “−” sign in the tensor: Lµν
e → −Le for each particle →
anti-particle.
6.2. ELECTRON-MUON SCATTERING 99
– Tr(γ α γ β γ µ γ ν ) = 4 g αβ g µν − g αµ g βν + g αν g βµ
Using the first rule we can write out the tensor as a sum of traces,
Lµν µ ν
e = Tr [(6 pC + m) γ (6 pA + m) γ ]
The last two terms vanish because they contain an odd number of γ-matrices. For the
second term (‘case 2’) we find
Tr [mγ µ mγ ν ] = m2 Tr [γ µ γ ν ] = 4 m2 g µν . (6.30)
Finally, for the first term (‘case 1’) we have
Tr [6 pC γ µ 6 pA γ ν ] ≡ Tr γ α pC,α γ µ γ β pA,β γ ν
= T r γ α γ µ γ β γ ν pC,α pA,β
(6.31)
= 4 g αµ g βν − g αβ g µν + g αν g βµ pC,α pA,β
The expression for the muon tensor is obtained with the substitution (pA , pC , me ) →
(pB , pD , mµ ),
Lµν = 4 pµD pνB + pνD pµB + m2µ − pD · pB g µν
µ (6.33)
To compute the contraction of the two tensors, which appears in the amplitude, we just
write everything out
µ ν ν µ
Lµν 2
µν
µ
· 4 pDµ pBν + pDν pBµ + m2µ − pD · pB gµν
e Lµν = 4 pC pA + pC pA + me − pC · pA g
= 16 (pC · pD ) (pA · pB ) + (pC · pB ) (pA · pD ) − (pC · pA ) (pD · pB ) + (pC · pA ) m2µ
= 32 (pA · pB ) (pC · pD ) + (pA · pD ) (pC · pB ) − m2e (pD · pB ) − m2µ (pA · pC ) + 2m2e m2µ
Combining everything we obtain for the square of the unpolarized amplitude for electron-
muon scattering
e4 h
|M|2 = 8 (pC · pD ) (pA · pB ) +
q4
i
(pC · pB ) (pA · pD ) − m2e (pD · pB ) − m2µ (pA · pC ) + 2m2e m2µ (6.34)
100 LECTURE 6. SPIN-1/2 ELECTRODYNAMICS
We now consider the ultra-relativistic limit and ignore the rest masses of the particles.
The amplitude squared then becomes
e4 h i
|M|2 ' 8 4 (pC · pD ) (pA · pB ) + (pC · pB ) (pA · pD ) (6.35)
q
Furthermore, we define the Mandelstam variables
s ≡ (pA + pB )2 = p2A + p2B + 2 (pA · pB ) ' 2 (pA · pB )
t ≡ (pD − pB )2 ≡ q 2 ' −2 (pD · pB ) (6.36)
u ≡ (pA − pD )2 ' −2 (pA · pD )
where the approximation on the right follows in the ultra-relativistic limit (m ≈ 0).
From energy-momentum conservation (pµA + pµB = pµC + pµD ) we have
(pA + pB )2 = (pC + pD )2 pA · pB = pC · pD
(pD − pB )2 = (pC − pA )2 =⇒ pD · p B = pC · pA (6.37)
(pA − pD )2 = (pB − pC )2 pA · p D = pB · pC
which gives
1 1 1 2
(pA · pB ) (pC · pD ) = s s= s (6.38)
2 2
4
1 1 1
(pA · pD ) (pC · pB ) = − u − u = u2 (6.39)
2 2 4
(6.40)
Inserting this in the amplitude, we find
s2 + u 2
2 4
|M| ' 2 e (6.41)
t2
Finally, as we did for the spinless scattering in Lecture 4, consider again the scattering
process in the centre-of-momentum system. The four-vectors can then be written as
pµA = (|pA |, pA ) pµB = (|pA |, −pA )
pµC = (|pC |, pC ) pµD = (|pC |, −pC )
2 4 + (1 + cos θ)2
|M| ' 8 e4 . (6.43)
(1 − cos θ)2
Inserting this in the expression for the differential cross-section (which we obtained after
integrating over the final state momenta in exercise 2.4) we find
dσ 1 1 2 α2 4 + (1 + cos θ)2
= |M| ' (6.44)
dΩ c.m. 64π 2 s 2s (1 − cos θ)2
with α ≡ e2 /4π.
µ− µ− pC
ν
ieuD γ uB pA
θ
−iM = 2 q2
q
pD pB
ieuC γ µ uA
e− e−
Figure 6.3: e− µ− → e− µ− scattering. Left: the Feynman diagram. Right: definition of
scattering angle in C.M. frame.
M e− (p0A )e+ (p0B ) → µ+ (p0C )µ− (p0D ) = M e− (p0A )µ− (−p0C ) → e− (−p0B )µ− (p0D )
(6.45)
102 LECTURE 6. SPIN-1/2 ELECTRODYNAMICS
⇒
µ− (pB ) µ− (pD ) e− (pC ) µ− (pD )
⇒
e− (pA ) µ+ (−pB )
p0A = pA
p0B = −pC ⇐
p0C = −pB
p0D = pD
e+ (−pC ) µ− (pD )
Figure 6.4: Illustration of crossing. Use the anti-particle interpretation of a particle with the
4-momentum reversed in order to related the Matrix element of the “crossed” reaction to the
original one.
In other words, we can use the original computation of the amplitude provide that we
relabel the momenta as follows:
pA = p0A pB = −p0C pC = −p0B pD = p0D
Consequently, the Mandelstam variables of the ’original’ particle diagram are
s ≡ (pA + pB )2 = (p0A − p0C )2 ≡ t0
t ≡ (pD − pB )2 = (p0C + p0D )2 = s0 (6.46)
u ≡ (pA − pD )2 = (p0A − p0D )2 = u0
Using the result in Eq. (6.41) the amplitude squared for the two processes are then
2 s2 + u2
|M|e− µ− →e− µ− = 2 e4 ”t-channel”: q2 = t
t2
2 u02 + t02
|M|e− e+ →µ− µ+ = 2 e4 ”s-channel”:
s02 q2 = s
It is customary to label these as the t-channel and the s-channel process, because we
have q 2 = t and q 2 = s, respectively.
We can express the momenta in the centre-of-momentum frame in terms of an initial
momentum p and a scattering angle θ, where θ is now the angle between the incoming
6.3. CROSSING: THE PROCESS E− E+ → µ− µ+ 103
e− (p0A ) and the outgoing µ− (p0C ). The expressions for u0 , s0 and t0 are identical to those
in (6.42).
We immediately get for the matrix element:
2 t02 + u02
4 4 2
|M|c.m.
=2e = e 1 + cos θ (6.47)
s02
The differential cross-section becomes
dσ α2
1 + cos2 θ
= (6.48)
dΩ 4s
Finally, to calculate the total cross-section for the process we integrate over the az-
imuthal angle φ and the polar angle θ:
4π α2
σ e+ e− → µ+ µ− = (6.49)
3 s
Note that the ‘shape’ of the angular distribution does not depend on the available energy,
but that the total cross-section scales as 1/s: the higher the cms energy, the smaller the
cross-section. If you look back to our original formulation of the golden rule, you’ll find
that the 1/s dependence comes from the density of the incoming waves. The faster the
relative velocity, the shorter the particles are in each others vicinity!
Figure 6.5 shows a comparison of the kinematic factors in the differential cross-section
of the t-channel process e− µ− → e− µ− and the s-channel process e− e+ → µ− µ+ for
spin-0 and spin- 12 leptons. For the t-channel process the difference is only visible in the
very backward region, while in the s-channel process there is a constant offset.
t-channel s-channel
α 2 dΩ
α 2 dΩ
4s dσ
4s dσ
3
103
spin-1/2
102 spin-0 2
10
1
1
0
-1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1
cos(θ) cos(θ)
Figure 6.5: Leading order QED differential cross-section dσ/dΩ divided by α2 /4s as function
of cos θ for the t-channel process e− µ− → e− µ− (left) and the s-channel process e− e+ → µ− µ+
(right) in the ultra-relativistic limit (me = mµ = 0).
Figure 6.6 shows a table copied from Halzen and Martin with the kinematic factors
for important leading order QED processes. These processes are related by crossing.
The interference terms follow via crossing procedure as well, provided that you add up
amplitudes (not amplitudes squared).
104 LECTURE 6. SPIN-1/2 ELECTRODYNAMICS
Figure 6.6: Leading order QED processes and their relations via crossing. From Halzen and
Martin, “Quarks and Leptons”.
In the computation of the e− µ− above we have seen only a subset of the Feynman
rules for QED. As an example of things we missed, consider the annihilation process
e+ e− → γγ. (Draw it!) To compute the cross-section for this process we need more
Feynman rules, namely those for the electron propagator and those for external photon
lines. We now briefly summarize the rules for QED. You can find these in more detail in
the textbooks, e.g. in appendix D of Griffiths and on the inside of the cover of Halzen
and Martin, or Thomson.
spin-0: nothing
incoming particle: u
outgoing particle: u
spin- 12 :
incoming anti-particle: v (6.50)
outgoing anti-particle: v
incoming: µ
spin-1:
outgoing: ∗µ
We have seen the photon polarization vectors in Lecture 3. Both the spin- 12 and spin-1
external lines carry also an index for the helicity. In calculations for cross-sections or
decays in which we measure the spin, we need explicit forms of the Dirac spinors and
the photon polarization vectors. However, often we sum over all incoming and outgoing
spins (’spin averaging’) and we can use the completeness relations.
For the internal lines (the propagators) we have
i
spin-0:
q2 − m2
i(6 q + m)
spin- 12 :
q 2 − m2 (6.51)
µν
−ig
massless:
q2
spin-1: µν µ ν 2
masssive: i [−g + q q /m ]
q 2 − m2
Exercises
Exercise 6.1 (The Gordon decomposition)
A spinless electron can interact with Aµ only via its charge; the coupling is proportional
to (pf + pi )µ . An electron with spin, on the other hand, can also interact with the
magnetic field via its magnetic moment. This coupling involves the factor iσ µν (pf − pi ).
The relation between the Dirac current and the Klein-Gordon current can be studied as
follows:
106 LECTURE 6. SPIN-1/2 ELECTRODYNAMICS
Exercise 6.2
Can you easily obtain the cross-section of the process e+ e− → e+ e− from the result of
e+ e− → µ+ µ− ? If yes: give the result, if no: why not?
2
(c) Use the principle of crossing to find |M| for e+ e− → π + π −
(Note the extra minus sign that appears from the 3rd crossing rule.)
(d) Determine the differential cross-section dσ/dΩ for e+ e− → π + π − in the centre-of-
momentum of the e+ e− -system.
Lecture 7
In 1896 Henri Becquerel studied the effect of fluorescence, which he thought was related
to X-rays that had been discovered by Wilhelm Röntgen. To test his hypothesis he
wrapped a photographic plate in black paper and placed various phosphorescent salts
on it. All results were negative until he used uranium salts. These affected photographic
plates even when put in the dark, such that the effects clearly had nothing to do with
fluorescence. Henri Becquerel had discovered natural radioactivity, and thereby the weak
interaction.
We know now that the most nuclear decays are the result of the transition of a neutron
to an electron, a proton and an anti-neutrino,
p+
n
e−
νe
or in a formula,
n → p + e− + ν e . (7.1)
A ‘free’ neutron has a lifetime of about 15 minutes, but the lifetime of various weakly
decaying isotopes spans a very wide range.
107
108 LECTURE 7. THE WEAK INTERACTION
The listed decay modes are the dominant decay modes. Other decay modes exist, but
they contribute marginally to the total decay width. As we have seen before, the lifetime
of a particle is inversely proportional to the total decay width,
hbar
τ = . (7.2)
Γ
We have also seen that the decay width to a particular final state is proportional to the
matrix element squared. For example, for the two-body decay A → B + C we had (in
particle A’s rest frame)
|M|2
Z
pB
Γ(A → B + C) = dΦ = |M|2 (7.3)
2EA 8πm2A
• The µ+ is a lepton and therefore does not couple to the strong interaction. It can-
not decay to an electron and photon, as the electromagnetic interaction conserves
lepton flavour. Its dominant decay is via the weak interaction to an electron and
neutrinos.
Considerations like these explain the gross features in the hierarchy of lifetimes. How-
ever, as you can also judge from the wide range in lifetimes of particles that decay
weakly, kinematic effects must be important as well.
Besides the fact proper that the weak interaction unlike the electromagnetic and strong
interaction does not ‘honour’ the quantum numbers for quark and lepton flavour, the
weak interaction is special in at least two more ways:
• it violates parity symmetry P . Until 1956, when the parity violating aspects of
the weak interaction were demonstrated, physicists were convinced that at least
at the level of fundamental interactions our world was left-right symmetric;
• in the quark sector, it even violates CP symmetry. That means, because of CP T
invariance, that it also violates T (time-reversal) symmetry. As we shall see, the
existence of a third quark family was predicted from the observation that neutral
Kaon decays exhibit CP violation.
1. e2 = 4πα is replaced by GF
e− e− 2. 1/q 2 is removed
110 LECTURE 7. THE WEAK INTERACTION
“Introduction to High Energy Physics”, 3rd edition, appendix D) that for the decay
n → pe− ν e
• S, P and T interactions imply that the helicity of the e− and the ν e have the same
sign;
• V and A interactions imply that they have opposite sign.
Fermi had assumed that the weak interaction was of the V type. In a number of
experiments performed in the late fifties it was established that the weak interaction
was a combination of V and A. Before we look at that in more detail, we need to discuss
the concept of parity.
7.3 Parity
Parity, or (space) inversion, is the operation that multiplies all spatial coordinates by
−1, so x → −x. It is closely related to reflection in a mirror: the parity operation
is identical to a reflection in a plane through the origin, followed by a rotation under
180 degrees around an axis through the origin perpendicular to the mirror. Therefore,
for systems that are rotation and translation invariant, the two are equivalent. When
illustrating parity violation in pictures, we usually use an image with a reflection in a
mirror. Yet, when formulating the effect of parity in a physics theory, we work with
space inversion.
Now, consider a process φi → φf for some initial state i and final state f . The relation
between i and f given by an operator that describes the time evolution,
φf = Ûf i φi (7.8)
(We can look at the process at any time scale. So Û can just be a continuous function
of time.) Denoting the parity operation by P̂ we can also consider the mirror process,
characterised by φ0i = P̂ φi and φ0f = P̂ φf . We define the process to be ‘symmetric under
parity’ when it does not make any difference whether we first transform φi to its mirror
image and then look at its time evolution, φi → φ0i → φ0f , or first wait for the system
to evolve and then reflect it, φi → φf → φ0f . Or, in terms more common in quantum
mechanics, the process is symmetric under parity when P̂ and Û commute,
[Û , P̂ ] = 0 (7.9)
Because for small times t we have Û (t) = e−iHt/~ ≈ 1 − iHt/~, it follows that such P̂
also commutes with the Hamiltonian. This definition of a symmetry is not limited to
mirror symmetry, but holds for any operator: if an operator Q̂ commutes with H then
it is called a symmetry operator.
If P̂ and H commute, then they have a common set of eigenvectors. If we consider
eigenvectors with energy E that are not degenerate (that is, there is no other state with
112 LECTURE 7. THE WEAK INTERACTION
equal energy) then this immediately implies that these states have definite parity: they
are eigenstates of the parity operator and there is an observable property (a quantum
number) associated with the parity operation.
If we apply the parity operator twice, then we put the system back in its original state.
Consequently, if p is the eigenvalue of our state under the parity operator, then p2 = 1.
(Strictly speaking, the system would be in the same state even if we had changed the
wave function by an arbitrary phase. However, for simplicity we will not deal with the
minor complications that this introduces.) Therefore, the eigenvalue is either +1 or −1.
We call such states states of even and odd parity respectively.
Until 1956 all the known laws of physics were invariant under inversion symmetry. At
the scale of elementary particles our world was perfectly left-right symmetric. This
symmetry was well tested for the electromagnetic and strong interaction and it was
generally assumed that it held for the weak interaction as well.
Since all all our leptons, mesons, and baryons are characterised by different masses (e.g.
by different eigenvalues of the total Hamiltonian) they all have definite parity: they are
either odd or even under the parity operation. (You will find their quantum numbers for
parity in the PDG.) These facts, definite parity for all (stable and meta-stable) particles,
and parity conservation in known interactions, is exactly what lead in the early fifties
to what was called the ‘theta-tau puzzle’.
The θ and the τ were charged particles with strangeness one that decayed through the
weak interaction to two and three pions respectively,
θ+ → π + + π 0
(7.10)
τ + → π + + 2π 0 or 2π + + π −
The pions were all known to have parity −1. Then, assuming parity to be conserved
in these processes the theta had even parity and the tau odd parity. However, what
was truly strange is that the theta and tau were otherwise seemingly identical particles:
they had the same mass and same lifetime.
After verifying that there had never been any experimental tests of parity conservation
in the weak interaction, Lee and Yang hypothesized in 1956, that the tau and theta
were actually the same particle, and that the weak interaction was responsible for the
apparent violation of parity. They also proposed a number of experiments that could
establish parity violation in weak decays directly. Within half a year two of these
experiments were performed (Wu et al. (1957), Garmin, Lederman and Weinrich (1957))
and the parity violating character of the weak interaction was firmly established.
7.4. COVARIANCE OF THE WAVE EQUATIONS UNDER PARITY 113
Hence βψ 0 (−x, t) does satisfy the Dirac equation. Consequently, one choice for the
parity operation is
ψ 0 (x0 , t) = βψ(−x0 , t) (7.18)
Again, we could insert any constant phase factor in the transformation. By convention,
we choose that factor to be one.
We now look at the solutions to the Dirac equation in the Pauli-Dirac representation.
In this representation, we have for the β = γ 0 matrix:
0 11 0
γ = (7.19)
0 −11
Consequently, the parity operator has opposite sign for the positive and negative energy
solutions. In other words, fermions and anti-fermions have opposite parity. With our
choice of the phase of the parity transformation, fermions have positive parity and anti-
fermions have negative parity.
What does this mean for the currents in the interactions? Under the parity operation
we find
S: ψψ → ψγ 0 γ 0 ψ = ψψ Scalar
Which combination of V and A appears in the weak interaction was established with a
famous experiment by Goldhaber.
152
Eu + e− −→ 152
Sm∗ + νe
direction of travel: ←− −→
−1/2 −1 +1/2
spin configuration A: ◦ ⇐ ⇐= ⇒
+1/2 +1 −1/2
spin configuration B: ◦ ⇒ =⇒ ⇐
For neutron decay, the measured vector and axial vector couplings are CV = 1.000 ±
0.003, CA = 1.260 ± 0.002
The Fermi theory has a 4-point interaction, unlike the Yukawa theory: there is no
propagator to ‘transmit’ the interaction from the lepton current to the hadron current.
However, we know now that forces are carried by bosons:
• the weak interaction is carried by the massive W , Z bosons. For vector boson
with non-zero mass M the propagator is
−i (g µν − q µ q ν /M 2 )
.
M 2 − q2
W
g g
G g2
strength: ∼ √F
2
∼ 8MW2
It is an experimental fact that the strength of the coupling of the weak interaction, the
coupling constant “g”, is identical for quarks and leptons of all flavours. For leptons
this is sometimes called ‘lepton-universality’.
e2
How “weak” is the weak interaction? For the electromagnetic coupling we have α = 4π ≈
g 2
1/137. It turns out that the weak coupling is equal to αw = 4π ≈ 1/29. We see that at
low energies, the weak interaction is ‘weak’ compared to the electromagnetic interaction
not because the coupling is small, but because the propagator mass is large! At high
energies q 2 & MW2
the weak interaction is comparable in strength to the electromagnetic
interaction.
7.7. MUON DECAY 117
e−(p’) νµ (k)
µ−(p) µ−
νe(k) (p)
W e−(p’)
νµ (k’) νe(−k’)
Figure 7.1: Muon decay: left: Labelling of the momenta, right: Feynman diagram. For the
spinor of the outgoing anti-particle we use uνe (−k 0 ) = vνe (k 0 ).
Using the Feynman rules we can write for the matrix element:
g 1 1 g 1
M = √ u(k) γ µ 1 − γ 5 √ u(p0 ) γµ 1 − γ 5 v(k 0 )
u(p)
2
2 |{z} 2 |{z} MW 2 | {z } 2 | {z }
outgoing νµ incoming µ |{z} outgoing e outgoing ν e
propagator
(7.25)
Next we square the matrix element and sum over the spin states, just like we did for
e+ e− → µ+ µ− . Then we use again Casimir’s tric, as well as the completeness relations,
to convert the sum over spins into a trace. The result is:
2 2
2 1 X 2 1 g
· Tr γ µ 1 − γ 5 (6 p 0 + me ) γ ν 1 − γ 5 6 k 0
|M| = |M| = 2
2 Spin 2 8MW
· Tr γµ 1 − γ 5 6 k γν 1 − γ 5 (6 p + mµ )
G g2
Now we use some more trace theorems (see below) and also √F
2
= 8MW2 to find the result:
2
|M| = 64 G2F (k · p0 ) (k 0 · p) (7.26)
Intermezzo: Trace theorems used (see also Halzen & Martin p 261):
Tr γ µ 1 − γ 5 6 a γ ν 1 − γ 5 6 b · Tr γµ 1 − γ 5 6 c γν 1 − γ 5 6 d = 256 (a · c) (b · d)
118 LECTURE 7. THE WEAK INTERACTION
E = muon energy
E0 = electron energy
ω0 = electron neutrino energy
ω = muon neutrino energy
First we evaluate the expression for the matrix element. Working in the rest frame of
the muon and ignoring the mass of the electron and the neutrinos, we find (Eq. 7.48 in
exercise 7.3),
2
|M| = 64 G2F (k · p0 ) (k 0 · p) = 32 G2F m2 − 2mω 0 mω 0
(7.29)
where m is the muon mass. Inserting this in the expression for the differential decay
width, we obtain
1 2 16G2F
(m2 − 2mω 0 mω 0 dQ
dΓ = |M| dQ = (7.30)
2E m
where we used that E = m in the muon rest frame. To obtain the total decay width we
must integrate over the phase space,
16G2F
Z Z
1 2
(m2 − 2mω 0 mω 0 dQ
Γ= |M| dQ = (7.31)
2E m
The integrand only depends on the neutrino energy ω 0 . So, let us first perform the
integral in dQ over the other energies and momenta:
d3 p0 d3 k0 d3 k
Z Z
1
dQ = δ (m − E 0 − ω 0 − ω) δ 3 (p0 + k0 + k)
other 8 (2π)5 E 0 ω0 ω
d3 p0 d3 k0
Z
1
= δ (m − E 0 − ω 0 − ω)
8 (2π)5 E 0ω0ω
where θ is the angle between the electron and the electron neutrino. We choose the
z-axis along k0 , the direction of the electron neutrino. From the equation for ω we
derive:
−2E 0 ω 0 sin θ −ω dω
dω = √ dθ ⇔ dθ = (7.33)
2 E 02 + ω 02 + 2E 0 ω 0 cos θ E 0 ω 0 sin θ
| {z }
ω
so that we get:
G2 m
Z
Γ= F 3 (m − 2ω 0 ) ω 0 dω 0 dE 0 (7.37)
(2π)
Before we do the integral over ω 0 we have to determine the limits:
νe e−
• maximum electron neutrino energy:
ω 0 = 12 m νµ
• minimum electron neutrino energy:
ω 0 = 12 m − E 0 νµ e−
νe
Therefore, we obtain for the distribution of the electron energy in the muon rest frame
Z 1m
G2F m G2F m2 02 E0
dΓ 2
0 0 0
= (m − 2ω ) ω dω = E 3−4 . (7.38)
dE 0 (2π)3 12 m−E 0 12π 3 m
VOLUME 14, NUMBER 12 PHYSICAL REVIE%' LETTERS 22 MARCH 1965
of internal consistency of the data over various mus for fabricating the chambers; J. Williams
momentum and angular regions lead to larger for computer programming; and C. Carlson,
uncertainty and a preliminary result of p =0.747 S. Herzka, B. Palatnick, and S. Stein for gen-
+ 0.005. eral assistance in the experiment.
120 %'e wish to thank Dr. G. Sutter for help in LECTURE 7. THE WEAK INTERACTION
the early phases of the experiment; F. Sippach *Work supported in part by the U. S. Office of Naval
for the design of the electronic system', G. Dore- Research under Contract No. Nonr-266(72}.
X
LA
CU
I-
IO 20 30 40 50
POS I TRON MOME NTU M M eV/c
FIG. 4. {a}Experimental points for magnetic-field settings, normalized to the overlap region. The solid line is
the theoretical spectrum for p= 0.75. The Michel spectrum, 24
1 GF m 5
Γ =≡zjl2x —
p(x)dx =12x +p[(32/3)x —8x ])dx, (7.39)
where x is the positron momentum
τ 192 π 3
divided by its maximum value, has been corrected for internal radiation,
bremsstrahlung, and ionization loss. (b) The deviation of experimental points from the best-fit theoretical curve
for p= 0.747, showing typical experimental errors for four points. Curves for p= 0.737 and 0.757 are shown for
The measurement of the muon lifetime is the standard method to determine the coupling
comparison.
constant of
452 the weak interaction. The muon lifetime has been measured to be τ =
2.19703 ± 0.00004 µs. From this we derive for the Fermi coupling constant GF =
(1.16639 ± 0.00002) · 10−5 GeV−2 .
The strong and electromagnetic interaction do not couple to currents that connect lep-
tons or quarks of different flavour. These interactions conserve the type of lepton or
quark at the interaction vertex.
This is different for the weak interaction. As the W is charged, it necessarily couples to
a current that contains two particles that differ by one unit in charge. For the quarks
7.8. QUARK MIXING 121
and leptons in the standard model, the Feynman diagrams for the interactions are
νe e− νµ µ− ντ τ−
g g g
W W W
u d c s t b
g g g
W W W
Leptons and quarks are usually ordered in three ‘generations’ to show how the weak
interaction couples:
νe νµ ντ u c t
Leptons: Quarks: (7.40)
e µ τ d s b
If the only couplings of the W are those shown in the Feynman diagrams, then the
lightest hadrons with a strange quark (such as the K − which is a sū bound state) would
be stable. However, K − mesons do decay, for instance to a muon and a muon neutrino:
u
?? g µ−
−
K W νµ
s
This decay looks a lot like that of the π − , which is a dū bound state:
u
g g µ−
π−
W νµ
d
Experimentally the K − decay is found to have a much smaller decay width than the
pion decay.
In 1963 Nicola Cabibbo provided a solution that explained most available data on strange
hadron decay by presenting the d quark in the current that couples to the W as a linear
combination of a d quark and an s quark:
d → d0 = d cos θc + s sin θc
(7.41)
s → s0 = −d sin θc + s cos θc
where θc is a ‘mixing’ angle now known as the Cabibbo angle. In matrix representation
the mixing can be written as
d0
cos θc sin θc d
= (7.42)
s0 − sin θc cos θc s
122 LECTURE 7. THE WEAK INTERACTION
u d u d0 u d u s
g g g cos θc g sin θc
⇒ = +
W W W W
Due to the mixing the amplitudes for pion and kaon decay contain factors cos θc and
sin θc :
1. Pion decay
u gcos θ g µ−
− −
π → µ νµ π−
Γπ− ∝ G2F cos2 θc W νµ
d
2. Kaon decay
u gsin θ g µ−
− − −
K → µ νµ K
ΓK − ∝ G2F sin2 θc W νµ
s
A proper calculation gives for the ratio of the decay rates
3 2
Γ (K − ) m2K − m2µ
mπ
−
≈ tan2 θc · (7.43)
Γ (π ) mK m2π − m2µ
From the experimental result on the lifetime ratios, the Cabibbo angle is then found to
be
θC ≈ 13.0◦ (7.44)
Even though Cabibbo’s theory explained strange decays, it did not quite get everything
right. The proposed quark mixing would allow neutral kaons (sd¯ combinations) to decay
to muons, via the amplitude represented by this Feynman diagram:
µ+
d W+
K 0 u νµ
−
s W
µ−
According to Cabibbo’s calculation this decay should have an appreciable rate, but it
was never found! An explanation was provided by Glashow, Iliopoulis and Maiani in
1970: They hypothesised the existence of the charm (c) quark, contributing with a
diagram
7.8. QUARK MIXING 123
µ+
d +
W
K 0 c νµ
−
s W
µ−
The up and charm quark amplitudes have opposite sign, which leads to a nearly van-
ishing decay rate. This mechanism, which is now known as the GIM mechanism, was
the first well-motivated prediction for a fourth quark. The charm quark was discovered
3 years later.
Including charm quarks, the couplings for the first two generations are:
u d c s u s c d
g cos θc g cos θc g sin θc g sin θc
W W W W
| {z } | {z }
Cabibbo “favoured00 decay Cabibbo “suppressed00 decay
The flavour eigenstates u, d, s, c arethe mass eigenstates of the total Hamiltonian
u c
describing quarks. The states , are the eigenstates of the weak interaction
d0 s0
Hamiltonian, which affects the decay of the particles. By convention mixing is presented
for ’down’ quarks, but in fact that choice is arbitrary: We could also consider the mixing
matrix to mix the u and c quarks.
Of course, the story of quarks did not stop with the discovery of the charm quark. In
1964 Cronin and Fitch had shown in experiments that CP symmetry is violated in
neutral kaon decays. Kobayashi and Maskawa found a solution in 1973: They extended
Cabibbo’s picture of quark mixing with a third family of quarks,
d0
Vud Vus Vub d
s0 = Vcd Vcs Vcb s (7.45)
b0 Vtd Vts Vtb b
| {z }
CKM matrix
The mixing matrix VCKM is a 3 × 3 unitary matrix. This matrix is not uniquely defined
since the phases of the quark field can be chosen arbitrarily. If the phases are ’‘absorbed’
in the quark fields, the matrix can be parametrized by four real parameters, which are
usually chosen to be three mixing angles between the quark generations θ12 , θ13 , θ23 ,
and one complex phase δ,
s13 e−iδ
c12 c13 s12 s13
VCKM = −s12 c23 − c12 s23 s13 eiδ c12 c23 − s12 s23 s13 eiδ s23 c13 (7.46)
iδ iδ
s12 s23 − c12 c23 s13 e −c12 s23 − s12 c23 s13 e c23 c13
where sij = sin θij and cij = cos θij .
Kobayashi and Maskawa realized that the fact elements of VCKM can have a non-trivial
complex phase — i.e. a phase that can not be removed by redefining the phase of the
quark fields — leads to CP violation in charged current decays. For CP violation to
occur this way, at least three generations of quarks are required. The bottom and top
quark were eventually discovered in 1977 and 1994, respectively.
In case neutrino particles have a non-zero mass, mixing occurs in the lepton sector as
well. Just like the down-type quarks were chosen to describe mixing in the quark sector,
the neutrinos are chosen for the lepton sector:
νe U11 U12 U13 ν1
νµ = U21 U22 U23 ν2 (7.47)
ντ U31 U32 U33 ν3
| {z }
PMNS-matrix
Exercises
Exercise 7.1 (Helicity versus chirality. See also H&M exercise 5.15)
(a) Write out the chirality operator γ 5 in the Dirac-Pauli representation.
(b) The helicity operator is defined as λ = 21 Σ · p̂, where p̂ is a unit vector along the
momentum and Σ is
σ 0
Σ= .
0 σ
Show that in the ultra-relativistic limit (E m) the helicity operator and the
chirality operator have the same effect on a spinor solution, i.e.
χ(s) χ(s)
5 5
γ ψ=γ σ·p ≈ 2λ σ·p = 2λ ψ
E+m
χ(s) E+m
χ(s)
Due to the fact that the quarks in the pion are not free particles we cannot just apply
the Dirac formalism for free particle waves. However, we know that the interaction is
transmitted by a W − and therefore the coupling must be of the type: V or A. (Also,
the matrix element must be a Lorentz scalar.) It turns out the decay amplitude has the
form:
GF
M = √ (q µ fπ ) u(p)γµ 1 − γ 5 v(k)
2
where pµ and k µ are the 4-momenta of the muon and the neutrino respectively, and q is
the 4-momentum carried by the W boson. fπ is called the decay constant.
126 LECTURE 7. THE WEAK INTERACTION
(c) Can the pion also decay to an electron and an electron-neutrino? Write down the
Matrix element for this decay.
Would you expect the decay width of the decay to electrons to be larger, smaller,
or similar to the decay width to the muon and muon-neutrino?
Base your argument on the available phase space in each of the two cases.
The measured lifetime of the pion is τπ = 2.6 · 10−8 s which means that fπ ≈ mπ . An
interesting observation is to compare the decay width to the muon and to the electron:
2 2
Γ(π − → e− ν e ) m2π − m2e
me
= ≈ 1.2 · 10−4 !!
Γ(π − → µ− ν µ ) mµ m2π − m2µ
(d) Can you give a reason why the decay rate into an electron and an electron-neutrino
is strongly suppressed in comparison to the decay to a muon and a muon-neutrino.
Consider the spin of the pion, the handedness of the W coupling and the helicity
of the leptons involved.
2k · p0 = m2 + 2p · k 0
In the next two lectures we discuss the theory of the electroweak interaction, the so-
called “Glashow-Salam-Weinberg model”. This theory can be formulated starting from
the principle of local gauge invariance.
8.1 Symmetries
Symmetries play a fundamental role in particle physics. A symmetry is always related
to a quantity that is fundamentally unobservable. In general one can distinguish1 four
types of symmetries:
• permutation symmetries: These lead to Bose-Einstein statistics for particles with
integer spin (bosons) and to Fermi-Dirac statistics for particles with half integer
spin (fermions). The unobservable is the absolute identity of a particle;
• continuous space-time symmetries: translation, rotation, acceleration, etc. The re-
lated unobservables are respectively: absolute position in space, absolute direction
and the equivalence between gravity and acceleration;
• discrete symmetries: space inversion, time inversion, charge inversion. The unob-
servables are absolute left/right handedness, the direction of time and an absolute
definition of the sign of charge;
• unitary symmetries or internal symmetries, also called ‘gauge invariance’: These
are the symmetries discussed in this lecture. As an example of an unobservable
quantity think of the complex phase of a wave function in quantum mechanics.
The relation between symmetries and conservation laws is expressed in a fundamental
theorem by Emmy Noether: each continuous symmetry transformation under which the
Lagrangian is invariant in form leads to a conservation law. Invariances under external
1
T.D. Lee: “Particle Physics and Introduction to Field Theory”
127
128 LECTURE 8. LOCAL GAUGE INVARIANCE
operations as time and space translation lead to conservation of energy and momentum,
and invariance under rotation to conservation of angular momentum. Invariances under
internal operations, like the shift of the complex phase of wave functions, lead to con-
served currents, or more specific, conservation of charge. We discuss the application of
Noether’s theorem to phase transformations in section 8.4.
In the standard model the elementary interactions of the quarks and leptons (electro-
magnetic, weak and strong) are all the result of gauge symmetries. Starting from a
Lagrangian that describes free quarks and leptons, the interactions can be constructed
by requiring the Lagrangian to be symmetric under particular transformations. The
idea of local gauge invariance will be discussed in this lecture and will be applied in the
unified electroweak theory in the next lecture.
In classical mechanics the equations of motion can be derived using the variational prin-
ciple of Hamilton. This principle states that the action integral S should be stationary
under arbitrary variations of the so-called generalized coordinates qi . For a pedagogical
discussion of the principle of least action read the Feynman lectures, Vol.2, chapter 19.
Generalized coordinates are coordinates that correspond to the actual degrees of freedom
of a system. As an example, consider a swinging pendulum in two dimensions. We could
describe the movement of the weight of the pendulum in terms of both its horizontal
coordinate x and its vertical coordinate y. However, only one of those is independent
since the length of the pendulum is fixed. Therefore, we say that the movement of the
pendulum can be described by one ‘generalized’ coordinate. We could choose x or y,
but also the angle of the pendulum with the vertical axis (usually called the amplitude).
We denote generalized coordinates with the symbol q and call the evolution of q with
time a trajectory or path.
The Lagrangian of the system can be defined as the kinetic energy minus the potential
energy,
L(q, q̇, t) = T (q̇) − V (q) , (8.1)
where the potential energy only depends on q (and eventually t) and the kinetic energy
only on the generalized velocity q̇ = dq/dt. We denote the action (or ‘action integral’)
of a path that starts at t1 and ends at t2 with
Z t1
S(q) = L(q, q̇, t) dt . (8.2)
t0
Hamilton’s principle states that the actual trajectory q(t) followed by the system is
the trajectory q(t) that minimizes the action. (It is said that the action is ‘stationary’
around this trajectory.) This is equivalent to requiring that for a given point q, q̇ on this
8.3. LAGRANGIAN DENSITY FOR FIELDS 129
trajectory, the change in the action following from a small deviation δq, δ q̇ is zero:
You will show in exercise 8.1 that for each of the coordinates qi , this leads to the so-called
Euler-Lagrange equation of motion
d ∂L ∂L
= . (8.4)
dt ∂ q˙i ∂qi
∂L ∂L
ṗi = with pi = , (8.5)
∂qi ∂ q˙i
and the equation of motion can also be written in the form of Hamilton’s equations,
∂H ∂H
ṗi = − and q˙i = . (8.7)
∂qi ∂pi
Finally, the classical system can be quantized by imposing the fundamental postulate
of quantum mechanics,
[qi , pj ] = i~δij . (8.8)
Following the principle of least action we obtain the Euler-Lagrange equation for the
fields:
∂L ∂L
= ∂µ (8.11)
∂φ (x) ∂ (∂µ φ (x))
(If at this point you are confused about the position of Lorentz indices on the right-
hand-side, then remember that what is meant is
∂L ∂ ∂L ∂ ∂L ∂ ∂L ∂ ∂L
= + + + (8.12)
∂φ (x) ∂x0 ∂(∂φ/∂x0 ) ∂x1 ∂(∂φ/∂x1 ) ∂x2 ∂(∂φ/∂x2 ) ∂x3 ∂(∂φ/∂x3 )
You could also use upper indices as long as you are consistent.)
To create a Lorentz covariant theory, we require the Lagrangian to be a Lorentz scalar.
(This also means that in the expression for L above the ‘loose’ Lorentz indices must
somehow be contracted with others.) This requirement imposes certain conditions on
the Lorentz transformation properties of the fields. (We have not discussed these in
detail. See textbooks.) Furthermore, although we consider complex fields, we always
require the Lagrangian to be real.
In quantum field theory, the coordinates φ become operators that obey the standard
quantum mechanical commutation relation with their associated generalized momenta.
The commutation relations lead to creation and annihilation operators. The wave func-
tions that we have considered before can be viewed as single particle excitations that
occur when the creation operators of the field act on the vacuum. For the discussions
here we do not need field theory. What is important to know is that field theory tells us
that, given a Lagrangian, we can find a set of Feynman rules that can be used to draw
diagrams and compute amplitudes.
Now consider the following Lagrangian for a complex scalar field:
You will show in an exercise that the equation of motion corresponding to this La-
grangian is the Klein-Gordon equation. Because the field is complex, it has two separate
components. We could choose these to be the real and imaginary part of the field, such
that φ = φ1 + iφ2 with φ1,2 real. It is easy to see what the Lagrangian looks like and
what the equations of motion become. However, rather than choosing φ1 and φ2 we can
also choose φ and φ∗ to represent the ‘independent’ components of the field.
A similar argument can be made for the bi-spinor and the adjoint bi-spinor in the
Lagrangian of the Dirac field. The latter is given by
L = ψ (iγ µ ∂µ − m) ψ (8.14)
and its equation of motion (treating ψ and ψ as independent components of the field)
is the Dirac equation.
8.4. GLOBAL PHASE INVARIANCE AND NOETHER’S THEOREM 131
φi → φi + i (x) (8.15)
where in the second step we have used the Euler-Lagrange equation to remove ∂L/∂φ.
Consequently, if the Lagrangian is insensitive to the transformation, then the quantity
X ∂L
jµ = i (8.17)
i
∂(∂µ φi )
is a conserved current.
Let’s now apply this to the complex scalar field for a (small) U (1) phase translation.
The two independent field components are φ and φ∗ . Under the phase translation these
change as
φ → φeiα ≈ φ(1 + iα)
(8.18)
φ∗ → φ∗ e−iα ≈ φ∗ (1 − iα)
Consequently, we have φ = iαφ and φ∗ = −iαφ∗ . Inserting these into the expression
for the Noether current, Eq. (8.17), we find
j µ = α i φ(∂ µ φ∗ ) − φ∗ (∂ µ φ)
(8.19)
Since α is an arbitrary constant, we omit it from the current. We have obtained exactly
the current that we constructed for the Klein-Gordon wave in Lecture 1.
If we make the replacement ψ(x) → eiα ψ(x) the expectation value of the observable
remains the same. We say that we cannot measure the absolute phase of the wave
function. (We can only measure relative phases between wave functions in interference
experiments.)
However, this only holds for a phase that is constant in space and time. Are we allowed
to choose a different phase convention on, say, the moon and on earth, for a wave
function ψ(x)? In other words, can we choose a phase that depends on space-time,
In general, we cannot do this without breaking the symmetry. The problem is that the
Lagrangian density L (ψ(x), ∂µ ψ(x)) depends on both the fields ψ(x) and the derivatives
∂µ ψ(x). The derivative term yields:
where Aµ is a new field and q is (for now) an arbitrary constant. Second, we require
that the field Aµ transforms as
1
Aµ (x) → A0µ (x) = Aµ (x) − ∂µ α(x) . (8.24)
q
By inserting the expression for A in the covariant derivative, we find that it just trans-
forms with the local phase α(x):
0 0 iα(x) 1
Dµ ψ(x) → Dµ ψ (x) = e ∂µ ψ(x) + i∂µ α(x)ψ(x) + iqAµ (x)ψ(x) − iq ∂µ α(x)ψ(x)
q
= eiα(x) Dµ ψ(x) (8.25)
As a consequence, terms in the derivative that look like ψ ∗ Dµ ψ are phase invariant.
With the substitution ∂µ → Dµ the Klein-Gordon and Dirac Lagrangians (and any
other real Lagrangian that we can construct with 2nd order terms from a complex field
and its derivatives) satisfy the local phase symmetry.
8.6. APPLICATION TO THE DIRAC LAGRANGIAN 133
L = ψ (iγ µ Dµ − m) ψ
= ψ (iγ µ ∂µ − m) ψ − qAµ ψγ µ ψ (8.26)
≡ Lfree + Lint ,
where we defined the interaction term
Lint = −J µ Aµ (8.27)
Therefore, the requirement of local U (1) invariance automatically implies that the pho-
ton is massless. This actually holds for other gauge symmetries as well. In chapters 11
134 LECTURE 8. LOCAL GAUGE INVARIANCE
to 14 we discuss how masses of vector bosons can be generated in the Higgs mechanism
by ‘breaking’ the symmetry.
You may wonder why we put so much emphasis on the principle of local gauge invariance.
After all, it looks like all we have done is find a different way of arriving at the equations
of motions of electrodynamics: is it really so attractive to formulate QED as a symmetry?
The reason that local gauge symmetries are so important is because of what is called
‘renormalizability’. By way of the Feynman rules, the Lagrangian encodes the infor-
mation to compute scattering and decay processes to arbitrary order. However, if you
compute anything beyond leading order you will quickly find that the result is not finite.
This can be solved by a number of different techniques, called collectively ‘renormaliza-
tion’. It was shown by ’t Hooft and Veldman in the early seventies that only Lagrangians
with interaction terms generated by local gauge symmetries are renormalizable. In other
words, if we want to have a theory in which we can compute something, then we cannot
have any other interactions than those derived from internal symmetries.
Have a careful look at what is written here: The doublet ψ is a 2-component column
vector with a Dirac spinor for each component. Each of the entries in the matrix in the
Lagrangian is again a 4x4 matrix.
Note that we have taken the two components to have identical mass m. Because they
have identical mass and no charge the nucleons are indistinguishable. Therefore, we
8.7. YANG-MILLS THEORY 135
As for the U (1) symmetry, we now try to promote the global symmetry to a local
symmetry. The strategy is similar to that for U (1), but because the group is non-abelian,
the implementation is more complicated. The first step is to make the parameters α
depend on space time. To simplify the notation we define the gauge transformation as
follows,
We have again, as in the case of QED, that the derivative transforms non-trivially
such that the Lagrangian is not phase invariant. To restore phase invariance, we intro-
duce the 2 × 2 covariant derivative
where g is a (so far arbitrary) coupling constant and Bµ a gauge field. In spinor space
the latter is a 2 × 2 unitary matrix with determinant 1. It is customary to parametrize
it in terms of three new real vector fields b1 , b2 and b3 ,
1 1X k k 1 b3 b1 − ib2
Bµ = τ · bµ = τ bµ = . (8.44)
2 2 2 b1 + ib2 −b3
k
We call the fields bi the gauge fields of the SU (2) symmetry. We need three fields rather
than one, because SU (2) has three generators.
In terms of the covariant derivative the Lagrangian is
Dµ0 ψ 0 = ∂µ + igBµ0 ψ 0
Since this expression must hold for all values of the field ψ, we can omit the field ψ from
this expression. If we subsequently multiply both sides of the equation on the right by
G−1 we find for the transformation of the gauge field
i
Bµ0 = GBµ G−1 + (∂µ G) G−1 . (8.50)
g
Although this looks rather complicated we can again try to interpret this by comparing
to the case of electromagnetism. For Gem = eiα(x) we have
i
A0µ = Gem Aµ G−1
em + (∂µ Gem ) G−1
em
q
1
= A µ − ∂µ α (8.51)
q
which is exactly the transformation rule that we had before.
We see that for an SU (2) symmetry the transformation of the gauge field Bµ involves
both a rotation and a gradient. The gradient term was already present in QED. The
rotation term is new. It arises due to the non-commutativity of the elements of SU (2).
If we write out the gauge field transformation formula in the components of the real
vector fields
0 1
bkµ = bkµ − klm αl bm − ∂µ αk (8.52)
g
we can see that there is a coupling between the different components of the field. We call
this the self-coupling. (To derive this start from Eq. 8.50, consider infinitesimal small
α(x) and use the commutation relation of the SU(2) generators, [τi , τj ] = 2ijk τk .)
The effect of the self-coupling becomes clear if one considers the kinetic term of the
SU (2) gauge field. Analogous to the QED case, the three new fields require their own
free Lagrangian, which we write as
1 X µν 1
Lfree
b =− Fl Fµν,l = − F µν · Fµν . (8.53)
4 l 4
Mass terms like m2 bν bν are again excluded by gauge invariance: as for the U (1) symme-
try, the gauge fields must be massless. However, while for the photon the field tensor in
the kinetic term was given by F µν = ∂ µ Aν − ∂ ν Aµ , this form does not work here because
it would break the symmetry. Rather, the individual components of the field tensor are
given by
Flµν = ∂ ν bµl − ∂ µ bνl + g jkl bµj bνk (8.54)
or in vector notation
F µν = ∂ µ bν − ∂ ν bµ − g bµ × bν (8.55)
138 LECTURE 8. LOCAL GAUGE INVARIANCE
As a consequence of the last term the Lagrangian contains contributions with 2, 3 and
4 factors of the b-field. These couplings are respectively referred to as bilinear, trilinear
and quadrilinear couplings. In QED there is only the bilinear photon propagator term.
In the SU (2) theory there are self interactions by a 3-gauge boson vertex and a 4 gauge
boson vertex.
p
Summarizing, we started from the free Lagrangian for a doublet ψ = of two
n
fields with equal mass,
Lfψree = ψ (iγ µ ∂µ − m) ψ
This Lagrangian has a global SU (2) symmetry. We then hypothesized a local SU (2)
phase invariance which we could implement by making the replacement ∂µ → Dµ =
∂µ + igBµ with Bµ = 21 τ · bµ . The full Lagrangian of the theory (which is called the
Yang-Mills theory) is then given by
1
LSU (2) = ψ(iγ µ Dµ − m)ψ − F µν · Fµν
4
1 (8.56)
= ψ(iγ ∂µ − m)ψ − gJ µ bµ − F µν · Fµν
µ
4
≡ Lfree
ψ + L interaction
+ L free
b
where we now absorbed the coupling constant g in the definition of the conserved cur-
rent,
g
J µ = ψγ µ τ ψ (8.57)
2
Comparing this to the QED Lagrangian
1
LU (1) = Lfψree − Aµ · J µ − F µν Fµν (8.58)
4
(with the electromagnetic current J µ = qψγ µ ψ), we see that instead of one field, we now
have three new fields. Furthermore, the kinetic term is more complicated and gives rise
to self-coupling vertices with three and four b-field lines.
copy of the Yang-Mills theory for SU (2) symmetry. The mediators of the force are the
eight massless gluons, corresponding to the 8 generators of the fundamental representa-
tion of SU (3), namely
0 1 0 0 −i 0 1 0 0
λ1 = 1 0 0 λ2 = i 0 0 λ3 = 0 −1 0
0 0 0 0 0 0 0 0 0
0 0 1 0 0 −i 0 0 0
λ4 = 0 0 0 λ5 = 0 0 0 λ6 = 0 0 1
1 0 0 i 0 0 0 1 0
0 0 0 1 0 0
1
λ7 = 0 0 −i λ8 = √ 0 1 0
0 i 0 3 0 0 −2
In this case, too, the Lagrangian contains self-coupling terms for the gauge fields. The
strong interaction is discusses in the Particle Physics II course.
The isospin symmetry in the proton-neutron is a flavour symmetry. Extended to the
system of all other hadrons it is essentially just the symmetry between u and d quarks.
We know that such a symmetry only exists if we ignore electromagnetic interactions,
and the small difference in mass between the u and the d quark (or the proton and the
neutron). Since the symmetry is not exact, we call it an approximate symmetry.
Although the Yang-Mills isospin theory is of no real use to the proton-neutron system, it
turns out to be exactly what is needed to describe the weak interactions. For historical
reasons the local SU (2) symmetry applied to the Lagrangian of Dirac fermion doublets,
discussed in the next lecture, is sometimes called ’weak isospin’. It should certainly not
be confused with the u − d flavour symmetry. In contrast with the flavour symmetry,
gauge symmetries are exact symmetries of the Lagrangian.
But if in addition the scale, or the unit of measure, for f changes by a factor (1 + S µ δxµ )
between x and x + δx, then the value of f becomes:
δf = (∂ µ + S µ ) f δxµ (8.61)
Exercises
Exercise 8.1 (Classical Euler-Lagrange equations)
Use Hamilton’s principle in Eq. (8.3) to derive the Euler-Lagrange equation of motion
(Eq. (8.4)).
Hint: Write down the change in the action δS for an arbitrary small change in the
trajectory (δq, δ q̇). Now use integration by parts to replace δ q̇ by δq. The bound-
ary condition (the fact that we know where the trajectory starts and ends) makes that
δq(t0 ) = δq(t1 ) = 0. The Euler-Lagrange equations then follow from the fact that for the
trajectory that minimizes the action, the change in the action should be zero, independent
of δq.
(b) Show the same for a complex scaler field starting from the Lagrangian
∗ 2 ∗
Lfree µ
KG = (∂µ φ) (∂ φ) − m φ φ (8.64)
taking φ and φ∗ as the (two) independent fields. (Alternatively, you can take the
real and imaginary part of φ. Note that you obtain two equations of motion, one
for φ and one for φ∗ .)
(c) Show that the Euler-Lagrange equations for the Lagrangian
Lfree µ
Dirac = iψγµ ∂ ψ − mψψ (8.65)
leads to the Dirac equations for ψ and for ψ. Note again that you need to consider
ψ and ψ as independent fields.
(d) Show that the Lagrangian
1 1
L = LEM = − (∂ µ Aν − ∂ ν Aµ ) (∂µ Aν − ∂ν Aµ ) − j µ Aµ = − F µν Fµν − j µ Aµ
4 4
(8.66)
leads to the Maxwell equations:
∂µ (∂ µ Aν − ∂ ν Aµ ) = j ν (8.67)
(a) (i) Start with the Lagrange density for a complex Klein-Gordon field
(b) (i) Start with the Lagrange density for a Dirac field
Electroweak Theory
In the previous lecture we have seen how imposing a local gauge symmetry requires a
modification of the free Lagrangian in such a way that a theory with interactions is
obtained. We studied two symmetries, namely
• local U (1) gauge invariance:
ψ (iγ µ Dµ − m) ψ = ψ (iγ µ ∂µ − m) − q ψ γ µ ψ Aµ (9.1)
| {z }
Jµ
For the U (1) symmetry we can identify the Aµ field as the photon. The Feynman rules
for QED, as we discussed them in previous lectures, follow automatically.
Yang and Mills implemented the SU (2) local gauge symmetry, hoping that they could
derive the strong interaction from proton-neutron isospin symmetry. Although that did
not work, we now show that the SU (2) gauge symmetry is still useful, but then to explain
the weak interaction. (The strong interaction follows from SU (3) gauge invariance.)
Before we continue with SU (2), we make a small modification to the interaction terms
above. First consider the U (1) symmetry. Every fermion field has its own charge.
Within the Standard Model we cannot explain why the charge of an up quark is two-
thirds of the charge of an electron. This is why the symbol q appears in the interaction
term above: it is a dimensionless parameter that signifies the strength of the interaction
and it can be different for different fields.
At this point it is customary to introduce a charge operator Q which acts as the generator
of the U (1) symmetry group for electromagnetic interactions. It appears in the field
transformation rule as
ψ 0 = eiα(x)Q ψ . (9.3)
143
144 LECTURE 9. ELECTROWEAK THEORY
where we now use the same coupling gEM for all fields. The fields ψ are eigenstates of
the charge operator Q with an eigenvalue equal to the charge in units of the positron
charge. (Why we do this will become clear later.)
A similar strategy is taken for the isospin symmetry. Rather than τ as the generator
we consider an operator T which for unit isospin charge is given by T = τ /2. It enters
into the douplet transformation rule as
Ψ0 = eiα(x)T Ψ (9.6)
JTµ = Ψ γ µ T Ψ , (9.7)
When, in the following sections, we consider SU (2) symmetry to generate the weak
interaction, the coupling constant g is taken to be the same for all douplets, but the
physical fields (which are eigenstates of T3 ) each have their own value of ‘weak-isospin
charge’, the eigenvalue for T3 . In the Standard Model, this eigenvalue is always ±1/2.
The coupling constants gEM and g in the interaction terms are dimensionless. For
the electromagnetic interaction the coupling is related to the unit charge and the fine-
structure constant as √
e
gEM = √ = 4πα . (9.9)
0 ~c
In our system of units 0 ~c = 1. Therefore, in the following we substitute gEM = e, as
you will find in most textbooks.
We define for any Dirac field ψ the left- and right-handed chiral projections,
1 1
1 − γ5 ψ 1 + γ5 ψ .
ψL ≡ and ψR ≡ (9.10)
2 2
9.1. SU(2) SYMMETRY FOR LEFT-HANDED DOUPLETS 145
As we have seen in Lecture 7, for particles with E m these correspond to the negative
and positive helicity states, respectively. Using the fact that (see exercise 9.1)
ψ γ µ ψ = ψ L γ µ ψL + ψ R γ µ ψR (9.11)
The mass terms ‘mix’ the left- and right-handed components. That is incovenient for
what we are going to do next. Therefore, in the following we consider only massless
fields and deal with non-zero mass later.
Let us now introduce the following doublets for the left-handed chirality states of the
leptons and quarks in the first family:
νL uL
ΨL = and ΨL = (9.14)
eL dL
We call these “weak isospin” doublets. Again, Ψ is not a Dirac spinor, but a doublet of
Dirac spinors. Consider the Lagrangian for the electron and neutrino and verify that it
can be written as (c.f. Eq. 8.35)
L = eR iγ µ ∂µ eR + ν R iγ µ ∂µ ν R +
(9.15)
µ
iγ ∂µ 0 νL
+ (ν L eL )
0 iγ µ ∂µ eL
Now it comes: We impose the SU (2) gauge symmetry on the left-handed doublets only.
That is, we require that the Lagrangian be invariant for local rotations of the doublet.
To do this we need to ignore that the two components of a doublet have different charge,
a problems that we will clearly need to deal with later. As in the Yang-Mills theory, we
also need to ignore that they have different mass, which is another motivation for only
considering massless fields.
The fact that we only impose the gauge symmetry on left-handed states leads to a weak
interaction that is completely left-right asymmetric. This is why it is referred to as
maximal violation of parity.
To construct the weak SU (2)L theory1 we start again with the free Dirac Lagrangian
and we impose SU (2) symmetry on the weak isospin doublets:
Lf ree = ΨL iγ µ ∂µ ΨL (9.16)
1
The subscript L is used to indicate that we only consider SU (2) transformations of the left-handed
doublet.
146 LECTURE 9. ELECTROWEAK THEORY
The generators τ1 and τ2 mix the components of a doublet, while τ3 does not. We define
the fields W ± as
1
Wµ± ≡ √ b1µ ∓ i b2µ
(9.21)
2
The ± index on the W refers to the electric charge. However, at this point we have not
yet shown that these fields are indeed electrically charged: That would require us to
look at the coupling of the W fields to the photon, which we will not do as part of these
lectures. As an alternative, we now show that these W fields couple to charge-lowering
and charge-raising currents. Charge conservation at each Feynman diagram vertex then
implies the charge of the gauge boson.
We define the charged current term of the interaction Lagrangian as
with
τ1 τ2
J 1µ = ΨL γ µ ΨL J 2µ = ΨL γ µ ΨL (9.23)
2 2
2
Note that in terms of physics strong and weak isospin have nothing to do with one another. It is
just that we use the same math!
9.2. THE CHARGED CURRENT 147
As you will show in exercise 9.2 we can rewrite the charged current Lagrangian as
LCC = −g Wµ+ J +µ − g Wµ− J −µ (9.24)
with
1
J µ,± = √ ΨL γ µ τ ± ΨL (9.25)
2
and τ ± 1
= 2 (τ1 ± iτ2 ), or in our representation
+ 0 1 − 0 0
τ = and τ = . (9.26)
0 0 1 0
The leptonic currents can then be written as
1 1
J +µ = √ ν L γ µ eL and J −µ = √ eL γ µ νL (9.27)
2 2
or written out with the left-handed projection operators:
1 1 1
J +µ = √ ν 1 + γ5 γµ 1 − γ5 e
(9.28)
2 2 2
and similar for J −µ . Verify for yourself that
1 + γ 5 γ µ 1 − γ 5 = 2γ µ 1 − γ 5
(9.29)
such that we can rewrite the leptonic charge raising current as
1
J +µ = √ ν γ µ 1 − γ 5 e
(9.30)
2 2
and the leptonic charge lowering current as
1
J −µ = √ e γ µ 1 − γ 5 ν
. (9.31)
2 2
Remembering that a vector interaction has an operator γ µ in the current and an axial
vector interaction a term γ µ γ 5 , we recognize in the charged weak interaction the famous
“V-A” interaction. The story for the quark doublet is identical. Drawn as diagrams,
the charged currents then look as follows:
νe u
Charge raising: W+ W+
e− d
νe u
Charge lowering: W− W−
e− d
148 LECTURE 9. ELECTROWEAK THEORY
The third component of the weak isospin gauge field leads to a neutral current interac-
tion,
Lint = −g b3µ J3µ (9.32)
with b3µ the third gauge boson (another real vector field) and the conserved current given
by
τ3
J3µ = ΨL γ µ ΨL . (9.33)
2
It is now tempting to identify this third component as the Z 0 boson and simply add
the electromagnetic interaction term that we had previously constructed with a U (1)
symmetry with the electromagnetic charge operator Q as generator.
However, this is not a valid way to extend the symmetry of the Lagrangian: the left-
handed doublets that we have constructed are not eigenfunctions of Q since they mix
fields with different charge. Therefore, our SU (2)L invariant Lagrangian cannot be
symmetric under a transformation with Q as generator.
The solution is to start from another U (1) gauge symmetry, called ‘weak hypercharge’.
We denote its generator with the symbol Y and require that it commutes with the
SUL (2) generators. The different members of the isospin multiplet then by construction
obtain the same value of hypercharge.
We denote the combined symmetry by SU (2)L ⊗ U (1)Y . Under this symmetry a left-
handed doublet transform as
where T = τ /2 are the SU (2) generators and Y is the generator for U (1)Y . At the
same time, the right-handed components of the fields in the doublet transform only
under hypercharge,
ΨR → Ψ0R = eiβ(x)Y ΨR . (9.35)
The conserved current corresponding to the U (1)Y symmetry is
JYµ = Ψ γ µ Y Ψ . (9.36)
The Lagrangian following from local SU (2)L ⊗ U (1)Y symmetry takes the form (see e.g.
Halzen and Martin, Chapter 13)
g0 µ
LEW = Lf ree − g JTµ · bµ − J aµ , (9.37)
2 Y
where aµ is the gauge field corresponding to U (1)Y and g 0 /2 is its coupling strength.
The factor 2 appears just because of a convention.
9.3. THE NEUTRAL CURRENT 149
for the τ3 /2 generator.) The right-handed electron is a singlet under SU (2)L and has
T3 = 0. Given a coupling constant e, the observed electromagnetic charge of the elec-
tron is −1. Therefore, the hypercharge of the right-handed electron is −2 while the
hypercharge of the left-handed electron and neutrino are both −1. (The latter two must
be equal, since the SU (2)L doublet is a singlet under U (1)Y .)
Expressing the interactions terms of the bµ3 and aµ fields in the Lagrangian above in
terms of the physical fields, we find
g0 µ JYµ
µ 3 µ 0
−gJ3 bµ − JY aµ = − g sin θW J3 + g cos θW Aµ
2 2
JYµ
µ 0
− g cos θW J3 − g sin θW Zµ
2
µ
≡ −eJEM Aµ − gZ JNµ C Zµ (9.41)
where in the last line we defined the currents and coupling constants associated to the
physical fields.
A direct consequence of Eq. (9.40) is that also the currents are related, namely by
µ 1
JEM = J3µ + JYµ . (9.42)
2
This relation implies that there are only two independent parameters. In particular, the
weak mixing angle is related to the SU (2)L and U (1)Y coupling constants by
g 0 /g = tan θW (9.44)
g0
µ µ µ
− g cos θW J3 − sin θW · 2 (JEM − J3 ) Zµ
2
= ...
e
J3µ − sin2 θW JEM
µ
= − Zµ
cos θW sin θW
Consequently, expressed in J3 and JEM the neutral current is
It is customary to write the term on the right in terms of a vector and an axial vector
coupling such that
µ 1 f
JNC,f = ψ γµ CV − CAf γ 5 ψ (9.51)
2
which implies that
CVf = T3f − 2Qf sin2 θW
(9.52)
CAf = T3f
Alternatively, we can write the current in terms of left- and right-handed fields as
1 f f µ f f
JNµ C,f = CL ψ L γ ψL + CRf ψ R γ µ ψRf (9.53)
2
with the left- and right-handed couplings given by
As stated before the values of the charge of the different fermion fields is not predicted.
Table 9.1 lists the quantum numbers and resulting vector and axial-vector couplings
for all fermions in the Standard Model. The model can be experimentally tested by
measuring these couplings in different processes.
Table 9.1: Gauge interaction quantum numbers and corresponding vector and axial vector
couplings for the fermions in the Standard Model.
Finally, expressed in terms of the left- and right-handed couplings, the Feynman rule
corresponding to Zψψ vertex becomes
f
g 1
Z 0 −i γ µ CVf − CAf γ 5
cos θW 2
f
GF g2
ρ√ = (9.56)
2 8MZ2 cos2 θW
The parameter ρ specifies the relative strength between the charged and neutral current
weak interactions. Comparing the two expressions, we have
2
MW
ρ= (9.57)
Mz2 cos2 θW
The masses of the W ± and Z 0 can be precisely measured, for instance by reconstructing
a ‘two-jet’ invariant mass distribution in high-energy e+ e− collisions: Provided that
9.5. THE MASS OF THE W AND Z BOSONS 153
the collision energy is large enough, the di-jet mass will show mass peaks for ‘on-shell’
produced W ± and Z 0 . The most precise measurement of the four-point coupling for
the charged current comes from the measurement of the muon lifetime. The ratio
of the charged and neutral current couplings was first measured by the Gargamelle
experiment, which exploited an intense neutrino beam to measure the cross-section for
a neutral current process νµ + nucleus → νµ + hadrons and a charged current process
νµ + nucleus → µ + hadrons.
Upon combination of measurements for the couplings an the masses it is found that the
experimental value for ρ is 1 within small uncertainties. This is actually a prediction
of the Higgs mechanism. In the Higgs mechanism the mass generated for the W and Z
are respectively
1 1
q
MW = v g and MZ = v g 2 + g 0 2 , (9.58)
2 2
where v is the so-called vacuum expectation value of the Higgs field. With g 0 /g = tan θW
we find that ρ = 1. Therefore, in the Standard Model the masses of the massive vectors
bosons are related by
MW = MZ cos θW . (9.59)
The best fit of the Standard Model to all experimental data gives approximately
sin2 θW = 0.231 (9.60)
s√
2 e
MW = = 80.4 GeV (9.61)
8GF sin θW
MZ = MW (gz /g) = MW /cos θ = 91.2 GeV (9.62)
Summary
We have introduced a local gauge symmetry SU (2)L ⊗ U (1)Y to obtain a Lagrangian
for electroweak interactions,
g0 µ
µ
− g JL · bµ + JY · aµ (9.63)
2
The coupling constants g and g 0 are free parameters. We can also take e and sin2 θW .
The electromagnetic and neutral weak currents are then given by:
µ 1
JEM = J3µ + JYµ
2
JYµ
JNµ C = J3µ − sin 2 µ
θW JEM 2
= cos θW J3µ 2
− sin θW
2
and the interaction term in the Lagrangian becomes:
µ e µ
− eJEM · Aµ + J · Zµ (9.64)
cos θW sin θW N C
154 LECTURE 9. ELECTROWEAK THEORY
Exercises
Exercise 9.1 (Currents for left and right-handed chirality)
We define the chiral projection operators as PL ≡ 12 (1 − γ 5 ) and PR ≡ 12 (1 + γ 5 ) =
1 − PL and the left- and right-handed chirality bi-spinor states as ψL ≡ PL ψ and
ψR ≡ PR ψ.
(a) Show that the left-handed adjoint spinor, defined as ψ L ≡ (ψL ), is given by
ψ L = ψPR
ψ γ µ ψ = ψ L γ µ ψL + ψ R γ µ ψR
ψ ψ = ψ R ψL + ψ L ψR
157
158 LECTURE 10. THE PROCESS e− e+ → µ− µ+
In the ultra-relativistic limit, we have x → 1 and the correspondence between the chiral
states and helicity states is obtained.
As a consequence, for any process that would violate helicity conservation in the ultra-
relativistic limit, such as the π + → e+ νe decay via the weak interaction, a helicity
suppressing factor (1 − x) appears in the amplitude. Simply said, this is because the
interaction can only couple to a ‘fraction’ (1 − x) of the lepton wave function.
You have shown in the previous lecture that for a vector coupling one can decompose
the current in the vertex factor as follows
ψγ µ ψ = ψ R γ µ ψR + ψ L γ µ ψL (10.5)
This means that a right-handed state only couples to a right-handed state, and a left-
handed state only to a left-handed state. This results holds equally well for an axial
vector coupling (γ µ γ 5 ). It is graphically illustrated in Fig. 10.1. Note that ‘crossing’ a
particle flips its chirality.
R R L L R L
or
L R
Figure 10.1: Helicity conservation in vector and axial-vector couplings. left: A right-handed
incoming electron scatters into a right-handed outgoing electron and vice versa in a vector
or axial vector interaction . right: In the crossed reaction the energy and momentum of one
electron is reversed: i.e. in the e+ e− pair production a right-handed electron and a left-handed
positron (or vice versa) are produced. This is the consequence of a spin=1 force carrier. (In
all diagrams time increases from left to right.)
For scalar couplings the situation is exactly opposite, as its decomposition would read
ψψ = ψ R ψL + ψ L ψR . (10.6)
As we shall see in the remainder of this chapter, conservation of helicity has interesting
consequences for the e− e+ → µ− µ+ process as well. For example, it allows us to under-
stand the angular dependence of a the polarized cross-sections without going in detail
through the kinematics.
µ+
θ −
e+ e
µ−
e+ µ− e+ µ−
Mγ : MZ :
γ Z
− + −
e µ e µ+
Figure 10.2: Leading order Feynman diagrams contributing to e− e+ → µ− µ+ .
In lecture 6 we considered this process in QED. At leading order there was only one
contribution to the amplitude, namely via an intermediate photon. In the electroweak
theory also the amplitude with an intermediate Z 0 boson contributes. The corresponding
Feynman diagrams are shown in Fig. 10.2.
Once we have computed the relevant amplitudes, the differential cross-section follows
as usual from the golden rule,
dσ(e− e+ → µ− µ+ ) 1 1 pf 2
= 2
|M| (10.7)
dΩ 64π s pi
where the invariant amplitude is the sum of the photon and Z 0 contributions.
In Lecture 6 we computed the spin-averaged amplitude via a rather lengthy procedure,
involving Casimir’s track and the trace theorems. Because it is actually a nice illustration
of the concept of helicity conservation, we will here follow a different approach.
Consider first only the matrix element of the photon contribution (evaluated using the
Feynman rules, see e.g. appendix C),
gµν
Mγ = −e2 ψ m γ µ ψm · 2 · ψ e γ ν ψe
(10.8)
q
where the subscript ‘m’ referes to the muon and the subscript ‘e’ to the electron. We
now decompose the spinors in left- and right-handed chirality states, as we did in lecture
160 LECTURE 10. THE PROCESS e− e+ → µ− µ+
11,
ψ m γ µ ψm ψ Lm γ µ ψLm + ψ Rm γ µ ψRm
=
ψ e γµ ψe = ψ Le γµ ψLe + ψ Re γµ ψRe .
e2
ψ Lm γ µ ψLm + ψ Rm γ µ ψrm ·
Mγ = −
s (10.9)
ψ Le γµ ψLe + ψ Re γµ ψRe
dσ unpolarized
1 dσ − + dσ − +
eL eR → µ− +
eL eR → µ− +
= L µR + R µL +
dΩ 4 dΩ dΩ
(10.10)
dσ − + dσ − +
eR eL → µ− + − +
L µR + e e → µR µL
dΩ dΩ R L
where we average over the incoming spins and sum over the final state spins. Note that
e+
R ≡ ψ Le etc.
The amplitudes for each of the spinconfigurations can be computed using the expressions
for the spinors in Eq. 5.83 and 5.84. We start with the outgoing µ− µ+ current, for the
configuration in which the µ− has positive helicity. In terms of the spinors, the transition
current then takes the form
† 0 µ
J µ (µ− + µ
L µR ) = u↑ (pµ− ) γ v↓ (pµ+ ) = [u↑ (pµ− )] γ γ v↓ (pµ+ ) (10.11)
Before we fill in the expressions from Eq. 5.83 and 5.84 we need to worry about the
momentum vectors. We denote the polar angles of the µ− with (θ, φ). Since the µ+
is going exactly in the opposite direction, its polar angles are (π − θ, φ + π). For this
transformation, we have for the various factors that appear in the anti-particle spinors
φ→φ+π
θ →π−θ
sin(θ/2) → cos(θ/2)
cos(θ/2) → sin(θ/2)
eiφ → −eiφ
10.2. THE CROSS-SECTION OF e− e+ → µ− µ+ 161
where the factor E is the product of the two wave function normalization factors and
we have used that u† = (u∗ )T .
The next step is to work out this result for all of the four components of γ µ . For the
zeroeth component we use that (γ 0 )2 = 11 and find J 0 = 0. For the spatial components
we use that γ 0 γ k = αk with the three matrices α given by
0 0 0 1 0 0 0 −i 0 0 1 0
α2 = 0 0 i 0 α3 = 0 0 0 −1 (10.13)
0 0 1 0
α1 =
0 1 0 0 0 −i 0 0 1 0 0 0
1 0 0 0 i 0 0 0 0 −1 0 0
Inserting these in the expression above, gives for the spatial components of the current
J k (µ− +
L µR ) = E 2e−iφ sin2 (θ/2) − 2eiφ cos2 (θ/2),
2ie−iφ sin2 (θ/2) + 2ieiφ cos2 (θ/2), (10.14)
4 cos(θ/2) sin(θ/2)
We can re-use this expression to obtain the electron current for this spin configuration
as well. The ingoing and outgoing particles are reversed, but as shown in Eq.6.22 that
is just a matter of taking the complex conjugate. The electron and position travel along
the z-axis, so θ = 0. The value of φ is arbitrary: it will just lead to an overall phase
which drops out when we take the amplitude squared. We choose φ = π since that
gives the correct correspondence for earlier expressions for spinors with spin along the
z-axis. Inserting these values for the polar angles and taking the complex conjugate to
change the order of the ingoing and outgoing spinors, the electron transition current
then becomes
J µ (e− + µ
L eR ) = v↓ γ u↑ = E 0, 2, 2i, 0 (10.15)
Taking the metric into account, we then obtain for the matrix element
e2 µ ν e2
M = − Je gµν Jm = 8E 2 eiφ cos2 (θ/2) = e2 eiφ (1 + cos θ) (10.16)
s s
where we used that s = 4E 2 . The factor eiφ drops out when we take the square of the
amplitude.
The other three amplitudes can be obtained in a similar way. However, because it
is more instructive, we will show you an alternative method to obtain the angular
162 LECTURE 10. THE PROCESS e− e+ → µ− µ+
dependence, which just uses quantum mechanics. Let us look in more detail at the
helicity dependence (H&M §6.6):
Initial state:
µ+
In the center of mass frame, scattering proceeds from an initial state with JZ = +1 or
−1 along axis ẑ into a final state with JZ0 = +1 or −1 along axis ẑ 0 . Since the interaction
proceeds via a photon with spin J = 1 the amplitude for scattering over an angle θ is
given by the rotation matrices1
where Jy is the y component of the angular momentum operator (which is also the
generator for rotations around the y-axis). The coefficients djm,m0 are sometimes called
‘Wigner d-matrices’. Computing them for the spin-1 system is not so hard (see e.g.
exercise 10.2, or H&M exercise 2.6) and gives
1
d1+1,+1 (θ) = d1−1,−1 (θ) = (1 + cos θ)
2 (10.19)
1
d1+1,−1 (θ) = d1−1,+1 (θ) = (1 − cos θ)
2
1
See H&M§2.2: X
e−iθJ2 |j mi = djm m0 (θ) |j m0 i (10.17)
m0
10.2.2 Z 0 contribution
Having written the total cross-section as a sum of polarized amplitudes we are ready to
include the contribution from the Z 0 boson amplitude. Using the Feynman rules (see
e.g. appendix C) we find for the invariant amplitude
g2 µ m m 5
gµν − qµ qν /MZ2 ν e e 5
MZ = − ψ m γ C V − C A γ ψm · 2
· ψ e γ C V − C A γ ψe
4 cos2 θw q 2 − MZ
(10.23)
√
We can simplify the Z 0 propagator if we ignore the lepton masses (m` s). In that
case the Dirac equation becomes:
ψ e (i∂µ γ µ − m) = 0 ⇒ ψ e (γ µ pµ,e ) = 0 (10.24)
Since pe = 21 q we also have:
1
ψ (γ µ qµ ) = 0 (10.25)
2 e
As a result the qµ qν term in the propagator vanishes and we obtain for the matrix
element
−g 2 1 µ m m 5
e e 5
MZ = 2
· ψ m γ C V − C A γ ψ m ψ e γµ C V − C A γ ψe .
4 cos2 θw q 2 − MZ
(10.26)
We have shown in the previous lecture how the (CVe − CAe γ 5 ) terms can be written in
terms of left- and right-handed couplings. Defining
CR ≡ CV − CA and CL ≡ CV + CA (10.27)
164 LECTURE 10. THE PROCESS e− e+ → µ− µ+
one finds
CV − CA γ 5 ψ = CR ψR + CL ψL .
(10.28)
Comparing this the expression to Eq. (10.9) we realize that we can obtain the polarized
cross-sections directly from the results obtained for the QED process. For two of the
four contributions we then obtain
dσ α2
e− + − +
(1 + cos θ)2 · |1 + r CLm CLe |2
e → µL R =
µ
dΩ γ,Z L R 4s
(10.30)
dσ α2
e− + − +
(1 − cos θ)2 · |1 + r CRm CLe |2
e → µR L =
µ
dΩ γ,Z L R 4s
with the relative contribution of the Z 0 and γ parameterized as
g2 1 s
r= . (10.31)
e 4 cos θw s − Mz2
2 2
The other two helicity configuration follow using the relation in Eq. (10.21) and replacing
CL by CR etc . Using the relation between the coupling constants
GF g2 g2
√ = 2
= . (10.32)
2 8MW 8MZ2 cos2 θw
we can also write r in terms of GF as
√
2GF MZ2 s
r= . (10.33)
e2 s − MZ2
The propagator for the massive vectors bosons has a ‘pole’ at the boson mass: it becomes
infinitely large for an ‘on-shell’ (p2 = m2 ) boson. As you can readily see from the
expression above,
√ this would lead to an infinite cross-section when we tune the beam
energies to s = MZ . The problem is that the propagator does not take into account
the finite decay width of the Z 0 . The Z 0 boson is not a stable particle and hence the ’on-
shell’ Z 0 is actually something with a rather broad mass distribution. We can account
for the width by replacing the mass in the propagator with
i
MZ → MZ − ΓZ (10.34)
2
10.2. THE CROSS-SECTION OF e− e+ → µ− µ+ 165
where Γ is the total decay width of ‘on-shell’ (i.e. not virtual) Z 0 -bosons.
A heuristic explanation (Halzen and Martin, §2.10) is as follows. The decay of an
unstable particle follows the exponential law
where |ψ(0)| is the probability (density) at t = 0 and 1/Γ is the lifetime. Therefore, the
time-dependence√of the wave function, which already involves the rest mass, must also
include a factor e−Γt/2 , or
Consequently, with the substitution above we can ‘correct’ the propagator mass for the
finite decay width. The lineshape that results from such a propagator is usually called
a (spin-1) Breit-Wigner.
With the replacement MZ → MZ −iΓZ /2 in the propagator, the expression for r becomes
√
2GF MZ2 s
r= 2 2
e s − MZ − i Γ2Z
√ (10.37)
2GF MZ2 s
=
Γ2
e2 s − Mz2 − 4Z + iMZ ΓZ
√
Note that r is now complex and that the phase of r depends on s.
Assuming “lepton universality”, we have CVe = CVµ and CAe = CAµ . The expressions for
A0 and A1 then become
2
A0 = 1 + 2 Re(r) CV2 + |r|2 CV2 + CA2
A1 = 4 Re(r) CA2 + 8|r|2 CV2 CA2
In the Standard Model the coefficients for leptons are CA = − 21 and CV = − 21 +2 sin2 θw .
Finally, we consider the total cross-section. In the integrated cross-section the term
proportional to cos θ vanishes and for the other term we use
Z
16π
(1 + cos2 θ) dΩ = (10.40)
3
Consequently, the total cross-section is
4πα2
σ(s) = A0 (s) (10.41)
3s
To summarize, on the amplitude level there are two diagrams that contribute:
e− µ− e− µ−
Mγ : γ MZ :
Z (10.42)
e+ µ+ e+ µ+
dσ
[Z, Z] = Z · Z ∝ |r|2
dΩ
dσ
[γZ] = γ · Z ∝ Re (r)
dΩ
dσ
[γ, γ] = γ · γ ∝1
dΩ
the expression for the differential cross-section is
dσ dσ dσ dσ
= [γ, γ] + [Z, Z] + [γ, Z] (10.43)
dΩ dΩ dΩ dΩ
with
dσ α2
1 + cos2 θ
[γ, γ] =
dΩ 4s
dσ α2 2 h e 2 2 2
i
CV + CAe 2 CVf + CAf 1 + cos2 θ + 8CVe CVf CAe CAf cos θ
[Z, Z] = |r|
dΩ 4s
dσ α2 h i
Re|r| CVe CVf 1 + cos2 θ + 2CAe CAf cos θ
[γ, Z] =
dΩ 4s
10.3. NEAR THE RESONANCE 167
Exactly at the resonance, this gives for the total cross-section to the final state f f¯:
+ − ¯ G2F s0 MZ2 e 2 e2
2
f f2
σ(e e → f f )|s=s0 = CV + CA CV + CA (10.48)
6π Γ2Z
For quark-antiquark final states (f = q) we need to take into account that there are
three distinct colour configurations, namely blue-anti-blue, red-anti-red and green-anti-
green. Therefore, for a quark-anti-quark pair, the cross-section involves another factor
Nc = 3,
+ − G2F s0 MZ2 e 2 e2
q2 q2
σ(e e → q q̄)|s=s0 = Nc · C V + CA CV + CA (10.49)
6π Γ2Z
Figure 10.3 shows the measured cross-section in hadronic final states as function of the
collision energy.
√
At collision energies well above the typical QCD binding energy ( s 2mπ ), the q q̄
state is observed as two ‘jets’, collimated showers of light mesons. The ratio between
the hadronic and leptonic event yields at the Z 0 resonance,
σ(e+ e− → hadrons)
Rl = (10.50)
σ(e+ e− → µ+ µ− )
provides an important test of the standard model, as shown in Fig. 10.4.
168 LECTURE 10. THE PROCESS e− e+ → µ− µ+
σhad [nb]
0
σ
40
ALEPH
DELPHI
L3
OPAL
30
ΓZ
20
measurements, error bars
increased by factor 10
10 σ from fit
QED unfolded
MZ
86 88 90 92 94
Ecm [GeV]
Figure 10.3: left: The Z 0 -lineshape: the cross-section for e+ e− → hadrons as a function of
√
s. right: Same but now near the resonance. The dashed line represents the leading order
computation,while the continuous gray line includes higher order corrections.
0.022
68% CL
0.018
0,l
mt
Afb
αs
0.014 mH
+−
ll
+ −
ee
+ −
µµ
+ −
ττ
0.01
20.6 20.7 20.8 20.9
Rl
combined in plots with SLD results
Figure 10.4: left: Tests of the standard model. The leptonic Af b vs. Rl . The contours show
the measurements while the arrows show the dependency on Standard Model parameters.
right: Determination of the vector and axial vector couplings.
A direct consequence of the photon-Z 0 interference is that the angular distribution is not
symmetric. Figure 10.5 shows the cos θ distribution observed at the Jade experiment,
which operated at the PETRA collider in Hamburg. The beam energy in this experiment
was not yet sufficient to directly produce Z 0 bosons. Still, the effect of the interference
was clearly visible long before the direct discovery of the Z 0 resonance.
At the peak and ignoring the pure photon exchange (because it is negligibly small), the
10.5. THE Z 0 DECAY WIDTH AND THE NUMBER OF LIGHT NEUTRINOS 169
√
Figure 10.5: Angular distribution for e+ e− → µ+ µ− for 25 < s < 37 GeV at the JADE
experiment. θ is the angle between the outgoing µ+ and the incoming e+ . The curves show
fits to the data p(1 + cos2 θ) + q cos θ (full curve) and p(1 + cos2 θ) (dashed curve). (Source:
JADE collaboration, Phys. Lett. Vol108B, p140-144, 1982.)
dσ 8
∝ 1 + cos2 θ + AFB cos θ , (10.51)
d cos θ 3
3 2CVf CAf
A0,f
FB = Ae Af with Af = . (10.52)
4 CV2 + CA2
Using the Feynman rules we can also compute the Z 0 → f f¯ decay width, represented
by the diagram f
f
You cannot easily do this computation yourself, since we have not discussed the external
line for the Z 0 in this course. (The computation needs to take into account the three
170 LECTURE 10. THE PROCESS e− e+ → µ− µ+
polarization states of the massive vector boson.) The result of the computation is
1 1 2
Γ Z → ff = M
16π MZ
g2 Mz f 2 f2
= C V + C A (10.53)
48π cos2 θw
GF MZ3 f 2 2
= √ CV + CAf
6 2 π
For quark-antiquark final states (f = q) we again need to multiply by the colour factor
Nc = 3,
GF MZ3 q 2 q2
Γ(Z → qq) = √ CV + CA · NC . (10.54)
6 2 π
The total decay width of the Z 0 is the sum of all partial widths to all accessible final
states,
where Nν is the number of neutrino species, which is equal to three in the standard
model.
Using all available data to extract information on the couplings we can compute the
decay widths to all final states within the standard model,
1
Γee ≈ Γµµ ≈ Γτ τ = 84 MeV CV ≈ 0 CA = −
2
1 1
Γνν = 167 MeV CV = CA =
2 2
1
Γuu ≈ Γcc = 276 MeV CV ≈ 0.19 CA =
2
1
Γdd ≈ Γss ≈ Γbb = 360 MeV CV ≈ −0.35 CA = −
2
√
A measurement of the lineshape (the cross-section as function of s) gives for the total
decay width of the Z 0 ,
ΓZ ≈ 2490 MeV
So, even though we cannot see the neutrino contribution, we can estimate the number
of neutrinos from the total width of the Z 0 . The result is
ΓZ − 3Γl − Γhad
Nν = = 2.984 ± 0.008 . (10.56)
Γνν
Figure 10.6 shows the predicted lineshape for different values of Nν . This results put
strong constraints on extra generations: if there is a fourth generation, then either it
has a very heavy neutrino, or its neutrino does not couple to the Z 0 . In either case, this
generation would be very different from the known generations of quarks and leptons.
10.5. THE Z 0 DECAY WIDTH AND THE NUMBER OF LIGHT NEUTRINOS 171
Exercises
Exercise 10.1 (Z 0 production and decay)
(a) Derive the expression for Re(κ) in Eq. (10.46).
(b) Calculate the relative contribution of the Z 0 -exchange and the γ exchange to the
cross-section at the Z 0 peak. Use sin2 θW = 0.23, Mz = 91GeV and ΓZ = 2.5GeV .
(c) Show also that at the peak
12π Γe Γµ
σpeak (e− e+ → µ− µ+ ) ≈ (10.57)
Mz2 Γ2Z
(d) Why does the top quark not contribute to the decay width of the Z 0 ?
(e) Calculate the value of Rl = Γhad /Γlep at the resonance s = s0 . Ignore the masses
of the fermions, as we did in the lecture. You may also ignore the contribution of
the photon, as it is very small at the resonance.
√
(f ) The actual line shape of the Z 0 -boson is not a pure Breit Wigner: at the high s
side of the peak the cross-section is higher then expected from the formula derived
in the lectures. Can you think of a reason why this would be the case?
(g) The number of light neutrino generations is determined from the “invisible width”
of the Z 0 -boson as follows:
ΓZ − 3Γl − Γhad
Nν =
Γν
Can you think of another way to determine the decay rate of Z 0 → ν ν̄ directly?
Do you think this method is more precise or less precise?
Hint: Write Sx in terms of the spin ladder operators and use that all states are normal-
ized to 1.
Lecture 11
Symmetry breaking
After a review of the shortcomings of the model of electroweak interactions in the Stan-
dard Model, in this section we study the consequences of spontaneous symmetry break-
ing of (gauge) symmetries. We will do this in three steps of increasing complexity and
focus on the principles of how symmetry breaking can be used to obtain massive gauge
bosons by working out in full detail the breaking of a local U(1) gauge invariant model
(QED) and give the photon a mass.
In the theory of Quantum ElectroDynamics (QED) the requirement of local gauge in-
0
variance, i.e. the invariance of the Lagrangian under the transformation φ → eiα(x) φ
plays a fundamental rôle. Invariance was achieved by replacing the partial derivative by
a covariant derivative, ∂µ → Dµ = ∂µ − ieAµ and the introduction of a new vector field
0
A with very specific transformation properties: Aµ → Aµ + 1e ∂µ α. This Lagrangian for
a free particle then changed to:
1
LQED = Lfree + Lint − Fµν F µν ,
4
which not only ’explained’ the presence of a vector field in nature (the photon), but also
automatically yields an interaction term Lint = eJ µ Aµ between the vector field and the
particle as explained in detail in the lectures on the electroweak model. Under these
173
174 LECTURE 11. SYMMETRY BREAKING
1 2 1 1 1 1
mγ Aµ Aµ = m2γ (Aµ + ∂µ α)(Aµ + ∂ µ α) 6= m2γ Aµ Aµ
2 2 e e 2
The example using only U(1) and the mass of the photon might sounds strange as the
photon is actually massless, but a similar argument holds in the electroweak model for
the W and Z bosons, particles that we know are massive and make the weak force only
present at very small distances.
Just like in QED, invariance under local gauge transformations in the electroweak model
~ µ + ig 0 1 Y Bµ
requires introducing a covariant derivative of the form Dµ = ∂µ + ig 12 ~τ · W 2
introducing a weak current, J weak and a different transformation for isospin singlets and
doublets. A mass term for a fermion in the Lagrangian would be of the form −mf ψ̄ψ,
but such terms in the Lagrangian are not allowed as they are not gauge invariant. This
is clear when we decompose the expression in helicity states:
3] Violating unitarity
To keep the theory renormalizable, we need a very high degree of symmetry (local
gauge invariance) in the model. Dropping the requirement of the local SU(2)L × U(1)Y
gauge invariance is therefore not a wise decision. Fortunately there is a way out of this
situation:
Introduce a new field with a very specific potential that keeps the full Lagrangian
invariant under SU(2)L × U(1)Y , but will make the vacuum not invariant under
this symmetry. We will explore this idea, spontaneous symmetry breaking of a
local gauge invariant theory (or Higgs mechanism), in detail in this section.
The Higgs mechanism: - Solves all the above problems
- Introduces a fundamental scalar → the Higgs boson !
L = T(kinetic) − V(potential)
In general, the Lagrangian for a real scalar particle (φ) is given by:
We can interpret the particle spectrum of the theory when studying the Lagrangian
under small perturbations. In expression (11.1), the constant (potential) term is for
most purposes of no importance as it does not appear in the equation of motion, the
term linear in the field has no direct interpretation (and should not be present as we will
explain later), the quadratic term in the fields represents the mass of the field/particle
and higher order terms describe interaction terms.
To describe the main idea of symmetry breaking we start with a simple model for a real
scalar field φ (or a theory to which we add a new field φ), with a specific potential term:
1
L = (∂µ φ)2 − V(φ)
2
1 1 1
= (∂µ φ)2 − µ2 φ2 − λφ4 (11.2)
2 2 4
As before, to investigate the particle spectrum in the theory, we have to look at small
perturbations around this minimum. To do this it is more natural to introduce a field
η (simply a shift of the φ field) that is centered at the vacuum: η = φ − v.
1 1
Potential term: V(η) = + µ2 (η + v)2 + λ(η + v)4
2 4
1 1
= λv 2 η 2 + λvη 3 + λη 4 − λv 4 ,
4 4
where we used µ2 = −λv 2 from equation (11.3). Although the Lagrangian is still
symmetric in φ, the perturbations around the minimum are not symmetric in η, i.e.
V(−η) 6= V(η). Neglecting the irrelevant 41 λv 4 constant term and neglecting terms or
order η 2 we have as Lagrangian:
1 1 1
Full Lagrangian: L(η) = (∂µ η)(∂ µ η) − λv 2 η 2 − λvη 3 − λη 4 − λv 4
2 4 4
1
= (∂µ η)(∂ µ η) − λv 2 η 2
2
From section 11.2 we see that this describes the kinematics for a massive scalar particle:
1 2 √ p
mη = λv 2 → mη = 2λv 2 = −2µ2 Note: mη > 0.
2
178 LECTURE 11. SYMMETRY BREAKING
At first glance, adding a V (φ) term as in equation (11.2) to the Lagrangian implies
adding a particle with imaginary mass with a four-point self-interaction. However,
when examining the particle spectrum using perturbations around the vacuum, we see
that it actually describes a massive scalar particle (real, positive mass) with three-
and four-point self-interactions. Although the Lagrangian retains its original symmetry
(symmetric in φ), the vacuum is not symmetric in the field η: spontaneous symmetry
breaking. Note that we have added a single degree of freedom to the theory: a scalar
particle.
Note that the Lagrangian is invariant under a U(1) global symmetry, i.e. under φ0 → eiα φ
since φ0∗ φ0 → φ∗ φe−iα e+iα = φ∗ φ.
There are again two distinct cases: µ2 > 0 and µ2 < 0. As in the previous section, we
investigate the particle spectrum by studying the Lagrangian under small perturbations
around the vacuum.
11.4.1 µ2 > 0
V(Φ)
This situation simply describes two massive scalar par-
ticles, each with a mass µ with additional interactions:
1 1 1 1
L(φ1 , φ2 ) = (∂µ φ1 )2 − µ2 φ21 + (∂µ φ2 )2 − µ2 φ22
φ2 |2 {z 2 } |2 {z 2 }
particle φ1 , mass µ particle φ2 , mass µ
φ1
+ interaction terms
11.4. BREAKING A GLOBAL SYMMETRY 179
11.4.2 µ2 < 0
Neglecting the constant and higher order terms, the full Lagrangian can be written as:
1 1
L(η, ξ) = (∂µ η)2 − (λv 2 )η 2 + (∂µ ξ)2 + 0 · ξ 2 + higher order terms
2
| {z } 2
| {z }
massive scalar particle η massless scalar particle ξ
Unlike the η-field, describing radial excitations, there is no ’force’ acting on oscillations
along the ξ-field. This is a direct consequence of the U(1) symmetry of the Lagrangian
and the massless particle ξ is the so-called Goldstone boson.
Goldstone theorem:
For each broken generator of the original symmetry group, i.e. for each generator that
connects the vacuum states one massless spin-zero particle will appear.
In this section we will take the final step and study what happens if we break a local
gauge invariant theory. As promised in the introduction, we will explore its consequences
using a local U(1) gauge invariant theory we know (QED). As we will see, this will allow
to add a mass-term for the gauge boson (the photon).
Local U(1) gauge invariance is the requirement that the Lagrangian is invariant under
0
φ → eiα(x) φ. From the lectures on electroweak theory we know that this can be achieved
by switching to a covariant derivative with a special transformation rule for the vector
field. In QED:
The local U(1) gauge invariant Lagrangian for a complex scalar field is then given by:
1
L = (Dµ φ)† (Dµ φ) − Fµν F µν − V (φ)
4
The term 41 Fµν F µν is the kinetic term for the gauge field (photon) and V (φ) is the extra
term in the Lagrangian we have seen before: V (φ∗ φ) = µ2 (φ∗ φ) + λ(φ∗ φ)2 .
11.5. BREAKING A LOCAL GAUGE INVARIANT SYMMETRY: THE HIGGS MECHANISM181
The situation µ2 > 0: we have a vacuum at 00 . The exact symmetry of the La-
grangian is preserved in the vacuum: we have QED with a massless photon and two
massive scalar particles φ1 and φ2 each with a mass µ.
In the situation µ2 < 0 we have an infinite number of vacua, each satisfying φ21 + φ22 =
−µ2 /λ = v 2 . The particle spectrum is obtained by studying the Lagrangian under small
oscillations using the same procedure as for the continuous global symmetry from section
(11.4.2). Because of local gauge invariance some important differences appear. Extra
terms will appear in the kinetic part of the Lagrangian through the covariant derivatives.
Using again the shifted fields η and ξ we define the vacuum as φ0 = √12 [(v + η) + iξ].
Potential term: V (η, ξ) = λv 2 η 2 , up to second order in the fields. See section 11.4.2.
1 1 1 1
L(η, ξ) = (∂µ η)2 − λv 2 η 2 + (∂µ ξ)2 − Fµν F µν + e2 v 2 A2µ − evAµ (∂ µ ξ) +int.-terms
|2 {z } |2 {z } |4 {z 2 } | {z }
η-particle ξ-particle photon field ?
(11.5)
At first glance: massive η, massless ξ (as before) and also a mass term for the photon.
However, the Lagrangian also contains strange terms that we cannot easily interpret:
−evAµ (∂ µ ξ). This prevents making an easy interpretation.
Looking at the terms involving the ξ-field, we see that we can rewrite them as:
2
1 1 1 2 2 1 1 2 2 0 2
(∂µ ξ)2 − evAµ (∂µ ξ) + e2 v 2 A2µ = e v Aµ − (∂µ ξ) = e v (Aµ )
2 2 2 ev 2
182 LECTURE 11. SYMMETRY BREAKING
This specific choice, i.e. taking α = −ξ/v, is called the unitary gauge. Of course, when
choosing this gauge (phase of rotation α) the field φ changes accordingly (see first part
of section 11.1 and dropping terms of O(ξ 2 , η 2 , ξη) ):
0 1 1 1
φ → e−i ξ/v
φ = e−i ξ/v
√ (v + η + iξ) = e−i ξ/v
√ (v + η)e+i ξ/v
= √ (v + h)
2 2 2
Here we have introduced the real h-field. When writing down the full Lagrangian in
this specific gauge, we will see that all terms involving the ξ-field will disappear and
that the additional degree of freedom will appear as the mass term for the gauge boson
associated to the broken symmetry.
1 1 2 2 2 1 1
= (∂µ h)2 − λv 2 h2 + e v Aµ + e2 vA2µ h + e2 A2µ h2 − λvh3 − λh4
|2 {z } |2 {z } | {z2 } | {z 4 }
massive scalar gauge field (γ) interaction Higgs Higgs self-
particle h with mass and gauge fields interactions
Expanding the terms in the Lagrangian associated to the vector field we see that we do
not only get terms proportional to A2µ , i.e. a mass term for the gauge field (photon), but
also automatically terms that describe the interaction of the Higgs field with the gauge
field. These interactions, related to the mass of the gauge boson, are a consequence of
the Higgs mechanism.
In our model, QED with a massive photon, when expanding 12 e2 A2µ (v + h)2 we get:
11.5. BREAKING A LOCAL GAUGE INVARIANT SYMMETRY: THE HIGGS MECHANISM183
γ
2] e 2
vA2µ h: photon-Higgs three-point interaction h
γ
3] 12 e2 A2µ h2 : photon-Higgs four-point interaction h
h
γ
We added a complex scalar field (2 degrees of freedom) to our existing theory and broke
the original symmetry by using a ’strange’ potential that yielded a large number of
vacua. The additional degrees of freedom appear in the theory as a mass term for the
gauge boson connected to the broken symmetry (mγ ) and a massive scalar particle (mh ).
Exercises
(b) Show that in this model the Higgs boson can decay into two photons and that the
coupling h → γγ is proportional to mγ .
(c) Draw all Feynman vertices that are present in this model and show that Higgs
three-point (self-)coupling, or h → hh, is proportional to mh .
184 LECTURE 11. SYMMETRY BREAKING
(d) Higgs boson properties: how can you see from the Lagrangian that the Higgs boson
is a scalar (spin 0) particle ? What defines the ’charge’ of the Higgs boson ?
Terms ∝ φ6 are allowed since they introduce additional interactions that are not can-
celled by gauge boson interactions, making the model non-renormalizable. Just ignore
this little detail for the moment and compute the ’prediction’ for the Higgs boson mass.
2
(c) Use V (φ† φ) = µ2 φ2 − λφ4 + 34 δφ6 , with µ2 < 0, λ > 0 and δ = − 2λ
µ2
.
q
Show that mh (new) = 32 mh (old), with ’old’: mh for the normal Higgs potential.
Lecture 12
In this section we will apply the idea of spontaneous symmetry breaking from section
11 to the model of electroweak interactions. With a specific choice of parameters we
can obtain massive Z and W bosons while keeping the photon massless.
To break the SU(2)L × U(1)Y symmetry we follow the ingredients of the Higgs mecha-
nism:
1) Add an isospin doublet:
φ+ 1 φ1 + iφ2
φ= =√
φ0 2 φ3 + iφ4
Since we would like the Lagrangian to retain all its symmetries, we can only add
SU(2)L × U(1)Y multiplets. Here we add a left-handed doublet (like the electron
neutrino doublet) with weak Isospin 21 . The electric charges of the upper and lower
component of the doublet are chosen to ensure that the hypercharge Y=+1. This
requirement is vital for reasons that will become more evident later.
2) Add a potential V(φ) for the field that will break (spontaneously) the symmetry:
2
V (φ) = µ2 (φ† φ) + λ(φ† φ) , with µ2 < 0
185
186 LECTURE 12. THE HIGGS MECHANISM IN THE STANDARD MODEL
1 ~ µ + ig 0 1 Y Bµ
Dµ = ∂µ + ig ~τ · W
2 2
3) Choose a vacuum:
We have seen that any choice of the vacuum that breaks a symmetry will generate a
mass for the corresponding gauge boson. The vacuum we choose has φ1 =φ2 =φ4 =0
and φ3 = v:
1 0
Vacuum =φ0 = √
2 v+h
How do we check if the symmetries associated to the gauge bosons are broken ? In-
variance implies that eiαZ φ0 = φ0 , with Z the associated ’rotation’. Under infinitesimal
rotations this means (1 + iαZ)φ0 = φ0 → Zφ0 = 0.
0 1 1 0 1 v+h
SU(2)L : τ1 φ0 = √ = +√ 6= 0 → broken
1 0 2 v+h 2 0
0 −i 1 0 i v+h
τ2 φ0 = √ = −√ 6= 0 → broken
i 0 2 v+h 2 0
1 0 1 0 1 0
τ3 φ0 = √ = −√ 6= 0 → broken
0 −1 2 v+h 2 v+h
1 0 1 0
U(1)Y : Y φ0 = Y φ0 √ = +√ 6= 0 → broken
2 v+h 2 v+h
12.3. SCALAR PART OF THE LAGRANGIAN: GAUGE BOSON MASS TERMS187
This means that all 4 gauge bosons (W1 , W2 , W3 and B) acquire a mass through the
Higgs mechanism. In the lecture on electroweak theory we have seen that the W1 and
W2 fields mix to form the charged W + and W − bosons and that the W3 and B field will
mix to form the neutral Z-boson and photon.
W1 W2 W3 B
| {z } | {z }
W+ and W− bosons Z-boson and γ
When computing the masses of these mixed physical states in the next sections, we
will see that one of these combinations (the photon) remains massless. Looking at the
symmetries we can already predict this is the case. For the photon to remain massless
the U(1)EM symmetry should leave the vacuum invariant. And indeed:
1 1 0 1 0
U(1)EM : Qφ0 = (τ3 + Y )φ0 = √ = 0 → unbroken
2 0 0 2 v+h
It is not so strange that U(1)EM is conserved as the vacuum is neutral and we have:
0
φ0 → eiαQφ0 φ0 = φ0
To obtain the masses for the gauge bosons we will only need to study the scalar part of
the Lagrangian:
Lscalar = (Dµ φ)† (Dµ φ) − V (φ) (12.1)
The V (φ) term will again give the mass term for the Higgs boson and the Higgs self-
interactions. The (Dµ φ)† (Dµ φ) terms:
1 1
~ µ + ig Y Bµ √
0 1 0
Dµ φ = ∂µ + ig ~τ · W
2 2 2 v+h
will give rise to the masses of the gauge bosons (and the interaction of the gauge bosons
with the Higgs boson) since, as we discussed in section 11.5.4, working out the (v + h)2 -
terms from equation (12.1) will give us three terms:
188 LECTURE 12. THE HIGGS MECHANISM IN THE STANDARD MODEL
We can then also easily compute (Dµ φ)† : (Dµ φ)† = − √iv8 g(W1 + iW2 ) , (−gW3 + g Yφ0 Bµ )
0
and we get the following expression for the kinetic part of the Lagrangian:
µ † 1 2h 2 2 2 0 2
i
(D φ) (Dµ φ) = v g (W1 + W2 ) + (−gW3 + g Yφ0 Bµ ) (12.2)
8
Before we can interpret this we need to rewrite this in terms of W+ , W− , Z and γ since
that are the gauge bosons that are observed in nature.
When discussing the charged current interaction on SU(2)L doublets we saw that the
charge raising and lowering operators connecting the members of isospin doublets were
τ+ and τ− , linear combinations of τ1 and τ2 and that each had an associated gauge
boson: the W+ and W− .
+ ν
W
1 0 1
τ+ = (τ1 + iτ2 ) =
2 0 0 e−
1 0 0 − ν
τ− = (τ1 − iτ2 ) = W
2 1 0
e−
12.3. SCALAR PART OF THE LAGRANGIAN: GAUGE BOSON MASS TERMS189
Looking at the terms involving W1 and W2 in the Lagrangian in equation (12.2), we see
that:
2 2
g 2 (W12 + W22 ) = g 2 (W + + W − ) or, alternatively, 2g 2 W + W − (12.3)
g2 −gg 0 Yφ0
0 2 W3
(−gW3 + g Yφ0 Bµ ) = (W3 , Bµ )
−gg 0 Yφ0 g02 Bµ
When looking at this expression there are some important things to note, especially
related to the role of the hypercharge of the vacuum, Yφ0 :
2 If Yφ0 = ± 1, the determinant of the mixing matrix vanishes and one of the
combinations will be massless (the coefficient for that gauge field squared is 0). In
our choice of vacuum we have Yφ0 = +1 (see Exercise 4 why that is a good idea).
In the rest of our discussion we will drop the term Yφ0 and simply use its value of
1.
The two eigenvalues and eigenvectors are given by [see Exercise 3]:
eigenvalue eigenvector
g0
1 1
λ=0 → p = (g 0 W3 + gBµ ) = Aµ photon(γ)
g
p
2 0 2 2 0 2
g +g g +g
2 1 g 1
λ = (g 2 + g 0 ) → p 0 =p (gW3 − g 0 Bµ ) = Zµ Z-boson (Z)
2
g +g 0 2 −g 2
g +g 0 2
2
(−gW3 + g 0 Yφ0 Bµ )2 = (g 2 + g 0 )Zµ2 + 0 · A2µ (12.4)
190 LECTURE 12. THE HIGGS MECHANISM IN THE STANDARD MODEL
Finally, by combining equation (12.3) and (12.4) we can rewrite the Lagrangian from
equation (12.2) in terms of the physical gauge bosons:
1
(Dµ φ)† (Dµ φ) = v 2 [g 2 (W + )2 + g 2 (W − )2 + (g 2 + g 0 )Zµ2 + 0 · A2µ ]
2
(12.5)
8
As a general mass term for a massive gauge boson V has the form 21 MV2 Vµ2 , from equation
(12.5) we see that:
1
MW + = MW − = vg
2 q
1
MZ = v (g 2 + g 0 2 )
2
Although since g and g 0 are free parameters, the SM makes no absolute predictions for
MW and MZ , it has been possible to set a lower limit before the W - and Z-boson were
discovered (see Exercise 2). The measured values are MW = 80.4 GeV and MZ = 91.2
GeV.
Although there is no absolute prediction for the mass of the W- and Z-boson, there is
a clear prediction on the ratio between the two masses. From discussions in QED we
know the photon couples to charge, which allowed us to relate e, g and g 0 (see Exercise
3):
e = g sin(θW ) = g 0 cos(θW ) (12.6)
In this expression θW is the Weinberg angle, often used to describe the mixing of the
W3 and Bµ -fields to form the physical Z boson and photon. From equation (12.6) we
see that g 0 /g = tan(θW ) and therefore:
1
MW vg
= 1 p2 = cos(θW )
MZ v g 2 + g02
2
12.5. MASS OF THE HIGGS BOSON 191
Similar to the Z boson we have now a mass for the photon: 21 Mγ2 = 0, so:
Mγ = 0. (12.7)
Although v is known (v ≈ 246 GeV, see below), since λ is a free parameter, the mass of
the Higgs boson is not predicted in the Standard Model.
νµ
Extra: how do we know v ?: µ GF
s νe
2 G
g GF 1 Fermi: ∝ √F e
Muon decay: 2
= √ →v= √ 2
8MW 2 2GF νµ
µ g
1 −5
We used MW = 2
vg.Given GF = 1.166 · 10 , we νe
g2 W g
see that v = 246 GeV. This energy scale is known EW: ∝ 8MW2 e
Exercises
Exercise 12.1 (Higgs - Vector boson couplings)
In the lecture notes we focussed on the masses of the gauge bosons, i.e. part 1) when
expanding the ((v + h)2 )-terms as discussed in Section 11.5.4 and 12.3. Looking now at
the terms in the Lagrangian that describe the interaction between the gauge fields and
the Higgs field, show that the four vertex factors describing the interaction between the
Higgs boson and gauge bosons: hWW, hhWW, hZZ, hhZZ are given by:
192 LECTURE 12. THE HIGGS MECHANISM IN THE STANDARD MODEL
MV2 MV2
3-point: 2i v
g µν and 4-point: 2i v2
g µν , with (V = W,Z).
Note: A vertex factor is obtained by multiplying the term involving the interacting
fields in the Lagrangian by a factor i and a factor n! for n identical particles in the
vertex.
(c) bonus: Imagine that we would have chosen Yφ0 = −1. What, in that scenario,
0
0 0
would be the (mass-)eigenvectors Aµ and Zµ , the ’photon’ and ’Z-boson’ ? In such
a model, what would be their masses ? Compare them to those in the Standard
Model.
(a) Looking only at the part involving Wµ3 and Bµ show that:
gg 0
Y 1 2 02 Y
Dµ = ∂µ + iAµ p T3 + + iZµ p g T3 − g
g 02 + g 2 2 g 02 + g 2 2
(b) Make also a final interpretation step for the Aµ part and show that:
gg 0 Y
p = e and T3 + = Q, the electric charge.
g 02 + g 2 2
12.5. MASS OF THE HIGGS BOSON 193
(c) bonus: Imagine that we would have chosen Yφ0 = −1. Show explicitly that in
0
that case the photon does not couple to the electric charge.
In this section we discuss how fermions acquire a mass and use our knowledge on the
Higgs coupling to fermions and gauge bosons to predict how the Higgs boson decays as a
function of its mass. Even though the Higgs boson has been discovered, we also discuss
what theoretical information we have on the mass of the Higgs boson as it reveals the
impact on the Higgs boson at higher energy scales (evolution of the universe).
In section 11 we saw that terms like 12 Bµ B µ and mψ̄ψ were not gauge invariant. Since
these terms are not allowed in the Lagrangian, both gauge bosons and fermions are
massless. In the previous section we have seen how the Higgs mechanism can be used
to accommodate massive gauge bosons in our theory while keeping the local gauge in-
variance. As we will now see, the Higgs mechanism can also give fermions a mass: ’twee
vliegen in een klap’.
A term like −mψ̄ψ = −m[ψ¯L ψR + ψ¯R ψL ], i.e. a decomposition in chiral states (see
exercise 1). Such a term in the Lagrangian is not gauge invariant since the left handed
fermions form an isospin doublet (for example νe ) and the right handed fermions
L
form isospin singlets like eR . They transform differently under SU(2)L × U(1)Y .
0 ~ ~
left handed doublet = χL → χL = χL eiW ·T +iαY
0
right handed singlet = ψR → ψR = ψR eiαY
195
196 LECTURE 13. FERMION MASSES, HIGGS DECAY AND LIMITS ON MH
This means that the term is not invariant under all SU(2)L × U(1)Y ’rotations’.
If we could make a term in the Lagrangian that is a singlet under SU(2)L and U(1)Y
, it would remain invariant. This can be done using the complex (Higgs) doublet we
introduced in the previous section. It can be shown that the Higgs has exactly the right
quantum numbers to form an SU (2)L and U (1)Y singlet in the vertex: −λf ψ̄L φψR ,
where λf is a so-called Yukawa coupling.
When we write out this term we’ll see that this does not only describe an interaction
between the Higgs field and fermion, but that the fermions will acquire a finite mass if
the φ-doublet has a non-zero expectation value. This is the case as φ0 = √12 v +0 h as
before.
1 0 ν
Le = −λe √ (ν̄, ē)L eR + ēR (0, v + h)
2 v + h e L
λe (v + h)
= − √ [ēL eR + ēR eL ]
2
λe (v + h)
= − √ ēe
2
λe v λ
= − √ ēe − √e hēe
2 2
| {z } | {z }
electron mass term electron-higgs interaction
λe v λ
me = √ √e ∝ me
2 2
A few side-remarks:
13.1. FERMION MASSES 197
√ m
1) The Yukawa coupling is often expressed as λf = 2 vf and the coupling of the
λ m
fermion to the Higgs field is √f2 = vf , so proportional to the mass of the fermion.
2) The mass of the electron is not predicted since λe is a free parameter. In that
sense the Higgs mechanism does not say anything about the electron mass itself.
The fermion mass term Ldown = λf ψ̄L φψR (leaving out the hermitian conjugate term
ψ̄R φ̄ψL for clarity) only gives mass to ’down’ type fermions, i.e. only to one of the isospin
doublet components. To give the neutrino a mass and give mass to the ’up’ type quarks
(u, c, t), we need another term in the Lagrangian. Luckily it is possible to compose a
new term in the Lagrangian, using again the complex (Higgs) doublet in combination
with the fermion fields, that is gauge invariant under SU(2)L × U(1)Y and gives a mass
to the up-type quarks. The mass-term for the up-type fermions takes the form:
c ∗ 1 (v + h)
φ̃ = −iτ2 φ = − √ (13.2)
2 0
As we will discuss now, this is not the whole story. If we look more closely we’ll see that
we can construct more fermion-mass-type terms in the Lagrangian that cannot easily
be interpreted. Getting rid of these terms is at the origin of quark mixing.
198 LECTURE 13. FERMION MASSES, HIGGS DECAY AND LIMITS ON MH
This section will discuss in full detail the consequences of all possible allowed quark
’mass-like’ terms and study the link between the Yukawa couplings and quark mixing in
the Standard Model: the difference between mass eigenstates and flavour eigenstates.
If we focus on the part of the SM Lagrangian that describes the dynamics of spinor
(fermion) fields ψ, the kinetic terms, we see that:
We saw that using the Higgs field φ we could construct terms in the Lagrangian of the
form given in equation (13.1). For up and down type fermions (leaving out the hermitian
conjugate term) that would allow us to write for example:
where the strength of the interactions between the Higgs and the fermions, the so-called
Yukawa couplings, had again to be added by hand.
13.2. YUKAWA COUPLINGS AND THE ORIGIN OF QUARK MIXING 199
This looks straightforward, but there is an additional complication when you realize
that in the most general realization the Λ’s are matrices. This will introduce mixing
between different flavours as we will see a little bit later. In the most general case, again
leaving out the h.c., the expression for the fermion masses is written as:
where the last term is the mass term for the charged leptons. The matrices Yijd , Yiju and
Yijl are arbitrary complex matrices that connect the flavour eigenstate since also terms
like Yuc will appear. These terms have no easy interpretation:
v v v
− LYukawa = ... + −Λdd √ u¯I uI −Λus √ u¯I sI −Λss √ s¯I sI + . . . (13.4)
2 2 2
| {z } | {z } | {z }
mass−term down quark ?? mass−term strange quark
To interpret the fields in the theory as physical particles, the fields in our model should
have a well-defined mass. This is not the case in equation (13.4). If we write out all
Yukawa terms in the Lagrangian we realize that it is possible to re-write them in terms
of mixed fields that do have a well-defined mass. These states are the physical particles
in the theory
Since this is the crucial part of flavour physics, we spell out the term Yijd QILi φ dIRj
explicitly and forget about the other 2 terms in expression (13.3):
φ+
Yijd QILi φ dIRj = Yijd (up-type down-type)IiL (down-type)IRj =
φ
+ + +
I φ I φ I φ
Y11 (u d)L 0 Y12 (u d)L 0 Y13 (u d)L 0
φ+ φ+ φ+
I
dR
φ φ φ
Y21 (c s)IL 0 Y 22 (c s) I
L 0 Y 23 (c s) I
L 0
· sIR
φ φ φ
+ + + bIR
φ φ φ
Y31 (t b)IL Y32 (t b)IL Y33 (t b)IL
φ0 φ0 φ0
After symmetry breaking we get the following mass terms for the fermion fields:
−Lquarks d I I u I
Yukawa = Yij QLi φ dRj + Yij QLi φ̃ uRj
I
v v
= Yijd dILi √ dIRj + Yiju uILi √ uIRj + ...
2 2
d I I u I I
= Mij dLi dRj + Mij uLi uRj + , (13.5)
200 LECTURE 13. FERMION MASSES, HIGGS DECAY AND LIMITS ON MH
where we omitted the corresponding interaction terms of the fermion fields to the Higgs
field, q̄qh(x) and the hermitian conjugate terms. Note that the d’s and u’s in equation
(13.5) still each represent the three down-type and up-type quarks respectively, so the
’mixed’-terms are still there. To obtain mass eigenstates, i.e. states with proper mass
terms, we should diagonalize the matrices M d and M u . We do this with unitary matrices
V d as follows:
d
Mdiag = VLd M d VRd†
u
Mdiag = VLu M d VRu†
Using the requirement that the matrices V are unitary (VLd† VLd = 1) and leaving out
again the hermitian conjugate terms the Lagrangian can now be expressed as follows:
−Lquarks I d I I u I
Yukawa = dLi Mij dRj + uLi Mij uRj + ...
= dILi VLd† VLd Mijd VRd† VRd dIRj + uILi VLu† VLu Miju VRu† VRu uIRj + ...
= dLi (Mijd )diag dRj + uLi (Miju )diag uRj + ...,
where in the last line the matrices V have been absorbed in the quark states. Note that
the up-type and down-type fields are now no longer the interaction states uI and dI , but
are now ’simply’ u and d. A bit more explicit, we now have the following quark mass
eigenstates:
The interaction terms are obtained by imposing gauge invariance by replacing the partial
derivative by the covariant derivate
W-bosons:
i
= iQILi γµ ∂ µ + gWiµ τi QILi
Lkinetic, weak (QL )
2
u I
I µ i µ
= i(u d)iL γµ ∂ + gWi τi
2 d iL
g g
= iuIiL γµ ∂ µ uIiL + idIiL γµ ∂ µ dIiL − √ uIiL γµ W −µ dIiL − √ dIiL γµ W +µ uIiL + ...
2 2
If we now express the Lagrangian in terms of the quark mass eigenstates d, u instead
of the weak interaction eigenstates dI , uI , the ’price’ to pay is that the quark mixing
between families (i.e. the off-diagonal elements) appear in the charged current inter-
action as each of the interaction fields is now replaced by a combination of the mass
eigenstates:
g g
Lkinetic, cc (QL ) = √ uIiL γµ W −µ dIiL + √ dIiL γµ W +µ uIiL + ...
2 2
g g
= √ uiL (VLu VL )ij γµ W −µ diL + √ diL (VLd VLu† )ij γµ W +µ uiL + ...
d†
2 2
The combination of matrices (VLd VLu† )ij , a unitary 3×3 matrix is known under the short-
hand notation VCKM , the famous Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix.
By convention, the interaction eigenstates and the mass eigenstates are chosen to be
equal for the up-type quarks, whereas the down-type quarks are chosen to be rotated,
going from the interaction basis to the mass basis:
uIi = uj
dIi = VCKM dj
or explicitly:
dI Vud Vus Vub d
sI = Vcd Vcs Vcb s (13.7)
bI Vtd Vts Vtb b
We should note here that in principle a similar matrix exists that connects the lep-
ton flavour and mass eigenstates. In this case, contrary to the quarks, the down-type
interaction doublet-states (charged leptons) are chosen to be the same as the mass eigen-
states. The rotation between mass and interaction eigenstates is in the neutrino sector.
This matrix is known as the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix and
has a completely different structure than the one for quarks. Just like for the CMS
matrix, the origin of the observed patterns are completely unknown. A last thing to
remember: neutrino interaction eigenstates are known as νe , νµ and ντ , whereas the
physical particles, the mass eigenstates, are ν1 , ν2 and ν3 .
Now that we have derived the coupling of fermions and gauge bosons to the Higgs field,
we can look in more detail at the decay of the Higgs boson.
The general expression for the two-body decay rate:
dΓ |M|2
= |pf | S, (13.8)
dΩ 32π 2 s
with M the matrix element, |pf | the momentum of the√produced particles and√S = n!1
for n identical particles. In a two-body decay we have s = mh and |pf | = 12 β s (see
exercise 2). Since the Higgs boson is a scalar particle, the Matrix element takes a simple
form:
imf
−iM = ū(p1 ) v(−p2 ) h
v
−imf
iM† = v̄(−p2 ) u(p1 )
v
13.3. HIGGS BOSON DECAY 203
Since there are no polarizations for the scalar Higgs boson, computing the Matrix ele-
ment squared is ’easy’:
m 2 X
2 f
M = (v̄)s2 (−p2 )us1 (p1 )(ū)s1 (p1 )vs2 (−p2)
v s1 ,s2
m 2 X
f
X
= us1 (p1 )(ū)s1 (p1 ) v̄s2 (−p2 )vs2 (−p2 )
v s1 s1
m 2
f
= Tr ((6 p1 + mf )(− 6 p2 − mf ))
v
m 2
− Tr(6 p1 6 p2 ) − m2f Tr(1))
f
=
v
m 2
f
−4p1 · p2 − 4m2f
=
v
use: s = (p1 − p2 )2 = p21 + p22 − 2p1 · p2 and since p21 = p22 = m2f
and s = m2h we have m2h = 2m2f − 2p1 · p2
m 2
f
2m2h − 8m2f
=
v s
m 2
f 2 2
4m2f
= 2mh β , with β = 1 − 2
v mh
Decay rate:
√ √
Starting from equation (13.8) and using M2 (above), |pf | = 21 β s, S=1 and s = mh
we get:
dΓ |M|2 Nc mh mf 2 3
= |p f | S = β
dΩ 32π 2 s 32π 2 v
R
Doing the angular integration dΩ = 4π we finally end up with:
Nc 2
Γ(h → f f¯) = m mh βf3 .
8πv 2 f
The decay ratio to gauge bosons is a bit more tricky, but is explained in great detail in
Exercise 5.
204 LECTURE 13. FERMION MASSES, HIGGS DECAY AND LIMITS ON MH
h √ 4m2f
Γ(h → f f¯) = Nc
8πv 2
m2f mh 1−x , with x = m2h
g2
√
h Γ(h → V V ) = 64πMW2 m3h SV V (1 − x + 34 x2 ) 1 − x
4MV2
, with x = m2h
and SW W,ZZ = 1, 12 .
The decay of the Higgs boson to two off-shell gauge bosons is given by:
3M 4 0
Γ(h → V V ∗ ) = 32π2Vv4 mh δV R(x) , with
h 0 0
δW = 1, δZ = 12 − 9 sin θW + 27 sin4 θW
7 10 2 40
, with
3(1−8x+20x2 ) 3x−1
1−x
R(x) = √
4x−1
acos 2x3/2 − 2x (2 − 13x + 47x2 )
− 23 (1 − 6x + 4x2 ) ln(x)
Since the coupling of the Higgs boson to gauge bosons is so much larger than that to
fermions, the Higgs boson decays to off-shell gauge bosons even though MV ∗ + MV <
2MV . The increase in coupling ’wins’ from the Breit-Wigner suppression. For example:
at mh = 140 GeV, the h → W W ∗ is already larger than h → bb̄.
γ 2
α2 4 (f )
m3h 2
P
h Γ(h → γγ) = 256π 3 v 2 3 f Nc ef − 7
γ
γ , where ef is the fermion’s electromagnetic charge.
h Note: - WW contribution ≈ 5 times top contribution
γ
- Some computation also gives h → γZ
2
αs2
95 7Nf αs
Γ(h → gluons) = m3 1+ − + ...
72π 3 v 2 h 4 6 π
h
Note: - The QCD higher order terms are large.
- Reading the diagram from right to left you see the dominant
production mechanism of the Higgs boson at the LHC.
13.4. THEORETICAL BOUNDS ON THE MASS OF THE HIGGS BOSON 205
2
10
τ τ
+ -
and Section 13.3.2 we can compute the rel- 10
--
ative branching fractions for the decay of a cc
Although the Higgs mass is not predicted within the minimal SM, there are theoretical
upper and lower bounds on the mass of the Higgs boson if we assume there is no new
physics between the electroweak scale and some higher scale called Λ. In this section
we present a quick sketch of the various arguments and present the obtained limits.
As the Higgs boson mass is now known to quite some precision this section might
feel strange and unnecessary to revisit. Since similar arguments are used to obtain
theoretical limits on the mass of hypothetical particles that are predicted in models
that go beyond the Standard Model it is good to understand the various elements that
enter in such a discussion.
13.4.1 Unitarity
In the absence of a scalar field the amplitude for elastic scattering of longitudinally
polarised massive gauge bosons (e.g. WL+ WL− → WL+ WL− ) diverges quadratically with
the centre-of-mass energy when calculated in perturbation theory and at an energy of
1.2 TeV this process violates unitarity. In the Standard Model, the Higgs boson plays
an important role in the cancellation of these high-energy divergences. Once diagrams
involving a scalar particle (the Higgs boson) are introduced in the gauge boson scattering
mentioned above, these divergences are no longer present and the theory remains unitary
and renormalizable. Focusing on solving these divergences alone also yields most of the
Higgs bosons properties. This cancellation only works however if the Higgs boson is not
too heavy. By requiring that perturbation theory remains valid an upper limit on the
Higgs mass can be extracted. With the requirement of unitarity and using all (coupled)
206 LECTURE 13. FERMION MASSES, HIGGS DECAY AND LIMITS ON MH
This number comes from an analysis that uses a partial wave decomposition for the
matrix element M, i.e.:
l=∞
dσ 1 X
= M2 , with M = 16π (2l + 1)Pl (cos θ)al ,
dΩ 64πs l=0
where Pl are Legendre polynomials and al are spin-l partial waves. Since (WL+ WL− +
ZL + ZL + HH)2 is well behaved, it must respect unitarity, i.e. |ai | < 1 or |Re(ai )| ≤ 0.5.
As the largest amplitude is given by:
GF m2h 3
amax
0 =− √ ·
4π 2 2
This limit is soft, i.e. it means that for Higgs boson masses > 700 GeV perturbation
theory breaks down.
In this section, the running of the Higgs self-coupling λ with the renormalisation scale
µ is used to put both a theoretical upper and a lower limit on the mass of the Higgs
boson as a function of the energy scale Λ.
Similar to the gauge coupling constants, the coupling λ ’runs’ with energy.
dλ
= βλ , where t = ln(Q2 ).
dt
13.4. THEORETICAL BOUNDS ON THE MASS OF THE HIGGS BOSON 207
Although these evolution functions (called β-functions) have been calculated for all SM
couplings up to two loops, to focus on the physics, we sketch the arguments to obtain
these mass limits by using only the one-loop results. At one-loop the quartic coupling
runs with the renormalisation scale as:
dλ 3 2 1 2 1 4 0
≡ βλ = λ + λht − ht + B(g, g ) (13.9)
dt 4π 2 2 4
, where ht is the top-Higgs Yukawa coupling as given in equation (13.1). The dominant
terms in the expression are the terms involving the Higgs self-coupling λ and the top
quark Yukawa coupling ht . The contribution from the gauge bosons is small and explic-
itly given by B(g, g 0 ) = − 18 λ(3g 2 + g 02 ) + 64
1
(3g 4 + 2g 2 g 02 + g 04 ). The terms involving the
mass of the Higgs boson, top quark and gauge bosons can be understood from looking
in more detail at the effective coupling at higher energy scales, where contributions from
higher order diagrams enter:
= + + + + ...
This expression allows to evaluate the value of λ(Λ) relative to the coupling at a reference
scale which is taken to be λ(v).
If we study the β-function in 2 special regimes: λ g, g 0 , ht or λ g, g 0 , ht , we’ll see
that we can set both a lower and an upper limit on the mass of the Higgs boson as a
function of the energy-scale cut-off in our theory (Λ):
For large values of λ (heavy Higgs boson since m2h = 2λv 2 ) and neglecting the effects
from gauge interactions and the top quark, the evolution of λ is given by the dominant
term in equation (13.9) that can be easily solved for λ(Λ):
dλ 3 2 λ(v)
= λ ⇒ λ(Λ) = 3λ(v)
(13.10)
dt 4π 2 Λ2
1− 4π 2
ln v2
Note:
208 LECTURE 13. FERMION MASSES, HIGGS DECAY AND LIMITS ON MH
Λ2
3λ(v) 2 /3λ(v)
ln = 1 → At a scale Λ = ve2π λ(Λ) is infinite.
4π 2 v2
downwards’, i.e. find λ(v) for which λ(Λ) = ∞ (the Landau pole) we find:
s
4π 2 8π 2 v 2
λmax (v) = 2 ⇒ mh < 2 (13.11)
3 ln Λv2 3 ln Λv2
For Λ=1016 GeV the upper limit on the Higgs mass is 160 GeV/c2 . This limit gets less
restrictive as Λ decreases. The upper limit on the Higgs mass as a function of Λ from a
computation that uses the two-loop β function and takes into account the contributions
from top-quark and gauge couplings is shown in the Figure at the end of Section 13.4.4.
For small λ (light Higgs boson since m2h = 2λv 2 ), a lower limit on the Higgs mass is
found by the requirement that the minimum of the potential be lower than that of the
unbroken theory and that the electroweak vacuum is stable. In equation (13.9) it is
clear that for small λ the dominant contribution comes from the top quark through the
Yukawa coupling (−h4t ).
1 4 3 4 2 02 2
βλ = −3ht + (2g + (g + g ) )
16π 2 16
3 4 4 4
= 2MW + MZ − 4mt
16π 2 v 4
< 0.
13.4. THEORETICAL BOUNDS ON THE MASS OF THE HIGGS BOSON 209
The requirement that λ remains positive up to a scale Λ, such that the Higgs vac-
uum is the global minimum below some cut-off scale, puts a lower limit on λ(v) and
therefore on the Higgs mass:
Λ2
dλ
= βλ → λ(Λ) − λ(v) = βλ ln and require λ(Λ) > 0.
dt v2
Λ2
λ(v) > βλ ln and λmin (v) → (mmin 2
h ) > 2λ
min
(v)v 2 , so
v2
2
Λ
m2h > 2
2v βλ ln
v2
3
(mmin 2 4
+ MZ4 − 4m4t
h ) = 2 2
2MW
8π v
Λ2
> −493 ln
v2
Note: This result makes no sense, but is meant to describe the logic. If we go to the
2-loop beta-function we get a new limit: mh > 130 − 140 GeV if Λ = 1019 GeV. A
detailed evaluation taking into account these considerations has been performed. The
region of excluded Higgs masses as a function of the scale Λ from this analysis is also
shown in the Figure at the end of Section 13.4.4 by the lower excluded region.
210 LECTURE 13. FERMION MASSES, HIGGS DECAY AND LIMITS ON MH
800
M H (GeV/c 2 )
In the Figure on the right the theoretically
mt = 175 GeV/c2
allowed range of Higgs masses is shown as a 600
function of Λ.
400
For a small window of Higgs masses around
Landau pole
160 GeV/c2 the Standard Model is valid up
to the Planck scale (∼ 1019 GeV). For other 200 Allowed
values of the Higgs mass the Standard Model
Vacuum instability
is only an effective theory at low energy and
0
new physics has to set in at some scale Λ. 10
3 6
10 10
9 12
10 10
15
10
18
Λ (GeV)
The electroweak gauge sector of the SM is described by only three independent pa-
rameters: g, g 0 and v. The predictions for electroweak observables, are often presented
using three (related) variables that are known to high precision: GF , MZ and αQED .
To obtain predictions to a precision better than the experimental uncertainties (often
at the per mill level) higher order loop corrections have to be computed. These higher
order radiative corrections contain, among others, contributions from the mass of the
top quark and the Higgs boson. Via the precision measurements one is sensitive to these
small contributions and thereby to the masses of these particles.
Radiative corrections
Apart from the mass of the W -boson, there are more measurements that provide sensi-
tivity to the mass of the Higgs boson. A summary of the measurements of several SM
measurements is given in the left plot of Figures 13.1.
While the corrections connected to the top quark behave as m2t , the sensitivity to the
mass of the Higgs boson is unfortunately only logarithmic (∼ ln mh ):
The results from a global fit to the electroweak data with only the Higgs mass as a free
parameter is shown in the right plot of Figure 13.1. The plot shows the ∆χ2 distribution
as a function of mh . The green band indicates the remaining theoretical uncertainty
in the fit. The result of the fit suggested a rather light Higgs boson and it could be
10
∆ χ2
AUG 11
MZ G fitter
Tevatron 95% CL
SM
LEP 95% CL
0.1
AUG 11
ΓZ G fitter SM
0.1
9 3σ
σ0had -1.7
0
Rlep -1.0
0,l
8
AFB -0.9
Al(LEP) 0.2
7
Al(SLD) -2.0
sin2Θlept (Q ) -0.7
eff FB
0,c
6
AFB 0.9
0,b
AFB 2.5
Ac
5
-0.1
Ab
2σ
0.6
0
4
Rc 0.1
0
Rb -0.8
3 Theory uncertainty
(5)
∆αhad(M2) -0.1
Z
Fit including theory errors
MW -1.3
ΓW
2 Fit excluding theory errors
0.2
1σ
mc -0.0
1
mb -0.0
mt 0.3
0
-3 -2 -1 0 1 2 3 50 100 150 200 250 300
(O - Omeas) / σmeas
fit
MH [GeV]
Figure 13.1: Status of various SM measurements (left) and the ∆χ2 distribution as a function
of mh from a global fit with only mh as a free parameter (right). Before the discovery.
summarised by the central value with its one standard deviation and the one-sided (95%
CL) upper limit:
mh = 95+30
−24
+74
−43 GeV/c2 and mh < 162 GeV/c2 (at 95% CL).
In July 2012 the ATLAS and CMS experiments at the Large Hadron Collider at CERN
announced the discovery of the Higgs boson. We will discuss the details of the search
for the Higgs boson and its discovery in a separate lecture, but we since we cannot have
212 LECTURE 13. FERMION MASSES, HIGGS DECAY AND LIMITS ON MH
a lecture note on the Higgs boson without proof of its discovery I include here 4 plots
that were in the discovery paper of the ATLAS experiment.
Events/5 GeV
Data ATLAS
(*)
25 Background ZZ
Σ weights / 2 GeV
(*)
ATLAS Data S/B Weighted H→ZZ →4l
100 Background Z+jets, tt
Sig+Bkg Fit (mH=126.5 GeV)
Bkg (4th order polynomial)
Signal (m =125 GeV)
80 20 H
Syst.Unc.
60
1 0σ (*)
10-1 1σ
H → WW → lν lν
s = 7 TeV: ∫ Ldt = 4.7 fb-1
-2 2σ
10 s = 8 TeV: ∫ Ldt = 5.8 fb-1
10-3 3σ
10-4
H → γγ
s = 7 TeV: ∫ Ldt = 4.8 fb
-1
4σ
10-5 s = 8 TeV: ∫ Ldt = 5.9 fb
-1
10-6 (*)
10-7
5σ H → ZZ → 4l
s = 7 TeV: ∫ Ldt = 4.8 fb-1
10-8 s = 8 TeV: ∫ Ldt = 5.8 fb-1
10-9 6σ
10-10
Combined
10-11 s = 7 TeV: ∫ Ldt = 4.6 - 4.8 fb
-1
µ = 1.4 ± 0.3
110 115 120 125 130 135 140 145 150 s = 8 TeV: ∫ Ldt = 5.8 - 5.9 fb
-1
mH [GeV]
-1 0 1
Signal strength (µ)
Figure 13.2: Plots from the Higgs discovery paper from ATLAS. Two-photon invariant mass
distribution (top left), the 4-lepton invariant mass distribution (top right), the p-value as a
function of the Higgs mass (bottom left) and the measurement of the coupling strength of the
Higgs boson to gauge bosons and fermions (bottom right).
All results on the Higgs boson from the ATLAS and CMS experiments at the LHC can
be found on these locations:
ATLAS: https://twiki.cern.ch/twiki/bin/view/AtlasPublic/HiggsPublicResults
CMS: http://cms.web.cern.ch/org/cms-higgs-results
13.5. EXPERIMENTAL LIMITS ON THE MASS OF THE HIGGS BOSON 213
Exercises
Exercise 13.1
Show that ūu = (ūL uR + ūR uL )
Exercise 13.2
Show that in a two body decay (a heavy particle M decaying into two particles with
mass m) the momentum of the decay particles can be written as:
√
s √ 4m2
|pf | = β, with β = 1 − x and x =
2 M2
Exercise 13.3
Higgs decay into fermions for mh = 100 GeV.
Use mb = 4.5 GeV, mτ = 1.8 GeV, mc = 1.25 GeV
(a) Compute Γ(H→ bb̄).
(b) Compute Γ(H→ all) assuming only decay into the three heaviest fermions.
(c) What is the lifetime of the Higgs boson. Compare it to that of the Z boson.
(c) Show that the matrix element can finally be written as:
g2 3 2 4M2V
M2 = m4
(1 − x + x ), with x =
4M2W h 4 m2h
(d) Show that the Higgs decay into vector bosons can be written as:
g 2 SV V 3 3 2 √
Γ(h → V V ) = m (1 − x + x ) 1 − x,
64πM2W h 4
4MV2
with x = m2h
and SWW,ZZ = 1, 21 .
Although the Higgs mechanism cures many of the problems in the Standard Model,
there are also several ’problems’ associated to the Higgs mechanism. We will explore
these problems in this section and very briefly discuss the properties of non-SM Higgs
bosons.
Since the Higgs field occupies all of space, the non-zero vacuum expectation value of
the Higgs field (v) will contribute to the vacuum energy, i.e. it will contribute to the
cosmological constant in Einstein’s equations: Λ = 8πG
c4
N
ρvac .
Note that we cannot simply redefine Vmin to be 0, or any arbitrary number since quantum
corrections will always yield a value like the one (order of magnitude) given above. The
Higgs mass is unknown, but since we have a lower limit on the (Standard Model) Higgs
215
216LECTURE 14. PROBLEMS WITH THE HIGGS MECHANISM AND HIGGS SEARCHES
boson mass from direct searches at LEP (mh > 114.4 GeV/c2 ) we can compute the
contribution of the Higgs field to ρvac .
1 2 2
ρHiggs
vac = m v
8 h
1
> 1 · 108 GeV4 and since GeV =
r
> 1 · 108 GeV/r3 (energy density)
Ωm ≈ 30% and ΩΛ ≈ 70% ∼ 10−46 GeV4 → empty space is really quite empty.
In the electroweak theory of the SM, loop corrections are small. In the loops the inte-
gration is done over momenta up to a cut-off value Λ.
mh = mbare
h +∆mh
ferm.
+∆mgauge
h +∆mHiggs
h +... The corrections from the fermions (mainly
from the top quark) are large. Expressed in
terms of the loop-momentum cut-off Λ given
h
t
h h
W/Z
h h
h
h
by:
and +
2 top 3
= − 2 λ2t Λ2
t h
∆m h
W/Z
8π
14.2. HIGGS BOSONS IN MODELS BEYOND THE SM (SUSY) 217
Most popular theoretical solution to the hierarchy problem is the concept of Supersym-
metry, where for every fermion/boson there is a boson/fermion as partner. For example,
the top and stop (supersymmetric bosonic partner of the top quark) contributions (al-
most) cancel. The quadratic divergences have disappeared and we are left with
2 2 2 Λ
∆mh ∝ (mf − mS ) ln .
mS
SM: Add 4 degrees of freedom → 3 massive gauge bosons → 1 Higgs boson (h)
SUSY: Add 8 degrees of freedom → 3 massive gauge bosons → 5 Higgs boson (h, H, A, H+ , H− )
v2
parameters:tan(β) = v1
and MA .
218LECTURE 14. PROBLEMS WITH THE HIGGS MECHANISM AND HIGGS SEARCHES
With the new parameters, all couplings to gauge bosons and fermions change:
SUSY SM
ghV V = ghV V sin(β − α)
Exercises
Exercise 14.1 (b-tagging at LEP)
A Higgs boson of 100 GeV decays at LEP: given a lifetime of a B mesons of roughly 1.6
picoseconds, what distance does it travel in the detector before decaying ? What is the
most likely decay distance ?
(a) Why is there a ’dip’ in te fraction of Higgs bosons that decays to 2 Z bosons
(between 160 and 180 GeV)?
(b) How many events H → ZZ → e+ e− µ+ µ− muons are produced in 1 fb−1 of data for
mh = 140, 160, 180 and 200 GeV ? The expected number of evets is the product
of the luminosity and the cross-section: N = L · σ
On the LHC slides, one of the LHC experiments shows its expectation for an analysis
aimed at trying to find the Higgs boson in the channel with 2 electrons and 2 muons.
We concentrate on mh =140 GeV.
14.2. HIGGS BOSONS IN MODELS BEYOND THE SM (SUSY) 219
(c) What is the fraction of events in which all 4 leptons have been well reconstructed
in the detector ? What is the single (high-energy) lepton detection efficiency ?
Name reasons why not all leptons are detected.
We do a counting experiment using the two bins around the expected Higgs boson mass
(we assume for the moment that the background is extremely well known and does not
fluctuate). In a counting experiment a Poisson distribution describes the probabilities
to observe x events when λ are expected:
λx e−λ
P(x|λ) =
x!
(d) Does this experiment expect to be able to discover the mh =140 GeV hypothesis
after 9.3 fb−1 .
(e) Imagine the data points was the actual measurement after 9.3 fb−1 . Can this
experiment claim to have discovered the Higgs boson at mh =140 GeV?
220LECTURE 14. PROBLEMS WITH THE HIGGS MECHANISM AND HIGGS SEARCHES
Appendix A
This appendix lists some properties of the operators αi and β in the Dirac Hamiltonian:
∂
~
Eψ = i ψ = −i~
α · ∇ + βm ψ
∂t
1. αi and β are hermitian.
They have real eigenvalues because the operators E and p~ are hermitian. (Think
µ
of a plane wave equation: ψ = N e−ipµ x .)
2. T r(αi ) = T r(β) = 0.
Since αi β = −βαi , we have also: αi β 2 = −βαi β. Since β 2 = 1, this implies:
αi = −βαi β and therefore T r(αi ) = −T r(βαi β) = −T r(αi β 2 ) = −T r(αi ), where
we used that T r(A · B) = T r(B · A).
3. The eigenvalues of αi and β are ±1.
To find the eigenvalues bring αi , β to diagonal form and since (αi )2 = 1, the square
of the diagonal elements are 1. Therefore the eigenvalues are ±1. The same is
true for β.
4. The dimension of αi and β matrices is even.
The T r(αi ) = 0. Make αi diagonal with a unitary rotation: U αi U −1 . Then, using
again T r(AB) = T r(BA), we find: T r(U αi U −1 ) = T r(αi U −1 U ) = T r(αi ). Since
U αi U −1 has only +1 and −1 on the diagonal (see 3.) we have: T r(U αi U −1 ) =
j(+1) + (n − j)(−1) = 0. Therefore j = n − j or n = 2j. In other words: n is
even.
221
222 APPENDIX A. SOME PROPERTIES OF DIRAC MATRICES αI AND β
Appendix B
223
224 APPENDIX B. PROPAGATORS AND VERTEX FACTORS
in the limit m = 0. (The propagator for spin-1 particles with non-zero mass is called
the Proca propagator and we will see it in Lecture 9.)
Interactions in the Lagrangian density take the form of terms that have fields of different
types. You cannot just write down arbitrary terms: the requirement that they be
Lorentz invariant restricts their form. For instance, the interaction term for the photon
with the spin-0 field that we have seen above, takes the form
where g is the coupling constant (e for the EM interaction) and j µ is the current for the
complex scalar field
jµ ≡ i [φ∗ (∂µ φ) − (∂µ φ)∗ φ] . (B.6)
It is important that φ and φ∗ are different fields here: in our discussion above one would
be the incoming field and the other the outgoing field. The interaction term then has
three fields and leads to a vertex with three lines. Labeling φ by A and φ∗ by C, we
now express the current in momentum space,
µ
jAC = i [φ∗C (pµA φA ) + (pµC φ∗C ) φA ] (B.7)
(Note how the plus sign appears: it is because pµ ↔ i~∂ µ .) The rule to obtain the vertex
factor is now to omit all the fields from the interaction term and write down what is
left, multiplied by i,
vertex factor = ig(pµA + pµC ) (B.8)
This summarized the recipe to get the vertex factors and propagators. It directly leads
to the ‘Feynman rules’ for a particular Lagrangian. Using the rules we can now compute
the amplitude M for any Feynman diagram.
Appendix C
Arrange the left-handed projections of the lepton and quark fields in doublets
νL uL
ΨL = or ΨL = (C.2)
eL dL
Ignore their masses (or choose them equal within the doublet). Now consider that the
Lagrangian remains invariant under
U (1)Y :
ψ → ψ 0 = eiY β(x) ψ (C.3)
SU (2)L :
ΨL → Ψ0L = exp [iY α
~ (x) · ~τ ] ΨL (C.4)
To keep the Lagrangian invariant compensating gauge fields must be introduced. These
transform simultaneously with the Dirac spinors in the doublet:
U (1)Y : hypercharge field aµ
Y
∂µ → Dµ = ∂µ + ig 0 aµ (C.5)
2
SU (2)L : weak isospin fields b1µ , b2µ , b3µ (only couple to left-handed doublet):
Ignoring the kinetic and self-coupling terms of the gauge fields, the Lagrangian becomes
g0
L = Lfree − i JYµ aµ − ig J~Lµ · ~bµ (C.7)
2
225
226 APPENDIX C. SUMMARY OF ELECTROWEAK THEORY
For the generators of SU (2) we choose the Pauli spin matrices. The first field in a
left-handed doublet has T3 = +1/2 and the second field T3 = −1/2. By construction
the right-handed projections are singlets under SU (3)L and therefore have T3 = 0.
The physical gauge fields (connecting the particle fields) become
“charged currents”
b1µ ∓ ib2µ
Wµ± = √ (C.8)
2
“neutral currents”
Zµ = −aµ sin θw + b3µ cos θw
(C.9)
Aµ = aµ cos θw + b3µ sin θw
The Higgs mechanism takes care that 3 out of 4 gauge bosons get mass. For the field
Aµ (the photon) to be massless, we need
g0
tan θw = (C.10)
g
The coupling of the massless field becomes proportional to a charge
Q = T3 + 12 Y (C.11)
MW = MZ cos θ (C.12)
with
e = g sin θw gz = g/ cos θw tan θw = g 0 /g (C.14)
and
CVf = T3f − 2Qf sin2 θw CAf = T3f (C.15)
The relevant quantum numbers for our fields are
227
f Q TL3 TR3
u, c, t + 32 + 12 0
d, s, b − 13 − 21 0
νe , νµ , ντ 0 + 12 –
e− , µ− , τ − −1 − 21 0
Till now we have ignored that the weak interaction mixes the quark fields. Inserting the
CKM matrix we get for the charged currents,
g
LC.C. = − i √ Vud ψ u γ µ 12 1 − γ 5 ψd Wµ+
2
g (C.16)
∗
ψ d γ µ 12 1 − γ 5 ψu Wµ−
− i √ Vud
2
228 APPENDIX C. SUMMARY OF ELECTROWEAK THEORY
The Feynman rules for the vertex factors are then as follows
f f
f f
ν d
g
W+ −i √ γ µ 12 (1 − γ 5 ) W+ −i √g2 Vud γ µ 12 (1 − γ 5 )
2
`− u