§8.8CP Violation in the Decay

Part III Bettini pp. 338–342 · ~20 min read

  • ε′
  • η₊₋ and η₀₀
  • double ratio
  • weak phase vs strong phase

Two conditions have to hold at once before a decay asymmetry can exist at all, and arranging for both is why this measurement took thirty years.

🎯 Why this matters

A non-zero ε′ moved CP violation out of the states and into the interaction itself. Mixing alone could have been an accident of how two particles happen to overlap; an asymmetry between decay amplitudes cannot be.

§8.5 established CP violation in the mixing: the propagating states are not CP eigenstates, and the impurity is ε\varepsilon. This section asks the independent question — do the decay amplitudes themselves differ between a process and its CP conjugate?

The answer is yes, by a further factor of a thousand, and extracting it took from 1970 to 2002.

Two conditions, and both are necessary

Start with the general statement, because it governs §8.9 and §8.10 as well. Suppose a decay proceeds through two amplitudes. Each carries a weak phase, which flips sign under CP, and a strong phase from final-state interactions, which does not:

Af=A1ei(δS1+δW1)+A2ei(δS2+δW2),Aˉfˉ=A1ei(δS1δW1)+A2ei(δS2δW2)A_f = |A_1|e^{i(\htmlClass{t-s}{\delta_{S1}}+\htmlClass{t-w}{\delta_{W1}})} + |A_2|e^{i(\htmlClass{t-s}{\delta_{S2}}+\htmlClass{t-w}{\delta_{W2}})}, \qquad \bar A_{\bar f} = |A_1|e^{i(\htmlClass{t-s}{\delta_{S1}}-\htmlClass{t-w}{\delta_{W1}})} + |A_2|e^{i(\htmlClass{t-s}{\delta_{S2}}-\htmlClass{t-w}{\delta_{W2}})}
(8.80)

Bettini p. 342. The only difference between the two lines is the sign of the weak phases. Everything follows from that.

Every symbol, one at a time

Hover or tap a symbol above — it lights up in the equation and its meaning, units and type appear here.

Subtract the two rates and almost everything cancels:

Af2Aˉfˉ2=4A1A2sin(δS1δS2)sin(δW1δW2)|A_f|^2 - |\bar A_{\bar f}|^2 = \htmlClass{t-pre}{-4\,|A_1||A_2|}\, \htmlClass{t-ss}{\sin(\delta_{S1}-\delta_{S2})}\,\htmlClass{t-sw}{\sin(\delta_{W1}-\delta_{W2})}
(8.63)

The master formula for a direct CP asymmetry. Two sines multiplied together — and the whole of experimental CP violation is the problem of making neither of them zero.

Every symbol, one at a time

Hover or tap a symbol above — it lights up in the equation and its meaning, units and type appear here.

A product of two sines. Set either to zero and there is no asymmetry, however large the other is. That is the section’s central claim, and it is worth playing with rather than reading:

A rate asymmetry needs two phase differences, not one

ReImA₁AĀ
|A|²
0.8911
|Ā|²
2.5076
asymmetry (|A|² − |Ā|²) / (|A|² + |Ā|²)
-0.4756

Both present. Flipping the weak phase reflects A₂ about the direction of A₁, and because the two are not collinear the reflected chain closes somewhere else. The rate difference is −4|A₁||A₂| sin(δ_S1−δ_S2) sin(δ_W1−δ_W2) = -1.6164, matching |A|² − |Ā|² = -1.6164.

Two amplitudes reaching one final state. Under CP the weak phases flip, which reflects the second phasor about the direction of the first — so if they are collinear (no strong phase difference) nothing changes, and if the weak phases are equal there is nothing to reflect. Drag either slider to zero.

⚙️ Engineer’s bridge — a phase is only observable against a reference, and CP has to move one and not the other

The reason both conditions are needed is not an accident of the algebra. It is the same reason a single phasor has no measurable phase at all.

An overall phase on an amplitude is unobservable — eiθA2=A2|e^{i\theta}A|^2 = |A|^2 for any θ\theta. To see a phase you need two amplitudes and their interference, which measures only the difference. That is the first condition: one amplitude gives you nothing.

Now ask what CP does. It flips the weak phases and leaves the strong ones. So the CP conjugate rate differs only if flipping the weak phases changes the relative angle between the two phasors — and if the two phasors are collinear to begin with, reflecting one about the other’s direction puts it back where it was. Hence the second condition.

An engineer will recognise the structure as a lock-in detection with a reference channel that must be in quadrature. A lock-in multiplies your signal by a reference and low-passes; if the signal and reference are exactly in phase the output is insensitive to small phase shifts, and if they are 90° apart it is maximally sensitive. The strong phase is that quadrature: it is nature’s reference channel, and sin(δS1δS2)\sin(\delta_{S1}-\delta_{S2}) is exactly the quadrature factor.

Where it breaks: the lock-in analogy assumes you supply the reference. Here nature supplies it, and you cannot calculate it. A strong phase you cannot compute sits directly in front of the quantity you want. ε\varepsilon' is proportional to sin(δ2δ0)\sin(\delta_2-\delta_0), and that number comes from pion scattering data, not from the CKM matrix. Direct CP violation is therefore never a clean CKM measurement in the way sin2β\sin2\beta was (§8.6) — which is exactly why the chapter’s cleanest result came from interference rather than from decay.

The kaon case: two isospin amplitudes

For K2πK\to2\pi the two amplitudes are supplied by isospin. The two pions are in a spatially symmetric state, so their isospin wave function must be symmetric too: I=0I = 0 or I=2I = 2, never 1. Call the corresponding weak amplitudes A0A_0 and A2A_2. Then (Eq. 8.67)

A(K0π+π)=13(eiδ2A2+2eiδ0A0),A(K0π0π0)=13(2eiδ2A2eiδ0A0)A(K^0\to\pi^+\pi^-) = \tfrac{1}{\sqrt3}\left(\htmlClass{t-a2}{e^{i\delta_2}A_2} + \htmlClass{t-cg}{\sqrt2}\,\htmlClass{t-a0}{e^{i\delta_0}A_0}\right), \qquad A(K^0\to\pi^0\pi^0) = \tfrac{1}{\sqrt3}\left(\htmlClass{t-cg}{\sqrt2}\,\htmlClass{t-a2}{e^{i\delta_2}A_2} - \htmlClass{t-a0}{e^{i\delta_0}A_0}\right)
(8.67)

The two physical decays written in terms of two isospin amplitudes. This is the change of basis that supplies the 'two paths to one final state' the master formula above demands.

Every symbol, one at a time

Hover or tap a symbol above — it lights up in the equation and its meaning, units and type appear here.

and the conjugate amplitudes are the same with A0,2A0,2A_{0,2} \to A_{0,2}^* — which is CPT, not an assumption about the decay.

Experimentally A0A_0 dominates: A2/A01/22|A_2|/|A_0| \approx 1/22, the ΔI = 1/2 rule. That single number can be checked against the measured rate ratio, and doing so turns up something:

🔢 Worked example — the K_S rate ratio, from Eq. (8.67) and two numbers

Square the two amplitudes above, keeping the interference term:

A(π+π)2A(π0π0)2=r2+2+22rcos(δ2δ0)2r2+122rcos(δ2δ0),rA2/A0\frac{|A(\pi^+\pi^-)|^2}{|A(\pi^0\pi^0)|^2} = \frac{r^2 + 2 + 2\sqrt2\,r\cos(\delta_2-\delta_0)} {2r^2 + 1 - 2\sqrt2\,r\cos(\delta_2-\delta_0)}, \qquad r \equiv |A_2|/|A_0|

With r=1/22r = 1/22 and δ2δ0=48.3°\delta_2-\delta_0 = 48.3°, and multiplying by the phase-space ratio p(π+π)/p(π0π0)=0.9855p^*(\pi^+\pi^-)/p^*(\pi^0\pi^0) = 0.9855:

Γ(KSπ+π)Γ(KSπ0π0)=2.24\frac{\Gamma(K_S\to\pi^+\pi^-)}{\Gamma(K_S\to\pi^0\pi^0)} = \mathbf{2.24}

Now the measurement. The PDG branching ratios are 69.20%69.20\,\% and 30.69%30.69\,\%, so the ratio is

0.69200.3069=2.255±0.005\frac{0.6920}{0.3069} = \mathbf{2.255} \pm 0.005

Agreement to 0.7 % — which simultaneously confirms the ΔI = 1/2 rule and the strong phase difference, neither of which was used to obtain the branching ratios. Note that pure I=0I=0 alone would give only 1.97; the 14 % excess is the A2A_2 admixture.

Erratum — Eq. (8.68) prints 2.55 where 2.255 belongs

The book gives

Γ(KSπ+π)/Γ(KSπ0π0)=2.55±0.005\Gamma(K_S\to\pi^+\pi^-)/\Gamma(K_S\to\pi^0\pi^0) = 2.55 \pm 0.005

Four independent checks say this is a transposed digit for 2.255:

  1. The PDG branching ratios the book cites give 2.2548. 0.6920/0.30690.6920/0.3069.
  2. The quoted uncertainty fits 2.255 and not 2.55. Propagating ±0.0005\pm0.0005 on each branching ratio gives ±0.004\pm0.004 — the printed ±0.005\pm0.005. On a value of 2.55 that would be a 0.2 % measurement of nothing in particular.
  3. The book’s own Eq. (8.67) predicts 2.24, using the book’s own A2/A01/22|A_2/A_0| \approx 1/22 and δ2δ0=48.3°\delta_2-\delta_0 = 48.3°. That is 0.7 % from 2.255 and 12 % from 2.55.
  4. 2.55 would require A2/A0=1/10.8|A_2/A_0| = 1/10.8, contradicting the sentence immediately before the equation.

Nothing downstream uses the number, so no conclusion changes. But it is quoted precisely as evidence that the final state is nearly pure I=0I = 0, and at 2.55 that evidence points the wrong way.

ε′, and why it is so hard

With the Wu–Yang convention (A0A_0 real and positive), the two observable amplitude ratios η₊₋ and η₀₀ come out as, together with the direct-violation parameter ε′ ,

η+=ε+ε,η00=ε2ε,ε=i2ImA2A0ei(δ2δ0)\htmlClass{t-pm}{\eta_{+-} = \varepsilon + \varepsilon'}, \qquad \htmlClass{t-00}{\eta_{00} = \varepsilon - 2\varepsilon'}, \qquad \varepsilon' = \frac{i}{\sqrt2}\frac{\htmlClass{t-im}{\mathrm{Im}\,A_2}}{\htmlClass{t-a0}{A_0}}\htmlClass{t-str}{e^{i(\delta_2-\delta_0)}}
(8.71)

The two observable ratios, and the parameter that separates CP violation in the decay from CP violation in the mixing. All three conditions of the master formula are visible in the last expression.

Every symbol, one at a time

Hover or tap a symbol above — it lights up in the equation and its meaning, units and type appear here.

Read ε\varepsilon' carefully — it is zero unless A2A_2 is both non-zero and non-real. Non-zero is the ΔI = 1/2 rule being imperfect; non-real is a weak phase. And the ei(δ2δ0)e^{i(\delta_2-\delta_0)} is the strong phase that makes the real part observable. All three conditions of the bridge, in one expression.

The trouble is the size. If CP violation lived only in the mixing then η+=η00\eta_{+-} = \eta_{00} exactly, so the search is for a difference between two nearly equal numbers:

η00/η+=0.9950±0.0007,ϕ00ϕ+=(0.01±0.07)°|\eta_{00}/\eta_{+-}| = 0.9950 \pm 0.0007, \qquad \phi_{00} - \phi_{+-} = (0.01 \pm 0.07)°

💡 What this really says — why the observable is a double ratio and not anything simpler

The natural thing to measure would be Γ(K0π+π)\Gamma(K^0\to\pi^+\pi^-) against Γ(Kˉ0π+π)\Gamma(\bar K^0\to\pi^+\pi^-) — a decay and its conjugate, differing by 2Reε2\,\mathrm{Re}\,\varepsilon'. You cannot: by the time anything decays, the propagating states are KSK_S and KLK_L, not K0K^0 and Kˉ0\bar K^0.

So the observable has to be built from KSK_S and KLK_L rates. Take the ratio of π+π\pi^+\pi^- to π0π0\pi^0\pi^0 for each, then the ratio of those — the double ratio :

Γ(KLπ+π)Γ(KSπ+π)/Γ(KLπ0π0)Γ(KSπ0π0)=η+η0021+6Re(ε/ε)\htmlClass{t-r1}{\frac{\Gamma(K_L\to\pi^+\pi^-)}{\Gamma(K_S\to\pi^+\pi^-)}} \Big/ \htmlClass{t-r2}{\frac{\Gamma(K_L\to\pi^0\pi^0)}{\Gamma(K_S\to\pi^0\pi^0)}} = \left|\frac{\eta_{+-}}{\eta_{00}}\right|^2 \simeq \htmlClass{t-res}{1 + 6\,\mathrm{Re}(\varepsilon'/\varepsilon)}
(8.74)

The double ratio — the most systematics-driven observable in the book. Every element of its design exists so that something cancels.

Every symbol, one at a time

Hover or tap a symbol above — it lights up in the equation and its meaning, units and type appear here.

Every source of systematic error that is common to a pair cancels in that pair’s ratio, and everything common to the two pairs cancels again in the double ratio. Incident kaon flux, detector acceptance for a given topology, trigger efficiency — all gone.

That is why the experiments were built the way §8.8 describes: simultaneous detection of both final states so the flux cancels, and KLK_L and KSK_S beams present at the same time with matched energy spectra and matched decay distributions so the acceptance cancels. Every element of the design exists to make a cancellation exact.

The engineering instinct is the same one behind a Wheatstone bridge or a four-terminal resistance measurement: when the quantity you want is a small difference between two large ones, do not measure the two and subtract. Build a configuration in which everything you do not want cancels before the measurement, and read out only the residual.

Re ⁣(εε)=16[1η00η+2]=16[1Γ(KLπ0π0)/Γ(KSπ0π0)Γ(KLπ+π)/Γ(KSπ+π)]\mathrm{Re}\!\left(\frac{\varepsilon'}{\varepsilon}\right) = \frac{1}{6}\left[1 - \left|\frac{\eta_{00}}{\eta_{+-}}\right|^2\right] = \frac{1}{6}\left[1 - \htmlClass{t-d}{\frac{\Gamma(K_L\to\pi^0\pi^0)/\Gamma(K_S\to\pi^0\pi^0)} {\Gamma(K_L\to\pi^+\pi^-)/\Gamma(K_S\to\pi^+\pi^-)}}\right]
(8.77)

Bettini p. 341. Four decay rates, arranged so that everything except direct CP violation cancels.

Every symbol, one at a time

Hover or tap a symbol above — it lights up in the equation and its meaning, units and type appear here.

from the measured double ratio to ε′

eta_ratio, eps = 0.9950, 2.232e-3
d = 1 - eta_ratio**2
print("the double ratio and what it gives:")
print(f"  |eta_00 / eta_+-| = {eta_ratio} +- 0.0007")
print(f"  1 - |eta_00/eta_+-|^2 = {d:.3e}")
print(f"  Re(eps'/eps) = (1/6) x {d:.3e} = {d/6:.3e}")
print( "  book (8.78): (1.66 +- 0.23)e-03      -- exact agreement")

print("\nso how big is direct CP violation, absolutely?")
print(f"  |eps|          = {eps:.3e}   (violation in the MIXING)")
print(f"  Re(eps'/eps)   = {d/6:.3e}")
print(f"  |eps'|         ~ {eps*d/6:.2e}   (violation in the DECAY)")

print("\nthree levels, each about a thousand times smaller than the last:")
print( "  CP-conserving amplitude              1")
print(f"  mixing violation      |eps|       {eps:.1e}")
print(f"  decay violation       |eps'|      {eps*d/6:.1e}")
prints
the double ratio and what it gives:
|eta_00 / eta_+-| = 0.995 +- 0.0007
1 - |eta_00/eta_+-|^2 = 9.975e-03
Re(eps'/eps) = (1/6) x 9.975e-03 = 1.662e-03
book (8.78): (1.66 +- 0.23)e-03      -- exact agreement

so how big is direct CP violation, absolutely?
|eps|          = 2.232e-03   (violation in the MIXING)
Re(eps'/eps)   = 1.662e-03
|eps'|         ~ 3.71e-06   (violation in the DECAY)

three levels, each about a thousand times smaller than the last:
CP-conserving amplitude              1
mixing violation      |eps|       2.2e-03
decay violation       |eps'|      3.7e-06

Thirty years, and a standoff

The measurement is worth telling as history because it shows what a 10310^{-3} effect on top of a 10310^{-3} effect actually costs.

Bettini §8.8. The 1993 pair disagreed at 3.5σ and the question stayed open for nine more years.
experimentyearresult
NA31 (CERN)19932.30 ± 0.653.5σ from zero — a claim of discovery, on 428 000 KLπ0π0K_L\to\pi^0\pi^0 decays
E731 (FNAL)1993−0.74 ± 0.56similar statistics, consistent with zero. The two disagreed with each other at 3.5σ
KTeV (FNAL)20032.07 ± 0.28107\sim10^7 π0π0\pi^0\pi^0 decays and much better systematics
NA48 (CERN)20021.47 ± 0.22the same, independently — and the effect was settled

the 1993 standoff, quantified

import numpy as np
m = [('NA31', 2.30, 0.65), ('E731', -0.74, 0.56), ('KTeV', 2.07, 0.28), ('NA48', 1.47, 0.22)]

print("NA31 vs E731, 1993:")
d, sd = 2.30 - (-0.74), np.hypot(0.65, 0.56)
print(f"  2.30 +- 0.65  against  -0.74 +- 0.56")
print(f"  difference {d:.2f} +- {sd:.3f}  ->  {d/sd:.1f} sigma apart")
print( "  one claimed a 3.5 sigma discovery, the other was consistent with zero.")

print("\nKTeV vs NA48, the next generation:")
d2, sd2 = 2.07 - 1.47, np.hypot(0.28, 0.22)
print(f"  2.07 +- 0.28  against  1.47 +- 0.22")
print(f"  difference {d2:.2f} +- {sd2:.3f}  ->  {d2/sd2:.1f} sigma apart")
print( '  the book says "the two values agree"; 1.7 sigma is marginal but fair.')

wsum = sum(v/s**2 for _, v, s in m); w = sum(1/s**2 for _, v, s in m)
avg, err = wsum/w, 1/np.sqrt(w)
chi2 = sum(((v-avg)/s)**2 for _, v, s in m)
print(f"\nnaive weighted average of all four: {avg:.2f} +- {err:.2f}")
print(f"  chi^2 = {chi2:.1f} for 3 dof, so the four are NOT mutually consistent;")
print( "  the PDG inflates the error for exactly this reason, reaching")
print( "  book (8.78): 1.66 +- 0.23")
prints
NA31 vs E731, 1993:
2.30 +- 0.65  against  -0.74 +- 0.56
difference 3.04 +- 0.858  ->  3.5 sigma apart
one claimed a 3.5 sigma discovery, the other was consistent with zero.

KTeV vs NA48, the next generation:
2.07 +- 0.28  against  1.47 +- 0.22
difference 0.60 +- 0.356  ->  1.7 sigma apart
the book says "the two values agree"; 1.7 sigma is marginal but fair.

naive weighted average of all four: 1.54 +- 0.16
chi^2 = 21.6 for 3 dof, so the four are NOT mutually consistent;
the PDG inflates the error for exactly this reason, reaching
book (8.78): 1.66 +- 0.23

Three experimental obstacles made it that hard, and all three are named in the book:

  • KLπ0π0K_L\to\pi^0\pi^0 has BR =103= 10^{-3}, and sits under KLπ0π0π0K_L\to\pi^0\pi^0\pi^0, which is 200 times more frequent and gives the same photons plus two more. Losing two photons out of six turns the background into the signal.
  • The two kaons have wildly different decay lengths. At 110 GeV, λL=3.4\lambda_L = 3.4 km against λS=6\lambda_S = 6 m — yet the two decay distributions inside the detector must be made as similar as possible, or the acceptance does not cancel.
  • Everything must be simultaneous. Both final states at once so the flux cancels; both beams at once with matched spectra so the acceptance does.

Erratum — Eq. (8.81) is missing a factor of 2

The book prints

Af2Aˉfˉ2=2A1A2sin(δS1δS2)sin(δW1δW2)|A_f|^2 - |\bar A_{\bar f}|^2 = -2\,|A_1||A_2|\sin(\delta_{S1}-\delta_{S2})\sin(\delta_{W1}-\delta_{W2})

The coefficient is −4. It is pure algebra from Eq. (8.80): the cross term in z1+z22|z_1+z_2|^2 already carries a factor 2, and

cos(X+Y)cos(XY)=2sinXsinY\cos(X+Y) - \cos(X-Y) = -2\sin X\sin Y

supplies another. 2×2=42 \times 2 = 4.

Checked numerically over 20 000 random choices of the two magnitudes and four phases: with 4-4 the maximum residual is 2×10142\times10^{-14}; with 2-2 as printed it is 1717 — i.e. the printed relation simply fails.

Nothing in the argument depends on it. The point of Eq. (8.81) is the product of two sines, and that is right as printed; only the overall coefficient is wrong, and the equation is never used numerically.

Nine years of not knowing — Re(ε′/ε) × 10³, four experiments

zero — no direct CP violation−2−10+1+2+3Re(ε′/ε) × 10³NA31, 1993E731, 1993NA48, 2002KTeV, 20033.5σ apart

Supplied — the table above has the four numbers and the disagreement only becomes obvious when they share an axis. Two experiments of comparable statistics, published the same year, disagreeing at 3.5σ: one claiming a discovery, the other consistent with no direct CP violation at all. Neither was wrong about its own data; the gap was systematics, which is what a 10⁻³ effect measured on top of a 10⁻³ effect costs. The next generation, with ~10⁷ π⁰π⁰ decays instead of 4×10⁵, agrees with itself to 1.7σ and excludes zero decisively — and note that the settled value sits closer to NA31 than to E731, so the discovery claim was right and the disagreement was still real.

🔑 If you remember only three things

  • It is a part in a thousand of a part in a thousand. ε′ sits on top of an effect that was itself barely visible, and the smallness is the whole difficulty.

  • Two experiments disagreed for a decade and neither was careless. That is what a double ratio’s systematics produce when the signal is this small.

  • The kaon happens to have the structure the question needs. Two amplitudes with different strong phases are required, and not every system offers them.

Where this goes next

CP violation in the decay is now established in the kaon at Re(ε/ε)=1.66×103\mathrm{Re}(\varepsilon'/\varepsilon) = 1.66\times10^{-3} — direct violation being a thousand times smaller again than the indirect violation of §8.5.

The B mesons make it much easier. Their masses open hundreds of channels, so individual branching ratios are 105\sim10^{-5} and asymmetries of several per cent are within reach of 10810^8 pairs — no double ratio required, because the rates of two charge-conjugate decays can simply be compared.

§8.9 applies all of this to charm, the only up-type system available, where CP violation was not seen until 2019 and the observable is a difference of two asymmetries — the double-ratio trick one level up. §8.10 then does the charged B, where there is no mixing at all and the payoff is the unitarity-triangle angle γ\gamma from tree diagrams alone.

Check yourself — CP violation in the decay

0/6 answered · 0 correct

  1. 1.A rate asymmetry between a decay and its CP conjugate requires two phase differences. Why is one not enough?

  2. 2.ε′ ∝ Im A₂ · e^{i(δ₂−δ₀)}. What does each factor require?

  3. 3.Why is the observable a double ratio of four decay rates rather than a single comparison of a decay with its conjugate?

  4. 4.Eq. (8.68) gives Γ(K_S→π⁺π⁻)/Γ(K_S→π⁰π⁰) = 2.55 ± 0.005. What is wrong?

  5. 5.In 1993 NA31 reported (2.30 ± 0.65) × 10⁻³ and E731 reported (−0.74 ± 0.56) × 10⁻³. What should one have concluded?

  6. 6.Why is K_L → π⁰π⁰ so much harder than K_L → π⁺π⁻?

Study aid derived from A. Bettini, Introduction to Elementary Particle Physics, 3rd ed., Cambridge University Press 2024 — published Open Access under CC-BY-NC 4.0, DOI 10.1017/9781009440745. Not the book: an independently written interactive companion, figures redrawn.