Alpha Decay Fundamentals
Alpha decay is the emission of a helium-4 nucleus ($^4_2$He or alpha particle) from an unstable parent nucleus. It is the dominant decay mode for heavy nuclei with $Z > 82$:
The Q-value (kinetic energy released) is:
Alpha decay is energetically possible when $Q_\alpha > 0$. The alpha particle carries most of the kinetic energy due to momentum conservation:
The Geiger-Nuttall Law
Geiger and Nuttall (1911) empirically discovered that alpha-decay half-lives span an enormous range (from nanoseconds to billions of years) but correlate strongly with the decay energy:
where $a$ and $b$ are constants that depend on the parent nucleus charge. A change in $Q_\alpha$ by a factor of 2 can change the half-life by 20 orders of magnitude! This extraordinary sensitivity is explained by quantum tunneling.
Gamow Theory of Alpha Decay
George Gamow (1928) provided the first successful quantum mechanical explanation of alpha decay using the WKB approximation for tunneling through the Coulomb barrier.
The Coulomb Barrier
Inside the nucleus ($r < R$), the alpha particle is bound by the nuclear potential. Outside, it sees the repulsive Coulomb barrier:
The barrier height at the nuclear surface is typically 25-30 MeV, while the alpha particle energy is only 4-9 MeV. Classically, escape is impossible.
WKB Tunneling Probability
The tunneling probability through the barrier using the WKB approximation is:
where $r_{\text{tp}} = 2(Z-2)e^2/(4\pi\epsilon_0 Q_\alpha)$ is the classical turning point and $\mu$ is the reduced mass. Evaluating the integral analytically:
where $\rho = R/r_{\text{tp}}$. The decay rate is then:
where $f$ is the assault frequency (how often the alpha particle hits the barrier).
Python Simulation: Gamow Tunneling
Calculates the Gamow tunneling probability and resulting half-lives as a function of alpha decay energy, demonstrating the Geiger-Nuttall relationship.
Gamow Tunneling Probability
PythonTunneling probability and half-life vs alpha-decay energy for Z=82 daughter
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Fortran Implementation
Calculates alpha-decay half-lives for specific nuclei using the Gamow tunneling formula.
Alpha Decay Half-Life Calculator
FortranComputes Gamow tunneling half-lives for various alpha-emitting nuclei
Click Run to execute the Fortran code
Code will be compiled with gfortran and executed on the server
Derivation: Gamow Theory Step by Step
We present the complete derivation of the Gamow tunneling probability, starting from the WKB approximation applied to the Coulomb barrier.
Step 1: The Potential Model
Model the alpha-nucleus system with a nuclear potential inside $r < R$ and a Coulomb potential outside:
where $B = 2(Z-2)k_e$ with $k_e = e^2/(4\pi\epsilon_0) = 1.44$ MeV$\cdot$fm, and $R$is the nuclear radius at which the alpha particle separates. The Coulomb barrier height at the nuclear surface is:
For $^{238}$U ($Z = 92$), $R \approx 9.3$ fm, giving $V_B \approx 28$ MeV, while $Q_\alpha = 4.27$ MeV. The alpha particle must tunnel through ~24 MeV of barrier.
Step 2: WKB Tunneling Integral
The WKB tunneling probability for a particle of energy $E = Q_\alpha$ through a barrier $V(r)$ is:
where $b = B/E$ is the classical outer turning point ($V(b) = E$) and$\mu = m_\alpha m_D/(m_\alpha + m_D)$ is the reduced mass. Substituting $V(r) = B/r$:
Step 3: Evaluating the Integral
Substituting $r = b\sin^2\theta$, so $dr = 2b\sin\theta\cos\theta\,d\theta$:
where $\sin^2\theta_0 = R/b = \rho$. Evaluating:
Therefore the Gamow factor is:
Step 4: Limiting Cases
For thin barriers ($\rho = R/b \ll 1$, i.e., $E \ll V_B$):
The first term is the Sommerfeld parameter $\pi\eta$ where$\eta = Z_\alpha Z_D e^2/(4\pi\epsilon_0 \hbar v)$:
This shows that $\log T \propto -1/\sqrt{E}$ at low energies, which is the essence of the Geiger-Nuttall relation.
Derivation: Geiger-Nuttall Law
The Geiger-Nuttall law relates the decay constant $\lambda$ to the alpha particle energy. Starting from $\lambda = f \cdot T$ and taking logarithms:
From Gamow Factor to Geiger-Nuttall
In the thin-barrier limit, the Gamow exponent is dominated by the Sommerfeld parameter:
Writing $E = Q_\alpha \cdot (A-4)/A \approx Q_\alpha$ and using $v = \sqrt{2Q_\alpha/\mu}$:
where:
For a given element ($Z_D$ fixed), the log of the half-life is linear in$1/\sqrt{Q_\alpha}$. This explains the extraordinary sensitivity: a factor of 2 change in $Q_\alpha$ corresponds to ~20 orders of magnitude change in half-life for heavy nuclei.
Alpha Decay Systematics
Preformation Factor
The simple Gamow model assumes the alpha particle exists as a preformed cluster inside the nucleus. The decay rate is more accurately written as:
where $P_0$ is the preformation probability. For even-even nuclei,$P_0 \approx 0.01\text{-}0.1$. The preformation factor depends on nuclear structure and can be estimated from the overlap of the parent wave function with the alpha-daughter configuration:
The preformation probability is largest for nuclei just above doubly-magic daughters (e.g., alpha emitters above $^{208}$Pb) and accounts for the favored/unfavored alpha decay classification.
Hindrance Factors
The hindrance factor measures the reduction in decay rate relative to an unhindered transition:
Typical values: HF $\approx 1$ for favored (ground-state to ground-state) decays of even-even nuclei; HF $\approx 4\text{-}10$ for decays involving single-particle state changes; HF $\approx 100\text{-}1000$ for transitions requiring rearrangement of multiple nucleon orbits.
Alpha Spectroscopy and Fine Structure
High-resolution alpha spectroscopy reveals that alpha decay does not always populate the ground state of the daughter nucleus. The alpha spectrum shows multiple discrete lines corresponding to transitions to different excited states.
Fine Structure
For a parent decaying to the $i$-th excited state of the daughter at excitation energy $E_i^*$:
The branching ratio to state $i$ relative to the ground state is:
Since the tunneling probability is extremely sensitive to energy, even small differences in $Q_i$ lead to large changes in branching ratios. The ground-state transition typically dominates unless angular momentum selection rules suppress it.
Angular Momentum Selection Rules
For alpha decay from parent state $J_P^\pi$ to daughter state $J_D^\pi$, the orbital angular momentum $\ell$ carried by the alpha particle must satisfy:
Non-zero $\ell$ adds a centrifugal barrier $\hbar^2 \ell(\ell+1)/(2\mu r^2)$ to the Coulomb barrier, further suppressing the tunneling probability. Each unit of angular momentum reduces the decay rate by roughly an order of magnitude.
Cluster Radioactivity
In 1984, Rose and Jones observed the emission of $^{14}$C from $^{223}$Ra, confirming the prediction of cluster radioactivity by Sandulescu, Poenaru, and Greiner (1980). This is a generalization of alpha decay where heavier clusters are emitted.
Observed Cluster Emissions
Clusters observed in radioactive decay include:
| Parent | Cluster | Daughter | Branch. ratio |
|---|---|---|---|
| $^{223}$Ra | $^{14}$C | $^{209}$Pb | $\sim 10^{-10}$ |
| $^{231}$Pa | $^{24}$Ne | $^{207}$Tl | $\sim 10^{-11}$ |
| $^{236}$Pu | $^{28}$Mg | $^{208}$Pb | $\sim 10^{-14}$ |
| $^{242}$Cm | $^{34}$Si | $^{208}$Pb | $\sim 10^{-16}$ |
Theoretical Framework
The decay rate for cluster emission follows the same Gamow tunneling formalism:
The tunneling probability is much lower due to the larger charge$Z_{\rm cluster} > 2$ and mass:
A striking feature is that nearly all observed cluster emissions produce daughters near$^{208}$Pb (doubly magic), reflecting the enhanced stability and preformation probability associated with the Z=82, N=126 shell closures. The branching ratios relative to alpha decay are extremely small ($10^{-10}$ to $10^{-16}$).