Radioactive Decay

Radioactive decay (also known as nuclear decay, radioactivity, or nuclear disintegration) is the spontaneous, stochastic process by which an unstable atomic nucleus (parent radionuclide) loses excess nuclear binding energy by emitting ionizing particles or electromagnetic radiation, transforming into a more stable daughter nuclide (or a lower-energy nuclear isomer). Discovered in uranium salts by Henri Becquerel in 1896, extended to polonium and radium by Marie Skłodowska-Curie and Pierre Curie in 1898, and formulated quantitatively as a first-order statistical transmutation law by Ernest Rutherford and Frederick Soddy in 1902, radioactive decay is governed jointly by the strong nuclear force, electromagnetic Coulomb repulsion, and the weak nuclear interaction.

Within nuclear physics, isotope geochemistry, and medical physics, radioactive decay provides both a window into quantum-mechanical nuclear structure and an invariant physical chronometer. Because nuclear binding energies (mega-electronvolts, MeV) exceed chemical valence bond energies (electronvolts, eV) by six orders of magnitude, the decay constant λ of nearly all alpha, beta, and gamma transitions is immune to ambient temperature, pressure, and chemical bonding, enabling radiometric geochronology across geological and cosmological timescales via characteristic half-lives.

First-Order Rutherford–Soddy Poisson Law and Activity Units (Bq vs. Ci)

At the level of an individual unstable nucleus, radioactive disintegration is a memoryless Poisson quantum process: the probability dP that a given nucleus decays within an infinitesimal time interval dt is constant and equal to λ dt, where λ (in s⁻¹) is the isotope-specific decay constant, independent of the age of the nucleus. For a macroscopic ensemble containing N(t) undecayed parent nuclei at time t, the mean rate of disintegration obeys the first-order linear differential equation:

The radioactive activity (A) of a sample is defined as the number of nuclear disintegrations per unit time:

In the International System of Units (SI), activity is measured in the becquerel (Bq), defined as 1 Bq = 1 disintegration s⁻¹, whereas the historical non-SI unit, the curie (Ci) (originally the activity of 1 gram of ²²⁶Ra), is fixed at 1 Ci = 3.7 × 10¹⁰ Bq. The inverse of the decay constant is the mean lifetime τ = 1 / λ, while the time required for N(t) or A(t) to decay to half its initial value is the half-life t_1/2 = ln(2) / λ = τ ln(2) ≈ 0.693147 τ. When a single parent nuclide can decay via k competing branching channels (e.g., ⁴⁰K decaying via β⁻ to ⁴⁰Ca with branching fraction 89.28% and via electron capture to ⁴⁰Ar with 10.72%), the partial decay constants add linearly: λ_tot = Σ_i λ_i.

Relativistic Mass Defect and Reaction Q-Value Conservation

Every spontaneous radioactive decay must be exothermic (Q > 0), converting rest-mass defect Δ*m* between the parent nucleus and the decay products into the total relativistic translational kinetic energy and photon energy of the products according to Albert Einstein's mass–energy equivalence:

where atomic masses are conventionally expressed in unified atomic mass units (1 u = 931.49410242 MeV/c²).

Alpha Decay, Gamow Quantum Tunneling, and the Geiger–Nuttall Law

Alpha (α) decay occurs predominantly in heavy nuclei beyond lead (Z > 82, such as ²³⁸U, ²²⁶Ra, ²³⁹Pu, and ²¹⁰Po), where short-range attractive strong nucleon-nucleon forces (saturating as A) are increasingly offset by long-range electrostatic Coulomb repulsion among protons (scaling as Z(Z − 1)/A^1/3). Because the ⁴₂He alpha particle is a doubly magic nucleus (Z = 2, N = 2) with an exceptionally high binding energy of 28.296 MeV (7.07 MeV per nucleon), pre-clustering an α-particle inside a heavy parent nucleus releases positive Q_α (typically 4 to 9 MeV).By conservation of linear momentum, the emitted α-particle carries away a discrete translational kinetic energy K_α = Q_α (A − 4)/A (~98% of Q_α), while the heavy daughter nucleus recoils with K_recoil = Q_α (4/A).

Classically, a 4.27 MeV α-particle inside ²³⁸U cannot escape because the surrounding electrostatic Coulomb barrier peaks at V_C ≈ 26–30 MeV at the nuclear radius R ≈ 9 fm. In 1928, George Gamow (and independently Ronald Gurney and Edward Condon) solved this paradox as the first triumph of quantum mechanics applied to the nucleus: the α-particle's wave function penetrates the classically forbidden barrier (r_in < r < r_out) via quantum tunneling with WKB transmission probability T ≈ exp(−2*G*), where the Gamow factor G scales as Z_d Q_α^−1/2:

This exponential WKB tunneling sensitivity T ≈ exp(−2*G*) explains the Geiger–Nuttall law (1911): doubling Q_α from 4.01 MeV in ²³²Th to 8.95 MeV in ²¹²Po causes the alpha half-life to plummet by 24 orders of magnitude, from 1.40 × 10¹⁰ years down to 0.299 microseconds.

Weak-Interaction Beta Decay (β⁻, β⁺, EC) and Neutrino Conservation

Nuclides that lie off the valley of beta stability (where the neutron-to-proton ratio N/Z increases from 1.0 in light Z ≤ 20 nuclei to ~1.54 in ²⁰⁸Pb to counterbalance proton Coulomb repulsion) undergo isobaric beta transitions (constant mass number A) mediated by the charged-current weak nuclear interaction (W^± vector bosons) formulated by Enrico Fermi in 1934:

• Beta-Minus (β⁻) Decay (Neutron-Rich Nuclei, e.g., ¹⁴C → ¹⁴N, ³H → ³He, ⁶⁰Co → ⁶⁰Ni): A down quark inside a bound neutron converts into an up quark (d → u + W⁻ → u + e⁻ + ν̄_e), increasing atomic number by one (Z → Z + 1).
• Beta-Plus (β⁺) Positron Decay (Proton-Rich Nuclei, e.g., ¹⁸F → ¹⁸O, ¹¹C → ¹¹B): An up quark inside a bound proton converts into a down quark (u → d + W⁺ → d + e⁺ + ν_e), decreasing atomic number by one (Z → Z − 1). Because the neutral atomic mass balance loses two electron masses, β⁺ decay is energetically possible only if the neutral atomic mass difference exceeds 2*m*_ec² = 1.022 MeV.
• Orbital Electron Capture (EC, e.g., ⁴⁰K → ⁴⁰Ar, ⁷Be → ⁷Li): A proton captures an inner-shell (90% K-shell 1*s*) atomic electron (p + e⁻ → n + ν_e), emitting a monoenergetic neutrino and leaving a K-shell vacancy that relaxes via characteristic daughter X-rays or Auger electrons.

Because β⁻ and β⁺ decays are three-body disintegrations (daughter nucleus + e^± + neutrino), the emitted beta particles exhibit a continuous energy spectrum from 0 up to the endpoint energy Q_β—a spectrum that led Wolfgang Pauli in 1930 to postulate the electrically neutral neutrino to rescue conservation of energy, linear momentum, and spin-1/2 angular momentum. Transitions are classified by nuclear spin and parity change into Fermi (Δ*I* = 0, spins antiparallel) and Gamow–Teller (Δ*I* = 0, ±1, spins parallel) allowed vs. forbidden transitions.

Gamma De-Excitation, Isomeric States, and Internal Conversion

Following α or β decay, the daughter nucleus is frequently left in an excited nuclear shell state (E^* ~ 10 keV to 5 MeV). De-excitation to the ground state occurs within picoseconds via emission of a gamma-ray photon (γ) of multipolarity Eλ or Mλ, or—when a large spin difference (Δ*I* ≥ 3–4) delays electromagnetic transition for seconds to hours—via a metastable nuclear isomer such as technetium-99m (^99mTc, t_1/2 = 6.0067 h, E_γ = 140.5 keV), the most widely used diagnostic radiopharmaceutical in single-photon emission computed tomography (SPECT). Alternatively, the excited nucleus can transfer its excitation energy directly to a K- or L-shell orbital electron via internal conversion (IC), ejecting a monoenergetic conversion electron of energy K_e = E_γ − B_K/L.

Secular vs. Transient Equilibrium and Isochron Radiometric Dating

Heavy actinide nuclei decay through multi-step natural decay series (A = 4*n* Thorium series from ²³²Th → ²⁰⁸Pb; A = 4*n* + 2 Uranium/Radium series from ²³⁸U → ²⁰⁶Pb; A = 4*n* + 3 Actinium series from ²³⁵U → ²⁰⁷Pb). For a parent N₁ (decay constant λ₁) decaying to a radioactive daughter N₂ (decay constant λ₂), the Bateman equation (Harry Bateman, 1910) gives:

When the parent half-life is vastly longer than the daughter half-life (λ₁ << λ₂, e.g., ²³⁸U [t_1/2 = 4.468 Gyr] → ²³⁴Th [t_1/2 = 24.1 d], or ²²⁶Ra [1,600 yr] → ²²²Rn [3.82 d]), after ~7 daughter half-lives the daughter's activity matches the parent's activity: A₂ ≈ A₁ (secular equilibrium, where N₂/N₁ = λ₁/λ₂ = t_1/2,2/t_1/2,1). When λ₁ < λ₂ are comparable (such as the medical generator ⁹⁹Mo [66.0 h] → ^99mTc [6.01 h]), the system reaches transient equilibrium where A₂/A₁ = λ₂/(λ₂ − λ₁) > 1.

In radiometric geochronology, accumulating stable daughter nuclei D(t) = D₀ + N(t)[exp(λ*t*) − 1] in closed mineral lattices—such as U–Pb concordia dating in zircon (ZrSiO₄, which incorporates U⁴⁺ while rejecting initial Pb²⁺ at crystallization), ⁸⁷Rb–⁸⁷Sr isochron dating, ⁴⁰K–⁴⁰Ar / ⁴⁰Ar–³⁹Ar dating, and ¹⁴C radiocarbon dating (t_1/2 = 5,700 ± 30 yr)—dates archeological artifacts, mass extinctions, and the 4.567-billion-year age of the Solar System and Earth (Clair Patterson, 1956).

Can temperature, pressure, or chemical reactions speed up or slow down radioactive decay?

For alpha, beta-minus, beta-plus, and gamma decays, no terrestrial temperature, pressure, or chemical bond change alters the decay rate, because thermal and chemical energies (~0.01 to 10 eV) are a million times smaller than nuclear energy levels (~10⁵ to 10⁷ eV). The sole minor exception is orbital electron capture (EC) (and low-energy internal conversion) in light nuclei such as ⁷Be: because EC depends on the overlap of inner-shell atomic electrons with the nucleus (‖ψ_1s(0)‖²), extreme diamond-anvil-cell compression or full ionization can shift the ⁷Be EC decay rate by ~0.1% to 1.5%.

Why do alpha particles come out at discrete energies, whereas beta particles have a continuous energy spectrum?

Alpha decay (P → D + α) is a two-body breakup: conservation of total energy (Q) and linear momentum (p_α = −p_D) uniquely fixes the translational kinetic energy of the alpha particle to a single discrete value (K_α = Q(A − 4)/A) for each nuclear state transition. Beta decay (P → D + e⁻ + ν̄_e) is a three-body breakup: the fixed decay energy Q_β is shared continuously in all proportions between the emitted electron and the virtually undetectable antineutrino, producing a smooth bell-shaped beta energy distribution up to the maximum endpoint E_max = Q_β.

How does uranium-lead (U–Pb) zircon dating achieve internal cross-checking of rock ages?

Zircon crystals (ZrSiO₄) incorporate trace uranium (U⁴⁺ substituting for Zr⁴⁺) when crystallizing from magma but strongly exclude lead (Pb²⁺), setting initial radiogenic lead near zero. Crucially, natural uranium contains two independent radioactive isotopes that decay at different rates to two different stable lead isotopes: ²³⁸U → ²⁰⁶Pb (t_1/2 = 4.468 billion years) and ²³⁵U → ²⁰⁷Pb (t_1/2 = 703.8 million years). Plotting ²⁰⁶Pb/²³⁸U against ²⁰⁷Pb/²³⁵U defines the Wetherill concordia curve, providing a built-in internal check that detects and corrects for any post-crystallization lead loss.

1. Rutherford, E., & Soddy, F. (1902). "The cause and nature of radioactivity." Philosophical Magazine, Series 6, 4(21), 370–396. https://doi.org/10.1080/14786440209462856
2. Gamow, G. (1928). "Zur Quantentheorie des Atomkernes." Zeitschrift für Physik, 51(3–4), 204–212. https://doi.org/10.1007/BF01343196
3. Fermi, E. (1934). "Versuch einer Theorie der β-Strahlen. I." Zeitschrift für Physik, 88(3–4), 161–177. https://doi.org/10.1007/BF01351864
4. Patterson, C. (1956). "Age of meteorites and the earth." Geochimica et Cosmochimica Acta, 10(4), 230–237. https://doi.org/10.1016/0016-7037(56)90036-990036-9)