Case Study 2: Complete Analysis of the ${}^{238}\text{U}$ Decay Chain — From Uranium to Lead

Introduction

The ${}^{238}\text{U}$ decay chain — also called the uranium series or the $4n + 2$ series — is the longest of the four natural radioactive decay series. Starting from uranium-238 ($t_{1/2} = 4.468 \times 10^9\,\text{yr}$, comparable to the age of the Earth), it proceeds through 14 decay steps involving 8 alpha decays and 6 beta-minus decays before terminating at the stable nucleus ${}^{206}\text{Pb}$. Along the way, it passes through radium-226 (historically the first isolated radioactive element, by Marie Curie), radon-222 (the noble gas responsible for the dominant source of natural radiation exposure), and polonium-210 (a notoriously toxic alpha emitter).

This case study demonstrates a complete Option B analysis, integrating decay physics (Part III), nuclear structure (Part II), radiometric dating (Part VI), and environmental radiation (Part VI).

1. The Complete Chain

1.1 Chain Table

The complete ${}^{238}\text{U}$ decay series:

The net reaction is:

$${}^{238}\text{U} \to {}^{206}\text{Pb} + 8\alpha + 6e^- + 6\bar\nu_e + Q_{\text{total}}$$

The total energy released: $Q_{\text{total}} = 51.7\,\text{MeV}$, distributed among alpha particle kinetic energies, beta electron and neutrino energies, and gamma-ray energies from de-excitation.

1.2 Geiger-Nuttall Systematics

The alpha decays in the chain span a remarkable range of half-lives — from $4.5 \times 10^9\,\text{years}$ (${}^{238}\text{U}$) to $164\,\mu\text{s}$ (${}^{214}\text{Po}$) — a range of $10^{23}$ in half-life, corresponding to $Q_\alpha$ values from $4.27$ to $7.83\,\text{MeV}$. This extraordinary sensitivity of the half-life to the $Q$-value is the hallmark of quantum tunneling through the Coulomb barrier.

The Geiger-Nuttall relation for even-even alpha emitters in this chain:

$$\log_{10}(t_{1/2}/\text{s}) = a Z / \sqrt{Q_\alpha / \text{MeV}} + b$$

Plotting $\log_{10} t_{1/2}$ versus $Q_\alpha^{-1/2}$ for the eight alpha emitters yields a nearly linear relationship, confirming the WKB tunneling model of Chapter 13. The slope is determined by the Coulomb barrier height and width; the intercept incorporates the nuclear radius, the assault frequency, and the preformation probability of the alpha particle.

The general trend is clear: higher $Q_\alpha$ corresponds to shorter half-life. The scatter around a single line arises because the Geiger-Nuttall relation is approximate — it neglects details such as angular momentum barriers, nuclear deformation, and alpha preformation factors that differ between nuclei.

1.3 Beta Decay Transitions

The six beta-minus decays in the chain adjust the neutron-to-proton ratio as the chain moves through the chart of nuclides:

The large $\log ft$ values (all $> 6$) indicate that none of these are allowed transitions. The parity changes ($\pi_i \cdot \pi_f = -1$) and the angular momentum changes ($\Delta I = 0, 1$) are consistent with first-forbidden selection rules. The ${}^{210}\text{Pb}$ decay ($\log ft \approx 11.5$) is particularly slow — a first-forbidden unique transition ($\Delta I = 2$, parity change) with an exceptionally small $Q$-value of only $63.5\,\text{keV}$.

2. Bateman Equation Solutions

2.1 The Differential Equations

For a chain of $n$ nuclides with no branching:

$$\frac{dN_1}{dt} = -\lambda_1 N_1$$

$$\frac{dN_i}{dt} = \lambda_{i-1} N_{i-1} - \lambda_i N_i \quad (i = 2, 3, \ldots, n)$$

$$\frac{dN_n}{dt} = \lambda_{n-1} N_{n-1} \quad \text{(stable end product)}$$

with initial conditions $N_1(0) = N_0$, $N_i(0) = 0$ for $i > 1$ (starting with pure ${}^{238}\text{U}$).

2.2 Analytical Solution for the First Three Members

The Bateman solution for the $i$-th member (assuming distinct decay constants) is:

$$N_i(t) = N_0 \prod_{j=1}^{i-1} \lambda_j \sum_{k=1}^{i} \frac{e^{-\lambda_k t}}{\prod_{m=1, m \neq k}^{i} (\lambda_m - \lambda_k)}$$

For the first three members (${}^{238}\text{U}$, ${}^{234}\text{Th}$, ${}^{234}\text{Pa}$):

${}^{238}\text{U}$:

$$N_1(t) = N_0 e^{-\lambda_1 t} \approx N_0 \quad \text{(for } t \ll 4.5 \times 10^9\,\text{yr}\text{)}$$

${}^{234}\text{Th}$:

$$N_2(t) = N_0 \frac{\lambda_1}{\lambda_2 - \lambda_1}\left(e^{-\lambda_1 t} - e^{-\lambda_2 t}\right)$$

Since $\lambda_2 \gg \lambda_1$ ($t_{1/2,2} = 24.1\,\text{d} \ll t_{1/2,1} = 4.5 \times 10^9\,\text{yr}$), after a few half-lives of ${}^{234}\text{Th}$:

$$N_2 \approx N_0 \frac{\lambda_1}{\lambda_2} = N_0 \frac{t_{1/2,2}}{t_{1/2,1}}$$

This is secular equilibrium: $A_2 = \lambda_2 N_2 = \lambda_1 N_1 = A_1$.

${}^{234}\text{Pa}$: Similarly reaches secular equilibrium with $A_3 = A_2 = A_1$.

2.3 Numerical Solution

The full 14-step chain with branching is solved numerically using the decay_chains.py toolkit module. The ODE system is stiff (decay constants span 23 orders of magnitude), requiring an implicit solver such as scipy.integrate.solve_ivp with the BDF or Radau method.

Key features of the solution:

1. Short-lived daughters build up quickly. ${}^{234}\text{Th}$ reaches secular equilibrium within $\sim 4$ months (about 5 half-lives of 24.1 days).

2. ${}^{234}\text{U}$ builds up slowly. With $t_{1/2} = 2.455 \times 10^5\,\text{yr}$, it takes $\sim 10^6\,\text{yr}$ to approach secular equilibrium.

3. ${}^{226}\text{Ra}$ and below are all in secular equilibrium with ${}^{238}\text{U}$ in natural uranium ores that have been undisturbed for $> 10^6\,\text{yr}$.

4. The total activity of a sample in secular equilibrium is 14 times the activity of the parent (since there are 14 radioactive members, each with $A_i = A_1$).

5. ${}^{214}\text{Po}$ ($t_{1/2} = 164\,\mu\text{s}$) never accumulates to measurable amounts — its activity equals $A_1$, but its number of atoms is negligible.

2.4 Activity Plot Description

A plot of activity versus time for the first $10^7$ years (starting with pure ${}^{238}\text{U}$) shows:

• ${}^{238}\text{U}$ activity: essentially constant (barely decayed)
• ${}^{234}\text{Th}$ and ${}^{234}\text{Pa}$: jump to $A_1$ within days/hours (invisible on this timescale)
• ${}^{234}\text{U}$: exponential growth toward $A_1$ with time constant $\sim t_{1/2}/\ln 2 \approx 3.5 \times 10^5\,\text{yr}$
• ${}^{230}\text{Th}$: delayed growth, approaching $A_1$ on a timescale of $\sim 5 \times 10^5\,\text{yr}$
• ${}^{226}\text{Ra}$ and below: even more delayed, reaching equilibrium after $\sim 10^4$--$10^5\,\text{yr}$

3. Equilibrium Analysis

3.1 Secular Equilibrium

Since $t_{1/2}({}^{238}\text{U}) = 4.468 \times 10^9\,\text{yr}$ is vastly longer than all daughter half-lives, every parent-daughter pair in the chain satisfies the secular equilibrium condition ($t_{1/2,\text{parent}} \gg t_{1/2,\text{daughter}}$) with respect to ${}^{238}\text{U}$.

At secular equilibrium:

$$A_1 = A_2 = A_3 = \cdots = A_{14}$$

where $A_i = \lambda_i N_i$ is the activity (decays per second) of the $i$-th member. This means:

$$\lambda_1 N_1 = \lambda_2 N_2 = \cdots = \lambda_{14} N_{14}$$

or equivalently:

$$N_i = N_1 \frac{\lambda_1}{\lambda_i} = N_1 \frac{t_{1/2,i}}{t_{1/2,1}}$$

Numerical example: For 1 gram of natural uranium (99.27% ${}^{238}\text{U}$):

$$N_1 = \frac{0.9927 \times 6.022 \times 10^{23}}{238} = 2.51 \times 10^{21}\,\text{atoms}$$

$$A_1 = \lambda_1 N_1 = \frac{0.693}{4.468 \times 10^9 \times 3.156 \times 10^7} \times 2.51 \times 10^{21} = 1.23 \times 10^4\,\text{Bq}$$

This is $12,300\,\text{Bq} = 0.33\,\mu\text{Ci}$ — the specific activity of natural uranium.

At secular equilibrium, each daughter also contributes $12,300\,\text{Bq}$, so the total activity of the chain is $14 \times 12,300 = 172,000\,\text{Bq}$ per gram of uranium.

3.2 Time to Equilibrium

The time for the slowest daughter (${}^{234}\text{U}$, $t_{1/2} = 2.455 \times 10^5\,\text{yr}$) to reach 99% of its equilibrium activity is:

$$t_{99\%} = \frac{\ln 100}{\lambda_2} = \frac{4.605 \times t_{1/2}}{\ln 2} = 4.605 \times 2.455 \times 10^5 / 0.693 \approx 1.63 \times 10^6\,\text{yr}$$

For the entire chain to be within 1% of secular equilibrium requires about $1.6\,\text{million years}$ — a blink compared to the $4.5\,\text{billion year}$ half-life of ${}^{238}\text{U}$.

4. Radiometric Dating: The U-Pb System

4.1 The Age Equation

The fundamental equation for U-Pb dating from the ${}^{238}\text{U}$ chain is:

$${}^{206}\text{Pb}^* = {}^{238}\text{U}\left(e^{\lambda_{238} t} - 1\right)$$

where ${}^{206}\text{Pb}^*$ is the radiogenic lead (produced by uranium decay) and $t$ is the age. Including initial (non-radiogenic) lead:

$$\frac{{}^{206}\text{Pb}}{{}^{204}\text{Pb}} = \left(\frac{{}^{206}\text{Pb}}{{}^{204}\text{Pb}}\right)_0 + \frac{{}^{238}\text{U}}{{}^{204}\text{Pb}}\left(e^{\lambda_{238} t} - 1\right)$$

where ${}^{204}\text{Pb}$ is the non-radiogenic reference isotope (not produced in any decay chain).

4.2 The Concordia Diagram

A more powerful technique uses both uranium decay chains simultaneously. Defining:

$$x = \frac{{}^{207}\text{Pb}^*}{{}^{235}\text{U}} = e^{\lambda_{235} t} - 1$$

$$y = \frac{{}^{206}\text{Pb}^*}{{}^{238}\text{U}} = e^{\lambda_{238} t} - 1$$

The locus of $(x, y)$ for different ages $t$ defines the concordia curve. A sample that has remained a closed system (no lead loss or uranium gain) since its formation at time $t$ plots on the concordia at the point corresponding to $t$.

Lead loss: If a sample experienced lead loss at some time $t_1 < t_0$ (where $t_0$ is the formation age), the measured $(x, y)$ point falls below the concordia, on a line connecting the concordia points at $t_0$ and $t_1$. This discordia line intersects the concordia at two points, giving both the formation age and the lead-loss event age.

4.3 Application: Age of the Earth

The age of the Earth was first determined by Clair Patterson in 1956 using U-Pb data from meteorites. The key insight: meteorites formed from the same solar nebula as Earth and have remained closed systems since formation. Patterson measured lead isotope ratios in iron meteorites (no uranium, so only primordial lead) and stony meteorites (with uranium, so primordial + radiogenic lead). The intercept of the isochron gives:

$$t_{\text{Earth}} = 4.55 \pm 0.07\,\text{Gyr}$$

This result, derived entirely from the nuclear physics of the ${}^{238}\text{U}$ and ${}^{235}\text{U}$ decay chains, is one of the most important measurements in the history of science. Modern U-Pb measurements on the oldest zircon crystals give $t = 4.404 \pm 0.008\,\text{Gyr}$ for the oldest terrestrial material.

4.4 Sources of Uncertainty

1. Decay constant uncertainty: $\lambda_{238} = (1.55125 \pm 0.00017) \times 10^{-10}\,\text{yr}^{-1}$ — known to $0.011\%$, the most precisely determined of all radiometric constants.
2. Common lead correction: Subtracting the initial lead requires knowledge of the primordial lead isotopic composition (from iron meteorites or model-dependent estimates).
3. Open system behavior: Lead loss (e.g., from metamorphic events) or uranium mobility can shift samples off the concordia.
4. Intermediate daughter disequilibrium: If ${}^{234}\text{U}$ is not in secular equilibrium with ${}^{238}\text{U}$ (common in low-temperature aqueous environments), the simple age equation requires correction. This is important for samples younger than $\sim 10^6$ years.

5. Environmental Radiation: The Radon Problem

5.1 Radon-222 in the Decay Chain

${}^{222}\text{Rn}$ (Step 7 in the chain) is a noble gas with $t_{1/2} = 3.823\,\text{d}$. Unlike all other members of the chain (which are solids), radon can escape from rock and soil into the atmosphere and accumulate in enclosed spaces — particularly basements and ground floors of buildings.

Once inhaled, radon itself is relatively harmless (it is exhaled before decaying). The danger comes from its solid daughters: ${}^{218}\text{Po}$, ${}^{214}\text{Pb}$, ${}^{214}\text{Bi}$, and ${}^{214}\text{Po}$. These "radon daughters" or "radon progeny" are electrically charged when formed (the alpha decay of radon produces a charged polonium ion) and attach to aerosol particles in the lung, where they continue to decay. The critical dose comes from the alpha emissions of ${}^{218}\text{Po}$ ($6.115\,\text{MeV}$) and ${}^{214}\text{Po}$ ($7.833\,\text{MeV}$), which deposit their energy in a thin layer of bronchial epithelial tissue.

5.2 Dose Calculation

The average indoor radon concentration in the United States is approximately $48\,\text{Bq/m}^3$ ($1.3\,\text{pCi/L}$). The EPA action level is $148\,\text{Bq/m}^3$ ($4\,\text{pCi/L}$).

The effective dose from radon exposure can be estimated:

1. Equilibrium factor $F$: In typical indoor air, radon daughters are not in full equilibrium with radon (some plate out on surfaces). The equilibrium factor $F \approx 0.4$ represents the ratio of actual daughter concentration to equilibrium concentration.

2. Dose conversion factor: The ICRP recommends a dose coefficient of $9\,\text{nSv/(Bq}\cdot\text{h/m}^3\text{)}$ for radon progeny, including the equilibrium factor.

3. Annual dose: For 7000 hours/year indoor occupancy at $48\,\text{Bq/m}^3$:

$$H = 48 \times 0.4 \times 9 \times 10^{-9} \times 7000 \approx 1.2\,\text{mSv/yr}$$

This is the dominant component of natural background radiation exposure for most people — about half the total natural background of $\sim 2.4\,\text{mSv/yr}$.

5.3 The Nuclear Physics Behind Radon Risk

The radon risk is a direct consequence of nuclear physics:

1. Why radon is a gas: $Z = 86$ gives radon the electron configuration of a noble gas. This is a chemical property, but it determines the transport pathway.

2. Why ${}^{222}\text{Rn}$ is the dangerous isotope: Its 3.82-day half-life is long enough to escape from soil and accumulate indoors, but short enough that it (and its daughters) produce significant activity. The other radon isotopes — ${}^{220}\text{Rn}$ ($t_{1/2} = 55.6\,\text{s}$, thoron series) and ${}^{219}\text{Rn}$ ($t_{1/2} = 3.96\,\text{s}$, actinium series) — are too short-lived to accumulate significantly.

3. Why the daughters are dangerous: The alpha energies (6.1 and 7.8 MeV) deposit $\sim 100\,\text{keV/}\mu\text{m}$ in tissue — high linear energy transfer (LET) radiation that causes dense ionization tracks, leading to double-strand DNA breaks and elevated cancer risk.

6. Connection to Reactor Physics

6.1 The Neptunium Production Chain

When ${}^{238}\text{U}$ captures a neutron in a reactor, it initiates a chain that produces fissile plutonium:

$${}^{238}\text{U}(n,\gamma){}^{239}\text{U} \xrightarrow[\beta^-]{23.5\,\text{min}} {}^{239}\text{Np} \xrightarrow[\beta^-]{2.357\,\text{d}} {}^{239}\text{Pu}$$

${}^{239}\text{Pu}$ is fissile and contributes significantly to power production in thermal reactors (about 40% of energy in a typical PWR comes from plutonium fission by end of fuel cycle). It is also the primary material for plutonium-based nuclear weapons — the connection between peaceful and military nuclear technology that drives nonproliferation policy.

6.2 ${}^{238}\text{U}$ as Fertile Material

${}^{238}\text{U}$ is not fissile (it does not sustain a chain reaction with thermal neutrons), but it is fertile: it can be converted to fissile ${}^{239}\text{Pu}$ by neutron capture. In a breeder reactor, this conversion can be made efficient enough that the reactor produces more fissile material than it consumes — a concept that, if deployed at scale, would extend uranium fuel resources by a factor of $\sim 60$.

The fast fission cross section of ${}^{238}\text{U}$ ($\sigma_f \approx 0.3\,\text{b}$ at $E > 1\,\text{MeV}$) also contributes the fast fission factor $\epsilon \approx 1.03$ in the four-factor formula.

7. Synthesis

The ${}^{238}\text{U}$ decay chain is a microcosm of nuclear physics:

Nuclear structure (Parts I-II): The alternation of alpha and beta decays traces a path through the chart of nuclides governed by the interplay between nuclear binding (which favors certain $N/Z$ ratios) and the energetics of alpha and beta decay (which depend on mass differences between parent and daughter). The $Q_\alpha$ values are set by the binding energy landscape; the Geiger-Nuttall systematics encode the quantum tunneling physics of the Coulomb barrier.

Decay physics (Part III): Every major decay mode is represented: alpha emission (8 steps), beta-minus (6 steps), and the competition between them at branch points. The half-life range of $10^{23}$ (from $4.5 \times 10^9\,\text{yr}$ to $164\,\mu\text{s}$) is explained by the exponential sensitivity of tunneling to the barrier height and width.

Applications (Part VI): The chain connects to radiometric dating (the age of the Earth), environmental radiation (radon exposure), reactor physics (fertile-to-fissile conversion), and nuclear forensics (isotopic signatures of uranium enrichment and reactor operations).

Astrophysics (Part V): ${}^{238}\text{U}$ itself was produced by the r-process, likely in a neutron star merger. The fact that any ${}^{238}\text{U}$ remains on Earth (its half-life is comparable to the age of the solar system) constrains the timing and yield of r-process nucleosynthesis in the early Milky Way.

The thread connecting all of this is the $Q$-value — the mass difference between parent and daughter that determines whether a decay is allowed, how fast it proceeds, and how much energy it releases. The $Q$-values are ultimately set by the nuclear binding energies, which are determined by the nuclear force. From the nuclear force to the age of the Earth, from femtometers to billions of years, the ${}^{238}\text{U}$ decay chain embodies the extraordinary reach of nuclear physics.