Chapter 2 — Nuclear Properties: Size, Shape, Mass, Spin, and Moments

"The nucleus is not merely small — it is precisely small, and the precision with which we know its properties is one of the triumphs of twentieth-century experimental physics."

In Chapter 1, we discovered the nucleus and built the chart of nuclides as our organizing framework. We measured the binding energy curve and saw that the nuclear density is approximately constant — a remarkable fact that demands explanation. Now we must go further. A nucleus is not just a point with a mass. It has a finite size, a definite shape, an angular momentum, and electromagnetic moments that encode its internal structure. This chapter is a systematic survey of these measurable properties — the empirical foundation on which every nuclear model in this book will be tested.

The approach throughout is the same: we will see what experiments measure, write down the simplest theoretical expectation, compare the two, and pay close attention to where they agree and where they disagree. The disagreements are not failures — they are the signals of deeper physics.

We organize the chapter around six measurable properties: size (Section 2.1), mass (Section 2.2), spin and parity (Section 2.3), magnetic dipole moments (Section 2.4), electric quadrupole moments (Section 2.5), and isospin (Section 2.6). Each section follows the same pattern: physical motivation, mathematical setup, key results, comparison to data, and interpretation. By the end, you will have a working vocabulary of nuclear observables and — equally important — a calibrated sense of when simple models work and when they fail.

2.1 Nuclear Sizes and Charge Distributions

2.1.1 The Radius Formula and Density Saturation

The most fundamental fact about nuclear sizes was already apparent from Rutherford's scattering experiments and was quantified definitively by Robert Hofstadter's electron scattering measurements at Stanford in the 1950s (for which he shared the 1961 Nobel Prize): the nuclear radius scales as

$$R = r_0 A^{1/3}$$

where $A$ is the mass number and $r_0$ is a constant. This formula is the direct consequence of nuclear density saturation — the nuclear interior has approximately the same density for all nuclei from ${}^{12}\text{C}$ to ${}^{208}\text{Pb}$ and beyond. To see why, assume the nucleus is a uniform sphere of radius $R$ and density $\rho_0$:

$$A \cdot m_N = \rho_0 \cdot \frac{4}{3}\pi R^3$$

where $m_N$ is the nucleon mass. Solving for $R$:

$$R = \left(\frac{3 m_N}{4\pi \rho_0}\right)^{1/3} A^{1/3} = r_0 A^{1/3}$$

The constant $r_0$ therefore encodes the saturation density:

$$\rho_0 = \frac{3 m_N}{4\pi r_0^3}$$

Inserting $r_0 \approx 1.2 \text{ fm}$ gives $\rho_0 \approx 0.17 \text{ nucleons/fm}^3$, or equivalently about $2.3 \times 10^{17} \text{ kg/m}^3$ — roughly $2 \times 10^{14}$ times the density of water. This is the same density found inside neutron stars (Chapter 25), not because nuclear matter is a neutron star, but because the underlying saturation mechanism — the short-range repulsive core of the nuclear force — operates in both systems.

The value of $r_0$ depends slightly on what one means by "radius" and which experimental method is used. The following table summarizes the situation:

The differences are not contradictions — they reflect the fact that "the nuclear radius" is not a single number. The half-density radius $R_{1/2}$ and the rms radius $R_\text{rms}$ are related by the shape of the density profile, and the matter radius includes the neutron distribution while the charge radius does not. Throughout this book, we will use $r_0 \approx 1.2$ fm for the half-density radius unless stated otherwise, and $r_0 \approx 0.95$ fm when referring to the rms charge radius.

2.1.2 What Does "Radius" Mean? Different Probes, Different Answers

A nucleus is not a billiard ball with a sharp edge. The nuclear density distribution is better described by a Fermi (Woods-Saxon) function:

$$\rho(r) = \frac{\rho_0}{1 + \exp\left(\frac{r - R_{1/2}}{a}\right)}$$

where $R_{1/2}$ is the half-density radius (the radius at which $\rho = \rho_0/2$) and $a \approx 0.54$ fm is the diffuseness parameter characterizing the surface thickness. The surface thickness $t = 4a \ln 3 \approx 2.4$ fm is essentially constant across the periodic table. This is a second remarkable fact: all nuclei have the same surface thickness.

Different experimental probes measure different aspects of this distribution:

Elastic electron scattering measures the charge radius $R_\text{ch}$, because electrons interact via the electromagnetic force and see only the proton distribution (plus meson-exchange currents and the proton's own charge radius). The differential cross section is:

$$\frac{d\sigma}{d\Omega} = \left(\frac{d\sigma}{d\Omega}\right)_\text{Mott} |F(q)|^2$$

where the form factor $F(q)$ is the Fourier transform of the charge distribution:

$$F(q) = \frac{1}{Ze}\int \rho_\text{ch}(\mathbf{r})\, e^{i\mathbf{q}\cdot\mathbf{r}}\, d^3r$$

Here $q = |\mathbf{q}|$ is the momentum transfer. For a uniform sphere of radius $R$, the form factor is:

$$F(q) = \frac{3[\sin(qR) - qR\cos(qR)]}{(qR)^3}$$

which has zeros (diffraction minima) at $qR \approx 4.49, 7.73, \ldots$ The positions of these minima directly determine $R$. For ${}^{208}\text{Pb}$, electron scattering gives $R_\text{ch} = 5.501 \pm 0.005$ fm, consistent with $r_0 \approx 1.21$ fm when we account for the finite proton charge radius ($r_p \approx 0.84$ fm):

$$R_\text{ch}^2 = R_\text{point}^2 + r_p^2 + \frac{N}{Z}r_n^2 + \text{spin-orbit corrections}$$

where $R_\text{point}$ is the radius of the point-proton distribution and $r_n^2 \approx -0.12$ fm$^2$ accounts for the neutron's internal charge distribution.

Hadronic probes — proton elastic scattering, pion scattering, antiprotonic atoms — measure the matter radius $R_m$, which includes both protons and neutrons. For a neutron-rich nucleus, the matter radius exceeds the charge radius because the neutron distribution extends further. The difference defines the neutron skin thickness:

$$\Delta r_{np} = R_n - R_p$$

For ${}^{208}\text{Pb}$, the neutron skin has been the subject of intense experimental effort. The PREX-2 experiment at Jefferson Lab used parity-violating electron scattering — exploiting the fact that the weak force couples preferentially to neutrons via $Z^0$ exchange — to measure $\Delta r_{np} = 0.283 \pm 0.071$ fm. This result, published in 2021, has profound implications for the equation of state of neutron-rich matter and hence for neutron star radii (Chapter 25).

Muonic atoms provide a complementary determination of charge radii. A muon in a 1$s$ orbit around a heavy nucleus has a Bohr radius $a_\mu = a_0 (m_e / m_\mu) \approx a_0/207$, which for heavy nuclei falls inside the nuclear charge distribution. The muonic x-ray transition energies are sensitive to the nuclear charge radius at the level of $10^{-3}$ fm. It was the muonic hydrogen measurement of the proton charge radius by the PSI group in 2010 that triggered the "proton radius puzzle" — a 5$\sigma$ discrepancy with electronic measurements that took a decade to resolve (and whose resolution favored the smaller muonic value, $r_p = 0.841$ fm).

Isotope shifts in optical spectroscopy measure changes in the mean-square charge radius $\delta\langle r^2 \rangle$ between isotopes. Modern laser spectroscopy at facilities like ISOLDE (CERN) and TRIUMF can measure charge radii of nuclei far from stability, including short-lived isotopes produced at radioactive beam facilities. The odd-even staggering of charge radii — where the radius jumps at certain neutron numbers — provides direct evidence for shell structure (Chapter 6).

The isotope shift in an optical transition has two components: the mass shift (due to the change in reduced mass when neutrons are added) and the field shift (due to the change in the nuclear charge distribution). For heavy nuclei, the field shift dominates and is proportional to $\delta\langle r^2 \rangle$:

$$\delta\nu_\text{field} \approx \frac{2\pi Z e^2}{3} |\psi(0)|^2 \delta\langle r^2 \rangle$$

where $|\psi(0)|^2$ is the electron density at the nucleus. The collinear laser spectroscopy technique, in which a laser beam overlaps with a fast ion beam, achieves precisions of $\delta\langle r^2 \rangle \sim 10^{-3}$ fm$^2$. This has enabled systematic mapping of charge radii across long isotope chains, from the proton drip line to far beyond the valley of stability. The copper isotopes ($Z = 29$), for example, have been measured from ${}^{57}\text{Cu}$ to ${}^{78}\text{Cu}$ — a span of 21 neutrons — revealing the onset of deformation near $N = 40$ and the restoration of sphericity near the $N = 50$ shell closure.

💡 Concept Check: If ${}^{40}\text{Ca}$ has a charge radius of 3.478 fm and ${}^{48}\text{Ca}$ has a charge radius of 3.477 fm, what does this near-equality tell you? (Answer: Despite adding 8 neutrons, the proton distribution barely changes. The extra neutrons form a neutron skin. The shell closure at $N = 28$ keeps the proton core compact.)

⚠️ Subtlety: The near-equality of the ${}^{40}\text{Ca}$ and ${}^{48}\text{Ca}$ charge radii is actually somewhat surprising — most mass models predicted that adding 8 neutrons to a doubly magic nucleus would increase the charge radius slightly through the proton-neutron interaction. The measurement constrains nuclear forces in a neutron-rich environment and has become a benchmark for ab initio nuclear structure calculations.

2.1.3 Numerical Examples with Real Data

Let us compute the expected radii and compare to experiment.

Example 2.1: For ${}^{40}\text{Ca}$ ($Z = 20$, $N = 20$): $$R = 1.21 \times 40^{1/3} = 1.21 \times 3.420 = 4.14 \text{ fm}$$

The experimental charge radius is $R_\text{ch} = 3.478$ fm. The point-proton radius is approximately $R_p \approx \sqrt{R_\text{ch}^2 - r_p^2} = \sqrt{3.478^2 - 0.84^2} \approx 3.38$ fm. The discrepancy between 4.14 fm and 3.38 fm tells us that $r_0 = 1.21$ fm is the half-density radius of the matter distribution, while the charge radius is systematically smaller because the Fermi surface profile reduces the rms radius relative to the uniform-sphere value. For a Fermi distribution, the rms radius is:

$$R_\text{rms} = \sqrt{\frac{3}{5}}\, R_{1/2}\sqrt{1 + \frac{7\pi^2 a^2}{3 R_{1/2}^2}}$$

which for $R_{1/2} = 1.21 \times 40^{1/3} = 4.14$ fm and $a = 0.54$ fm gives $R_\text{rms} = 3.37$ fm — in excellent agreement with the measured value.

Example 2.2: For ${}^{208}\text{Pb}$ ($Z = 82$, $N = 126$): $$R_{1/2} = 1.21 \times 208^{1/3} = 1.21 \times 5.925 = 7.17 \text{ fm}$$ $$R_\text{rms} = \sqrt{\frac{3}{5}} \times 7.17 \times \sqrt{1 + \frac{7\pi^2(0.54)^2}{3(7.17)^2}} = 5.52 \text{ fm}$$

Experimental value: $R_\text{ch} = 5.501$ fm. The agreement to 0.3% is a testament to the quality of the $R = r_0 A^{1/3}$ parameterization with the correct density profile.

Let us also derive the rms radius formula from the Fermi distribution to make the connection explicit. The mean-square radius of a charge distribution $\rho(r)$ is:

$$\langle r^2 \rangle = \frac{\int_0^\infty r^2 \rho(r) 4\pi r^2 \, dr}{\int_0^\infty \rho(r) 4\pi r^2 \, dr}$$

For the Fermi distribution, this integral can be evaluated analytically in the limit $a \ll R_{1/2}$ (thin-skin approximation) using the Sommerfeld expansion:

$$\langle r^2 \rangle = \frac{3}{5}R_{1/2}^2 + \frac{7}{5}\pi^2 a^2$$

from which:

$$R_\text{rms} = \sqrt{\langle r^2 \rangle} = \sqrt{\frac{3}{5}}R_{1/2}\sqrt{1 + \frac{7\pi^2 a^2}{3 R_{1/2}^2}}$$

The correction term $7\pi^2 a^2/(3R_{1/2}^2)$ is small for heavy nuclei (about 3% for ${}^{208}\text{Pb}$) but significant for light nuclei (about 10% for ${}^{12}\text{C}$). This is why the simple formula $R_\text{rms} = \sqrt{3/5}\, r_0 A^{1/3}$ without the diffuseness correction works poorly for light nuclei — the surface is a larger fraction of the total volume.

2.2 Nuclear Masses and Binding Energies

2.2.1 Mass-Energy Equivalence and Nuclear Binding

The mass of a nucleus is not the sum of its constituent masses. The mass defect

$$\Delta M = Z m_p + N m_n - M(A,Z)$$

represents the mass equivalent of the binding energy:

$$B(A,Z) = \Delta M \cdot c^2 = [Z m_p + N m_n - M(A,Z)] c^2$$

In practice, nuclear masses are tabulated as the atomic mass $M_\text{atom}$, which includes the electron masses but (to high accuracy) also the electron binding energies. The atomic mass excess is defined as:

$$\Delta = M_\text{atom} - A \cdot u$$

where $u = 931.494\,102\,42(28)$ MeV/$c^2$ is the unified atomic mass unit, defined as $1/12$ the mass of a ${}^{12}\text{C}$ atom. Mass excesses are more convenient than full masses because they are small numbers (typically $-100$ to $+100$ MeV), making differences — which are the physically important quantities — easy to compute.

The connection between atomic masses and nuclear binding energies deserves a careful derivation. Starting from the atomic mass $M(A,Z)$ (which includes $Z$ electrons):

$$B(A,Z) = \left[Z M({}^{1}\text{H}) + N m_n - M(A,Z)\right] c^2$$

In terms of mass excesses, since $M = A \cdot u + \Delta/c^2$:

$$B(A,Z) = \left[Z \Delta_H + N \Delta_n - \Delta(A,Z)\right]$$

where $\Delta_H = 7288.971$ keV is the mass excess of the hydrogen atom and $\Delta_n = 8071.317$ keV is the mass excess of the neutron. This formula is the workhorse for computing binding energies from the AME tables.

Example 2.3: For ${}^{4}\text{He}$: $$B = 2(7288.971) + 2(8071.317) - 2424.916 = 30\,295.66 \text{ keV} = 28.296 \text{ MeV}$$ $$B/A = 7.074 \text{ MeV}$$

This is the binding energy of the alpha particle — one of the most tightly bound light nuclei, a fact that has consequences for alpha decay (Chapter 13) and stellar nucleosynthesis (Chapter 22).

2.2.2 Mass Spectrometry: From Aston to Penning Traps

The history of nuclear mass measurements is a history of increasing precision, and each order of magnitude in precision has opened new physics.

Magnetic deflection mass spectrometers (Aston, 1919; Dempster, 1935; Nier, 1940s) achieve $\delta m / m \sim 10^{-6}$. The basic principle is simple: ions with charge $q$ and mass $m$ moving at velocity $v$ perpendicular to a uniform magnetic field $B$ follow circular orbits with radius $\rho = mv/(qB)$. By measuring $\rho$ (or equivalently the position on a photographic plate or detector), one determines $m/q$. These instruments established the existence of isotopes and provided the first mass table.

Doublet measurements improved precision to $\delta m / m \sim 10^{-8}$ by comparing nearly degenerate mass multiplets — for example, comparing ${}^{12}\text{C}{}^{1}\text{H}_4$ with ${}^{16}\text{O}$ (both nominally 16 u). The mass difference is small and can be measured very precisely.

Penning traps represent the modern gold standard. A Penning trap confines a single ion using a superposition of a uniform magnetic field $B$ (providing radial confinement via the cyclotron motion) and a quadrupole electrostatic field (providing axial confinement). The ion executes three independent oscillations:

1. Axial oscillation at frequency $\nu_z = \frac{1}{2\pi}\sqrt{\frac{qU_0}{md^2}}$, where $U_0$ is the trapping voltage and $d$ is a characteristic trap dimension.

2. Modified cyclotron motion at frequency $\nu_+$.

3. Magnetron drift at frequency $\nu_-$.

The free-space cyclotron frequency is recovered via the invariance theorem:

$$\nu_c = \frac{1}{2\pi}\frac{qB}{m} = \sqrt{\nu_+^2 + \nu_z^2 + \nu_-^2}$$

or, in practice, via the relation $\nu_c = \nu_+ + \nu_-$. By measuring $\nu_c$ for the ion of interest and a reference ion of known mass, the mass ratio is determined:

$$\frac{m_\text{ion}}{m_\text{ref}} = \frac{\nu_{c,\text{ref}}}{\nu_{c,\text{ion}}}$$

Modern Penning traps (ISOLTRAP at CERN, JYFLTRAP at Jyv\"askyl\"a, CPT at Argonne, LEBIT at MSU/FRIB) achieve $\delta m / m \sim 10^{-9}$ to $10^{-11}$, corresponding to absolute uncertainties of less than 1 keV even for heavy nuclei. The JYFLTRAP measurement of the ${}^{76}\text{Ge}$-${}^{76}\text{Se}$ $Q$-value for double beta decay — critical for determining the neutrino mass from neutrinoless double beta decay experiments — has an uncertainty of $0.8$ keV.

This precision matters because many important quantities are differences of large masses:

• Q-values for reactions and decays: $Q = (M_\text{initial} - M_\text{final})c^2$
• Nucleon separation energies: $S_n = B(A,Z) - B(A-1,Z)$, $S_p = B(A,Z) - B(A,Z-1)$
• Pairing energies: $\Delta_n = \frac{1}{2}[S_n(A,Z) - S_n(A-1,Z)]$ (for even $A$)

A 10 keV error in the mass of a single nuclide can produce a 10 keV error in a Q-value, which can change a predicted half-life for a beta decay near the endpoint by an order of magnitude, or shift the predicted path of the astrophysical r-process (Chapter 23).

2.2.3 The Atomic Mass Evaluation

The Atomic Mass Evaluation (AME) is the community standard for nuclear masses, maintained by the collaboration of M. Wang, W.J. Huang, F.G. Kondev, G. Audi, and S. Naimi. The most recent published evaluation is AME2020, which includes:

• 2,457 experimentally measured masses (out of ~3,400 known nuclides)
• Estimated masses for ~900 additional nuclides from systematic trends
• All data available from the IAEA Nuclear Data Services

The AME performs a least-squares adjustment over all input data — direct mass measurements, decay Q-values, reaction Q-values — to produce a self-consistent mass table. This is a nontrivial global fit because many masses are connected by multiple independent measurements. For example, the mass of ${}^{56}\text{Ni}$ is constrained by Penning trap measurements, by the ${}^{56}\text{Co}(\beta^+){}^{56}\text{Ni}$ Q-value, and by the ${}^{58}\text{Ni}(p,t){}^{56}\text{Ni}$ reaction Q-value, among others.

Selected mass excesses from AME2020 (in keV):

Separation energy systematics reveal shell structure with crystal clarity. A plot of $S_{2n}$ (two-neutron separation energy) vs. $N$ shows sudden drops at $N = 28, 50, 82, 126$ — the magic numbers. We will derive the theoretical origin of these numbers in Chapter 6, but the experimental signature is unambiguous: at a magic number, the last two neutrons are bound by $\sim 4$–6 MeV less than the preceding pair. This is the mass-table fingerprint of quantum shell closure.

📊 Data Point: The two-neutron separation energy of ${}^{132}\text{Sn}$ ($N = 82$) is $S_{2n} = 12.59$ MeV. For ${}^{134}\text{Sn}$ ($N = 84$), it drops to $S_{2n} = 7.49$ MeV — a decrease of more than 5 MeV in a single step. This is one of the sharpest shell closures observed.

2.2.4 Separation Energies and the Drip Lines

The neutron drip line is the boundary in the chart of nuclides beyond which adding another neutron produces an unbound system: $S_n \leq 0$. Similarly, the proton drip line is defined by $S_p \leq 0$. Beyond the drip lines, nuclei exist only as short-lived resonances that decay by nucleon emission on timescales of $\sim 10^{-22}$ s.

The neutron drip line has been experimentally reached only for light elements ($Z \leq 10$), most recently with the observation of ${}^{31}\text{F}$ at RIKEN in 2019. For heavier elements, the neutron drip line is predicted by mass models but lies far from the valley of stability. The proton drip line is known up to about $Z = 83$ because fewer neutrons are needed to reach the boundary — the Coulomb force assists in unbinding the last proton.

The asymmetry between the two drip lines is a direct consequence of the Coulomb interaction. For a nucleus at the neutron drip line, adding one more neutron costs no Coulomb energy — only the nuclear surface and asymmetry energies matter. The neutron drip line lies far from stability (at very high $N/Z$ ratios), making the relevant nuclei difficult to produce and short-lived. For the proton drip line, the Coulomb barrier provides an additional energy cost for each added proton, and the drip line occurs closer to stability. Furthermore, proton-unbound nuclei can still have measurable half-lives ($\sim$ms to $\sim$s) because the proton must tunnel through the Coulomb barrier to escape — the same tunneling physics that governs alpha decay (Chapter 13).

Two-proton radioactivity is a particularly exotic decay mode observed near the proton drip line. In nuclei like ${}^{45}\text{Fe}$, ${}^{48}\text{Ni}$, and ${}^{54}\text{Zn}$, the one-proton separation energy is positive ($S_p > 0$, so single-proton emission is energetically forbidden) but the two-proton separation energy is negative ($S_{2p} < 0$), and the nucleus decays by emitting two protons simultaneously. This is a genuinely quantum-mechanical process — the two protons are emitted in a correlated pair, not sequentially — and was first observed at GSI in 2002.

🔗 Cross-Reference: We will use these separation energy systematics extensively in Chapter 4 when fitting the semi-empirical mass formula, and in Chapter 23 when tracing the r-process nucleosynthesis path along the neutron drip line.

💡 Concept Check: Why do we plot $S_{2n}$ rather than $S_n$ when looking for shell closure signatures? (Answer: The odd-even staggering of $S_n$ — caused by the pairing force, which makes even-$N$ nuclei more tightly bound — obscures the shell-closure signal. Plotting $S_{2n}$ eliminates the pairing oscillation and reveals the shell gaps cleanly.)

2.3 Nuclear Spin and Parity

2.3.1 Definitions and Notation

Every nucleus in its ground state has a definite total angular momentum $J$ and parity $\pi$, written together as $J^\pi$. The total angular momentum is the vector sum of the orbital angular momenta and intrinsic spins of all $A$ nucleons:

$$\mathbf{J} = \sum_{i=1}^{A} (\boldsymbol{\ell}_i + \mathbf{s}_i)$$

Since angular momentum is quantized, $J$ takes integer or half-integer values, and the parity $\pi = \pm 1$ is the eigenvalue of the parity operator $\hat{P}$ acting on the nuclear wave function:

$$\hat{P} \Psi = \pi \Psi = \prod_{i=1}^{A} (-1)^{\ell_i} \Psi$$

The nuclear spin $J$ (often loosely called "the nuclear spin," though strictly it is the total angular momentum quantum number) is always:

• Integer for even-$A$ nuclei (even number of fermions → boson-like composite)
• Half-integer for odd-$A$ nuclei (odd number of fermions → fermion-like composite)

2.3.2 Empirical Rules

The following empirical rules capture the dominant pattern:

1. Even-even nuclei ($Z$ even, $N$ even): The ground state is always $J^\pi = 0^+$. This is one of the most robust empirical facts in nuclear physics — there is no known exception among the ~800 even-even nuclei with measured ground-state spins. The pairing interaction (Chapter 7) drives nucleons to couple in pairs with $J = 0$.

2. Odd-$A$ nuclei (one odd, one even): The ground-state spin and parity are determined by the single unpaired nucleon. If the unpaired nucleon occupies a shell-model orbit with quantum numbers $(\ell, j)$, then $J^\pi = j^{(-1)^\ell}$.

3. Odd-odd nuclei ($Z$ odd, $N$ odd): The ground-state spin results from coupling the unpaired proton (with $j_p$) and unpaired neutron (with $j_n$). The Nordheim rules give guidance: if $j_p = \ell_p \pm 1/2$ and $j_n = \ell_n \pm 1/2$, then the strong rule (both $j = \ell + 1/2$ or both $j = \ell - 1/2$) predicts $J = |j_p - j_n|$, while the weak rule (one of each) gives $J = |j_p - j_n|$ or $j_p + j_n$ — and the weak rule frequently fails.

Example 2.3: ${}^{17}\text{O}$ ($Z = 8$, $N = 9$). Oxygen-17 has an even number of protons (all paired) and one unpaired neutron. The 9th neutron fills the $1d_{5/2}$ orbit: $\ell = 2$, $j = 5/2$, parity $= (-1)^2 = +$. The predicted ground state is $J^\pi = 5/2^+$. The experimental value is $5/2^+$. ✓

Example 2.4: ${}^{207}\text{Pb}$ ($Z = 82$, $N = 125$). Lead-207 has a doubly magic proton core ($Z = 82$) and one neutron hole in the $N = 126$ shell closure. The 126th neutron orbit is $1i_{13/2}$: missing one neutron from this filled shell is equivalent to one neutron hole with $\ell = 6$, $j = 13/2$, parity $= (-1)^6 = +$. Wait — that would give $13/2^+$. But the measured ground state is $J^\pi = 1/2^-$, corresponding to the $3p_{1/2}$ orbit ($\ell = 1$, $j = 1/2$). This tells us that the last neutron below the $N = 126$ gap is not in $1i_{13/2}$ but in $3p_{1/2}$ — the shell-model ordering near $N = 126$ places $3p_{1/2}$ just below the gap. The point is not that the rule fails, but that applying it correctly requires knowing the shell-model level ordering, which we will derive in Chapter 6.

2.3.3 The Pairing Interaction and $0^+$ Ground States

The universal $J^\pi = 0^+$ ground state for even-even nuclei demands explanation. The key is the pairing interaction — a short-range, attractive component of the nuclear force that preferentially couples pairs of like nucleons (two protons or two neutrons) in time-reversed orbits to total angular momentum $J = 0$. This is analogous to Cooper pairing in superconductors, and the analogy is not superficial: the nuclear pairing gap $\Delta \approx 12/\sqrt{A}$ MeV behaves much like the BCS gap in condensed-matter physics (Chapter 7).

The pairing energy manifests in several observable ways:

• Odd-even mass staggering: Even-$N$ (or even-$Z$) nuclei are systematically more bound than their odd-$N$ (odd-$Z$) neighbors, by an amount $\delta \approx 12/\sqrt{A}$ MeV. This is visible as a zigzag pattern in plots of $B(A,Z)$ vs. $N$ at fixed $Z$.

• The $E(2^+_1)$ energy gap: In even-even nuclei, the first excited state (typically $2^+$) lies at an energy that reflects the cost of breaking a pair. Near magic numbers, $E(2^+_1) \sim 1$–4 MeV; midshell, collective effects lower it to $\sim 80$–300 keV.

• Moments of inertia: The nuclear moment of inertia is smaller than the rigid-body value, because pairing correlations reduce the response to rotation — the system is partly "superfluid."

2.3.4 How Spin and Parity Are Measured

Nuclear spins and parities are determined from multiple experimental techniques, each with its own strengths:

• Atomic hyperfine structure: The nuclear spin $I$ (using the spectroscopic notation $I$ instead of $J$) splits atomic energy levels into $2I + 1$ components (for $J_\text{atom} > I$) or $2J_\text{atom} + 1$ components (for $I > J_\text{atom}$). Counting the number of hyperfine components gives $I$. The ratio of the hyperfine splitting between two atomic levels gives the nuclear magnetic moment (Section 2.4).

• Nuclear reaction angular distributions: The angular distribution of products from reactions like $(d,p)$ or $(e,e')$ encodes the orbital angular momentum $\ell$ transferred, which constrains $J^\pi$. In a $(d,p)$ stripping reaction, the angular distribution of the outgoing proton peaks at an angle determined by $\ell$, allowing model-independent $\ell$ assignments. Combined with the known spins of the target and projectile, this gives $J^\pi$ for the final state.

• Gamma-ray angular correlations and polarization: The multipolarity and mixing ratios of gamma transitions between states constrain the spins and parities of both states (Chapter 9). Modern gamma-ray tracking arrays (GRETINA at FRIB, AGATA in Europe) measure both the energy and the direction of gamma rays with sufficient precision to determine spin-parity assignments for states populated in in-beam spectroscopy.

• Beta-decay $\log ft$ values: Allowed beta decays ($\Delta J = 0, 1$; no parity change) have $\log ft \approx 3$–6, while forbidden transitions have progressively larger $\log ft$ values, constraining the spin-parity change (Chapter 14). A measured $\log ft = 3.5$ essentially proves an allowed Gamow-Teller transition ($\Delta J = 0, 1$, no parity change), providing a strong constraint on the spin-parity of the daughter state.

• Coulomb excitation: The excitation of nuclear states by the electric field of a passing heavy ion is highly selective: $E2$ Coulomb excitation populates $2^+$ states from $0^+$ ground states, $E1$ excitation populates $1^-$ states, and so on. The multipolarity selectivity makes this a clean tool for spin-parity determination.

⚠️ Common Misconception: Students often confuse nuclear spin $I$ (or $J$) with the spin quantum number $s = 1/2$ of a single nucleon. The nuclear spin is the total angular momentum of the entire nucleus, which can be 0, 1/2, 1, 3/2, ... up to values exceeding 30 in some isomeric states (such as the $J^\pi = 35/2^-$ isomer in ${}^{179}\text{Hf}$, which has a half-life of 25 days). The single-nucleon spin $s = 1/2$ is just one ingredient.

2.4 Magnetic Dipole Moments

2.4.1 The Nuclear Magneton

A nucleon with angular momentum $\mathbf{j} = \boldsymbol{\ell} + \mathbf{s}$ has a magnetic moment arising from both its orbital motion (for protons) and its intrinsic spin (for both protons and neutrons):

$$\hat{\boldsymbol{\mu}} = g_\ell \,\mu_N \hat{\boldsymbol{\ell}} + g_s \,\mu_N \hat{\mathbf{s}}$$

where $\mu_N = e\hbar/(2m_p) = 5.050\,783\,7461(15) \times 10^{-27}$ J/T is the nuclear magneton — the natural unit for nuclear magnetic moments, analogous to the Bohr magneton for electrons but smaller by a factor of $m_e/m_p \approx 1/1836$.

The $g$-factors for a free nucleon are:

The proton orbital $g$-factor is 1 because the proton carries one unit of charge. The neutron orbital $g$-factor is 0 because the neutron is uncharged — it generates no magnetic moment from its orbital motion. The spin $g$-factors are the famous anomalous magnetic moments, which were among the first indications that nucleons are composite particles (Chapter 31). For a structureless Dirac particle, $g_s$ would be exactly 2; the departures arise from the quark substructure and the pion cloud.

The fact that the neutron has a large anomalous magnetic moment ($\mu_n = g_s^n \mu_N / 2 = -1.913\,\mu_N$) despite being electrically neutral is itself remarkable. A truly elementary neutral particle would have $\mu = 0$. The nonzero neutron moment proves that the neutron has internal structure — a distribution of charged quarks ($udd$) whose orbital and spin angular momenta produce a net magnetic moment. The negative sign means the magnetic moment points opposite to the spin, consistent with the outer region of the neutron being dominated by negatively charged $d$ quarks.

The proton magnetic moment is $\mu_p = g_s^p \mu_N / 2 = +2.793\,\mu_N$. The sum $\mu_p + \mu_n = 0.880\,\mu_N$ (the isoscalar moment) and the difference $\mu_p - \mu_n = 4.706\,\mu_N$ (the isovector moment) are important diagnostic quantities that constrain quark models and have implications for nuclear magnetic moments of mirror nuclei.

2.4.2 Schmidt Values: The Single-Particle Prediction

For an odd-$A$ nucleus, the magnetic moment is determined (in the single-particle picture) by the unpaired nucleon. If this nucleon has angular momentum quantum numbers $\ell$ and $j$, the magnetic moment (defined as the expectation value of $\hat{\mu}_z$ in the state $m_j = j$) is:

$$\mu = \langle j, m_j = j | \hat{\mu}_z | j, m_j = j \rangle$$

Using the projection theorem (a consequence of the Wigner-Eckart theorem applied to vector operators), the expectation value of any vector operator $\hat{V}$ in a state $|j, m\rangle$ is proportional to $\langle \hat{J}\rangle$:

$$\langle j, m | \hat{V}_z | j, m \rangle = \frac{\langle j, m | \hat{\mathbf{V}} \cdot \hat{\mathbf{J}} | j, m \rangle}{j(j+1)} \cdot m$$

Applying this to the magnetic moment operator $\hat{\boldsymbol{\mu}} = g_\ell \mu_N \hat{\boldsymbol{\ell}} + g_s \mu_N \hat{\mathbf{s}}$ and evaluating in the stretched state $m_j = j$:

$$\mu = g_j \, j \, \mu_N$$

where the Landé $g$-factor is derived from $\langle \hat{\boldsymbol{\mu}} \cdot \hat{\mathbf{J}} \rangle = g_\ell \mu_N \langle \hat{\boldsymbol{\ell}} \cdot \hat{\mathbf{J}} \rangle + g_s \mu_N \langle \hat{\mathbf{s}} \cdot \hat{\mathbf{J}} \rangle$. Using $\hat{\mathbf{J}} = \hat{\boldsymbol{\ell}} + \hat{\mathbf{s}}$ to write $\hat{\boldsymbol{\ell}} \cdot \hat{\mathbf{J}} = (\hat{\mathbf{J}}^2 + \hat{\boldsymbol{\ell}}^2 - \hat{\mathbf{s}}^2)/2$ and similarly for $\hat{\mathbf{s}} \cdot \hat{\mathbf{J}}$:

$$g_j = g_\ell \frac{j(j+1) + \ell(\ell+1) - s(s+1)}{2j(j+1)} + g_s \frac{j(j+1) + s(s+1) - \ell(\ell+1)}{2j(j+1)}$$

This is the nuclear analog of the atomic Landé formula, but with the free-nucleon $g$-factors instead of $g_\ell = 1$ and $g_s = 2$. The key difference from the atomic case is the large anomalous moments of the nucleon: $g_s^p = 5.586 \neq 2$ and $g_s^n = -3.826 \neq 0$.

For $s = 1/2$, the two cases $j = \ell + 1/2$ and $j = \ell - 1/2$ give the Schmidt values:

Case 1: $j = \ell + 1/2$ $$\mu_\text{Schmidt} = \left(j - \frac{1}{2} + g_s/2\right) \mu_N$$

which simplifies to: $$\mu_\text{Schmidt} = \left(\ell + g_s/2\right) \mu_N \quad \text{(proton, } j = \ell + 1/2\text{)}$$

Case 2: $j = \ell - 1/2$ $$\mu_\text{Schmidt} = \frac{j}{j+1}\left(j + \frac{3}{2} - g_s/2\right) \mu_N$$

Example 2.5: ${}^{17}\text{O}$ ground state $5/2^+$, single neutron in $1d_{5/2}$ ($\ell = 2$, $j = 5/2$, $j = \ell + 1/2$):

$$\mu_\text{Schmidt} = \left(2 + \frac{-3.826}{2}\right)\mu_N = (2 - 1.913)\mu_N = 0.087\,\mu_N$$

Experimental value: $\mu = -1.894\,\mu_N$.

This is a spectacular failure — the Schmidt value has the wrong sign! This disagreement is one of the classic puzzles of nuclear physics and tells us that the single-particle picture is incomplete. Core polarization, meson exchange currents, and configuration mixing all contribute to the discrepancy.

Example 2.6: ${}^{209}\text{Bi}$ ground state $9/2^-$, single proton in $1h_{9/2}$ ($\ell = 5$, $j = 9/2$, $j = \ell - 1/2$):

$$g_j = g_\ell \frac{j(j+1) + \ell(\ell+1) - 3/4}{2j(j+1)} + g_s \frac{j(j+1) + 3/4 - \ell(\ell+1)}{2j(j+1)}$$

$$= 1 \cdot \frac{(9/2)(11/2) + 30 - 3/4}{2(9/2)(11/2)} + 5.586 \cdot \frac{(9/2)(11/2) + 3/4 - 30}{2(9/2)(11/2)}$$

$$= 1 \cdot \frac{99/4 + 30 - 3/4}{99/2} + 5.586 \cdot \frac{99/4 + 3/4 - 30}{99/2}$$

$$= 1 \cdot \frac{54}{49.5} + 5.586 \cdot \frac{-5}{49.5} = 1.091 - 0.564 = 0.527$$

Therefore:

$$\mu_\text{Schmidt} = g_j \cdot j \cdot \mu_N = 0.527 \times 4.5 \times \mu_N = 2.37\,\mu_N$$

(Using the simplified formula: $\mu_\text{Schmidt} = \frac{9/2}{11/2}(6 - 2.793)\mu_N = \frac{9}{11}(3.207)\mu_N = 2.624\,\mu_N$; the small numerical difference arises from rounding in the intermediate steps.)

Experimental value: $\mu = +4.111\,\mu_N$. Again, significant discrepancy — the Schmidt value is in the right ballpark and has the right sign, but underestimates the moment by about 1.5 $\mu_N$. The systematic underprediction for $j = \ell - 1/2$ states and overprediction for $j = \ell + 1/2$ states is a universal feature of the Schmidt model.

Example 2.7: Let us work through one more case that shows good agreement. ${}^{3}\text{H}$ (tritium), ground state $1/2^+$, with the unpaired proton in $1s_{1/2}$ ($\ell = 0$, $j = 1/2$, $j = \ell + 1/2$):

$$\mu_\text{Schmidt} = \left(0 + \frac{5.586}{2}\right)\mu_N = 2.793\,\mu_N$$

Experimental value: $\mu = +2.979\,\mu_N$. The agreement to about 6% reflects the fact that the $A = 3$ system has relatively few configurations available for mixing, making the single-particle picture a reasonable (though not exact) starting point.

2.4.3 The Schmidt Diagrams

When the measured magnetic moments of all odd-proton and odd-neutron nuclei are plotted against $j$, they cluster between the two Schmidt lines ($j = \ell + 1/2$ and $j = \ell - 1/2$) but rarely fall on them. This is the Schmidt diagram. The systematic pattern is:

• Moments for $j = \ell + 1/2$ states tend to lie below the upper Schmidt line.
• Moments for $j = \ell - 1/2$ states tend to lie above the lower Schmidt line.
• The data cluster at roughly 60–70% of the way from the center toward the Schmidt lines.

The physical reasons for the departure from Schmidt values include:

1. Configuration mixing: The ground state is not a pure single-particle state but includes admixtures of other configurations, each contributing to the magnetic moment.

2. Core polarization: The unpaired nucleon polarizes the even-even core, inducing a small magnetic moment in the core that adds to (or subtracts from) the single-particle moment.

3. Meson exchange currents: The nucleon-nucleon interaction is mediated by meson exchange, and the exchanged mesons carry charge and current, modifying the effective magnetic moments.

4. Quenching of $g_s$: Empirically, replacing the free-nucleon $g_s$ by an "effective" value $g_s^\text{eff} \approx 0.7 \, g_s^\text{free}$ brings the Schmidt values much closer to experiment. This quenching is a manifestation of the above effects combined.

📊 Data Point: The magnetic moment of ${}^{41}\text{Sc}$ ($7/2^-$, single proton in $1f_{7/2}$) is $\mu = +5.431\,\mu_N$. The Schmidt value is $\mu_\text{Schmidt} = 5.793\,\mu_N$ for $j = \ell + 1/2$. The agreement to ~6% is unusually good and reflects the fact that ${}^{41}\text{Sc}$ is a single proton outside the doubly magic ${}^{40}\text{Ca}$ core, making the single-particle picture a good approximation.

2.5 Electric Quadrupole Moments

2.5.1 Definition and Physical Meaning

The electric quadrupole moment $Q$ describes the departure of the nuclear charge distribution from spherical symmetry. For a nuclear state with angular momentum $J$ and maximum projection $m = J$, the spectroscopic quadrupole moment is:

$$eQ = \langle J, m = J | \hat{Q}_{20} | J, m = J \rangle = \langle J, m = J | \sum_i e_i(3z_i^2 - r_i^2) | J, m = J \rangle$$

where $e_i = e$ for protons and $e_i = 0$ for neutrons (in the simplest picture), and the sum runs over all protons. The factor of $e$ on the left side is conventional so that $Q$ has dimensions of area (typically reported in units of fm$^2$ or barns = 100 fm$^2$).

The sign convention is:

• $Q > 0$: prolate deformation (elongated along the symmetry axis, like a rugby ball)
• $Q < 0$: oblate deformation (flattened along the symmetry axis, like a doorknob)
• $Q = 0$: spherical, or $J = 0$ or $J = 1/2$ (for which $Q = 0$ trivially by angular momentum selection rules, regardless of deformation)

This last point is critical: a measured $Q = 0$ for a state with $J = 0$ or $J = 1/2$ tells you nothing about the intrinsic shape of the nucleus. The quadrupole moment vanishes because the wave function is isotropic in the lab frame, not because the nucleus is spherical. To determine the intrinsic shape, one must use other probes (rotational bands, Coulomb excitation — see Chapter 8).

2.5.2 Single-Particle Quadrupole Moments

For a single nucleon in an orbit with angular momentum $j$, the single-particle quadrupole moment is:

$$Q_\text{sp} = -e_\text{eff} \langle r^2 \rangle \frac{2j - 1}{2(j + 1)}$$

where $\langle r^2 \rangle$ is the mean-square radius of the orbit and $e_\text{eff}$ is the effective charge (equal to $e$ for a proton in the simplest approximation). Note the sign: for a single proton, $Q_\text{sp} < 0$ (oblate), because a single particle in a circular orbit creates a toroidal charge distribution concentrated in the equatorial plane.

For a single-particle estimate with $\langle r^2 \rangle \approx (3/5)R^2 = (3/5)(r_0 A^{1/3})^2$:

$$Q_\text{sp} \approx -\frac{3}{5}r_0^2 A^{2/3} \frac{2j - 1}{2(j+1)} \text{ fm}^2$$

Example 2.7: ${}^{209}\text{Bi}$ ($9/2^-$, single proton in $1h_{9/2}$):

$$Q_\text{sp} = -\frac{3}{5}(1.2)^2(209)^{2/3}\frac{2(9/2) - 1}{2(9/2 + 1)} = -0.864 \times 35.23 \times \frac{8}{11} = -22.1 \text{ fm}^2 = -0.22 \text{ b}$$

Experimental value: $Q = -0.516(15)$ b $= -51.6$ fm$^2$.

The measured quadrupole moment is more than twice as large (in magnitude) as the single-particle estimate. This enhancement is the hallmark of core polarization — the unpaired proton polarizes the even-even core, inducing a collective quadrupole deformation that amplifies the moment. This effect, first explained by Bohr and Mottelson, was a key insight leading to the unified model of nuclear structure (Chapter 8).

2.5.3 Intrinsic vs. Spectroscopic Quadrupole Moments

For nuclei that are deformed in their intrinsic (body-fixed) frame, the intrinsic quadrupole moment $Q_0$ is related to the measured spectroscopic moment $Q$ by:

$$Q = \frac{3K^2 - J(J+1)}{(J+1)(2J+3)} Q_0$$

where $K$ is the projection of $J$ on the nuclear symmetry axis. For the ground-state rotational band with $K = J$:

$$Q = \frac{3J^2 - J(J+1)}{(J+1)(2J+3)} Q_0 = \frac{J(2J-1)}{(J+1)(2J+3)} Q_0$$

For large $J$, $Q \to Q_0$. For $J = 0$ or $J = 1/2$, $Q = 0$ regardless of $Q_0$.

The intrinsic quadrupole moment is related to the deformation parameter $\beta_2$ by:

$$Q_0 = \frac{3}{\sqrt{5\pi}} Z R_0^2 \beta_2 (1 + 0.36 \beta_2 + \ldots)$$

where $R_0 = r_0 A^{1/3}$ and $\beta_2$ parameterizes the departure from sphericity in the expansion of the nuclear surface:

$$R(\theta, \phi) = R_0 \left[1 + \beta_2 Y_2^0(\theta)\right]$$

Typical deformation parameters range from $|\beta_2| \lesssim 0.05$ (near-spherical, e.g., nuclei near magic numbers) to $|\beta_2| \sim 0.3$–0.4 (well-deformed, e.g., rare-earth and actinide nuclei) to $|\beta_2| \sim 0.6$ (superdeformed, observed in high-spin states).

2.5.4 Systematic Trends

The quadrupole moment systematics across the chart of nuclides reveal a striking pattern:

• Near closed shells (magic numbers), $Q$ is small and consistent with single-particle values. ${}^{208}\text{Pb}$ is the paradigmatic case: $Q = 0$ (because $J^\pi = 0^+$), and the low-lying excitation spectrum shows no rotational bands. It is as spherical as a nucleus can be.

• Between closed shells (midshell), $|Q|$ can be enormous — up to $\sim 8$ b for rare-earth and actinide nuclei. For example, ${}^{176}\text{Lu}$ ($J^\pi = 7^-$) has $Q = +8.0$ b, corresponding to an intrinsic deformation $\beta_2 \approx 0.30$. The nucleus is roughly 30% longer than it is wide.

• The transition from spherical to deformed is often abrupt. In the samarium isotopes, $Q$ is near zero for ${}^{144}\text{Sm}$ ($N = 82$, magic) but jumps to $Q_0 \sim 6$ b by ${}^{154}\text{Sm}$ ($N = 92$), just 10 neutrons beyond the shell closure. This phase transition from spherical to deformed shapes is one of the most dramatic phenomena in nuclear structure (Chapter 8).

• ${}^{208}\text{Pb}$ is the paradigmatic spherical nucleus. With $J^\pi = 0^+$, its quadrupole moment is trivially zero. But the key evidence for sphericity is that its low-lying spectrum shows no rotational band structure. The first excited state is $3^-$ at 2.614 MeV (an octupole vibration), and the first $2^+$ state lies at 4.085 MeV — vastly higher than in deformed nuclei (where $E(2^+_1) \sim 50$–100 keV). The high $2^+$ energy reflects the stiffness of the doubly magic ${}^{208}\text{Pb}$ against quadrupole deformation.

The following data table summarizes measured spectroscopic quadrupole moments for a selection of nuclei spanning the chart, illustrating the range from near-spherical to strongly deformed:

The deuteron (${}^{2}\text{H}$) deserves special mention. Its quadrupole moment $Q = +0.00286$ b is small but nonzero. Since the deuteron is a $J = 1$ system, a nonzero $Q$ is allowed, and its positive sign indicates a slight prolate deformation. This tiny quadrupole moment was one of the first pieces of evidence that the nuclear force contains a tensor component — a force that depends on the orientation of the nucleon spins relative to the line connecting them. A purely central force would give $Q = 0$ for the deuteron. We will derive this in detail in Chapter 3.

💡 Key Insight: The electric quadrupole moment is perhaps the single most informative observable for nuclear structure. It bridges the single-particle picture (Chapter 6) and the collective picture (Chapter 8), and its failure in the single-particle model was one of the key motivations for the Bohr-Mottelson unified model.

2.6 Isospin

2.6.1 Motivation and Definition

The proton and neutron have nearly identical masses ($m_n - m_p = 1.293$ MeV/$c^2$, or 0.14%), similar sizes, and — crucially — the nuclear (strong) force between two protons, two neutrons, or a proton-neutron pair is essentially the same (after correcting for the Coulomb interaction). This charge independence of the nuclear force was recognized by Heisenberg in 1932, who proposed treating the proton and neutron as two states of a single particle, the nucleon, distinguished by a new quantum number called isospin.

Formally, isospin $\mathbf{T}$ is defined in exact analogy to ordinary spin:

• Each nucleon has isospin $t = 1/2$.
• The proton is the $t_3 = +1/2$ state: $|p\rangle = |t = 1/2, t_3 = +1/2\rangle$.
• The neutron is the $t_3 = -1/2$ state: $|n\rangle = |t = 1/2, t_3 = -1/2\rangle$.

(Note: the sign convention varies in the literature. The convention $t_3(p) = +1/2$ is used by most nuclear physics texts, including Krane and Heyde. The particle physics convention reverses the signs.)

For a nucleus with $Z$ protons and $N$ neutrons, the third component of total isospin is:

$$T_3 = \frac{1}{2}(Z - N)$$

The total isospin $T$ satisfies $T \geq |T_3|$. For ground states, the empirical rule is:

$$T = |T_3| = \frac{1}{2}|N - Z|$$

This means the ground state has the minimum value of $T$ consistent with the given $Z$ and $N$ — a consequence of the isospin-dependent part of the nuclear force being more attractive in the $T = 0$ channel than in the $T = 1$ channel (for $N = Z$ nuclei).

Example 2.9: Isospin assignments for selected nuclei:

Note that ${}^{14}\text{C}$, ${}^{14}\text{N}$ (in its $T = 1$ excited state at 2.31 MeV), and ${}^{14}\text{O}$ form a $T = 1$ isospin triplet. The ground state of ${}^{14}\text{N}$ has $T = 0$ because the $T = 0$ state is lower — the nuclear force is more attractive in the proton-neutron $T = 0$ channel. The energy of the $T = 1$ state in ${}^{14}\text{N}$ matches the ground-state energies of ${}^{14}\text{C}$ and ${}^{14}\text{O}$ after correcting for the Coulomb displacement. This is a textbook example of isobaric analog states.

2.6.2 Isospin Symmetry and Its Breaking

If the nuclear force were exactly charge-independent and the Coulomb force did not exist, isospin would be an exact quantum number — nuclear energy levels would form degenerate multiplets labeled by $T$, with $(2T + 1)$ members at the same energy. In reality, isospin symmetry is broken by:

1. The Coulomb interaction: This is the dominant source of isospin-symmetry breaking, contributing energy differences of order $\sim Z e^2 / R \sim$ several MeV for heavy nuclei.

2. The proton-neutron mass difference: $m_n - m_p = 1.293$ MeV/$c^2$, which contributes to isospin breaking at the percent level.

3. Charge-symmetry-breaking (CSB) and charge-independence-breaking (CIB) components of the nuclear force: These are smaller corrections arising from the $u$-$d$ quark mass difference and electromagnetic effects within QCD. They contribute at the level of $\sim 100$–300 keV.

Despite these symmetry-breaking effects, isospin remains an approximately good quantum number, especially for light nuclei. The evidence comes from the systematic study of isobaric analog states (IAS) — states in different nuclei that are members of the same isospin multiplet. The energies of IAS members can be related by the isobaric multiplet mass equation (IMME):

$$M(A, T, T_3) = a(A, T) + b(A, T) T_3 + c(A, T) T_3^2$$

where $a$, $b$, $c$ are constants for a given multiplet. The physical interpretation of these coefficients is illuminating:

• $a$: The charge-independent average mass. This is the "nuclear physics" part — what the mass would be if the Coulomb force did not exist and protons and neutrons were identical.
• $b$: The linear Coulomb displacement. This is proportional to $T_3$ because the Coulomb energy increases linearly with the number of protons (approximately). It also includes the $n$-$p$ mass difference contribution.
• $c$: The quadratic Coulomb term. This arises from the fact that the Coulomb energy goes as $Z(Z-1)/R \propto (A/2 + T_3)(A/2 + T_3 - 1)/R$, which has a quadratic term in $T_3$.

The quadratic form works remarkably well — cubic ($d T_3^3$) and quartic terms are typically less than 10 keV, demonstrating that isospin breaking is dominantly a one- and two-body (Coulomb) effect. The precision of the IMME has been tested for hundreds of isospin multiplets, and the handful of cases where the cubic term might be nonzero (notably the $A = 32$, $T = 2$ quintet) are the subject of active experimental investigation at facilities like TRIUMF and GANIL.

Numerical example: For the $A = 14$, $T = 1$ triplet (${}^{14}\text{C}$, ${}^{14}\text{N}^*$, ${}^{14}\text{O}$), the measured mass excesses of the IAS members give: - $a = 2863$ keV (the charge-independent mass excess, close to $\Delta({}^{14}\text{N})$) - $b = -3880$ keV (the Coulomb displacement — negative because adding protons costs Coulomb energy, and our $T_3$ convention has $T_3 > 0$ for proton excess) - $c = 235$ keV (the quadratic correction)

2.6.3 Mirror Nuclei: A Precision Test

Mirror nuclei are pairs of nuclei obtained by interchanging $Z$ and $N$: if one nucleus has ($Z$, $N$), its mirror has ($N$, $Z$). Since $A = Z + N$ is the same, mirror nuclei are isobars. If isospin symmetry were exact, mirror nuclei would have identical energy level spectra (after correcting for the Coulomb energy).

The comparison of mirror pairs is one of the cleanest tests of isospin symmetry:

The $A = 3$ system: ${}^{3}\text{H}$ and ${}^{3}\text{He}$

These are the simplest mirror pair. Both have $J^\pi = 1/2^+$ and $T = 1/2$. Their binding energies are:

$$B({}^{3}\text{H}) = 8.482 \text{ MeV}, \quad B({}^{3}\text{He}) = 7.718 \text{ MeV}$$

The difference, $\Delta B = 0.764$ MeV, is almost entirely accounted for by the Coulomb energy difference between a system with one proton (${}^{3}\text{H}$) and two protons (${}^{3}\text{He}$). A simple estimate using a uniform sphere of radius $R = r_0 A^{1/3} = 1.2 \times 3^{1/3} = 1.73$ fm gives:

$$\Delta E_C = \frac{3}{5}\frac{e^2}{4\pi\epsilon_0}\frac{Z'(Z'-1) - Z(Z-1)}{R} = \frac{3}{5}\frac{1.44 \text{ MeV}\cdot\text{fm}}{1.73 \text{ fm}}(2 \times 1 - 1 \times 0) = 0.998 \text{ MeV}$$

The overestimate (0.998 vs. 0.764 MeV) reflects the fact that the uniform-sphere model overestimates the Coulomb energy; using a more realistic charge distribution reduces the Coulomb estimate and brings it into better agreement. The remaining discrepancy after a careful Coulomb correction — about 60–80 keV — is attributed to the charge-symmetry-breaking component of the nuclear force and the $n$-$p$ mass difference.

The $A = 11$ system: ${}^{11}\text{B}$ and ${}^{11}\text{C}$

Both have ground-state $J^\pi = 3/2^-$ and $T = 1/2$. Their excitation spectra are nearly identical up to a uniform Coulomb shift — the energies of the first several excited states agree to better than 100 keV after the Coulomb correction. This is a stringent test of isospin symmetry: the nuclear force produces (almost) the same level scheme in both nuclei, and the deviations are consistent with the known sources of isospin breaking.

Example 2.10: For the mirror pair ${}^{41}\text{Sc}$ ($Z = 21$, $N = 20$) and ${}^{41}\text{Ca}$ ($Z = 20$, $N = 21$), both have $J^\pi = 7/2^-$ and $T = 1/2$. These are particularly clean mirror nuclei because the odd nucleon (a proton in ${}^{41}\text{Sc}$, a neutron in ${}^{41}\text{Ca}$) sits outside a doubly magic ${}^{40}\text{Ca}$ core.

The Coulomb energy difference in the uniform-sphere approximation is:

$$\Delta E_C \approx \frac{3}{5}\frac{e^2}{4\pi\epsilon_0 R}[Z'(Z'-1) - Z(Z-1)] = \frac{3}{5}\frac{1.44}{1.2 \times 41^{1/3}}(21 \times 20 - 20 \times 19) = \frac{3}{5}\frac{1.44}{4.13}(420 - 380)$$

$$= 0.6 \times 0.349 \times 40 = 8.37 \text{ MeV}$$

The measured mass difference $M({}^{41}\text{Sc}) - M({}^{41}\text{Ca}) = 7.278$ MeV. The ~1 MeV discrepancy comes from two sources: (1) the oversimplified uniform-sphere Coulomb model overestimates the Coulomb energy because the charge distribution is diffuse, not uniform; (2) the exchange Coulomb energy (the quantum mechanical correction from antisymmetrization) is not included. A proper Hartree-Fock calculation with realistic charge distributions gives a Coulomb energy difference of $\sim 6.76$ MeV, which is closer but still underestimates the measured difference by about 0.5 MeV. This residual discrepancy is the Nolen-Schiffer anomaly — a persistent puzzle that has been partially attributed to charge-symmetry-breaking components of the nuclear force (see Exercises, Problem 30).

2.6.4 Isospin Selection Rules

Isospin conservation by the strong interaction leads to powerful selection rules:

• Nuclear reactions: The strong-interaction Hamiltonian commutes with $\hat{T}^2$ and $\hat{T}_3$, so nuclear reactions conserve $T$ (approximately) and $T_3$ (exactly, since $T_3$ is related to charge conservation). For example, the reaction ${}^{14}\text{N}(\alpha, p){}^{17}\text{O}$ is allowed only for specific isospin channels.

• Gamma-ray transitions: The electromagnetic interaction has $\Delta T = 0$ and $\Delta T = 1$ components. The $E1$ operator is predominantly isovector ($\Delta T = 0, 1$), leading to the isospin selection rule that $E1$ transitions between states of the same isospin ($\Delta T = 0$) are hindered in self-conjugate ($N = Z$) nuclei.

• Beta decay: The weak interaction changes $T_3$ by $\pm 1$ (since it changes a neutron to a proton or vice versa). Fermi beta decay has $\Delta T = 0$ (the isospin analog transition), while Gamow-Teller decay can have $\Delta T = 0$ or $\pm 1$.

Example 2.11: The reaction $d + d \to {}^{4}\text{He} + \pi^0$ is forbidden by isospin conservation, even though it conserves charge, baryon number, and energy-momentum. The reasoning: the deuteron has $T = 0$, so the initial state has $T_\text{initial} = 0$. The alpha particle has $T = 0$ and $\pi^0$ is the $T_3 = 0$ member of the $T = 1$ pion triplet. Thus $T_\text{final} = 1 \neq T_\text{initial} = 0$. This reaction has never been observed despite extensive searches, providing strong experimental evidence for isospin conservation in the strong interaction.

By contrast, $p + p \to d + \pi^+$ is allowed: the initial $pp$ state has $T = 1$, $T_3 = 1$; the final state has $T(d) = 0$, $T(\pi^+) = 1$ (with $T_3 = 1$), and the total isospin $T = 1$, $T_3 = 1$ is conserved.

🔗 Forward Reference: Isospin will return as a central organizing principle in Chapter 3 (charge independence of the nuclear force), Chapter 14 (beta decay selection rules), and Chapter 31 (the Standard Model, where isospin connects to the fundamental $u$-$d$ quark symmetry).

📊 The Big Picture: Isospin is the most powerful approximate symmetry in nuclear physics after rotational invariance. It organizes the nuclear level scheme into multiplets, provides selection rules for reactions and decays, and connects directly to the quark-level symmetry $SU(2)_\text{flavor}$ between $u$ and $d$ quarks. The precision with which isospin symmetry holds — and the pattern of its breaking — constrains our understanding of the nuclear force at the most fundamental level.

2.7 Spaced Review: Connections to Chapter 1

Before proceeding, let us connect the properties surveyed in this chapter to concepts from Chapter 1. These connections reinforce the key ideas and help build the web of understanding that characterizes expert knowledge of nuclear physics.

1. Nuclear density saturation (Ch 1) is the physical basis for the radius formula $R = r_0 A^{1/3}$ (Section 2.1). The constant density explains why $r_0$ is universal. In Chapter 3, we will see that saturation arises from the balance between the attractive nuclear force at medium range ($\sim 1$ fm) and the repulsive core at short range ($< 0.5$ fm).

2. The binding energy curve $B/A$ vs. $A$ (Ch 1) is now understood more precisely through mass excesses and separation energies (Section 2.2). The local structure in $B/A$ — the peaks at $A = 4, 12, 16, 28, 56$ — will be explained by shell structure, which is signaled by the separation energy systematics. The smooth trend is captured by the semi-empirical mass formula (Chapter 4); the deviations are the domain of the shell model (Chapter 6).

3. The chart of nuclides (Ch 1) organizes all the properties surveyed here. Nuclear radii, spins, moments, and masses all show systematic trends when plotted on the $N$-$Z$ plane. The magic numbers appear as horizontal and vertical lines of enhanced stability; the deformed regions appear as islands of large quadrupole moments between magic numbers.

4. Stable vs. unstable nuclei (Ch 1): The drip lines (Section 2.2.4) define the boundaries of nuclear existence. The valley of stability is the region where both $S_n > 0$ and $S_p > 0$. The total number of bound nuclides is estimated at $\sim 7,000$, of which only about 3,400 have been experimentally observed as of 2024.

5. ${}^{208}\text{Pb}$ as a benchmark (introduced Ch 1): In this chapter, ${}^{208}\text{Pb}$ appears as the paradigmatic spherical, doubly magic nucleus. Its charge radius (5.501 fm) benchmarks our radius formula. Its ground-state spin ($0^+$) is the universal even-even prediction. Its quadrupole moment is zero (by angular momentum selection rules), and its neighbors' moments are consistent with single-particle values plus core polarization. It will return in virtually every chapter of this book.

2.8 Progressive Project: Nuclear Data Analysis Toolkit — Checkpoint 2

nuclear_data_access.py

In this chapter's contribution to the Nuclear Data Analysis Toolkit, we build nuclear_data_access.py — a module that:

1. Stores a curated subset of AME2020 mass data for ~50 key nuclei spanning the chart of nuclides.
2. Plots nuclear radii vs. $A^{1/3}$, demonstrating the linearity that encodes density saturation.
3. Plots two-neutron separation energies $S_{2n}$ for tin isotopes, revealing the magic number signature at $N = 82$.

The code uses only numpy and matplotlib — no external nuclear data libraries required. Run it with:

python nuclear_data_access.py

The script produces two publication-quality figures:

• Figure 1: Charge radii vs. $A^{1/3}$ with a linear fit extracting $r_0$.
• Figure 2: $S_{2n}$ vs. $N$ for Sn isotopes, with the shell closure marked.

See the code/ directory for the full implementation and code/project-checkpoint.md for the project status update.

✅ Self-Check: After running the code, verify that: (1) the fitted $r_0$ falls between 0.94 and 0.97 fm for the rms charge radius slope (equivalently, $r_0 \approx 1.21$ fm for the half-density radius), and (2) the $S_{2n}$ plot shows a clear drop at $N = 82$.

2.9 Summary and What's Next

This chapter has established the empirical landscape that nuclear theory must explain:

• Nuclear sizes scale as $R = r_0 A^{1/3}$ with a diffuse surface, reflecting saturation of the nuclear density.
• Nuclear masses are known with sub-keV precision from Penning traps, and the separation energy systematics reveal shell closures.
• Nuclear spins follow from the coupling of nucleon angular momenta, with the pairing interaction driving all even-even ground states to $0^+$.
• Magnetic moments fall between but not on the Schmidt lines, signaling that the single-particle picture needs corrections.
• Quadrupole moments range from near-zero (spherical, near magic numbers) to several barns (strongly deformed, midshell), providing the most direct evidence for nuclear deformation.
• Isospin is an approximate symmetry of the nuclear force, tested precisely through mirror nuclei.

In Chapter 3, we turn to the force responsible for all of this: the nuclear force that holds protons and neutrons together. We will begin with the simplest bound nucleus — the deuteron — and build up to the modern understanding of nuclear interactions based on meson exchange and chiral effective field theory. The deuteron, with its small but nonzero quadrupole moment (Section 2.5), will provide our first evidence for the tensor component of the nuclear force.

The chain of reasoning in Part I now becomes clear: Chapter 1 established the landscape (the chart of nuclides, the binding energy curve). This chapter has catalogued the properties that any theory must explain. Chapter 3 will provide the force, Chapter 4 the simplest model (the semi-empirical mass formula), and Chapter 5 the mathematical tools. With all five chapters complete, we will be prepared for the detailed nuclear structure physics of Part II.

🔗 Chapter 2 at a Glance: Nuclear radii ($R = r_0 A^{1/3}$) encode density saturation. Masses (Penning traps, AME2020) reveal shell closures through separation energies. Spins ($0^+$ for all even-even nuclei) reflect pairing. Magnetic moments (Schmidt values ± corrections) probe single-particle structure. Quadrupole moments (small near magic numbers, large midshell) signal deformation. Isospin (tested by mirror nuclei) is an approximate symmetry broken primarily by the Coulomb force.

Chapter 2 of "Nuclear Physics: From Nuclear Forces to Neutron Stars." Suggested next: Chapter 3 — The Nuclear Force: What Holds the Nucleus Together.