Chapter 1 — The Discovery of the Nucleus: From Rutherford to the Chart of Nuclides

"It was quite the most incredible event that ever happened to me in my life. It was almost as incredible as if you had fired a 15-inch shell at a piece of tissue paper and it came back and hit you." — Ernest Rutherford, on the back-scattering of alpha particles (1911)

Chapter Overview

In 1909, two young physicists — Hans Geiger and Ernest Marsden — fired alpha particles at a thin gold foil and observed something that should not have happened. Most alpha particles sailed through as expected, but a few bounced nearly straight back. That observation, and Rutherford's brilliant analysis of it, revealed that atoms contain a tiny, massive core: the nucleus.

This chapter traces the path from that discovery to the modern organizing framework of nuclear physics: the chart of nuclides, a map of every known nucleus classified by its proton and neutron numbers. Along the way, we will:

• Derive the Rutherford scattering cross section from first principles and compare it to the original Geiger-Marsden data.
• Learn the nuclear notation ($A$, $Z$, $N$) that labels every nuclide and understand the classification of isotopes, isotones, and isobars.
• Measure nuclear sizes using electron scattering and discover the remarkable fact that nuclear density is nearly constant — a property called density saturation.
• Compute binding energies from measured masses and construct the most important plot in nuclear physics: the binding energy per nucleon curve — the single graph that explains why both fission and fusion release energy.

Every concept introduced here recurs throughout this book. The chart of nuclides is the map; the binding energy curve is the legend; and the nucleus itself — a tiny, dense, quantum-mechanical system of strongly interacting protons and neutrons — is the territory we will spend forty chapters exploring.

🏃 Fast Track: If you have seen Rutherford scattering before, skim Sections 1.1–1.2 and begin at Section 1.3 (nuclear notation). The essential new material is in Sections 1.5 (nuclear sizes) and 1.6 (binding energy).

🔬 Deep Dive: The full Rutherford derivation (Section 1.2) rewards careful study — the relationship between impact parameter, scattering angle, and cross section is a pattern you will encounter repeatedly in scattering theory (Chapter 10).

1.1 The Road to the Nucleus

1.1.1 The State of the Atom in 1909

By the first decade of the twentieth century, physicists knew that atoms contained negatively charged electrons (Thomson, 1897) and that the atom as a whole was electrically neutral. The dominant model was J.J. Thomson's "plum pudding" — electrons embedded in a diffuse sphere of positive charge, like raisins in a pudding. This model made a clear prediction: alpha particles (helium nuclei, charge $+2e$, mass $\approx 4\,\text{u}$) fired at a thin foil should experience only small deflections as they passed through the smeared-out positive charge.

📜 Historical Context: Alpha particles had been identified by Rutherford himself in 1899 as a form of radiation emitted by uranium. By 1908, he and Thomas Royds had conclusively shown that alpha particles are helium nuclei by trapping them and observing helium spectral lines.

1.1.2 The Geiger-Marsden Experiment

In Rutherford's laboratory at the University of Manchester, Geiger and Marsden set up a deceptively simple experiment. A collimated beam of alpha particles from a radioactive source struck a thin gold foil (thickness $\sim 0.4\,\mu\text{m}$, roughly 1000 atomic layers). A zinc sulfide screen — which produced a visible scintillation when struck by a single alpha particle — was placed at various angles $\theta$ around the foil.

The overwhelming majority of alpha particles passed straight through or scattered at small angles. But a small fraction — about 1 in 8000 — scattered at angles greater than 90°. Some bounced nearly straight back ($\theta \approx 180°$).

The plum pudding model could not explain this. A diffuse charge distribution produces only small-angle scattering; the probability of a large deflection drops exponentially with angle. Rutherford recognized that only a concentrated point-like charge could produce such large-angle scattering through the Coulomb force.

💡 Intuition: Imagine rolling a marble across a table scattered with cotton balls versus one that has a few steel ball bearings hidden under a cloth. The cotton balls barely deflect the marble; the ball bearings send it careening off at large angles. The few large-angle scatters reveal the concentrated masses even though most of the table is empty.

1.2 Rutherford Scattering: The Derivation

We now derive the Rutherford scattering formula — the differential cross section for Coulomb scattering of a charged particle by a point-like nucleus. This derivation is a masterclass in classical scattering theory.

1.2.1 Setup: The Scattering Geometry

Consider an alpha particle (charge $z_1 e = 2e$, mass $m_\alpha$) approaching a gold nucleus (charge $z_2 e = 79e$, mass $M_\text{Au}$) with kinetic energy $T$ in the laboratory frame. Since $M_\text{Au} \gg m_\alpha$, we treat the gold nucleus as stationary (we will correct for recoil shortly).

The alpha particle approaches along a straight line offset from the nucleus by the impact parameter $b$ — the perpendicular distance between the initial trajectory and a parallel line through the nucleus. The Coulomb repulsion between the two charges deflects the alpha particle through a scattering angle $\theta$.

The Coulomb potential energy at separation $r$ is:

$$U(r) = \frac{k \, z_1 z_2 e^2}{r} = \frac{z_1 z_2 e^2}{4\pi\epsilon_0 \, r}$$

1.2.2 Relating Impact Parameter to Scattering Angle

This is a central-force problem. The trajectory is a hyperbola (repulsive Coulomb), and the relationship between $b$ and $\theta$ can be derived from conservation of energy and angular momentum.

Conservation of energy: At infinite separation, all energy is kinetic: $E = T = \frac{1}{2}\mu v_0^2$ where $\mu$ is the reduced mass (approximately $m_\alpha$ for a heavy target) and $v_0$ is the initial speed. At the distance of closest approach $r_\text{min}$:

$$T = \frac{1}{2}\mu v_\text{min}^2 + \frac{k z_1 z_2 e^2}{r_\text{min}}$$

Conservation of angular momentum: $L = \mu v_0 b = \mu v_\text{min} r_\text{min}$.

The orbit equation: For the Coulomb potential, the classical orbit of the projectile is:

$$\frac{1}{r} = \frac{1}{a}\left(-1 + \frac{1}{\varepsilon}\cos(\phi - \phi_0)\right)$$

where $a = \frac{k z_1 z_2 e^2}{2T}$ is half the distance of closest approach for a head-on collision, and $\varepsilon = \sqrt{1 + (b/a)^2}$ is the eccentricity of the hyperbola.

The key geometric result — derived by analyzing the asymptotic angles of the hyperbola — is:

$$\boxed{b = a \cot\frac{\theta}{2}}$$

where

$$a \equiv \frac{z_1 z_2 e^2}{4\pi\epsilon_0 \cdot 2T} = \frac{k z_1 z_2 e^2}{2T}$$

is the half-distance of closest approach for a head-on collision ($b = 0$), also sometimes written $d/2$ where $d$ is the full distance of closest approach.

🔄 Check Your Understanding: For 5.5 MeV alpha particles on gold ($z_1 = 2$, $z_2 = 79$), compute $a$. You should get $a \approx 20.7\,\text{fm}$. Verify: $a = (1.44\,\text{MeV}\cdot\text{fm})(2)(79)/(2 \times 5.5\,\text{MeV}) = 20.7\,\text{fm}$, using $ke^2 = 1.44\,\text{MeV}\cdot\text{fm}$.

Derivation of $b = a\cot(\theta/2)$:

We derive this from conservation laws. The angular momentum is $L = \mu v_0 b$. The total energy is $E = \frac{1}{2}\mu v_0^2 = T$. We define the Coulomb parameter $a = kz_1z_2e^2/(2T)$.

For a repulsive Coulomb orbit, the particle travels along a hyperbola. At asymptotically large distances, the particle approaches along one asymptote and recedes along the other. The angle between the two asymptotes (measured from the focus, i.e., the scattering center) is $\pi - \theta$, where $\theta$ is the scattering angle. From the geometry of the hyperbola:

$$\cos\left(\frac{\pi - \theta}{2}\right) = \frac{1}{\varepsilon}$$

which gives $\sin(\theta/2) = 1/\varepsilon$. Since $\varepsilon = \sqrt{1 + (b/a)^2}$:

$$\sin^2\frac{\theta}{2} = \frac{1}{1 + (b/a)^2}$$

Solving for $b$:

$$\frac{b^2}{a^2} = \frac{1}{\sin^2(\theta/2)} - 1 = \frac{\cos^2(\theta/2)}{\sin^2(\theta/2)} = \cot^2\frac{\theta}{2}$$

Therefore $b = a\cot(\theta/2)$, as claimed. $\square$

1.2.3 From Impact Parameter to Cross Section

The differential cross section $d\sigma/d\Omega$ is defined such that the number of particles scattered into solid angle $d\Omega$ at angle $\theta$ is proportional to $d\sigma/d\Omega$. Physically, it connects the microscopic scattering dynamics to the macroscopic counting rate.

For azimuthally symmetric scattering (no $\phi$ dependence), particles hitting an annular ring of impact parameters between $b$ and $b + db$ scatter into angles between $\theta$ and $\theta + d\theta$. The area of the annular ring is:

$$d\sigma = 2\pi b \, |db|$$

The solid angle element is:

$$d\Omega = 2\pi \sin\theta \, |d\theta|$$

Therefore:

$$\frac{d\sigma}{d\Omega} = \frac{b}{\sin\theta} \left|\frac{db}{d\theta}\right|$$

From $b = a\cot(\theta/2)$:

$$\frac{db}{d\theta} = -\frac{a}{2\sin^2(\theta/2)}$$

Substituting:

$$\frac{d\sigma}{d\Omega} = \frac{a\cot(\theta/2)}{\sin\theta} \cdot \frac{a}{2\sin^2(\theta/2)}$$

Using $\sin\theta = 2\sin(\theta/2)\cos(\theta/2)$:

$$\frac{d\sigma}{d\Omega} = \frac{a^2 \cos(\theta/2)/\sin(\theta/2)}{2\sin(\theta/2)\cos(\theta/2)} \cdot \frac{1}{2\sin^2(\theta/2)}$$

$$\boxed{\frac{d\sigma}{d\Omega}\bigg|_{\text{Ruth}} = \left(\frac{a}{2}\right)^2 \frac{1}{\sin^4(\theta/2)} = \left(\frac{z_1 z_2 e^2}{4 \cdot 4\pi\epsilon_0 T}\right)^2 \frac{1}{\sin^4(\theta/2)}}$$

This is the Rutherford scattering formula. Often written compactly as:

$$\frac{d\sigma}{d\Omega}\bigg|_{\text{Ruth}} = \left(\frac{a}{2\sin^2(\theta/2)}\right)^2$$

1.2.4 Features of the Rutherford Cross Section

Several features are worth noting:

1. $\sin^{-4}(\theta/2)$ divergence: The cross section diverges as $\theta \to 0$ (forward scattering). This is physically reasonable — distant encounters produce small deflections, and the Coulomb force is long-ranged. In practice, atomic electron screening cuts off the divergence at very small angles.

2. $T^{-2}$ dependence: Higher-energy projectiles are harder to deflect, so the cross section decreases as the square of the energy.

3. $(z_1 z_2)^2$ dependence: Heavier nuclei (larger $Z$) scatter more strongly. This is why Rutherford used gold — a high-$Z$ target maximizes the signal.

4. No dependence on nuclear size: The formula assumes a point charge. Deviations from the Rutherford formula at large angles (small impact parameters) reveal that the nucleus has finite extent — this is how we measure nuclear radii (Section 1.5).

📊 Real-World Application: Rutherford Back-Scattering (RBS) is used today in materials science to determine the composition and thickness of thin films. A beam of MeV helium ions is directed at a sample, and the energy spectrum of back-scattered ions reveals which elements are present and at what depth. The technique is directly descended from Geiger and Marsden's original experiment.

1.2.5 Numerical Example

For 5.5 MeV alpha particles on gold ($z_1 = 2$, $z_2 = 79$):

$$a = \frac{(1.44\,\text{MeV}\cdot\text{fm})(2)(79)}{2(5.5\,\text{MeV})} = 20.7\,\text{fm}$$

At $\theta = 60°$:

$$\frac{d\sigma}{d\Omega} = \left(\frac{20.7}{2}\right)^2 \cdot \frac{1}{\sin^4(30°)} = (10.35)^2 \cdot \frac{1}{(0.5)^4} = 107.1 \times 16 = 1714\,\text{fm}^2/\text{sr}$$

Converting: $1714\,\text{fm}^2 = 17.14\,\text{b}$ (barns), since $1\,\text{b} = 100\,\text{fm}^2$.

At $\theta = 150°$:

$$\frac{d\sigma}{d\Omega} = (10.35)^2 \cdot \frac{1}{\sin^4(75°)} = 107.1 \cdot \frac{1}{(0.9659)^4} = 107.1 \cdot \frac{1}{0.8706} = 123\,\text{fm}^2/\text{sr} = 1.23\,\text{b/sr}$$

The ratio of cross sections at $60°$ and $150°$ is $1714/123 \approx 14$, reflecting the strong forward peaking.

⚠️ Common Pitfall: Students often confuse the half-distance of closest approach $a = kz_1z_2e^2/(2T)$ with the full distance of closest approach for a head-on collision $d_0 = 2a = kz_1z_2e^2/T$. Be careful about factors of 2. In this book we use $a = d_0/2$.

1.3 Chadwick's Discovery of the Neutron (1932)

1.3.1 The Missing Piece

By 1920, Rutherford had established that the nucleus contains protons ($Z$ of them, matching the atomic number). But nuclei were too heavy: helium had $Z = 2$ but mass number $A = 4$. Where did the extra mass come from?

Rutherford hypothesized in his 1920 Bakerian Lecture that there might exist an electrically neutral particle of approximately the proton's mass — a "neutron." For twelve years, the neutron eluded detection. Being electrically neutral, it does not ionize matter and leaves no track in cloud chambers or photographic emulsions.

1.3.2 Beryllium Radiation

In 1930, Walter Bothe and Herbert Becker observed that when alpha particles from polonium struck beryllium, a highly penetrating radiation was emitted. They assumed it was gamma radiation (high-energy photons). In 1932, Irène and Frédéric Joliot-Curie showed that this "beryllium radiation" could knock protons out of paraffin wax with considerable energy.

1.3.3 Chadwick's Analysis

James Chadwick, working at the Cavendish Laboratory in Cambridge, realized the critical inconsistency. If the beryllium radiation were gamma rays, conservation of energy and momentum required the photon energy to be implausibly large — about 55 MeV. But if the radiation consisted of neutral particles of approximately the proton's mass, the kinematics worked out perfectly.

The reaction is:

$${}^{9}\text{Be} + {}^{4}\text{He} \to {}^{12}\text{C} + n$$

Chadwick measured the recoil energies of both protons and nitrogen nuclei struck by the mystery radiation. From the kinematics of elastic collisions, he deduced:

• Recoil proton energy: $T_p \leq 5.7\,\text{MeV}$
• Recoil nitrogen energy: $T_N \leq 1.4\,\text{MeV}$

For a neutral particle of mass $m_n$ colliding elastically with a stationary target of mass $M$, the maximum recoil energy is:

$$T_\text{max} = \frac{4 m_n M}{(m_n + M)^2} T_n$$

where $T_n$ is the neutron kinetic energy. Taking the ratio for proton and nitrogen recoils eliminates $T_n$:

$$\frac{T_p}{T_N} = \frac{4m_n m_p/(m_n+m_p)^2}{4m_n m_N/(m_n+m_N)^2} = \frac{m_p(m_n+m_N)^2}{m_N(m_n+m_p)^2}$$

Inserting $T_p/T_N = 5.7/1.4 \approx 4.07$, $m_p = 1\,\text{u}$, $m_N = 14\,\text{u}$, and solving for $m_n$ gives $m_n \approx 1.15\,\text{u}$, close to the proton mass. Chadwick refined this to $m_n = 1.0 \pm 0.1\,\text{u}$, conclusively establishing the existence of the neutron.

📜 Historical Context: Chadwick's paper "The Existence of a Neutron" appeared in the Proceedings of the Royal Society A in 1932. It is a model of concise scientific writing — only 4 pages. He received the Nobel Prize in Physics in 1935. The modern value of the neutron mass is $m_n = 939.565\,\text{MeV}/c^2 = 1.008665\,\text{u}$.

🔄 Check Your Understanding: Verify Chadwick's result. With $T_p = 5.7\,\text{MeV}$ and $T_N = 1.4\,\text{MeV}$, confirm that $m_n \approx 1.15\,\text{u}$. (Hint: set up the ratio equation and solve the resulting quadratic in $m_n$.)

1.4 Nuclear Notation and the Classification of Nuclides

1.4.1 Standard Notation

A nuclide is a specific nuclear species characterized by its atomic number $Z$ (number of protons), neutron number $N$ (number of neutrons), and mass number $A = Z + N$ (total number of nucleons).

The standard notation is:

$${}^{A}_{Z}\text{X}_N$$

where X is the chemical symbol. Since $Z$ is uniquely determined by the chemical symbol and $N = A - Z$, the most common shorthand is simply ${}^{A}\text{X}$. For example:

1.4.2 Isotopes, Isotones, and Isobars

Nuclides are classified by which quantum number they share:

• Isotopes: Same $Z$ (same element), different $N$. Example: ${}^{12}\text{C}$, ${}^{13}\text{C}$, ${}^{14}\text{C}$ all have $Z = 6$.
• Isotones: Same $N$, different $Z$. Example: ${}^{13}\text{C}$ ($N = 7$) and ${}^{14}\text{N}$ ($N = 7$).
• Isobars: Same $A$, different $Z$. Example: ${}^{14}\text{C}$ ($Z = 6$) and ${}^{14}\text{N}$ ($Z = 7$).
• Isomers: Same $Z$ and $A$ but different energy states. Example: ${}^{99m}\text{Tc}$ (metastable excited state) and ${}^{99}\text{Tc}$ (ground state).

1.4.3 How Many Nuclides Exist?

As of the most recent Atomic Mass Evaluation (AME2020), approximately 3300 nuclides have been experimentally observed. Of these, only about 252 are stable (they do not decay on any experimentally accessible timescale) or 286 if we include those with half-lives exceeding the age of the universe. The rest are radioactive, decaying through alpha emission, beta decay, gamma emission, or more exotic modes (proton emission, neutron emission, fission).

The ratio is striking: for every stable nuclide, roughly 12 unstable ones are known. Nuclear instability is the rule; stability is the exception. Understanding why certain combinations of $Z$ and $N$ are stable is one of the central goals of nuclear physics.

🔗 Connection: The question of nuclear stability connects to astrophysics: the elements heavier than iron in your body were synthesized in supernovae and neutron star mergers. In 2017, gravitational waves from colliding neutron stars (event GW170817) were detected by LIGO and Virgo, followed by electromagnetic observations that confirmed heavy-element nucleosynthesis in the merger ejecta. We will tell this story in full in Chapter 23.

1.5 The Chart of Nuclides

1.5.1 A Map of All Known Nuclei

The chart of nuclides (sometimes called the Segrè chart, after Emilio Segrè) is the organizing framework for nuclear physics. It plots each known nuclide as a point (or colored square) on a grid with neutron number $N$ on the horizontal axis and proton number $Z$ on the vertical axis.

Each square is typically color-coded by decay mode: - Black: Stable - Blue: $\beta^-$ decay (neutron-rich side) - Red/Pink: $\beta^+$ decay or electron capture (proton-rich side) - Yellow: Alpha decay (heavy nuclei) - Green: Spontaneous fission - Orange: Proton emission

1.5.2 The Valley of Stability

The stable nuclides form a narrow band called the valley of stability. For light nuclei ($A \lesssim 40$), the valley follows $N \approx Z$. For heavier nuclei, the valley curves toward $N > Z$ — heavy nuclei require a neutron excess to remain stable. This is because the short-range nuclear (strong) force acts equally between all nucleon pairs ($pp$, $nn$, $pn$), but the long-range Coulomb repulsion acts only between proton pairs. Adding extra neutrons dilutes the Coulomb repulsion without weakening the nuclear attraction. (We will quantify this with the semi-empirical mass formula in Chapter 3.)

1.5.3 Drip Lines and the Limits of Nuclear Existence

Moving away from the valley of stability in either direction (adding protons or neutrons), nuclei become increasingly unstable. At the extreme, we reach the drip lines: the boundaries beyond which it is energetically favorable for the nucleus to spontaneously emit a proton or neutron.

• Proton drip line: The boundary beyond which the last proton is unbound. Experimentally well-established up to about $Z \sim 90$.
• Neutron drip line: The boundary beyond which the last neutron is unbound. Experimentally established only up to about $Z \sim 10$ (oxygen), because producing extremely neutron-rich nuclei requires rare-isotope beam facilities.

📊 Real-World Application: Modern rare-isotope facilities — FRIB (Facility for Rare Isotope Beams) at Michigan State University, RIBF at RIKEN in Japan, FAIR at GSI in Germany — are designed specifically to produce and study nuclei near and beyond the drip lines. Understanding these exotic nuclei tests our models of nuclear structure in extreme conditions.

1.5.4 Magic Numbers

Certain proton or neutron numbers produce nuclei of exceptional stability. These magic numbers are:

$$2, \quad 8, \quad 20, \quad 28, \quad 50, \quad 82, \quad 126$$

Doubly magic nuclei — magic in both $Z$ and $N$ — are particularly stable. The most prominent examples are:

${}^{208}\text{Pb}$ will serve as one of our anchor examples throughout this book. Its exceptional stability — it is the heaviest stable nucleus with a doubly magic configuration — makes it an ideal testing ground for nuclear models. We will return to it in Chapter 6 (shell model), Chapter 11 (collective excitations), and Chapter 15 (nuclear reactions).

🔄 Check Your Understanding: On a chart of nuclides, where would you find the magic numbers? They appear as horizontal lines ($Z$ magic) and vertical lines ($N$ magic). Nuclei along these lines tend to have more stable isotopes/isotones than average.

1.5.5 The N-Z Pattern: Why Heavy Nuclei Need More Neutrons

The valley of stability departs from $N = Z$ in a systematic way that reflects the competition between the nuclear force and the Coulomb force. Let us trace the pattern quantitatively:

The physical explanation is straightforward. The nuclear (strong) force acts between all nucleon pairs — $pp$, $nn$, and $np$ — with nearly equal strength (charge independence, Chapter 2). But the Coulomb repulsion acts only between proton pairs. As the proton number $Z$ increases, the total Coulomb energy grows roughly as $Z^2$ (more precisely, as $Z(Z-1)/2$ for $Z$ protons). To compensate this growing repulsion without sacrificing nuclear attraction, the nucleus adds extra neutrons. The neutrons contribute attractive nuclear force without adding Coulomb repulsion, effectively diluting the electrostatic cost.

This has a profound consequence: the heaviest elements cannot exist with equal numbers of protons and neutrons. A hypothetical nucleus with $Z = N = 92$ ($A = 184$) would be catastrophically unstable — the Coulomb repulsion of 92 protons in a volume of radius $\sim 7\,\text{fm}$ would tear it apart. Uranium exists only because it has 146 neutrons providing additional binding while contributing no Coulomb repulsion.

Pairing and the even-odd effect. The chart of nuclides also reveals that even-$Z$, even-$N$ nuclei (even-even) are more abundant and more stable than odd-$A$ or odd-odd nuclei. Of the approximately 252 stable nuclides: 166 are even-even, 53 are even-$Z$/odd-$N$, 29 are odd-$Z$/even-$N$, and only 4 are odd-odd (${}^{2}\text{H}$, ${}^{6}\text{Li}$, ${}^{10}\text{B}$, ${}^{14}\text{N}$). This pattern reflects the pairing interaction — the tendency of nucleons to form spin-zero pairs, lowering the total energy. Each pair of identical nucleons (two protons or two neutrons) in the same orbital with opposite spin projections gains a pairing energy of roughly $\delta \approx 12/\sqrt{A}\,\text{MeV}$. Even-even nuclei benefit from full pairing in both the proton and neutron systems; odd-odd nuclei have one unpaired proton and one unpaired neutron, making them less stable.

Isobars and beta stability. For a fixed mass number $A$, the most stable isobar is the one closest to the bottom of the valley — the optimal $Z$ for that $A$. Isobars with too many neutrons undergo $\beta^-$ decay ($n \to p + e^- + \bar{\nu}_e$), moving toward higher $Z$. Isobars with too many protons undergo $\beta^+$ decay ($p \to n + e^+ + \nu_e$) or electron capture, moving toward lower $Z$. For odd-$A$ nuclei, there is typically one stable isobar. For even-$A$ nuclei, the pairing interaction creates two separate parabolas (one for even-even, one for odd-odd), and there can be two or even three stable isobars.

1.6 Nuclear Sizes: Electron Scattering and Density Saturation

1.6.1 Why Electrons?

To measure the size of a nucleus, we need a probe that interacts with it primarily through a known force. Electrons are ideal: they are point-like (no internal structure) and interact with the nuclear charge distribution through the well-understood electromagnetic force. At the energies used ($100$–$500\,\text{MeV}$), the electron's de Broglie wavelength is:

$$\lambda = \frac{\hbar c}{pc} \approx \frac{197\,\text{MeV}\cdot\text{fm}}{200\,\text{MeV}} \approx 1\,\text{fm}$$

This is comparable to nuclear dimensions, making electrons excellent probes of nuclear structure.

1.6.2 The Nuclear Charge Distribution: Hofstadter's Experiments

Robert Hofstadter (Nobel Prize 1961) pioneered elastic electron-nucleus scattering at Stanford in the 1950s, using the Mark III linear electron accelerator to produce beams of 100–550 MeV electrons. His program fundamentally changed our understanding of nuclear structure by providing the first model-independent measurements of the charge distribution inside nuclei.

The method. An electron scattering from a point-like nucleus obeys the Mott cross section — the relativistic generalization of the Rutherford formula that accounts for the electron's spin:

$$\frac{d\sigma}{d\Omega}\bigg|_{\text{Mott}} = \frac{d\sigma}{d\Omega}\bigg|_{\text{Ruth}} \cdot \left(1 - \beta^2 \sin^2\frac{\theta}{2}\right)$$

where $\beta = v/c$ for the electron. For a nucleus with a finite charge distribution $\rho_{\text{ch}}(r)$, the measured cross section deviates from the Mott prediction by the square of the charge form factor $F(q)$:

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

where $q$ is the magnitude of the three-momentum transfer. For elastic scattering, $q = (4\pi/\lambda)\sin(\theta/2) = 2k\sin(\theta/2)$. The form factor is the Fourier transform of the charge distribution:

$$F(q) = \frac{4\pi}{Ze} \int_0^{\infty} \rho_{\text{ch}}(r) \, \frac{\sin(qr)}{qr} \, r^2 \, dr$$

By measuring $d\sigma/d\Omega$ over a wide range of angles (and therefore a wide range of $q$), one maps out $|F(q)|^2$ and can invert the Fourier transform to obtain $\rho_{\text{ch}}(r)$. In practice, one fits the data to a parameterized density distribution and adjusts the parameters.

Key results. Hofstadter's group measured scattering from targets ranging from ${}^{12}\text{C}$ to ${}^{209}\text{Bi}$. The form factors showed striking diffraction patterns — minima at specific momentum transfers, analogous to single-slit diffraction in optics. The position of the first diffraction minimum is inversely related to the nuclear radius, confirming $R \propto A^{1/3}$. For heavy nuclei, the diffraction minima are closely spaced and deeply pronounced, reflecting a sharp surface. For light nuclei, the pattern is smoother, indicating a more diffuse charge distribution.

For calcium-40, Hofstadter measured the first diffraction minimum at $q \approx 1.85\,\text{fm}^{-1}$. From the relation $qR \approx 4.49$ (the first zero of $j_1(x)/x$, the form factor for a uniform sphere), this gives $R \approx 2.43 \cdot A^{1/3} \approx 3.5\,\text{fm}$ — but the actual Fermi distribution fit yields more precise parameters accounting for the surface diffuseness.

Beyond Hofstadter: modern electron scattering. Facilities such as the Mainz Microtron (MAMI) in Germany and Jefferson Lab (JLab) in Virginia have pushed electron scattering to unprecedented precision. The PREX and CREX experiments at JLab used parity-violating electron scattering (measuring the interference between photon and $Z^0$ exchange) to determine the neutron skin thickness of ${}^{208}\text{Pb}$ and ${}^{48}\text{Ca}$ — the difference between the neutron and proton rms radii. The PREX-2 result (Adhikari et al., 2021, Physical Review Letters 126:172502) found a neutron skin of $\Delta r_{np} = 0.283 \pm 0.071\,\text{fm}$ for ${}^{208}\text{Pb}$, a result with profound implications for the equation of state of neutron-rich matter and the structure of neutron stars (Chapter 22).

The measured angular distributions — which deviate from the Rutherford/Mott formula at large angles — can be inverted to extract the nuclear charge density $\rho_{\text{ch}}(r)$.

The result, for nuclei with $A \gtrsim 20$, is well described by the Fermi (or Woods-Saxon) distribution:

$$\rho_{\text{ch}}(r) = \frac{\rho_0}{1 + \exp\left(\frac{r - R}{a}\right)}$$

where: - $\rho_0 \approx 0.07\,e/\text{fm}^3$ is the central charge density (approximately $0.16\,\text{nucleons}/\text{fm}^3$ for the total nucleon density) - $R$ is the half-density radius (where $\rho = \rho_0/2$) - $a \approx 0.54\,\text{fm}$ is the skin thickness parameter (the 90%-to-10% density falloff occurs over $\approx 4.4a \approx 2.4\,\text{fm}$)

1.6.3 The Nuclear Radius Formula

The half-density radius follows the empirical law:

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

with $r_0 = 1.21 \pm 0.02\,\text{fm}$ for the charge radius. This $A^{1/3}$ scaling has a profound implication.

Numerical examples with measured charge radii:

Note: The root-mean-square (rms) charge radius $\langle r^2 \rangle^{1/2}$ differs from the half-density radius $R$. For the Fermi distribution, $\langle r^2 \rangle^{1/2} \approx \sqrt{(3/5)} R \sqrt{1 + (7/3)(\pi a/R)^2}$. The measured rms charge radii are well tabulated (see Angeli and Marinova, 2013).

1.6.4 Nuclear Density Saturation

The volume of a nucleus scales as:

$$V = \frac{4}{3}\pi R^3 = \frac{4}{3}\pi r_0^3 A$$

Therefore the nuclear density is:

$$\rho_{\text{nuc}} = \frac{A}{V} = \frac{A}{(4/3)\pi r_0^3 A} = \frac{3}{4\pi r_0^3} \approx \frac{3}{4\pi (1.21)^3} \approx 0.135\,\text{nucleons/fm}^3$$

The $A$ cancels. Nuclear density is approximately the same for all nuclei — from carbon to uranium. This is nuclear density saturation, and it tells us something deep about the nuclear force: it is short-ranged and saturates. Each nucleon interacts primarily with its nearest neighbors, not with every other nucleon in the nucleus. (If the force were long-ranged like gravity, heavier nuclei would be denser.)

💡 Intuition: Nuclear matter behaves somewhat like an incompressible liquid. Adding more nucleons increases the volume proportionally, just as adding more water to a container increases the volume without changing the density. This is the physical basis for the liquid drop model of the nucleus (Chapter 3).

In more familiar units: $\rho_{\text{nuc}} \approx 2.3 \times 10^{17}\,\text{kg/m}^3$. A teaspoon of nuclear matter would weigh about 1 billion tonnes. Neutron stars achieve densities comparable to and exceeding nuclear density — we will explore this in Chapter 22.

1.6.5 Quantifying Nuclear Density: A Closer Look

The central nuclear density $\rho_0$ is one of the most fundamental constants in nuclear physics. Let us examine the numbers more carefully with several nuclei.

Density from the Fermi distribution. The normalization condition for the charge density is:

$$\int_0^{\infty} \rho_{\text{ch}}(r) \, 4\pi r^2 \, dr = Z e$$

For the total nucleon density (protons plus neutrons), the normalization gives:

$$\int_0^{\infty} \rho_{\text{nuc}}(r) \, 4\pi r^2 \, dr = A$$

For a uniform sphere of radius $R$, $\rho_0 = 3A/(4\pi R^3) = 3/(4\pi r_0^3)$. Using $r_0 = 1.21\,\text{fm}$:

$$\rho_0 = \frac{3}{4\pi (1.21)^3} = \frac{3}{7.42} = 0.135\,\text{nucleons/fm}^3$$

However, the actual Fermi distribution has a diffuse surface, which slightly modifies this value. When Hofstadter's data are fit to the Fermi form for a variety of nuclei, the extracted central densities are:

The central densities are all within $\pm 5\%$ of the value $\rho_0 \approx 0.16\,\text{nucleons/fm}^3$, confirming saturation across a wide range of nuclei.

Why does density saturate? The saturation of nuclear density is a direct consequence of the saturation property of the nuclear force — the fact that each nucleon interacts strongly with only a limited number of nearest neighbors (roughly 4–6), not with all $A - 1$ other nucleons. This is analogous to molecules in a liquid, where each molecule interacts with its nearest neighbors but not with distant ones. If the nuclear force were long-ranged (like gravity), every nucleon would attract every other, the nucleus would contract as $A$ increases, and the density would grow — which is not what is observed. The short range ($\sim 1$–$2\,\text{fm}$) and the repulsive core at very short distances ($< 0.5\,\text{fm}$) prevent the nucleons from collapsing into an ever-denser state.

This saturation property will be central to the liquid drop model (Chapter 3), where the nucleus is treated as a droplet of incompressible nuclear fluid with a well-defined surface tension. It also determines the nuclear equation of state, which governs the structure of neutron stars (Chapter 22) and the dynamics of core-collapse supernovae (Chapter 23).

1.7 Nuclear Binding Energy

1.7.1 The Mass Defect

The mass of a nucleus is less than the sum of the masses of its constituent nucleons. This "missing mass" is the mass defect $\Delta$:

$$\Delta = \left[Z m_p + N m_n\right] - M(A,Z)$$

where $M(A,Z)$ is the nuclear mass. By Einstein's mass-energy equivalence, the mass defect corresponds to the binding energy:

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

The binding energy is the energy required to completely disassemble the nucleus into $Z$ free protons and $N$ free neutrons. A larger binding energy means a more tightly bound (more stable) nucleus.

⚠️ Common Pitfall — The Threshold Concept of This Chapter: Binding energy is NOT "energy stored inside the nucleus." It is the opposite. A nucleus with a large binding energy has less mass than its constituent parts. You must add energy equal to $B$ to break it apart. Think of it as the depth of a potential well: the deeper the well, the more energy is required to climb out.

1.7.2 Working with Atomic Masses

In practice, masses are tabulated for neutral atoms, not bare nuclei. The atomic mass $M_{\text{atom}}$ includes the mass of $Z$ electrons plus the electronic binding energy (which is small compared to nuclear binding energies). Using atomic masses:

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

where $M({}^{1}\text{H}) = 1.007825\,\text{u}$ is the atomic mass of hydrogen. The electron masses cancel (hydrogen atom has one electron; the atom ${}^A_Z\text{X}$ has $Z$ electrons). The small electronic binding energy difference ($\sim\text{keV}$) is negligible compared to nuclear binding energies ($\sim\text{MeV}$).

The atomic mass unit is defined as $1\,\text{u} = 1/12$ the mass of ${}^{12}\text{C}$: $$1\,\text{u} = 931.494\,\text{MeV}/c^2$$

Fundamental masses:

1.7.3 Numerical Examples

Example 1: ${}^{4}\text{He}$ (alpha particle)

$$B({}^{4}\text{He}) = \left[2(1.007825) + 2(1.008665) - 4.002603\right] \times 931.494\,\text{MeV}$$ $$= \left[2.015650 + 2.017330 - 4.002603\right] \times 931.494$$ $$= 0.030377 \times 931.494 = 28.296\,\text{MeV}$$

So $B/A = 28.296/4 = 7.074\,\text{MeV/nucleon}$.

Example 2: ${}^{56}\text{Fe}$

$$B({}^{56}\text{Fe}) = \left[26(1.007825) + 30(1.008665) - 55.934936\right] \times 931.494$$ $$= \left[26.203450 + 30.259950 - 55.934936\right] \times 931.494$$ $$= 0.528464 \times 931.494 = 492.26\,\text{MeV}$$

So $B/A = 492.26/56 = 8.790\,\text{MeV/nucleon}$.

Example 3: ${}^{208}\text{Pb}$

$$B({}^{208}\text{Pb}) = \left[82(1.007825) + 126(1.008665) - 207.976652\right] \times 931.494$$ $$= \left[82.641650 + 127.091790 - 207.976652\right] \times 931.494$$ $$= 1.756788 \times 931.494 = 1636.43\,\text{MeV}$$

So $B/A = 1636.43/208 = 7.868\,\text{MeV/nucleon}$.

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

$$B({}^{238}\text{U}) = \left[92(1.007825) + 146(1.008665) - 238.050788\right] \times 931.494$$ $$= \left[92.719900 + 147.265090 - 238.050788\right] \times 931.494$$ $$= 1.934202 \times 931.494 = 1801.69\,\text{MeV}$$

So $B/A = 1801.69/238 = 7.570\,\text{MeV/nucleon}$.

🔄 Check Your Understanding: Verify the ${}^{12}\text{C}$ binding energy. With $M({}^{12}\text{C}) = 12.000000\,\text{u}$ exactly (by definition), you should find $B = 92.16\,\text{MeV}$, giving $B/A = 7.680\,\text{MeV/nucleon}$.

Example 5: ${}^{16}\text{O}$ (doubly magic)

$$B({}^{16}\text{O}) = \left[8(1.007825) + 8(1.008665) - 15.994915\right] \times 931.494$$ $$= \left[8.062600 + 8.069320 - 15.994915\right] \times 931.494$$ $$= 0.137005 \times 931.494 = 127.62\,\text{MeV}$$

So $B/A = 127.62/16 = 7.976\,\text{MeV/nucleon}$. Notice this is higher than the smooth trend for nuclei of this mass — the double magic numbers ($Z = 8$, $N = 8$) confer extra binding.

Example 6: ${}^{40}\text{Ca}$ (doubly magic)

$$B({}^{40}\text{Ca}) = \left[20(1.007825) + 20(1.008665) - 39.962591\right] \times 931.494$$ $$= \left[20.156500 + 20.173300 - 39.962591\right] \times 931.494$$ $$= 0.367209 \times 931.494 = 342.05\,\text{MeV}$$

So $B/A = 342.05/40 = 8.551\,\text{MeV/nucleon}$. Again, the double magic character ($Z = N = 20$) gives ${}^{40}\text{Ca}$ extra binding relative to its neighbors.

Example 7: ${}^{120}\text{Sn}$ (magic $Z = 50$)

$$B({}^{120}\text{Sn}) = \left[50(1.007825) + 70(1.008665) - 119.902199\right] \times 931.494$$ $$= \left[50.391250 + 70.606550 - 119.902199\right] \times 931.494$$ $$= 1.095601 \times 931.494 = 1020.55\,\text{MeV}$$

So $B/A = 1020.55/120 = 8.505\,\text{MeV/nucleon}$. Tin, with its magic proton number $Z = 50$, has the largest number of stable isotopes of any element (ten), a direct consequence of the extra stability conferred by the closed proton shell.

📊 Summary Table of Binding Energies:

Nuclide	$A$	$B$ (MeV)	$B/A$ (MeV/nucleon)	Notes
${}^{4}\text{He}$	4	28.30	7.074	Doubly magic
${}^{12}\text{C}$	12	92.16	7.680	Triple-$\alpha$ product
${}^{16}\text{O}$	16	127.62	7.976	Doubly magic
${}^{40}\text{Ca}$	40	342.05	8.551	Doubly magic
${}^{56}\text{Fe}$	56	492.26	8.790	Near peak
${}^{62}\text{Ni}$	62	545.26	8.795	Highest $B/A$
${}^{120}\text{Sn}$	120	1020.55	8.505	Magic $Z$
${}^{208}\text{Pb}$	208	1636.43	7.868	Doubly magic
${}^{238}\text{U}$	238	1801.69	7.570	Heaviest primordial

1.7.4 The Binding Energy Per Nucleon Curve

The most important plot in nuclear physics is $B/A$ versus $A$. It reveals the essential systematics of nuclear stability:

1. Rapid rise for light nuclei: $B/A$ increases steeply from deuterium ($1.112\,\text{MeV}$) through ${}^{4}\text{He}$ ($7.074\,\text{MeV}$) to about $A \sim 12$.

2. Broad maximum near $A \approx 56$–$62$: The peak is at ${}^{62}\text{Ni}$ ($8.795\,\text{MeV/nucleon}$) — the most tightly bound nucleus per nucleon. ${}^{58}\text{Fe}$ ($8.792$) and ${}^{56}\text{Fe}$ ($8.790$) are close behind. (It is a common misconception that ${}^{56}\text{Fe}$ is the most tightly bound; ${}^{62}\text{Ni}$ holds that distinction.)

3. Gradual decline for heavy nuclei: $B/A$ decreases slowly from the peak, reaching $7.570\,\text{MeV}$ at ${}^{238}\text{U}$.

4. Local peaks at magic numbers and $\alpha$-particle clusters: ${}^{4}\text{He}$, ${}^{12}\text{C}$, ${}^{16}\text{O}$, ${}^{40}\text{Ca}$ all sit above the smooth trend.

💡 Intuition — Why the Curve Has This Shape:

Rising part (light nuclei): As nucleons are added, each new nucleon interacts with more partners, increasing the total binding. But the nuclear force saturates — each nucleon has only a limited number of nearest neighbors.

Falling part (heavy nuclei): As the nucleus grows, the fraction of nucleons on the surface (with fewer neighbors) decreases slowly, but Coulomb repulsion between protons grows as $Z^2$ and eventually dominates, reducing $B/A$.

The interplay between the attractive nuclear force (saturating, short-ranged) and the repulsive Coulomb force (long-ranged, growing with $Z$) creates the peak and determines the boundary between fusion and fission.

1.7.5 Why Fission and Fusion Both Release Energy

The binding energy curve immediately explains the two processes by which nuclear energy can be released:

Fusion (light nuclei): Combining two light nuclei (to the left of the peak) into a heavier nucleus moves up the $B/A$ curve. The product is more tightly bound per nucleon than the reactants, so it has less mass. The mass difference is released as kinetic energy.

Example: ${}^{2}\text{H} + {}^{3}\text{H} \to {}^{4}\text{He} + n$

• Reactants: $B/A = 1.112$ (D) and $2.827$ (T)
• Product: $B/A = 7.074$ (He-4)
• Energy released: $\approx 17.6\,\text{MeV}$ (one of the most energetic fusion reactions)

Fission (heavy nuclei): Splitting a heavy nucleus (to the right of the peak) into two medium-mass fragments moves up the $B/A$ curve from the heavy side. Again, the products are more tightly bound, and energy is released.

Example: ${}^{235}\text{U} + n \to \text{fission fragments} + \text{neutrons}$

• ${}^{235}\text{U}$: $B/A \approx 7.59\,\text{MeV}$
• Typical fragments (e.g., ${}^{141}\text{Ba}$, ${}^{92}\text{Kr}$): $B/A \approx 8.3$–$8.5\,\text{MeV}$
• Energy released: $\approx 200\,\text{MeV}$ per fission event

⚠️ Common Pitfall: Students sometimes think "fission releases energy because heavy nuclei have lots of binding energy." This is backwards. Heavy nuclei do have large total $B$, but their $B/A$ is lower than the fragments. Energy is released because the products are more tightly bound per nucleon, not because binding energy is "released" from the parent.

🔗 Connection: The binding energy curve connects nuclear physics to both power generation and astrophysics. Fission powers nuclear reactors (${}^{235}\text{U}$) and was harnessed in the first nuclear weapons. Fusion powers every star in the universe — including the Sun, which fuses hydrogen into helium via the pp chain (Chapter 21). Every day in every hospital, nuclear physics also saves lives: PET (positron emission tomography) imaging uses the annihilation of positrons emitted by ${}^{18}\text{F}$-FDG, a radioactive fluorine-labeled glucose analog, to map metabolic activity in the body. The physics of positron emission starts with the concepts in this chapter and is developed in Chapter 5.

1.7.6 Separation Energies

A closely related quantity is the nucleon separation energy — the energy required to remove a single nucleon from the nucleus:

Neutron separation energy: $$S_n(A,Z) = B(A,Z) - B(A-1,Z) = \left[M(A-1,Z) + m_n - M(A,Z)\right]c^2$$

Proton separation energy: $$S_p(A,Z) = B(A,Z) - B(A-1,Z-1) = \left[M(A-1,Z-1) + M({}^{1}\text{H}) - M(A,Z)\right]c^2$$

Separation energies are the nuclear analogs of ionization energies in atomic physics. Discontinuities in $S_n$ and $S_p$ at magic numbers provide the most direct evidence for nuclear shell structure (Chapter 6).

Example — Neutron separation energies around $N = 126$ (Lead):

For ${}^{208}\text{Pb}$ ($N = 126$): - $S_n({}^{208}\text{Pb}) = B(208,82) - B(207,82) = 1636.43 - 1629.06 = 7.37\,\text{MeV}$ - $S_n({}^{209}\text{Pb}) = B(209,82) - B(208,82) = 1640.25 - 1636.43 = 3.82\,\text{MeV}$

The dramatic drop from $7.37$ to $3.82\,\text{MeV}$ when crossing $N = 126$ is direct evidence of the magic number. The 126th neutron is the last one to fill a shell; the 127th starts a new shell and is much less tightly bound.

Example — Proton separation energies around $Z = 50$ (Tin):

The proton magic number $Z = 50$ manifests in proton separation energies of tin and its neighbors:

• $S_p({}^{120}\text{Sn}) = B(120,50) - B(119,49) = 1020.55 - 1010.28 = 10.27\,\text{MeV}$ (adding the 50th proton)
• $S_p({}^{121}\text{Sb}) = B(121,51) - B(120,50) = 1026.33 - 1020.55 = 5.78\,\text{MeV}$ (adding the 51st proton)

The separation energy drops by nearly a factor of two when we cross the $Z = 50$ shell closure — the 51st proton occupies a higher-energy orbit outside the closed shell and is far less tightly bound.

Example — Neutron separation energies around $N = 82$ (Barium):

• $S_n({}^{138}\text{Ba}) = B(138,56) - B(137,56) = 1158.30 - 1149.69 = 8.61\,\text{MeV}$ ($N = 82$, closing the shell)
• $S_n({}^{139}\text{Ba}) = B(139,56) - B(138,56) = 1163.29 - 1158.30 = 4.99\,\text{MeV}$ ($N = 83$, opening the next shell)

Again, the signature is unmistakable: a drop of $\sim 3.6\,\text{MeV}$ at the shell closure.

💡 Insight — Separation Energies and the Atomic Analogy: Separation energies in nuclear physics play the same role as ionization energies in atomic physics. The ionization energy of a noble gas atom (closed electron shell) is much larger than that of the following alkali metal atom. Similarly, the nucleon separation energy of a magic nucleus is much larger than that of the nucleus one nucleon beyond the closed shell. This analogy was one of the original motivations for the nuclear shell model (Mayer and Jensen, 1949 — Chapter 6).

Systematic behavior. Plotting $S_n$ or $S_p$ across an isotopic or isotonic chain reveals not only the shell closures but also:

• Pairing staggering: $S_n$ zigzags between even-$N$ and odd-$N$ isotopes. Even-$N$ nuclei have higher $S_n$ because removing a neutron from a paired system costs extra energy (you must break a pair). The magnitude of the staggering, about $1$–$2\,\text{MeV}$ for medium-mass nuclei, is a direct measure of the pairing interaction strength.

• Smooth trends between shells: Between magic numbers, $S_n$ decreases gradually as neutrons fill a shell — each additional neutron in a partially filled shell is slightly less bound than the previous one, because the Fermi energy rises within the shell.

• Drip line prediction: The nuclear drip line is defined by $S_n = 0$ (neutron drip line) or $S_p = 0$ (proton drip line). Beyond the drip line, the last nucleon is unbound and the nucleus decays by nucleon emission on a timescale of $\sim 10^{-22}$ s.

Separation energies will recur throughout this book: in the shell model (Chapter 6) as the primary evidence for magic numbers, in decay energetics (Chapters 12–15) where the $Q$-value of a decay is a difference of separation energies, and in nuclear reactions (Chapters 17–20) where separation energies determine thresholds and final-state energies.

1.8 Project Checkpoint: Nuclear Data Analysis Toolkit

This chapter's contribution to the progressive project — the Nuclear Data Analysis Toolkit — consists of two Python scripts.

Script 1: binding_energy_curve.py

This script plots the binding energy per nucleon ($B/A$) versus mass number ($A$) for stable and long-lived nuclides. It uses a curated set of nuclear data values embedded directly in the code (from the 2020 Atomic Mass Evaluation, AME2020). The script:

• Plots $B/A$ vs. $A$ for approximately 50 representative nuclides spanning the full range of $A$
• Marks the peak near ${}^{62}\text{Ni}$ / ${}^{56}\text{Fe}$
• Shades the "fusion releases energy" region (light nuclei, to the left of the peak)
• Shades the "fission releases energy" region (heavy nuclei, to the right of the peak)
• Labels several notable nuclides: ${}^{2}\text{H}$, ${}^{4}\text{He}$, ${}^{12}\text{C}$, ${}^{56}\text{Fe}$, ${}^{62}\text{Ni}$, ${}^{208}\text{Pb}$, ${}^{238}\text{U}$

Script 2: nuclide_chart.py

This script creates a basic chart of nuclides visualization — a scatter plot of $N$ (horizontal axis) versus $Z$ (vertical axis), with points colored by stability or decay mode. It uses a representative dataset of known nuclides and illustrates:

• The valley of stability ($N \approx Z$ line for reference)
• The neutron-rich and proton-rich sides
• Magic number lines

See the code/ directory for the complete, runnable scripts and code/project-checkpoint.md for documentation.

💻 Computational Note: Both scripts require only numpy and matplotlib. Run them with: python binding_energy_curve.py and python nuclide_chart.py.

Chapter Summary

1. Rutherford scattering — the scattering of charged particles by a Coulomb potential — revealed the nuclear atom. The differential cross section $d\sigma/d\Omega = (a/2)^2 / \sin^4(\theta/2)$ is derived from the relationship $b = a\cot(\theta/2)$ and predicts the $\sin^{-4}(\theta/2)$ angular dependence confirmed by Geiger and Marsden.

2. Chadwick's discovery of the neutron (1932) completed the nuclear model. The nucleus contains $Z$ protons and $N$ neutrons, with mass number $A = Z + N$.

3. Nuclear notation (${}^{A}_{Z}\text{X}_{N}$) classifies nuclides. Isotopes share $Z$; isotones share $N$; isobars share $A$.

4. The chart of nuclides maps all ~3300 known nuclides on an $N$–$Z$ grid. The valley of stability curves toward $N > Z$ for heavy nuclei. Magic numbers ($2, 8, 20, 28, 50, 82, 126$) mark nuclei of exceptional stability.

5. Nuclear sizes from electron scattering follow $R = r_0 A^{1/3}$ with $r_0 \approx 1.21\,\text{fm}$. Because $V \propto A$, nuclear density is approximately constant: density saturation ($\rho \approx 0.16\,\text{nucleons/fm}^3$).

6. Binding energy $B = [Zm_p + Nm_n - M(A,Z)]c^2$ measures nuclear stability. The binding energy per nucleon $B/A$ peaks near $A \approx 62$ (${}^{62}\text{Ni}$) and decreases for both lighter and heavier nuclei.

7. Fusion (combining light nuclei) and fission (splitting heavy nuclei) both release energy because the products have higher $B/A$ — they are more tightly bound per nucleon than the reactants. This is the central insight of the binding energy curve.

8. Separation energies ($S_n$, $S_p$) — the energy to remove one nucleon — show sharp discontinuities at magic numbers, providing direct evidence for nuclear shell structure.

What's Next

In Chapter 2, we tackle the force that holds the nucleus together: the nuclear force. How can protons — all positively charged and crammed into a volume a few femtometers across — resist the enormous Coulomb repulsion? The answer involves a force far stronger than electromagnetism, but with a range of only $\sim 1\text{–}2\,\text{fm}$. We will examine the deuteron (the simplest bound nuclear system), nucleon-nucleon scattering data, and the basic properties of the nuclear force: its short range, its charge independence, and its spin dependence.