Unit 2: Atomic Nucleus: Composition, Binding Energy & Structural Models

The atomic nucleus is an ultra-dense, quantum many-body system bound by the strong nuclear interaction. It consists of two constituent fermionic species collectively termed nucleons:

• Protons ($p$): Baryons with positive charge $+e = +1.60217663 \times 10^{-19}\text{ C}$, rest mass $m_p = 1.0072764666\text{ u} = 938.272\text{ MeV}/c^2$, intrinsic spin $s = 1/2$, and isospin projection $T_3 = +1/2$.
• Neutrons ($n$): Electrically neutral baryons ($q = 0$), rest mass $m_n = 1.0086649158\text{ u} = 939.565\text{ MeV}/c^2$, intrinsic spin $s = 1/2$, and isospin projection $T_3 = -1/2$.

The composition of any nuclide $^{A}_{Z}\text{X}_N$ is defined by:

• Atomic number $Z$: Number of protons.
• Neutron number $N$: Number of neutrons.
• Mass number $A = Z + N$: Total nucleon count.

Nuclear Dimensions and Radius Scaling

Because nucleons are subject to the Pauli exclusion principle and the strong nuclear force exhibits a hard repulsive core at inter-nucleon separations $r < 0.5\text{ fm}$ alongside short-range saturation at $r \approx 1.0 - 1.5\text{ fm}$, nuclear matter is essentially incompressible.

Consequently, the volume of a nucleus $V$ is directly proportional to the total number of nucleons $A$:

The nuclear radius $R$ is parameterized by the foundational empirical scaling law:

where $R_0$ is the nuclear radius parameter.

• High-energy electron scattering measurements (Hofstadter, 1950s) sensitive to the nuclear charge distribution yield:

• Nuclear potential measurements (neutron scattering, alpha decay barriers) yield slightly larger values:

Nuclear Density and Incompressibility

Using $R_0 = 1.25\text{ fm} = 1.25 \times 10^{-15}\text{ m}$ and average nucleon mass $m \approx 1.67 \times 10^{-27}\text{ kg}$:

Notice that mass number $A$ cancels out completely!

The nuclear matter nucleon number density $\rho_0$ is:

This incredible density—over 200 trillion times denser than liquid water—is uniform across all nuclei from helium to uranium, confirming that nuclear matter behaves like an incompressible quantum liquid drop.

Woods-Saxon Radial Charge Distribution

High-energy elastic electron scattering demonstrates that nuclei do not possess sharp, hard-sphere boundaries. Instead, the charge density profile $\rho(r)$ is modeled accurately by the Woods-Saxon (two-parameter Fermi) distribution:

where:

• $\rho_0$ is the central interior nuclear density.
• $c$ is the half-density radius ($c \approx 1.07 A^{1/3}\text{ fm}$), where $\rho(c) = 0.5 \rho_0$.
• $a$ is the surface diffuseness parameter ($a \approx 0.54\text{ fm}$), governing the rate of density falloff.
• The skin thickness $t_{90-10}$, defined as the distance over which the density drops from $90\%$ to $10\%$ of $\rho_0$, is universally:

across virtually all stable nuclei.

One of the most striking findings of precise mass spectrometry (Aston, Bainbridge) is that the precise mass of any bound nucleus $M(A,Z)$ is strictly less than the sum of the rest masses of its constituent free protons and neutrons.

The Mass Defect $\Delta m$

The difference between the total mass of the individual constituent nucleons at infinite separation and the actual bound atomic mass is defined as the mass defect $\Delta m$:

where:

• $m(^1\text{H}) = 1.007825032\text{ u}$ is the atomic mass of neutral hydrogen-1 (accounting for electron mass $m_e$).
• $m_n = 1.008664916\text{ u}$ is the free neutron mass.
• $M(A,Z)$ is the neutral atomic mass of the nuclide.
• $1\text{ u} \equiv \frac{1}{12} m(^{12}\text{C}) = 1.66053906660 \times 10^{-27}\text{ kg} \equiv 931.49410242\text{ MeV}/c^2$.

Total Nuclear Binding Energy ($B$)

According to Einstein's mass-energy equivalence principle ($E = mc^2$), this missing mass was liberated as binding energy during the nucleosynthetic coalescence of the nucleus:

Binding Energy per Nucleon ($B/A$) Curve Morphology

To compare the relative thermodynamic stability of different nuclear species, we define the binding energy per nucleon $B/A$:

``` B/A (MeV / Nucleon) 9 ▲ ⁵⁶Fe (8.79 MeV) ⁶²Ni (8.79 MeV) │ ▲ 8 │ ⁴He / \_____________ │ ▲ / \_______ 7 │ / \ / \_____ ²³⁸U (7.57 MeV) │ / \____/ \ 6 │ / 5 │ / EXOTHERMIC FUSION EXOTHERMIC FISSION 4 │ / 3 │/ 2 │ ²H (1.11 MeV) 1 │ 0 └──┴───────────┴─────────────────────────┴────────► Mass Number A 0 56 238 ```

The universal curve of binding energy per nucleon reveals fundamental properties of nuclear forces:

1. Low Mass Region ($A < 20$):

The curve rises steeply from $1.11\text{ MeV/nucleon}$ for deuterium ($^{2}\text{H}$) up to $\approx 8\text{ MeV}$. Prominent periodic spikes appear at $^{4}\text{He}$ ($7.07\text{ MeV}$), $^{8}\text{Be}$, $^{12}\text{C}$ ($7.68\text{ MeV}$), and $^{16}\text{O}$ ($7.98\text{ MeV}$). These peaks reflect the exceptional stability of tightly bound $\alpha$-conjugate nuclei ($Z = N = 2n$).

2. The Global Peak at $A \approx 56 - 62$:

The maximum of the curve occurs in the iron-nickel group:

• $^{56}_{26}\text{Fe}$: $B/A = 8.790\text{ MeV/nucleon}$
• $^{62}_{28}\text{Ni}$: $B/A = 8.7946\text{ MeV/nucleon}$ (highest absolute binding energy per nucleon of all known nuclides)
• $^{58}_{26}\text{Fe}$: $B/A = 8.792\text{ MeV/nucleon}$

These nuclides represent the thermodynamically most stable nuclear states in the cosmos ("iron peak" in stellar nucleosynthesis).

3. High Mass Plateau and Decline ($A > 60$):

Beyond iron, the binding energy per nucleon gently and monotonically declines from $\sim 8.8\text{ MeV}$ down to $7.57\text{ MeV/nucleon}$ for uranium-238 ($^{238}\text{U}$). This decline is driven directly by long-range Coulomb electrostatic repulsion between protons, which scales quadratically with $Z^2$ and overcomes the short-range strong nuclear attraction.

4. Thermodynamic Implications for Nuclear Energy:

• Nuclear Fusion: Coalescing light nuclei ($A < 56$, e.g., $^2\text{H} + ^3\text{H} \to ^4\text{He} + n$) moves upward along the curve toward higher $B/A$, releasing millions of electron volts per reaction.
• Nuclear Fission: Splitting heavy nuclei ($A > 200$, e.g., $n + ^{235}\text{U} \to \text{Ba} + \text{Kr} + 3n$) breaks a less tightly bound nucleus ($7.6\text{ MeV/nucleon}$) into two intermediate fragments ($8.5\text{ MeV/nucleon}$), releasing $\approx 0.9\text{ MeV/nucleon} \times 235 \approx 200\text{ MeV}$ of net kinetic energy per fission.

If we plot all known nuclides on a Segrè chart (the chart of nuclides) with proton number $Z$ on the ordinate and neutron number $N$ on the abscissa, stable nuclei form a narrow ribbon termed the valley of beta-stability.

``` Z (Proton Number) 92 ▲ / N = Z line │ / 82 │ Pb / │ / Valley of Beta Stability │ / (N > Z for heavy nuclei) │ /• 40 │ Zr /•• │ /•• 20 │ Ca /••• (N ≈ Z for light nuclei) │ /•••• │ /•••• └────────────┴──────────────────────────► N (Neutron Number) 0 20 40 82 126 ```

The $N/Z$ Ratio Trajectory

1. Light Stable Nuclei ($Z \le 20$):

For light nuclei, the valley follows the line of symmetry $N = Z$ ($N/Z \approx 1.0$). Examples include $^{4}_{2}\text{He}$, $^{12}_{6}\text{C}$, $^{14}_{7}\text{N}$, $^{16}_{8}\text{O}$, $^{20}_{10}\text{Ne}$, $^{28}_{14}\text{Si}$, and $^{40}_{20}\text{Ca}$. This symmetry arises from the Pauli exclusion principle, which penalizes unequal filling of neutron and proton quantum energy levels.

2. Heavy Stable Nuclei ($Z > 20$):

As $Z$ increases, Coulomb repulsion among protons scales as $Z(Z-1)/R \propto Z^2 / A^{1/3}$. To maintain binding against this disruptive force, the nucleus must incorporate an increasing excess of neutrons, which supply attractive strong nuclear force without adding disruptive electrostatic charge. Consequently, the $N/Z$ ratio increases monotonically:

• Iron-56: $^{56}_{26}\text{Fe}_{30} \implies N/Z = 30/26 \approx 1.15$
• Silver-107: $^{107}_{47}\text{Ag}_{60} \implies N/Z = 60/47 \approx 1.28$
• Lead-208: $^{208}_{82}\text{Pb}_{126} \implies N/Z = 126/82 \approx 1.54$
• Uranium-238: $^{238}_{92}\text{U}_{146} \implies N/Z = 146/92 \approx 1.59$

Nuclear Decay Modes Relative to the Valley of Stability

Nuclides lying off the central floor of the valley undergo spontaneous radioactive transitions to reach stability:

• Neutron-Rich Nuclides (Lying below/to the right of the valley):

Excess neutrons result in negative beta decay ($\beta^-$) where a neutron converts to a proton:

For extreme neutron excess, prompt neutron emission ($S_n \le 0$) defines the neutron dripline.

• Proton-Rich Nuclides (Lying above/to the left of the valley):

Excess protons result in positron emission ($\beta^+$) or electron capture (EC):

For extreme proton excess, proton emission ($S_p \le 0$) defines the proton dripline.

• Superheavy Unstable Nuclides ($Z > 82, A > 209$):

Coulomb repulsion is so massive that beta decay cannot stabilize the nucleus; alpha decay ($\Delta Z = -2, \Delta A = -4$) and spontaneous fission become dominant. The heaviest completely stable nuclide is lead-208 ($^{208}_{82}\text{Pb}$). Bismuth-209 is very weakly alpha active ($T_{1/2} = 2.01 \times 10^{19}\text{ yr}$).

Even-Odd Nuclear Stability Systematics

The pairing interaction in nuclear forces produces dramatic patterns in stable isotope abundances:

| Nucleus Type | $Z$ | $N$ | Number of Stable Nuclides | Average Stable Isotopes per Element | Examples | | :--- | :--- | :--- | :--- | :--- | :--- | | Even-Even | Even | Even | 166 | $\sim 5.5$ | $^{12}\text{C}, ^{16}\text{O}, ^{56}\text{Fe}, ^{208}\text{Pb}$ | | Even-Odd | Even | Odd | 57 | $\sim 1.4$ | $^{13}\text{C}, ^{57}\text{Fe}, ^{207}\text{Pb}$ | | Odd-Even | Odd | Even | 53 | $\sim 1.3$ | $^{19}\text{F}, ^{23}\text{Na}, ^{63}\text{Cu}$ | | Odd-Odd | Odd | Odd | 9 (only 4 stable light) | $< 0.1$ | $^{2}\text{H}, ^{6}\text{Li}, ^{10}\text{B}, ^{14}\text{N}$ ($^{50}\text{V}, ^{138}\text{La}, ^{180m}\text{Ta}$) |

Over $60\%$ of all stable nuclides in the universe are even-even, while only four completely stable odd-odd nuclides exist in nature ($^{2}\text{H}, ^{6}\text{Li}, ^{10}\text{B}, ^{14}\text{N}$). This disparity provides direct evidence that nucleons of identical species couple into anti-parallel pairs ($J^\pi = 0^+$) with enhanced binding energy.

Nuclear species (nuclides) are classified into specific families based on relationships between their proton count $Z$, neutron count $N$, and total mass number $A$.

1. Isotopes ($\Delta Z = 0$)

Nuclides possessing the same atomic number $Z$ (identical chemical identity) but differing neutron numbers $N$ and mass numbers $A$.

• Examples:
• Hydrogen isotopes: $^{1}_{1}\text{H}$ (protium), $^{2}_{1}\text{H}$ (deuterium), $^{3}_{1}\text{H}$ (tritium).
• Uranium isotopes: $^{234}_{92}\text{U}$, $^{235}_{92}\text{U}$, $^{238}_{92}\text{U}$.
• Because their electron configurations are identical, isotopes exhibit nearly identical chemical properties (governed by the same electron shell structure), but display distinct nuclear properties (spins, magnetic moments, cross sections, decay half-lives).

2. Isobars ($\Delta A = 0$)

Nuclides possessing the same total mass number $A$ but differing atomic numbers $Z$ and neutron numbers $N$.

• Examples:
• $A = 14$: $^{14}_{6}\text{C}$, $^{14}_{7}\text{N}$, $^{14}_{8}\text{O}$.
• $A = 40$: $^{40}_{18}\text{Ar}$, $^{40}_{19}\text{K}$, $^{40}_{20}\text{Ca}$.
• Isobars possess completely different chemical identities. In any isobaric chain, Mattauch's Isobar Rule states that two neighboring stable isobars ($\Delta Z = 1$) cannot exist; one must beta decay into the other.

3. Isotones ($\Delta N = 0$)

Nuclides possessing the same neutron number $N$ but differing atomic numbers $Z$ and mass numbers $A$.

• Examples ($N = 82$, magic neutron shell):
• $^{138}_{56}\text{Ba}_{82}$, $^{139}_{57}\text{La}_{82}$, $^{140}_{58}\text{Ce}_{82}$, $^{141}_{59}\text{Pr}_{82}$, $^{142}_{60}\text{Nd}_{82}$.
• Isotones are invaluable in nuclear structure research because variations in binding energy and excitation spectra reflect changes in proton configuration against an identical neutron core.

4. Nuclear Isomers

Nuclides possessing the identical $Z$ and identical $A$ that exist in different, long-lived metastable nuclear energy states (isomeric states, denoted with an "m", e.g., $^{99m}\text{Tc}$, $^{180m}\text{Ta}$).

• A metastable excited state arises when the angular momentum difference $\Delta I$ between the excited isomer and the lower state is large (high multipolarity, e.g., $E4$ or $M4$) and the transition energy $\Delta E$ is low, severely suppressing the probability of gamma electromagnetic de-excitation (isomeric transition, IT).
• Example: Technetium-99m ($^{99m}_{43}\text{Tc}$):

The $142.6\text{ keV}$ state has spin-parity $1/2^-$ while the ground state $^{99}\text{Tc}$ has $9/2^+$. The large spin change $\Delta I = 4$ produces an exceptionally long half-life of $T_{1/2} = 6.01\text{ hours}$, making it the premier radiotracer in nuclear medicine.

• Tantalum-180m ($^{180m}_{73}\text{Ta}$, $I^\pi = 9^-$) has a half-life exceeding $4.5 \times 10^{16}\text{ years}$—longer than the age of the universe!

5. Mirror Nuclei

A pair of isobaric nuclei ($A_1 = A_2$) where the proton count of one equals the neutron count of the other:

• Examples: $(^3_1\text{H}, ^3_2\text{He})$, $(^7_3\text{Li}, ^7_4\text{Be})$, $(^{11}_5\text{B}, ^{11}_6\text{C})$, $(^{15}_7\text{N}, ^{15}_8\text{O})$.
• Because the strong nuclear interaction is charge-symmetric (the $p-p$, $n-n$, and $p-n$ strong potentials are identical in the same quantum state), the binding energy difference between mirror nuclei is due almost exclusively to the difference in Coulomb electrostatic self-energy:

Measuring this mass difference enables precise experimental determination of the nuclear charge radius $R_0$.

In 1935, Carl Friedrich von Weizsäcker developed the Liquid Drop Model, analogizing the atomic nucleus to a macroscopic droplet of incompressible charged liquid. Nucleons interact through short-range, saturating attractive strong forces analogous to van der Waals forces in a liquid, counteracted by long-range electrostatic repulsion between protons.

The resulting Semi-Empirical Mass Formula (SEMF) provides an analytical formulation for total binding energy $B(A,Z)$:

``` SEMF Binding Energy Components B(A,Z) = + Volume Term [+ a_v * A]

• Surface Term [- a_s * A^(2/3)]
• Coulomb Term [- a_c * Z(Z-1) / A^(1/3)]
• Asymmetry Term [- a_a * (A - 2Z)^2 / A]

+ Pairing Term [+ δ(A,Z)] ```

Physical Derivation of the Five Terms

1. Volume Energy Term ($+a_v A$):

Because the strong nuclear force is short-ranged and exhibits saturation, each interior nucleon interacts only with its immediate nearest neighbors ($\sim 12$ nucleons), contributing a constant binding energy independent of total drop size:

2. Surface Energy Term ($-a_s A^{2/3}$):

Nucleons residing on the nuclear surface have fewer neighboring nucleons than those in the interior, reducing the net binding. By analogy with surface tension in liquids, this deficit is proportional to the surface area of the sphere:

3. Coulomb Repulsion Term ($-a_c Z(Z-1)/A^{1/3}$):

The $Z$ protons uniformly distributed throughout the nuclear sphere repel each other electrostatically. The classical self-energy of a uniformly charged sphere of radius $R$ is:

Correcting for quantum self-interaction (a proton does not repel itself, yielding $Z(Z-1)$ pairs):

4. Asymmetry Energy Term ($-a_a (A - 2Z)^2 / A$):

Quantum mechanics (Pauli exclusion principle) dictates that protons and neutrons fill separate fermion energy wells. For a fixed total nucleon count $A$, the state of lowest total kinetic energy occurs when $N = Z = A/2$. Any neutron excess $(N - Z) = (A - 2Z)$ forces nucleons into higher unoccupied quantum states. In a Fermi gas approximation, expanding the energy difference in powers of $(N - Z)$:

5. Pairing Energy Term ($\delta(A,Z)$):

Due to the spin-orbit pairing force, identical nucleons with antiparallel spins form pairs ($J^\pi = 0^+$) with enhanced binding energy. The pairing term is:

where $a_p \approx 11.18\text{ MeV}$ (or alternatively parameterized as $a_p / A^{3/4}$).

Prediction of the Most Stable Isobar ($Z_{\text{stable}}$)

For a constant mass number $A$, the binding energy is a quadratic parabola in $Z$:

Evaluating the derivative:

Neglecting small terms, we solve for $Z_{\text{stable}}$:

• For light nuclei ($A \to 0$): $Z_{\text{stable}} \to A/2$ ($N = Z$).
• For heavy nuclei ($A = 208$): $Z_{\text{stable}} = \frac{208}{2 + 0.015(208)^{2/3}} = \frac{208}{2 + 0.015(35.1)} = \frac{208}{2.527} \approx 82.3 \implies Z = 82$ (Lead-208), in remarkable agreement with experiment!

While the Liquid Drop Model explains macroscopic binding energies and fission, it completely fails to explain quantum micro-structure:

1. Exceptional binding energy spikes at specific proton or neutron counts:

1. Nuclei with magic numbers of nucleons (such as $^{4}_{2}\text{He}_2$, $^{16}_{8}\text{O}_8$, $^{40}_{20}\text{Ca}_{20}$, $^{48}_{20}\text{Ca}_{28}$, $^{208}_{82}\text{Pb}_{126}$) are "doubly magic" and possess zero electric quadrupole moments, exceptionally high first-excited state energies, and extraordinarily small neutron capture cross sections.
2. Ground-state nuclear spins and parities ($J^\pi$) and magnetic dipole moments.

The Shell Model Potential and Mayer-Jensen Spin-Orbit Coupling

In 1949, Maria Goeppert Mayer and J. Hans D. Jensen independently recognized that nucleons move in a mean central potential created by all other nucleons, accompanied by an extraordinarily strong spin-orbit interaction ($\vec{l}\cdot\vec{s}$ coupling):

Mathematical Derivation of Spin-Orbit Level Splitting

For a nucleon with orbital angular momentum $\vec{l}$ and spin $\vec{s}$ ($s = 1/2$), the total single-particle angular momentum is:

Squaring $\vec{j}$:

Evaluating the quantum expectation values:

With $s = 1/2 \implies s(s+1) = 3/4$:

1. For state with $j = l + 1/2$:

1. For state with $j = l - 1/2$:

The energy splitting between the two sub-states is:

Crucially, in the nuclear interaction (unlike atomic atomic fine structure), the spin-orbit potential $V_{ls}(r)$ is negative. Therefore, the state with higher total angular momentum $j = l + 1/2$ is pushed downward in energy, while $j = l - 1/2$ is pushed upward!

``` Harmonic Oscillator With Spin-Orbit Splitting Magic Shell Levels (V_ls < 0) Closure ──────────────────────────────────────────────────────────────── 1f (l=3) ────────── 1f_5/2 (6 states) ─────── \ \──────── 1f_7/2 (8 states) ─────── ► [28] MAGIC ──────────────────────────────────────────────────────────────── 1d (l=2) ────────── 1d_3/2 (4 states) ─────── ► [20] MAGIC 2s (l=0) ────────── 2s_1/2 (2 states) 1d (l=2) ────────── 1d_5/2 (6 states) ──────────────────────────────────────────────────────────────── 1p (l=1) ────────── 1p_1/2 (2 states) ─────── ► [8] MAGIC ────────── 1p_3/2 (4 states) ──────────────────────────────────────────────────────────────── 1s (l=0) ────────── 1s_1/2 (2 states) ─────── ► [2] MAGIC ```

For large orbital angular momentum $l$ ($l = 3$ for $f$-states, $l = 4$ for $g$-states), this downward shift is so massive that the $j = l + 1/2$ level drops completely across the oscillator gap into the major shell below. This accounts for every observed magic number:

• $1s_{1/2}$ (capacity 2) $\implies \mathbf{2}$
• $1p_{3/2}, 1p_{1/2}$ (capacity $4 + 2 = 6$) $\implies 2 + 6 = \mathbf{8}$
• $1d_{5/2}, 2s_{1/2}, 1d_{3/2}$ (capacity $6 + 2 + 4 = 12$) $\implies 8 + 12 = \mathbf{20}$
• Intruder state $1f_{7/2}$ pushed down (capacity 8) $\implies 20 + 8 = \mathbf{28}$
• $2p_{3/2}, 1f_{5/2}, 2p_{1/2}$ plus intruder $1g_{9/2}$ (capacity 22) $\implies 28 + 22 = \mathbf{50}$
• $2d_{5/2}, 1g_{7/2}, 1h_{11/2}$, etc. $\implies \mathbf{82}$
• Intruder $1i_{13/2}$, etc. $\implies \mathbf{126}$

The spherical Shell Model succeeds near magic numbers, but fails in the mid-shell regions ($150 < A < 190$ and $A > 220$). In these regions, nuclei exhibit:

• Electric quadrupole moments $Q_0$ up to 30 times larger than the single-particle limit.
• Rotational band excitation spectra following exact $E_J \propto J(J+1)$ sequences.
• Strongly enhanced collective electric quadrupole transition rates $B(E2)$.

The Bohr-Mottelson Unified Collective Model (1953)

Aage Bohr, Ben Mottelson, and James Rainwater resolved this disparity by modeling the nucleus as a deformed, non-spherical liquid drop whose collective surface vibrations and rotations couple to individual valence nucleon orbits.

The deformed nuclear surface is parameterized via spherical harmonics $Y_{\lambda \mu}(\theta, \phi)$:

For axially symmetric quadrupole deformations ($\lambda = 2, \mu = 0$), this simplifies in terms of the deformation parameter $\beta$:

• $\beta > 0$: Prolate spheroid (cigar-shaped, major axis along symmetry axis; predominant in nature).
• $\beta < 0$: Oblate spheroid (doorknob / disk-shaped).

Rotational Energy Spectra

For an even-even deformed nucleus with ground state $J^\pi = 0^+$, quantum rotation of the collective core about an axis perpendicular to the symmetry axis yields kinetic rotational energy:

where $\mathcal{I}$ is the effective moment of inertia of the deformed nucleus. The excitation energy ratios in a pure rotational band follow:

Experimental spectra for deformed nuclei (e.g., $^{160}\text{Gd}, ^{174}\text{Yb}, ^{238}\text{U}$) match this $3.33$ ratio precisely, confirming collective rotation.

The Fermi Gas Model of Nuclear Matter

For high-energy nuclear reactions and statistical level densities, the nucleus is treated as a degenerate quantum gas of non-interacting fermions (protons and neutrons) confined within a spherical volume $V = \frac{4}{3}\pi R^3$.

The number of spatial momentum states in phase space volume $d^3r \, d^3p$ is:

where the factor 2 accounts for spin degeneracy ($s_z = \pm 1/2$). Integrating from $p = 0$ to the Fermi momentum $p_F$:

Substituting nuclear matter density $\rho = N / V$:

For symmetric nuclear matter ($N = Z = A/2$) with $\rho_0 \approx 0.17\text{ fm}^{-3}$:

The Fermi Energy $E_F$ of a nucleon is:

The average kinetic energy per nucleon in the ground state is:

Adding the average nucleon binding energy ($B/A \approx 8\text{ MeV}$), the total depth of the nuclear potential well is:

In 1932, immediately following Chadwick's discovery of the neutron, Werner Heisenberg proposed that the proton and neutron are not fundamentally distinct particles, but two different charge states of a single underlying physical entity: the nucleon.

The Isospin Vector Formalism

Heisenberg introduced the concept of isotopic spin (or isospin, denoted by vector $\vec{T}$ or $\vec{I}$), formulating a mathematical analogy with quantum mechanical spin-1/2: The nucleon possesses total isospin $T = 1/2$. In an abstract three-dimensional "isospace":

• Proton state: Third component projection $T_3 = +1/2$ (or $-1/2$ by nuclear convention):

• Neutron state: Third component projection $T_3 = -1/2$ (or $+1/2$):

For a nucleus composed of $Z$ protons and $N$ neutrons, the total third isospin component $T_3$ is strictly fixed by its composition:

The total isospin quantum number $T$ can take any integer or half-integer value satisfying:

Charge Symmetry and Charge Independence

High-energy scattering experiments demonstrate that the strong nuclear force obeys two profound symmetries:

1. Charge Symmetry: The proton-proton ($p-p$) strong potential equals the neutron-neutron ($n-n$) strong potential in the identical spatial and spin state ($V_{pp} = V_{nn}$).

2. Charge Independence: The proton-neutron ($p-n$) strong interaction equals the $p-p$ and $n-n$ interactions in identical $T = 1$ states ($V_{pp} = V_{nn} = V_{pn}|_{T=1}$).

Mathematically, the nuclear Hamiltonian $\hat{H}_{\text{strong}}$ commutes with the total isospin operator:

Consequently, states with identical total isospin $T$ but differing projections $T_3$ form an isobaric multiplet (or isospin analog states) with nearly identical nuclear structure, differing in energy solely due to electromagnetic Coulomb perturbation:

Eugene Wigner expanded this into Wigner Supermultiplet Theory, classifying nuclear energy levels according to $SU(4)$ symmetry combining space, spin, and isospin degrees of freedom.

Comparative Systematics of Semi-Empirical Mass Formula (SEMF) Parameter Sets

Different empirical fits to experimental nuclear masses yield slightly varying coefficient sets (all in MeV):

| Parameter Set | Volume $a_v$ | Surface $a_s$ | Coulomb $a_c$ | Asymmetry $a_a$ | Pairing $a_p$ | RMS Mass Error | | :--- | :--- | :--- | :--- | :--- | :--- | :--- | | von Weizsäcker (1935) | $15.75$ | $17.80$ | $0.711$ | $23.70$ | $11.18$ | $\sim 2.5\text{ MeV}$ | | Green (1954) | $15.75$ | $17.80$ | $0.710$ | $23.69$ | $12.00$ | $\sim 2.1\text{ MeV}$ | | Myers & Swiatecki (1966) | $15.68$ | $18.56$ | $0.717$ | $28.07$ | $11.00$ | $\sim 1.8\text{ MeV}$ | | Möller et al. (FRDM, 1995) | $16.25$ | $19.00$ | $0.730$ | $32.00$ | $11.50$ | $\sim 0.67\text{ MeV}$ |

Modern Finite-Range Droplet Models (FRDM) incorporate microscopic Strutinsky shell corrections and macroscopic nuclear surface diffuseness, predicting nuclear binding energies for thousands of exotic isotopes out to the neutron and proton driplines.