Mechanisms of Nuclear Decay: Alpha, Beta & Gamma Transitions

1. Kinematics and Q-Value of Alpha Disintegration

In spontaneous alpha decay, a parent nucleus $^A_Z\text{X}$ emits an alpha particle ($^4_2\text{He}$ nucleus) and transmutes into a daughter nucleus $^{A-4}_{Z-2}\text{Y}$:

By conservation of relativistic energy in the rest frame of the parent ($M_P c^2 = M_D c^2 + T_D + M_\alpha c^2 + T_\alpha$), the decay energy or $Q_\alpha$-value is:

where $M$ denotes neutral atomic masses (atomic electron binding energies cancel to high precision). Alpha emission is energetically permissible if and only if $Q_\alpha > 0$. Using the SEMF, $Q_\alpha$ becomes positive for mass numbers $A \gtrsim 150$, but observable half-lives ($T_{1/2} < 10^{20}\text{ years}$) require $A \gtrsim 208$ ($Q_\alpha \gtrsim 4 - 9\text{ MeV}$).

2. Kinetic Energy Sharing & Daughter Recoil

By linear momentum conservation in the rest frame of the parent ($p_D = p_\alpha = p$):

Substituting into $Q_\alpha = T_\alpha + T_D$:

For heavy nuclei ($A \sim 200 - 240$), the emitted alpha carries approximately $98\%$ of the total available energy, while the recoiling daughter carries $\approx 2\%$ (typically $80 - 120\text{ keV}$). This recoil energy is sufficient to displace the daughter atom from its crystal lattice, creating radiation damage tracks in minerals.

3. Alpha Spectrum Fine Structure & Hindrance Factors

If the alpha transition populates the ground state of the daughter nucleus, a monoenergetic group of alpha particles is emitted. However, if the daughter nucleus possesses low-lying excited states, the decay branches into multiple discrete energy lines, known as alpha fine structure. The partial decay width to each excited level depends critically on:

• Energy suppression: Lower $Q_i = Q_0 - E^*_i$ drastically thickens the Coulomb barrier, exponentially depressing transition probability.
• Orbital angular momentum: If the parent and daughter states differ in spin $\vec{I}_i \to \vec{I}_f$, the alpha particle must carry orbital angular momentum $\vec{L} = \vec{I}_i - \vec{I}_f$ ($|I_i - I_f| \le L \le I_i + I_f$), adding a centrifugal barrier $V_\ell(r) = \frac{\hbar^2 \ell(\ell+1)}{2\mu r^2}$.
• Hindrance Factor ($HF$): The ratio of the theoretical barrier penetration half-life to the experimental partial half-life. Transitions between spherical states with $\Delta L = 0$ have $HF \approx 1$ (favored decays), whereas transitions requiring intrinsic nucleon rearrangement or high $\ell$ have $HF \sim 10^2 - 10^4$ (hindered decays).

1. The Geiger-Nuttall Empirical Systematics

In 1911, Hans Geiger and John Mitchell Nuttall discovered an astonishing empirical regularity connecting the decay constant $\lambda$ of alpha emitters to the range $R_\alpha$ (or kinetic energy $E_\alpha$) of the emitted particles:

A modest factor of two change in alpha particle energy ($Q_\alpha \approx 4\text{ MeV}$ for $^{238}\text{U}$ to $Q_\alpha \approx 8.8\text{ MeV}$ for $^{212}\text{Po}$) causes the half-life to plunge over 24 orders of magnitude—from $4.47 \times 10^9\text{ years}$ down to $0.3\text{ \mu s}$. Classical mechanics could offer no explanation: the Coulomb barrier height between an alpha particle ($Z_1=2$) and daughter nucleus ($Z_2=Z-2$) at nuclear contact radius $R \approx 8\text{ fm}$ is:

An alpha particle with $E_\alpha \approx 4 - 8\text{ MeV}$ lacks more than $20\text{ MeV}$ of energy required to surmount the barrier classically.

2. George Gamow's Quantum Tunneling Derivation (1928)

George Gamow (and independently Ronald Gurney and Edward Condon) resolved this puzzle by recognizing alpha decay as the quantum mechanical tunneling of a pre-formed alpha particle through the classically forbidden Coulomb potential barrier.

The decay constant $\lambda$ is modeled as the product of two probabilities:

1. Assault Frequency ($f$): The frequency with which the alpha particle inside the nuclear well strikes the barrier. If the particle has internal velocity $v_0 \approx \sqrt{2 E_0 / m_\alpha} \approx 2 \times 10^7\text{ m/s}$ in a nuclear well of diameter $2R \approx 1.5 \times 10^{-14}\text{ m}$:
   $$f = \frac{v_0}{2R} \approx 10^{21} - 10^{22}\text{ s}^{-1}$$
2. Tunneling Transmission Probability ($P$): In the WKB (Wentzel-Kramers-Brillouin) approximation, the transmission through a potential barrier $V(r)$ from the nuclear surface $R$ to the classical turning point $b$ ($V(b) = Q_\alpha$) is:
   $$P = \exp\left( -2 \int_{R}^{b} k(r) dr \right) = \exp\left( -2 \int_{R}^{b} \sqrt{\frac{2\mu}{\hbar^2} [V_C(r) - Q_\alpha]} dr \right)$$

   where $\mu = \frac{M_D M_\alpha}{M_D + M_\alpha}$ is the reduced mass and $V_C(r) = \frac{2(Z - 2)e^2}{4\pi\varepsilon_0 r}$. At the turning point, $b = \frac{2(Z - 2)e^2}{4\pi\varepsilon_0 Q_\alpha}$.

3. Analytical Evaluation of the Gamow Factor

Evaluating the WKB integral by setting $r = b \cos^2\theta$ yields:

For heavy nuclei, $b \gg R$, so $\arccos(\sqrt{R/b}) \approx \frac{\pi}{2} - \sqrt{R/b}$. Thus:

This reproduces the Geiger-Nuttall law from first principles: the exponent depends inversely on $\sqrt{Q_\alpha}$, explaining how minor variations in alpha energy produce astronomical differences in nuclear half-lives.

1. The Three Modes of Nuclear Beta Decay

Beta decay is a weak interaction process in which an isobaric nucleon changes its isospin state ($\Delta Z = \pm 1$, constant $A$):

1. Beta-Minus ($\beta^-$) Decay: Occurs in neutron-rich nuclei. A bound neutron transforms into a proton, emitting an electron ($e^-$) and an electron antineutrino ($\bar{\nu}_e$):
   $$^A_Z\text{X} \longrightarrow {^A_{Z+1}\text{Y}} + e^- + \bar{\nu}_e \quad (n \to p + e^- + \bar{\nu}_e)$$
   $$Q_{\beta^-} = [M(A, Z) - M(A, Z+1)] c^2$$
2. Beta-Plus ($\beta^+$) Decay (Positron Emission): Occurs in proton-rich nuclei. A bound proton transforms into a neutron, emitting a positron ($e^+$) and an electron neutrino ($\nu_e$):
   $$^A_Z\text{X} \longrightarrow {^A_{Z-1}\text{Y}} + e^+ + \nu_e \quad (p \to n + e^+ + \nu_e)$$

   Accounting for the two electron rest masses in neutral atomic mass accounting ($M_P - M_D - 2m_e$):

   $$Q_{\beta^+} = [M(A, Z) - M(A, Z-1) - 2m_e] c^2 = [M(A, Z) - M(A, Z-1)] c^2 - 1.022\text{ MeV}$$

   Positron decay requires a minimum mass threshold difference of $2 m_e c^2 = 1.022\text{ MeV}$.

3. Orbital Electron Capture ($EC$): The nucleus captures an inner atomic orbital electron (typically $K$-shell):
   $$^A_Z\text{X} + e^- \longrightarrow {^A_{Z-1}\text{Y}} + \nu_e \quad (p + e^- \to n + \nu_e)$$
   $$Q_{EC} = [M(A, Z) - M(A, Z-1)] c^2 - B_n$$

   where $B_n$ is the atomic binding energy of the captured electron. Electron capture competes with $\beta^+$ and is the sole decay mode when $[M(A,Z) - M(A,Z-1)]c^2 < 1.022\text{ MeV}$.

2. The Continuous Beta Spectrum Crisis & Pauli's Neutrino Hypothesis

In 1914, James Chadwick demonstrated that whereas alpha particles are emitted with sharp discrete energies, beta electrons exhibit a continuous kinetic energy distribution extending from zero up to a well-defined maximum end-point energy $E_{\text{max}} = Q_\beta$. If beta decay were a two-body transition ($^A_Z\text{X} \to {^A_{Z+1}\text{Y}} + e^-$), conservation of energy and linear momentum would require the electron to possess a unique discrete energy:

The continuous spectrum, combined with apparent violations of angular momentum (e.g., $^{14}_6\text{C}(0^+) \to {^{14}_7\text{N}}(1^+) + e^-(1/2)$ has non-conserved half-integer spin), led Niels Bohr to suggest that energy might only be conserved statistically. To preserve strict conservation laws, Wolfgang Pauli (1930) proposed his "desperate remedy": an undetectable, neutral, spin-$1/2$ fermion with vanishingly small rest mass emitted simultaneously with the electron. Enrico Fermi named this elusive particle the neutrino ($\nu$). In three-body decay:

The three-body phase space partitions energy continuously between $T_e$ and $T_\nu$, exactly resolving the anomaly.

1. Fermi's Quantum Formulation of Beta Disintegration

In 1934, Enrico Fermi formulated the quantum theory of beta decay using time-dependent perturbation theory (Fermi's Golden Rule #2). The transition rate per unit energy interval is:

where $V_{fi} = g \int [\psi_f^* \phi_e^*(\vec{r}) \phi_\nu^*(\vec{r})] \mathcal{O}_{\text{weak}} \psi_i d^3r$, $g \approx 1.4 \times 10^{-62}\text{ J}\cdot\text{m}^3$ ($G_F \approx 1.166 \times 10^{-5}\text{ GeV}^{-2}$) is Fermi's weak coupling constant, and $\rho(E_0)$ is the density of accessible two-particle continuum states.

Because the electron de Broglie wavelength ($\lambda_e \sim 1000\text{ fm}$) is vastly larger than the nuclear radius ($R \sim 5\text{ fm}$), the lepton wave functions can be approximated as plane waves and expanded as $e^{i \vec{k}\cdot\vec{r}} \approx 1 + i \vec{k}\cdot\vec{r} + \dots$. Retaining the leading term ($\ell = 0$) defines allowed transitions.

2. Theoretical Electron Energy Spectrum & The Fermi Function

The statistical phase space factor combined with the Coulomb correction of the daughter nucleus yields the differential electron momentum distribution:

where $E_0 = Q_\beta + m_e c^2$ is the total endpoint energy ($E_e^2 = p_e^2 c^2 + m_e^2 c^4$), and $F(Z', p_e)$ is the relativistic Fermi Function representing Coulomb distortion of the electron wave function by the daughter nucleus ($Z' = +Z_D$ for $e^-$, $Z' = -Z_D$ for $e^+$):

3. The Fermi-Kurie Plot

Linearizing the spectral distribution provides a sensitive test of the theory and precise measurement of the decay endpoint:

Plotting $K(p_e)$ against electron total energy $E_e$ yields a perfect straight line for allowed transitions whose intercept on the horizontal axis determines $E_0$. A non-zero neutrino mass $m_\nu$ would produce a vertical downturn with infinite slope at the extreme endpoint $E_e = E_0 - m_\nu c^2$. Current tritium beta decay experiments (KATRIN) establish an upper bound $m_\nu < 0.45\text{ eV}/c^2$.

4. Classification & Selection Rules: Fermi vs Gamow-Teller

In allowed transitions ($\ell = 0$), the emitted leptons carry zero orbital angular momentum. Their intrinsic spins ($s_e = 1/2, s_\nu = 1/2$) can couple in two distinct orientations:

Superallowed $0^+ \to 0^+$ pure Fermi decays (e.g., $^{14}\text{O} \to {^{14m}\text{N}}$) have matrix elements $|M_F|^2 = 2$ governed strictly by isospin symmetry, allowing high-precision determination of the vector coupling constant $G_V$ and the Cabibbo-Kobayashi-Maskawa matrix element $V_{ud}$.

1. Electromagnetic Multipole Radiations

Following alpha or beta decay, the daughter nucleus is typically left in an excited quantum state. It de-excites to the ground state by emitting a gamma-ray photon ($\gamma$). Photons are spin-$1$ bosons; because a photon has intrinsic spin $1$, monoenergetic transitions between two spin-zero states ($0^+ \to 0^+$) via single photon emission are strictly forbidden ($\gamma$ carries at least $L = 1\hbar$ angular momentum).

Gamma transitions are classified by the multipole order $L$ (dipole $L=1$, quadrupole $L=2$, octupole $L=3$) and electromagnetic character (Electric $EL$ or Magnetic $ML$):

The parity selection rules dictate:

2. Weisskopf Single-Particle Estimates

Victor Weisskopf derived standard reference single-particle transition rates $\lambda(EL)$ and $\lambda(ML)$ assuming a single proton transitions between single-particle shell model orbitals within a sphere of radius $R = R_0 A^{1/3}$ ($E_\gamma$ in MeV, $A$ mass number):

Key physical conclusions from Weisskopf rates:

• For a given multipole order $L$, electric transitions are typically two orders of magnitude faster than magnetic: $\lambda(EL) / \lambda(ML) \sim 100$.
• Each unit increase in multipole order $L$ suppresses the decay rate by a factor of roughly $10^5 - 10^6$ due to the small nuclear size parameter $(k R)^2 \approx (E_\gamma R / \hbar c)^2 \ll 1$.
• Consequently, transitions proceed predominantly via the lowest allowed multipole order ($L_{\text{min}} = |I_i - I_f|$), with $M1/E2$ mixing commonly observed when $L=1$ and $L=2$ compete.

3. Nuclear Isomerism & Metastable States

When an excited nuclear state requires a high multipole transition ($L \ge 3$ or $4$) combined with a low transition energy ($E_\gamma \lesssim 100\text{ keV}$), the transition probability becomes exceedingly small. The excited state exhibits a remarkably long lifetime (seconds, hours, or even years), termed a nuclear isomer (denoted with an 'm', e.g., $^{99m}_{43}\text{Tc}$ with $T_{1/2} = 6.01\text{ h}$, decaying via an $M4$ transition to the ground state). The longest known isomer is $^{180m}_{73}\text{Ta}$ ($9^-$ state), with a half-life exceeding $10^{15}\text{ years}$, exceeding the age of the universe.

1. Internal Conversion (IC)

Internal conversion is an electromagnetic de-excitation mechanism that competes directly with gamma-ray photon emission. Instead of emitting a photon, the excited nucleus interacts directly via the near-field Coulomb interaction with an atomic inner-shell electron ($K, L, M$), ejecting the electron into the continuum:

where $E^*$ is the nuclear excitation energy and $B_e$ is the electron atomic binding energy. Unlike beta decay electrons, internal conversion electrons are strictly monoenergetic. Vacancies created in inner atomic shells subsequently trigger the emission of characteristic X-rays or Auger electrons.

The internal conversion coefficient ($\alpha$) is defined as the branching ratio:

Internal conversion dominates under three conditions:

• High atomic number: $\alpha \propto Z^3$, because inner atomic electrons spend greater time inside the nucleus.
• Low transition energy: $\alpha \propto E_\gamma^{-(L + 5/2)}$.
• High multipolarity: $\alpha$ increases dramatically with multipole order $L$. For $0^+ \to 0^+$ transitions (e.g., $^{16}\text{O}^*$ at $6.05\text{ MeV}$, $^{72}\text{Ge}$), single photon emission is forbidden ($\lambda_\gamma = 0$), forcing decay to occur 100% via internal conversion ($\alpha = \infty$) or $e^+e^-$ pair conversion.

2. Four Primary Photon Interactions with Matter

As gamma-ray photons traverse matter, they do not lose energy continuously; instead, they undergo catastrophic single-interaction scattering or absorption events characterized by a linear attenuation coefficient $\mu(E_\gamma)$:

1. Photoelectric Absorption ($\sigma_{\text{pe}} \propto Z^{4-5} / E_\gamma^{3.5}$): The photon is completely absorbed by a bound atomic electron, which is ejected with kinetic energy $T_e = E_\gamma - B_K$. Dominates at low energies ($E_\gamma < 0.1\text{ MeV}$) and in high-$Z$ absorbers (e.g., Lead, $Z=82$).
2. Compton Scattering ($\sigma_C \propto Z / E_\gamma$): Inelastic scattering of the photon from a quasi-free electron. By energy-momentum conservation, the scattered photon energy is given by the Compton formula:
   $$E_\gamma' = \frac{E_\gamma}{1 + \frac{E_\gamma}{m_e c^2}(1 - \cos\theta)}$$
   The scattered electron recoils with kinetic energy $T_e = E_\gamma - E_\gamma'$, reaching its maximum at backscattering ($\theta = 180^\circ$), defining the sharp Compton edge in gamma spectroscopy:
   $$T_{\text{max}} = E_\gamma \left[ \frac{2 E_\gamma / m_e c^2}{1 + 2 E_\gamma / m_e c^2} \right]$$
   Compton scattering dominates in the intermediate energy regime ($0.5\text{ MeV} < E_\gamma < 5\text{ MeV}$).
3. Pair Production ($\sigma_{\text{pp}} \propto Z^2 \ln(E_\gamma)$): In the strong Coulomb field of an atomic nucleus, a high-energy photon transforms into an electron-positron pair: $\gamma \to e^- + e^+$. Requires a strict threshold energy $E_{\text{threshold}} = 2 m_e c^2 = 1.022\text{ MeV}$. Dominates at high energies ($E_\gamma > 5 - 10\text{ MeV}$). Subsequent positron annihilation produces two collinear $511\text{ keV}$ annihilation photons.
4. Photonuclear Reactions ($\gamma, n$): For photon energies exceeding the nuclear neutron separation energy ($E_\gamma > S_n \sim 7 - 10\text{ MeV}$), the photon excites the Giant Dipole Resonance (GDR), ejecting a neutron.