Chemistry Labs

Grade 12

Half-life and the radioactive decay law

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 −dNdt-\frac{\mathrm{d}N}{\mathrm{d}t} of a radionuclide is directly proportional to the number of radioactive nuclei N(t)N(t) present: −dNdt=λN-\frac{\mathrm{d}N}{\mathrm{d}t} = \lambda N, where λ\lambda is the characteristic decay constant (s−1\text{s}^{-1}). Integrating yields N(t)=N0e−λt=N02−t/t1/2N(t) = N_0 \mathrm{e}^{-\lambda t} = N_0 2^{-t / t_{1/2}}. The half-life t1/2=ln⁡2λ≈0.69315λt_{1/2} = \frac{\ln 2}{\lambda} \approx \frac{0.69315}{\lambda} is the time required for half the initial nuclei to disintegrate. The activity A(t)=λN(t)=A0e−λtA(t) = \lambda N(t) = A_0 \mathrm{e}^{-\lambda t} quantifies disintegrations per second, measured in becquerels (1 Bq=1 disintegration/s1\text{ Bq} = 1\text{ disintegration/s}) or curies (1 Ci=3.7×1010 Bq1\text{ Ci} = 3.7 \times 10^{10}\text{ Bq}).

N(t)=N0e−λt=N0 2−t/t1/2,A(t)=λN(t),t1/2=ln⁡2λ,τ=1λ≈1.4427 t1/2N(t) = N_0 \mathrm{e}^{-\lambda t} = N_0 \, 2^{-t / t_{1/2}}, \qquad A(t) = \lambda N(t), \qquad t_{1/2} = \frac{\ln 2}{\lambda}, \qquad \tau = \frac{1}{\lambda} \approx 1.4427 \, t_{1/2}
Representative radionuclides and their half-lives
NuclideHalf-lifeDecay modePractical application
X614X26214C\ce{^{14}_6C}5730 yearsβ−\beta^-Radiocarbon dating of organic archaeological artifacts
X53131X2532131I\ce{^{131}_{53}I}8.02 daysβ−,γ\beta^-, \gammaThyroid gland targeted radionuclide therapy and imaging
X4399mX243299mTc\ce{^{99m}_{43}Tc}6.01 hoursIsomeric transition (γ,140.5 keV\gamma, 140.5\text{ keV})Clinical diagnostic SPECT medical imaging
X2760X227260Co\ce{^{60}_{27}Co}5.27 yearsβ−,γ\beta^-, \gammaIndustrial food irradiation and external beam radiotherapy
X55137X2552137Cs\ce{^{137}_{55}Cs}30.1 yearsβ−,γ\beta^-, \gammaFission product tracer and laboratory calibration source
X88226X2882226Ra\ce{^{226}_{88}Ra}approx. 1600 yearsα,γ\alpha, \gammaHistorical Curie activity benchmark and brachytherapy needle
X92238X2922238U\ce{^{238}_{92}U}4.47 billion yearsα\alphaPrimordial geochronology and nuclear reactor fuel

Example: Radiocarbon dating and decay fraction calculations

An organic wooden artifact recovered from an archaeological site displays a 14C^{14}\text{C} activity equal to 25%25\% of modern wood (t1/2=5730 yearst_{1/2} = 5730\text{ 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%0.1\% (1/10001/1000) of its initial level.

Solution

Fraction remaining: N/N0=0.25=(1/2)2N/N_0 = 0.25 = (1/2)^2, so exactly n=2n = 2 half-lives have passed. Age of artifact: t=2×5730=11460 yearst = 2 \times 5730 = 11460\text{ years} (approx. 11.511.5 millennia). For a fraction of 1/10001/1000: (1/2)n=10−3  ⟹  n=ln⁡1000ln⁡2≈6.90780.69315≈9.97 half-lives(1/2)^n = 10^{-3} \implies n = \frac{\ln 1000}{\ln 2} \approx \frac{6.9078}{0.69315} \approx 9.97\text{ half-lives} (approx. 10 half-lives, or ≈57100 years\approx 57100\text{ years}).

Example: Activity of a microgram sample of iodine-131

A carrier-free radiopharmaceutical vial contains 1.0 μg1.0\ \mu\text{g} of pure iodine-131 (X53131X2532131I\ce{^{131}_{53}I}, atomic mass M≈130.9 g/molM \approx 130.9\text{ g/mol}, t1/2=8.02 dayst_{1/2} = 8.02\text{ days}). Calculate the decay constant λ\lambda in s−1\text{s}^{-1}, the total number of radioactive nuclei NN, and the initial activity in becquerels (Bq) and gigabecquerels (GBq).

Solution

Half-life in seconds: t1/2=8.02×86400 s≈6.929×105 st_{1/2} = 8.02 \times 86400\text{ s} \approx 6.929 \times 10^5\text{ s}. Decay constant: λ=ln⁡2t1/2=0.693156.929×105≈1.000×10−6 s−1\lambda = \frac{\ln 2}{t_{1/2}} = \frac{0.69315}{6.929 \times 10^5} \approx 1.000 \times 10^{-6}\text{ s}^{-1}. Number of nuclei: N=1.0×10−6 g130.9 g/mol×6.022×1023 mol−1≈4.60×1015 atomsN = \frac{1.0 \times 10^{-6}\text{ g}}{130.9\text{ g/mol}} \times 6.022 \times 10^{23}\text{ mol}^{-1} \approx 4.60 \times 10^{15}\text{ atoms}. Activity: A=λN=1.000×10−6 s−1×4.60×1015≈4.60×109 Bq=4.60 GBqA = \lambda N = 1.000 \times 10^{-6}\text{ s}^{-1} \times 4.60 \times 10^{15} \approx 4.60 \times 10^9\text{ Bq} = 4.60\text{ GBq} (approx. 124 mCi124\text{ mCi}).

UndergraduateUniversity tier: Bateman equations and decay equilibria

dN2dt=λ1N1−λ2N2  ⟹  N2(t)=λ1λ2−λ1N1,0(e−λ1t−e−λ2t),tmax⁡=ln⁡(λ2/λ1)λ2−λ1\frac{\mathrm{d}N_2}{\mathrm{d}t} = \lambda_1 N_1 - \lambda_2 N_2 \implies N_2(t) = \frac{\lambda_1}{\lambda_2 - \lambda_1} N_{1,0} \left( \mathrm{e}^{-\lambda_1 t} - \mathrm{e}^{-\lambda_2 t} \right), \qquad t_{\max} = \frac{\ln(\lambda_2 / \lambda_1)}{\lambda_2 - \lambda_1}

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\lambda_1 \ll \lambda_2, such as X226X22226Ra\ce{^{226}Ra} with t1/2=1600 yt_{1/2} = 1600\text{ y} decaying into X222X22222Rn\ce{^{222}Rn} with t1/2=3.82 dt_{1/2} = 3.82\text{ d}). After several daughter half-lives, e−λ2t→0\mathrm{e}^{-\lambda_2 t} \to 0 and the daughter activity reaches exact parity with the parent: A2(t)→A1(t)A_2(t) \to A_1(t). 2. Transient equilibrium develops when the parent is longer-lived than the daughter but has a comparable lifespan (λ1<λ2\lambda_1 < \lambda_2), exemplified by the X99X2299Mo→X99mX2299mTc\ce{^{99}Mo} \to \ce{^{99m}Tc} medical generator (t1/2=65.94 ht_{1/2} = 65.94\text{ h} vs 6.01 h6.01\text{ h}). Daughter activity surpasses parent activity and attains a distinct peak at tmax⁡=ln⁡(λ2/λ1)λ2−λ1≈22.8 ht_{\max} = \frac{\ln(\lambda_2 / \lambda_1)}{\lambda_2 - \lambda_1} \approx 22.8\text{ h} (approx. 23 hours), after which both isotopes decay synchronously governed by the parent half-life: A2(t)→λ2λ2−λ1A1(t)≈1.10 A1(t)A_2(t) \to \frac{\lambda_2}{\lambda_2 - \lambda_1} A_1(t) \approx 1.10 \, A_1(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\lambda_{\text{eff}} = \lambda_{\text{phys}} + \lambda_{\text{bio}}), yielding the reciprocal harmonic relation for effective half-life: 1Teff=1Tphys+1Tbio\frac{1}{T_{\text{eff}}} = \frac{1}{T_{\text{phys}}} + \frac{1}{T_{\text{bio}}}. If a nuclide disintegrates via competing paths (branching decay with partial constants λa\lambda_a and λb\lambda_b), the total decay constant is λ=λa+λb\lambda = \lambda_a + \lambda_b, and the branching ratio for mode aa is fa=λa/λf_a = \lambda_a / \lambda.

AdvancedAdvanced physical depth: Poisson statistics and quantum tunneling

Because individual radioactive disintegrations are independent Poisson events occurring in enormous populations, counting NN disintegrations carries an irreducible statistical uncertainty with standard deviation σ=N\sigma = \sqrt{N} and relative uncertainty σ/N=1/N\sigma / N = 1/\sqrt{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 second10^{21}\text{ times per second}, where the exponential barrier transmission probability P∝exp⁡[−2πη]P \propto \exp[-2\pi\eta] explains the empirical Geiger-Nuttall law, linking tiny increases in alpha particle kinetic energy to enormous multi-order drops in half-life.

References

  • Introductory Nuclear Physics · K. S. Krane, 1987
  • Radiochemistry and Nuclear Chemistry, 4th Edition · G. Choppin, J.-O. Liljenzin, J. Rydberg, C. Ekberg, 2013
  • Atmospheric Helium Three and Radiocarbon from Cosmic Radiation · W. F. Libby, 1946