Case Study 1 — The Deuteron: The Simplest Nucleus Teaches the Deepest Lessons

Introduction

The deuteron is the hydrogen atom of nuclear physics — the simplest system in which the fundamental interaction operates, and therefore the system from which we can extract the most unambiguous information about the force. But where the hydrogen atom reveals the clean, elegant structure of the Coulomb interaction, the deuteron reveals a force that is messy, complicated, and full of surprises. Every one of the deuteron's measured properties carries a lesson about nuclear physics.

This case study examines what we can learn from the deuteron's five key observables: binding energy, spin-parity, magnetic moment, quadrupole moment, and matter radius.

Lesson 1: The Binding Energy — A Nuclear Force Just Barely Strong Enough

The deuteron binding energy, $B_d = 2.22457$ MeV, is measured with extraordinary precision by the photodisintegration threshold:

$$\gamma + d \to p + n$$

The minimum photon energy required is exactly $B_d$ (in the center of mass; the laboratory threshold is slightly higher due to recoil). Chadwick and Goldhaber performed this measurement in 1934, just two years after the discovery of the neutron.

What makes $B_d$ remarkable is its smallness. The typical depth of the nucleon-nucleon potential well is $V_0 \sim 40$--$50$ MeV. The deuteron's binding energy is only about 5% of this depth. This means the deuteron is barely bound — it exists only because the nuclear force is just barely strong enough to create a bound state.

The consequences are far-reaching:

No excited states. In the square well model, the first excited state would require a potential depth roughly nine times larger than the minimum for any bound state. With $V_0/V_0^{\min} \approx 1.45$, the deuteron is far from this threshold.

Large spatial extent. The wavefunction decays as $e^{-\kappa r}$ outside the potential well, with $\kappa = 0.232$ fm$^{-1}$. The characteristic decay length $1/\kappa \approx 4.3$ fm is more than twice the range of the nuclear force. Roughly 70% of the deuteron's probability density lies outside the well. The deuteron is, in a sense, more "empty space" than "nucleus."

Fine-tuning. If the nuclear force were about 2% weaker, the deuteron would not exist. Without the deuteron, there would be no stepping stone from hydrogen to heavier elements in Big Bang nucleosynthesis. The existence of the chemical elements — and therefore of chemistry, biology, and us — depends on the nuclear force being within a narrow window of strength. This observation, sometimes called the "fine-tuning" of the nuclear force, is one of the data points in discussions of the anthropic principle. Nuclear physics provides the tightest anthropic constraints known.

Lesson 2: Spin and Parity — The Force Depends on Spin

The deuteron has $J^\pi = 1^+$. Since the neutron and proton each have spin 1/2, their total spin can be $S = 0$ (singlet) or $S = 1$ (triplet). The positive parity requires $L$ to be even: $L = 0, 2, 4, \ldots$ The quantum numbers $J = 1$ with $S = 1$ and $L = 0$ (the $^3S_1$ state) are the dominant configuration.

The critical fact is that the $S = 0$ (singlet) neutron-proton state is not bound. The nuclear force in the $^1S_0$ channel is almost — but not quite — strong enough to bind. The difference between the singlet and triplet channels is a direct measurement of the spin dependence of the nuclear force. Quantitatively:

• In the triplet channel, the potential depth (in a square well model) is about $V_0 \approx 35$ MeV.
• In the singlet channel, $V_0 \approx 25$ MeV — strong enough to produce a near-threshold virtual state (with $a_s = -23.7$ fm) but not a true bound state.

The spin dependence arises from the $\boldsymbol{\sigma}_1 \cdot \boldsymbol{\sigma}_2$ term in the nuclear potential. In the triplet state, $\langle \boldsymbol{\sigma}_1 \cdot \boldsymbol{\sigma}_2 \rangle = +1$, while in the singlet state, $\langle \boldsymbol{\sigma}_1 \cdot \boldsymbol{\sigma}_2 \rangle = -3$. A spin-spin interaction of the form $V_\sigma(r)\,\boldsymbol{\sigma}_1 \cdot \boldsymbol{\sigma}_2$ with $V_\sigma < 0$ (attractive) would make the triplet state more attractive, consistent with the data.

Lesson 3: The Magnetic Moment — Almost Right, But Not Quite

If the deuteron were a pure $^3S_1$ state (two spin-1/2 particles with $L = 0$, coupled to $S = J = 1$), its magnetic moment would be:

$$\mu_d^{(S)} = \mu_p + \mu_n = 2.7928 + (-1.9130) = 0.8798 \; \mu_N$$

The measured value is $\mu_d = 0.8574 \; \mu_N$, which differs from the $S$-state prediction by about 2.5%.

The discrepancy has two sources:

1. $D$-state admixture. The tensor force mixes $L = 0$ and $L = 2$ components. The magnetic moment of the $^3D_1$ state differs from the $^3S_1$ state because the orbital angular momentum contributes. A careful Clebsch-Gordan analysis gives, for the $D$-state contribution to $\mu$:

$$\mu_d = (1 - P_D)\mu_d^{(S)} + P_D \left[ \mu_d^{(D)} \right]$$

where $P_D$ is the $D$-state probability. The magnetic moment is sensitive to $P_D$ and provides an independent determination, though not a precise one because of the second effect.

2. Meson exchange currents (MEC). The nucleon magnetic moments are measured for free nucleons. Inside the deuteron, virtual mesons are exchanged between the nucleons, and these mesons carry electric current. The resulting "meson exchange currents" modify the electromagnetic properties of the bound system. The MEC correction to $\mu_d$ is of order 1--2% and has the same sign as the $D$-state correction, making it impossible to extract $P_D$ from $\mu_d$ alone without a model for the exchange currents.

This is a recurring lesson in nuclear physics: even the simplest observable in the simplest nucleus is affected by the full complexity of the nuclear force, including its meson exchange structure.

Lesson 4: The Quadrupole Moment — The Smoking Gun for the Tensor Force

The deuteron's electric quadrupole moment was measured by Kellogg, Rabi, Ramsey, and Zacharias in 1939 using molecular beam magnetic resonance on HD molecules. The result, $Q_d = 0.2860 \pm 0.0015$ fm$^2$, was a shock: it proved that the nuclear force is not purely central.

A central force (depending only on $r$, $\boldsymbol{\sigma}_1 \cdot \boldsymbol{\sigma}_2$, and $\boldsymbol{\tau}_1 \cdot \boldsymbol{\tau}_2$, but not on the direction of $\hat{r}$ relative to the spins) cannot mix different $L$ states. The $^3S_1$ state would remain pure, $Q = 0$. The nonzero $Q_d$ demands a non-central component — specifically, the tensor force with operator $S_{12}$.

The tensor operator $S_{12}$ has a characteristic angular structure:

$$S_{12} = 3(\boldsymbol{\sigma}_1 \cdot \hat{r})(\boldsymbol{\sigma}_2 \cdot \hat{r}) - \boldsymbol{\sigma}_1 \cdot \boldsymbol{\sigma}_2$$

which depends on the orientation of the spin vectors relative to the interparticle axis. In the $S = 1$ state, $S_{12}$ is nonzero and has nonzero matrix elements between $L = 0$ and $L = 2$ states with the same $J$:

$$\langle ^3D_1 | S_{12} | ^3S_1 \rangle \neq 0$$

The one-pion exchange potential (OPEP) naturally produces a tensor force of precisely this form. In fact, the OPEP tensor component dominates the $S$-$D$ mixing in the deuteron. The quadrupole moment is thus direct evidence for pion exchange.

The positive sign of $Q_d$ corresponds to a prolate deformation: the charge distribution is elongated along the spin axis. Physically, the tensor force from pion exchange is most attractive when the two nucleon spins are aligned along the line connecting them, stretching the deuteron along that axis.

Modern calculations with realistic potentials (Argonne $v_{18}$, chiral N$^3$LO) predict $Q_d$ in excellent agreement with experiment, provided both the $D$-state wavefunction and meson exchange current corrections are properly included. The $D$-state probability from these calculations is $P_D \approx 4$--$7$%, depending on the potential.

Lesson 5: The Matter Radius — A Loosely Bound Object

The deuteron's RMS charge radius is $r_{\text{ch}} = 2.128$ fm, and its matter radius is $r_m = 1.97$ fm. Both are substantially larger than the range of the nuclear force ($\sim 2$ fm) and much larger than the charge radii of the proton ($0.84$ fm) and neutron ($\sim 0$ fm net, but $\langle r^2 \rangle_n = -0.116$ fm$^2$).

The large radius is a direct consequence of the small binding energy. In the asymptotic region ($r > R$), the wavefunction is:

$$u(r) \propto e^{-\kappa r} \quad \text{with} \quad \kappa = \sqrt{2\mu B_d}/\hbar = 0.232 \text{ fm}^{-1}$$

The RMS radius is dominated by this exponential tail and scales as $\sim 1/(2\kappa) \approx 2.2$ fm. The deuteron is, in effect, two nucleons orbiting each other at a distance much larger than the range of their mutual attraction — like two people holding hands at arm's length while the "glue" connecting them has a reach of only a few centimeters.

This has practical implications. The deuteron is easily broken apart by photons of energy just above 2.225 MeV (photodisintegration). Its large size makes it a useful probe in nuclear reactions — deuteron stripping reactions ($d + A \to p + (A+1)$, where the neutron is captured) are among the most important tools for studying nuclear structure, precisely because the loosely bound deuteron readily donates one of its nucleons.

Synthesis: What the Deuteron Tells Us

No single nucleus teaches us more about the nuclear force. The deuteron is where the program of understanding nuclear interactions begins, and every advance — from Yukawa's meson theory to chiral EFT — is tested first against these five numbers.

The Deuteron as a Probe of Nuclear Structure

The deuteron is not only interesting in its own right — it is one of the most important experimental probes for studying other nuclei. Several techniques exploit its unique properties:

Deuteron stripping reactions ($d + A \to p + (A+1)^*$ or $d + A \to n + (A+1)^*$) transfer one nucleon from the loosely bound deuteron to the target nucleus. The angular distribution of the outgoing nucleon reveals the orbital angular momentum of the transferred nucleon, directly measuring the single-particle structure of the final nucleus. This technique, pioneered in the 1950s and 1960s, was instrumental in establishing the shell model (Chapter 6).

Deuteron breakup reactions ($d + A \to p + n + A^*$) are sensitive to the nuclear force between the two nucleons of the deuteron and the target, providing complementary information to elastic scattering.

Tensor analyzing powers in deuteron-nucleus scattering probe the tensor force directly, using the deuteron's spin-1 nature as a built-in tensor probe.

The deuteron thus serves a dual role in nuclear physics: it is both the simplest nuclear bound state and one of the most versatile experimental tools for studying nuclear structure.

The Anthropic Significance of the Deuteron

The deuteron occupies a unique position in discussions of fine-tuning in physics. Its marginal binding — $B_d/V_0 \approx 6\%$ — means that a small change in the strength of the nuclear force would either unbind it (eliminating the stepping stone to heavier elements in Big Bang nucleosynthesis) or produce additional bound states in the $nn$ and $pp$ channels (dramatically changing the landscape of nuclear physics and chemistry). Detailed calculations by Barr and Khan (2007) and others have shown that the window of nuclear force strengths compatible with a universe containing complex nuclei is narrow, perhaps 10--20% of the actual value. Whether this narrow window is a coincidence, a consequence of a deeper principle, or evidence for a multiverse is a question that lies outside nuclear physics proper — but the data come from the five numbers in the table above.

Discussion Questions

1. The deuteron has isospin $T = 0$. This means the $np$ pair is in a symmetric isospin state (analogous to the spin triplet). The $T = 1$ states include $pp$, $nn$, and the $T = 1$ component of $np$. Why does no $T = 1$ bound state exist? Is this a consequence of the nuclear force being weaker in the $T = 1$ channel, or of the Pauli principle, or both?

2. The deuteron's $D$-state probability is model-dependent: different potentials that fit the same scattering data give different values of $P_D$ (ranging from 4% to 7%). Why? What does this tell you about the relationship between observables and wavefunction components?

3. If the nuclear force had no tensor component, the deuteron would be a pure $S$-state with $Q_d = 0$. How would this affect nuclear structure more broadly? (Hint: consider the role of the tensor force in the shell model and in nuclear matter saturation.)