Case Study 21.2: Fermi's Golden Rule in Nuclear and Particle Physics

How a formula derived for atomic transitions became the central computational tool of nuclear physics and high-energy particle physics

From Atoms to Nuclei: The Universality of Fermi's Golden Rule

Fermi's golden rule, $\Gamma = (2\pi/\hbar)|V_{fi}|^2\rho(E_f)$, was derived in the context of atomic transitions — an electron changing energy levels in a hydrogen atom. But the formula's power lies in its universality: it applies to any quantum system making a transition from one state to another under the influence of a weak perturbation, provided there is a continuum (or near-continuum) of final states.

In nuclear and particle physics, the "perturbation" is the fundamental interaction (strong, electromagnetic, or weak force), the "initial state" is the unstable particle or nucleus, and the "final states" are the decay products or scattering final states. The same three ingredients — matrix element, density of states, and $2\pi/\hbar$ — determine every decay rate and scattering cross section.

This case study traces how Fermi's golden rule has been applied to three landmark problems in nuclear and particle physics, each illustrating a different aspect of the formula.

Application 1: Beta Decay and Fermi's Theory of the Weak Interaction

The Problem

In 1934, Enrico Fermi proposed the first theory of beta decay — the process in which a neutron inside a nucleus transforms into a proton, emitting an electron and an antineutrino:

$$n \to p + e^- + \bar{\nu}_e$$

This was, in modern terms, a theory of the weak nuclear force. Fermi treated the interaction as a point-like "four-fermion" coupling with a coupling constant $G_F$ (now called the Fermi constant):

$$\hat{V} = \frac{G_F}{\sqrt{2}}\,(\bar{\psi}_p \gamma^\mu \psi_n)(\bar{\psi}_e \gamma_\mu \psi_\nu) + \text{h.c.}$$

(We write this in the modern notation for clarity; Fermi's original formulation did not use Dirac spinors.) The matrix element is:

$$V_{fi} = \frac{G_F}{\sqrt{2}}\,\langle p, e, \bar{\nu}|\hat{V}|n\rangle$$

Applying Fermi's Golden Rule

The transition rate is:

$$\Gamma = \frac{2\pi}{\hbar}\,|V_{fi}|^2\,\rho(E_f)$$

Here $\rho(E_f)$ is the density of final states for a three-body decay (proton + electron + antineutrino). The proton's recoil kinetic energy is negligible (it is 1,836 times heavier than the electron), so the phase space is essentially that of the electron and antineutrino sharing the available energy $Q = (m_n - m_p)c^2 = 1.293$ MeV.

The density of final states for the electron (with momentum between $p_e$ and $p_e + dp_e$) and antineutrino is:

$$\rho(E_f)\,dE_f = \frac{4\pi p_e^2\,dp_e}{(2\pi\hbar)^3}\,\frac{4\pi p_\nu^2}{(2\pi\hbar)^3}\,\frac{V^2}{dp_\nu/dE_\nu}$$

After integration over the antineutrino momentum (fixed by energy conservation: $E_\nu = Q - T_e$, where $T_e$ is the electron kinetic energy), Fermi obtained the electron energy spectrum:

$$\frac{d\Gamma}{dT_e} \propto |V_{fi}|^2\, p_e\, E_e\, (Q - T_e)^2\, F(Z, E_e)$$

where $F(Z, E_e)$ is the Fermi function accounting for the Coulomb attraction between the emitted electron and the daughter nucleus. The total decay rate is obtained by integrating over $T_e$ from 0 to $Q$:

$$\Gamma = \frac{G_F^2\,|M_{fi}|^2}{2\pi^3\hbar^7 c^6}\, f(Q, Z)$$

where $f(Q, Z)$ is the dimensionless Fermi integral (a function of the $Q$-value and the nuclear charge $Z$) and $|M_{fi}|^2$ is the nuclear matrix element.

The $ft$-Value and Nuclear Structure

The combination $ft_{1/2}$, where $t_{1/2} = \ln 2/\Gamma$ is the half-life and $f$ is the Fermi integral, is a measure of the nuclear matrix element alone:

$$ft_{1/2} = \frac{\ln 2 \cdot 2\pi^3\hbar^7 c^6}{G_F^2\,|M_{fi}|^2}$$

For superallowed Fermi transitions ($0^+ \to 0^+$), the nuclear matrix element is fixed by isospin symmetry: $|M_{fi}|^2 = 2$ (for $T = 1 \to T = 1$ transitions). The measured $ft$-values of 14 superallowed beta decays are all consistent with $ft = 3072.27 \pm 0.72$ s, providing the most precise determination of the Fermi constant:

$$G_F = (1.16638 \pm 0.00001) \times 10^{-5}\,\text{GeV}^{-2}$$

This is Fermi's golden rule at its most powerful: the same formula that describes an atom absorbing a photon also determines the fundamental coupling constant of the weak force.

📊 By the Numbers: The free neutron decay has $t_{1/2} = 611.6 \pm 0.6$ s (about 10 minutes). The $Q$-value is $1.293$ MeV. The $ft$-value is approximately 1,130 s — a "mixed" transition (part Fermi, part Gamow-Teller) with $|M_F|^2 + 3|M_{GT}|^2 = |1|^2 + 3|g_A|^2$, where $g_A \approx 1.276$ is the axial coupling constant.

Application 2: Nuclear Gamma Decay and the Weisskopf Estimates

Single-Particle Transition Rates

When a nucleus transitions from an excited state to a lower state by emitting a gamma ray, the physics is almost identical to atomic photon emission — the same electric dipole, magnetic dipole, and electric quadrupole matrix elements, the same selection rules, the same $\omega^3$ (or $\omega^5$ for E2) scaling. The differences are:

1. The "particles" are protons and neutrons rather than electrons.
2. The relevant length scale is the nuclear radius $R \approx 1.2\,A^{1/3}$ fm rather than $a_0 = 0.53 \times 10^5$ fm.
3. The photon energies are MeV rather than eV, so the selection rules involve nuclear spins and parities.

Victor Weisskopf derived order-of-magnitude estimates for nuclear gamma transition rates using a single-particle model. Applying Fermi's golden rule with the dipole matrix element $|\vec{d}| \sim eR$ (a single proton moving across the nucleus):

Electric dipole (E1): $$\Gamma_{\text{E1}}^W = \frac{4\alpha}{3}\left(\frac{E_\gamma}{\hbar c}\right)^3 \left(\frac{3}{3+l}\right)^2 R^2\, c \approx 1.0 \times 10^{14}\, A^{2/3}\, E_\gamma^3\,\text{s}^{-1}$$

(with $E_\gamma$ in MeV). For a typical $E_\gamma = 1$ MeV transition in a medium-mass nucleus ($A = 100$):

$$\Gamma_{\text{E1}}^W \approx 2.2 \times 10^{15}\,\text{s}^{-1}, \qquad \tau \approx 5 \times 10^{-16}\,\text{s}$$

Electric quadrupole (E2): $$\Gamma_{\text{E2}}^W \approx 7.3 \times 10^{7}\, A^{4/3}\, E_\gamma^5\,\text{s}^{-1}$$

For $E_\gamma = 1$ MeV, $A = 100$: $\tau \approx 6 \times 10^{-12}$ s.

Magnetic dipole (M1): $$\Gamma_{\text{M1}}^W \approx 3.1 \times 10^{13}\, E_\gamma^3\,\text{s}^{-1}$$

For $E_\gamma = 1$ MeV: $\tau \approx 3 \times 10^{-14}$ s.

Selection Rules for Nuclear Transitions

The same angular momentum and parity selection rules from atomic physics apply to nuclear transitions, with nuclear spin $I$ replacing atomic angular momentum $J$:

Each step up in multipole order suppresses the rate by a factor of roughly $(kR)^2 \sim (E_\gamma R/\hbar c)^2 \sim 10^{-7}$ for $E_\gamma = 1$ MeV and $R = 6$ fm. This factor is much smaller than the atomic analogue $(a_0/\lambda)^2 \sim 10^{-7}$ — it happens to be about the same numerically, because the nuclear radius is to MeV photons what the atomic radius is to eV photons.

Nuclear Isomers: When Selection Rules Create Long-Lived States

When the only available decay path requires a high multipole order (e.g., $\Delta I = 4$ requiring E4 or M4), the decay rate is suppressed by many orders of magnitude, creating nuclear isomers — excited nuclear states with lifetimes ranging from nanoseconds to years.

The most famous isomer is ${}^{180m}\text{Ta}$, an excited state of tantalum-180 with spin $I^\pi = 9^-$. Its ground state has $I^\pi = 1^+$, requiring a $\Delta I = 8$ transition — which would need an M8 or E8 multipole. The suppression is so extreme that ${}^{180m}\text{Ta}$ has a measured lower bound on its lifetime of $> 10^{15}$ years — longer than the age of the universe. It is the only naturally occurring nuclear isomer.

This is Fermi's golden rule working in reverse: when the matrix element is tiny (high-multipole transitions), the transition rate is correspondingly tiny, and "unstable" states become effectively stable.

Application 3: Particle Decay Rates and the CKM Matrix

Weak Decays of Hadrons

In modern particle physics, the weak decays of quarks and leptons are described by Fermi's golden rule with the perturbation provided by the $W$ boson exchange:

$$\hat{V} = \frac{g^2}{8M_W^2}\,\bar{u}_f\,\gamma^\mu(1-\gamma^5)\,u_i \cdot \bar{u}_\ell\,\gamma_\mu(1-\gamma^5)\,v_\nu$$

where $g$ is the weak coupling constant and $M_W = 80.4$ GeV/$c^2$ is the $W$ boson mass. At low energies ($E \ll M_W c^2$), this reduces to Fermi's four-fermion interaction with $G_F/\sqrt{2} = g^2/(8M_W^2)$.

The matrix element for a quark transition $u_i \to u_j$ (where $i, j$ label quark flavors) contains the Cabibbo-Kobayashi-Maskawa (CKM) matrix element $V_{ij}$:

$$|V_{fi}|^2 \propto |V_{ij}|^2$$

The CKM matrix encodes the mixing between quark mass eigenstates and weak interaction eigenstates:

$$V_{\text{CKM}} = \begin{pmatrix} V_{ud} & V_{us} & V_{ub} \\ V_{cd} & V_{cs} & V_{cb} \\ V_{td} & V_{ts} & V_{tb} \end{pmatrix} \approx \begin{pmatrix} 0.974 & 0.227 & 0.004 \\ 0.226 & 0.973 & 0.041 \\ 0.008 & 0.040 & 0.999 \end{pmatrix}$$

The off-diagonal elements are small, explaining why flavor-changing transitions are suppressed. This is directly measurable through Fermi's golden rule: the decay rate of a particle involving a $u_i \to u_j$ transition is proportional to $|V_{ij}|^2$.

Example: Pion Decay

The charged pion $\pi^+$ (a bound state of $u$ and $\bar{d}$ quarks) decays via the weak interaction:

$$\pi^+ \to \mu^+ + \nu_\mu$$

Fermi's golden rule gives the decay rate:

$$\Gamma = \frac{G_F^2}{8\pi}\,|V_{ud}|^2\, f_\pi^2\, m_\mu^2\left(1 - \frac{m_\mu^2}{m_\pi^2}\right)^2 m_\pi$$

where $f_\pi \approx 130$ MeV is the pion decay constant (encoding the quark-level matrix element) and the factor $m_\mu^2$ reflects the helicity suppression of the decay (the outgoing lepton must be in the "wrong" helicity state, which costs a factor of $m_\ell^2$).

This gives $\tau_{\pi^+} = 1/\Gamma = 2.60 \times 10^{-8}$ s, in excellent agreement with the measured value of $2.6033 \times 10^{-8}$ s.

The helicity suppression predicts that the ratio of $\pi \to e\nu$ to $\pi \to \mu\nu$ rates should be:

$$\frac{\Gamma(\pi \to e\nu)}{\Gamma(\pi \to \mu\nu)} = \frac{m_e^2}{m_\mu^2}\,\frac{(m_\pi^2 - m_e^2)^2}{(m_\pi^2 - m_\mu^2)^2} = 1.28 \times 10^{-4}$$

The experimental value is $(1.230 \pm 0.004) \times 10^{-4}$, a triumph of the Standard Model. The small discrepancy is understood as a radiative correction — a higher-order effect in the perturbative expansion that Fermi's golden rule naturally invites us to compute.

Example: $B$-Meson Mixing

One of the most spectacular applications of Fermi's golden rule in modern particle physics is $B^0$-$\bar{B}^0$ mixing. The $B^0$ meson (a $\bar{b}d$ bound state) can oscillate into its antiparticle $\bar{B}^0$ ($b\bar{d}$) through a second-order weak process involving virtual $W$ bosons and top quarks.

The mixing amplitude is given by a "box diagram" — a second-order application of Fermi's golden rule, with the sum over intermediate states including virtual top quarks:

$$\langle \bar{B}^0|\hat{H}_{\text{eff}}|B^0\rangle \propto G_F^2 M_W^2 (V_{tb}^* V_{td})^2\, f_{B}^2\, m_B\, S(m_t^2/M_W^2)$$

where $S$ is the Inami-Lim function (the loop integral over virtual top-quark momenta). The oscillation frequency is:

$$\Delta m_d = \frac{G_F^2 M_W^2}{6\pi^2}\, |V_{tb}^* V_{td}|^2\, \eta_B\, S(x_t)\, f_B^2\, B_B\, m_B$$

The measured $\Delta m_d = 3.334 \times 10^{-13}$ GeV provides a determination of $|V_{td}| \approx 0.008$, one of the least-well-known CKM elements. The entire analysis rests on Fermi's golden rule — extended to second order in the weak interaction, but the same fundamental formula.

The Deeper Lesson: Why One Formula Rules Them All

Fermi's golden rule is not a formula specific to atomic physics. It is a direct consequence of the structure of quantum mechanics itself:

1. Perturbation theory works because interactions are weak. The electromagnetic coupling $\alpha \approx 1/137$ ensures that QED perturbation theory converges rapidly. The Fermi constant $G_F \approx 10^{-5}\,\text{GeV}^{-2}$ ensures that weak decays are amenable to lowest-order calculations.

2. Phase space (density of states) is universal. Whether the final state is a photon, an electron-antineutrino pair, or a quark-antiquark jet, the phase space is determined by the same relativistic kinematics.

3. The matrix element encodes all the physics. Different interactions (electromagnetic, weak, strong) enter through different matrix elements $V_{fi}$, but the framework for computing the rate is always the same.

From hydrogen spectral lines to gravitational wave detectors, from nuclear reactors to the Large Hadron Collider, Fermi's golden rule is the formula that connects quantum amplitudes to measurable rates. It earned its name.

"Golden Rule No. 2" — Fermi's 1950 lecture notes. (Golden Rule No. 1, concerning the Born approximation, is essentially the same formula applied to scattering — Chapter 22.)