Case Study 2 — $^{208}$Pb: The Perfect Doubly Magic Nucleus

"If you want to understand nuclear structure, study lead-208. If you still don't understand, study it harder." — Attributed to a senior experimentalist at Oak Ridge (possibly apocryphal, but widely quoted)

The Status of $^{208}$Pb in Nuclear Physics

$^{208}$Pb occupies a singular position in nuclear physics. It is the heaviest stable nucleus that is doubly magic — with $Z = 82$ protons filling through the sixth proton shell closure and $N = 126$ neutrons filling through the seventh neutron shell closure. It is the end point of three natural radioactive decay chains (the thorium, uranium, and actinium series), which means that a substantial fraction of the lead in the Earth's crust is $^{208}$Pb produced by radioactive decay over geological time. It is the single most studied nucleus in nuclear structure physics, with thousands of published measurements of its properties. And it provides the most stringent, most precise, and most successful test of the nuclear shell model.

This case study examines the properties of $^{208}$Pb through the lens of the shell model, demonstrating how a theoretical framework built from a mean-field potential and spin-orbit coupling explains — quantitatively — the observed behavior of a real nucleus.

Shell Structure

The proton configuration of $^{208}$Pb fills all orbits through the $Z = 82$ shell closure:

$$\pi: \underbrace{(1s_{1/2})^2}_{2} \; \underbrace{(1p_{3/2})^4 (1p_{1/2})^2}_{8} \; \underbrace{(1d_{5/2})^6 (2s_{1/2})^2 (1d_{3/2})^4}_{20}$$

$$\underbrace{(1f_{7/2})^8}_{28} \; \underbrace{(2p_{3/2})^4 (1f_{5/2})^6 (2p_{1/2})^2 (1g_{9/2})^{10}}_{50}$$

$$\underbrace{(1g_{7/2})^8 (2d_{5/2})^6 (2d_{3/2})^4 (3s_{1/2})^2 (1h_{11/2})^{12}}_{82}$$

The neutron configuration fills through $N = 126$:

$$\nu: [\text{same as protons through } N = 82]$$

$$\underbrace{(2f_{7/2})^8 (1h_{9/2})^{10} (3p_{3/2})^4 (2f_{5/2})^6 (3p_{1/2})^2 (1i_{13/2})^{14}}_{126}$$

Every single-particle orbit up to the respective Fermi surface is completely filled. There are no unpaired nucleons, no partially occupied orbits, no valence particles. The nucleus is a perfect closed-shell system.

Ground-State Properties

Spin and Parity: $J^{\pi} = 0^+$

Since all orbits are filled, all nucleon angular momenta pair to zero. The ground-state spin is $J = 0$, and the parity is $\pi = +1$ (the product of parities of all occupied orbits in pairs gives $+1$). This is confirmed experimentally.

Magnetic Dipole Moment: $\mu = 0$

A $J = 0$ state has no orientation in space, so all electromagnetic multipole moments vanish. The measured magnetic moment is zero. This is a necessary consequence of the $0^+$ ground state, not a separate prediction.

Electric Quadrupole Moment: $Q = 0$

A $J = 0$ state also has zero quadrupole moment. This is again required by angular momentum selection rules. But the spectroscopic quadrupole moment of the $2^+_1$ excited state is small compared to mid-shell nuclei, confirming that $^{208}$Pb is very close to spherical.

Binding Energy

The experimental binding energy of $^{208}$Pb is:

$$B_{\text{exp}} = 1636.43 \text{ MeV}$$

The semi-empirical mass formula (with standard Bethe-Weizsacker coefficients from Chapter 4) predicts approximately 1629 MeV. The positive residual of $\delta B \approx +7$ MeV represents the extra binding from the simultaneous closure of both proton and neutron shells. This is among the largest shell corrections in the chart of nuclides.

The binding energy per nucleon is $B/A = 7.868$ MeV, which sits near the peak of the $B/A$ curve (the maximum is at $^{56}$Fe with $B/A = 8.790$ MeV). $^{208}$Pb is less bound per nucleon than iron because the increasing Coulomb repulsion in heavy nuclei reduces $B/A$. But among nuclei with $A > 150$, $^{208}$Pb is remarkably tightly bound — the shell effect partially compensates the Coulomb penalty.

Separation Energies: Direct Evidence for Shell Gaps

The one-neutron separation energy provides the most direct measurement of the shell gap:

$$S_n(^{208}\text{Pb}) = B(^{208}\text{Pb}) - B(^{207}\text{Pb}) = 7.37 \text{ MeV}$$

$$S_n(^{209}\text{Pb}) = B(^{209}\text{Pb}) - B(^{208}\text{Pb}) = 3.94 \text{ MeV}$$

The drop of 3.43 MeV in $S_n$ upon crossing $N = 126$ directly measures the gap between the last occupied neutron orbit ($3p_{1/2}$ at $N = 126$) and the first unoccupied orbit ($2g_{9/2}$ at $N = 127$). Similarly, for protons:

$$S_p(^{208}\text{Pb}) = B(^{208}\text{Pb}) - B(^{207}\text{Tl}) = 8.01 \text{ MeV}$$

$$S_p(^{209}\text{Bi}) = B(^{209}\text{Bi}) - B(^{208}\text{Pb}) = 3.80 \text{ MeV}$$

The proton shell gap (drop of 4.21 MeV) is even larger. Both gaps are among the largest in the nuclear chart, reflecting the exceptional robustness of the $Z = 82$ and $N = 126$ closures.

Single-Particle States: Particle and Hole Spectra

The most detailed test of the shell model comes from the single-particle states in nuclei adjacent to $^{208}$Pb. If $^{208}$Pb truly has a closed-shell structure, then:

  • $^{209}$Pb ($N = 127$) should show single-neutron states above the $N = 126$ gap
  • $^{207}$Pb ($N = 125$) should show single-neutron-hole states in the $N = 82$-$126$ shell
  • $^{209}$Bi ($Z = 83$) should show single-proton states above the $Z = 82$ gap
  • $^{207}$Tl ($Z = 81$) should show single-proton-hole states in the $Z = 50$-$82$ shell

Neutron Particle States: $^{209}$Pb

The low-lying spectrum of $^{209}$Pb:

$E_x$ (MeV) $J^{\pi}$ Shell model orbit
0 $9/2^+$ $2g_{9/2}$
0.779 $11/2^+$ $1i_{11/2}$
1.423 $15/2^-$ $1j_{15/2}$
1.567 $5/2^+$ $3d_{5/2}$
2.032 $1/2^+$ $4s_{1/2}$
2.150 $3/2^+$ $3d_{3/2}$
2.491 $7/2^+$ $2g_{7/2}$

These states map directly onto the single-neutron orbits above $N = 126$. The $2g_{9/2}$ ground state confirms the shell model filling order. The energies give the single-particle energies in the next major shell.

Neutron Hole States: $^{207}$Pb

The low-lying spectrum of $^{207}$Pb:

$E_x$ (MeV) $J^{\pi}$ Shell model orbit (hole)
0 $1/2^-$ $3p_{1/2}$
0.570 $5/2^-$ $2f_{5/2}$
0.898 $3/2^-$ $3p_{3/2}$
1.633 $13/2^+$ $1i_{13/2}$
2.340 $7/2^-$ $2f_{7/2}$

These are the orbits of the $N = 82$-$126$ shell, observed as hole states. The $3p_{1/2}$ ground state is the last orbit to fill before $N = 126$, consistent with the shell model level ordering.

The energy spacings between these hole states give the single-particle energy differences within the $N = 82$-$126$ shell. For example, the $3p_{1/2}$-$2f_{5/2}$ gap is 0.570 MeV, and the $3p_{1/2}$-$1i_{13/2}$ gap is 1.633 MeV.

Proton States: $^{209}$Bi and $^{207}$Tl

Similar spectra exist for the proton particle ($^{209}$Bi) and proton hole ($^{207}$Tl) nuclei. The $^{209}$Bi ground state is $9/2^-$ ($1h_{9/2}$), confirming the first proton orbit above $Z = 82$. The $^{207}$Tl ground state is $1/2^+$ ($3s_{1/2}$), the last proton orbit below $Z = 82$.

The First Excited State: $E(3^-_1) = 2.615$ MeV

The first excited state of $^{208}$Pb is at 2.615 MeV with $J^{\pi} = 3^-$. This is a remarkably high first excitation energy — by comparison, the first excited state of $^{238}$U (a well-deformed mid-shell nucleus) is at 0.045 MeV, lower by a factor of 58.

The $3^-$ state is predominantly a collective octupole vibration — a coherent superposition of many one-particle-one-hole excitations across the shell gap. The fact that even this collective mode sits at 2.6 MeV testifies to the large shell gaps in $^{208}$Pb. In nuclei where the shell gaps are small, the first $3^-$ state typically appears below 1 MeV.

The first $2^+$ state — the more common first excitation in non-magic even-even nuclei — occurs at 4.085 MeV in $^{208}$Pb. This is one of the highest $E(2^+_1)$ values known, and it demonstrates the extreme resistance of $^{208}$Pb to quadrupole deformation.

Charge Radius and Nuclear Shape

Electron scattering experiments at Saclay, Mainz, and Jefferson Lab have measured the charge distribution of $^{208}$Pb with extraordinary precision. The charge radius is:

$$R_{\text{ch}} = 5.5012 \pm 0.0013 \text{ fm}$$

The charge distribution is well described by a Fermi function with half-density radius $c = 6.624$ fm and diffuseness $a = 0.549$ fm. The higher-order moments of the distribution confirm that $^{208}$Pb is spherical to an excellent approximation.

The Neutron Skin: Connecting Structure to Astrophysics

Because $^{208}$Pb has 44 more neutrons than protons, the neutron distribution extends slightly beyond the proton distribution. The difference $R_n - R_p$ — the neutron skin thickness — is a quantity of profound importance because it directly constrains the symmetry energy of nuclear matter and, through the equation of state, the structure of neutron stars.

The PREX (Lead Radius Experiment) at Jefferson Lab measured the neutron skin using parity-violating electron scattering — a technique that exploits the fact that the weak neutral current couples preferentially to neutrons. The PREX-2 result (2021) is:

$$R_n - R_p = 0.283 \pm 0.071 \text{ fm}$$

This measurement connects the shell model of $^{208}$Pb directly to the mass-radius relationship of neutron stars (Chapter 25). A larger neutron skin implies a stiffer symmetry energy, which in turn predicts larger neutron star radii. The PREX result is consistent with neutron star radius measurements from NICER and gravitational wave observations from LIGO/Virgo.

$^{208}$Pb as a Calibration Standard

Because its properties are so well understood theoretically, $^{208}$Pb serves as a calibration standard for:

  1. Nuclear energy density functionals. Every Skyrme, Gogny, or relativistic mean-field parametrization is fitted to reproduce the binding energy, charge radius, and single-particle spectrum of $^{208}$Pb.

  2. Ab initio nuclear theory. The ability to compute $^{208}$Pb from the bare NN interaction is a benchmark for modern many-body methods. As of the mid-2020s, coupled-cluster and in-medium similarity renormalization group calculations can reproduce $^{208}$Pb properties to within a few percent — a triumph of nuclear theory.

  3. Neutrino physics. $^{208}$Pb is a detector material for supernova neutrinos (the HALO detector at SNOLAB) because its well-understood nuclear structure allows precise predictions of neutrino cross sections.

  4. Radiation shielding. Lead is the standard shielding material for gamma radiation, and $^{208}$Pb's large neutron capture cross section minimum (a consequence of the $N = 126$ closure) influences the isotopic choice for shielding applications.

Nuclear Reactions and $^{208}$Pb

The magic character of $^{208}$Pb profoundly influences nuclear reaction physics. In neutron capture reactions, $^{208}$Pb has one of the smallest thermal capture cross sections of any nucleus: $\sigma_{n,\gamma} \approx 0.49$ mb. For comparison, neighboring non-magic nuclei like $^{197}$Au have cross sections nearly $10^5$ times larger. This makes lead an effective neutron moderator and reflector material in some reactor designs, particularly lead-bismuth eutectic cooled reactors.

In proton-nucleus and heavy-ion reactions, $^{208}$Pb is used as a target for precision measurements precisely because its well-understood structure provides a clean theoretical baseline. Coulomb excitation experiments with $^{208}$Pb targets have been used to measure electromagnetic transition rates in projectile nuclei, and elastic scattering from $^{208}$Pb is a standard probe for optical model parametrizations.

$^{208}$Pb in the Cosmos

The shell closure at $N = 126$ has dramatic consequences for nucleosynthesis. In the slow neutron capture process (s-process) that operates in asymptotic giant branch (AGB) stars, the nuclear reaction flow encounters a bottleneck at $^{208}$Pb because its small neutron capture cross section slows the capture chain. Material accumulates at $A = 208$, producing the well-known s-process abundance peak near lead and bismuth. The solar system abundance of $^{208}$Pb is enhanced by a factor of approximately 3 relative to neighboring isotopes — a direct imprint of the $N = 126$ magic number on the composition of the cosmos.

In the rapid neutron capture process (r-process), which produces about half of all elements heavier than iron, the magic number $N = 126$ creates a "waiting point" where the r-process path stalls because beta-decay lifetimes are long and neutron capture rates are suppressed. The r-process abundance peak near $A \approx 195$ (shifted down from $A = 208$ because the r-process path runs through neutron-rich nuclei with $N = 126$ but $Z < 82$) is another cosmic fingerprint of the nuclear shell model.

Summary: Why $^{208}$Pb Matters

$^{208}$Pb is not merely a textbook example — it is the single most important nucleus for validating nuclear structure theory. Its doubly magic character makes it amenable to the cleanest possible theoretical treatment; its experimental accessibility (it is stable, abundant, and heavy enough for precision measurements) makes it the most thoroughly measured nucleus in existence; and its connections to astrophysics, neutrino physics, and applied nuclear science ensure its continued relevance at the forefront of research.

Every property of $^{208}$Pb — its $0^+$ ground state, its high first excited state, its spherical shape, its separation energy discontinuities, the single-particle spectra of its neighbors, its small neutron capture cross section, and its prominence in cosmic abundance patterns — confirms the shell model. If the shell model is the foundation of nuclear structure, then $^{208}$Pb is the cornerstone of that foundation.

Discussion Questions

  1. The SEMF residual for $^{208}$Pb is about +7 MeV. Given that the total binding energy is 1636 MeV, this represents only 0.4% of the total. Yet it has enormous physical consequences. Explain why a "small" correction to the binding energy can determine whether a nucleus is stable or radioactive.

  2. The neutron skin of $^{208}$Pb constrains the equation of state of neutron-rich matter. Why is it difficult to measure $R_n$ using electromagnetic probes (which couple primarily to charge), and how does parity-violating electron scattering circumvent this problem?

  3. $^{208}$Pb is the end point of three natural decay chains, but it is not the nucleus with the highest $B/A$. Explain why the decay chains terminate at $^{208}$Pb rather than continuing to a nucleus with higher binding energy per nucleon.

  4. Compare the single-particle spectra of $^{209}$Pb and $^{207}$Pb. Why does the particle-hole symmetry work so well — that is, why are the hole states in $^{207}$Pb well-described as single holes in the $^{208}$Pb core? Under what circumstances would this description break down?