Nuclear Chemistry Fundamentals and Radiochemistry Applications

Positron emission tomography (PET) relies on $\beta^+$-emitting nuclides — $^{18}$F ($t_{1/2}=109.8$ min), $^{11}$C (20.4 min), $^{68}$Ga (68 min) — attached to biologically active molecules such as fluorodeoxyglucose. Each positron annihilates with an electron within ~1 mm of emission, producing two collinear 511 keV photons detected in coincidence to localize metabolic activity, making PET central to oncology staging, cardiac viability assessment, and neurodegeneration research.

Single-photon emission computed tomography (SPECT) uses gamma-emitting nuclides directly, most commonly $^{99m}$Tc (140 keV gamma, $t_{1/2}=6.01$ h), produced on-demand from a $^{99}$Mo/$^{99m}$Tc generator ("technetium cow") exploiting the 66-hour half-life of the parent $^{99}$Mo to supply fresh daughter activity to hospitals for a week at a time — a direct clinical application of the same transient-equilibrium mathematics used in Example 3.

Targeted radiotherapy uses alpha emitters (e.g., $^{223}$Ra for bone metastases, exploiting radium's calcium-mimicking chemistry) and beta emitters (e.g., $^{131}$I for thyroid ablation, $^{177}$Lu-labeled peptides for neuroendocrine tumors) to deliver cytotoxic radiation with sub-millimeter precision, converting nuclear decay energy directly into localized biological damage while sparing surrounding tissue.

Radiometric dating spans timescales from $^{14}$C (5730 yr half-life, for archaeological and recent-geological samples up to ~50,000 years) to $^{238}$U/$^{206}$Pb, $^{40}$K/$^{40}$Ar, and $^{87}$Rb/$^{87}$Sr systems (half-lives of $10^9$–$10^{11}$ yr) for dating the oldest rocks and meteorites, using exactly the exponential decay and secular-equilibrium relations derived above.

Neutron activation analysis (NAA) irradiates a sample with a neutron flux, converting stable isotopes into radioactive ones via $(n,\gamma)$ capture, then identifies trace elements from the characteristic gamma-ray energies and half-lives of the induced activities — a technique sensitive to parts-per-billion concentrations, used in forensic science, archaeology (provenancing artifacts), and quality control of high-purity materials, all resting on the same decay-mode identification logic used throughout this article.

Nuclear power and weapons rest directly on the fission energetics derived above: $^{235}$U and $^{239}$Pu undergo neutron-induced fission releasing ~200 MeV per event, sustained as a chain reaction moderated (in reactors) to controlled criticality, while fusion of deuterium and tritium — releasing energy for the opposite reason on the binding-energy curve — remains the target of ongoing magnetic and inertial confinement fusion research (ITER, National Ignition Facility) as a prospective clean energy source.

Ask the tutor: How does the shell model's spin–orbit-split magic numbers quantitatively predict the half-lives of superheavy elements near the island of stability?

Every nucleus is a bound state of $Z$ protons and $N$ neutrons held together by the strong nuclear force, a short-range (~1–2 fm) attraction that acts almost equally between any pair of nucleons, competing against the long-range Coulomb repulsion that acts only between protons. Because the strong force saturates — each nucleon binds tightly only to its immediate neighbors — the binding energy per nucleon rises quickly with mass number $A$ for light nuclei (more neighbors to bond to) but then falls slowly for heavy nuclei as Coulomb repulsion, growing as $Z^2$, increasingly destabilizes the whole assembly. This produces the famous binding-energy curve that peaks near $^{56}$Fe/$^{62}$Ni and explains, in one picture, both why light nuclei fuse and heavy nuclei fission — both processes move nuclides toward the energetic floor of the curve. Superimposed on this smooth trend are quantum shell effects: nuclei with "magic numbers" of protons or neutrons (2, 8, 20, 28, 50, 82, 126) are anomalously stable, exactly analogous to the extra stability of noble-gas electron configurations, because nucleons fill discrete shell-model orbitals (analogous to atomic $1s, 1p, 1d2s,\ldots$ levels but generated by a nuclear mean-field potential plus strong spin–orbit coupling). Doubly magic nuclei such as $^{16}$O, $^{40}$Ca, and $^{208}$Pb are unusually inert; this shell structure is also why superheavy element synthesis targets a hypothesized "island of stability" around $Z=114$–126, $N\approx184$. The chart of nuclides — protons on one axis, neutrons on the other — shows a "belt of stability" that starts along the $N=Z$ diagonal for light elements and curves toward increasing neutron excess for heavy elements. This curvature is a direct, quantitative consequence of the Coulomb term growing as $Z^2$ (repulsion between every pair of protons) while the strong-force contribution grows only linearly with $A$; extra neutrons add binding without adding Coulomb repulsion, so heavy stable nuclei need proportionally more of them. Any nuclide off this belt decays by whichever mode moves it back toward the valley floor: neutron-rich nuclides undergo $\beta^-$ decay (a neutron converts to a proton), proton-rich (neutron-poor) nuclides undergo $\beta^+$ decay or electron capture (a proton converts to a neutron), and very heavy nuclides shed mass and charge via $\alpha$ decay or spontaneous fission.

Alpha decay is a quantum tunneling phenomenon: the alpha particle is preformed inside the heavy nucleus and tunnels through the Coulomb barrier that classically would confine it. Because the barrier height and width depend sensitively on the decay energy $Q_\alpha$, the Geiger–Nuttall relation links $\log t_{1/2}$ almost linearly to $Q_\alpha^{-1/2}$: small increases in decay energy produce enormous decreases in half-life, spanning from microseconds to billions of years across the actinide and superheavy series. Beta decay, by contrast, proceeds via the weak interaction, converting a neutron into a proton (or vice versa) and emitting an electron/antineutrino pair (or positron/neutrino pair); electron capture achieves the same proton-to-neutron conversion by absorbing an atomic (usually K-shell) electron, and is often accompanied by characteristic X-ray or Auger-electron emission as the resulting vacancy is filled. Gamma emission is not a separate decay mode in the sense of changing $Z$ or $N$; it is simply the electromagnetic de-excitation of a nucleus left in an excited state after alpha or beta decay, exactly analogous to fluorescence but with MeV-scale nuclear level spacings instead of eV-scale electronic ones.

The binding energy of a nucleus, via the semi-empirical (liquid-drop) mass formula, is

where: $a_V \approx 15.8$ MeV (volume term, bulk saturation of the strong force), $a_S \approx 18.3$ MeV (surface term, correcting for under-bonded surface nucleons), $a_C \approx 0.714$ MeV (Coulomb repulsion among $Z(Z-1)/2$ proton pairs, scaled by $A^{-1/3}\propto 1/R$), $a_A \approx 23.2$ MeV (asymmetry term penalizing $N\neq Z$, a Pauli-exclusion effect from filling separate proton/neutron Fermi seas), and $\delta(A,Z)$ the pairing term (+12/$\sqrt{A}$ MeV for even–even, 0 for odd-$A$, $-12/\sqrt{A}$ MeV for odd–odd nuclei), reflecting the extra binding from paired nucleon spins — a nuclear analogue of electron pairing energy in Crystal Field Theory and d-Orbital Splitting.

Reading: every term has a physical origin already familiar from bonding intuition — a "cohesive energy" term, a "surface tension" penalty, an electrostatic repulsion, an exchange/exclusion penalty for asymmetric filling, and a pairing stabilization — and their competition, not any single term, sets both the location of the stability valley and the energy released in fission or fusion. The radioactive decay law follows from treating decay as a memoryless, first-order stochastic process:

where: $N(t)$ is the number of undecayed parent nuclei at time $t$, $\lambda$ is the decay constant (probability of decay per nucleus per unit time, fixed by nuclear structure), and $t_{1/2}$ is the half-life.

Reading: because $\lambda$ depends only on nuclear structure (not on temperature, pressure, or chemical bonding — with rare, tiny exceptions for electron-capture nuclides sensitive to electron density at the nucleus), radioactive decay is an essentially perfect clock, and activity $A = \lambda N$ decreases with exactly the same exponential law as $N$ itself. The energy released in any nuclear transformation follows directly from the mass deficit:

where: $Q$ is the energy released (shared as kinetic energy of products and any emitted photons), and masses are atomic (or nuclear, consistently) rest masses in energy units (1 u = 931.494 MeV/$c^2$).

Reading: a decay or reaction is energetically allowed if and only if $Q > 0$; this single mass-bookkeeping criterion decides whether alpha decay, beta decay, fission, or fusion can occur spontaneously for a given nuclide, and its magnitude sets the kinetic energy spectrum of the emitted particles.

Example 1 (foundational): Half-life arithmetic for ¹⁸F in a PET tracer. Fluorine-18 has $t_{1/2} = 109.8$ min. A radiopharmacy ships a dose with activity $A_0 = 370$ MBq at synthesis. If transport and scanning together take 219.6 min (exactly 2 half-lives), find the activity delivered to the patient. Since decay is first-order,

With $t/t_{1/2} = 219.6/109.8 = 2.000$,

Now take a non-round, asymmetric check: suppose instead the elapsed time is 165 min (not a whole number of half-lives). Then

Limiting check: as $t \to 0$, $A \to A_0$ exactly (no decay yet), and as $t \to \infty$, $A \to 0$, both consistent with the exponential form; at $t = t_{1/2}$ exactly, $A = A_0/2 = 185$ MBq, confirming the formula reduces correctly to the definition of half-life.

Example 2 (applied): Decay-energy calculation for α decay of ²³⁸Pu, a radioisotope thermoelectric generator (RTG) fuel. $^{238}$Pu decays by alpha emission to $^{234}$U. Given atomic masses $m(^{238}\text{Pu}) = 238.049560$ u, $m(^{234}\text{U}) = 234.040952$ u, $m(^{4}\text{He}) = 4.002602$ u, find $Q_\alpha$ and the power output per gram assuming all decay energy converts to heat.

(matches the accepted value of 5.593 MeV for $^{238}$Pu, of which the alpha particle carries off ≈5.49 MeV and the recoiling $^{234}$U nucleus the small remainder, by momentum conservation). For power output, first find the decay constant from $t_{1/2} = 87.7$ yr $= 2.766\times10^{9}$ s:

Number of atoms per gram: $N = \dfrac{6.022\times10^{23}}{238.05} = 2.530\times10^{21}$ atoms/g. Activity per gram:

Power per gram (energy per decay converted to watts):

This is close to the accepted specific power of $^{238}$PuO$_2$ RTG fuel (~0.5 W/g accounting for the oxide's added mass and some energy loss), confirming the calculation.

Limiting check: if $Q_\alpha \to 0$ (mass of parent exactly equal to sum of products), $P \to 0$ regardless of $\lambda$ — no energy release means no power, as physically required, and this is the boundary between energetically allowed and forbidden alpha decay.

Example 3 (exam/competition-level): Secular equilibrium and activity ratio in the ²³⁸U decay chain used for U–Pb dating. In old rocks (age $\gg t_{1/2}$ of intermediate daughters), $^{238}$U ($t_{1/2} = 4.468\times10^9$ yr) sits in secular equilibrium with its much shorter-lived daughter $^{234}$Th ($t_{1/2}=24.10$ days) within the decay chain. Derive the equilibrium activity ratio and use it to find, for a rock of unknown age, the fraction of an initially pure $^{238}$U sample that has converted to stable $^{206}$Pb after $t = 1.80$ Gyr (asymmetric, non-round age).

Setting up the Bateman relation for two members: the daughter's population obeys

In secular equilibrium ($\lambda_2 \gg \lambda_1$, and after times long compared with $1/\lambda_2$ but short compared with $1/\lambda_1$), $dN_2/dt \to 0$ effectively on the timescale of interest, giving

so every daughter in a long secular-equilibrium chain has the same activity as the parent, even though $N_2 \ll N_1$ in absolute atom count — this is why a single old rock shows measurable activity from many short-lived daughters despite tiny standing populations of each.

Finding the age from the parent decay law: because the overall chain from $^{238}$U to stable $^{206}$Pb is rate-limited by the parent's enormous half-life, the number of Pb-206 atoms produced simply tracks U-238 depletion:

With $\lambda_1 = \ln 2 / (4.468\times10^9\ \text{yr}) = 1.5513\times10^{-10}\ \text{yr}^{-1}$ and $t = 1.80\times10^9$ yr,

So a 1.80 Gyr-old rock retains 75.6% of its original $^{238}$U and has accumulated $^{206}$Pb equal to 24.4% of the original uranium atom count — precisely the measurable $^{206}\text{Pb}/^{238}\text{U}$ ratio ($=0.2437/0.7563 = 0.3223$) that a geochronologist inverts, via this exact exponential relation, to read off rock age.

Limiting check: as $t\to 0$, the Pb fraction $\to 0$ (no time to decay), and as $t \to \infty$, the fraction $\to 1$ (complete conversion), both physically required; at $t = t_{1/2}$ exactly the fraction is precisely 0.5, recovering the definition of half-life as an internal consistency check on the Bateman-derived formula.

Marie Curie's laboratory notebooks are still measurably radioactive today, over a century later, a fact that is itself a small, poignant demonstration of the very physics she discovered: the half-life of $^{226}$Ra (1600 years) guarantees her belongings will remain hazardous for many human lifetimes. Her insistence on quantitative, reproducible measurement of activity — rather than qualitative description of "rays" — transformed radioactivity from a physical curiosity (discovered by Becquerel in uranium salts) into a genuine branch of chemistry, with its own separations, its own periodic patterns (the decay series, the transuranics later isolated by Seaborg and coworkers), and its own analytical applications that now underpin the entire framework connected to Actinide Chemistry and the 5f Series, since essentially every isotope of every actinide is understood, produced, and identified through exactly the decay-mode and half-life logic she pioneered.