Case Study 1: From QCD to Nuclear Forces — The Most Complex Emergence Problem in Physics

The Problem

We know the fundamental theory of the strong interaction — quantum chromodynamics — and we can write its Lagrangian in a single line. We also know that nuclear forces bind protons and neutrons into the roughly 3,300 nuclides on the chart of nuclides, with properties measured to extraordinary precision. The problem is connecting the two: deriving the nuclear force, quantitatively and with controlled uncertainties, from the QCD Lagrangian.

This is arguably the most complex emergence problem in physics. The electromagnetic analogue — deriving the van der Waals force between neutral atoms from QED — can be done perturbatively on the back of an envelope. The gravitational analogue — deriving the Newtonian gravitational force from general relativity — is a straightforward weak-field limit. But deriving the nuclear force from QCD requires navigating a strongly coupled regime where no analytical methods work, and the hierarchy of scales (from quarks at $\sim 10^{-16}$ m to nuclei at $\sim 10^{-14}$ m) demands a chain of effective theories.

The Multi-Scale Chain

Scale 1: The QCD Lagrangian ($E > 2$ GeV, $r < 0.1$ fm)

At the highest energies, QCD is perturbative. Quarks and gluons are the explicit degrees of freedom. Deep inelastic scattering experiments have verified QCD predictions at the percent level in this regime. The PDFs of the nucleon are well-determined. This is the "easy" part of QCD.

Scale 2: The Confinement Scale ($E \sim 0.2$--$1$ GeV, $r \sim 0.2$--$1$ fm)

As the energy decreases below about 1--2 GeV, the coupling $\alpha_s$ grows and perturbation theory fails. This is where quarks are confined into hadrons — protons, neutrons, pions, and the other members of the hadron zoo. Lattice QCD operates in this regime, discretizing spacetime and computing path integrals numerically.

Milestone achievement: The BMW collaboration's calculation of the light hadron spectrum (2008 and subsequent updates). Starting from the QCD Lagrangian with physical quark masses, they computed the masses of the proton, neutron, and other light hadrons to within 1--3% of experimental values. No parameters were adjusted after fixing the quark masses and the overall scale — the hadron masses are genuine predictions.

Remaining challenge: Computing multi-baryon systems. The signal-to-noise ratio in lattice QCD degrades exponentially with baryon number: $\text{S/N} \propto e^{-A(m_N - 3m_\pi/2) T}$, where $A$ is the baryon number and $T$ is the Euclidean time separation. For the deuteron ($A = 2$), the noise grows roughly as $e^{-0.35T}$ per fm of temporal extent, making precise energy determinations extremely expensive.

Scale 3: The Chiral EFT Bridge ($E \sim 50$--$500$ MeV, $r \sim 0.4$--$4$ fm)

Chiral effective field theory provides the bridge between QCD and nuclear forces. The key insight is that at energies well below $\Lambda_\chi \sim 1$ GeV, the relevant degrees of freedom are no longer quarks and gluons but nucleons and pions. QCD constrains this effective theory through its symmetries — most importantly, the approximate chiral symmetry of the light quark sector.

The program: 1. Write the most general Lagrangian for nucleons and pions consistent with the symmetries of QCD (Lorentz invariance, chiral symmetry, parity, time reversal, etc.) 2. Organize the terms by chiral order (powers of $Q/\Lambda_\chi$) 3. Determine the low-energy constants (LECs) by fitting to experimental data ($NN$ scattering, $\pi N$ scattering, few-body observables) 4. Make predictions for nuclear structure using these forces as input to many-body calculations

Milestone achievement: The LENPIC (Low Energy Nuclear Physics International Collaboration) project has constructed chiral $NN$ potentials at N$^3$LO and N$^4$LO that reproduce the world $NN$ scattering database with $\chi^2/\text{datum} \approx 1$, and include consistent three-nucleon forces at N$^2$LO. These forces, used as input to ab initio many-body calculations, reproduce the ground-state energies and spectra of nuclei up to the calcium region.

Remaining challenge: Convergence of the chiral expansion. At the momenta relevant for nuclear binding ($Q \sim m_\pi$), the expansion parameter is $Q/\Lambda_\chi \sim 0.14$, suggesting rapid convergence. But certain observables — particularly the $^1S_0$ scattering length and the deuteron binding energy — are "fine-tuned" (unnaturally large), and the convergence pattern for these quantities is less clean.

Scale 4: Many-Body Nuclear Structure ($E \sim 1$--$20$ MeV, $r \sim 2$--$10$ fm)

Given nuclear forces from chiral EFT, the nuclear many-body problem must be solved. Modern ab initio methods include:

• Quantum Monte Carlo (QMC): Exact for $A \leq 12$ with local potentials (Argonne + Illinois)
• No-Core Shell Model (NCSM): Exact diagonalization for $A \leq 16$
• Coupled Cluster (CC): Polynomial-scaling method, applied to $^{16}$O, $^{40}$Ca, $^{48}$Ca, $^{56}$Ni, and beyond
• In-Medium Similarity Renormalization Group (IM-SRG): Applied to medium-mass nuclei up to $^{100}$Sn
• Self-Consistent Green's Function (SCGF): Accessed open-shell nuclei in the oxygen and calcium chains

Milestone achievement: Ab initio calculations with chiral EFT forces have reproduced the ground-state energies, charge radii, and excitation spectra of nuclei throughout the oxygen, calcium, and nickel isotopic chains, with quantified theoretical uncertainties. The prediction of the oxygen drip line at $^{24}$O — driven by three-nucleon forces — and the subsequent experimental confirmation was a landmark validation.

Remaining challenge: Heavy nuclei ($A > 100$), nuclear reactions, and the nuclear equation of state at supra-saturation density remain at the frontier.

A Concrete Example: The Triton Binding Energy

The triton ($^3$H) provides a concrete illustration of the emergence chain. Its binding energy, $B(^3$H$) = 8.482$ MeV, is one of the most precisely measured quantities in nuclear physics.

From QCD (lattice): At unphysically heavy pion masses ($m_\pi \sim 800$ MeV), the NPLQCD collaboration has computed the binding energy of the three-nucleon system, obtaining qualitative agreement with the expected trend. At the physical pion mass, the signal-to-noise problem makes this computation extremely expensive — no controlled result exists as of the mid-2020s.

From chiral EFT + exact few-body: Using chiral N$^3$LO two-body forces alone, the triton is underbound: $B(^3$H$)_{\text{NN only}} \approx 7.6$ MeV, about 10% below experiment. Adding the chiral N$^2$LO three-nucleon force (with the two parameters $c_D$ and $c_E$ fitted to reproduce $B(^3$H$)$ and the $nd$ doublet scattering length) gives the experimental binding energy by construction. But the critical test is predictive: with these parameters fixed, the $^4$He binding energy, the $A = 6$--$12$ spectra, the oxygen drip line, and the calcium isotope chain are all predicted — and the predictions agree with experiment at the few-percent level.

From phenomenological forces: The Argonne $v_{18}$ two-body potential gives $B(^3$H$) = 7.62$ MeV (underbound by 10%). Adding the Illinois-7 three-body force and solving the Faddeev equations exactly yields $B(^3$H$) = 8.48 \pm 0.01$ MeV — in precise agreement with experiment. This demonstrates that the physics of the three-nucleon force is correctly captured by both the phenomenological and chiral EFT frameworks.

The triton is the simplest system where the three-nucleon force is essential. Its binding energy serves as a calibration point for the entire nuclear structure program from chiral EFT.

The State of the Art: How Close Are We?

As of the mid-2020s, the chain from QCD to nuclear structure is:

The weakest link is the direct lattice QCD calculation of nuclear forces at physical quark masses. If this link could be made quantitative, the low-energy constants of chiral EFT could be determined from first principles rather than from scattering data — closing the chain from QCD to nuclei entirely.

A Second Concrete Example: Nuclear Matter Saturation

Perhaps the most dramatic test of the emergence chain is the saturation of nuclear matter — the observation that the interior of all heavy nuclei has approximately the same density ($\rho_0 \approx 0.16$ fm$^{-3}$) and the same binding energy per nucleon ($B/A \approx 16$ MeV).

With two-body forces only: When the nuclear matter equation of state is computed using any realistic two-nucleon potential (Argonne $v_{18}$, CD-Bonn, Nijmegen II, or chiral N$^3$LO), the saturation point falls on the "Coester band" — a narrow strip in the $(\rho_0, E/A)$ plane that runs from $\rho \sim 0.2$ fm$^{-3}$, $E/A \sim -12$ MeV to $\rho \sim 0.5$ fm$^{-3}$, $E/A \sim -25$ MeV. All potentials lie on this band, but none pass through the empirical saturation point. Nuclear matter is overbound and oversaturated with two-body forces alone.

With three-nucleon forces: Including chiral three-nucleon forces at N$^2$LO shifts the saturation point toward the empirical region. The repulsive component of the 3NF at high density (arising from the $c_E$ contact term) provides the additional pressure needed to prevent collapse. Modern calculations with chiral NN + 3NF give $\rho_0 \approx 0.16 \pm 0.01$ fm$^{-3}$ and $E/A \approx -15.5 \pm 1.0$ MeV — consistent with the empirical values.

This result has profound implications for neutron star physics (Chapter 25): the same three-nucleon forces that fix nuclear matter saturation also stiffen the equation of state of neutron-rich matter, supporting neutron star masses above $2 M_\odot$ — consistent with the observed pulsars PSR J0348+0432 ($M = 2.01 \pm 0.04 \, M_\odot$) and PSR J0740+6620 ($M = 2.08 \pm 0.07 \, M_\odot$).

Lessons

1. Known laws do not mean solved problems. QCD is the correct theory of the strong interaction, confirmed by decades of high-energy experiments. But deriving nuclear physics from QCD is a problem of extraordinary computational and conceptual difficulty, precisely because the coupling is strong at nuclear energy scales.

2. Effective theories are essential, not merely convenient. Chiral EFT is not a compromise or an admission of failure — it is the correct description of QCD at low energies, built on rigorous symmetry principles. The low-energy constants encode the effects of short-distance QCD in a systematically improvable way.

3. The emergence hierarchy is real. Nuclear physics is "derived from" QCD in the same sense that chemistry is derived from quantum electrodynamics — true in principle, but the emergent phenomena (nuclear binding, shell structure, collective motion) cannot be anticipated or understood without analysis at the appropriate level.

4. Computation is the new microscope. Lattice QCD, ab initio many-body methods, and the algorithms that drive them are as essential to modern nuclear physics as accelerators and detectors. The interplay between numerical computation and physical insight defines the frontier.

5. The chain will close. The trajectory is clear: lattice QCD will eventually compute the low-energy constants of chiral EFT with controlled uncertainties, eliminating the need to fit them to scattering data. When this happens — perhaps in the 2030s or 2040s with exascale computing and algorithmic advances — the chain from QCD to nuclear structure will be complete. We will be able to say, rigorously, that the chart of nuclides follows from six quark masses and one coupling constant.

Discussion Questions

1. The nuclear force is a residual strong force, analogous to the van der Waals force. But the van der Waals force can be computed perturbatively, while the nuclear force cannot. What specific property of QCD makes the nuclear emergence problem so much harder than the atomic one?

2. Should the chiral EFT program, which fits its low-energy constants to experimental $NN$ scattering data, be considered a "derivation from QCD" or a "parametrization of data"? What would constitute a genuine derivation?

3. The oxygen drip line prediction from chiral EFT three-nucleon forces was confirmed experimentally. Is this a test of QCD, of chiral EFT, or of the many-body calculation? Can these be separated?

4. If lattice QCD computing power increased by a factor of $10^6$ (roughly 20 years of Moore's law), what nuclear physics calculations would become possible that are currently intractable? Prioritize the most important ones.