Experiment: CP Violation

Contents
  1. Historical context and the prediction under test
  2. Apparatus
  3. Procedure
  4. Observations and data
  5. Interpretation
  6. Direct CP violation in the kaon system
  7. CP violation beyond the kaons
  8. Cosmological relevance and open questions
  9. Primary references

Tests Equation (105.46), Equation (105.47) and Phenomenon 103.31. Assuming Theorems 105.34, 105.36 and 105.38.

After parity fell in 1957 (Experiment: Parity Violation) the combined operation \(CP\) was expected to be exact, and the neutral kaon system was built around that expectation: the long-lived state was supposed to be \(CP\)-odd and therefore forbidden to decay to two pions. This chapter records the two-pion decay that Christenson, Cronin, Fitch and Turlay found at the Brookhaven alternating-gradient synchrotron in a branching ratio of about \(2\times 10^{-3}\) [Christenson:1964], and everything measured since: direct \(CP\) violation in kaon decay [Fanti:1999] [AlaviHarati:1999], the large asymmetries of the \(B\) mesons [Aubert:2001] [Abe:2001], and the small one in charm [Aaij:2019]. It sits after the weak-interaction and flavour chapters (Weak Interactions and Flavour Physics and Neutrinos) because it is the measurement that forced a third quark generation [Kobayashi:1973], and beside Discrete Symmetries and CPT, whose \(CPT\) theorem is what makes \(CP\) violation and \(T\) violation the same statement.

The division of labour with those two chapters is worth stating at the outset, because it decides what is and is not repeated here. The two-state formalism of the neutral kaon, the Wigner–Weisskopf elimination of the decay continuum, the definition of \(\epsilon\) and \(\epsilon'\), the double ratio, the Jarlskog invariant and the Sakharov conditions are all derived in Sections 105.3 and 105.6; the flavour structure that ties them to the quark mixing matrix is Section 103.5 and Section 103.4.3. What follows is the experimental side: what the beams and the detectors physically were, how each number was extracted, what background had to be excluded, and what the numbers are with their uncertainties in SI units. Where a result of the theory chapters is needed it is cited by its equation number rather than re-derived, and where the experimental argument requires a derivation of its own — the purity of the beam, the regeneration background, the geometry of the spectrometer, the extraction of \(\sin2\beta\) from a vertex separation — it is carried out here.

Its closing sections take the result out of particle physics and into cosmology: \(CP\) violation is one of Sakharov's three conditions for a baryon asymmetry [Sakharov:1967], and the honest arithmetic is that the Cabibbo–Kobayashi–Maskawa phase is many orders of magnitude too small to produce the asymmetry observed in the sky [Gavela:1994] — which is stated here as an open problem, not papered over.

Remark 111.1 (Units in this chapter, and in its sources).

Every equation below carries \(\hbar\) and \(c\) explicitly and every number is quoted in SI, with the electronvolt — an accepted non-SI unit — used for energies and its joule equivalent given where the comparison matters. The primary literature of this subject does not do this: it works in units in which \(\hbar\) and \(c\) are set to unity, so that a mass, a momentum and an energy are all quoted in \(\mathrm{MeV}\) and a lifetime in \(/\mathrm{MeV}\). To map a formula from those papers onto one here, restore one factor of \(c\) for every unit of momentum written as an energy, one factor of \(c^{2}\) for every mass written as an energy, and one factor of \(\hbar\) for every inverse energy written as a time. Two conversions recur and are worth fixing once: \(1\,\mathrm{MeV}=1.602176634\times 10^{-13}\,\mathrm{J}\) exactly, since the elementary charge is exact in the present SI [Mohr:2025], and \(\hbar=6.582119569\times 10^{-16}\,\mathrm{eV}\,\mathrm{s}\), so a width of \(1\,\mathrm{MeV}\) is a lifetime of \(6.58\times 10^{-22}\,\mathrm{s}\). No derivation in this chapter is carried out in those units.

Historical context and the prediction under test

The neutral kaon system

The experiment of 1964 was possible only because an unusual physical system had been identified nine years earlier, and the identification was made on paper before any of it was seen.

Gell-Mann and Pais asked what happens to a neutral particle that is not its own antiparticle when the interaction that produces it conserves a quantum number that the interaction destroying it does not [GellMann:1955]. The strong interaction makes the neutral kaon in a state of definite strangeness — \(\ket{K^{0}}\) with \(S=+1\), \(\ket{\bar{K}^{0}}\) with \(S=-1\) — and conserves that label, so the two are distinct particles. The weak interaction does not conserve it, and therefore connects them. The states of definite mass and lifetime are then not the strangeness eigenstates at all but the combinations that diagonalize the weak effective Hamiltonian, which under an exact \(CP\) symmetry are the \(CP\) eigenstates \(\ket{K_{1}}\) and \(\ket{K_{2}}\) of Equation (105.46).

The consequence Gell-Mann and Pais drew was a prediction of a new particle. Two pions have \(CP=+1\) and three neutral pions \(CP=-1\) (Equation (105.47)), so if \(CP\) is exact the \(CP\)-even combination may decay to \(\pi\pi\) and the \(CP\)-odd one may not. Since three pions barely fit inside the kaon mass, the second must be long-lived. Numerically, with the masses of Table 111.5, the energy released into \(3\pi^{0}\) is

\begin{equation}\tag{111.1} Q_{3\pi^{0}}=m_{K^{0}}c^{2}-3m_{\pi^{0}}c^{2} =497.611\,\mathrm{MeV}-404.930\,\mathrm{MeV} =92.68\,\mathrm{MeV}\ec \end{equation}

against \(Q_{\pi\pi}=218.47\,\mathrm{MeV}\) for \(\pi^{+}\pi^{-}\). The non-relativistic phase-space volume for \(n\) bodies sharing a kinetic energy release \(Q\) grows as \(Q^{(3n-5)/2}\): the momentum delta function leaves \(3(n-1)\) independent components, the energy constraint restricts them to a sphere of radius proportional to \(\sqrt{Q}\) in that space, and the surface measure of that sphere divided by the gradient of the energy supplies \(R^{3(n-1)-1}/R\propto Q^{(3n-5)/2}\). Two bodies therefore give \(Q^{1/2}\) and three bodies \(Q^{2}\), so the three-body channel is far more strongly throttled near its own threshold, and it has less than half the energy release to work with. That is what produces the observed factor of several hundred in lifetime. Lande, Booth, Impeduglia, Lederman and Chinowsky found the long-lived partner at the Brookhaven Cosmotron the following year [Lande:1956], with a lifetime some six hundred times that of the known short-lived \(V\) particle. The modern values,

\begin{equation}\tag{111.2} \tau_{S}=8.954(4)\times 10^{-11}\,\mathrm{s}\ec\qquad \tau_{L}=5.116(21)\times 10^{-8}\,\mathrm{s}\ec \end{equation}

give a ratio \(\tau_{L}/\tau_{S}=571\) [Navas:2024].

The second consequence is the one that makes the system an instrument rather than a curiosity. Because the two mass eigenstates are superpositions of the two strangeness eigenstates, a beam prepared as \(K^{0}\) evolves into a mixture whose strangeness content oscillates, and the oscillation frequency is the mass difference divided by \(\hbar\). This is Equation (105.61). The measured splitting is

\begin{equation}\tag{111.3} \Delta m_{K}c^{2}=3.484(6)\times 10^{-6}\,\mathrm{eV} =5.582\times 10^{-25}\,\mathrm{J}\ec\qquad \frac{\Delta m_{K}c^{2}}{\hbar}=5.293\times 10^{9}\,/\mathrm{s}\ec \end{equation}

corresponding to a mass difference of \(\Delta m_{K}=6.211\times 10^{-42}\,\mathrm{kg}\), which is \(\Delta m_{K}/m_{K}=7.0\times10^{-15}\) of the kaon's own mass (Equations (105.60) and (105.63)). No spectrometer resolves a fractional mass difference of \(10^{-14}\); an interferometer whose fringe is a phase accumulated over a flight path of metres does, and that is what the neutral kaon is. Everything in this chapter rests on that fact: the quantities measured are interference terms, not energies.

Remark 111.2 (Why the long-lived state was expected to be pure).

It is worth being precise about the strength of the pre-1964 expectation, because it is what makes the measurement decisive rather than merely interesting. If \(CP\) commutes with the full Hamiltonian then \(\ket{K_{2}}\) is an exact eigenstate of it and the amplitude for \(K_{2}\to\pi\pi\) is exactly zero — not small, not suppressed by a coupling, but forbidden by a superselection of the same kind that forbids a \(0^{+}\to0^{-}\) electromagnetic transition. There was no parameter to adjust. A single observed event of the right kind, above background and above every mundane explanation, falsifies the hypothesis; and that is what Section 111.4.1 reports.

Regeneration

The mundane explanation that had to be excluded is regeneration, and because it is simultaneously the cleanest demonstration of quantum superposition in the subject and the principal background of the 1964 experiment, it is derived here rather than assumed.

Pais and Piccioni observed that a long-lived beam traversing matter emerges with a short-lived component [Pais:1955]. The reason is that the strong interaction sees strangeness and not \(CP\): a \(\bar{K}^{0}\) can convert a nucleon into a hyperon, a \(K^{0}\) cannot, so the two forward scattering amplitudes differ. A medium that attenuates or phase-shifts the two strangeness components differently is, in the \(CP\) basis, a medium that rotates \(\ket{K_{L}}\) into \(\ket{K_{S}}\).

Proposition 111.3 (Coherent regeneration in a thin slab).

Let a \(K_{L}\) beam of momentum \(p\) traverse a slab of thickness \(\ell\) containing \(N\) scattering centres per unit volume, with complex forward scattering amplitudes \(f\) for \(K^{0}\) and \(\bar{f}\) for \(\bar{K}^{0}\), both carrying the SI dimension of length. If \(\ell\) is short compared with the \(K_{S}\) decay length, the emerging state contains a coherent \(K_{S}\) amplitude

\begin{equation}\tag{111.4} \rho=\ii\tan\frac{\Delta\varphi}{2}\approx\frac{\ii\Delta\varphi}{2} =\ii\pi N\left(f-\bar{f}\right)\ell\,\frac{\hbar}{p}\ec \end{equation}

relative to the surviving \(K_{L}\) amplitude. It vanishes when the medium cannot distinguish the two strangeness states. Rests on Equations (105.46) and (105.64).

Derivation. Derives Proposition 111.3. A wave of reduced wavelength \(\hbar/p\) propagating through a medium of \(N\) centres per unit volume with forward amplitude \(f\) acquires, per unit length, the extra phase \(2\pi Nf\hbar/p\): this is the optical relation \(k^{2}(n-1)=2\pi Nf\) with \(k=p/\hbar\), so that the phase per unit length is \(k(n-1)=2\pi Nf/k=2\pi Nf\hbar/p\), and it is the only piece of scattering theory needed. Both sides of the optical relation carry the SI dimension \(/\mathrm{m}^{2}\), since \([N]=/\mathrm{m}^{3}\) and \([f]=\mathrm{m}\). Over the slab the two strangeness components therefore acquire

\[ \varphi=2\pi N f\ell\,\frac{\hbar}{p}\ec\qquad \bar{\varphi}=2\pi N\bar{f}\ell\,\frac{\hbar}{p}\ec\qquad \Delta\varphi=\varphi-\bar{\varphi}\ep \]

Write the incident state, neglecting the \(10^{-3}\) admixture \(\epsilon\), as \(\ket{K_{L}}=\left(\ket{K^{0}}-\ket{\bar{K}^{0}}\right)/\sqrt{2}\) from Equation (105.46). Applying the two phases and factoring out the common part,

\begin{align*} \ket{\psi}&=\frac{1}{\sqrt{2}} \left(\ee^{\ii\varphi}\ket{K^{0}} -\ee^{\ii\bar{\varphi}}\ket{\bar{K}^{0}}\right)\\ &=\frac{\ee^{\ii(\varphi+\bar{\varphi})/2}}{\sqrt{2}} \left(\ee^{\ii\Delta\varphi/2}\ket{K^{0}} -\ee^{-\ii\Delta\varphi/2}\ket{\bar{K}^{0}}\right)\\ &=\ee^{\ii(\varphi+\bar{\varphi})/2} \left[\cos\frac{\Delta\varphi}{2}\ket{K_{L}} +\ii\sin\frac{\Delta\varphi}{2}\ket{K_{S}}\right]\ec \end{align*}

where the last line uses \(\ket{K_{S}}=\left(\ket{K^{0}}+\ket{\bar{K}^{0}}\right)/\sqrt{2}\). The ratio of the two amplitudes is Equation (111.4), and it is proportional to \(f-\bar{f}\), hence zero for a medium blind to strangeness. Retaining the \(K_{S}\) decay and the \(K_{S}\)–\(K_{L}\) phase difference over the slab replaces the linear \(\ell\) by a bounded oscillating factor of the kind written in Equation (105.64); for the thin, dilute media of interest here that factor reduces to \(\ell\) and the difference is immaterial. Only \(\abs{\rho}\) is used below, so the overall phase convention of the thick-slab formula does not enter.

Two features of Equation (111.4) decide the design of the experiment. The first is that \(\rho\) is proportional to \(N\), so a dilute medium regenerates weakly and a vacuum not at all. The second, less obvious, is that the regenerated amplitude is coherent with the incident beam: it is forward, it carries the beam's momentum direction, and its decay to \(\pi^{+}\pi^{-}\) is therefore collinear with the beam — exactly like the signal. Coherent regeneration cannot be removed by an angular cut. Only its incoherent relatives can: diffractive regeneration off a single nucleus, and inelastic production, both leave the reconstructed parent momentum at a measurable angle to the beam, and those are what the collinearity requirement of Section 111.3 eliminates. The coherent part has to be made small, and made small by construction; that is what the helium bag is for, and Proposition 111.5 puts a number on it.

CP as the symmetry that was supposed to survive

The parity experiments of 1957 (Experiment: Parity Violation) showed that the weak interaction distinguishes left from right maximally: in the beta decay of polarized cobalt the electrons come out preferentially opposite to the nuclear spin [Wu:1957], and in the pion–muon–electron chain the muon is fully polarized [Garwin:1957]. Charge conjugation fares no better, and for the same reason — applied to the left-handed neutrino it produces a left-handed antineutrino, which does not exist.

Landau's response was to propose that the true symmetry is the product [Landau:1957]. It is an economical suggestion: \(CP\) maps the left-handed neutrino onto the right-handed antineutrino, which does exist, so maximal violation of \(C\) and of \(P\) separately is exactly compatible with an exact \(CP\). Together with the \(CPT\) theorem of local Lorentz-invariant quantum field theory [Luders:1957] — proved as Theorem 105.34 — the proposal has three testable consequences, and it is worth listing them because the experiment tested the first and the other two survived:

  1. \(K_{L}\to\pi\pi\) is strictly forbidden, the long-lived state being the \(CP\)-odd combination and the two-pion state \(CP\)-even.

  2. A particle and its antiparticle have exactly equal masses and exactly equal total lifetimes (Theorem 105.36, Theorem 105.38), whether or not \(CP\) holds; these follow from \(CPT\) alone.

  3. Any \(CP\) violation is, given \(CPT\), a \(T\) violation of the same size, so a direct comparison of a process with its motion reverse must show it too — which is the programme of Section 111.7.4.

The observable and what a measurement must deliver

The prediction under test is item (i), and stating it as an experimental requirement fixes the whole design.

A neutral kaon beam is prepared far enough downstream of its production target that the \(K_{S}\) component has decayed to an utterly negligible level. Any \(\pi^{+}\pi^{-}\) pair observed to emerge from the surviving beam must then satisfy two independent conditions if it is a genuine two-body decay of a beam particle: the invariant mass of the pair, computed under the pion hypothesis, must equal the kaon mass; and the vector sum of the two pion momenta — the reconstructed parent momentum — must be collinear with the beam. A measurement must therefore deliver

Each of the five is met by a specific piece of the apparatus, and the next two sections take them in that order.

Apparatus

The Brookhaven beam and the decay volume

The source was the alternating-gradient synchrotron at Brookhaven, delivering protons of \(30\,\mathrm{GeV}\) — a kinetic energy of \(4.8\times 10^{-9}\,\mathrm{J}\) per proton — onto an internal target. A neutral beam was taken off at \(30\,^\circ\) to the circulating protons, passed through a collimator, swept clear of charged particles by a magnet, and re-collimated; the detector viewed the beam \(17.4\,\mathrm{m}\) (57 feet in the units of the original report) from the target [Christenson:1964]. The momentum spectrum of the surviving neutral kaons peaked near \(1.1\,\mathrm{GeV}/c\), which is the value used for every kinematic estimate below.

That distance is the first design requirement, and it is met with an enormous margin.

Proposition 111.4 (The beam is pure long-lived kaons).

At \(p=1.1\,\mathrm{GeV}/c\) the \(K_{S}\) decay length is \(5.93\times 10^{-2}\,\mathrm{m}\) and the \(K_{L}\) decay length is \(33.9\,\mathrm{m}\). A flight path of \(17.4\,\mathrm{m}\) is \(293\) \(K_{S}\) decay lengths and \(0.513\) \(K_{L}\) decay lengths. Rests on Equation (111.2).

Derivation. Derives Proposition 111.4. With \(m_{K}c^{2}=497.611\,\mathrm{MeV}\) and \(pc=1.1\,\mathrm{GeV}\), the total energy is

\[ E=\sqrt{p^{2}c^{2}+m_{K}^{2}c^{4}} =1.2073\,\mathrm{GeV}=1.9343\times 10^{-10}\,\mathrm{J}\ec \]

so \(\gamma=E/(m_{K}c^{2})=2.426\), \(\beta=pc/E=0.9111\) and \(\gamma\beta=pc/(m_{K}c^{2})=2.211\). The laboratory decay length of a state of proper lifetime \(\tau\) is \(\gamma\beta c\tau\); with Equation (111.2),

\begin{align*} \lambda_{S}&=2.211\times2.99792458\times 10^{8}\,\mathrm{m}/\mathrm{s} \times8.954\times 10^{-11}\,\mathrm{s}=5.934\times 10^{-2}\,\mathrm{m}\ec\\ \lambda_{L}&=2.211\times2.99792458\times 10^{8}\,\mathrm{m}/\mathrm{s} \times5.116\times 10^{-8}\,\mathrm{s}=33.90\,\mathrm{m}\ep \end{align*}

Hence \(17.4\,\mathrm{m}/\lambda_{S}=293\) and \(17.4\,\mathrm{m}/\lambda_{L}=0.513\). The surviving \(K_{S}\) fraction is \(\ee^{-293}\), a number smaller than \(10^{-127}\) and therefore not a background of any kind; the surviving \(K_{L}\) fraction is \(\ee^{-0.513}=0.599\), so more than half the long-lived beam is still available where the apparatus looks at it. The flight path is close to the optimum: shortening it would not increase the \(K_{L}\) flux much and would begin to admit \(K_{S}\), lengthening it would cost \(K_{L}\) flux exponentially for no gain.

The decay volume itself was a bag of helium at atmospheric pressure. The reason is Proposition 111.3: coherent regeneration scales as the number density of scatterers times the difference of their forward amplitudes, and both factors are minimized by a light, dilute gas. A true vacuum over a volume of several cubic metres was impractical in 1964; helium was the practical approximation to one, and the following estimate is what justifies the choice.

Proposition 111.5 (Helium suppresses the regeneration background).

Over a decay path of \(2\,\mathrm{m}\) of helium at atmospheric pressure, the coherently regenerated \(K_{S}\) amplitude is of order \(3\times10^{-5}\), that is about one per cent of the signal amplitude \(\abs{\eta_{+-}}=2.2\times10^{-3}\); in air it would have been several times larger. Rests on Equation (111.4).

Derivation. Derives Proposition 111.5. Insert into Equation (111.4). The reduced wavelength of the beam is

\[ \frac{\hbar}{p} =\frac{1.054571817\times 10^{-34}\,\mathrm{J}\,\mathrm{s}} {5.879\times 10^{-19}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}} =1.794\times 10^{-16}\,\mathrm{m}\ec \]

where \(p=1.1\,\mathrm{GeV}/c\) has been converted using \(1\,\mathrm{eV}=1.602176634\times 10^{-19}\,\mathrm{J}\) [Mohr:2025]. Helium at \(101.3\,\mathrm{kPa}\) and room temperature, taken as \(293\,\mathrm{K}\), has \(N=P/(k_{\mathrm{B}}T) =2.50\times 10^{25}\,/\mathrm{m}^{3}\) atoms; the familiar \(2.69\times 10^{25}\,/\mathrm{m}^{3}\) is the same pressure at \(273.15\,\mathrm{K}\). The amplitude difference \(\abs{f-\bar{f}}\) is not calculable from first principles here, but it cannot exceed the geometric scale of the nucleus, so taking \(\abs{f-\bar{f}}\sim10^{-15}\,\mathrm{m}\) is a generous bound. With \(\ell=2\,\mathrm{m}\),

\[ \abs{\rho}\sim\pi\times2.50\times 10^{25}\,/\mathrm{m}^{3} \times10^{-15}\,\mathrm{m}\times2\,\mathrm{m} \times1.794\times 10^{-16}\,\mathrm{m}=3\times10^{-5}\ep \]

Air is denser in nucleons by a factor \(1.20\,\mathrm{kg}/\mathrm{m}^{3}/0.166\,\mathrm{kg}/\mathrm{m}^{3}=7.2\) — both densities at \(293\,\mathrm{K}\) and \(101.3\,\mathrm{kPa}\), since setting a room-temperature air density against the tabulated \(273.15\,\mathrm{K}\) helium value of \(0.179\,\mathrm{kg}/\mathrm{m}^{3}\) would understate the ratio as \(6.7\) — and its nuclei are heavier, so \(\abs{f-\bar{f}}\) is larger by roughly \((14/4)^{2/3}\approx2.3\) per atom while the atom density rises by about a factor of two; the product is an increase of four or five. The regenerated amplitude adds coherently to \(\eta_{+-}\), so a regeneration amplitude that is one per cent of \(\eta_{+-}\) perturbs the measured rate by about two per cent, while one that is five per cent of it perturbs the rate by about ten per cent — the difference between a negligible correction and a systematic error comparable with the statistical one. This is an order-of-magnitude estimate and was treated as such: the experiment also measured the regeneration directly by inserting known material, which is the calibration described in Section 111.3.3.

The two-arm spectrometer

The detector was two symmetric arms viewing the helium bag, each consisting of a bending magnet with a spark chamber before it and a spark chamber after it, followed by a water Cherenkov counter and a scintillation counter [Christenson:1964]. Each arm measured one charged particle: its direction from the two spark chambers, its momentum from the deflection between them, and its identity, loosely, from the Cherenkov signal. The two arms together therefore delivered the two three-vectors \(\vect{p}_{1}\) and \(\vect{p}_{2}\) from which every quantity in the analysis is built.

The arms sat at about \(22\,^\circ\) to the beam on either side, and that angle is not arbitrary.

Proposition 111.6 (Why the arms sit near 22 degrees).

A \(K^{0}\) of momentum \(1.1\,\mathrm{GeV}/c\) decaying to \(\pi^{+}\pi^{-}\) sends both pions to the same laboratory angle \(20.5^\circ\), on opposite sides of the beam, when they are emitted transversely to the beam in the kaon rest frame. That symmetric configuration is the one a pair of symmetric arms can accept with both arms at the same angle, and it is what fixes that angle. A single pion can reach as far as \(41.9^\circ\) from the beam, attained at \(\cos\theta^{*}=-0.909\), but only when its partner has collapsed to \(4.9^\circ\). Rests on Proposition 111.4.

Derivation. Derives Proposition 111.6. In the kaon rest frame the two pions are back to back with momentum

\[ p^{*}c=\sqrt{\frac{m_{K}^{2}c^{4}}{4}-m_{\pi}^{2}c^{4}} =\sqrt{\left(248.806\,\mathrm{MeV}\right)^{2} -\left(139.570\,\mathrm{MeV}\right)^{2}} =205.97\,\mathrm{MeV}\ec \]

that is \(p^{*}=1.101\times 10^{-19}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}\), and each carries energy \(E^{*}=m_{K}c^{2}/2=248.806\,\mathrm{MeV}\). Boosting to the laboratory with \(\gamma=2.426\), \(\beta=0.9111\) from Proposition 111.4, a pion emitted at \(90^\circ\) in the rest frame keeps its transverse momentum \(p_{\perp}=p^{*}\) and acquires

\[ p_{\parallel}=\gamma\beta\frac{E^{*}}{c} =2.211\times248.806\,\mathrm{MeV}/c=550.1\,\mathrm{MeV}/c\ec \]

so its laboratory angle is \(\arctan\left(205.97/550.1\right)=20.5^\circ\). Its partner, emitted transversely in the opposite sense, arrives at the same \(20.5^\circ\) on the other side of the beam: this is the only configuration in which the two pions make equal angles with the beam, and it is therefore the configuration a symmetric two-arm spectrometer is built around.

The maximum laboratory angle of a single pion is a different number and is worth separating out. With \(u=\cos\theta^{*}\) the laboratory angle satisfies

\[ \tan\theta=\frac{p^{*}\sqrt{1-u^{2}}} {\gamma\left(p^{*}u+\beta E^{*}/c\right)} =\frac{1}{\gamma}\,\frac{\sqrt{1-u^{2}}}{u+\beta/\beta^{*}}\ec \qquad \beta^{*}=\frac{p^{*}c}{E^{*}}=0.8278\ep \]

Differentiating the last fraction gives \(-u(u+\beta/\beta^{*})-(1-u^{2})=0\), that is \(u=-\beta^{*}/\beta=-0.909\); a maximum exists precisely because \(\beta=0.9111\) exceeds \(\beta^{*}\), so that the origin lies outside the ellipse of laboratory momenta and the tangent from it touches at \(u=-0.909\) rather than at the transverse point. There \(\tan\theta=0.8968\) and \(\theta=41.9^\circ\), while the partner at \(u=+0.909\) is thrown forward to \(\tan\theta=0.0857\), that is \(4.9^\circ\). Two arms placed at \(22\,^\circ\) therefore sit on the symmetric configuration, accept a wide band of rest-frame angles about it at once, and keep clear of the beam; the wide-angle pairs are lost, but they are lost together with their forward partners and cost only acceptance, not bias.

The two remaining pieces of the apparatus are there to reject the common decays. The Cherenkov counters used water, of refractive index \(n=1.33\), whose threshold velocity is \(\beta_{\mathrm{th}}=1/n=0.752\); for a pion this is a momentum of

\begin{equation}\tag{111.5} p_{\mathrm{th}}=\frac{m_{\pi}c\,\beta_{\mathrm{th}}} {\sqrt{1-\beta_{\mathrm{th}}^{2}}} =159\,\mathrm{MeV}/c\ec \end{equation}

so essentially every pion in the accepted momentum band radiates while the counters remain useful in the trigger. The scintillators supplied the fast timing that defined the coincidence.

The resolutions matter more than any of this, because the entire argument of the experiment is a two-dimensional separation. The spectrometer reconstructed the invariant mass of the pair to about \(4\,\mathrm{MeV}/c^{2}\) and the direction of the reconstructed parent to a few milliradians; the signal bin used in the analysis, \(\cos\theta>0.9999\), corresponds to \(\theta<14.1\,\mathrm{mrad}\) (Table 111.1).

QuantityValueNote
Proton energy (AGS)\(30\,\mathrm{GeV}\)$4.8\times 10^{-9}\,\mathrm{J}$
Neutral-beam production angle\(30\,^\circ\)to circulating protons
Target-to-detector flight path\(17.4\,\mathrm{m}\)57 feet as published
Beam momentum at the peak$1.1\,\mathrm{GeV}/c$$5.879\times 10^{-19}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}$
Lorentz factor of the beam$\gamma=2.426$derived
$K_{S}$ decay length\(5.93\times 10^{-2}\,\mathrm{m}\)derived
$K_{L}$ decay length\(33.9\,\mathrm{m}\)derived
Flight path in $K_{S}$ decay lengths$293$derived
Decay volumehelium, \(101.3\,\mathrm{kPa}\)bag, not a vacuum vessel
Spectrometer arm angle\(22\,^\circ\)each side of the beam
Symmetric-pair laboratory angle\(20.5^\circ\)derived
Maximum single-pion angle\(41.9^\circ\)derived, at $\cos\theta^{*}=-0.909$
Rest-frame pion momentum$206\,\mathrm{MeV}/c$derived
Cherenkov radiatorwater, $n=1.33$pion threshold $159\,\mathrm{MeV}/c$
Invariant-mass resolution$\approx4\,\mathrm{MeV}/c^{2}$pion hypothesis
Signal angular bin$\cos\theta>0.9999$$\theta<14.1\,\mathrm{mrad}$
Parameters of the Christenson–Cronin–Fitch–Turlay apparatus and beam [Christenson:1964], with the derived kinematic quantities of Propositions 111.4 and 111.6. Imperial figures in the original report are converted; the \(17.4\,\mathrm{m}\) flight path is the 57 feet of the published layout. Quantities marked “derived” are computed in this chapter from the beam momentum and the particle masses of Table 111.5, not taken from the source.

Procedure

Trigger and event recording

An event was accepted when the scintillation and Cherenkov counters of both arms fired in coincidence, indicating one charged particle in each arm within the resolving time of the counters. The coincidence pulsed the spark chambers to a few tens of kilovolts, so that the ionization trails left by the two particles broke down into visible sparks, and triggered a camera which photographed the chambers [Christenson:1964]. The recorded object was therefore a photograph, and the analysis a measurement of positions on film, followed by a reconstruction.

This is worth stating plainly because it bounds the sample size. The experiment could not record a large fraction of the beam; it recorded about twenty-two thousand charged decays in total, and the signal it found is forty-five events. Everything about the argument — the use of a two-dimensional discriminant rather than a single cut, the reliance on sidebands for the background, the insistence on a control sample — follows from having a sample of that size and no possibility of enlarging it quickly.

Reconstruction in the mass-angle plane

From the two measured momenta the analysis formed two numbers. The first is the invariant mass of the pair under the hypothesis that both particles are pions,

\begin{equation}\tag{111.6} m_{\pi\pi}^{2}c^{4}=2m_{\pi}^{2}c^{4} +2\left(E_{1}E_{2}-\vect{p}_{1}\cdot\vect{p}_{2}c^{2}\right)\ec \qquad E_{i}=\sqrt{p_{i}^{2}c^{2}+m_{\pi}^{2}c^{4}}\ep \end{equation}

The second is the angle \(\theta\) between the reconstructed parent momentum \(\vect{p}_{1}+\vect{p}_{2}\) and the beam direction, reported as \(\cos\theta\).

The discriminating power lies in requiring both at once, and the reason is a kinematic identity rather than a statistical convention.

Proposition 111.7 (A two-body decay is a point in the mass-angle plane).

If a beam particle of mass \(m_{K}\) decays to exactly two pions, then \(m_{\pi\pi}=m_{K}\) and \(\cos\theta=1\) identically. If it decays to three or more bodies of which two are detected, then \(m_{\pi\pi}<m_{K}\), and \(\cos\theta<1\) unless the undetected momentum happens to vanish. Rests on Equation (111.6).

Derivation. Derives Proposition 111.7. For a two-body decay, four-momentum conservation gives \(p_{K}^{\mu}=p_{1}^{\mu}+p_{2}^{\mu}\), so the invariant mass of the pair is the parent mass and its three-momentum is the parent's, which points along the beam. Both statements are exact and both are independent of the decay angle, which is what makes the pair \((m_{\pi\pi},\cos\theta)\) a point rather than a distribution.

For a decay \(K\to\pi\pi X\) with \(X\) undetected of four-momentum \(q^{\mu}\) and mass \(m_{X}\ge0\),

\[ m_{\pi\pi}^{2}c^{4}=\left(p_{K}-q\right)^{2}c^{2} =m_{K}^{2}c^{4}+m_{X}^{2}c^{4}-2p_{K}\!\cdot\!q\,c^{2}\ep \]

Here a four-momentum carries the SI dimension of momentum, so that \(p^{2}=m^{2}c^{2}\) and every term above is an energy squared. In the kaon rest frame \(p_{K}\cdot q\,c^{2}=m_{K}c^{2}E_{X}^{*}\) with \(E_{X}^{*}\ge m_{X}c^{2}\), and \(E_{X}^{*}\) is bounded below by the kinematic minimum for the given final state, so \(m_{\pi\pi}<m_{K}\) strictly. The missing three-momentum likewise tilts \(\vect{p}_{1}+\vect{p}_{2}\) away from the beam, by an angle of order \(q_{\perp}/(p_{1}+p_{2})\); for the dominant \(K_{L}\) modes — \(\pi^{\pm}e^{\mp}\nu\), \(\pi^{\pm}\mu^{\mp}\nu\) and \(\pi^{+}\pi^{-}\pi^{0}\) — the missing transverse momentum is of order \(100\,\mathrm{MeV}/c\) against a total of order \(1\,\mathrm{GeV}/c\), giving angles of order \(100\,\mathrm{mrad}\), seven times the signal bin. The two variables therefore fail together for background and hold together for signal, and it is their joint behaviour, not the size of the excess, that carries the argument.

For the semileptonic modes there is a second, independent handicap. Because the electron or muon is assigned the pion mass in Equation (111.6), its energy is overestimated, and the reconstructed mass is pushed further from \(m_{K}\) rather than towards it — so the misassignment cannot manufacture a signal, only remove background from the peak.

Backgrounds and the control samples

Three backgrounds had to be dealt with, and each was measured rather than modelled.

Three-body decays under the peak. The overwhelming majority of the recorded events are \(\pi e\nu\), \(\pi\mu\nu\) and \(3\pi\), which populate the whole \((m_{\pi\pi},\cos\theta)\) plane. Their density in the signal region was fixed empirically from mass sidebands: the \(\cos\theta\) distribution was formed in the mass intervals immediately below and above the kaon mass, where no two-body signal can exist, and interpolated across. The interpolation is safe because the background varies smoothly in \(m_{\pi\pi}\) on the scale of the mass resolution while the signal does not.

Incoherent regeneration. Diffractive and inelastic regeneration in the residual material produce a \(K_{S}\) whose momentum is not collinear with the beam, and are therefore removed by the same angular requirement that defines the signal, as Proposition 111.7 makes explicit.

Coherent regeneration. This one is collinear and cannot be cut away, so it was measured. Slabs of material were deliberately placed in the beam, producing a large and unmistakable coherent \(K_{S}\to\pi^{+}\pi^{-}\) signal at \(m_{\pi\pi}=m_{K}\) and \(\cos\theta=1\); that sample served two purposes at once. It calibrated the mass and angular resolution of the spectrometer on real two-body decays of the right mass and the right momentum spectrum — the ideal control, since it is the signal in every respect except its origin — and it fixed the constant of proportionality in Equation (111.4) for the materials involved, from which the residual regeneration in helium follows by scaling. The observed forward peak did not scale with the amount of material in the beam, which is the direct experimental statement that it is not regeneration.

Finally, the normalization. The denominator of the published ratio is the number of all charged decays of the same beam recorded in the same apparatus. Taking a ratio to that sample cancels the beam flux, the live time, the film-scanning efficiency and a large part of the geometric acceptance, which is why the result could be quoted to twenty per cent from forty-five events.

Observations and data

The forward peak at the kaon mass

In the mass interval containing the kaon mass, and only there, the \(\cos\theta\) distribution shows a sharp excess in the single bin \(\cos\theta>0.9999\). In the two adjacent mass intervals, ten \(\mathrm{MeV}\) below and above, the same bin shows nothing above the smooth background. The excess amounts to \(45\pm9\) events, and the sample of all charged decays recorded in the same apparatus numbers \(22\,700\) [Christenson:1964].

Phenomenon 111.8 (The long-lived neutral kaon decays to two pions).

A neutral kaon beam allowed to propagate far enough that every short-lived component has decayed still produces two-charged-pion events whose reconstructed invariant mass is the kaon mass and whose reconstructed parent momentum is collinear with the beam: an excess of \(45\pm9\) events over an estimated background of about twelve in the forward bin, corresponding to

\begin{equation}\tag{111.7} \frac{\Gamma\left(K_{L}\to\pi^{+}\pi^{-}\right)} {\Gamma\left(K_{L}\to\text{all charged}\right)} =\frac{45\pm9}{22\,700}=(2.0\pm0.4)\times10^{-3} \end{equation}

[Christenson:1964]. The decay is forbidden if \(CP\) is a symmetry. Rests on Equations (105.46) and (105.47).

Derivation. Derives Phenomenon 111.8. The two pions from a spinless parent are in a state of zero relative orbital angular momentum, \(L=0\). Space inversion multiplies such a state by the product of the two intrinsic parities and by \((-1)^{L}\); the pion is pseudoscalar, so \(P\ket{\pi^{+}\pi^{-}}=(-1)^{2}(-1)^{L}\ket{\pi^{+}\pi^{-}} =+\ket{\pi^{+}\pi^{-}}\). Charge conjugation exchanges the two pions, which for relative angular momentum \(L\) contributes the same factor \((-1)^{L}=+1\). Hence

\[ CP\ket{\pi^{+}\pi^{-}}=+\ket{\pi^{+}\pi^{-}}\ec \]

the final state is \(CP\)-even. Before 1964 the long-lived kaon was identified with the \(CP\)-odd combination \(\ket{K_{2}}=\left(\ket{K^{0}}-\ket{\bar{K}^{0}}\right)/\sqrt{2}\), for which \(CP\ket{K_{2}}=-\ket{K_{2}}\), precisely because that assignment forbids the two-pion channel and so explains the observed factor of some six hundred between the two lifetimes. An interaction commuting with \(CP\) cannot connect states of opposite \(CP\) eigenvalue, so Equation (111.7) establishes that the weak interaction does not commute with \(CP\). What the observation does not by itself decide is whether the violation lies in the mixing that defines the long-lived state or in the decay amplitude; that is the question of Section 111.5.1.

The 1964 numbers and their modern successors are collected in Table 111.2. One comparison in that table deserves comment, because it is the only place where the original measurement can be checked against a modern one and the two are not quoted in the same variable. The published ratio has “all charged decays” in its denominator, whereas the modern branching fraction has all decays. The charged modes of the \(K_{L}\) — \(\pi^{\pm}e^{\mp}\nu\), \(\pi^{\pm}\mu^{\mp}\nu\), \(\pi^{+}\pi^{-}\pi^{0}\) and \(\pi^{+}\pi^{-}\) itself — account for \(0.8033\pm0.0014\) of the total width [Navas:2024], the quoted uncertainty being the quadrature combination of the uncertainties on those four branching fractions. It is negligible beside the twenty per cent statistical error on the 1964 ratio, so

\begin{equation}\tag{111.8} \mathcal{B}\left(K_{L}\to\pi^{+}\pi^{-}\right)_{1964} =(2.0\pm0.4)\times10^{-3}\times0.8033 =(1.6\pm0.3)\times10^{-3}\ec \end{equation}

against the present \((1.967\pm0.010)\times10^{-3}\) [Navas:2024]: agreement at \(1.2\) standard deviations. The original measurement was correct, and it was correct to twenty per cent on forty-five events.

QuantityValueYear
Charged decays recorded$22\,700$1964
Excess in the forward mass peak$45\pm9$ events1964
Estimated background in that bin$\approx12$ events1964
$\Gamma(K_{L}\to\pi^{+}\pi^{-})/ \Gamma(K_{L}\to\text{charged})$$(2.0\pm0.4)\times10^{-3}$1964
Converted branching fraction$(1.6\pm0.3)\times10^{-3}$1964
$\mathcal{B}(K_{L}\to\pi^{+}\pi^{-})$$(1.967\pm0.010)\times10^{-3}$2024
$\mathcal{B}(K_{L}\to\pi^{0}\pi^{0})$$(8.64\pm0.06)\times10^{-4}$2024
Charged fraction of the $K_{L}$ width$0.8033\pm0.0014$2024
The 1964 measurement and its modern counterparts. The first block is what Christenson, Cronin, Fitch and Turlay reported [Christenson:1964]; the second is the current world average [Navas:2024]. The conversion between the two normalizations is Equation (111.8). Uncertainties are as quoted by the sources; the 1964 error is dominated by the counting statistics of the forty-five-event excess.

The branching ratio and the parameter epsilon

The branching fraction is not the natural variable. What the theory predicts is an amplitude ratio, because the two-pion final state is reached from the \(CP\)-even admixture in the long-lived state, and the natural comparison is with the same final state reached from the short-lived one. The definitions are Equation (105.65):

\begin{equation}\tag{111.9} \abs{\eta_{+-}}^{2} =\frac{\Gamma\left(K_{L}\to\pi^{+}\pi^{-}\right)} {\Gamma\left(K_{S}\to\pi^{+}\pi^{-}\right)} =\frac{\mathcal{B}\left(K_{L}\to\pi^{+}\pi^{-}\right)\tau_{S}} {\mathcal{B}\left(K_{S}\to\pi^{+}\pi^{-}\right)\tau_{L}}\ec \end{equation}

the second form following because a branching fraction is a partial width divided by a total width and the total widths are \(\hbar/\tau\). Inserting \(\mathcal{B}(K_{L}\to\pi^{+}\pi^{-})=1.967\times10^{-3}\), \(\mathcal{B}(K_{S}\to\pi^{+}\pi^{-})=0.6920\) and the lifetimes of Equation (111.2) gives

\[ \abs{\eta_{+-}}^{2} =\frac{1.967\times10^{-3}\times8.954\times 10^{-11}\,\mathrm{s}} {0.6920\times5.116\times 10^{-8}\,\mathrm{s}}=4.98\times10^{-6}\ec \]

that is \(\abs{\eta_{+-}}=2.23\times10^{-3}\), which is the measured value in Table 111.3. The arithmetic is worth seeing once: the tiny branching fraction and the enormous lifetime ratio very nearly cancel, and what survives is a number of order \(10^{-3}\) that has nothing to do with either.

The phase of \(\eta_{+-}\) is measured independently, by observing the interference between the \(K_{S}\) and \(K_{L}\) amplitudes in a beam containing both, and it comes out at \(43.51\pm0.05^\circ\). That value is predicted, not fitted: Equation (105.59) gives \(\arg\epsilon=\arctan\left[2\Delta m_{K}c^{2}\tau_{S}/\hbar\right]\), and with Equation (111.3) and Equation (111.2),

\begin{equation}\tag{111.10} \arctan\frac{2\times3.484\times 10^{-6}\,\mathrm{eV}} {7.3510\times 10^{-6}\,\mathrm{eV}}=43.47^\circ\ec \end{equation}

where \(\hbar/\tau_{S}=7.3510\times 10^{-6}\,\mathrm{eV}\) is the \(K_{S}\) full width. Equivalently, dividing numerator and denominator by two, \(\arctan\left[\Delta m_{K}c^{2}/(\hbar/2\tau_{S})\right]\) with the half width \(3.6755\times 10^{-6}\,\mathrm{eV}\): the factor of two in Equation (105.59) is the one that turns the half width into the full width, and must not be applied twice.

The agreement of a measured phase with a number built out of a mass difference and a lifetime is one of the sharper consistency checks in the subject, and it is what licenses the identification of \(\eta_{+-}\) with the mixing parameter \(\epsilon\) of Equation (105.54) up to the small direct term of Section 111.6.

QuantityValueSI equivalentYear
$m_{K^{0}}c^{2}$\(497.611(13)\,\mathrm{MeV}\)$7.9726\times 10^{-11}\,\mathrm{J}$2024
$m_{K^{0}}$$8.8707\times 10^{-28}\,\mathrm{kg}$2024
$\tau_{S}$\(8.954(4)\times 10^{-11}\,\mathrm{s}\)2024
$\tau_{L}$\(5.116(21)\times 10^{-8}\,\mathrm{s}\)2024
$\Delta m_{K}c^{2}$\(3.484(6)\times 10^{-6}\,\mathrm{eV}\)$5.582\times 10^{-25}\,\mathrm{J}$2024
$\Delta m_{K}c^{2}/\hbar$$5.293\times 10^{9}\,/\mathrm{s}$2024
$\abs{\eta_{+-}}$$(2.232\pm0.011)\times10^{-3}$dimensionless2024
$\arg\eta_{+-}$\(43.51\pm0.05^\circ\)\(0.7594\,\mathrm{rad}\)2024
$\abs{\epsilon}$$(2.228\pm0.011)\times10^{-3}$dimensionless2024
$\delta_{L}$$(3.32\pm0.06)\times10^{-3}$dimensionless2024
$\mathrm{Re}(\epsilon'/\epsilon)$$(1.66\pm0.23)\times10^{-3}$dimensionless2024
Measured parameters of the neutral kaon system [Navas:2024], with SI equivalents where the customary unit is not SI. The mass values are those of the machine-readable Particle Data Group release; the mass difference, the lifetimes and the $CP$ parameters are from the same review. The last column gives the year of the compilation, not of the first measurement; the years of first measurement are in Table 111.4.

Charge asymmetry in semileptonic decay

The two-pion decay establishes that \(CP\) is violated; it does not by itself give a sign. The observable that does is the charge asymmetry of the semileptonic decays of the same long-lived state, and it deserves attention out of all proportion to its size.

Phenomenon 111.9 (Charge asymmetry in semileptonic kaon decay).

The long-lived kaon decays slightly more often to \(\pi^{-}e^{+}\nu_{e}\) than to \(\pi^{+}e^{-}\bar{\nu}_{e}\):

\begin{equation}\tag{111.11} \delta_{L}=\frac{\Gamma(\pi^{-}e^{+}\nu_{e}) -\Gamma(\pi^{+}e^{-}\bar{\nu}_{e})} {\Gamma(\pi^{-}e^{+}\nu_{e})+\Gamma(\pi^{+}e^{-}\bar{\nu}_{e})} =(3.32\pm0.06)\times10^{-3} \end{equation}

[Navas:2024]. Small as it is, this is an absolute, convention-free distinction between matter and antimatter: which sign of lepton charge is the commoner one can be communicated to a correspondent with whom we share no sample and no convention. Rests on Equation (105.56) and Theorem 105.34.

Derivation. Derives Phenomenon 111.9. Let the long-lived state be \(\ket{K_{L}}\propto(1+\epsilon)\ket{K^{0}} -(1-\epsilon)\ket{\bar{K}^{0}}\) with \(\abs{\epsilon}\ll1\), which is Equation (105.56). The selection rule \(\Delta S=\Delta Q\) allows \(K^{0}\to\pi^{-}e^{+}\nu_{e}\) and \(\bar{K}^{0}\to\pi^{+}e^{-}\bar{\nu}_{e}\) but not the charge-exchanged pair, and \(CPT\) makes the two allowed amplitudes equal in magnitude. The two rates are therefore proportional to \(\abs{1+\epsilon}^{2}\) and \(\abs{1-\epsilon}^{2}\), so

\[ \delta_{L}=\frac{\abs{1+\epsilon}^{2}-\abs{1-\epsilon}^{2}} {\abs{1+\epsilon}^{2}+\abs{1-\epsilon}^{2}} =\frac{2\,\mathrm{Re}\,\epsilon}{1+\abs{\epsilon}^{2}} \approx2\,\mathrm{Re}\,\epsilon\ep \]

With \(\abs{\epsilon}=2.228\times10^{-3}\) and its phase \(43.5^\circ\) [Navas:2024], \(2\,\mathrm{Re}\,\epsilon=2\times2.228\times10^{-3} \times\cos43.5^\circ=3.23\times10^{-3}\), against the measured \((3.32\pm0.06)\times10^{-3}\) of Equation (111.11). The asymmetry is thus the same \(\epsilon\) that Phenomenon 111.8 exhibits, seen in an entirely different channel — which is what makes the mixing interpretation testable rather than merely available.

The convention-free character of the statement is worth spelling out, since it is the one experimental result in this book that changes what can be said about the words used to describe the world. Before 1964 the labels “matter” and “antimatter” were conventional: \(C\) violation and \(P\) violation alone leave open the possibility that what we call matter is simply the mirror image of what a distant correspondent calls matter, because a mirror reflection can be compensated by an exchange of particles for antiparticles. Equation (111.11) closes that possibility. The instruction “prepare a long-lived neutral kaon beam, count the charge of the emitted lepton, and call the commoner sign positive” can be transmitted in words, needs no shared sample, and gives the same answer everywhere.

The measured CP-violating observables

The remaining sections of the chapter analyse one measurement at a time. Table 111.4 collects them all first, with the year and the source of each, so that the numbers can be read against one another; Table 111.5 gives the masses used throughout, and Table 111.6 the three mixing frequencies in SI. Every value is quoted as its source quotes it, with the uncertainties separated into statistical and systematic where the source separates them.

SystemObservableValueYearSource
$K$$\Gamma_{\pi\pi}/\Gamma_{\mathrm{charged}}$$(2.0\pm0.4)\times10^{-3}$1964[Christenson:1964]
$K$$\abs{\epsilon}$$(2.228\pm0.011)\times10^{-3}$2024[Navas:2024]
$K$$\delta_{L}$$(3.32\pm0.06)\times10^{-3}$2024[Navas:2024]
$K$$\mathrm{Re}(\epsilon'/\epsilon)$, NA31$(3.3\pm1.1)\times10^{-3}$1988[Burkhardt:1988]
$K$$\mathrm{Re}(\epsilon'/\epsilon)$, NA48$(18.5\pm4.5\pm5.8)\times10^{-4}$1999[Fanti:1999]
$K$$\mathrm{Re}(\epsilon'/\epsilon)$, KTeV$(28.0\pm3.0\pm2.8)\times10^{-4}$1999[AlaviHarati:1999]
$K$$\mathrm{Re}(\epsilon'/\epsilon)$, average$(1.66\pm0.23)\times10^{-3}$2024[Navas:2024]
$K$$A_{T}$, CPLEAR$(6.6\pm1.3\pm1.0)\times10^{-3}$1998[Angelopoulos:1998]
$B$$\sin2\beta$, BaBar$0.59\pm0.14\pm0.05$2001[Aubert:2001]
$B$$\sin2\beta$, Belle$0.99\pm0.14\pm0.06$2001[Abe:2001]
$B$$\sin2\beta$, average$0.708\pm0.011$2024[Navas:2024]
$B$$\Delta S_{T}^{+}$, BaBar$-1.37\pm0.14\pm0.06$2012[Lees:2012]
$B$$A_{CP}(B^{0}\to K^{+}\pi^{-})$$-0.080\pm0.007\pm0.003$2013[Aaij:2013]
$B_{s}$$A_{CP}(B_{s}^{0}\to K^{-}\pi^{+})$$0.27\pm0.04\pm0.01$2013[Aaij:2013]
$D$$\Delta A_{CP}$$(-15.4\pm2.9)\times10^{-4}$2019[Aaij:2019]
$n$$\abs{d_{n}}$ (\(90\,\mathrm{\%}\) CL)$<2.9\times 10^{-47}\,\mathrm{C}\,\mathrm{m}$2020[Abel:2020]
CKM$J$$\left(3.08^{+0.15}_{-0.13}\right)\times10^{-5}$2024[Navas:2024] [Charles:2005]
Cosmos$\eta$$(6.12\pm0.04)\times10^{-10}$2020[Aghanim:2020]
Measured $CP$- and $T$-violating observables, in the order they are discussed. Where two uncertainties are given the first is statistical and the second systematic. All the quantities are dimensionless except the neutron electric dipole moment. The world-average entries are the compilations of the Particle Data Group [Navas:2024], which absorb later data than the discovery papers in the same block.
Particle$mc^{2}$$mc^{2}$ (J)$m$ (kg)
$K^{0}$\(497.611(13)\,\mathrm{MeV}\)$7.9726\times 10^{-11}$$8.8707\times 10^{-28}$
$K^{\pm}$\(493.677(15)\,\mathrm{MeV}\)$7.9096\times 10^{-11}$$8.8006\times 10^{-28}$
$\pi^{\pm}$\(139.57039(18)\,\mathrm{MeV}\)$2.2362\times 10^{-11}$$2.4881\times 10^{-28}$
$\pi^{0}$\(134.9768(5)\,\mathrm{MeV}\)$2.1626\times 10^{-11}$$2.4062\times 10^{-28}$
$J/\psi$\(3.096900(6)\,\mathrm{GeV}\)$4.9618\times 10^{-10}$$5.5207\times 10^{-27}$
Masses used in the kinematic derivations of this chapter, taken from the machine-readable Particle Data Group release [Navas:2024], with the SI conversion made using the exact elementary charge and speed of light of the present SI [Mohr:2025]. The energy equivalent is the value the source quotes; the mass in kilograms is that value divided by $c^{2}$.
System$\Delta m\,c^{2}$$\Delta m\,c^{2}$ (J)$\Delta m\,c^{2}/\hbar$$\Delta m/m$
$K^{0}$\(3.484(6)\times 10^{-6}\,\mathrm{eV}\)$5.582\times 10^{-25}$$5.293\times 10^{9}\,/\mathrm{s}$$7.0\times 10^{-15}$
$B^{0}$\(3.334(13)\times 10^{-4}\,\mathrm{eV}\)$5.341\times 10^{-23}$$5.065\times 10^{11}\,/\mathrm{s}$$6.3\times 10^{-14}$
$B_{s}^{0}$\(1.1693(4)\times 10^{-2}\,\mathrm{eV}\)$1.8734\times 10^{-21}$$1.7765\times 10^{13}\,/\mathrm{s}$$2.2\times 10^{-12}$
Mixing frequencies of the three neutral-meson systems in which $CP$ violation has been established, expressed as an energy, as a frequency and as a mass [Navas:2024]. The oscillation frequency is $\Delta m\,c^{2}/\hbar$ throughout, with $\hbar=1.054571817\times 10^{-34}\,\mathrm{J}\,\mathrm{s}$ [Mohr:2025]. The fractional column is $\Delta m/m$, and it is the reason the neutral meson is the most sensitive interferometer available.

Interpretation

Indirect versus direct violation

The 1964 result establishes that the weak interaction does not commute with \(CP\). It does not say where the violation sits, and there are two places available.

It may sit in the state. If the long-lived eigenvector of the effective Hamiltonian is not the \(CP\)-odd combination \(\ket{K_{2}}\) but carries a small admixture \(\epsilon\ket{K_{1}}\) of the \(CP\)-even one, then the two-pion decay proceeds entirely from that admixture, using a decay amplitude that itself respects \(CP\) perfectly. This is \(CP\) violation in the mixing, and it is a single complex number \(\epsilon\) for the whole system.

Or it may sit in the decay. If the amplitudes for \(K^{0}\to\pi\pi\) and \(\bar{K}^{0}\to\pi\pi\) differ in more than the strong rescattering phase they are obliged to share, then even an exact \(CP\) eigenstate would decay asymmetrically. This is \(CP\) violation in the decay, or direct \(CP\) violation, and by Equation (105.72) it is measured by a second complex number \(\epsilon'\).

The distinction is not academic, and for thirty-five years it was the open question of the subject, because Wolfenstein had immediately supplied a hypothesis in which the second number is exactly zero [Wolfenstein:1964]. Suppose there exists a new interaction, enormously weaker than the weak interaction, which changes strangeness by two units and violates \(CP\) maximally. Acting once, it contributes to the off-diagonal element \(M_{12}\) of the mass matrix (Equation (105.50)) and so to \(\epsilon\); it never contributes to a decay, because the decays are \(\Delta S=1\) transitions. Such a “superweak” interaction reproduces every kaon measurement then available with one free parameter and predicts

\begin{equation}\tag{111.12} \epsilon'=0\ec\qquad \eta_{+-}=\eta_{00}=\epsilon\ec\qquad \arg\epsilon=\arctan\frac{2\Delta m_{K}c^{2}\tau_{S}}{\hbar}\ec \end{equation}

and the last of these is the phase relation verified numerically in Equation (111.10). That agreement is worth pausing on: the superweak hypothesis was not a desperate one, it fitted the phase correctly, and the experiment that killed it had to measure a difference between \(\eta_{00}\) and \(\eta_{+-}\) of a few parts in a thousand of a quantity that is itself a few parts in a thousand. That is a measurement at the \(10^{-6}\) level in the ratio of rates, and it is why it took thirty-five years. It is Section 111.6.

The CKM matrix and the Jarlskog invariant

The alternative to a new interaction is that the violation is already present in the interaction we have, and this is what Kobayashi and Maskawa proposed [Kobayashi:1973]. Their argument is arithmetic and is proved as Theorem 105.25: a unitary \(N\times N\) quark mixing matrix, after all the rephasings of the quark fields have been used up, retains \(N(N-1)/2\) real angles and \((N-1)(N-2)/2\) phases. For \(N=2\) — the Cabibbo case [Cabibbo:1963], one angle and no phase — the matrix can be made real, and a real Lagrangian conserves \(CP\). For \(N=3\) there is exactly one phase, and it cannot be rotated away.

The status of this as a prediction should not be softened. In 1973 three quarks were known. The measurement being explained was a branching ratio of \(2\times10^{-3}\) in a kaon. The inference was that three more quarks exist. The charm quark appeared the following year (Section 103.2.3), the bottom quark in 1977 and the top quark in 1995, and the Nobel award of 2008 to Kobayashi and Maskawa recognized an argument, not a discovery of a particle.

Two further pieces of apparatus, in the theoretical sense, make the phase measurable. The first is a parametrization in which the observed hierarchy of the mixing angles is manifest: Wolfenstein's expansion in \(\lambda=\sin\theta_{C}\) [Wolfenstein:1983], whose four parameters, measured, are Equation (103.24). The second is a rephasing-invariant measure, since the phase \(\delta\) itself depends on the convention chosen for the quark fields and is therefore not an observable. Jarlskog's quantity [Jarlskog:1985], defined at Definition 105.26, is the imaginary part of a quartet of matrix elements and is invariant; its measured value is Equation (105.78),

\begin{equation}\tag{111.13} J=\left(3.08^{+0.15}_{-0.13}\right)\times10^{-5}\ep \end{equation}

Geometrically \(J/2\) is the area of every unitarity triangle (Proposition 103.30), so the single statement “\(CP\) is violated” is the single statement “the triangle is not degenerate”. Numerically, \(J\) is small not because the phase is small — it is of order unity, \(\delta=65.70^\circ\) by Equation (103.25) — but because the mixing angles are. That distinction returns with consequences in Section 111.8.2.

Direct CP violation in the kaon system

NA31 and the first evidence

The observable that separates the two hypotheses of Section 111.5.1 is the double ratio of the four two-pion widths, and the reason it and not something simpler is that every quantity common to the four channels cancels in it: the kaon flux, the solid angle, most of the reconstruction efficiency. What does not cancel is the difference between reconstructing two charged tracks and reconstructing four photons, and that difference is the systematic error of every measurement in this section.

The first experiment to see a non-zero value was NA31 at the CERN Super Proton Synchrotron [Burkhardt:1988]. Protons of \(450\,\mathrm{GeV}\) struck a beryllium target; the resulting neutral beam was allowed to run down a long evacuated decay region for the \(K_{L}\) measurement, and for the \(K_{S}\) measurement a second target mounted on a movable train was brought close to the detector, so that the same apparatus saw \(K_{S}\) decays distributed over the same fiducial region. Alternating the two configurations, rather than running them at once, is what makes the acceptance cancel.

NA31 carried no magnetic spectrometer. Its detector was a liquid-argon and lead sampling calorimeter for the photons of \(\pi^{0}\pi^{0}\), an iron–scintillator hadron calorimeter, and wire chambers for the charged mode; the charged pion momenta were never measured individually, and the charged decays were identified from the calorimetric energy and the transverse position. The result,

\begin{equation}\tag{111.14} \mathrm{Re}\frac{\epsilon'}{\epsilon} =(3.3\pm1.1)\times10^{-3}\ec \end{equation}

is three standard deviations from zero and was the first evidence for direct \(CP\) violation.

It was not accepted, and the reason is that the contemporaneous Fermilab experiment E731, using a different technique — two parallel \(K_{L}\) beams with a regenerator placed alternately in one of them, so that \(K_{S}\) and \(K_{L}\) decays were recorded simultaneously — obtained a value consistent with zero, of order \(7\times10^{-4}\) with an uncertainty of similar size. The two results were about two standard deviations apart, which is not a crisis in itself; what made the decade that followed uncomfortable is that they disagreed about a qualitative question, the existence or not of \(\epsilon'\), and therefore about whether the superweak hypothesis had been excluded. Neither collaboration could settle it with the data in hand, and both built successors designed around the other's technique. That is the correct response to an irreconcilable pair of measurements, and it is worth recording as such.

Remark 111.10 (A source this book cannot cite).

The E731 result quoted above is described but not cited, because the bibliography of this treatise contains no entry for it. The number is given to an order of magnitude and its role in the argument is historical rather than evidential — the argument turns on the two 1999 measurements, which are cited — but the reader should know that this one sentence rests on an uncited source and treat it accordingly.

NA48 and KTeV

Both successors reported in 1999, and both found a non-zero value.

NA48, at CERN, produced \(K_{L}\) and \(K_{S}\) simultaneously rather than alternately [Fanti:1999]. The \(K_{L}\) beam came from a target \(126\,\mathrm{m}\) upstream of the decay region. A small fraction of the protons that missed that target were deflected by a bent silicon crystal onto a second target \(6\,\mathrm{m}\) upstream, producing a \(K_{S}\) beam nearly collinear with the first. The two decay populations were then recorded in the same detector at the same time, in the same beam conditions, with the acceptance difference reduced to the small difference in the decay distributions. Distinguishing a \(K_{S}\) decay from a \(K_{L}\) decay was done by timing: an array of scintillator counters — the tagger — in the proton beam heading for the \(K_{S}\) target recorded the passage of each proton to about \(200\,\mathrm{ps}\), and a decay coincident with a tagged proton is a \(K_{S}\) decay.

The detector was a magnetic spectrometer of four drift chambers around a dipole magnet, for \(\pi^{+}\pi^{-}\), together with a liquid-krypton electromagnetic calorimeter of about \(10\,\mathrm{m}^{3}\) and some thirteen thousand cells, twenty-seven radiation lengths deep, for the four photons of \(\pi^{0}\pi^{0}\). The krypton calorimeter is the critical instrument: the neutral mode has no tracks, the kaon's decay point must be reconstructed from the photon energies and positions alone by imposing the \(\pi^{0}\) mass, and the energy resolution therefore propagates directly into the fiducial volume definition. The 1999 result was

\begin{equation}\tag{111.15} \mathrm{Re}\frac{\epsilon'}{\epsilon} =(18.5\pm4.5\pm5.8)\times10^{-4}\ec \end{equation}

the second uncertainty systematic and larger than the first.

KTeV, at Fermilab, was the successor of E731 and kept its technique [AlaviHarati:1999]. Protons of \(800\,\mathrm{GeV}\) on a beryllium-oxide target produced two parallel neutral beams; a regenerator — a block of scintillator, itself instrumented so that inelastic interactions in it could be vetoed — sat in one of the two and was moved from one beam to the other about once a minute, so that the regenerated \(K_{S}\) sample and the \(K_{L}\) sample were recorded simultaneously and each beam served in turn as each. The detector was a drift-chamber spectrometer around an analysing magnet and an electromagnetic calorimeter of \(3100\) pure caesium-iodide crystals, whose energy resolution below one per cent for the photons of interest is what made the neutral mode competitive. The result was

\begin{equation}\tag{111.16} \mathrm{Re}\frac{\epsilon'}{\epsilon} =(28.0\pm3.0\pm2.8)\times10^{-4}\ep \end{equation}

The two 1999 numbers differ by \(9.5\times10^{-4}\) against a combined uncertainty of \(8.4\times10^{-4}\), that is by \(1.1\) standard deviations, which is unremarkable. Their distances from zero are not alike: combining each experiment's two errors in quadrature, KTeV is \(6.8\) standard deviations from zero and NA48 is \(2.5\). It is KTeV that established the effect in 1999; NA48 confirmed it, at a significance that on its own would not have settled the question, and the two together left no room for zero. With the later data of both experiments the world average is Equation (105.74), \(\mathrm{Re}(\epsilon'/\epsilon)=(1.66\pm0.23)\times10^{-3}\).

Phenomenon 111.11 (Direct CP violation in kaon decay).

\(CP\) violation occurs not only in the mixing that defines the long-lived state but in the decay amplitudes themselves. The double ratio of the four two-pion widths,

\begin{equation}\tag{111.17} \frac{\Gamma(K_{L}\to\pi^{0}\pi^{0})\,\Gamma(K_{S}\to\pi^{+}\pi^{-})} {\Gamma(K_{L}\to\pi^{+}\pi^{-})\,\Gamma(K_{S}\to\pi^{0}\pi^{0})} =1-6\,\mathrm{Re}\left(\epsilon'/\epsilon\right)\ec \end{equation}

differs from unity. Two experiments with independent systematics found \(\mathrm{Re}(\epsilon'/\epsilon)\) of order \(10^{-3}\) and non-zero [Fanti:1999] [AlaviHarati:1999], and the world average is \((1.66\pm0.23)\times10^{-3}\) [Navas:2024]. A superweak \(\Delta S=2\) interaction, which reproduces \(\epsilon\) exactly, predicts \(\epsilon'=0\) and is thereby excluded. Rests on Equation (105.71), Equation (105.72) and Proposition 105.39.

Derivation. Derives Phenomenon 111.11. Two steps: the double ratio in terms of \(\epsilon'/\epsilon\), and the reason \(\epsilon'\) requires two amplitudes with different strong phases.

The double ratio. Two pions in a relative \(S\) wave from a spinless parent have a symmetric isospin wavefunction, so only \(I=0\) and \(I=2\) are available; write the two amplitudes as \(A_{I}\ee^{\ii\delta_{I}}\) with \(\delta_{I}\) the real strong rescattering phase, as in Equation (105.67). Projecting the physical final states onto the isospin basis with the Clebsch–Gordan coefficients Equation (105.70) gives Equation (105.71),

\[ \eta_{+-}=\epsilon+\epsilon'\ec\qquad \eta_{00}=\epsilon-2\epsilon'\ec \]

the factor \(-2\) being nothing but the ratio of the weight of \(A_{2}\) in \(\ket{\pi^{0}\pi^{0}}\) to its weight in \(\ket{\pi^{+}\pi^{-}}\). The four widths of Equation (111.17) are \(\abs{\eta_{00}}^{2}\) over \(\abs{\eta_{+-}}^{2}\), since each \(\Gamma(K_{L}\to f)\) divided by \(\Gamma(K_{S}\to f)\) is the squared modulus of the corresponding amplitude ratio and everything else cancels. Hence

\[ \abs{\frac{\eta_{00}}{\eta_{+-}}}^{2} =\abs{\frac{1-2\epsilon'/\epsilon}{1+\epsilon'/\epsilon}}^{2} =\abs{1-3\frac{\epsilon'}{\epsilon} +O\!\left(\left(\epsilon'/\epsilon\right)^{2}\right)}^{2} =1-6\,\mathrm{Re}\frac{\epsilon'}{\epsilon} +O\!\left(10^{-6}\right)\ec \]

which is Equation (111.17). With the measured average the double ratio is \(1-0.00996\), so the experiments are looking for a one per cent departure from unity in a ratio of four separately measured rates — and the factor of six in front, which triples the naive \(-2\epsilon'/\epsilon\), is the only thing that makes it feasible at all.

Why a non-zero value needs two phases. From Equation (105.72),

\[ \epsilon'=\frac{\ii}{\sqrt{2}}\, \ee^{\ii(\delta_{2}-\delta_{0})}\, \frac{\mathrm{Re}\,A_{2}}{\mathrm{Re}\,A_{0}} \left[\frac{\mathrm{Im}\,A_{2}}{\mathrm{Re}\,A_{2}} -\frac{\mathrm{Im}\,A_{0}}{\mathrm{Re}\,A_{0}}\right]\ec \]

and the bracket is a difference. A phase common to both isospin amplitudes is a rephasing of the kaon state and carries no physics; it cancels in the bracket, leaving \(\epsilon'=0\). So direct \(CP\) violation requires at least two amplitudes with different weak phases. It requires more than that: the prefactor carries \(\ee^{\ii(\delta_{2}-\delta_{0})}\), and if the two strong phases were equal that factor would be real, \(\epsilon'\) would be purely imaginary relative to \(\epsilon\), and \(\mathrm{Re}(\epsilon'/\epsilon)\) — the observable of Equation (111.17) — would vanish even with \(\epsilon'\ne0\). Both a weak-phase difference and a strong-phase difference are necessary. This is the general statement proved at Proposition 105.39, and it is used again without change for the charm asymmetry in Section 111.7.3.

Remark 111.12 (Why the measured value is so small).

That \(\epsilon'/\epsilon\) is of order \(10^{-3}\) rather than of order unity is not a statement about the \(CP\)-violating phase. The prefactor of \(\epsilon'\) contains \(\mathrm{Re}\,A_{2}/\mathrm{Re}\,A_{0}\), and the \(\Delta I=\tfrac12\) rule — the empirical observation that \(\abs{A_{0}/A_{2}}\approx22\) in kaon decay, itself an unexplained strong-interaction fact — suppresses that ratio by a factor of twenty-two while simultaneously enhancing \(\epsilon\). Roughly five hundred of the three orders of magnitude come from that one number.

The Standard Model prediction and what it is worth

The measurement excludes the superweak hypothesis by itself, and that conclusion needs no theory. Comparing the measured value with the Cabibbo–Kobayashi–Maskawa prediction is a different and much weaker statement, and the reason should be stated plainly rather than left to a citation.

The short-distance part of the calculation is under control: the \(\Delta S=1\) effective Hamiltonian is a sum of local four-quark operators with coefficients computed perturbatively, and the phase enters only through the CKM combination \(\mathrm{Im}\left(V_{td}V_{ts}^{*}\right)\). The long-distance part is not: what is needed is the matrix element of each operator between a kaon and a two-pion state at the physical kinematics, including the final-state rescattering that produces \(\delta_{0}\) and \(\delta_{2}\), and there is no expansion parameter in which that is small. Lattice calculations of these matrix elements now exist and give a value compatible with Equation (105.74) with an uncertainty several times larger, dominated by the \(I=0\) channel where the rescattering is strongest; those calculations are not cited here because this treatise's bibliography carries no entry for them, and the reader should not take the agreement on the authority of this sentence.

What is established without any of that machinery is the qualitative result: \(\epsilon'\ne0\), hence \(CP\) is violated in a \(\Delta S=1\) decay amplitude, hence no interaction acting only in the mass matrix can account for the phenomenon. That is the content of Phenomenon 111.11 and it is secure.

Derivation pending.

The short-distance box-diagram computation of the \({[}\Delta S=2{]}\) mixing amplitude, giving the mass difference and epsilon in terms of the CKM elements and the quark masses — the step that converts the measured epsilon into a constraint on the unitarity triangle. It is an appendix-scale calculation requiring the loop functions of the box with internal charm and top, and it belongs in Appendix A rather than in an experiment chapter.

CP violation beyond the kaons

The B factories

If one phase in one matrix accounts for all of it, then the kaon's part in a thousand and an effect of order unity somewhere else are the same physics. The Standard Model says where to look: in \(B^{0}\to J/\psi\,K^{0}_{S}\), where the interference is between decay after mixing and decay without mixing, and no small ratio of amplitudes suppresses it.

The machines. Two asymmetric-energy electron–positron colliders were built for this measurement, PEP-II at SLAC and KEKB at KEK. Both run at a centre-of-mass energy of \(10.58\,\mathrm{GeV}\), the mass of the \(\Upsilon(4S)\), which decays to a \(B^{0}\bar{B}^{0}\) pair with almost no energy to spare — with \(m_{B^{0}}c^{2}=5279.66\,\mathrm{MeV}\) each \(B\) has a momentum \(p=\sqrt{(m_{\Upsilon}c^{2}/2)^{2}-(m_{B}c^{2})^{2}}/c =331\,\mathrm{MeV}/c\) in the \(\Upsilon\) frame. Asymmetric energies are the whole point. PEP-II collided \(9.0\,\mathrm{GeV}\) electrons with \(3.1\,\mathrm{GeV}\) positrons, KEKB \(8.0\,\mathrm{GeV}\) with \(3.5\,\mathrm{GeV}\), so that the centre of mass moves in the laboratory with

\begin{equation}\tag{111.18} \beta\gamma=\frac{E_{-}-E_{+}}{\sqrt{4E_{-}E_{+}}} =0.56\ \text{(PEP-II)}\ec\qquad 0.43\ \text{(KEKB)}\ep \end{equation}

The one line behind this: for head-on beams whose energies far exceed the electron mass, the total four-momentum has energy \(E_{-}+E_{+}\) and momentum \((E_{-}-E_{+})/c\) along the electron direction, so \(\beta=(E_{-}-E_{+})/(E_{-}+E_{+})\), while \((E_{-}+E_{+})^{2}-(E_{-}-E_{+})^{2}=4E_{-}E_{+}\) is the square of the invariant centre-of-mass energy, so \(\gamma=(E_{-}+E_{+})/\sqrt{4E_{-}E_{+}}\). The product is Equation (111.18), and \(\sqrt{4E_{-}E_{+}}\) is \(10.56\,\mathrm{GeV}\) for PEP-II and \(10.58\,\mathrm{GeV}\) for KEKB, as it must be.

The detectors. BaBar and Belle are conventional collider detectors of their generation, built around the interaction point inside a \(1.5\,\mathrm{T}\) solenoid: a silicon strip vertex detector of five layers immediately around the beam pipe, a gas drift chamber for momentum measurement, a Cherenkov device for separating kaons from pions — internally reflecting quartz bars in BaBar, aerogel threshold counters and time-of-flight in Belle — a caesium-iodide crystal calorimeter, and an instrumented flux return for muons. What the measurement demands of them is unusual, however: not energy resolution but position resolution along the beam axis, to better than \(100\,\mu\mathrm{m}\) per decay vertex.

Phenomenon 111.13 (A large CP asymmetry in B meson decay).

In \(\Upsilon(4S)\) decays to \(B^{0}\bar{B}^{0}\) pairs, the rate for \(B^{0}\to J/\psi\,K^{0}_{S}\) differs from that for \(\bar{B}^{0}\to J/\psi\,K^{0}_{S}\) by an amount that oscillates with the proper-time difference between the two \(B\) decays, with an amplitude of order unity; two experiments reported it simultaneously [Aubert:2001] [Abe:2001]. The parameter extracted is \(\sin2\beta=0.708\pm0.011\) [Navas:2024]. \(CP\) violation is therefore not intrinsically a small effect: in the kaon system it is a part in a thousand, in this \(B\) channel it is maximal in order of magnitude. Rests on Equation (105.87), Proposition 105.28 and Equation (105.82).

Derivation. Derives Phenomenon 111.13. Three things have to be shown: that the observable is a difference of decay times and not a decay time, that its amplitude is \(\sin2\beta\), and that the boost of Equation (111.18) makes the difference measurable.

The observable. The \(\Upsilon(4S)\) has \(J^{PC}=1^{--}\) and decays to two pseudoscalars in a \(P\) wave, so the pair is in the antisymmetric state Equation (105.87). Antisymmetry is preserved until one meson decays; therefore at every instant before the first decay there is exactly one \(B^{0}\) and exactly one \(\bar{B}^{0}\), and no oscillation is observable in the pair as a whole. The clock starts when one of them decays to a flavour-specific final state — a semileptonic decay, whose lepton charge tags the flavour, or a \(D^{(*)}\pi\) state — because that decay projects the partner onto the opposite flavour at that instant. The partner then oscillates freely until it decays to \(J/\psi K^{0}_{S}\), and the phase it has accumulated is \(\Delta m_{d}c^{2}\Delta t/\hbar\) with \(\Delta t=t_{CP}-t_{\mathrm{tag}}\), which may be of either sign. The asymmetry is therefore Equation (105.81),

\[ \mathcal{A}(\Delta t) =S_{f}\sin\frac{\Delta m_{d}c^{2}\Delta t}{\hbar} -C_{f}\cos\frac{\Delta m_{d}c^{2}\Delta t}{\hbar}\ec \]

odd in \(\Delta t\) for the \(CP\)-violating term, which is what makes the sign of \(\Delta t\) — and hence the vertex ordering — essential rather than incidental.

The amplitude. With \(\lambda=(q/p)\left(\bar{\mathcal{A}}_{f}/\mathcal{A}_{f}\right)\), the coefficients are \(S_{f}=2\mathrm{Im}\,\lambda/(1+\abs{\lambda}^{2})\) and \(C_{f}=(1-\abs{\lambda}^{2})/(1+\abs{\lambda}^{2})\) (Proposition 105.28). For \(f=J/\psi K^{0}_{S}\) the decay \(b\to c\bar{c}s\) proceeds through one dominant amplitude, so \(\abs{\lambda}=1\) and \(C_{f}=0\), and \(\lambda\) is a pure phase built from three CKM ratios: the \(B^{0}\) mixing, dominated by the top quark in the box, contributes \(q/p=V_{tb}^{*}V_{td}/\left(V_{tb}V_{td}^{*}\right)\); the decay amplitude ratio contributes \(V_{cb}V_{cs}^{*}/\left(V_{cb}^{*}V_{cs}\right)\); and the neutral kaon in the final state must itself be projected onto \(K^{0}_{S}\), which contributes \(V_{cd}V_{cs}^{*}/\left(V_{cd}^{*}V_{cs}\right)\). The product, with the \(CP\) eigenvalue \(\eta_{f}=-1\) of \(J/\psi K^{0}_{S}\), is

\[ \lambda=-\frac{V_{tb}^{*}V_{td}}{V_{tb}V_{td}^{*}}\, \frac{V_{cb}V_{cd}^{*}}{V_{cb}^{*}V_{cd}} =-\ee^{-2\ii\beta}\ec\qquad \beta=\arg\!\left(-\frac{V_{cd}V_{cb}^{*}}{V_{td}V_{tb}^{*}}\right)\ec \]

the factors of \(V_{cs}\) cancelling between the decay and the kaon. Since \(\abs{\lambda}=1\), \(S_{f}=\mathrm{Im}\,\lambda=\mathrm{Im}\left(-\ee^{-2\ii\beta}\right) =\sin2\beta\), which is Equation (105.82). Note what has and has not been assumed: unitarity of the mixing matrix, top dominance of the box, and a single decay amplitude. Nothing about hadronic matrix elements enters, which is why this channel measures an angle of the unitarity triangle to a per-cent theoretical accuracy while the kaon's \(\epsilon'\) does not.

The measurement of \(\Delta t\). A \(B\) meson lives \(\tau_{B}=1.519\times 10^{-12}\,\mathrm{s}\) and moves with the centre of mass, so the two decay vertices are separated along the beam axis by

\[ \Delta z=\beta\gamma c\,\Delta t\ec\qquad \avg{\abs{\Delta z}}\sim\beta\gamma c\tau_{B} =0.56\times2.998\times 10^{8}\,\mathrm{m}/\mathrm{s}\times1.519\times 10^{-12}\,\mathrm{s} =2.5\times 10^{-4}\,\mathrm{m} \]

at PEP-II, and \(1.9\times 10^{-4}\,\mathrm{m}\) at KEKB. Without the boost the two mesons would be nearly at rest and \(\Delta z\) would be a few micrometres, far below any silicon detector's resolution; with it, \(\Delta z\) is a quarter of a millimetre and a vertex resolution of about \(180\,\mu\mathrm{m}\) gives \(\sigma(\Delta t)=\sigma(\Delta z)/(\beta\gamma c) \approx1.1\times 10^{-12}\,\mathrm{s}\), comparable to \(\tau_{B}\) itself and to the oscillation period \(2\pi\hbar/(\Delta m_{d}c^{2})=1.24\times 10^{-11}\,\mathrm{s}\). The measurement is therefore possible but resolution-limited, which is why the fitted asymmetry must be convolved with a resolution function determined from the data themselves.

Dilution. The flavour tag is not perfect: a fraction \(w\) of tags is wrong, and each wrong tag contributes with the opposite sign, so the observed amplitude is \((1-2w)S_{f}\) and the statistical power of a sample of \(N\) events with tagging efficiency \(\varepsilon\) scales as \(\varepsilon(1-2w)^{2}N\). Both experiments used several tagging categories with different purities and fitted them simultaneously, reaching an effective efficiency \(\varepsilon(1-2w)^{2}\) of about \(0.3\). This is the dominant reason the first measurements, on samples of about thirty million \(B\bar{B}\) pairs each, carried statistical uncertainties of \(0.14\).

The two 2001 results were \(\sin2\beta=0.59\pm0.14\pm0.05\) from BaBar, on a sample of about thirty-two million \(B\bar{B}\) pairs [Aubert:2001], and \(0.99\pm0.14\pm0.06\) from Belle, on about thirty-one million [Abe:2001]. They differ by roughly two standard deviations, which is a useful reminder that first measurements scatter and that the agreement of two simultaneous results is not to be expected as a matter of course. Both are far from zero, which is the discovery; the average settled at \(0.708\pm0.011\) as the samples grew by a factor of ten [Navas:2024].

Remark 111.14 (The prediction this confirms).

The value was predicted before it was measured, from quantities having nothing to do with \(CP\) asymmetries. Taking the apex of the unitarity triangle from \(\abs{V_{ub}/V_{cb}}\), from the \(B^{0}\) and \(B_{s}^{0}\) oscillation frequencies and from \(\epsilon\) itself — Equation (103.24) — gives \(\beta=22.7^\circ\) and hence \(\sin2\beta=0.713\), against the measured \(0.708\pm0.011\). Three loop processes and one tree process, none of them an asymmetry, anticipate an asymmetry of order unity to better than one per cent. That is the strongest single test the Cabibbo–Kobayashi–Maskawa picture has passed.

Neutral $B_{s}$ mesons at the LHC

The \(B_{s}^{0}\) system, in which the light quark is strange rather than down, oscillates thirty-five times faster than the \(B^{0}\) (Table 111.6) and was for that reason the last of the neutral-meson systems to be brought under control.

The apparatus. LHCb is not a collider detector in the usual sense: it is a single-arm forward spectrometer, covering pseudorapidities from about \(2\) to \(5\), built on the observation that at the Large Hadron Collider the \(b\) and \(\bar{b}\) of a produced pair both go forward at small angles. Its elements, in order from the interaction point, are a silicon strip vertex locator whose sensors approach to within about \(8\,\mathrm{mm}\) of the beams, a tracking system around a warm dipole magnet of about \(4\,\mathrm{T}\,\mathrm{m}\) bending power, two ring-imaging Cherenkov detectors with gas radiators, which give kaon–pion separation over the momentum range from a few \(\mathrm{GeV}/c\) to about \(100\,\mathrm{GeV}/c\), calorimeters, and muon chambers. The two numbers that matter for flavour physics are an impact-parameter resolution of about \(20\,\mu\mathrm{m}\) and a decay-time resolution of about \(45\,\mathrm{fs}\) — the latter being what makes a \(B_{s}^{0}\) oscillation of period \(2\pi\hbar/(\Delta m_{s}c^{2})=3.5\times 10^{-13}\,\mathrm{s}\) resolvable at all.

That a hadron collider can do this at all rests on two things. The production cross section for \(b\bar{b}\) pairs is large, so statistics are never the limitation; and the \(b\) hadron flies a measurable distance — of order a centimetre at LHC momenta — before decaying, so a vertex detached from the primary interaction is by itself a powerful trigger against the enormous light-quark background.

The result. The first observation of \(CP\) violation in the decays of \(B_{s}^{0}\) mesons was of the direct kind, in the decay \(B_{s}^{0}\to K^{-}\pi^{+}\) [Aaij:2013]. The measured asymmetry,

\begin{equation}\tag{111.19} A_{CP}\left(B_{s}^{0}\to K^{-}\pi^{+}\right) =0.27\pm0.04\pm0.01\ec \end{equation}

is non-zero at \(6.5\) standard deviations, and the same analysis measured \(A_{CP}(B^{0}\to K^{+}\pi^{-})=-0.080\pm0.007\pm0.003\). Two features are worth noting. The asymmetry is of order a quarter — \(CP\) violation in a decay amplitude can be large when the interfering amplitudes are comparable, as tree and penguin are here. And the two asymmetries have opposite signs, which is a prediction of the Cabibbo–Kobayashi–Maskawa structure and not an accident: the ratio of the two, weighted by the branching fractions and lifetimes, is predicted to be \(-1\) by \(U\)-spin symmetry, the interchange of the \(d\) and \(s\) quarks. The reason in outline is that \(U\)-spin maps \(B^{0}\to K^{+}\pi^{-}\) onto \(B_{s}^{0}\to K^{-}\pi^{+}\) while leaving the strong phases and the moduli of the interfering tree and penguin amplitudes invariant and flipping the sign of the weak phase difference, so the two \(CP\) asymmetries have equal magnitude and opposite sign once each is normalized to its own width.

Derivation pending.

The \(U\)-spin relation between the direct \(CP\) asymmetries of \({[}B^{0}\to K^{+}\pi^{-}{]}\) and \({[}B_{s}^{0}\to K^{-}\pi^{+}{]}\) is stated here in outline only. Owed: the explicit tree-plus-penguin decomposition, the action of the \(d \leftrightarrow s\) interchange on each amplitude, and the width weighting that turns the amplitude statement into the ratio \(-1\), together with the size of the \(U\)-spin-breaking corrections. This belongs in Appendix A. No key in the bibliography covers the \(U\)-spin literature, so the claim carries no citation.

The null test. The observable usually associated with the \(B_{s}^{0}\) is not Equation (111.19) but the mixing phase \(\phi_{s}\), the analogue of \(2\beta\) for the strange-beauty system, measured from the time-dependent asymmetry in \(B_{s}^{0}\to J/\psi\,\phi\) and related channels. Its Standard Model value is small and sharply predicted, because in the Wolfenstein expansion the relevant combination \(V_{ts}V_{tb}^{*}\) is almost real; the world average, \(\phi_{s}\approx-0.05\,\mathrm{rad}\) with an uncertainty of about \(0.02\,\mathrm{rad}\) [Navas:2024], agrees with it. A small, precisely predicted quantity is the ideal place to look for a new short-distance contribution, since anything new entering the box would shift it; nothing has been seen, and the measurement is therefore a constraint rather than a discovery.

Charm

The up-type quarks were the last to yield, twenty years after the \(B\) mesons and fifty-five after the kaon.

Phenomenon 111.15 (CP violation in charm decay).

The difference between the time-integrated \(CP\) asymmetries of \(D^{0}\to K^{-}K^{+}\) and \(D^{0}\to\pi^{-}\pi^{+}\) is not zero,

\begin{equation}\tag{111.20} \Delta A_{CP}=A_{CP}\left(D^{0}\to K^{-}K^{+}\right) -A_{CP}\left(D^{0}\to\pi^{-}\pi^{+}\right) =(-15.4\pm2.9)\times10^{-4}\ec \end{equation}

established at \(5.3\) standard deviations [Aaij:2019]. With this the effect is observed in the strange, bottom and charm sectors alike; it has never been observed in any process that does not involve the weak interaction. Rests on Proposition 105.39.

Derivation. Derives Phenomenon 111.15. Why a difference is measured and not an asymmetry. A raw counting asymmetry between \(D^{0}\) and \(\bar{D}^{0}\) decays contains three contributions: the physics asymmetry, a production asymmetry, because a proton–proton collider does not produce charm and anticharm symmetrically, and a detection asymmetry, because the tagging pion of the \(D^{*+}\to D^{0}\pi^{+}\) chain is reconstructed with a different efficiency according to its charge. The last two are common to the two final states \(K^{-}K^{+}\) and \(\pi^{-}\pi^{+}\), both of which are self-conjugate, so they cancel in the difference Equation (111.20) to first order. This is the same design principle as the double ratio of Equation (111.17): build the observable so that everything not being measured appears in both terms.

What generates the asymmetry. A direct asymmetry requires two interfering amplitudes with both a weak and a strong phase difference, by Proposition 105.39. Here they are the tree amplitude \(c\to s\bar{s}u\) (or \(d\bar{d}u\)) and the loop, or “penguin”, amplitude in which the \(c\) converts to \(u\) via a \(W\) and a down-type quark. Write the total as \(\mathcal{A}=A_{1}\ee^{\ii(\varphi_{1}+\delta_{1})} +A_{2}\ee^{\ii(\varphi_{2}+\delta_{2})}\), with \(\varphi_{i}\) the weak phases, which change sign under \(CP\), and \(\delta_{i}\) the strong phases, which do not; the conjugate amplitude \(\bar{\mathcal{A}}\) is the same with \(\varphi_{i}\to-\varphi_{i}\). Then, with \(r=A_{2}/A_{1}\), \(\Delta\varphi=\varphi_{1}-\varphi_{2}\) and \(\Delta\delta=\delta_{1}-\delta_{2}\),

\begin{align*} \abs{\bar{\mathcal{A}}}^{2}-\abs{\mathcal{A}}^{2} &=2A_{1}A_{2}\left[\cos\left(\Delta\delta-\Delta\varphi\right) -\cos\left(\Delta\delta+\Delta\varphi\right)\right] =4A_{1}A_{2}\sin\Delta\delta\,\sin\Delta\varphi\ec\\ A_{CP}&=\frac{\abs{\bar{\mathcal{A}}}^{2}-\abs{\mathcal{A}}^{2}} {\abs{\bar{\mathcal{A}}}^{2}+\abs{\mathcal{A}}^{2}} =\frac{2r\sin\Delta\delta\,\sin\Delta\varphi} {1+2r\cos\Delta\delta\,\cos\Delta\varphi+r^{2}}\ep \end{align*}

Both sines must be non-zero. If the two amplitudes shared a weak phase it could be rotated away; if they shared a strong phase the two conjugate amplitudes would have equal moduli. For \(r\ll1\) this reduces to \(A_{CP}\approx2r\sin\Delta\delta\sin\Delta\varphi\), and with \(r\sim10^{-3}\) from the CKM factor \(\abs{V_{cb}V_{ub}^{*}/(V_{cs}V_{us}^{*})}\approx6\times10^{-4}\) and strong phases of order unity, an asymmetry of order \(10^{-3}\) is what the Standard Model leads one to expect — which is what was found.

Why this does not test the theory sharply. The expression above contains \(r\) and \(\Delta\delta\), and neither is calculable. The charm quark mass sits in the worst possible place: too heavy for chiral perturbation theory, too light for the heavy-quark expansion that makes \(B\) decays tractable, and the final states \(K\bar{K}\) and \(\pi\pi\) lie in a region dense with hadronic resonances that supply large and unpredictable strong phases. Estimates of the Standard Model value of Equation (111.20) span an order of magnitude and do not exclude the measured number at either end. The measurement therefore establishes the phenomenon in the up-type sector, which is a genuine addition to knowledge, and constrains rather than tests the Cabibbo–Kobayashi–Maskawa picture.

The apparatus and the sample. The measurement is LHCb's, using proton–proton collisions at centre-of-mass energies of \(7\), \(8\) and \(13\,\mathrm{TeV}\) and an integrated luminosity of about \(9\times 10^{43}\,/\mathrm{m}^{2}\). Two independent tagging methods were combined: the charge of the soft pion in \(D^{*+}\to D^{0}\pi^{+}\), and the charge of the muon in a semileptonic \(b\)-hadron decay producing the \(D^{0}\). The samples run to tens of millions of decays in each final state, and at that size the measurement is systematics-limited in practice — which is why the whole design is arranged so that the systematics cancel in a difference.

Direct T violation

Given the \(CPT\) theorem, every result above is also a statement about time reversal. But that inference assumes the theorem, whose hypotheses — locality, Lorentz invariance, the spin-statistics connection — are the deepest in the subject; a test of \(T\) that does not assume it is independent information. Two experiments have done it, and the difficulty in both is the same: a genuine \(T\) test compares a transition with its reverse, exchanging initial and final states, and conjugates nothing.

CPLEAR. The neutral kaon offers a transition and its exact reverse without any need to reverse momenta or spins, because both states are the same particle at rest: \(K^{0}\to\bar{K}^{0}\) and \(\bar{K}^{0}\to K^{0}\) are related by exchanging initial and final state, and nothing else [Kabir:1970]. The Kabir asymmetry Equation (105.84) must vanish identically if \(T\) is a symmetry.

The apparatus was built around the tagging problem, which is the whole difficulty: one must know the strangeness at production and at decay, in the same event. CPLEAR used the Low Energy Antiproton Ring at CERN, stopping antiprotons of about \(200\,\mathrm{MeV}/c\) in a gaseous hydrogen target at high pressure at a rate of order \(10^{6}\) per second. About two annihilations in a thousand give \(\bar{p}p\to K^{\mp}\pi^{\pm}K^{0}/\bar{K}^{0}\), and strangeness conservation in the strong production ties the sign of the accompanying charged kaon to the strangeness of the neutral one: the initial state is tagged by a charge measured in the same event. The final strangeness is tagged by the charge of the lepton in the semileptonic decay, using \(\Delta S=\Delta Q\) exactly as in Equation (111.11). The detector was a cylindrical assembly inside a solenoid of about \(0.44\,\mathrm{T}\) — tracking chambers, a threshold Cherenkov and scintillator sandwich for identifying the charged kaon, and a lead and gas-sampling electromagnetic calorimeter. The measured asymmetry Equation (105.86),

\begin{equation}\tag{111.21} A_{T}=(6.6\pm1.3\pm1.0)\times10^{-3}\ec \end{equation}

is the first direct observation of a time-reversal asymmetry [Angelopoulos:1998], and it agrees with the prediction \(A_{T}=4\,\mathrm{Re}\,\epsilon=6.46\times10^{-3}\) of Proposition 105.29 evaluated with \(\abs{\epsilon}=2.228\times10^{-3}\) and \(\arg\epsilon=43.5^\circ\).

BaBar. The cleaner measurement uses entanglement as a state-preparation device [Lees:2012]. From the antisymmetric \(B^{0}\bar{B}^{0}\) state Equation (105.87), observing one meson decay to a flavour-specific state projects the other onto a definite flavour at that instant, and observing a decay to \(J/\psi K^{0}_{S}\) projects it onto a definite \(CP\) label. The comparison is then available in both time orders: events in which the flavour tag comes first measure \(B^{0}\to B_{-}\), events in the opposite order measure \(B_{-}\to B^{0}\), and these are a process and its reverse with nothing conjugated. This is the entanglement machinery of Entanglement and Bell Tests used not to test locality but as an instrument. On a sample of \(468\) million \(B\bar{B}\) pairs BaBar found the \(T\)-odd parameters non-zero at fourteen standard deviations,

\begin{equation}\tag{111.22} \Delta S_{T}^{+}=-1.37\pm0.14\pm0.06\ec\qquad \Delta S_{T}^{-}=1.17\pm0.18\pm0.11\ec \end{equation}

while the \(CPT\)-odd parameters extracted from the same fit were consistent with zero [Lees:2012].

The \(CPT\) tests obtained as a by-product are worth recording, because they are among the sharpest tests of anything. The same formalism that yields Equation (111.21) bounds the \(CPT\)-violating parameter of the kaon mass matrix, and the resulting limit on the fractional mass difference between \(K^{0}\) and \(\bar{K}^{0}\) is of order \(10^{-18}\) — the most precise equality of a particle and antiparticle property in physics. The neutral-meson interferometer is what makes it possible, for the reason set out at Equation (105.63): the observable is an accumulated phase, not an energy. The details belong to Section 105.5.3.

Remark 111.16 (What these measurements do and do not assume).

Neither experiment is entirely free of assumptions, and both have been criticized on that ground. CPLEAR's final-state tag uses the semileptonic amplitudes, so a \(CPT\) violation in those amplitudes, or a failure of \(\Delta S=\Delta Q\), could in principle produce the observed asymmetry with no \(T\) violation in the mixing; the same data bound those parameters well below the measured \(A_{T}\), but the argument is internal. BaBar's mesons have a finite width, so the reversed transition occupies a different interval of proper time and is not an exact motion reverse; and the projection interpretation of the tag requires the tagging decay to be strictly flavour-specific. What is established, and what no earlier measurement established, is that the \(T\)-odd and \(CP\)-odd parameters are separately determined from one data set and are separately non-zero.

The unitarity triangle

The individual measurements above are of very different quantities in very different apparatus. What ties them together is one closure condition, the orthogonality of the first and third columns of the mixing matrix,

\begin{equation}\tag{111.23} V_{ud}V_{ub}^{*}+V_{cd}V_{cb}^{*}+V_{td}V_{tb}^{*}=0\ec \end{equation}

which is Equation (105.80). Three complex numbers summing to zero close a triangle in the complex plane; dividing by the middle term puts its base on the real axis from \(0\) to \(1\) and its apex at \((\bar{\rho},\bar{\eta})\). Every measurement in this chapter is a constraint on that apex, and they come in two kinds:

The measured angles are Equation (103.32), \(\alpha=85\pm5^\circ\), \(\beta=22.2\pm0.7^\circ\) and \(\gamma=66\pm3^\circ\); they sum to \(173^\circ\pm6^\circ\) against the identity \(\alpha+\beta+\gamma=\pi\), which is a consistency test the three measurements had no way of enforcing on one another. The global fits [Charles:2005] combine all of it, and the constraints cross in one region of the plane a few per cent across. The full account, with the fitted apex and every input, is Section 103.4.3 and Phenomenon 103.31; what belongs here is the experimental summary and the honest qualifications.

Those qualifications are three, and they are reported as they stand rather than smoothed over. The angle \(\alpha\) carries a \(5^\circ\) uncertainty dominated not by statistics but by the hadronic penguin amplitudes that must be subtracted from \(B\to\pi\pi\) — that is, by theory. The side \(\abs{V_{ub}/V_{cb}}\) inherits an unresolved tension: the inclusive determinations, which sum over all final states and use an operator product expansion, sit two to three standard deviations above the exclusive ones, which measure a single channel and need a lattice form factor, for both \(\abs{V_{ub}}\) and \(\abs{V_{cb}}\); choosing one or the other moves the apex by about one standard deviation. And the most precisely measured unitarity statement of all is not a triangle at all but the first-row relation \(\abs{V_{ud}}^{2}+\abs{V_{us}}^{2}+\abs{V_{ub}}^{2}=1\), which at present does not come out right at the \(10^{-3}\) level (Remark 103.32). The correct summary is that one phase accounts for all observed \(CP\) violation in the quark sector at the few-per-cent level, and that two specific tensions remain.

Cosmological relevance and open questions

The Sakharov conditions

The measurements of this chapter acquire a significance beyond particle physics because of an observation about the universe and an argument by Sakharov [Sakharov:1967].

Phenomenon 111.17 (The universe is made of matter).

No astrophysical concentration of antimatter has been found: neither the annihilation radiation that a matter–antimatter boundary would produce, nor antinuclei in the cosmic radiation beyond the rate expected from secondary production. The baryon-to-photon ratio inferred from the acoustic peaks of the cosmic microwave background is

\begin{equation}\tag{111.24} \eta=\frac{n_{B}-n_{\bar{B}}}{n_{\gamma}} =(6.12\pm0.04)\times10^{-10}\ec \end{equation}

in agreement with the value independently required by primordial nucleosynthesis [Aghanim:2020]. A universe that began symmetric under \(C\) and \(CP\), in thermal equilibrium and with conserved baryon number, could not have arrived at Equation (111.24). Rests on Theorem 105.41, Theorem 105.36 and Equation (105.105).

Derivation. Derives Phenomenon 111.17. The three requirements are proved as Theorem 105.41; what is carried out here is the arithmetic that turns them into a number, and the identification of which laboratory measurement supplies each.

Baryon number violation. Without it the initial value of \(n_{B}-n_{\bar{B}}\) is the final one, and a symmetric start stays symmetric. The Standard Model supplies it: the baryon current has an Adler–Bell–Jackiw anomaly, so its divergence is a total derivative of a topological density [tHooft:1976a], and configurations interpolating between vacua of different winding number change \(B\) by three units per generation. At zero temperature the process is suppressed by a barrier factor of order \(\ee^{-4\pi/\alpha_{W}}\) and is entirely unobservable; above the electroweak scale the barrier is crossed thermally by the sphaleron of Equation (105.104), and the rate is unsuppressed [Kuzmin:1985]. This condition is met.

\(C\) and \(CP\) violation. If \(C\) or \(CP\) were a symmetry, every baryon-number-violating process would be accompanied by its mirror at exactly the same rate and no net asymmetry could accumulate. This is the condition that the present chapter has measured, and it is met — that is the whole content of Phenomenon 111.8 and everything after it.

Departure from thermal equilibrium. In equilibrium the number density of a species depends only on its mass, chemical potential and the temperature, and \(CPT\) makes the masses of a particle and its antiparticle equal (Theorem 105.36); with \(B\) not conserved the corresponding chemical potential vanishes, so the two densities are equal whatever the microphysics. Any asymmetry generated is therefore erased as fast as it is made unless the expansion outruns the reactions.

The number. What must be explained is Equation (111.24), and that value is itself derived rather than observed directly: the acoustic peak structure of the microwave background measures \(\Omega_{b}h^{2}=0.02237\pm0.00015\) [Aghanim:2020], and dividing the resulting baryon number density by the blackbody photon density at \(T_{0}=2.7255\,\mathrm{K}\) gives \(\eta=2.74\times10^{-8}\,\Omega_{b}h^{2}\), which is Equation (105.105) and Equation (103.45). The independent determination from the light-element abundances (Evidence-Based Cosmology) agrees, which is what makes Equation (111.24) a measurement rather than a fit. The absence of antimatter concentrations is a separate observational statement: a matter–antimatter boundary anywhere within the observable universe would radiate annihilation gamma rays at a level long since excluded, and no search of the cosmic radiation has produced an antinucleus heavier than the antiproton, whose observed flux is consistent with secondary production on the interstellar medium.

Insufficiency of Standard Model CP violation

All three Sakharov conditions are formally present in the Standard Model. Two of them fail quantitatively, and the failures are independent of each other.

The phase is far too small, and it is the masses that make it so. \(CP\) violation in the quark sector is controlled not by \(J\) alone but by the rephasing-invariant combination Equation (105.79), which multiplies \(J\) by the six squared-mass differences and therefore carries the SI dimension \(\mathrm{J}^{12}\). An asymmetry generated at temperature \(T\) is a pure number, so it can depend on that invariant only through the dimensionless ratio Equation (105.106), obtained by dividing by \((k_{B}T)^{12}\). Evaluated with the running quark masses at \(k_{B}T\approx100\,\mathrm{GeV}\), where the sphaleron processes are active, this gives Equation (105.107),

\begin{equation}\tag{111.25} d_{CP}\approx2\times10^{-21}\ec \end{equation}

against the measured \(\eta=6.12\times10^{-10}\). The ratio is

\begin{equation}\tag{111.26} \frac{\eta}{d_{CP}} =\frac{6.12\times10^{-10}}{2\times10^{-21}} =3\times10^{11}\ec \end{equation}

a shortfall of eleven orders of magnitude, and the detailed calculations [Gavela:1994] do not improve it. The reason is structural rather than numerical. \(CP\) violation in this theory requires all six quark masses to be distinct (Proposition 105.27); at \(100\,\mathrm{GeV}\) five of the six are negligible against the temperature, so the effect is suppressed by twelve powers of small mass ratios. The phase itself is of order unity — \(\delta=65.70^\circ\) — and the masses kill it. This is the same arithmetic as Equation (103.46), and it is quoted here in the same form so that the two chapters cannot drift apart.

The third condition fails outright. A departure from equilibrium at the electroweak scale requires the transition from the symmetric to the broken phase to be first order, proceeding by the nucleation of bubbles whose walls sweep the plasma — and strongly first order at that, since the sphalerons must switch off fast enough after the transition not to erase what was just made. Lattice computations locate the endpoint of the first-order line at a scalar mass of about \(72\,\mathrm{GeV}/c^{2}\) [Kajantie:1996]; beyond it the change is a smooth crossover with no bubbles, no walls and no departure from equilibrium. The measured Higgs mass, Equation (104.66), is \(125.20 \pm 0.11\,\mathrm{GeV}/c^{2}\) [Navas:2024], far above the endpoint (Remark 104.65). No amount of \(CP\) violation would help.

Remark 111.18 (Stated as a failure, not as a puzzle solved).

Mechanisms that would generate Equation (111.24) exist in outline — heavy Majorana neutrinos decaying out of equilibrium, additional scalars driving a first-order transition — and not one of them has any experimental support. The shortfall itself is recorded in Section 133.5, which states that no confirmed mechanism supplies it, and Remark 103.45 makes the same point from the flavour side. What is measured is Equation (111.24); what is calculated is Equation (111.25); the eleven orders of magnitude between them are one of the few places where the Standard Model is known, from laboratory numbers plus one cosmological observation, to be incomplete.

Electric dipole moments and the strong CP problem

Every measurement so far in this chapter is of a \(CP\) violation that exists. The last one is of a \(CP\) violation that does not, and the null result is the more constraining of the two kinds.

Phenomenon 111.19 (The neutron has no electric dipole moment).

Repeated searches over seven decades find no permanent electric dipole moment of the neutron. The current bound is

\begin{equation}\tag{111.27} \abs{d_{n}}<2.9\times 10^{-47}\,\mathrm{C}\,\mathrm{m} \end{equation}

at \(90\,\mathrm{\%}\) confidence [Abel:2020], improving on [Baker:2006]. The same limit is customarily written \(1.8\times10^{-26}\) in units of the elementary charge times the centimetre, the conversion being \(1\,e\,\mathrm{cm}=1.602176634\times 10^{-21}\,\mathrm{C}\,\mathrm{m}\). This is a null result about a symmetry, not about a nuclear property. Rests on Theorem 105.30 and Equation (105.91).

Derivation. Derives Phenomenon 111.19. A neutron in a state of definite spin has no vector at its disposal other than its spin: by the projection theorem carried by the Wigner–Eckart theorem, the expectation value of any vector operator in a nondegenerate state of given angular momentum is proportional to \(\avg{\vect{S}}\). Hence \(\avg{\vect{d}}=d_{n}\avg{\vect{S}}/\abs{\avg{\vect{S}}}\) for some real number \(d_{n}\), and the interaction with a static electric field is

\[ H_{d}=-\avg{\vect{d}}\cdot\vect{E} \propto-d_{n}\,\avg{\vect{S}}\cdot\vect{E}\ep \]

Under space inversion \(\vect{E}\mapsto-\vect{E}\) while \(\vect{S}\mapsto\vect{S}\), so \(H_{d}\mapsto-H_{d}\); under time reversal \(\vect{E}\mapsto\vect{E}\) while \(\vect{S}\mapsto-\vect{S}\), so again \(H_{d}\mapsto-H_{d}\). A non-zero \(d_{n}\) therefore requires the violation of \(P\) and of \(T\) separately, and — given the \(CPT\) theorem of Discrete Symmetries and CPT — of \(CP\). The full statement, with the \(C\) properties made explicit, is Theorem 105.30.

The nondegeneracy clause is what makes the statement non-trivial: a state with an exact degeneracy of opposite parity, such as the hydrogen \(n=2\) level, shows a linear Stark effect with no symmetry violation at all, because there the dipole matrix element is off-diagonal and the argument above does not apply (Remark 105.31).

From the bound to a frequency. The experiment measures a frequency, not a dipole. Ultracold neutrons are stored in a cell in parallel magnetic and electric fields and their Larmor precession is measured by Ramsey's method of separated oscillatory fields [Ramsey:1950]; reversing \(\vect{E}\) and differencing isolates \(\Delta\nu=4d_{n}E/h\) (Equation (105.91)). At the Paul Scherrer Institute storage field \(E\approx1.1\times 10^{6}\,\mathrm{V}/\mathrm{m}\) a dipole equal to Equation (111.27) would shift the resonance by

\[ \Delta\nu=\frac{4\times2.9\times 10^{-47}\,\mathrm{C}\,\mathrm{m} \times1.1\times 10^{6}\,\mathrm{V}/\mathrm{m}} {6.62607015\times 10^{-34}\,\mathrm{J}\,\mathrm{s}}=1.9\times 10^{-7}\,\mathrm{Hz}\ec \]

against a Larmor frequency of about \(29\,\mathrm{Hz}\) at the working field of \(1\,\mu\mathrm{T}\): the bound is a frequency ratio measured to \(7\times10^{-9}\), twice, with the field reversed. It is not a theoretical estimate of any kind.

The number is worth converting once, because it is otherwise unimaginable. A dipole of \(2.9\times 10^{-47}\,\mathrm{C}\,\mathrm{m}\) corresponds to separating one elementary charge from its opposite by \(1.8\times 10^{-28}\,\mathrm{m}\), about \(10^{13}\) times smaller than the neutron itself, whose charge radius is of order \(0.8\,\mathrm{fm}\).

What the null result constrains. Two things, and they are of very different character. The first is new physics: the Cabibbo–Kobayashi–Maskawa phase generates a neutron dipole only at three loops with a severe chiral suppression, some six orders of magnitude below Equation (111.27).

Derivation pending.

The three-loop, chirally suppressed Cabibbo–Kobayashi–Maskawa contribution to the neutron electric dipole moment is quoted here as an order of magnitude and is not derived. Owed: the argument that the one- and two-loop quark diagrams cancel by the Glashow–Iliopoulos–Maiani mechanism because the rephasing invariant needs all six quark masses, and the resulting estimate of order \(3\times 10^{-53}\,\mathrm{C}\,\mathrm{m}\). This belongs in Appendix A. The bibliography carries no key for the electric-dipole calculation, so the estimate is uncited.

The entire experimental range is one in which the Standard Model predicts nothing observable. A non-zero result would be unambiguous, and the null result excludes new \(CP\)-violating sources over a wide range with no background to subtract.

The second is a parameter of the Standard Model itself. Quantum chromodynamics admits a term proportional to \(\theta\,\tr\left(G_{\mu\nu}\tilde{G}^{\mu\nu}\right)\) which is a total derivative, invisible in perturbation theory, and physical nonperturbatively because the gauge field has topologically distinct sectors. Only the combination \(\bar{\theta}=\theta+\arg\det M_{q}\) is invariant under the chiral rotations that move phase between the gluon and Yukawa sectors (Equation (105.111)), and \(\bar{\theta}\) violates \(P\) and \(T\), hence generates a neutron dipole. The chiral estimates [Baluni:1979] [Crewther:1979] give \(d_{n}\approx\bar{\theta} \times1.6\times 10^{-37}\text{–}6.4\times 10^{-37}\,\mathrm{C}\,\mathrm{m}\) (Equation (105.112)), so Equation (111.27) implies

\begin{equation}\tag{111.28} \abs{\bar{\theta}}\lesssim \frac{2.9\times 10^{-47}\,\mathrm{C}\,\mathrm{m}}{1.6\times 10^{-37}\,\mathrm{C}\,\mathrm{m}} \approx1.8\times10^{-10}\ec \end{equation}

that is \(\abs{\bar{\theta}}\lesssim10^{-10}\), with a factor-of-a-few uncertainty inherited from the chiral coefficient (Equation (105.113)).

This is the strong \(CP\) problem, and what makes it a problem is worth stating precisely, because it is unlike the other open questions in this chapter. \(\bar{\theta}\) is a free dimensionless parameter, periodic with period \(2\pi\), so there is no sense in which a small value sits near a boundary of its range; a value of order unity is the natural expectation for a number about which nothing else is known. It is measured to be smaller than \(10^{-10}\). No symmetry of the Standard Model requires it to vanish — the one candidate, a massless up quark, is excluded by the lattice determination \(m_{u}c^{2}=2.16 \pm 0.07\,\mathrm{MeV}\) [Navas:2024] — and no dynamical mechanism within the theory relaxes it. The problem is an observation, and it is recorded as such in Section 133.1.

The best-known response is to make \(\bar{\theta}\) dynamical: a global symmetry, spontaneously broken and itself anomalous, replaces it by a field whose instanton-generated potential is minimized at zero [Peccei:1977], at the cost of a light pseudo-Goldstone boson, the axion [Weinberg:1978] [Wilczek:1978]. No axion has been observed. The haloscope, helioscope and light-shining-through-walls searches report exclusions and nothing else, and the mechanism is a hypothesis with an attractive structure and no experimental support; it is listed as such in Section 133.7 and named here only to say so. What this chapter records as evidence is Equation (111.28): a free parameter that could have been of order unity is smaller than \(10^{-10}\), and nothing explains why.

Primary references

The literature of \(CP\) violation divides into papers that predicted a system, papers that measured something in it, and papers that supplied the framework in which the measurement means something. It is worth separating them, because the reader who goes to the sources will find that several of the most-cited ones do not contain what they are cited for.

The system. Gell-Mann and Pais [GellMann:1955] is the paper that invented the neutral kaon as a two-state problem and predicted the long-lived partner; it is a short paper and it is worth reading for the argument alone, which uses nothing beyond the superposition principle and a conservation law. Lande and collaborators [Lande:1956] reported the long-lived particle at the Brookhaven Cosmotron the following year. Pais and Piccioni [Pais:1955] predicted regeneration; the paper is the origin of the technique that produces \(K_{S}\) beams and of the background that Proposition 111.5 bounds.

The discovery. Christenson, Cronin, Fitch and Turlay [Christenson:1964] is two and a half pages of Physical Review Letters, and every quantitative statement about the 1964 apparatus in Sections 111.2 and 111.3 comes from it. The paper is unusual in modern eyes for how much of its argument is carried by one two-dimensional plot and how carefully it enumerates the alternatives it must exclude. Note that the published ratio has “all charged decay modes” in its denominator, not the total width; the conversion is Equation (111.8), and a reader comparing the 1964 number directly with a modern branching fraction will find a spurious twenty per cent discrepancy.

The frameworks. Landau [Landau:1957] proposed \(CP\) as the surviving symmetry after the parity results of Wu [Wu:1957] and Garwin [Garwin:1957]; Lüders [Luders:1957] proved the \(CPT\) theorem which is what makes \(CP\) violation and \(T\) violation the same statement. Wolfenstein [Wolfenstein:1964] supplied the superweak alternative in the same year as the discovery — a model that survived thirty-five years of measurement and was killed by \(\epsilon'\). Kobayashi and Maskawa [Kobayashi:1973] extended Cabibbo's angle [Cabibbo:1963] to three generations and showed a phase then becomes irremovable; the parametrization in which the hierarchy is manifest is Wolfenstein's second contribution [Wolfenstein:1983], and the convention-independent measure of the effect is Jarlskog's [Jarlskog:1985]. Kabir [Kabir:1970] is the source of the observable that CPLEAR measured.

The epsilon-prime programme. NA31 [Burkhardt:1988] reported the first evidence; the contemporaneous Fermilab result is described in Section 111.6.1 but not cited, since this treatise's bibliography carries no entry for it, and Remark 111.10 says so. The two 1999 papers that settled the question are Fanti and collaborators for NA48 [Fanti:1999] and Alavi-Harati and collaborators for KTeV [AlaviHarati:1999]. The world average is the Particle Data Group's [Navas:2024]; the lattice calculations to which it should be compared are not cited, for the same reason, and Section 111.6.3 says so.

Beyond the kaon. BaBar [Aubert:2001] and Belle [Abe:2001] reported the \(B^{0}\to J/\psi K^{0}_{S}\) asymmetry simultaneously in 2001; the two papers are worth reading side by side, since they disagree by two standard deviations while both establishing the effect. LHCb's first observation of \(CP\) violation in \(B_{s}^{0}\) decays [Aaij:2013] and in charm [Aaij:2019] complete the set of quark flavours. CPLEAR [Angelopoulos:1998] and BaBar [Lees:2012] are the two direct \(T\) tests. The global fit that combines all of it is CKMfitter [Charles:2005].

The cosmological side. Sakharov [Sakharov:1967] is a two-page note and states the three conditions in a form that has not needed amendment. 't Hooft [tHooft:1976a] showed that baryon number is anomalous in the electroweak theory; Kuzmin, Rubakov and Shaposhnikov [Kuzmin:1985] showed that the resulting rate is unsuppressed at high temperature. Gavela and collaborators [Gavela:1994] carried out the calculation that quantifies the shortfall, and Kajantie and collaborators [Kajantie:1996] showed by lattice computation that the electroweak transition is a crossover for any Higgs mass above about \(72\,\mathrm{GeV}/c^{2}\), which removes the third Sakharov condition independently. The measured baryon-to-photon ratio is Planck's [Aghanim:2020].

The null results. The current neutron electric dipole moment bound is from ultracold neutrons at the Paul Scherrer Institute [Abel:2020], improving the Institut Laue-Langevin result [Baker:2006]; the technique is Ramsey's [Ramsey:1950] and the argument for measuring the quantity at all was made by Purcell and Ramsey [Purcell:1950] at a time when parity was universally assumed. The chiral estimates relating the bound to \(\bar{\theta}\) are Baluni [Baluni:1979] and Crewther, Di Vecchia, Veneziano and Witten [Crewther:1979]. Peccei and Quinn [Peccei:1977] proposed the dynamical relaxation of \(\bar{\theta}\) and Weinberg [Weinberg:1978] and Wilczek [Wilczek:1978] pointed out the axion that follows; none of it has been observed.

The data. Every world-average number in this chapter is from the Particle Data Group's 2024 Review [Navas:2024], and the masses in Table 111.5 are from its machine-readable release, stored with this treatise as evidence. The SI conversions — the elementary charge, the speed of light, the Planck and reduced Planck constants — are CODATA 2022 [Mohr:2025], likewise stored; in the present SI the first three of those are exact, so the conversions in this chapter carry no uncertainty of their own.