Chapter 14: Addition of Angular Momentum — Clebsch-Gordan Coefficients and Coupled Representations

"The theory of angular momentum coupling is the grammar of atomic and nuclear physics." — Eugene Wigner

"Angular momentum is to quantum mechanics what Fourier analysis is to signal processing — a universal language that appears in every problem." — Adapted from J. J. Sakurai

You have spent Chapters 12 and 13 mastering the angular momentum algebra of a single degree of freedom: the eigenvalue spectrum of $\hat{J}^2$ and $\hat{J}_z$ from commutation relations alone, the beautiful matrix representations for any spin $j$, Pauli matrices, the Bloch sphere, Larmor precession. You know how to handle one angular momentum with precision and confidence.

Nature, however, rarely hands you just one angular momentum.

An electron in hydrogen has orbital angular momentum $\hat{\mathbf{L}}$ and spin $\hat{\mathbf{S}}$. A helium atom has two electron spins $\hat{\mathbf{S}}_1$ and $\hat{\mathbf{S}}_2$. A deuterium nucleus has proton spin, neutron spin, and orbital angular momentum of the relative motion. Every multi-particle system, every atom beyond hydrogen, every nucleus, every hadron — all of them require you to combine angular momenta.

This chapter teaches you how. The payoff is enormous: once you learn to couple angular momenta, you will understand the fine structure of hydrogen, the term symbols of multi-electron atoms (Chapter 16), selection rules for spectroscopic transitions, and the deep mathematical machinery that underpins nuclear and particle physics.

The mathematical framework we develop here — Clebsch-Gordan coefficients, the Wigner-Eckart theorem, irreducible tensor operators — is not merely a technical toolkit. It is the natural language of rotational symmetry, and rotational symmetry is the most fundamental symmetry in atomic, molecular, nuclear, and particle physics. Learning this language will unlock doors that remain permanently closed without it. The investment is substantial; the return is extraordinary.

🏃 Fast Track: If you need the results quickly for later chapters, focus on Sections 14.1–14.5 (the coupling procedure and Clebsch-Gordan coefficients) and Section 14.7 (the Wigner-Eckart theorem). Sections 14.6 and 14.8 provide applications and depth but can be deferred.

14.1 Why We Need to Add Angular Momenta

The Physical Situation

Consider an electron in the hydrogen atom. From Chapter 5, you know it has:

• Orbital angular momentum $\hat{\mathbf{L}}$ with quantum numbers $\ell$ and $m_\ell$, arising from its motion around the proton.
• Spin angular momentum $\hat{\mathbf{S}}$ with quantum number $s = 1/2$ and $m_s = \pm 1/2$, an intrinsic property discovered in Chapter 13.

These are two independent angular momenta living in two different Hilbert spaces. The orbital part acts on the spherical harmonics $Y_\ell^{m_\ell}(\theta, \phi)$; the spin part acts on the two-dimensional spinor space $\{|\!\uparrow\rangle, |\!\downarrow\rangle\}$. The full state lives in the tensor product:

$$|\psi\rangle \in \mathcal{H}_{\text{orbital}} \otimes \mathcal{H}_{\text{spin}}$$

In the absence of any coupling between $\hat{\mathbf{L}}$ and $\hat{\mathbf{S}}$, the uncoupled basis $|\ell, m_\ell\rangle \otimes |s, m_s\rangle = |\ell, m_\ell; s, m_s\rangle$ is perfectly adequate. The two angular momenta are independent, and $m_\ell$ and $m_s$ are separately good quantum numbers.

But the moment you introduce the spin-orbit interaction — a relativistic correction proportional to $\hat{\mathbf{L}} \cdot \hat{\mathbf{S}}$ — these individual projections are no longer conserved. Instead, the total angular momentum

$$\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}$$

is conserved, and the good quantum numbers are $J$ (total) and $M$ (its $z$-projection).

🔗 Connection: Recall from Chapter 5 that the hydrogen atom energy levels $E_n = -13.6\,\text{eV}/n^2$ depend only on the principal quantum number $n$. The spin-orbit coupling lifts this degeneracy, splitting each level into states labeled by $J$. This is the fine structure, which we will study quantitatively in Chapter 18 using the angular momentum coupling tools built in this chapter.

The Central Question

Given two angular momenta $\hat{\mathbf{J}}_1$ (with quantum number $j_1$) and $\hat{\mathbf{J}}_2$ (with quantum number $j_2$), we want to form the total angular momentum $\hat{\mathbf{J}} = \hat{\mathbf{J}}_1 + \hat{\mathbf{J}}_2$ and find:

1. What values of $J$ are allowed?
2. How do we express the coupled states $|J, M\rangle$ in terms of the uncoupled states $|j_1, m_1; j_2, m_2\rangle$?
3. What are the coefficients in this expansion?

Those coefficients are the Clebsch-Gordan (CG) coefficients, and they are the central mathematical objects of this chapter.

🔗 Spaced Review — Chapter 12: The eigenvalue spectrum of angular momentum. Recall that for any angular momentum $\hat{\mathbf{J}}$, we have $\hat{J}^2|j, m\rangle = j(j+1)\hbar^2|j, m\rangle$ and $\hat{J}_z|j, m\rangle = m\hbar|j, m\rangle$, where $j = 0, 1/2, 1, 3/2, \ldots$ and $m = -j, -j+1, \ldots, j$. The raising and lowering operators $\hat{J}_\pm|j, m\rangle = \hbar\sqrt{j(j+1) - m(m\pm 1)}|j, m\pm 1\rangle$ will be essential tools throughout this chapter.

Commutation Relations of the Total Angular Momentum

A critical mathematical fact underpins everything that follows. If $\hat{\mathbf{J}}_1$ and $\hat{\mathbf{J}}_2$ each satisfy the angular momentum commutation relations, and if they commute with each other (because they act on different Hilbert spaces), then the total $\hat{\mathbf{J}} = \hat{\mathbf{J}}_1 + \hat{\mathbf{J}}_2$ also satisfies the angular momentum commutation relations:

$$[\hat{J}_x, \hat{J}_y] = [\hat{J}_{1x} + \hat{J}_{2x}, \hat{J}_{1y} + \hat{J}_{2y}] = [\hat{J}_{1x}, \hat{J}_{1y}] + [\hat{J}_{2x}, \hat{J}_{2y}] = i\hbar\hat{J}_{1z} + i\hbar\hat{J}_{2z} = i\hbar\hat{J}_z$$

The cross-terms vanish because $[\hat{J}_{1x}, \hat{J}_{2y}] = 0$. This means $\hat{\mathbf{J}}$ is a perfectly valid angular momentum operator. Its eigenvalue spectrum must follow the standard pattern from Chapter 12: $\hat{J}^2|J, M\rangle = J(J+1)\hbar^2|J, M\rangle$ with $M = -J, -J+1, \ldots, J$.

However, knowing the commutation relations only tells us that $J$ must be a non-negative integer or half-integer. It does not tell us which values of $J$ actually appear when we combine specific $j_1$ and $j_2$. Answering that question is the content of the triangle rule (Section 14.3).

Let us also record which operators commute with which. This is essential for understanding the two bases:

$$[\hat{J}^2, \hat{J}_{1z}] \neq 0, \qquad [\hat{J}^2, \hat{J}_{2z}] \neq 0$$

This is why you cannot specify $J$ and $m_1$ simultaneously: the total angular momentum magnitude does not commute with the individual projections. To see this explicitly, note that:

$$\hat{J}^2 = \hat{J}_1^2 + \hat{J}_2^2 + 2\hat{\mathbf{J}}_1 \cdot \hat{\mathbf{J}}_2 = \hat{J}_1^2 + \hat{J}_2^2 + 2\hat{J}_{1z}\hat{J}_{2z} + \hat{J}_{1+}\hat{J}_{2-} + \hat{J}_{1-}\hat{J}_{2+}$$

The terms $\hat{J}_{1+}\hat{J}_{2-}$ and $\hat{J}_{1-}\hat{J}_{2+}$ raise one projection while lowering the other. They preserve $M = m_1 + m_2$ but change $m_1$ and $m_2$ individually, which is why $\hat{J}^2$ does not commute with $\hat{J}_{1z}$ or $\hat{J}_{2z}$.

14.2 The Problem: $|j_1, m_1\rangle \otimes |j_2, m_2\rangle \rightarrow |J, M\rangle$

Two Bases for the Same Space

The tensor product space $\mathcal{H}_{j_1} \otimes \mathcal{H}_{j_2}$ has dimension $(2j_1 + 1)(2j_2 + 1)$. We can describe any state in this space using either of two complete orthonormal bases:

Uncoupled basis: $\{|j_1, m_1; j_2, m_2\rangle\}$, where $m_1$ ranges from $-j_1$ to $j_1$ and $m_2$ ranges from $-j_2$ to $j_2$. These states are simultaneous eigenstates of $\hat{J}_1^2$, $\hat{J}_{1z}$, $\hat{J}_2^2$, and $\hat{J}_{2z}$.

Coupled basis: $\{|J, M\rangle\}$, where $J$ ranges over the allowed values (to be determined) and $M$ ranges from $-J$ to $J$. These states are simultaneous eigenstates of $\hat{J}_1^2$, $\hat{J}_2^2$, $\hat{J}^2$, and $\hat{J}_z$.

Notice that both bases share the operators $\hat{J}_1^2$ and $\hat{J}_2^2$ — the magnitudes of the individual angular momenta are always well-defined. What differs is whether we specify the individual projections ($m_1, m_2$) or the total magnitude and projection ($J, M$).

💡 Key Insight: You cannot simultaneously specify $m_1$, $m_2$, $J$, and $M$ (except in trivial cases) because $\hat{J}_{1z}$ and $\hat{J}^2$ do not commute. The choice between the uncoupled and coupled bases is a choice between knowing the individual projections and knowing the total angular momentum.

The Unitary Transformation

The two bases span the same space, so they are related by a unitary transformation:

$$|J, M\rangle = \sum_{m_1, m_2} |j_1, m_1; j_2, m_2\rangle \langle j_1, m_1; j_2, m_2 | J, M\rangle$$

The coefficients $\langle j_1, m_1; j_2, m_2 | J, M\rangle$ are the Clebsch-Gordan coefficients. They encode the complete relationship between the two bases.

The inverse transformation is:

$$|j_1, m_1; j_2, m_2\rangle = \sum_{J, M} |J, M\rangle \langle J, M | j_1, m_1; j_2, m_2\rangle$$

Since the transformation is unitary (a change of orthonormal basis), the CG coefficients satisfy:

$$\sum_{m_1, m_2} |\langle j_1, m_1; j_2, m_2 | J, M\rangle|^2 = 1$$

$$\sum_{J, M} |\langle j_1, m_1; j_2, m_2 | J, M\rangle|^2 = 1$$

These orthogonality relations will be essential for checking your calculations.

The $M = m_1 + m_2$ Selection Rule

One property follows immediately from the fact that $\hat{J}_z = \hat{J}_{1z} + \hat{J}_{2z}$:

$$\langle j_1, m_1; j_2, m_2 | J, M\rangle = 0 \quad \text{unless} \quad M = m_1 + m_2$$

This is because both sides must be eigenstates of $\hat{J}_z$ with the same eigenvalue. The total $z$-projection is always the sum of the individual projections. This constraint dramatically reduces the number of nonzero CG coefficients.

✅ Checkpoint: Before continuing, make sure you can explain why the uncoupled basis has exactly $(2j_1 + 1)(2j_2 + 1)$ states. Can you verify that the coupled basis has the same total number? (You will, once you learn which values of $J$ are allowed.)

14.3 Allowed Values of $J$: The Triangle Rule

Dimension Counting

What values of $J$ appear in the coupled basis? Consider the maximum possible $M$: the largest value of $m_1 + m_2$ is $j_1 + j_2$, achieved when $m_1 = j_1$ and $m_2 = j_2$. This state must belong to a multiplet with $J = j_1 + j_2$ (since $M_{\max} = J$ for any angular momentum multiplet).

Similarly, the minimum value of $J$ cannot be zero (in general) — there is a lower bound. The allowed values turn out to be:

$$J = j_1 + j_2, \quad j_1 + j_2 - 1, \quad j_1 + j_2 - 2, \quad \ldots, \quad |j_1 - j_2|$$

This is the triangle rule (or triangle inequality): $J$ takes all values from $|j_1 - j_2|$ to $j_1 + j_2$ in integer steps.

Proof by Dimension Count

The total number of states in the coupled basis must equal $(2j_1 + 1)(2j_2 + 1)$. Let us verify. For each allowed $J$, there are $2J + 1$ states (the $M$ values from $-J$ to $J$). So:

$$\sum_{J = |j_1 - j_2|}^{j_1 + j_2} (2J + 1) = (2j_1 + 1)(2j_2 + 1)$$

Let us verify this with $j_1 \geq j_2$ (without loss of generality). Set $J_{\min} = j_1 - j_2$ and $J_{\max} = j_1 + j_2$. The number of terms in the sum is $J_{\max} - J_{\min} + 1 = 2j_2 + 1$. Using the formula for the sum of an arithmetic series:

$$\sum_{J = j_1 - j_2}^{j_1 + j_2} (2J + 1) = (2j_2 + 1) \cdot \frac{(2(j_1 - j_2) + 1) + (2(j_1 + j_2) + 1)}{2}$$

$$= (2j_2 + 1) \cdot \frac{(2j_1 - 2j_2 + 1) + (2j_1 + 2j_2 + 1)}{2} = (2j_2 + 1)(2j_1 + 1)$$

The dimensions match. This is a necessary consistency check, and it works because the triangle rule is correct.

Visualizing the Coupling: $M$-Value Decomposition

Another way to derive the triangle rule uses what we might call the "$M$-counting" argument. In the uncoupled basis, the number of states with a given $M = m_1 + m_2$ equals the number of ways to choose $m_1$ and $m_2$ such that $m_1 + m_2 = M$, with $|m_1| \leq j_1$ and $|m_2| \leq j_2$.

For the maximum $M = j_1 + j_2$, there is exactly one uncoupled state: $|j_1, j_1; j_2, j_2\rangle$. This must equal $|J = j_1 + j_2, M = j_1 + j_2\rangle$.

For $M = j_1 + j_2 - 1$, there are exactly two uncoupled states: $|j_1, j_1 - 1; j_2, j_2\rangle$ and $|j_1, j_1; j_2, j_2 - 1\rangle$. In the coupled basis, one of these belongs to the $J = j_1 + j_2$ multiplet (the $M = J - 1$ state), and the other must belong to a new multiplet with $J = j_1 + j_2 - 1$.

Continuing this argument: at each step, one new multiplet begins, until we reach $M = j_1 - j_2$ (assuming $j_1 \geq j_2$), after which the number of uncoupled states starts decreasing. The result is exactly the triangle rule.

📊 By the Numbers — Coupling Examples:

$j_1$	$j_2$	Allowed $J$	Dim check
$1/2$	$1/2$	$1, 0$	$(2)(2) = (3)+(1) = 4$ ✓
$1$	$1/2$	$3/2, 1/2$	$(3)(2) = (4)+(2) = 6$ ✓
$1$	$1$	$2, 1, 0$	$(3)(3) = (5)+(3)+(1) = 9$ ✓
$2$	$1$	$3, 2, 1$	$(5)(3) = (7)+(5)+(3) = 15$ ✓
$3/2$	$1/2$	$2, 1$	$(4)(2) = (5)+(3) = 8$ ✓

⚠️ Common Misconception: The triangle rule does not mean $J = j_1 + j_2$ only. This is a common error. The vector model analogy — "adding" angular momentum vectors — would suggest a fixed resultant. But in quantum mechanics, the magnitudes and projections follow different rules. The total angular momentum magnitude $\sqrt{J(J+1)}\hbar$ is generally less than $\sqrt{j_1(j_1+1)}\hbar + \sqrt{j_2(j_2+1)}\hbar$, and $J$ ranges over all allowed values, not a single one.

The Vector Model: Useful Intuition, Dangerous Literalism

A classical analogy sometimes helps build intuition, even though it is ultimately misleading. Imagine two classical angular momentum vectors $\mathbf{J}_1$ and $\mathbf{J}_2$ of fixed lengths $\sqrt{j_1(j_1+1)}\hbar$ and $\sqrt{j_2(j_2+1)}\hbar$, respectively. Their resultant $\mathbf{J} = \mathbf{J}_1 + \mathbf{J}_2$ has a length that depends on the angle between them:

• When the vectors are parallel, $|\mathbf{J}| = |\mathbf{J}_1| + |\mathbf{J}_2|$ (maximum).
• When the vectors are antiparallel, $|\mathbf{J}| = ||\mathbf{J}_1| - |\mathbf{J}_2||$ (minimum).
• At intermediate angles, $|\mathbf{J}|$ takes intermediate values.

The quantum triangle rule $|j_1 - j_2| \leq J \leq j_1 + j_2$ captures this range, but with a crucial difference: quantum mechanically, the vectors are never fully parallel or antiparallel. The angular momentum components perpendicular to $\hat{z}$ are always uncertain (recall from Chapter 12 that $[\hat{J}_x, \hat{J}_y] = i\hbar\hat{J}_z$, so specifying $J_z$ leaves $J_x$ and $J_y$ uncertain). The "vector model" should be understood as a mnemonic for the triangle rule, not as a description of what the angular momentum vectors are actually doing.

In fact, you can prove that $\sqrt{J(J+1)}\hbar < \sqrt{j_1(j_1+1)}\hbar + \sqrt{j_2(j_2+1)}\hbar$ for all allowed $J$ — the total angular momentum magnitude is always less than the naive sum of magnitudes. This is because perfect alignment of quantum angular momentum vectors is impossible.

🔵 Historical Note: The vector model was developed by Arnold Sommerfeld and Alfred Lande in the early 1920s, before the modern quantum theory was established. Despite its classical flavor, it proved remarkably effective at organizing atomic spectroscopic data. The reason it works is precisely the content of this chapter: the addition of angular momenta obeys strict algebraic rules (the triangle rule, CG coefficients) that the vector model approximates pictorially.

Formal Proof of the Triangle Rule via Commutators

The triangle rule can be proven rigorously from the commutation relations alone. The key steps:

1. $M_{\max} = j_1 + j_2$: Since $\hat{J}_z = \hat{J}_{1z} + \hat{J}_{2z}$ and the maximum eigenvalues of $\hat{J}_{1z}$ and $\hat{J}_{2z}$ are $j_1\hbar$ and $j_2\hbar$, the maximum eigenvalue of $\hat{J}_z$ is $(j_1 + j_2)\hbar$. This state must belong to a $J = j_1 + j_2$ multiplet, because $M_{\max} = J$ within any multiplet.

2. $J_{\max} = j_1 + j_2$: We just established that $J = j_1 + j_2$ must be one of the allowed values. Can $J$ be larger? No, because that would require $M > j_1 + j_2$, which is impossible.

3. $J_{\min} = |j_1 - j_2|$: The number of uncoupled states with a given $M$ equals the number of ways to write $M = m_1 + m_2$ with $|m_1| \leq j_1$ and $|m_2| \leq j_2$. At each step in the $M$-counting argument, one new $J$ multiplet begins. The process terminates when we have exhausted all states, which happens at $J = |j_1 - j_2|$.

4. All integer steps: The step size between successive $J$ values must be 1 (integer), because each new $J$ multiplet that opens as $M$ decreases from $j_1 + j_2$ opens at one integer step below the previous one.

14.4 Clebsch-Gordan Coefficients: Definition and Calculation

Notation and Conventions

The Clebsch-Gordan coefficient $\langle j_1, m_1; j_2, m_2 | J, M\rangle$ is the overlap between an uncoupled state and a coupled state. Different textbooks use different notations:

We adopt the Condon-Shortley phase convention, which is the most widely used. Under this convention, all CG coefficients are real, and $\langle j_1, j_1; j_2, J - j_1 | J, J\rangle > 0$.

The Recursion Method

The standard method for calculating CG coefficients uses the lowering operator $\hat{J}_- = \hat{J}_{1-} + \hat{J}_{2-}$. The strategy:

1. Start at the top. The "stretched" state $|J = j_1 + j_2, M = j_1 + j_2\rangle$ is unique — it equals the uncoupled state $|j_1, j_1; j_2, j_2\rangle$:

$$|j_1 + j_2, j_1 + j_2\rangle = |j_1, j_1; j_2, j_2\rangle$$

This gives us one CG coefficient for free: $\langle j_1, j_1; j_2, j_2 | j_1 + j_2, j_1 + j_2\rangle = 1$.

2. Apply $\hat{J}_-$ to generate lower $M$ states. Acting with $\hat{J}_- = \hat{J}_{1-} + \hat{J}_{2-}$ on both sides:

$$\hat{J}_-|J, M\rangle = \hbar\sqrt{J(J+1) - M(M-1)}|J, M-1\rangle$$

$$(\hat{J}_{1-} + \hat{J}_{2-})|j_1, m_1; j_2, m_2\rangle = \hbar\sqrt{j_1(j_1+1) - m_1(m_1 - 1)}|j_1, m_1 - 1; j_2, m_2\rangle + \hbar\sqrt{j_2(j_2+1) - m_2(m_2-1)}|j_1, m_1; j_2, m_2 - 1\rangle$$

Matching coefficients generates all CG coefficients for the $J = j_1 + j_2$ multiplet.

3. Use orthogonality for lower $J$ values. The state $|J = j_1 + j_2 - 1, M = j_1 + j_2 - 1\rangle$ is the state in the two-dimensional $M = j_1 + j_2 - 1$ subspace that is orthogonal to the state $|j_1 + j_2, j_1 + j_2 - 1\rangle$ obtained in step 2. Find it by Gram-Schmidt orthogonalization, then apply $\hat{J}_-$ to fill out the $J = j_1 + j_2 - 1$ multiplet.

4. Repeat for each successively lower $J$ until $J = |j_1 - j_2|$.

The General Recursion Relation

Applying $\hat{J}_-$ to $|J, M\rangle = \sum_{m_1, m_2} \langle j_1, m_1; j_2, m_2 | J, M\rangle |j_1, m_1; j_2, m_2\rangle$ produces the recursion relation:

$$\sqrt{J(J+1) - M(M-1)} \langle j_1, m_1; j_2, m_2 | J, M-1\rangle$$

$$= \sqrt{j_1(j_1+1) - m_1(m_1+1)} \langle j_1, m_1+1; j_2, m_2 | J, M\rangle$$

$$+ \sqrt{j_2(j_2+1) - m_2(m_2+1)} \langle j_1, m_1; j_2, m_2+1 | J, M\rangle$$

where $M = m_1 + m_2 + 1$ on the left side and $M = (m_1 + 1) + m_2 = m_1 + (m_2 + 1)$ on the right side. This relation, combined with the normalization condition and starting values, determines all CG coefficients.

Wigner 3j Symbols

For applications in atomic and nuclear physics, it is often more convenient to use the Wigner 3j symbol, which has better symmetry properties than the CG coefficient:

$$\begin{pmatrix} j_1 & j_2 & J \\ m_1 & m_2 & M \end{pmatrix}$$

The relationship to CG coefficients is:

$$\langle j_1, m_1; j_2, m_2 | J, M\rangle = (-1)^{-j_1 + j_2 - M}\sqrt{2J + 1} \begin{pmatrix} j_1 & j_2 & J \\ m_1 & m_2 & -M \end{pmatrix}$$

The 3j symbol vanishes unless: - $m_1 + m_2 + M = 0$ (note the sign convention: $-M$ in the CG definition) - The triangle rule is satisfied: $|j_1 - j_2| \leq J \leq j_1 + j_2$ - Each $|m_i| \leq j_i$

The 3j symbol has elegant symmetry properties under permutation of columns and sign reversal of all $m$ values, which make it the preferred tool in advanced applications.

Key Properties of CG Coefficients

1. Reality: Under the Condon-Shortley convention, all CG coefficients are real.

2. Vanishing condition: $\langle j_1, m_1; j_2, m_2 | J, M\rangle = 0$ unless $M = m_1 + m_2$ and $|j_1 - j_2| \leq J \leq j_1 + j_2$.

3. Orthogonality (completeness over uncoupled states):

$$\sum_{m_1, m_2} \langle j_1, m_1; j_2, m_2 | J, M\rangle \langle j_1, m_1; j_2, m_2 | J', M'\rangle = \delta_{JJ'}\delta_{MM'}$$

4. Orthogonality (completeness over coupled states):

$$\sum_{J, M} \langle j_1, m_1; j_2, m_2 | J, M\rangle \langle j_1, m_1'; j_2, m_2' | J, M\rangle = \delta_{m_1 m_1'}\delta_{m_2 m_2'}$$

5. Symmetry relations:

$$\langle j_1, m_1; j_2, m_2 | J, M\rangle = (-1)^{j_1 + j_2 - J}\langle j_2, m_2; j_1, m_1 | J, M\rangle$$

$$\langle j_1, m_1; j_2, m_2 | J, M\rangle = (-1)^{j_1 + j_2 - J}\langle j_1, -m_1; j_2, -m_2 | J, -M\rangle$$

🔴 Warning: The phase factor $(-1)^{j_1 + j_2 - J}$ in the symmetry relations is the source of many sign errors. Always check your phase conventions when using CG coefficients from different sources. The Condon-Shortley convention is standard in most modern textbooks (Griffiths, Sakurai, Shankar), but older references and some nuclear physics literature may differ.

14.5 Worked Example: Two Spin-$\frac{1}{2}$ Particles — Singlet and Triplet

This is the simplest nontrivial example and one of the most important in all of physics. It describes the total spin state of two electrons, two nucleons, or any pair of spin-1/2 particles.

Setup

We have $j_1 = j_2 = 1/2$. The uncoupled basis has $(2)(2) = 4$ states:

$$|\!\uparrow\uparrow\rangle \equiv |{\textstyle\frac{1}{2}}, {\textstyle\frac{1}{2}}; {\textstyle\frac{1}{2}}, {\textstyle\frac{1}{2}}\rangle, \quad |\!\uparrow\downarrow\rangle, \quad |\!\downarrow\uparrow\rangle, \quad |\!\downarrow\downarrow\rangle$$

where $|\!\uparrow\downarrow\rangle \equiv |{\textstyle\frac{1}{2}}, {\textstyle\frac{1}{2}}; {\textstyle\frac{1}{2}}, -{\textstyle\frac{1}{2}}\rangle$, etc.

The triangle rule gives $J = 1/2 + 1/2 = 1$ and $J = 1/2 - 1/2 = 0$. So the allowed values are $J = 1$ (three states: the triplet) and $J = 0$ (one state: the singlet). Dimension check: $3 + 1 = 4$ ✓.

Building the Triplet ($J = 1$)

Step 1: The stretched state is unique:

$$|1, 1\rangle = |\!\uparrow\uparrow\rangle$$

CG coefficient: $\langle {\textstyle\frac{1}{2}}, {\textstyle\frac{1}{2}}; {\textstyle\frac{1}{2}}, {\textstyle\frac{1}{2}} | 1, 1\rangle = 1$.

Step 2: Apply $\hat{J}_- = \hat{S}_{1-} + \hat{S}_{2-}$ to both sides.

Left side: $$\hat{J}_-|1, 1\rangle = \hbar\sqrt{1(2) - 1(0)}|1, 0\rangle = \hbar\sqrt{2}|1, 0\rangle$$

Right side: $$(\hat{S}_{1-} + \hat{S}_{2-})|\!\uparrow\uparrow\rangle = \hbar\sqrt{\tfrac{1}{2}\cdot\tfrac{3}{2} - \tfrac{1}{2}(-\tfrac{1}{2})}|\!\downarrow\uparrow\rangle + \hbar\sqrt{\tfrac{1}{2}\cdot\tfrac{3}{2} - \tfrac{1}{2}(-\tfrac{1}{2})}|\!\uparrow\downarrow\rangle = \hbar(|\!\downarrow\uparrow\rangle + |\!\uparrow\downarrow\rangle)$$

Setting left = right:

$$\hbar\sqrt{2}|1, 0\rangle = \hbar(|\!\uparrow\downarrow\rangle + |\!\downarrow\uparrow\rangle)$$

$$\boxed{|1, 0\rangle = \frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle + |\!\downarrow\uparrow\rangle)}$$

Step 3: Apply $\hat{J}_-$ again:

$$\hat{J}_-|1, 0\rangle = \hbar\sqrt{2}|1, -1\rangle$$

Acting on the right side: $$(\hat{S}_{1-} + \hat{S}_{2-})\frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle + |\!\downarrow\uparrow\rangle) = \frac{1}{\sqrt{2}}(|\!\downarrow\downarrow\rangle + |\!\downarrow\downarrow\rangle) = \frac{2}{\sqrt{2}}|\!\downarrow\downarrow\rangle = \sqrt{2}\hbar|\!\downarrow\downarrow\rangle$$

So:

$$\boxed{|1, -1\rangle = |\!\downarrow\downarrow\rangle}$$

Building the Singlet ($J = 0$)

The singlet $|0, 0\rangle$ lives in the $M = 0$ subspace, which is spanned by $|\!\uparrow\downarrow\rangle$ and $|\!\downarrow\uparrow\rangle$. It must be orthogonal to $|1, 0\rangle = \frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle + |\!\downarrow\uparrow\rangle)$.

The orthogonal combination (with the Condon-Shortley phase convention) is:

$$\boxed{|0, 0\rangle = \frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle - |\!\downarrow\uparrow\rangle)}$$

Summary: The Complete CG Table for $\frac{1}{2} \otimes \frac{1}{2}$

💡 Key Insight: The triplet states are symmetric under particle exchange ($|\!\uparrow\downarrow\rangle + |\!\downarrow\uparrow\rangle$), while the singlet is antisymmetric ($|\!\uparrow\downarrow\rangle - |\!\downarrow\uparrow\rangle$). This is not a coincidence. The triplet has $S = 1$ (symmetric under exchange of two spin-1/2 particles), and the singlet has $S = 0$ (antisymmetric). This symmetry property has profound consequences for identical particles, which you will explore in Chapter 15.

🔗 Spaced Review — Chapter 11: If you studied tensor products in Chapter 11, you may recognize the singlet state $\frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle - |\!\downarrow\uparrow\rangle)$ as one of the four Bell states. The singlet is maximally entangled — measuring one spin instantly determines the other. This connection between angular momentum coupling and entanglement runs deep and will resurface in Chapter 24 (Bell inequalities).

🧪 Experiment: The singlet-triplet structure of two-electron spin states is directly observable. In helium, the singlet (parahelium, $S = 0$) and triplet (orthohelium, $S = 1$) states have dramatically different energy spectra because the spin symmetry determines the spatial symmetry through the Pauli exclusion principle (Chapter 15). Historically, parahelium and orthohelium were thought to be different elements!

Verification: $\hat{\mathbf{S}}^2$ Eigenvalues

We can verify our results by computing $\hat{\mathbf{S}}^2 = (\hat{\mathbf{S}}_1 + \hat{\mathbf{S}}_2)^2 = \hat{S}_1^2 + \hat{S}_2^2 + 2\hat{\mathbf{S}}_1 \cdot \hat{\mathbf{S}}_2$. Using $\hat{\mathbf{S}}_1 \cdot \hat{\mathbf{S}}_2 = \hat{S}_{1z}\hat{S}_{2z} + \frac{1}{2}(\hat{S}_{1+}\hat{S}_{2-} + \hat{S}_{1-}\hat{S}_{2+})$:

For the triplet $|1, 0\rangle$:

$$\hat{\mathbf{S}}_1 \cdot \hat{\mathbf{S}}_2 \frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle + |\!\downarrow\uparrow\rangle) = \frac{1}{\sqrt{2}}\left[-\frac{\hbar^2}{4}(|\!\uparrow\downarrow\rangle + |\!\downarrow\uparrow\rangle) + \frac{\hbar^2}{2}(|\!\downarrow\uparrow\rangle + |\!\uparrow\downarrow\rangle)\right] = \frac{\hbar^2}{4} \cdot \frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle + |\!\downarrow\uparrow\rangle)$$

So $\hat{\mathbf{S}}_1 \cdot \hat{\mathbf{S}}_2 = \frac{\hbar^2}{4}$ for the triplet, giving $\hat{S}^2 = \frac{3}{4}\hbar^2 + \frac{3}{4}\hbar^2 + 2 \cdot \frac{\hbar^2}{4} = 2\hbar^2 = 1(1+1)\hbar^2$. ✓

For the singlet $|0, 0\rangle$:

$$\hat{\mathbf{S}}_1 \cdot \hat{\mathbf{S}}_2 \frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle - |\!\downarrow\uparrow\rangle) = \frac{1}{\sqrt{2}}\left[-\frac{\hbar^2}{4}(|\!\uparrow\downarrow\rangle - |\!\downarrow\uparrow\rangle) + \frac{\hbar^2}{2}(|\!\downarrow\uparrow\rangle - |\!\uparrow\downarrow\rangle)\right] = -\frac{3\hbar^2}{4} \cdot \frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle - |\!\downarrow\uparrow\rangle)$$

So $\hat{\mathbf{S}}_1 \cdot \hat{\mathbf{S}}_2 = -\frac{3\hbar^2}{4}$ for the singlet, giving $\hat{S}^2 = \frac{3}{4}\hbar^2 + \frac{3}{4}\hbar^2 + 2 \cdot (-\frac{3\hbar^2}{4}) = 0 = 0(0+1)\hbar^2$. ✓

💡 Key Insight: A useful identity from this calculation: $\hat{\mathbf{S}}_1 \cdot \hat{\mathbf{S}}_2 = \frac{1}{2}[\hat{S}^2 - \hat{S}_1^2 - \hat{S}_2^2] = \frac{\hbar^2}{2}[S(S+1) - s_1(s_1+1) - s_2(s_2+1)]$. For two spin-1/2 particles, this gives $+\hbar^2/4$ for the triplet and $-3\hbar^2/4$ for the singlet. This identity appears constantly in spin-orbit coupling, exchange interactions, and NMR.

14.6 Worked Example: Spin-Orbit Coupling in Hydrogen

The Physical Problem

In hydrogen, the electron has orbital angular momentum $\hat{\mathbf{L}}$ (with quantum number $\ell$) and spin $\hat{\mathbf{S}}$ (with $s = 1/2$). The spin-orbit Hamiltonian is:

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

where $V(r)$ is the Coulomb potential. This interaction couples $\hat{\mathbf{L}}$ and $\hat{\mathbf{S}}$, and the natural basis is the coupled basis $|j, m_j\rangle$ where $\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}$.

🔗 Spaced Review — Chapter 5: The hydrogen atom quantum numbers $n, \ell, m_\ell$ came from solving the Schrodinger equation with the Coulomb potential. The energy depended only on $n$, with the $\ell$-degeneracy being a special feature of the $1/r$ potential. Spin-orbit coupling breaks this degeneracy, as we will see.

Allowed $J$ Values

For a given $\ell \geq 1$, coupling $\hat{\mathbf{L}}$ ($j_1 = \ell$) with $\hat{\mathbf{S}}$ ($j_2 = 1/2$) gives:

$$j = \ell + \tfrac{1}{2} \quad \text{or} \quad j = \ell - \tfrac{1}{2}$$

For $\ell = 0$, only $j = 1/2$ is possible (the triangle rule gives $|0 - 1/2| = 1/2$ to $0 + 1/2 = 1/2$).

Spectroscopic notation: These states are labeled $n\ell_j$. For the $n = 2$ level: - $\ell = 0$: $j = 1/2$ only → $2s_{1/2}$ - $\ell = 1$: $j = 3/2$ or $j = 1/2$ → $2p_{3/2}$ and $2p_{1/2}$

The spin-orbit coupling splits $2p$ into $2p_{3/2}$ and $2p_{1/2}$, with an energy difference that explains the fine structure of the hydrogen spectrum.

Coupled States for $\ell = 1$, $s = 1/2$

We couple $j_1 = 1$ and $j_2 = 1/2$ to get $J = 3/2$ and $J = 1/2$. The uncoupled basis has $3 \times 2 = 6$ states; the coupled basis has $4 + 2 = 6$ states.

The $J = 3/2$ multiplet:

The stretched state: $$|{\textstyle\frac{3}{2}}, {\textstyle\frac{3}{2}}\rangle = |1, 1\rangle|\!\uparrow\rangle$$

Applying $\hat{J}_-$: $$\hat{J}_-|{\textstyle\frac{3}{2}}, {\textstyle\frac{3}{2}}\rangle = \hbar\sqrt{3}|{\textstyle\frac{3}{2}}, {\textstyle\frac{1}{2}}\rangle$$

$$(\hat{L}_- + \hat{S}_-)|1, 1\rangle|\!\uparrow\rangle = \hbar\sqrt{2}|1, 0\rangle|\!\uparrow\rangle + \hbar|1, 1\rangle|\!\downarrow\rangle$$

Therefore: $$|{\textstyle\frac{3}{2}}, {\textstyle\frac{1}{2}}\rangle = \sqrt{\frac{2}{3}}|1, 0\rangle|\!\uparrow\rangle + \sqrt{\frac{1}{3}}|1, 1\rangle|\!\downarrow\rangle$$

Continuing with $\hat{J}_-$: $$|{\textstyle\frac{3}{2}}, -{\textstyle\frac{1}{2}}\rangle = \sqrt{\frac{1}{3}}|1, -1\rangle|\!\uparrow\rangle + \sqrt{\frac{2}{3}}|1, 0\rangle|\!\downarrow\rangle$$

$$|{\textstyle\frac{3}{2}}, -{\textstyle\frac{3}{2}}\rangle = |1, -1\rangle|\!\downarrow\rangle$$

The $J = 1/2$ multiplet (orthogonal to the above in each $M$ subspace):

$$|{\textstyle\frac{1}{2}}, {\textstyle\frac{1}{2}}\rangle = -\sqrt{\frac{1}{3}}|1, 0\rangle|\!\uparrow\rangle + \sqrt{\frac{2}{3}}|1, 1\rangle|\!\downarrow\rangle$$

$$|{\textstyle\frac{1}{2}}, -{\textstyle\frac{1}{2}}\rangle = -\sqrt{\frac{2}{3}}|1, -1\rangle|\!\uparrow\rangle + \sqrt{\frac{1}{3}}|1, 0\rangle|\!\downarrow\rangle$$

The CG Table for $1 \otimes \frac{1}{2}$

🔴 Warning: When reading CG tables, always note the sign conventions carefully. The negative signs in the $J = 1/2$ states arise from the Condon-Shortley phase convention and from orthogonality to the $J = 3/2$ states. These signs have physical consequences — they affect transition amplitudes and selection rules.

Energy Shifts from Spin-Orbit Coupling

Using the identity $\hat{\mathbf{L}} \cdot \hat{\mathbf{S}} = \frac{1}{2}(\hat{J}^2 - \hat{L}^2 - \hat{S}^2)$, the expectation value in a coupled state $|n, \ell, j, m_j\rangle$ is:

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

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

The energy splitting is proportional to these values. The $j = 3/2$ level is higher in energy, and the $j = 1/2$ level is lower, with a splitting ratio of $2:1$ weighted by the degeneracies $4:2$. This is the fine structure of hydrogen's $p$ levels.

📊 By the Numbers: For the $2p$ level in hydrogen, the spin-orbit splitting is approximately $4.5 \times 10^{-5}$ eV, corresponding to the famous Lamb shift experiments. The $2p_{3/2}$ state lies above $2p_{1/2}$ by this amount. This tiny energy difference was one of the first precision tests of quantum electrodynamics (QED) and earned Willis Lamb the 1955 Nobel Prize.

L-S Coupling vs. j-j Coupling: Two Philosophies of Angular Momentum Addition

When an atom has multiple electrons, each with orbital angular momentum $\hat{\boldsymbol{\ell}}_i$ and spin $\hat{\mathbf{s}}_i$, the order in which we couple these angular momenta matters physically — it determines which quantum numbers are approximately good.

L-S coupling (Russell-Saunders coupling) applies when the electrostatic (Coulomb) interaction between electrons is much stronger than the spin-orbit coupling. This is typical for light atoms ($Z \lesssim 30$). The procedure is:

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

States are labeled by term symbols $^{2S+1}L_J$. For example, the ground state of carbon has two $2p$ electrons coupling to $L = 1$, $S = 1$, $J = 0$, giving the term symbol $^3P_0$. The notation uses spectroscopic letters for $L$: $S, P, D, F, G, H, \ldots$ for $L = 0, 1, 2, 3, 4, 5, \ldots$

The physical justification for L-S coupling is that in light atoms, the residual Coulomb repulsion between electrons (the part not captured by the central field approximation) is typically $10^4$--$10^5$ cm$^{-1}$, while the spin-orbit coupling is only $10$--$10^2$ cm$^{-1}$. Since the electrostatic interaction depends on the spatial configuration of the electrons (which depends on $L$) but not on spin, the appropriate first step is to couple the spatial degrees of freedom together, then the spin degrees of freedom together, and only then couple $L$ to $S$.

j-j coupling applies when the spin-orbit coupling for individual electrons is comparable to or exceeds the residual Coulomb interaction. This occurs for heavy atoms ($Z \gtrsim 70$), where the spin-orbit interaction scales as $Z^4$ while the Coulomb repulsion scales more slowly. The procedure is:

1. For each electron, couple its own $\hat{\boldsymbol{\ell}}_i$ and $\hat{\mathbf{s}}_i$: $\hat{\mathbf{j}}_i = \hat{\boldsymbol{\ell}}_i + \hat{\mathbf{s}}_i$, with quantum number $j_i$.
2. Couple all individual $\hat{\mathbf{j}}_i$: $\hat{\mathbf{J}} = \sum_i \hat{\mathbf{j}}_i$.

States are labeled by the individual $j_i$ values: $(j_1, j_2, \ldots)_J$.

For most elements in the periodic table, neither pure L-S nor pure j-j coupling is exactly correct. The true coupling is intermediate, and the CG coefficient formalism provides the mathematical tool for transforming between the two limiting descriptions. The Wigner 6j symbol (Section 14.8) connects the L-S and j-j coupling bases.

⚖️ Interpretation: The distinction between L-S and j-j coupling is not just a mathematical convenience — it determines which quantum numbers are "good" (nearly conserved). In L-S coupling, $L$ and $S$ are approximately good quantum numbers because the electrostatic interaction conserves them. In j-j coupling, the individual $j_i$ are approximately good because the spin-orbit interaction conserves them. The choice between coupling schemes is a physical statement about the relative strengths of the interactions.

The Lande Interval Rule

In L-S coupling, the spin-orbit interaction produces a fine structure splitting that follows a simple pattern. For a term $^{2S+1}L_J$ with fixed $L$ and $S$, the energy splitting between adjacent $J$ levels satisfies:

$$E(J) - E(J - 1) = A \cdot J$$

where $A$ is the spin-orbit coupling constant (which depends on the radial wavefunction but not on $J$). This Lande interval rule is a direct consequence of the $\hat{\mathbf{L}} \cdot \hat{\mathbf{S}}$ form of the coupling and the eigenvalue formula:

$$\langle \hat{\mathbf{L}} \cdot \hat{\mathbf{S}} \rangle_J = \frac{1}{2}[J(J+1) - L(L+1) - S(S+1)]$$

Taking the difference between $J$ and $J-1$:

$$\langle \hat{\mathbf{L}} \cdot \hat{\mathbf{S}} \rangle_J - \langle \hat{\mathbf{L}} \cdot \hat{\mathbf{S}} \rangle_{J-1} = \frac{1}{2}[J(J+1) - (J-1)J] = J$$

The interval rule is a powerful experimental diagnostic: if observed energy splittings obey the rule, L-S coupling is a good description. Departures from the interval rule signal the breakdown of pure L-S coupling and the onset of intermediate or j-j coupling.

🧪 Experiment: The classic test of the Lande interval rule uses the fine structure of the calcium $4s4p \; ^3P$ term ($L = 1$, $S = 1$, $J = 0, 1, 2$). The measured splittings are $E(J=1) - E(J=0) = 52.1$ cm$^{-1}$ and $E(J=2) - E(J=1) = 105.0$ cm$^{-1}$. The ratio is $105.0/52.1 = 2.02$, very close to the predicted ratio of $2/1 = 2.0$. The small deviation arises from second-order spin-orbit effects.

14.7 The Wigner-Eckart Theorem

The Wigner-Eckart theorem is one of the most powerful results in quantum mechanics. It dramatically simplifies the calculation of matrix elements of operators that transform in a definite way under rotations.

Irreducible Tensor Operators: Definition

An irreducible tensor operator of rank $k$ is a set of $2k + 1$ operators $\hat{T}_q^{(k)}$ (with $q = -k, -k+1, \ldots, k$) that transform under rotations exactly like the angular momentum states $|k, q\rangle$. Formally, they satisfy the commutation relations:

$$[\hat{J}_z, \hat{T}_q^{(k)}] = q\hbar \hat{T}_q^{(k)}$$

$$[\hat{J}_\pm, \hat{T}_q^{(k)}] = \hbar\sqrt{k(k+1) - q(q \pm 1)} \hat{T}_{q\pm 1}^{(k)}$$

These commutation relations are the defining property of an irreducible tensor operator. They say that the components $\hat{T}_q^{(k)}$ transform among themselves under rotations in exactly the same way as the states $|k, q\rangle$. The operator "carries" angular momentum $k$ and projection $q$.

The language of irreducible tensor operators may seem abstract at first, but it provides a unified framework for classifying all physically relevant operators by their angular momentum content. Once you identify the rank of an operator, the Wigner-Eckart theorem immediately tells you everything about the $m$-dependence of its matrix elements. The key examples:

Rank 0 (scalar): A single operator $\hat{T}_0^{(0)}$ that commutes with all components of $\hat{\mathbf{J}}$. The Hamiltonian is a scalar operator if it has rotational symmetry.

Rank 1 (vector): Three operators $\hat{T}_{-1}^{(1)}, \hat{T}_0^{(1)}, \hat{T}_{+1}^{(1)}$ that transform like a vector under rotations. The position operator $\hat{\mathbf{r}}$ and the angular momentum $\hat{\mathbf{J}}$ itself are rank-1 tensor operators. The spherical components of $\hat{\mathbf{r}}$ are:

$$\hat{T}_0^{(1)} = \hat{z} = r\cos\theta, \qquad \hat{T}_{\pm 1}^{(1)} = \mp\frac{1}{\sqrt{2}}(\hat{x} \pm i\hat{y}) = \mp\frac{r}{\sqrt{2}}\sin\theta \, e^{\pm i\phi}$$

Rank 2 (quadrupole): Five operators that transform like the $\ell = 2$ spherical harmonics. The electric quadrupole moment operator is rank 2.

Statement of the Wigner-Eckart Theorem

The Wigner-Eckart theorem states that the matrix element of an irreducible tensor operator $\hat{T}_q^{(k)}$ between angular momentum eigenstates factors into a geometric part (a CG coefficient that depends on $m$, $q$, $M$) and a dynamical part (a reduced matrix element that is independent of all magnetic quantum numbers):

$$\langle j', m' | \hat{T}_q^{(k)} | j, m\rangle = \frac{\langle j' \| \hat{T}^{(k)} \| j \rangle}{\sqrt{2j' + 1}} \langle j, m; k, q | j', m'\rangle$$

Here: - The left side is the full matrix element we want to compute. - $\langle j, m; k, q | j', m'\rangle$ is a Clebsch-Gordan coefficient — we just learned how to calculate these! - $\langle j' \| \hat{T}^{(k)} \| j \rangle$ is the reduced matrix element (or double-bar matrix element). It depends on $j'$, $k$, and $j$, but not on $m'$, $q$, or $m$. It contains all the dynamical information — the physics of the specific operator and states.

💡 Key Insight: The Wigner-Eckart theorem says that the $m$-dependence of every matrix element of an irreducible tensor operator is completely determined by geometry (the CG coefficient). The physics — the hard part — is compressed into the single reduced matrix element $\langle j' \| T^{(k)} \| j \rangle$. Instead of calculating $(2j'+1)(2k+1)(2j+1)$ separate matrix elements, you need only one reduced matrix element. The CG coefficient handles all the rest.

Why It Works: The Intuition

The Wigner-Eckart theorem is a consequence of rotational symmetry. The operator $\hat{T}_q^{(k)}$ carries angular momentum $(k, q)$, just as a state $|k, q\rangle$ does. When we compute the matrix element $\langle j', m' | \hat{T}_q^{(k)} | j, m\rangle$, we are "coupling" the angular momentum $(j, m)$ of the initial state with the angular momentum $(k, q)$ of the operator to produce the angular momentum $(j', m')$ of the final state.

This is exactly the angular momentum addition problem! The geometric factor — how the angular momenta combine — is given by the CG coefficient. The dynamical factor — the strength of the coupling — is the reduced matrix element.

Convention Warning

🔴 Warning: There are two common conventions for the Wigner-Eckart theorem. The convention used here (and in Griffiths, Sakurai, and most modern physics texts) places $1/\sqrt{2j'+1}$ with the reduced matrix element. Some sources (particularly in nuclear physics, following Edmonds) write the theorem without this factor, absorbing it into the definition of the reduced matrix element. Always check which convention a source uses before combining results from different references.

The version in terms of the Wigner 3j symbol is:

$$\langle j', m' | \hat{T}_q^{(k)} | j, m\rangle = (-1)^{j'-m'} \begin{pmatrix} j' & k & j \\ -m' & q & m \end{pmatrix} \langle j' \| \hat{T}^{(k)} \| j \rangle$$

A Concrete Application: Dipole Matrix Elements

Consider the electric dipole operator $\hat{\mathbf{d}} = e\hat{\mathbf{r}}$, which is a rank-1 tensor operator. For hydrogen, we want the matrix elements:

$$\langle n', \ell', m' | er_q | n, \ell, m\rangle$$

where $r_q$ is the $q$-th spherical component of $\hat{\mathbf{r}}$ ($q = -1, 0, +1$).

By the Wigner-Eckart theorem:

$$\langle n', \ell', m' | er_q | n, \ell, m\rangle = \frac{\langle n', \ell' \| er \| n, \ell \rangle}{\sqrt{2\ell' + 1}} \langle \ell, m; 1, q | \ell', m'\rangle$$

The CG coefficient immediately gives us: 1. $m' = m + q$ (from the CG selection rule $M = m_1 + m_2$) 2. $|\ell - 1| \leq \ell' \leq \ell + 1$ (from the triangle rule with $k = 1$)

For the specific case of the electric dipole, there is an additional constraint from parity: $\ell'$ must differ from $\ell$ by an odd integer. Combined with the triangle rule, this gives:

$$\Delta\ell = \pm 1, \qquad \Delta m = 0, \pm 1$$

These are the electric dipole selection rules for atomic transitions. We have derived them with almost no work — the Wigner-Eckart theorem did the heavy lifting.

The reduced matrix element $\langle n', \ell' \| er \| n, \ell \rangle$ involves a radial integral:

$$\langle n', \ell' \| er \| n, \ell \rangle = e(-1)^{\ell'}\sqrt{(2\ell'+1)(2\ell+1)} \begin{pmatrix} \ell' & 1 & \ell \\ 0 & 0 & 0 \end{pmatrix} \int_0^\infty R_{n'\ell'}(r) \, r \, R_{n\ell}(r) \, r^2 \, dr$$

where $R_{n\ell}(r)$ are the hydrogen radial wavefunctions from Chapter 5. The radial integral contains the dynamical information; the 3j symbol contains the geometry.

✅ Checkpoint: Can you explain in words why the Wigner-Eckart theorem is useful? Here is the test: suppose you need the matrix elements $\langle 2, m' | \hat{T}_q^{(1)} | 3, m\rangle$ for all $m, m', q$. Without the theorem, you have $(5)(3)(7) = 105$ matrix elements to calculate. With the theorem, how many reduced matrix elements do you need? (Answer: just one — $\langle 2 \| T^{(1)} \| 3 \rangle$.)

14.8 Irreducible Tensor Operators and Selection Rules

Constructing Tensor Operators

Given two vector operators $\hat{\mathbf{A}}$ and $\hat{\mathbf{B}}$ (each a rank-1 tensor), we can construct tensor operators of ranks 0, 1, and 2 from their product — just as coupling $j_1 = 1$ with $j_2 = 1$ gives $J = 0, 1, 2$:

Rank 0 (scalar product): $$[\hat{\mathbf{A}} \otimes \hat{\mathbf{B}}]_0^{(0)} = -\frac{1}{\sqrt{3}}\hat{\mathbf{A}} \cdot \hat{\mathbf{B}}$$

Rank 1 (vector product): $$[\hat{\mathbf{A}} \otimes \hat{\mathbf{B}}]_q^{(1)} \propto (\hat{\mathbf{A}} \times \hat{\mathbf{B}})_q$$

Rank 2 (quadrupole): $$[\hat{\mathbf{A}} \otimes \hat{\mathbf{B}}]_q^{(2)} = \sum_{q_1, q_2} \langle 1, q_1; 1, q_2 | 2, q\rangle A_{q_1} B_{q_2}$$

This construction — building higher-rank tensors from lower-rank ones using CG coefficients — is ubiquitous in nuclear and atomic physics. It is the operator analog of coupling angular momenta.

Selection Rules from the Wigner-Eckart Theorem

The Wigner-Eckart theorem provides selection rules through the CG coefficient $\langle j, m; k, q | j', m'\rangle$. This coefficient vanishes unless:

1. $m$-selection rule: $m' = m + q$
2. Triangle rule: $|j - k| \leq j' \leq j + k$

For specific types of operators:

Electric dipole ($k = 1$): $\Delta j = 0, \pm 1$ (but $j = 0 \to j' = 0$ is forbidden since the CG coefficient vanishes), $\Delta m = 0, \pm 1$. Combined with parity, $\Delta \ell = \pm 1$.

Electric quadrupole ($k = 2$): $\Delta j = 0, \pm 1, \pm 2$ (but $j = 0 \to j' = 0$ and $j = 1/2 \to j' = 1/2$ are forbidden), $\Delta m = 0, \pm 1, \pm 2$. Combined with parity, $\Delta \ell = 0, \pm 2$.

Magnetic dipole ($k = 1$, no parity change): $\Delta j = 0, \pm 1$, $\Delta m = 0, \pm 1$, $\Delta \ell = 0$.

💡 Key Insight: Selection rules are not arbitrary catalog entries to be memorized. They are theorems that follow directly from the angular momentum properties of the relevant operator and the Wigner-Eckart theorem. If you know the rank of the operator and the parity, you can derive the selection rules from scratch.

The Projection Theorem (Special Case of Wigner-Eckart)

An important special case occurs for vector operators ($k = 1$) when $j' = j$. The Wigner-Eckart theorem then implies the projection theorem:

$$\langle j, m | \hat{\mathbf{V}} | j, m\rangle = \frac{\langle j, m | \hat{\mathbf{J}} \cdot \hat{\mathbf{V}} | j, m\rangle}{j(j+1)\hbar^2} \langle j, m | \hat{\mathbf{J}} | j, m\rangle$$

In words: the expectation value of any vector operator $\hat{\mathbf{V}}$ within a fixed-$j$ subspace is proportional to the expectation value of $\hat{\mathbf{J}}$. The proportionality constant involves $\hat{\mathbf{J}} \cdot \hat{\mathbf{V}}$.

This explains, for example, why the magnetic moment of an atom in a state $|j, m_j\rangle$ points along $\hat{\mathbf{J}}$, even though it is built from $\hat{\mathbf{L}}$ and $\hat{\mathbf{S}}$ which point in "different directions." The projection theorem forces:

$$\langle \hat{\boldsymbol{\mu}} \rangle = -g_j \mu_B \frac{\langle \hat{\mathbf{J}} \rangle}{\hbar}$$

where $g_j$ is the Lande g-factor:

$$g_j = 1 + \frac{j(j+1) + s(s+1) - \ell(\ell+1)}{2j(j+1)}$$

This result, which governs the anomalous Zeeman effect, follows directly from the Wigner-Eckart theorem applied to $\hat{\boldsymbol{\mu}} = -\mu_B(\hat{\mathbf{L}} + 2\hat{\mathbf{S}})/\hbar$.

🔗 Spaced Review — Chapter 13: Recall from Chapter 13 that the electron spin g-factor is $g_s \approx 2$ (the anomalous magnetic moment). The Lande g-factor $g_j$ interpolates between $g_\ell = 1$ (pure orbital) and $g_s \approx 2$ (pure spin), weighted by how much of each angular momentum contributes to the total $\hat{\mathbf{J}}$. When $j = \ell$ (pure orbital), $g_j = 1$. When $j = s$ (pure spin), $g_j = 2$.

Parity and Selection Rules

Selection rules from the Wigner-Eckart theorem constrain the angular momentum quantum numbers. But there is an additional constraint from parity — the behavior of the operator under spatial inversion $\hat{\mathbf{r}} \to -\hat{\mathbf{r}}$.

An operator has even parity if it is unchanged under inversion (like $\hat{r}^2$ or $\hat{\mathbf{L}}$), and odd parity if it changes sign (like $\hat{\mathbf{r}}$ or $\hat{\mathbf{p}}$). The parity of hydrogen states is $(-1)^\ell$, so:

• E1 (electric dipole): $\hat{\mathbf{r}}$ has odd parity → parity must change → $\Delta\ell = \pm 1, \pm 3, \ldots$ Combined with the triangle rule ($\Delta\ell = 0, \pm 1$ for a rank-1 operator), this gives $\Delta\ell = \pm 1$ only.

• M1 (magnetic dipole): $\hat{\mathbf{L}}$ and $\hat{\mathbf{S}}$ have even parity → parity must not change → $\Delta\ell = 0$. (Combined with the triangle rule, $\Delta\ell = 0$ is the only allowed value.)

• E2 (electric quadrupole): $r^2 Y_2^q$ has even parity → parity must not change → $\Delta\ell = 0, \pm 2$. (From the triangle rule with $k = 2$, $\Delta\ell = 0, \pm 1, \pm 2$; the parity constraint eliminates $\Delta\ell = \pm 1$.)

These results show how the angular momentum selection rules from the Wigner-Eckart theorem combine with parity to give the complete set of selection rules. Neither constraint alone is sufficient — you need both.

⚠️ Common Misconception: Students sometimes think the selection rule $\Delta\ell = \pm 1$ for E1 transitions comes from the angular momentum alone (the triangle rule with $k = 1$). This is incorrect — the triangle rule alone gives $\Delta\ell = 0, \pm 1$. The elimination of $\Delta\ell = 0$ requires the parity argument. Without parity, the selection rule would be weaker.

Wigner 6j Symbols (Preview)

When coupling three angular momenta — for example, two electrons each with spin and orbital angular momentum — the CG coefficients are not sufficient. Instead, we need the Wigner 6j symbol (or Racah W-coefficient), which describes the recoupling of three angular momenta:

$$\begin{Bmatrix} j_1 & j_2 & j_{12} \\ j_3 & j & j_{23} \end{Bmatrix}$$

The 6j symbol appears when comparing different coupling schemes — for example, L-S coupling (first couple all $\ell$'s, then all $s$'s, then couple $L$ and $S$) versus j-j coupling (first couple each $\ell_i$ with $s_i$ to get $j_i$, then couple the $j_i$'s).

We will not derive the 6j symbols here — they are covered in the advanced exercises and references — but you should know they exist and recognize their role as the next level of angular momentum coupling machinery after CG coefficients.

14.9 Summary and Project Checkpoint

What We Have Learned

This chapter has taken you from the simple question "how do two angular momenta combine?" to a powerful, general formalism that governs all matrix elements of operators with definite rotational properties. Along the way, we have encountered some of the most elegant mathematics in all of physics — mathematics that is not merely beautiful but indispensable for real calculations.

Let us recapitulate the logical arc. We began with two angular momenta living in a tensor product space, and asked how to build eigenstates of the total angular momentum from eigenstates of the individual momenta. The answer was the coupled basis, related to the uncoupled basis by the Clebsch-Gordan coefficients. We then used the Wigner-Eckart theorem to show that these same coefficients control the matrix elements of all operators with definite angular momentum properties. This deep connection between state coupling and operator matrix elements is what makes the formalism so powerful.

The main results:

1. The Coupling Problem: Two angular momenta $\hat{\mathbf{J}}_1$ (quantum number $j_1$) and $\hat{\mathbf{J}}_2$ ($j_2$) combine to form the total $\hat{\mathbf{J}} = \hat{\mathbf{J}}_1 + \hat{\mathbf{J}}_2$. The uncoupled basis $|j_1, m_1; j_2, m_2\rangle$ and coupled basis $|J, M\rangle$ are related by a unitary transformation.

2. The Triangle Rule: The allowed total angular momentum quantum numbers are $J = |j_1 - j_2|, |j_1 - j_2| + 1, \ldots, j_1 + j_2$, and $M = m_1 + m_2$.

3. Clebsch-Gordan Coefficients: $\langle j_1, m_1; j_2, m_2 | J, M\rangle$ are the matrix elements of the basis transformation. They are calculated by the recursion method (lowering operators + orthogonality).

4. The Wigner-Eckart Theorem: Matrix elements of irreducible tensor operators of rank $k$ factor as:

$$\langle j', m' | \hat{T}_q^{(k)} | j, m\rangle = \frac{\langle j' \| \hat{T}^{(k)} \| j \rangle}{\sqrt{2j' + 1}} \langle j, m; k, q | j', m'\rangle$$

5. Selection Rules: The CG coefficient in the Wigner-Eckart theorem immediately gives $\Delta m = q$ and $|j - k| \leq j' \leq j + k$. For electric dipole transitions ($k = 1$), $\Delta \ell = \pm 1$ and $\Delta m = 0, \pm 1$.

Uncoupled vs. Coupled: When to Use Each Basis

Looking Ahead

• Chapter 15 (Identical Particles): The singlet/triplet structure from Section 14.5 determines the symmetry of two-electron spin states, which in turn determines the spatial symmetry through the Pauli exclusion principle.
• Chapter 16 (Multi-Electron Atoms): CG coefficients and the coupling machinery are essential for constructing term symbols ($^{2S+1}L_J$) and understanding the periodic table.
• Chapter 18 (Fine Structure): The spin-orbit coupling from Section 14.6 will be treated quantitatively using degenerate perturbation theory, with the coupled basis as the natural framework.
• Chapter 21 (Time-Dependent Perturbation Theory): The selection rules from Section 14.8 determine which atomic transitions are allowed, connecting to spectroscopy and lasers.

Project Checkpoint: Toolkit v1.4 — The Clebsch-Gordan Module

Your Progressive Project checkpoint for this chapter builds the angular momentum coupling module for the Quantum Simulation Toolkit.

Required implementations:

1. clebsch_gordan(j1, m1, j2, m2, J, M) — Compute a single CG coefficient
2. coupled_basis(j1, j2) — Generate all coupled states as linear combinations of uncoupled states
3. uncoupled_to_coupled(j1, j2) — Build the full unitary transformation matrix
4. wigner_3j(j1, j2, j3, m1, m2, m3) — Compute Wigner 3j symbols

Verification tests: - Reproduce the $\frac{1}{2} \otimes \frac{1}{2}$ table from Section 14.5 - Reproduce the $1 \otimes \frac{1}{2}$ table from Section 14.6 - Verify orthogonality relations (Equations in Section 14.4) - Check unitarity of the transformation matrix

See code/project-checkpoint.py for the implementation and code/example-01-coupling.py for demonstration scripts.

Testing your implementation: Your CG calculator should reproduce every coefficient in the tables of Sections 14.5 and 14.6 to machine precision. As a stress test, try the coupling $j_1 = 5/2$, $j_2 = 3/2$ (a 24-dimensional space with $J = 1, 2, 3, 4$) and verify that the transformation matrix is unitary. A non-unitary matrix is a certain sign of a bug — CG coefficients always form a unitary transformation.

🔗 Connection: The CG coefficient calculator you build here will be used by the Chapter 16 term symbol engine (coupling electron angular momenta in multi-electron atoms), the Chapter 18 fine structure calculator (spin-orbit matrix elements), and the Chapter 21 transition probability module (selection rules and dipole matrix elements). Invest the time to get it right — it pays off throughout the rest of the textbook.

Chapter 14 Notation Summary