Nuclear Forces, Meson Exchange Theory & Nucleon-Nucleon Symmetries

1. Phenomenological Characteristics of the Strong Nuclear Interaction

The force binding protons and neutrons together within atomic nuclei is the strongest fundamental interaction in nature, yet it displays intricate characteristics that distinguish it sharply from the inverse-square gravitational and Coulomb forces:

1. Short Range ($R \approx 1\text{ to }2\text{ fm}$): The nuclear force is negligible at inter-nucleon separations $r > 3\text{ fm}$. Its effective range is governed by the mass of the lightest mediating exchange boson (the pion, $\lambdabar_\pi = \hbar/m_\pi c \approx 1.4\text{ fm}$).
2. Extremely High Strength: Within its attractive well ($r \sim 1\text{ fm}$), the nuclear interaction is approximately $100$ times stronger than the electrostatic Coulomb repulsion between two protons and $10^{38}$ times stronger than Newtonian gravity.
3. Nuclear Saturation: The binding energy per nucleon ($B/A$) is approximately constant across all stable nuclei with $A > 16$, hovering near $B/A \approx 8.0 \pm 0.5\text{ MeV/nucleon}$, and the nuclear interior density is constant ($\rho_0 \approx 0.16\text{ nucleons/fm}^3$). If every nucleon interacted with all other $A-1$ nucleons, $B$ would scale as $A(A-1) \approx A^2$. The linear scaling $B \propto A$ demonstrates that a nucleon interacts only with its immediate nearest neighbors—a property termed saturation.
4. Charge Independence: Once electromagnetic Coulomb interactions are subtracted, the nuclear force between two protons ($p$-$p$), two neutrons ($n$-$n$), and a proton-neutron pair ($n$-$p$) in the same orbital, spin, and spatial quantum state is identical to within $1\%$: $$V_{pp}^{\text{nucl}} \approx V_{nn}^{\text{nucl}} \approx V_{np}^{\text{nucl}} \quad (\text{for identical } L, S, T)$$
5. Charge Symmetry: The interaction between two protons is identical to the interaction between two neutrons ($V_{pp}^{\text{nucl}} = V_{nn}^{\text{nucl}}$). Charge symmetry is a specific subgroup of the broader charge independence.
6. Spin and Isospin Dependence: As demonstrated in the deuteron ($B = 2.22\text{ MeV}$ for $S = 1$, unbound for $S = 0$), the interaction is strongly dependent on the relative spin and isospin alignments.
7. Non-Central Tensor Force ($S_{12}$): The force depends on the spatial orientation of the nucleon spins relative to the inter-nucleon separation vector $\vec{r}$, producing $D$-state mixing and non-zero electric quadrupole moments.
8. Velocity-Dependent Spin-Orbit Force ($\vec{L}\cdot\vec{S}$): At short distances, relativistic meson exchange generates a powerful spin-orbit interaction $\vec{L}\cdot\vec{S}$, which is responsible for the single-particle shell model magic numbers.
9. Repulsive Hard Core ($r \le 0.4\text{ to }0.5\text{ fm}$): At extremely short distances, the potential transitions from strong attraction to an impenetrable repulsive barrier, preventing nuclear collapse.

1. Invariance Principles and Conservation Laws

The nuclear Hamiltonian $\hat{H}$ is governed by fundamental spacetime and internal symmetries:

• Translational and Rotational Invariance: Guarantees conservation of total linear momentum $\vec{P}$ and total angular momentum $\vec{J} = \vec{L} + \vec{S}$. The potential can depend only on the relative separation $r = |\vec{r}_1 - \vec{r}_2|$, relative momentum $\vec{p}$, and relative spins $\vec{s}_1, \vec{s}_2$.
• Space Inversion Invariance (Parity $\hat{P}$): Strong nuclear interactions strictly conserve parity: $[\hat{H}, \hat{P}] = 0$. Under parity, $\vec{r} \to -\vec{r}$ and $\vec{p} \to -\vec{p}$, while spin vectors are pseudovectors ($\vec{s} \to +\vec{s}$). Consequently, any potential term linear in $\vec{\sigma}\cdot\vec{r}$ or $\vec{\sigma}\cdot\vec{p}$ is parity-violating and strictly excluded from the strong nuclear force. Allowed terms must be scalars under coordinate inversion: $r^2$, $\vec{\sigma}_1\cdot\vec{\sigma}_2$, $\vec{L}\cdot\vec{S}$, and $S_{12}$.
• Time-Reversal Invariance ($\hat{T}$): Under time reversal, $t \to -t$, $\vec{r} \to +\vec{r}$, $\vec{p} \to -\vec{p}$, $\vec{s} \to -\vec{s}$, and $\vec{L} \to -\vec{L}$. Terms like $(\vec{\sigma}_1\times\vec{\sigma}_2)\cdot\vec{r}$ change sign under $\hat{T}$ and are forbidden in time-reversal invariant nuclear potentials.

2. The Generalized Pauli Principle

According to Werner Heisenberg's isospin concept, the proton and neutron are not two distinct particles, but two different charge states of a single entity: the nucleon with isospin $I = 1/2$.

Because nucleons are fermions with intrinsic spin $1/2$ and isospin $1/2$, the generalized Pauli exclusion principle mandates that the total wavefunction of any two-nucleon system must be antisymmetric under the simultaneous exchange of all spatial, spin, and isospin coordinates:

$$\hat{P}_{\text{total}} \Psi(1, 2) = \hat{P}_{\text{space}} \hat{P}_{\text{spin}} \hat{P}_{\text{isospin}} \Psi(1, 2) = -\Psi(1, 2)$$

Evaluating the eigenvalues of the permutation operators:

• Spatial Exchange ($\hat{P}_{\text{space}}$): $\hat{P}_r Y_{LM}(\hat{r}) = (-1)^L Y_{LM}(\hat{r})$. Symmetric for even $L$ ($+1$), antisymmetric for odd $L$ ($-1$).
• Spin Exchange ($\hat{P}_{\text{spin}}$): Triplet ($S=1$) is symmetric ($+1$); singlet ($S=0$) is antisymmetric ($-1$). Eigenvalue: $(-1)^{S+1}$.
• Isospin Exchange ($\hat{P}_{\text{isospin}}$): Triplet ($I=1$) is symmetric ($+1$); singlet ($I=0$) is antisymmetric ($-1$). Eigenvalue: $(-1)^{I+1}$.

Combining these permutation eigenvalues yields the fundamental selection rule:

$$(-1)^L \times (-1)^{S+1} \times (-1)^{I+1} = -1 \implies (-1)^{L + S + I} = -1$$

Thus, for any physically allowed two-nucleon state:

$$L + S + I = \text{odd integer}$$

Classification of allowed two-nucleon states:

• If $I = 0$ (like the deuteron): $L + S$ must be odd. If $S = 1$ (triplet), $L$ must be even ($L = 0, 2 \implies {}^3S_1, {}^3D_1$ allowed!). If $S = 0$ (singlet), $L$ must be odd ($L = 1 \implies {}^1P_1$ allowed).
• If $I = 1$ (two protons $pp$, two neutrons $nn$, or $np$ triplet): $L + S$ must be even. If $S = 0$ (singlet), $L$ must be even ($L = 0 \implies {}^1S_0$ allowed!). If $S = 1$ (triplet), $L$ must be odd ($L = 1 \implies {}^3P_{0, 1, 2}$ allowed).

1. The Isospin Operators and Algebra

In analogy to ordinary spin-$1/2$ angular momentum, isospin space is a 2D complex vector space spanned by the proton and neutron state vectors:

$$|p\rangle = \begin{pmatrix} 1 \\ 0 \end{pmatrix} \quad (I_3 = +1/2), \qquad |n\rangle = \begin{pmatrix} 0 \\ 1 \end{pmatrix} \quad (I_3 = -1/2)$$

The isospin vector operator $\vec{I} = \frac{1}{2}\vec{\tau}$ satisfies the standard $SU(2)$ Lie algebra:

$$[I_i, I_j] = i \epsilon_{ijk} I_k, \quad [\tau_i, \tau_j] = 2i \epsilon_{ijk} \tau_k$$

where $\tau_x, \tau_y, \tau_z$ are the Pauli matrices in isospin space:

$$\tau_x = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}, \quad \tau_y = \begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}, \quad \tau_z = \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}$$

The electric charge operator of a nucleon is related to isospin projection by:

$$\hat{Q} = e\left( I_3 + \frac{1}{2} \right) = e\left( \frac{\tau_z + 1}{2} \right)$$

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

$$I_3 = \sum_{i=1}^A I_{3,i} = \frac{Z - N}{2}$$

2. Two-Nucleon Isospin Multiplets

Coupling the isospins of two nucleons ($\vec{I} = \vec{I}_1 + \vec{I}_2$):

• Isospin Triplet ($I = 1$): Symmetrical isospin states: $$|1, +1\rangle = |pp\rangle \quad (Q = +2e)$$ $$|1, 0\rangle = \frac{1}{\sqrt{2}}(|pn\rangle + |np\rangle) \quad (Q = +1e)$$ $$|1, -1\rangle = |nn\rangle \quad (Q = 0)$$
• Isospin Singlet ($I = 0$): Antisymmetrical isospin state: $$|0, 0\rangle = \frac{1}{\sqrt{2}}(|pn\rangle - |np\rangle) \quad (Q = +1e, \text{ the deuteron})$$

The scalar product $\vec{\tau}_1\cdot\vec{\tau}_2$ has distinct eigenvalues:

$$\vec{I}^2 = (\vec{I}_1 + \vec{I}_2)^2 = \vec{I}_1^2 + \vec{I}_2^2 + 2\vec{I}_1\cdot\vec{I}_2 = \frac{3}{4} + \frac{3}{4} + \frac{1}{2}\vec{\tau}_1\cdot\vec{\tau}_2 = \frac{3}{2} + \frac{1}{2}\vec{\tau}_1\cdot\vec{\tau}_2$$ $$\vec{\tau}_1\cdot\vec{\tau}_2 = 2 I(I+1) - 3 = \begin{cases} +1, & I = 1 \text{ (triplet)} \\ -3, & I = 0 \text{ (singlet)} \end{cases}$$

Under exact charge independence, the strong nuclear Hamiltonian commutes with total isospin: $[\hat{H}_{\text{strong}}, \vec{I}] = 0$. Consequently, the three members of the $I = 1$ multiplet (${}^2\text{He}$, ${}^2\text{H}^*$, ${}^2n$) have identical strong interaction energy levels. The tiny observed mass splittings ($\sim \text{few MeV}$) are caused entirely by electromagnetic Coulomb interactions and the small quark mass difference ($m_d - m_u \approx 3\text{ MeV}$).

1. The Fundamental Theory: Quantum Chromodynamics (QCD)

At the deepest fundamental level, nucleons are not point particles, but composite hadrons containing valence quarks bound by massless gauge bosons called gluons:

$$|p\rangle = |uud\rangle, \qquad |n\rangle = |udd\rangle$$

The fundamental strong interaction is governed by Quantum Chromodynamics (QCD), an $SU(3)_C$ non-Abelian Yang-Mills gauge theory based on color charge (red, green, blue). Gluons themselves carry color charge and undergo self-interactions, leading to two defining features:

1. Asymptotic Freedom: At ultra-high energy scales ($Q^2 \to \infty$) or extremely short distances ($r \ll 0.1\text{ fm}$), the effective strong coupling constant $\alpha_s(Q^2) \to 0$, and quarks behave as quasi-free particles (Gross, Politzer, Wilczek, Nobel Prize 2004).
2. Color Confinement: At low energies / nuclear distance scales ($r \sim 1\text{ fm}$), $\alpha_s \sim \mathcal{O}(1)$. The color flux lines between quarks form a narrow flux tube (string) with constant string tension $\kappa \approx 1\text{ GeV/fm}$. Free color charges cannot be isolated: all observable asymptotic physical particles must be color-singlet states ($\varepsilon_{ijk}q_i q_j q_k$ for baryons, $\delta_{ij} q_i \bar{q}_j$ for mesons).

2. The Nuclear Force as a Residual Interaction

Because nucleons are color singlets (color neutral), there is no net long-range color field outside the nucleon boundary, exactly as neutral atoms have no net electric charge.

The nuclear force binding protons and neutrons is the residual strong interaction, mathematically and physically analogous to the van der Waals molecular force in chemistry:

$$\text{Color Interaction between Quarks (QCD)} \iff \text{Coulomb Interaction between Nucleus \& Electrons (QED)}$$ $$\text{Residual Nuclear Force between Nucleons} \iff \text{Van der Waals Dipole Force between Neutral Molecules}$$

Due to the spontaneous breaking of chiral symmetry in the QCD vacuum ($SU(2)_L \times SU(2)_R \to SU(2)_V$), the pion ($\pi$) emerges as the pseudo-Goldstone boson. Its exceptionally small mass ($m_\pi \approx 140\text{ MeV} \ll m_N \approx 940\text{ MeV}$) makes virtual pion exchange the dominant mechanism governing the long-range residual force between nucleons.

1. Relativistic Field Equation for a Massive Scalar Field

In 1935, Hideki Yukawa proposed that nuclear forces arise from the exchange of an unknown massive quantum field $\phi(\vec{r}, t)$. A massive scalar particle obeys the relativistic Klein-Gordon equation:

$$\left( \frac{1}{c^2}\frac{\partial^2}{\partial t^2} - \nabla^2 + \mu^2 \right) \phi(\vec{r}, t) = -\frac{g}{\varepsilon_0} \rho(\vec{r}, t)$$

where $\mu \equiv \frac{m c}{\hbar}$ is the inverse Compton wavelength of the mediating boson, $g$ is the strong coupling constant, and $\rho(\vec{r})$ is the nucleon source density.

For static interactions ($\frac{\partial\phi}{\partial t} = 0$) generated by a point-like nucleon source at the origin ($\rho(\vec{r}) = \delta(\vec{r})$):

$$\left( \nabla^2 - \mu^2 \right) \phi(\vec{r}) = 0 \quad (\text{for } r > 0)$$

In spherically symmetric coordinates:

$$\frac{1}{r}\frac{d^2}{dr^2}(r\phi) - \mu^2 \phi = 0 \implies \frac{d^2}{dr^2}(r\phi) = \mu^2 (r\phi)$$

The general solution is:

$$r\phi(r) = C_1 e^{-\mu r} + C_2 e^{+\mu r}$$

Requiring that the field vanish as $r \to \infty$ forces $C_2 = 0$. As $r \to 0$, the field must approach the unshielded Coulomb-like source $\lim_{r\to 0}\phi(r) = \frac{g}{4\pi r}$, which fixes $C_1 = \frac{g}{4\pi}$.

2. The Classic Yukawa Potential

The interaction potential energy $V(r) = -g \phi(r)$ between two nucleons separated by distance $r$ is:

$$V_{\text{Yukawa}}(r) = -\frac{g^2}{4\pi} \frac{e^{-\mu r}}{r} = -g^2 \frac{e^{-r/R}}{r}$$

where the range of the nuclear interaction $R$ is:

$$R = \frac{1}{\mu} = \frac{\hbar}{m c}$$

Using the experimentally known nuclear force range $R \approx 1.4\text{ fm}$:

$$m c^2 = \frac{\hbar c}{R} = \frac{197.3\text{ MeV}\cdot\text{fm}}{1.4\text{ fm}} \approx 141\text{ MeV}$$

Yukawa thus famously predicted the existence of a meson with mass $\approx 270\text{ }m_e$. The discovery of the pion in cosmic rays (Powell, Occhialini, Lattes, 1947) with mass $m_{\pi^\pm} = 139.6\text{ MeV}$ brilliantly confirmed Yukawa's meson exchange hypothesis (Nobel Prize 1949).

1. Phenomenological Exchange Potentials

To account for nuclear saturation and the anomalous backward peaking observed in high-energy $n$-$p$ scattering, Eugene Wigner, Ettore Majorana, James Bartlett, and Werner Heisenberg introduced four fundamental exchange operators acting on the two-nucleon wavefunction $\psi(\vec{r}_1, \vec{r}_2, \vec{s}_1, \vec{s}_2)$:

1. Wigner Force (Ordinary Non-Exchange Force, $\hat{P}_W = \hat{\mathbf{1}}$): Leaves spatial and spin coordinates unaltered: $$V_W(r) = V_0(r) \hat{\mathbf{1}}$$
2. Majorana Force (Space-Exchange Force, $\hat{P}_M \equiv \hat{P}_r$): Exchanges the spatial coordinates of the two nucleons ($\vec{r}_1 \leftrightarrow \vec{r}_2 \implies \vec{r} \to -\vec{r}$): $$\hat{P}_r \psi(\vec{r}) = \psi(-\vec{r}) = (-1)^L \psi(\vec{r})$$ Its eigenvalue is $+1$ for even $L$ ($S, D$ states) and $-1$ for odd $L$ ($P, F$ states). The Majorana interaction is $V_M(r) \hat{P}_r$.
3. Bartlett Force (Spin-Exchange Force, $\hat{P}_B \equiv \hat{P}_\sigma$): Exchanges the spin states of the two nucleons ($\vec{s}_1 \leftrightarrow \vec{s}_2$): $$\hat{P}_\sigma = \frac{1 + \vec{\sigma}_1\cdot\vec{\sigma}_2}{2}$$ Because $\vec{\sigma}_1\cdot\vec{\sigma}_2 = +1$ for triplet ($S=1$) and $-3$ for singlet ($S=0$): $$\hat{P}_\sigma |S=1\rangle = +|S=1\rangle, \qquad \hat{P}_\sigma |S=0\rangle = -|S=0\rangle$$ Its eigenvalue is $(-1)^{S+1}$. The Bartlett interaction is $V_B(r) \hat{P}_\sigma$.
4. Heisenberg Force (Space and Spin Exchange Force, $\hat{P}_H \equiv \hat{P}_r \hat{P}_\sigma$): Exchanges both spatial and spin coordinates simultaneously: $$\hat{P}_H = \hat{P}_r \hat{P}_\sigma = -\hat{P}_\tau$$ Its eigenvalue is $(-1)^L (-1)^{S+1} = -(-1)^{I+1}$. The Heisenberg interaction is $V_H(r) \hat{P}_H$.

2. The Composite Central Potential

The general central nucleon-nucleon potential is a linear combination of all four operators:

$$V_{\text{central}}(r) = V_W(r) + V_M(r) \hat{P}_r + V_B(r) \hat{P}_\sigma + V_H(r) \hat{P}_r \hat{P}_\sigma$$

Evaluating the eigenvalues in the four basic two-nucleon spin-parity channels:

State (${}^{2S+1}L_J$)	$L$ Parity	Spin $S$	$\hat{P}_r$	$\hat{P}_\sigma$	$\hat{P}_H$	Effective Potential $V_{\text{eff}}(r)$
${}^3S_1$ (Triplet Even)	Even ($L=0$)	1	$+1$	$+1$	$+1$	$V_W + V_M + V_B + V_H$
${}^1S_0$ (Singlet Even)	Even ($L=0$)	0	$+1$	$-1$	$-1$	$V_W + V_M - V_B - V_H$
${}^3P$ (Triplet Odd)	Odd ($L=1$)	1	$-1$	$+1$	$-1$	$V_W - V_M + V_B - V_H$
${}^1P$ (Singlet Odd)	Odd ($L=1$)	0	$-1$	$-1$	$+1$	$V_W - V_M - V_B + V_H$

Nuclear saturation requires that the nuclear force be repulsive in odd-$L$ states while attractive in even-$L$ states, which is achieved when the Majorana exchange component dominates: $V_M \gg V_W, V_B, V_H$ (Serber force mixture).

1. The One-Boson-Exchange (OBE) Model

A single pion exchange explains only the long-range tail ($r > 2\text{ fm}$) of the nuclear force. A complete microscopic description requires summing exchanges of all low-lying meson states:

1. One-Pion Exchange Potential (OPEP, $r > 2\text{ fm}$): Pions ($\pi^0, \pi^\pm$) are pseudoscalar mesons ($J^P = 0^-$) with isospin $I = 1$ and mass $m_\pi \approx 138\text{ MeV}$. OPEP produces the tensor force and long-range spin-spin interaction: $$V_{\text{OPEP}}(r) = \frac{g_{\pi NN}^2}{12\pi}\left( \frac{m_\pi}{M} \right)^2 (\vec{\tau}_1\cdot\vec{\tau}_2) \left[ (\vec{\sigma}_1\cdot\vec{\sigma}_2) + S_{12}\left(1 + \frac{3}{m_\pi r} + \frac{3}{(m_\pi r)^2}\right) \right] \frac{e^{-m_\pi r}}{r}$$
2. Scalar Meson Exchange ($\sigma / f_0(500)$, $0.8\text{ fm} < r < 2.0\text{ fm}$): A scalar meson ($J^P = 0^+$) with isospin $I = 0$ and broad mass $m_\sigma \approx 500\text{ to }600\text{ MeV}$ (physically representing two correlated pions in an $s$-wave state, $\pi\pi$). It provides the powerful intermediate-range central attraction responsible for binding atomic nuclei.
3. Vector Meson Exchange ($\rho, \omega$, $r < 0.8\text{ fm}$):
   • $\omega$ meson ($J^P = 1^-, I = 0, m_\omega \approx 782\text{ MeV}$): Yields strong short-range central repulsion (generating the repulsive hard core) and a large Thomas-like spin-orbit force $\vec{L}\cdot\vec{S}$.
   • $\rho$ meson ($J^P = 1^-, I = 1, m_\rho \approx 770\text{ MeV}$): Provides isospin-dependent vector exchange, modifying the tensor force at short distances with opposite sign to OPEP.

2. Modern High-Precision Realistic Potentials

State-of-the-art nucleon-nucleon potentials fit thousands of $p$-$p$ and $n$-$p$ scattering data points up to $350\text{ MeV}$ with $\chi^2/\text{datum} \approx 1.0$:

• Argonne $v_{18}$ (AV18): Contains 18 operator terms, including 14 charge-independent terms (central, tensor, spin-orbit, quadratic spin-orbit) plus 4 charge-dependent and charge-symmetry-breaking terms.
• CD-Bonn: A non-local relativistic One-Boson-Exchange potential based on covariant Feynman amplitudes.
• Chiral Effective Field Theory ($\chi$EFT): Modern approach based directly on the symmetries of QCD, expanding the nuclear force systematically in powers of $(Q/\Lambda_\chi)$ (Weinberg, van Kolck, Epelbaum, Machleidt). It naturally predicts three-nucleon forces ($3NF$) at next-to-next-to-leading order (NNLO).