Case Study 1: Spin-Orbit Coupling — Where Angular Momentum Addition Meets Atomic Physics

Overview

The fine structure of hydrogen — the small but measurable splitting of energy levels that depends on the total angular momentum quantum number $j$ — is one of the most celebrated results in atomic physics. It was the first place where the addition of angular momentum moved from abstract mathematics to observable spectral lines. This case study traces the full story: from the classical magnetic interaction that causes spin-orbit coupling, through the quantum mechanical treatment using the coupled basis machinery of Chapter 14, to the precision measurements that eventually revealed the limitations of even this refined picture and led to quantum electrodynamics (QED).

Part 1: The Classical Origin of Spin-Orbit Coupling

The Electron's View

Consider an electron orbiting a proton in hydrogen. In the electron's rest frame, the proton appears to orbit the electron. This moving positive charge creates a magnetic field at the electron's position:

$$\mathbf{B} = -\frac{\mathbf{v} \times \mathbf{E}}{c^2}$$

where $\mathbf{v}$ is the electron's velocity and $\mathbf{E}$ is the electric field from the proton. Since the Coulomb field is radial, $\mathbf{E} = \frac{e}{4\pi\epsilon_0 r^3}\mathbf{r}$, the magnetic field is:

$$\mathbf{B} = \frac{e}{4\pi\epsilon_0 m_e c^2 r^3}\mathbf{L}$$

where $\mathbf{L} = m_e \mathbf{r} \times \mathbf{v}$ is the orbital angular momentum.

The electron's spin magnetic moment interacts with this field:

$$H_{\text{SO}} = -\boldsymbol{\mu}_s \cdot \mathbf{B} = \frac{e^2}{4\pi\epsilon_0 m_e^2 c^2 r^3}\mathbf{S} \cdot \mathbf{L}$$

The Thomas Correction

There is a subtlety: the electron's rest frame is not an inertial frame (it is accelerating due to the Coulomb force). Llewellyn Thomas showed in 1926 that the non-inertial frame introduces an additional precession — the Thomas precession — which reduces the spin-orbit coupling by exactly a factor of 2. The correct expression is:

$$\hat{H}_{\text{SO}} = \frac{e^2}{8\pi\epsilon_0 m_e^2 c^2 r^3}\hat{\mathbf{L}} \cdot \hat{\mathbf{S}} = \frac{1}{2m_e^2 c^2}\frac{1}{r}\frac{dV}{dr}\hat{\mathbf{L}} \cdot \hat{\mathbf{S}}$$

🔵 Historical Note: Thomas's 1926 paper resolved a factor-of-2 discrepancy that had plagued the theory for months. When Goudsmit and Uhlenbeck proposed electron spin in 1925, their prediction for the fine structure splitting was off by a factor of 2 from experiment. Thomas showed that the discrepancy arose from neglecting the relativistic precession of the electron's frame, not from the spin hypothesis itself. This was one of the first successes of combining special relativity with quantum mechanics.

Part 2: The Quantum Mechanical Treatment

Why the Uncoupled Basis Fails

In the uncoupled basis $|n, \ell, m_\ell, s, m_s\rangle$, the spin-orbit Hamiltonian is:

$$\hat{H}_{\text{SO}} = A(r) \hat{\mathbf{L}} \cdot \hat{\mathbf{S}}$$

where $A(r) > 0$ for hydrogen. The operator $\hat{\mathbf{L}} \cdot \hat{\mathbf{S}}$ is not diagonal in the uncoupled basis. To see this, write:

$$\hat{\mathbf{L}} \cdot \hat{\mathbf{S}} = \hat{L}_z \hat{S}_z + \frac{1}{2}(\hat{L}_+ \hat{S}_- + \hat{L}_- \hat{S}_+)$$

The terms $\hat{L}_+ \hat{S}_-$ and $\hat{L}_- \hat{S}_+$ connect states with different $m_\ell$ and $m_s$ (while preserving $m_j = m_\ell + m_s$). For example:

$$\hat{L}_+ \hat{S}_- |1, 0; \tfrac{1}{2}, \tfrac{1}{2}\rangle = \sqrt{2}\hbar \cdot \hbar |1, 1; \tfrac{1}{2}, -\tfrac{1}{2}\rangle$$

This means $m_\ell$ and $m_s$ are not good quantum numbers in the presence of spin-orbit coupling. But $m_j = m_\ell + m_s$ is conserved, because $[\hat{H}_{\text{SO}}, \hat{J}_z] = 0$.

The Coupled Basis Solution

In the coupled basis $|n, \ell, j, m_j\rangle$, the operator $\hat{\mathbf{L}} \cdot \hat{\mathbf{S}}$ is diagonal — and its eigenvalue is trivially computed:

$$\hat{\mathbf{L}} \cdot \hat{\mathbf{S}} = \frac{1}{2}(\hat{J}^2 - \hat{L}^2 - \hat{S}^2)$$

$$\langle \hat{\mathbf{L}} \cdot \hat{\mathbf{S}} \rangle = \frac{\hbar^2}{2}[j(j+1) - \ell(\ell+1) - s(s+1)]$$

This is the payoff of Chapter 14's coupling machinery: what is a complicated off-diagonal matrix in the uncoupled basis becomes a single number in the coupled basis.

Complete Fine Structure of the $n = 2$ Level

The unperturbed $n = 2$ level of hydrogen has eight degenerate states:

• $2s$ ($\ell = 0$): 2 states (spin up and down)
• $2p$ ($\ell = 1$): 6 states ($m_\ell = -1, 0, 1$ each with spin up/down)

Spin-orbit coupling splits these as follows:

$2s_{1/2}$ ($\ell = 0$, $j = 1/2$): No spin-orbit shift, since $\ell = 0$ implies $\hat{\mathbf{L}} = 0$.

$$\Delta E = 0$$

$2p_{1/2}$ ($\ell = 1$, $j = 1/2$): $$\langle \hat{\mathbf{L}} \cdot \hat{\mathbf{S}} \rangle = \frac{\hbar^2}{2}\left[\frac{3}{4} - 2 - \frac{3}{4}\right] = -\hbar^2$$

$$\Delta E_{1/2} = -A_{21} \hbar^2$$

where $A_{21} = \langle A(r) \rangle_{n=2, \ell=1}$ is the radial expectation value.

$2p_{3/2}$ ($\ell = 1$, $j = 3/2$): $$\langle \hat{\mathbf{L}} \cdot \hat{\mathbf{S}} \rangle = \frac{\hbar^2}{2}\left[\frac{15}{4} - 2 - \frac{3}{4}\right] = \frac{\hbar^2}{2}$$

$$\Delta E_{3/2} = +\frac{1}{2}A_{21}\hbar^2$$

The splitting between $2p_{3/2}$ and $2p_{1/2}$ is:

$$\Delta E_{3/2} - \Delta E_{1/2} = \frac{3}{2}A_{21}\hbar^2$$

The Interval Rule

A general pattern emerges: for a given $\ell$, the splitting between adjacent $j$ levels satisfies the Lande interval rule:

$$E(j) - E(j-1) \propto j$$

This rule is a direct consequence of the $\hat{\mathbf{L}} \cdot \hat{\mathbf{S}}$ form of the interaction and the eigenvalue formula for the coupled basis. It provides an experimental test: if observed splittings do not follow the interval rule, additional interactions beyond spin-orbit coupling are present.

Part 3: From Hydrogen to Multi-Electron Atoms

L-S Coupling (Russell-Saunders Coupling)

For light atoms (roughly $Z \lesssim 30$), the residual Coulomb interaction between electrons is much stronger than spin-orbit coupling. The appropriate coupling scheme is L-S coupling:

1. Couple all orbital angular momenta: $\hat{\mathbf{L}} = \sum_i \hat{\boldsymbol{\ell}}_i$
2. Couple all spin angular momenta: $\hat{\mathbf{S}} = \sum_i \hat{\mathbf{s}}_i$
3. Couple $\hat{\mathbf{L}}$ and $\hat{\mathbf{S}}$ to form $\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}$

States are labeled by term symbols $^{2S+1}L_J$, where $L$ is written as a letter ($S, P, D, F, \ldots$ for $L = 0, 1, 2, 3, \ldots$).

Example: Carbon ($Z = 6$)

The ground configuration is $1s^2 2s^2 2p^2$. The two $2p$ electrons can couple to: - $L = 0, 1, 2$ and $S = 0, 1$ - Pauli exclusion restricts the allowed combinations - Hund's rules predict the ground term: $^3P_0$ (triplet $P$, $J = 0$)

The fine structure of the $^3P$ term produces three levels: - $^3P_0$ ($J = 0$): ground state - $^3P_1$ ($J = 1$): $16.4\,\text{cm}^{-1}$ above ground - $^3P_2$ ($J = 2$): $43.4\,\text{cm}^{-1}$ above ground

The ratio of splittings: $(43.4 - 16.4) / 16.4 = 1.65$, close to the Lande interval rule prediction of $2/1 = 2$. The deviation signals that higher-order effects (spin-spin, spin-other-orbit) are not negligible.

j-j Coupling

For heavy atoms ($Z \gtrsim 70$), the spin-orbit coupling becomes comparable to or stronger than the residual Coulomb interaction. The appropriate scheme is j-j coupling:

1. For each electron $i$, couple $\hat{\boldsymbol{\ell}}_i$ and $\hat{\mathbf{s}}_i$ to form $\hat{\mathbf{j}}_i$
2. Couple all $\hat{\mathbf{j}}_i$ to form $\hat{\mathbf{J}} = \sum_i \hat{\mathbf{j}}_i$

States are labeled by the individual $j_i$ values.

Example: Lead ($Z = 82$)

The ground configuration includes $6p^2$. In j-j coupling: - Each $6p$ electron has $\ell = 1$, $s = 1/2$, so $j = 1/2$ or $3/2$ - The ground state has both electrons in $j = 1/2$: $(j_1 = 1/2, j_2 = 1/2)_{J=0}$

The energy levels of lead's $6p^2$ configuration do not follow the Lande interval rule — a clear signature that L-S coupling has broken down and j-j coupling is more appropriate.

The Transition Region

Most elements fall between pure L-S and pure j-j coupling. The real coupling scheme is intermediate coupling, where neither set of quantum numbers is perfectly good. The CG coefficients provide the mathematical framework for transforming between the two limiting descriptions.

Part 4: Precision Measurements and the Lamb Shift

The Dirac Prediction

In 1928, Paul Dirac's relativistic quantum mechanics predicted that the hydrogen fine structure depends only on $j$ and $n$ — not on $\ell$ separately. Specifically, levels with the same $n$ and $j$ but different $\ell$ (such as $2s_{1/2}$ and $2p_{1/2}$) should be exactly degenerate.

The Lamb-Retherford Experiment

In 1947, Willis Lamb and Robert Retherford used microwave spectroscopy to show that the $2s_{1/2}$ level lies approximately $1057\,\text{MHz}$ (about $4.4 \times 10^{-6}\,\text{eV}$) above the $2p_{1/2}$ level. This tiny but nonzero energy difference — the Lamb shift — could not be explained by the Dirac equation alone.

📊 By the Numbers: The Lamb shift is tiny compared to the fine structure splitting ($\sim 10^4\,\text{MHz}$ for $2p_{3/2}$ vs. $2p_{1/2}$), which is tiny compared to the Bohr energy ($\sim 10^{15}\,\text{Hz}$ for the $n = 2$ level). Yet it was measurable, and its explanation required a fundamental advance in theoretical physics.

Quantum Electrodynamics

The Lamb shift arises from the electron's interaction with the quantum vacuum — virtual photon emission and reabsorption, and virtual electron-positron pair creation. These are quantum electrodynamic (QED) effects that go beyond the Dirac equation.

Hans Bethe, in a famous calculation performed on the train from the 1947 Shelter Island conference, computed the non-relativistic contribution to the Lamb shift and obtained $1040\,\text{MHz}$ — remarkably close to the measured $1057\,\text{MHz}$. The full relativistic QED calculation, carried out by Julian Schwinger, Richard Feynman, and Sin-Itiro Tomonaga, gave excellent agreement and established QED as the most precisely tested theory in physics.

Part 5: Modern Applications

Nuclear Magnetic Resonance (NMR) and MRI

Nuclear spin-orbit coupling — the coupling between nuclear orbital and spin angular momenta — determines nuclear energy levels and magnetic moments. The Lande g-factor and the projection theorem (both consequences of angular momentum coupling theory) are essential for predicting NMR frequencies. Every MRI machine in every hospital relies on the physics of angular momentum coupling.

Atomic Clocks

Modern atomic clocks, including the cesium-133 standard that defines the second, rely on transitions between hyperfine levels of atoms. Hyperfine structure arises from the coupling of the electron's total angular momentum $\hat{\mathbf{J}}$ with the nuclear spin $\hat{\mathbf{I}}$ to form $\hat{\mathbf{F}} = \hat{\mathbf{J}} + \hat{\mathbf{I}}$. This is a direct application of angular momentum addition, with CG coefficients determining the transition strengths.

Quantum Computing with Trapped Ions

Trapped-ion quantum computers (such as those built by IonQ, Quantinuum, and others) encode qubits in the hyperfine levels of ions like $^{171}\text{Yb}^+$ or $^{43}\text{Ca}^+$. The qubit states are specific $|F, m_F\rangle$ levels, and the gate operations are transitions whose strengths are governed by CG coefficients. Angular momentum coupling theory is quite literally the language in which these quantum computers are programmed.

Discussion Questions

1. The Lamb shift was measured in 1947, but the theoretical tools to explain it (QED) were developed almost simultaneously. Is it fair to say that experiment drove theory, or was theory already anticipating such effects? How does this compare to the historical sequence for the discovery of spin?

2. The transition from L-S to j-j coupling as atoms get heavier is a smooth crossover, not a sharp phase transition. What does this tell us about the nature of "good quantum numbers" — are they absolute properties of the system, or convenient approximations?

3. If we could turn off the spin-orbit interaction in all atoms, how would the periodic table change? Would chemistry be significantly different?

4. The Lande g-factor formula assumes pure L-S coupling. For atoms with intermediate coupling, $g_j$ deviates from the formula. How could you use measured g-factors to determine the degree of configuration mixing in a real atom?