Chapter 31 — The Standard Model and Nuclear Physics: Quarks, Gluons, and the Strong Force

"The theory of the strong interactions is both the simplest and the most difficult of the fundamental theories. Simplest because the Lagrangian can be written on a single line. Most difficult because that Lagrangian cannot be solved analytically at the energy scales that matter most for nuclear physics." — Frank Wilczek, Nobel Lecture (2004)

Throughout this textbook, we have treated the proton and neutron as the fundamental building blocks of nuclei. We modeled the nuclear force as an interaction between nucleons — first through Yukawa's meson exchange (Chapter 3), then through the shell model's mean field (Chapter 6), and ultimately through the chiral effective field theory framework that connects nuclear forces to the symmetries of the underlying theory. We have used the weak interaction to understand beta decay (Chapter 14) and the electromagnetic interaction to probe nuclear structure (Chapter 9).

But the proton and neutron are not fundamental particles. They are composite objects, built from quarks bound together by gluons — the carriers of the strong force described by quantum chromodynamics (QCD). The nuclear force that binds protons and neutrons into nuclei is not a fundamental interaction at all. It is a residual effect of the strong force between quarks, in precisely the same way that the van der Waals force between electrically neutral atoms is a residual effect of the electromagnetic force between charged particles.

This chapter connects nuclear physics to its foundations in particle physics. We begin with the quark model of hadrons, develop the essentials of QCD, explain how the nuclear force emerges from QCD at low energies, describe the modern tools (lattice QCD and chiral EFT) that bridge the gap, and conclude with the internal structure of the nucleon itself — including two puzzles that have driven experimental programs for decades.

31.1 The Quark Model of Hadrons

A Brief History: From the Eightfold Way to Quarks

By the early 1960s, particle physics faced an embarrassment of riches. Accelerator experiments had discovered dozens of strongly interacting particles — collectively called hadrons — with a bewildering variety of masses, spins, and decay patterns. The situation was sometimes compared to atomic spectroscopy before the periodic table: a catalogue of data without an organizing principle.

The breakthrough came from Murray Gell-Mann (independently, George Zweig) in 1964, who proposed that all hadrons are composite objects built from a small number of fundamental constituents called quarks. Gell-Mann's scheme grew out of his earlier "Eightfold Way" classification (1961), which organized hadrons into multiplets of the flavor symmetry group SU(3) — the same mathematical structure as the irreducible representations of three fundamental objects.

The Six Quark Flavors

We now know that nature provides six quark flavors, organized into three generations:

Each quark carries spin $1/2$ (they are fermions) and baryon number $B = 1/3$. Each quark has a corresponding antiquark ($\bar{u}$, $\bar{d}$, etc.) with opposite charge and baryon number.

For nuclear physics, the first generation — up and down quarks — dominates overwhelmingly. The proton and neutron are built almost entirely from $u$ and $d$ quarks. The strange quark is relevant for kaons, hyperons, and the strange-quark content of the nucleon (a subtle effect we discuss in Section 31.7). The heavier quarks ($c$, $b$, $t$) are important for particle physics but play essentially no direct role in nuclear structure at ordinary energies.

A crucial and initially puzzling fact: the $u$ and $d$ quark masses ($\sim 2$--$5$ MeV) are tiny compared to the proton mass ($938$ MeV). The quarks account for less than 2% of the proton's mass. Where does the rest come from? We will return to this question in Section 31.3 — the answer is one of the most profound results in modern physics.

Baryons: Three-Quark States

Baryons are hadrons with baryon number $B = 1$, composed of three quarks ($qqq$). The two most important baryons — the ones that constitute all nuclear matter — are:

$$p = uud \qquad (\text{proton: charge } +1, \text{ mass } 938.3 \text{ MeV}/c^2)$$ $$n = udd \qquad (\text{neutron: charge } 0, \text{ mass } 939.6 \text{ MeV}/c^2)$$

The electric charges follow directly from the quark charges:

$$Q_p = \frac{2}{3} + \frac{2}{3} - \frac{1}{3} = +1$$ $$Q_n = \frac{2}{3} - \frac{1}{3} - \frac{1}{3} = 0$$

This is not a coincidence — it is the quark model working.

The proton-neutron mass difference, $m_n - m_p = 1.293$ MeV/$c^2$, arises from the $d$-$u$ mass difference ($m_d - m_u \approx 3$ MeV/$c^2$) partially offset by the electromagnetic self-energy of the proton (the proton's charge increases its electromagnetic mass relative to the neutral neutron). The fact that the neutron is heavier than the proton — by just $0.14\%$ — is essential for the stability of hydrogen and thus for the existence of chemistry and biology as we know them. If the proton were heavier, it would beta-decay into a neutron, and stable hydrogen atoms could not exist.

Other baryons important for nuclear and particle physics include:

The $\Delta(1232)$ resonance deserves special attention for nuclear physics. As we discussed in Chapter 3, the virtual excitation of a nucleon to a $\Delta$ is the dominant mechanism generating the three-nucleon force. The $\Delta$ has the same quark content as the nucleon but with all three quark spins aligned ($J = 3/2$ rather than $J = 1/2$), and it is 293 MeV heavier than the nucleon. The $\Delta$-nucleon mass splitting provides the natural energy scale for three-nucleon force effects.

The $\Omega^-$ ($sss$) baryon holds a special place in history. Its existence was predicted by Gell-Mann's SU(3) classification before it was discovered experimentally in 1964 at Brookhaven National Laboratory — a triumph of the quark model comparable to Mendeleev's prediction of undiscovered elements.

Mesons: Quark-Antiquark States

Mesons are hadrons with baryon number $B = 0$, composed of a quark-antiquark pair ($q\bar{q}$). The mesons most relevant to nuclear physics are:

The pions are by far the most important mesons for nuclear physics. As the lightest hadrons, they mediate the longest-range component of the nuclear force. Their unusually low mass — much lighter than any other hadron — is not accidental. It is a consequence of the approximate chiral symmetry of QCD, which we discuss in Section 31.5. If the $u$ and $d$ quark masses were exactly zero, the pions would be exactly massless Goldstone bosons. Their small but nonzero mass ($\sim 140$ MeV) reflects the small but nonzero quark masses.

🔄 Check Your Understanding: Verify the electric charges of the $\pi^+$, $\pi^-$, and $\pi^0$ from their quark content. Why is the $\pi^0$ a superposition of $u\bar{u}$ and $d\bar{d}$ rather than a pure state?

The Color Quantum Number and the $\Delta^{++}$ Problem

The quark model faced an immediate crisis. The $\Delta^{++}$ baryon has quark content $uuu$ and spin $J = 3/2$. In the ground state ($L = 0$), all three quarks occupy the same spatial state with all three spins aligned. But quarks are fermions — they should obey the Pauli exclusion principle. How can three identical fermions occupy the same quantum state?

The resolution, proposed by O.W. Greenberg (1964) and by Han and Nambu (1965), is that quarks carry an additional quantum number called color charge. Each quark comes in three color variants: red ($r$), green ($g$), and blue ($b$). (The names are whimsical — color charge has nothing to do with visible light.) Antiquarks carry anti-colors ($\bar{r}$, $\bar{g}$, $\bar{b}$). With this new quantum number, the three $u$ quarks in the $\Delta^{++}$ can be distinguished by their colors ($u_r u_g u_b$), and the total wavefunction is antisymmetric in the color indices (as required by the Pauli principle).

The color quantum number is not merely a bookkeeping device. It is the charge of the strong force — the analogue of electric charge in electromagnetism. This takes us from the quark model to the full dynamical theory.

31.2 Quantum Chromodynamics: The Theory of the Strong Interaction

QCD as a Gauge Theory

Quantum chromodynamics (QCD) is the gauge theory of the strong interaction. It is based on the symmetry group SU(3)$_{\text{color}}$ — the group of special unitary $3 \times 3$ matrices — acting on the color quantum number of quarks. The structure closely parallels quantum electrodynamics (QED), but with a crucial difference: where QED has one type of charge (electric) and one gauge boson (the photon), QCD has three types of charge (color) and eight gauge bosons (the gluons).

The QCD Lagrangian density can be written compactly:

$$\mathcal{L}_{\text{QCD}} = \sum_f \bar{\psi}_f^a \left(i \gamma^\mu D_\mu - m_f\right)_{ab} \psi_f^b - \frac{1}{4} G_{\mu\nu}^A G^{A\mu\nu}$$

where: - $\psi_f^a$ is the quark field for flavor $f$ and color $a = 1, 2, 3$ (red, green, blue) - $D_\mu = \partial_\mu - i g_s T^A A_\mu^A$ is the covariant derivative, with $g_s$ the strong coupling constant and $T^A$ ($A = 1, \ldots, 8$) the SU(3) generators (Gell-Mann matrices divided by 2) - $A_\mu^A$ are the eight gluon fields - $G_{\mu\nu}^A = \partial_\mu A_\nu^A - \partial_\nu A_\mu^A + g_s f^{ABC} A_\mu^B A_\nu^C$ is the gluon field strength tensor - $f^{ABC}$ are the SU(3) structure constants - The sum runs over all quark flavors $f = u, d, s, c, b, t$

This Lagrangian, occupying just one line, contains the complete theory of the strong interaction at the fundamental level. Every nuclear force, every hadron mass, every quark-gluon interaction is encoded here. The difficulty is not in writing the theory — it is in solving it.

Gluons: The Gauge Bosons of QCD

The eight gluons are massless spin-1 particles that mediate the strong force between quarks. Each gluon carries one unit of color and one unit of anti-color — for example, a gluon might carry the color combination $r\bar{g}$ (red, anti-green). This is the fundamental difference between QCD and QED: the photon is electrically neutral and does not interact with itself, but gluons carry color charge and therefore interact with other gluons.

This gluon self-interaction, encoded in the $g_s f^{ABC} A_\mu^B A_\nu^C$ term of the field strength tensor, is responsible for the two most remarkable properties of QCD: confinement and asymptotic freedom.

The Strong Coupling Constant

The dimensionless coupling constant of QCD, analogous to the fine structure constant $\alpha = e^2/(4\pi\epsilon_0 \hbar c) \approx 1/137$ of QED, is:

$$\alpha_s = \frac{g_s^2}{4\pi}$$

Unlike $\alpha$, which is small enough that perturbation theory works beautifully in QED, the value of $\alpha_s$ depends dramatically on the energy scale (momentum transfer) at which it is measured. This running of the coupling constant is the key to understanding both asymptotic freedom and confinement.

At high energies (short distances), $\alpha_s$ is small — perturbation theory works:

$$\alpha_s(M_Z) \approx 0.118 \quad (\text{at } Q = M_Z = 91.2 \text{ GeV})$$

$$\alpha_s(M_\tau) \approx 0.33 \quad (\text{at } Q = M_\tau = 1.78 \text{ GeV})$$

At nuclear energy scales ($Q \sim 0.2$--$1$ GeV), $\alpha_s$ becomes of order 1, and perturbation theory fails completely. The running of $\alpha_s$ is governed by the renormalization group equation:

$$\frac{d\alpha_s}{d\ln Q^2} = -\frac{\beta_0}{2\pi} \alpha_s^2 - \frac{\beta_1}{4\pi^2} \alpha_s^3 - \cdots$$

where the leading coefficient is:

$$\beta_0 = 11 - \frac{2}{3} n_f$$

with $n_f$ the number of active quark flavors. For $n_f \leq 16$ (and nature has $n_f = 6$), $\beta_0 > 0$, which means $\alpha_s$ decreases at high energies — this is asymptotic freedom. The "11" comes from gluon self-interactions; the "$-2n_f/3$" comes from quark loop screening (analogous to the vacuum polarization that makes $\alpha$ increase at high energies in QED). The gluon contribution dominates and reverses the sign — a qualitative difference from QED.

At leading order, the running coupling is:

$$\alpha_s(Q^2) = \frac{4\pi}{\beta_0 \ln(Q^2/\Lambda_{\text{QCD}}^2)}$$

where $\Lambda_{\text{QCD}} \approx 200$--$300$ MeV is the fundamental scale of the strong interaction. When $Q$ approaches $\Lambda_{\text{QCD}}$ from above, $\alpha_s$ diverges — perturbation theory ceases to be valid, and the physics of confinement takes over.

31.3 Confinement: Why We Never See Free Quarks

The Experimental Fact

Despite decades of searching — in cosmic rays, in accelerator collisions, in moon rocks, in seawater, in magnetic monopole detectors — no one has ever observed a free quark. Every experiment that has looked for fractional electric charge ($\pm 1/3 e$ or $\pm 2/3 e$) in isolation has come up empty (the one claimed exception, Fairbank's magnetic levitation experiment in the 1980s, was not confirmed and is now regarded as a systematic artifact).

This is confinement: quarks and gluons are permanently bound inside hadrons. Only color-neutral (color-singlet) combinations can exist as free particles. The allowed combinations are:

• Baryons: $q q q$ (one quark of each color, forming a color singlet $\epsilon_{abc} q^a q^b q^c$)
• Antibaryons: $\bar{q}\bar{q}\bar{q}$
• Mesons: $q\bar{q}$ (color-anticolor combination $\delta_{a\bar{a}} q^a \bar{q}^{\bar{a}}$)
• Exotic hadrons: tetraquarks ($qq\bar{q}\bar{q}$), pentaquarks ($qqqq\bar{q}$), glueballs ($gg$, $ggg$) — predicted by QCD and now experimentally confirmed (LHCb has discovered several candidates since 2015)

The Color String Picture

The physical mechanism of confinement can be understood through the color flux tube (or "string") picture. In QED, the electric field lines between an electron and a positron spread out in all directions, producing a $1/r$ potential. In QCD, the gluon self-interaction causes the color field lines between a quark and an antiquark to be squeezed into a narrow tube — a flux tube with an approximately constant energy per unit length, $\sigma \approx 0.9$ GeV/fm $\approx 14$ tonnes.

The resulting potential between a quark and an antiquark at large separation is:

$$V_{q\bar{q}}(r) \approx -\frac{4}{3} \frac{\alpha_s}{r} + \sigma r$$

The first term is the short-distance Coulomb-like color interaction (the factor $4/3$ is the color Casimir for the fundamental representation). The second term, linear in $r$, is the confining string. As you try to pull a quark-antiquark pair apart, the energy stored in the string increases linearly. At a separation of about 1 fm, the energy stored in the string ($\sim 0.9$ GeV) exceeds the threshold for creating a new $q\bar{q}$ pair from the vacuum. The string "breaks," producing two mesons rather than a free quark — you have created new hadrons instead of liberating a quark. This is why high-energy collisions produce jets of hadrons rather than free quarks.

Confinement and the Origin of Mass

This picture leads to one of the most profound facts in physics: the mass of visible matter is overwhelmingly generated by the strong force, not by the Higgs mechanism.

The $u$ and $d$ quark masses are approximately 2 and 5 MeV respectively. Three quarks contribute roughly $10$ MeV to the proton mass of $938$ MeV — barely $1\%$. The remaining $\sim 99\%$ comes from the kinetic energy of the quarks (confined in a small volume, the uncertainty principle demands large momenta) and the energy stored in the gluon field. This is $E = mc^2$ in reverse: the mass of the proton is overwhelmingly generated by the energy of the strong force.

Lattice QCD calculations (Section 31.5) have confirmed this picture quantitatively. In the limit where the $u$ and $d$ quark masses are set to zero, the proton mass barely changes — it drops from $938$ MeV to roughly $860$ MeV. The strong force generates mass from pure energy.

💡 Key Insight: The Higgs mechanism gives mass to the quarks and leptons, but it accounts for less than 2% of the mass of ordinary matter. The strong interaction — through confinement and the energy of the gluon field — generates the overwhelming majority of the mass of protons, neutrons, and therefore of you, this textbook, and the Earth.

31.4 Asymptotic Freedom: Why Quarks Behave Freely at Short Distances

The Discovery

In 1973, David Gross, Frank Wilczek, and H. David Politzer independently discovered that non-abelian gauge theories (like QCD) have the remarkable property of asymptotic freedom: the coupling constant decreases at short distances (high momentum transfers). They received the Nobel Prize in Physics in 2004.

Asymptotic freedom means that at very high energies ($Q \gg \Lambda_{\text{QCD}}$), quarks interact weakly and behave almost like free particles. This is why deep inelastic scattering experiments at SLAC in the late 1960s — which probed the proton's interior with high-energy electrons — observed nearly free point-like constituents (called "partons" by Feynman before they were identified with quarks). The partons of Feynman and the quarks of Gell-Mann are the same objects, and asymptotic freedom explains why they appear free when probed at short wavelengths.

The physical picture is the opposite of what happens in QED. In QED, the vacuum around a bare electron is filled with virtual electron-positron pairs that screen the charge — the effective charge appears smaller at large distances and larger at short distances. This is vacuum polarization, and it causes the fine structure constant $\alpha$ to increase (very slowly) at higher energies: $\alpha(M_Z) \approx 1/128$ versus $\alpha(0) \approx 1/137$.

In QCD, the gluon self-interaction introduces an anti-screening effect that overwhelms the quark screening. Virtual gluon loops spread the color charge, making it appear larger at large distances and smaller at short distances. The net effect — anti-screening dominating over screening when $n_f < 33/2$ — gives asymptotic freedom.

Why This Matters for Nuclear Physics

Asymptotic freedom has two critical consequences for the program of this textbook:

1. Perturbative QCD works at high energies. At momentum transfers $Q \gtrsim 2$ GeV ($\alpha_s \lesssim 0.3$), perturbative calculations in QCD are reliable. This is the regime of deep inelastic scattering, jet production, and hard processes in heavy-ion collisions. We can calculate with controlled uncertainty.

2. Perturbative QCD fails at nuclear energy scales. At the energy scales relevant to nuclear binding ($Q \sim 200$--$500$ MeV, $\alpha_s \sim 1$), the coupling is strong and perturbation theory is useless. This is the regime of confinement, hadron formation, and the nuclear force. We cannot calculate the nuclear force directly from the QCD Lagrangian using the standard Feynman diagram approach. We need non-perturbative methods.

This is why nuclear physics, despite being "derived" from a known fundamental theory, is so difficult. The nuclear force lives in the worst possible regime of QCD — too strongly coupled for perturbation theory, too complicated for exact analytic solution. The two bridges across this gap are lattice QCD (Section 31.5) and effective field theory (Section 31.6).

🔄 Check Your Understanding: Why does the gluon self-interaction reverse the sign of the beta function compared to QED? What would happen to QCD if there were 20 quark flavors instead of 6?

31.5 From QCD to Nuclear Forces: The Emergence Problem

The Central Challenge

We now arrive at the central question of this chapter: how does the nuclear force — an interaction between color-neutral nucleons separated by $\sim 1$ fm — emerge from QCD, a theory of colored quarks and gluons interacting at distances $\lesssim 0.1$ fm?

This is an emergence problem. The nuclear force is not "contained" in the QCD Lagrangian in any obvious way, just as the van der Waals attraction between argon atoms is not apparent from the Coulomb's law between charged particles. In both cases, a residual interaction between composite, internally neutral objects arises from the underlying fundamental force between charged constituents.

The analogy to the van der Waals force is precise and illuminating:

But there is a crucial difference in difficulty. The van der Waals force can be calculated perturbatively from QED, because the fine structure constant is small ($\alpha \approx 1/137$). The nuclear force cannot be calculated perturbatively from QCD, because $\alpha_s \sim 1$ at the relevant scale. We need non-perturbative tools.

Lattice QCD: First Principles on a Supercomputer

Lattice QCD is the only known method for solving QCD from first principles in the non-perturbative regime. The idea, proposed by Kenneth Wilson in 1974, is conceptually simple: discretize spacetime on a four-dimensional Euclidean lattice with spacing $a$, replace the continuous quark and gluon fields with fields defined on the lattice sites and links, and compute path integrals numerically using Monte Carlo methods.

The lattice QCD partition function is:

$$Z = \int \mathcal{D}U \, \mathcal{D}\bar{\psi} \, \mathcal{D}\psi \; \exp\left(-S_G[U] - S_F[\bar{\psi}, \psi, U]\right)$$

where $S_G$ is the gauge (gluon) action and $S_F$ is the fermion (quark) action, both discretized on the lattice. The gauge fields live on the links between lattice sites (as SU(3) matrices $U_\mu(x)$), and the quark fields live on the sites. The path integral is evaluated stochastically — generating an ensemble of gauge field configurations weighted by $e^{-S}$, then averaging observables over this ensemble.

The computational cost is enormous. The key challenges are:

1. Lattice spacing $a$ must be small compared to the physical scales of interest ($a \lesssim 0.05$--$0.1$ fm), requiring large lattices.
2. The lattice volume $L^3 \times T$ must be large compared to the hadron size ($L \gtrsim 3$--$5$ fm), requiring $L/a \gtrsim 30$--$100$ lattice sites per dimension.
3. Light quark masses are computationally expensive. The cost of fermion calculations scales roughly as $1/m_q^2$. Early lattice calculations used unphysically heavy quarks ($m_\pi \sim 300$--$500$ MeV) and extrapolated; modern calculations work directly at the physical pion mass ($m_\pi = 135$--$140$ MeV).
4. Multi-nucleon systems require special techniques (the signal-to-noise problem grows exponentially with baryon number).

Despite these challenges, lattice QCD has achieved remarkable successes.

Hadron Masses from First Principles

The most celebrated result of lattice QCD is the ab initio calculation of the hadron mass spectrum. The BMW (Budapest-Marseille-Wuppertal) collaboration published a landmark calculation in 2008 (updated subsequently) in which the masses of the lightest hadrons — the pion, kaon, proton, neutron, $\Delta$, $\Sigma$, $\Xi$, $\Omega$, and others — were computed directly from the QCD Lagrangian with physical quark masses, at the percent level of accuracy.

The inputs to such a calculation are: the QCD Lagrangian, the quark masses (or equivalently, $m_\pi$, $m_K$, and one baryon mass to set the scale), and $\alpha_s$. Everything else — the entire hadron spectrum — is a prediction. The agreement with experiment is spectacular:

The proton-neutron mass difference — a tiny 1.3 MeV signal on a 940 MeV background — has been computed by the BMW collaboration including electromagnetic and $m_u - m_d$ effects, obtaining $m_n - m_p = 1.51 \pm 0.28$ MeV, consistent with experiment.

This is a triumph: QCD, with just a handful of input parameters, correctly predicts the masses of the particles that make up the visible universe.

Nuclear Forces from the Lattice

Computing the nuclear force from lattice QCD is far more difficult than computing hadron masses, because it requires extracting the interaction energy between two (or more) baryons — a signal that is exponentially smaller than the statistical noise for large baryon number.

The pioneering work of the HAL QCD (Hadrons to Atomic nuclei from Lattice QCD) collaboration in Japan has developed methods to extract nuclear potentials from lattice QCD by studying the Nambu-Bethe-Salpeter wave function of two-baryon systems. Their results reproduce the qualitative features of the nuclear force:

• Attraction at intermediate range ($\sim 1$ fm)
• Repulsion at short range ($\lesssim 0.5$ fm)
• Tensor force from one-pion exchange at long range

The NPLQCD (Nuclear Physics from Lattice QCD) collaboration has computed the binding energy of the deuteron and light nuclei, though at unphysically heavy pion masses ($m_\pi \sim 300$--$800$ MeV). At the physical pion mass, the binding energies of light nuclei remain extremely challenging.

A milestone was achieved in 2013 when lattice QCD calculations provided evidence for the existence (or near-existence) of the $H$-dibaryon — a hypothetical six-quark state ($uuddss$) predicted by Jaffe in 1977. The lattice results suggest a near-threshold state, but its fate at the physical quark masses remains uncertain.

The current frontier is computing nuclear binding energies with controlled uncertainties at physical quark masses. While we are not yet able to derive the nuclear chart from QCD, the trajectory is clear — lattice QCD is progressing toward this goal, and it validates the qualitative picture of the nuclear force as a residual strong interaction.

📜 Historical Context: Wilson's Vision

Kenneth Wilson proposed lattice gauge theory in 1974 and received the Nobel Prize in 1982 for his work on the renormalization group and critical phenomena. Wilson understood that the lattice was not merely a numerical tool but a conceptual framework: it showed that QCD could confine quarks, even if an analytical proof remained elusive (and remains so — proving confinement rigorously from the QCD Lagrangian is one of the seven Millennium Prize problems in mathematics, with a $1 million prize). The exponential growth of computing power has turned Wilson's conceptual tool into a precision instrument.

31.6 The Bridge: Chiral Effective Field Theory

Why Effective Field Theory?

Lattice QCD is the direct approach — solve QCD numerically. But it has severe limitations for nuclear physics: multi-nucleon systems with $A > 4$ are extremely expensive, nuclei at finite temperature or density are essentially inaccessible, and the physical insight that comes from analytic expressions is absent. Nuclear physicists need a framework that captures the essential physics of QCD at low energies without solving the full theory.

This is the role of effective field theory (EFT): a systematic framework for describing physics at a given energy scale by retaining only the degrees of freedom relevant at that scale, while encoding the effects of higher-energy physics in a set of coupling constants.

The EFT relevant for nuclear forces is chiral effective field theory (chiral EFT), which we introduced briefly in Chapter 3. Here we develop the conceptual foundations more fully.

Chiral Symmetry of QCD

In the limit where the $u$ and $d$ quark masses are set to zero ($m_u = m_d = 0$), the QCD Lagrangian has an additional symmetry called chiral symmetry: independent rotations of the left-handed and right-handed quark fields. This is the symmetry group SU(2)$_L \times$ SU(2)$_R$ (for two flavors).

If chiral symmetry were an exact symmetry of the ground state (the QCD vacuum), then for every hadron with a given parity, there would be a partner of equal mass and opposite parity. But the $\rho$ meson ($J^\pi = 1^-$, $m = 775$ MeV) and the $a_1$ meson ($J^\pi = 1^+$, $m = 1230$ MeV) are not degenerate. Chiral symmetry is spontaneously broken by the QCD vacuum: the ground state does not respect the symmetry of the Lagrangian.

Goldstone's theorem then requires the existence of massless spin-0 particles — one for each broken generator. For SU(2)$_L \times$ SU(2)$_R \to$ SU(2)$_V$ (isospin), there are three broken generators, giving three Goldstone bosons. These are the pions ($\pi^+$, $\pi^0$, $\pi^-$).

In reality, the quark masses are small but not zero: $m_u \approx 2$ MeV, $m_d \approx 5$ MeV. This gives the pions a small mass:

$$m_\pi^2 \approx \frac{(m_u + m_d) \langle \bar{q}q \rangle}{f_\pi^2}$$

where $\langle \bar{q}q \rangle \approx -(250 \text{ MeV})^3$ is the quark condensate (the order parameter of chiral symmetry breaking) and $f_\pi \approx 92$ MeV is the pion decay constant. The pions are pseudo-Goldstone bosons — light because the symmetry-breaking quark masses are small, but not exactly massless.

This explains one of the deepest facts in nuclear physics: the pion is light because the quarks are light. And because the pion is light, it mediates the longest-range component of the nuclear force.

Chiral Perturbation Theory and Nuclear Forces

Chiral perturbation theory (ChPT), developed by Weinberg, Gasser, and Leutwyler in the 1970s and 1980s, provides a systematic expansion of low-energy QCD observables in powers of:

$$\frac{Q}{\Lambda_\chi} \quad \text{and} \quad \frac{m_\pi}{\Lambda_\chi}$$

where $Q$ is a typical momentum, $m_\pi$ is the pion mass, and $\Lambda_\chi \sim 4\pi f_\pi \sim 1$ GeV is the chiral symmetry breaking scale.

Weinberg's crucial insight (1990--1991) was to apply chiral EFT to nuclear forces. In this framework, the nuclear potential is organized as a systematic expansion in powers of $(Q/\Lambda_\chi)^\nu$, where the chiral order $\nu$ determines the importance of each contribution:

Leading order (LO), $\nu = 0$: - One-pion exchange (OPEP) — the longest-range part of the nuclear force - Two contact terms (short-range nucleon-nucleon interactions)

Next-to-leading order (NLO), $\nu = 2$: - Two-pion exchange with leading vertex corrections - Seven additional contact terms

Next-to-next-to-leading order (N$^2$LO), $\nu = 3$: - Two-pion exchange with subleading vertices - First appearance of the three-nucleon force (the Fujita-Miyazawa mechanism and its chiral partners) - No new $NN$ contact terms at this order

N$^3$LO, $\nu = 4$: - Three-loop two-pion exchange - Leading four-nucleon force (first appearance) - 15 additional $NN$ contact terms - This is the current state of the art for precision nuclear forces

The power counting provides a hierarchy: each successive order is suppressed by $(Q/\Lambda_\chi)$ relative to the previous one. For typical nuclear momenta $Q \sim m_\pi \approx 140$ MeV and $\Lambda_\chi \sim 1$ GeV, the expansion parameter is $Q/\Lambda_\chi \sim 0.14$, suggesting reasonable convergence.

Key Achievements of Chiral EFT

The chiral EFT program, carried out by groups worldwide (Bochum-Bonn, the LENPIC collaboration, the Oak Ridge group, and others), has produced nuclear forces that:

1. Reproduce $NN$ scattering data with $\chi^2/\text{datum} \approx 1$ at N$^3$LO, rivaling phenomenological potentials like Argonne $v_{18}$.

2. Predict three-nucleon forces consistently. Unlike older approaches (Urbana, Illinois models) where the 3NF was an ad hoc addition, in chiral EFT the two-body and three-body forces emerge from the same Lagrangian at a definite order. The 3NF at N$^2$LO depends on only two new parameters ($c_D$ and $c_E$), determined from $A = 3$ observables.

3. Provide systematic uncertainty estimates. Because the expansion is organized in powers of $Q/\Lambda_\chi$, the truncation error at each order can be estimated — a major advantage over phenomenological models.

4. Connect to QCD through lattice calculations. The low-energy constants of chiral EFT can, in principle, be computed from lattice QCD, closing the loop from fundamental theory to nuclear structure.

5. Successfully describe light nuclei. Combined with ab initio many-body methods (quantum Monte Carlo, no-core shell model, coupled cluster), chiral EFT nuclear forces reproduce the ground-state energies and spectra of nuclei up to $A \sim 100$ with controlled uncertainties.

The chiral EFT framework is the modern realization of the program hinted at in Chapter 3: a principled, systematically improvable bridge between QCD and nuclear physics.

🔄 Check Your Understanding: Why do three-nucleon forces first appear at N$^2$LO in the chiral expansion, not at leading order? What determines the chiral order of a given Feynman diagram contribution to the nuclear force?

31.7 The Nucleon: Internal Structure

Having established that nucleons are composite objects, we now turn to what experiments have revealed about their internal structure. Three classes of observables provide complementary windows into the nucleon: form factors (the spatial distribution of charge and magnetization), parton distribution functions (the momentum distribution of quarks and gluons), and spin structure.

Electromagnetic Form Factors

The most direct probe of nucleon structure is elastic electron-nucleon scattering. Because the electron interacts via the well-understood electromagnetic force, it serves as a clean, calibrated probe of the nucleon's charge and magnetization distributions.

For a point-like proton, the elastic scattering cross section would be the Mott cross section (the relativistic generalization of Rutherford scattering):

$$\left(\frac{d\sigma}{d\Omega}\right)_{\text{Mott}} = \frac{\alpha^2}{4E^2\sin^4(\theta/2)} \cdot \frac{E'}{E} \cos^2\frac{\theta}{2}$$

where $E$ and $E'$ are the initial and final electron energies. For a nucleon with internal structure, the cross section is modified by two form factors:

$$\frac{d\sigma}{d\Omega} = \left(\frac{d\sigma}{d\Omega}\right)_{\text{Mott}} \cdot \left[ \frac{G_E^2(Q^2) + \tau G_M^2(Q^2)}{1 + \tau} + 2\tau G_M^2(Q^2) \tan^2\frac{\theta}{2} \right]$$

This is the Rosenbluth formula, where $Q^2 = -q^2$ is the four-momentum transfer squared, $\tau = Q^2/(4M_N^2)$, and $G_E(Q^2)$ and $G_M(Q^2)$ are the Sachs electric and magnetic form factors of the nucleon.

At $Q^2 = 0$:

$$G_E^p(0) = 1, \quad G_E^n(0) = 0$$ $$G_M^p(0) = \mu_p = 2.793 \; \mu_N, \quad G_M^n(0) = \mu_n = -1.913 \; \mu_N$$

These are simply the total charge and magnetic moment.

In the non-relativistic limit, the Fourier transform of $G_E(Q^2)$ gives the spatial distribution of charge:

$$\rho_E(r) = \frac{1}{(2\pi)^3} \int G_E(\mathbf{q}^2) \, e^{i\mathbf{q}\cdot\mathbf{r}} \, d^3q$$

and the RMS charge radius is related to the slope of $G_E$ at $Q^2 = 0$:

$$\langle r^2 \rangle_E = -6 \frac{dG_E}{dQ^2}\bigg|_{Q^2=0}$$

Robert Hofstadter pioneered electron scattering from nuclei at Stanford in the 1950s, revealing for the first time that the proton has a finite size — a charge radius of approximately 0.87 fm. He received the Nobel Prize in 1961 "for his pioneering studies of electron scattering in atomic nuclei and for his thereby achieved discoveries concerning the structure of the nucleons."

The Dipole Approximation

A striking experimental finding is that all four nucleon form factors are approximately described by a single dipole function:

$$G_D(Q^2) = \left(1 + \frac{Q^2}{0.71 \; \text{GeV}^2}\right)^{-2}$$

such that:

$$G_E^p(Q^2) \approx G_D(Q^2), \quad \frac{G_M^p(Q^2)}{\mu_p} \approx \frac{G_M^n(Q^2)}{\mu_n} \approx G_D(Q^2)$$

The Fourier transform of the dipole form factor corresponds to an exponential charge distribution $\rho(r) \propto e^{-r/a}$ with $a \approx 0.23$ fm, giving an RMS charge radius $r_p \approx 0.81$ fm. More precise measurements give $r_p \approx 0.84$--$0.88$ fm (the range reflects the proton radius puzzle, discussed in Section 31.9).

The neutron's electric form factor $G_E^n(Q^2)$ is not zero for $Q^2 > 0$, even though the neutron is electrically neutral ($G_E^n(0) = 0$). This means the neutron has a non-trivial charge distribution — positive charge near the center (from the $u$ quark core) and negative charge in the outer region (from the $d$ quark cloud and the pion cloud). The resulting charge radius squared is small and negative: $\langle r^2 \rangle_n = -0.1161 \pm 0.0022$ fm$^2$.

📜 Historical Context: Hofstadter's Revolution

Before Hofstadter's experiments, it was widely assumed that the proton was a point particle, like the electron. His demonstration that the proton has a finite size and a smooth charge distribution was one of the first direct pieces of evidence for substructure — a harbinger of the quark model that would arrive a decade later. Hofstadter's experiments used the Mark III linear accelerator at Stanford, which produced electron beams up to 550 MeV. Modern electron scattering facilities (Jefferson Lab, Mainz MAMI) reach much higher energies and have mapped the form factors with exquisite precision, revealing subtle deviations from the dipole approximation.

Parton Distribution Functions

At high momentum transfer ($Q^2 \gg 1$ GeV$^2$), the probe resolves individual quarks and gluons inside the nucleon. Deep inelastic scattering (DIS) experiments — in which the electron transfers enough energy to shatter the nucleon — provide information about the parton distribution functions (PDFs) $f_i(x, Q^2)$, which give the probability of finding a parton of type $i$ (quark, antiquark, or gluon) carrying a fraction $x$ of the nucleon's momentum, when probed at scale $Q^2$.

The key experimental observable is the structure function:

$$F_2(x, Q^2) = \sum_i e_i^2 \, x \, [f_i(x, Q^2) + \bar{f}_i(x, Q^2)]$$

where $e_i$ is the electric charge of quark flavor $i$ (in units of $e$) and the sum runs over all active quark flavors.

The PDFs satisfy a sum rule (momentum conservation):

$$\int_0^1 x \left[\sum_i (f_i(x) + \bar{f}_i(x)) + f_g(x)\right] dx = 1$$

where $f_g(x)$ is the gluon distribution. Experiments show that quarks carry only about 45--50% of the proton's momentum — the remaining 50--55% is carried by gluons. This was one of the early confirmations of the gluon as a physical, momentum-carrying constituent.

For nuclear physics, the EMC effect is a particularly important DIS result. In 1983, the European Muon Collaboration (EMC) at CERN discovered that the structure functions of nucleons bound in nuclei are modified compared to free nucleons:

$$\frac{F_2^A(x)}{A \cdot F_2^N(x)} \neq 1$$

The ratio shows characteristic deviations: shadowing at small $x$ (ratio $< 1$), anti-shadowing at $x \sim 0.1$ (ratio $> 1$), the EMC effect at $x \sim 0.3$--$0.7$ (ratio $< 1$), and Fermi motion at $x > 0.7$ (ratio $> 1$). The EMC effect in the intermediate-$x$ region — a suppression of about 10--15% for heavy nuclei — was unexpected and remains incompletely understood after four decades. It shows that the quark-gluon structure of a nucleon is modified by the nuclear medium, challenging the simple picture of nuclei as collections of independent nucleons.

31.8 The Proton Spin Puzzle

The Crisis

The proton has spin $1/2$. In the naive quark model, this spin comes from the quark spins: two quarks with spin up and one with spin down, giving a net $1/2$. If this picture were correct, the fraction of the proton's spin carried by the quarks would be:

$$\Delta \Sigma = \Delta u + \Delta d + \Delta s = 1$$

where $\Delta q$ is the fraction of the proton spin carried by quarks of flavor $q$.

In 1988, the European Muon Collaboration (EMC) at CERN measured the spin-dependent structure function of the proton using deep inelastic scattering of polarized muons from polarized proton targets. Their result was shocking:

$$\Delta \Sigma = 0.12 \pm 0.09 \pm 0.14$$

The quarks carry only about 12% of the proton's spin — essentially zero within uncertainties. This was the "proton spin crisis."

The Resolution: A Decomposition

Subsequent experiments (COMPASS at CERN, HERMES at DESY, STAR and PHENIX at RHIC, and Jefferson Lab) have refined the picture. The modern decomposition of the proton spin is:

$$\frac{1}{2} = \frac{1}{2}\Delta\Sigma + \Delta G + L_q + L_g$$

where: - $\frac{1}{2}\Delta\Sigma \approx 0.15$--$0.20$ — quark spin contributes about 30--40% (the modern value is larger than the original EMC result, though still far less than 100%) - $\Delta G$ — the gluon spin contribution. RHIC measurements indicate $\Delta G \approx 0.2$--$0.3$ for gluons carrying momentum fraction $x > 0.05$. This accounts for roughly 40--60% when extrapolated to all $x$. - $L_q + L_g$ — orbital angular momentum of quarks and gluons. This is the least well determined component. Generalized parton distributions (GPDs), accessible through deeply virtual Compton scattering (DVCS) and deeply virtual meson production, provide experimental access to orbital angular momentum.

The proton spin puzzle is not fully resolved, but the picture is now far more nuanced than the naive quark model. The proton's spin arises from a complex interplay of quark spin, gluon spin, and orbital motion — a result that would have been impossible to anticipate without experimental guidance.

The Electron-Ion Collider (EIC), currently under construction at Brookhaven National Laboratory with an expected start date around 2032, will make definitive measurements of the gluon spin contribution and the orbital angular momentum terms, closing this decades-old puzzle.

💡 Key Insight: The proton spin puzzle illustrates a recurring theme: the internal structure of the nucleon is far richer and more complex than the simple $qqq$ picture suggests. Quarks, antiquarks, gluons, and their orbital motion all contribute to the nucleon's quantum numbers. The "valence quark model" is a useful classification tool but a poor description of the nucleon's actual wavefunction.

31.9 The Proton Radius Puzzle

The Puzzle

The proton charge radius $r_p$ has been measured by two independent methods:

1. Electron scattering. Fitting the $Q^2$ dependence of $G_E^p(Q^2)$ yields the radius. The CODATA 2010 world average from this method (combined with hydrogen spectroscopy) was:

$$r_p(\text{electron}) = 0.8775 \pm 0.0051 \text{ fm}$$

2. Muonic hydrogen spectroscopy. In 2010, the CREMA collaboration at PSI (Switzerland) measured the Lamb shift in muonic hydrogen — a hydrogen atom where the electron is replaced by a muon. Because the muon is 207 times heavier than the electron, its Bohr radius is 207 times smaller, and it spends far more time inside the proton. This dramatically enhances the sensitivity to the proton radius. Their result was:

$$r_p(\text{muon}) = 0.84087 \pm 0.00039 \text{ fm}$$

The discrepancy — $0.84$ fm vs. $0.88$ fm — is about $5\sigma$ in statistical significance. This became known as the proton radius puzzle.

The Stakes

The discrepancy was taken extremely seriously because it could, in principle, signal:

• New physics: a previously unknown interaction that distinguishes electrons from muons (violating lepton universality)
• A problem with QED calculations: an error in the higher-order QED corrections used to extract the radius from spectroscopic data
• Experimental systematic errors: in the electron scattering or hydrogen spectroscopy measurements

Toward Resolution

The situation has evolved significantly since 2010. Several developments point toward the smaller, muonic hydrogen value being correct:

1. The PRad experiment at Jefferson Lab (2019) measured the proton charge radius via electron-proton scattering with a novel calorimetric technique, obtaining $r_p = 0.831 \pm 0.007_{\text{stat}} \pm 0.012_{\text{syst}}$ fm — consistent with the muonic hydrogen value.

2. Hydrogen spectroscopy measurements at the Max Planck Institute for Quantum Optics (Garching), using the $2S$--$4P$ transition, obtained $r_p = 0.8335 \pm 0.0095$ fm (2017), again consistent with the muonic value.

3. The York/Toronto group re-analyzed the world electron scattering data and found that a more careful treatment of the low-$Q^2$ extrapolation yields a smaller radius, compatible with the muonic result.

4. The CODATA 2018 recommended value shifted to $r_p = 0.8414 \pm 0.0019$ fm, significantly closer to the muonic hydrogen result.

The emerging consensus (as of the mid-2020s) is that the puzzle is largely resolved in favor of the smaller value, $r_p \approx 0.84$ fm. The earlier larger values from electron scattering appear to have been affected by systematic issues in the fitting of form factor data at low $Q^2$. However, some tension remains, and the ongoing experimental programs at Jefferson Lab, Mainz, and elsewhere continue to refine the measurements.

The proton radius puzzle, even if it turns out not to require new physics, has been enormously productive. It has driven improvements in:

• Precision QED calculations for atomic physics
• Experimental techniques for electron scattering at very low $Q^2$
• Understanding of two-photon exchange corrections
• Development of muonic atom spectroscopy as a precision tool

🔄 Check Your Understanding: Why is muonic hydrogen more sensitive to the proton charge radius than ordinary (electronic) hydrogen? What is the ratio of the muon's Bohr radius to the electron's, and what fraction of its time does the muon spend inside the proton?

31.10 Connecting Back: From Quarks to Nuclei

We can now give a complete, multi-scale answer to the question posed in Chapter 3: what holds the nucleus together?

Scale 1: Quarks and gluons ($r < 0.1$ fm). At the shortest distances, QCD is the fundamental theory. Quarks interact via gluon exchange, with a coupling $\alpha_s$ that is small at short distances (asymptotic freedom) but large at distances approaching 1 fm (confinement).

Scale 2: Nucleon formation ($r \sim 0.5$--$1$ fm). The strong force confines quarks and gluons into color-neutral hadrons. The proton ($uud$) and neutron ($udd$) emerge as the lightest baryons. Their mass ($\sim 939$ MeV) is overwhelmingly generated by the energy of the confined gluon field, not by the quark masses ($\sim 10$ MeV total). The spontaneous breaking of chiral symmetry produces the pion as a light pseudo-Goldstone boson.

Scale 3: The nuclear force ($r \sim 1$--$2$ fm). At distances larger than the nucleon radius, color is screened and the relevant degrees of freedom are nucleons and pions — not quarks and gluons. The nuclear force emerges as a residual effect of QCD: - At long range ($r > 2$ fm): one-pion exchange, determined by chiral symmetry - At intermediate range ($1 < r < 2$ fm): two-pion exchange, correlated meson exchange - At short range ($r < 1$ fm): overlap of nucleon quark distributions, encoded in contact terms

Chiral EFT provides the systematic framework connecting QCD's symmetries to the nuclear force. The expansion in powers of $Q/\Lambda_\chi$ organizes the nuclear force, including three-nucleon and higher many-body forces.

Scale 4: Nuclear structure ($r \sim 2$--$10$ fm). The nuclear force, combined with many-body quantum mechanics, produces the shell model, collective motion, nuclear binding, and the chart of nuclides. The nuclear many-body problem is solved using methods developed in Chapters 6--8.

Scale 5: Neutron stars ($r \sim 10$ km). At the macroscopic scale, nuclear forces (extrapolated to extreme densities) determine the equation of state and the structure of neutron stars, as discussed in Chapter 25.

This hierarchy — from quarks to nuclei to neutron stars, spanning more than 19 orders of magnitude in distance — is one of the most complete examples of emergence in all of physics. Each level has its own effective degrees of freedom and its own effective theory, yet the levels are connected by principled theoretical frameworks: QCD $\to$ chiral EFT $\to$ nuclear forces $\to$ nuclear structure $\to$ the equation of state.

31.11 The Quark-Gluon Plasma: Nuclear Matter at Extreme Temperatures

Before concluding, we mention briefly a state of matter where the distinction between nuclear physics and QCD disappears entirely. At temperatures above approximately $T_c \approx 155$ MeV ($\sim 1.8 \times 10^{12}$ K) — as computed by lattice QCD and confirmed experimentally — nuclear matter undergoes a phase transition to a quark-gluon plasma (QGP): a state in which quarks and gluons are deconfined and move freely over distances much larger than the nucleon size.

The QGP existed in the early universe during the first $\sim 10$ microseconds after the Big Bang. It has been recreated in the laboratory at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory and at the Large Hadron Collider (LHC) at CERN, by colliding heavy nuclei (typically gold or lead) at ultrarelativistic energies.

The surprising discovery at RHIC (announced in 2005) was that the QGP is not the weakly interacting gas of quarks and gluons that perturbative QCD had suggested, but rather a strongly coupled liquid — the "most perfect liquid" ever created, with a ratio of viscosity to entropy density near the conjectured lower bound from string theory ($\eta/s \geq \hbar/(4\pi k_B)$).

Experimental signatures of the QGP include:

• Jet quenching: high-energy quarks and gluons lose energy traversing the QGP, producing suppressed jet yields
• $J/\psi$ suppression: the $c\bar{c}$ bound state ($J/\psi$) is disrupted by color screening in the QGP
• Collective flow: the QGP exhibits elliptic and higher-order flow patterns, described by relativistic hydrodynamics
• Strangeness enhancement: strange quarks are produced more copiously in the QGP than in hadronic collisions

The study of the QGP connects nuclear physics directly to QCD and to the earliest moments of the universe. It is a reminder that the boundary between nuclear physics and particle physics is not a wall but a gradient.

Chapter Summary

• The quark model organizes all known hadrons as composite states of quarks ($qqq$ = baryons, $q\bar{q}$ = mesons). The proton is $uud$ and the neutron is $udd$. The up and down quark masses ($\sim 2$--$5$ MeV) account for less than 2% of the nucleon mass.

• Quantum chromodynamics (QCD) is the gauge theory of the strong interaction, based on SU(3) color symmetry. Quarks carry color charge and interact by exchanging eight massless gluons. Unlike the photon, gluons carry color and interact with each other.

• Confinement means quarks are never observed as free particles — only color-neutral combinations (hadrons) exist. The confining color flux tube stores about 0.9 GeV/fm. Over 98% of the proton's mass comes from the energy of the strong force, not from quark masses.

• Asymptotic freedom means the strong coupling $\alpha_s$ decreases at short distances (high energies). This allows perturbative QCD at high energies but makes low-energy (nuclear) QCD non-perturbative and analytically unsolvable.

• The nuclear force is a residual effect of QCD, analogous to the van der Waals force being a residual effect of electromagnetism. At nuclear distances ($\sim 1$ fm), QCD is non-perturbative, and the nuclear force emerges from the complex dynamics of colored quarks and gluons inside nucleons.

• Lattice QCD discretizes spacetime and computes QCD observables numerically. It has reproduced the hadron mass spectrum at the percent level and is making progress on nuclear forces, though multi-nucleon systems remain extremely challenging.

• Chiral EFT provides the systematic bridge between QCD symmetries and nuclear forces. The nuclear potential is organized as an expansion in $Q/\Lambda_\chi$, with two-nucleon, three-nucleon, and higher-body forces appearing at definite orders. Three-nucleon forces first appear at N$^2$LO.

• Nucleon form factors, measured by electron scattering, reveal the charge and magnetization distributions of the proton and neutron. The proton charge radius is $r_p \approx 0.84$ fm.

• The proton spin puzzle: quarks carry only $\sim 30$--$40\%$ of the proton spin. The remainder comes from gluon spin ($\sim 40$--$60\%$) and orbital angular momentum. The Electron-Ion Collider will make definitive measurements.

• The proton radius puzzle (muonic vs. electronic hydrogen) is largely resolved in favor of the smaller value ($r_p \approx 0.84$ fm), but continues to drive precision measurements and theoretical developments.

Next: Chapter 32 — Fundamental Symmetries Tested with Nuclei: Parity, CP, and Beyond