The mathematics of nuclear survival: exponential decay, activity in becquerels, parent-daughter Bateman chains, and the transient equilibrium of technetium generators.
IntuitionIntuition: a clock built from pure probability
Imagine tossing ten thousand coins simultaneously every minute and discarding every coin landing on heads. In the first minute approximately half disappear, leaving five thousand. In the second minute, half of what remains vanishes, leaving twenty-five hundred. No individual coin can predict when its turn will come, and a coin that has survived ten minutes is no more likely to land on heads than a newly minted one. Radioactive nuclei behave with this identical memoryless randomness: each unstable nucleus possesses a constant, unalterable probability of decaying in any given second, entirely independent of its past age.
Because macroscopic matter contains astronomical ensembles of nuclei—a single microgram of an isotope contains trillions of atoms—this microscopic randomness averages out into an extraordinarily smooth, predictable exponential law. The number of surviving nuclei halves across every fixed time window called the half-life. Whether tracking an imaging tracer decaying over hours or uranium charting terrestrial geology across billions of years, the exponential decay law furnishes science with an absolute, unyielding atomic chronometer.
Interactive decay spectrum. Inspect generic exponential halving (set 0), parent-daughter kinetics in the ⁹⁹Mo/⁹⁹ᵐTc medical generator peaking near 23 hours (set 1), and radiocarbon decay marks across millennia (set 2).
SchoolCore concepts: half-life, decay constant and activity
Definition: Radioactive decay law and activity
The rate of disintegration −dtdN of a radionuclide is directly proportional to the number of radioactive nuclei N(t) present: −dtdN=λN, where λ is the characteristic decay constant (s−1). Integrating yields N(t)=N0e−λt=N02−t/t1/2. The half-life t1/2=λln2≈λ0.69315 is the time required for half the initial nuclei to disintegrate. The activity A(t)=λN(t)=A0e−λt quantifies disintegrations per second, measured in becquerels (1 Bq=1 disintegration/s) or curies (1 Ci=3.7×1010 Bq).
Radiocarbon dating of organic archaeological artifacts
X53131X2532131I
8.02 days
β−,γ
Thyroid gland targeted radionuclide therapy and imaging
X4399mX243299mTc
6.01 hours
Isomeric transition (γ,140.5 keV)
Clinical diagnostic SPECT medical imaging
X2760X227260Co
5.27 years
β−,γ
Industrial food irradiation and external beam radiotherapy
X55137X2552137Cs
30.1 years
β−,γ
Fission product tracer and laboratory calibration source
X88226X2882226Ra
approx. 1600 years
α,γ
Historical Curie activity benchmark and brachytherapy needle
X92238X2922238U
4.47 billion years
α
Primordial geochronology and nuclear reactor fuel
Example: Radiocarbon dating and decay fraction calculations
An organic wooden artifact recovered from an archaeological site displays a 14C activity equal to 25% of modern wood (t1/2=5730 years). Determine how many half-lives have elapsed and compute the age of the artifact. In addition, calculate the number of half-lives required for the activity to drop below 0.1% (1/1000) of its initial level.
Solution
Fraction remaining: N/N0=0.25=(1/2)2, so exactly n=2 half-lives have passed.
Age of artifact: t=2×5730=11460 years (approx. 11.5 millennia).
For a fraction of 1/1000: (1/2)n=10−3⟹n=ln2ln1000≈0.693156.9078≈9.97 half-lives (approx. 10 half-lives, or ≈57100 years).
Example: Activity of a microgram sample of iodine-131
A carrier-free radiopharmaceutical vial contains 1.0μg of pure iodine-131 (X53131X2532131I, atomic mass M≈130.9 g/mol, t1/2=8.02 days). Calculate the decay constant λ in s−1, the total number of radioactive nuclei N, and the initial activity in becquerels (Bq) and gigabecquerels (GBq).
Solution
Half-life in seconds: t1/2=8.02×86400 s≈6.929×105 s.
Decay constant: λ=t1/2ln2=6.929×1050.69315≈1.000×10−6 s−1.
Number of nuclei: N=130.9 g/mol1.0×10−6 g×6.022×1023 mol−1≈4.60×1015 atoms.
Activity: A=λN=1.000×10−6 s−1×4.60×1015≈4.60×109 Bq=4.60 GBq (approx. 124 mCi).
UndergraduateUniversity tier: Bateman equations and decay equilibria
When a radioactive parent 1 disintegrates into a radioactive daughter 2, the daughter inventory obeys Harry Bateman’s coupled differential equations. Depending on the ratio of the half-lives, two primary equilibria govern radioactive families:
1. Secular equilibrium arises when the parent is enormously longer-lived than the daughter (λ1≪λ2, such as X226X22226Ra with t1/2=1600 y decaying into X222X22222Rn with t1/2=3.82 d). After several daughter half-lives, e−λ2t→0 and the daughter activity reaches exact parity with the parent: A2(t)→A1(t).
2. Transient equilibrium develops when the parent is longer-lived than the daughter but has a comparable lifespan (λ1<λ2), exemplified by the X99X2299Mo→X99mX2299mTc medical generator (t1/2=65.94 h vs 6.01 h). Daughter activity surpasses parent activity and attains a distinct peak at tmax=λ2−λ1ln(λ2/λ1)≈22.8 h (approx. 23 hours), after which both isotopes decay synchronously governed by the parent half-life: A2(t)→λ2−λ1λ2A1(t)≈1.10A1(t).
In radiopharmacology and environmental biophysics, administered radionuclides undergo simultaneous physical decay and biological clearance. The effective decay constant is the direct sum of physical and biological rates (λeff=λphys+λbio), yielding the reciprocal harmonic relation for effective half-life: Teff1=Tphys1+Tbio1. If a nuclide disintegrates via competing paths (branching decay with partial constants λa and λb), the total decay constant is λ=λa+λb, and the branching ratio for mode a is fa=λa/λ.
AdvancedAdvanced physical depth: Poisson statistics and quantum tunneling
Because individual radioactive disintegrations are independent Poisson events occurring in enormous populations, counting N disintegrations carries an irreducible statistical uncertainty with standard deviation σ=N and relative uncertainty σ/N=1/N. In low-background metrology or paleolithic radiocarbon dating, accumulating substantial total counts is the sole physical method to suppress counting noise. Moreover, nuclear decay constants span more than twenty orders of magnitude across the isotope chart. In 1928, George Gamow (along with Gurney and Condon) demonstrated that alpha half-lives originate from quantum mechanical tunneling: an alpha cluster strikes the nuclear Coulomb barrier roughly 1021 times per second, where the exponential barrier transmission probability P∝exp[−2πη] explains the empirical Geiger-Nuttall law, linking tiny increases in alpha particle kinetic energy to enormous multi-order drops in half-life.