Coupling of Angular Momenta and Total J

Take a carbon atom in its ground configuration 2p22p^2. Ask a simple question: how many distinct energy levels does it actually have? The answer is not one, and not even three. It is five, spread over an electron-volt of energy by electrostatic repulsion. These levels are split further by relativistic spin-orbit effects. Untangling that structure requires combining the orbital angular momenta of the two electrons with each other. It also requires separately combining their spins. The resulting totals may then need to be combined again. This is the problem of angular momentum coupling. Given two (or more) angular momentum vectors, each individually quantized, what are the allowed magnitudes and quantum numbers of their vector sum? The problem matters because almost every observable involving more than one source of angular momentum is governed by the same addition rule. This includes two electron spins, an electron spin and its orbital motion, two nuclear spins, and a photon's angular momentum with an atom's. Such observables appear throughout atomic, molecular, and nuclear physics. The rule is deceptively simple to state and deeply consequential. Two angular momenta j1j_1 and j2j_2 combine into a total JJ. That total can take only the values ∣j1−j2∣,∣j1−j2∣+1,…,j1+j2|j_1-j_2|, |j_1-j_2|+1, \ldots, j_1+j_2, each appearing exactly once. This single fact explains why term symbols look the way they do. It explains why the periodic table's spectroscopic notation works. It explains why NMR spectra split into multiplets. And it explains why fine structure in atomic spectra has the specific number of lines it does. The big picture is a change of basis. There are two natural but inequivalent ways to describe a two-angular-momentum system. The first is the uncoupled basis ∣j1m1⟩∣j2m2⟩|j_1 m_1\rangle|j_2 m_2\rangle, which diagonalizes J12,J1z,J22,J2zJ_1^2, J_{1z}, J_2^2, J_{2z}. The second is the coupled basis ∣j1j2;JM⟩|j_1 j_2; J M\rangle, which diagonalizes J12,J22,J2,JzJ_1^2, J_2^2, J^2, J_z instead. Both are complete bases for the same (2j1+1)(2j2+1)(2j_1+1)(2j_2+1)-dimensional space. The Clebsch-Gordan coefficients are simply the entries of the unitary matrix connecting them. Whether the coupled or uncoupled description is "the right one" depends entirely on which operators commute with the Hamiltonian. That, in turn, is exactly the question that spin-orbit coupling and fine structure force you to confront. This article builds the coupling rule from the commutation algebra of angular momentum operators. It works through the Clebsch-Gordan decomposition by explicit vector-counting and by ladder operators. It also shows how coupled angular momenta produce the term symbols and multiplicities used throughout atomic and molecular spectroscopy.

Key takeaways
  • Two angular momenta j1,j2j_1, j_2 combine into total JJ with allowed values J=∣j1−j2∣,…,j1+j2J = |j_1-j_2|, \ldots, j_1+j_2 in integer steps — the triangle rule.
  • The coupled basis ∣J,M⟩|J,M\rangle and uncoupled basis ∣m1,m2⟩|m_1,m_2\rangle span the same space; Clebsch-Gordan coefficients ⟨j1m1j2m2∣JM⟩\langle j_1 m_1 j_2 m_2|JM\rangle connect them.
  • Dimension is conserved: (2j1+1)(2j2+1)=∑J(2J+1)(2j_1+1)(2j_2+1) = \sum_J (2J+1), which is the bookkeeping check on every coupling table.
  • Which basis is physical depends on the Hamiltonian: if HH commutes with J12,J22,J2,JzJ_1^2,J_2^2,J^2,J_z but not with J1z,J2zJ_{1z},J_{2z} separately, the coupled basis diagonalizes HH.
  • Term symbols 2S+1LJ^{2S+1}L_J for atoms arise from coupling orbital angular momenta to get LL, spins to get SS, then LL and SS to get JJ (Russell-Saunders/LS coupling).
  • Key idea: angular momentum coupling is nothing but a unitary change of basis dictated by which operators actually commute with the physical Hamiltonian.

Where this lives today

Angular momentum coupling is not a historical curiosity confined to term-symbol bookkeeping — it is the operating grammar behind entire branches of modern spectroscopy and technology. In atomic and molecular spectroscopy, term symbols 2S+1LJ^{2S+1}L_J generated by LS coupling remain the standard language for cataloging electronic states of atoms and small molecules, and the selection rules on ΔJ=0,±1\Delta J = 0,\pm1 for electric dipole transitions follow directly from the coupled-basis structure and the Wigner-Eckart theorem. Fine structure, which splits terms like 3P^3P of carbon or the sodium D-lines (2P1/2^2P_{1/2} and 2P3/2^2P_{3/2}) into distinct JJ levels, is the direct successor topic — Spin-Orbit Coupling and Fine Structure is precisely the physical mechanism (relativistic coupling of orbital and spin magnetic moments) that makes the coupled JJ basis the energetically preferred one, splitting what coupling theory shows to be otherwise degenerate levels. In nuclear magnetic resonance and EPR spectroscopy, the coupling of nuclear spins to each other (scalar JJ-coupling in NMR, hyperfine coupling in EPR) produces the characteristic multiplet splitting patterns — doublets, triplets, quartets — used every day to determine molecular structure; the intensity ratios within these multiplets are Clebsch-Gordan-derived binomial coefficients. In lanthanide and actinide chemistry, where spin-orbit coupling is comparable to or exceeds inter-electron repulsion, the same coupling algebra (now in jjjj or intermediate coupling schemes) determines the magnetic and optical properties exploited in MRI contrast agents (gadolinium complexes), lanthanide-based luminescent probes, and single-molecule magnets. In quantum information science, coupling of two qubits' angular-momentum-like degrees of freedom (e.g., coupling nuclear and electron spins in nitrogen-vacancy centers in diamond) uses precisely this formalism to design controlled quantum gates. And in nuclear and particle physics, the same triangle rule governs the coupling of nucleon spins and isospins, angular momentum conservation in nuclear reactions, and the classification of hadron multiplets in the quark model — all descendants of the identical algebraic structure first worked out for two orbiting electrons.

Ask the tutor: How does adding spin-orbit coupling as a perturbation actually split the 3P^3P term of carbon into its J=0,1,2J=0,1,2 fine-structure levels, and why is J=0J=0 lowest in energy?

The physical picture

Classically, adding two angular momentum vectors is trivial: you add the vectors component by component, and the magnitude of the sum can be anything between ∣j1−j2∣|j_1 - j_2| and j1+j2j_1+j_2 depending on the relative orientation. Quantum mechanically the situation is both similar and strikingly different. The allowed magnitudes of the total do fall in that same classical range — but instead of a continuum of possible resultant lengths, only a discrete ladder of values is allowed, each differing from its neighbor by exactly one unit of ℏ\hbar. The reason is representation theory in disguise. Each individual angular momentum J1\mathbf{J}_1 and J2\mathbf{J}_2 furnishes an irreducible representation of the rotation group, of dimension 2j1+12j_1+1 and 2j2+12j_2+1 respectively. The combined system's Hilbert space, of dimension (2j1+1)(2j2+1)(2j_1+1)(2j_2+1), carries a reducible representation of the rotation group (rotations act on both particles simultaneously). Decomposing that reducible representation into irreducibles — each irreducible piece being a definite total angular momentum JJ — is exactly the Clebsch-Gordan decomposition. This is the same logic that appears whenever Ladder Operators and the Angular Momentum Algebra are used to build the ∣j,m⟩|j,m\rangle ladder for a single angular momentum; coupling simply asks how two such ladders combine. Physically, the uncoupled states ∣j1m1⟩∣j2m2⟩|j_1 m_1\rangle|j_2 m_2\rangle are eigenstates of J1zJ_{1z} and J2zJ_{2z} separately — they describe a situation where each angular momentum vector has a definite projection on the zz-axis, but the relative orientation of J1\mathbf{J}_1 and J2\mathbf{J}_2 is completely undetermined (equivalently, J1⋅J2\mathbf{J}_1\cdot\mathbf{J}_2 has no definite value). The coupled states ∣JM⟩|JM\rangle instead have a definite total angular momentum magnitude J(J+1)ℏ2J(J+1)\hbar^2 and a definite zz-projection MℏM\hbar of the sum, but the individual zz-projections m1,m2m_1, m_2 are now indefinite — the two vectors are precessing around their fixed resultant. Neither description is more "real" than the other; they are two orthonormal bases for the identical physical Hilbert space, and the correct choice is whichever one diagonalizes the actual Hamiltonian governing the system, exactly as prefigured by the Spin Angular Momentum and the Stern-Gerlach Experiment discussion of measurement bases.

Figure
Vector model for coupling j1=1j_1=1 and j2=12j_2=\tfrac12: the five allowed MM states group into J=3/2J=3/2 (4 states) and J=1/2J=1/2 (2 states)
Uncoupled basis |m₁,m₂⟩ — 3×2 = 6 states m₁=+1, m₂=+½M = 3/2 m₁=+1, m₂=−½M = 1/2 m₁=0, m₂=+½M = 1/2 m₁=0, m₂=−½M = −1/2 m₁=−1, m₂=+½M = −1/2 m₁=−1, m₂=−½M = −3/2 Coupled basis |J,M⟩ J = 3/2: M = 3/2, 1/2, −1/2, −3/2 J = 1/2: M = 1/2, −1/2 4 states (J=3/2) + 2 states (J=1/2) = 6 ✓ Triangle rule: |1 − ½| = ½ ≤ J ≤ 1 + ½ = 3/2 so J = 3/2 or J = 1/2 only — each M sector's degeneracy (2,2,1,1 for M=3/2,½,−½,−3/2) fixes the count.

The dimension count in that figure is not a coincidence — it is the master consistency check for any coupling problem, however large j1j_1 and j2j_2 are. Every valid Clebsch-Gordan decomposition must satisfy

(2j1+1)(2j2+1)=∑J=∣j1−j2∣j1+j2(2J+1)(2j_1+1)(2j_2+1) = \sum_{J=|j_1-j_2|}^{j_1+j_2} (2J+1)

where:

  • 2j1+12j_1+1, 2j2+12j_2+1 are the dimensions of the individual irreducible spaces,
  • the sum runs over the triangle-rule ladder of allowed total JJ values in integer steps,
  • 2J+12J+1 is the dimension (number of MM states) of each total-JJ multiplet.

Reading: the left side counts uncoupled basis states ∣m1⟩∣m2⟩|m_1\rangle|m_2\rangle; the right side counts coupled basis states ∣JM⟩|JM\rangle grouped into multiplets. Since both bases span the same space, the totals must match — and they always do, because the sum of consecutive odd (or the appropriate arithmetic) integers telescopes exactly. A second essential picture is how the maximum MM value pins down where the ladder starts. Since m1m_1 ranges over 2j1+12j_1+1 values and m2m_2 over 2j2+12j_2+1 values, the largest possible sum is Mmax=j1+j2M_{max} = j_1+j_2, achieved in only one way: m1=j1,m2=j2m_1 = j_1, m_2 = j_2. That unique state must belong to a multiplet with J=j1+j2J = j_1+j_2 (since no larger JJ can exist — nothing produces a bigger MM). This is the top rung of the tallest ladder. Descending one unit at a time and comparing degeneracies at each MM level to how many are already "used up" by higher JJ multiplets, reveals successively that J=j1+j2−1,j1+j2−2,…J = j_1+j_2-1, j_1+j_2-2,\ldots must each appear exactly once, terminating at J=∣j1−j2∣J=|j_1-j_2| — precisely the triangle rule.

Figure
Degeneracy counting for j₁=2, j₂=1: reading off J values from the MM-degeneracy staircase
j₁ = 2, j₂ = 1 → 5×3 = 15 uncoupled states Mdegeneracy g(M)new J introduced 31J=3 starts 22J=2 starts 13J=1 starts 03(none new — all used) −13 −22 −31 J=3: 7 states (M=3…−3) J=2: 5 states (M=2…−2) J=1: 3 states (M=1…−1) 7 + 5 + 3 = 15 ✓ (matches 5×3 above) Triangle rule: |2−1|=1 ≤ J ≤ 2+1=3, so J=1,2,3 — no J=0, consistent with M=0 needing no *new* multiplet beyond J=1.

The core ideas

The starting point is the definition of the total angular momentum operator as a vector sum acting on the tensor product Hilbert space:

J=J1⊗12+11⊗J2≡J1+J2\mathbf{J} = \mathbf{J}_1 \otimes \mathbb{1}_2 + \mathbb{1}_1 \otimes \mathbf{J}_2 \equiv \mathbf{J}_1 + \mathbf{J}_2

where:

  • J1\mathbf{J}_1 acts only on the first particle's/subsystem's degrees of freedom, J2\mathbf{J}_2 only on the second's,
  • the tensor product structure guarantees [J1,J2]=0[\mathbf{J}_1, \mathbf{J}_2] = 0 (they act on different factors),
  • J\mathbf{J} is then shown to satisfy the standard angular momentum commutation relations [Jx,Jy]=iℏJz[J_x,J_y]=i\hbar J_z (cyclic) purely as a consequence of J1\mathbf{J}_1 and J2\mathbf{J}_2 separately satisfying them.

Reading: because J1J_1 and J2J_2 commute with each other (different subsystems), their sum automatically obeys the same algebra as either one alone — this is why J2J^2 and JzJ_z can be simultaneously diagonalized just like for a single angular momentum, giving a legitimate ∣JM⟩|JM\rangle ladder. The triangle rule (also called the Clebsch-Gordan series) states the allowed total quantum numbers:

J=j1+j2,  j1+j2−1,  …,  ∣j1−j2∣J = j_1+j_2,\; j_1+j_2-1,\; \ldots,\; |j_1-j_2|

where:

  • j1,j2j_1, j_2 are the (fixed) individual angular momentum quantum numbers, integer or half-integer,
  • JJ steps down by exactly 1 from the maximum to the minimum — never skipping,
  • each allowed JJ value occurs exactly once in the decomposition (multiplicity 1), not more.

Reading: this is the master selection rule for what totals are even possible — before asking about probabilities or coefficients, this tells you the complete list of total-JJ multiplets present. The four commuting operators of the coupled scheme, and their simultaneous eigenstates, are:

J12∣j1j2;JM⟩=j1(j1+1)ℏ2∣j1j2;JM⟩,J22∣j1j2;JM⟩=j2(j2+1)ℏ2∣j1j2;JM⟩J_1^2|j_1 j_2;JM\rangle = j_1(j_1+1)\hbar^2|j_1j_2;JM\rangle,\quad J_2^2|j_1j_2;JM\rangle = j_2(j_2+1)\hbar^2|j_1j_2;JM\rangle
J2∣j1j2;JM⟩=J(J+1)ℏ2∣j1j2;JM⟩,Jz∣j1j2;JM⟩=Mℏ ∣j1j2;JM⟩J^2|j_1j_2;JM\rangle = J(J+1)\hbar^2|j_1j_2;JM\rangle,\quad J_z|j_1j_2;JM\rangle = M\hbar\,|j_1j_2;JM\rangle

where:

  • the first pair fixes the individual magnitudes j1,j2j_1, j_2 (these don't change under coupling — J12J_1^2 commutes with everything since [J12,J2]=0[\mathbf J_1^2,\mathbf J_2]=0),
  • J2=J12+J22+2J1⋅J2J^2 = J_1^2+J_2^2+2\mathbf J_1\cdot\mathbf J_2 is the total-squared-magnitude operator,
  • M=m1+m2M = m_1+m_2 always, since Jz=J1z+J2zJ_z = J_{1z}+J_{2z} commutes with everything in both bases.

Reading: {J12,J22,J2,Jz}\{J_1^2,J_2^2,J^2,J_z\} is a complete set of commuting observables — an equally valid alternative to {J12,J22,J1z,J2z}\{J_1^2,J_2^2,J_{1z},J_{2z}\} — and the coupled kets ∣j1j2;JM⟩|j_1j_2;JM\rangle are simultaneous eigenkets of the first set. The connection between the two complete bases is the Clebsch-Gordan expansion:

∣j1j2;JM⟩=∑m1,m2⟨j1m1 j2m2∣JM⟩ ∣j1m1⟩∣j2m2⟩|j_1 j_2; JM\rangle = \sum_{m_1,m_2} \langle j_1 m_1\, j_2 m_2 | JM\rangle\, |j_1 m_1\rangle|j_2 m_2\rangle

where:

  • ⟨j1m1 j2m2∣JM⟩\langle j_1 m_1\, j_2 m_2|JM\rangle is the Clebsch-Gordan coefficient — a real number (by convention) fixed by the algebra, tabulated in standard references,
  • the sum is constrained by m1+m2=Mm_1+m_2 = M (all other terms vanish identically),
  • the coefficients form a real orthogonal matrix, so the inverse transformation (uncoupled in terms of coupled) uses the same coefficients.

Reading: this equation is the precise statement of "coupled states are superpositions of uncoupled states with definite m1+m2=Mm_1+m_2=M" — the coefficients are what carries all the quantitative information about how the two angular momenta correlate in a state of definite total J,MJ,M.

Derivation ADeriving the triangle rule by counting states

Goal. Show from first principles, using only dimension-counting of MM values, that coupling j1j_1 and j2j_2 produces exactly one copy each of J=j1+j2,j1+j2−1,…,∣j1−j2∣J=j_1+j_2, j_1+j_2-1,\ldots,|j_1-j_2| and no others.

Step 1 — Enumerate the uncoupled states by MM. The uncoupled basis has (2j1+1)(2j2+1)(2j_1+1)(2j_2+1) states ∣m1⟩∣m2⟩|m_1\rangle|m_2\rangle, each with m1∈{−j1,…,j1}m_1\in\{-j_1,\ldots,j_1\} and m2∈{−j2,…,j2}m_2\in\{-j_2,\ldots,j_2\} in integer steps. Group them by total M=m1+m2M=m_1+m_2. Assume without loss of generality j1≥j2j_1\ge j_2. The degeneracy g(M)g(M) — the number of (m1,m2)(m_1,m_2) pairs giving a fixed MM — rises linearly from g(j1+j2)=1g(j_1+j_2)=1 up to a plateau of g(M)=2j2+1g(M) = 2j_2+1 for ∣M∣≤j1−j2|M|\le j_1-j_2, then falls symmetrically back to g(−j1−j2)=1g(-j_1-j_2)=1. This is just the discrete convolution of two "boxcar" ranges.

Step 2 — Identify the top multiplet. The unique state at M=j1+j2M=j_1+j_2 (only m1=j1,m2=j2m_1=j_1,m_2=j_2 gives this) cannot belong to any J>j1+j2J>j_1+j_2, since such JJ would require an even larger MM state, which doesn't exist in the space. So this state initiates a J=j1+j2J=j_1+j_2 multiplet — a full ladder of 2(j1+j2)+12(j_1+j_2)+1 states running from M=j1+j2M=j_1+j_2 down to M=−(j1+j2)M=-(j_1+j_2), built by repeatedly applying J−J_- (the Ladder Operators and the Angular Momentum Algebra machinery guarantees this ladder exists and doesn't truncate early since J2J^2 eigenvalue is fixed).

Step 3 — Subtract and repeat. At M=j1+j2−1M=j_1+j_2-1 the raw degeneracy is g(j1+j2−1)=2g(j_1+j_2-1)=2 (from (m1,m2)=(j1,j2−1)(m_1,m_2)=(j_1,j_2-1) and (j1−1,j2)(j_1-1,j_2)). One of these two states is already accounted for as the second rung of the J=j1+j2J=j_1+j_2 ladder from Step 2. The remaining one state must be the top of a new ladder, necessarily J=j1+j2−1J=j_1+j_2-1 (it cannot be any higher JJ, since a higher JJ ladder would already have shown up at a higher MM, and it cannot be lower without violating M≤JM\le J). This generates a full J=j1+j2−1J=j_1+j_2-1 multiplet of 2(j1+j2−1)+12(j_1+j_2-1)+1 states.

Step 4 — Iterate down to the plateau. Continue this bookkeeping: at each descending integer step in MM, the degeneracy g(M)g(M) increases by exactly 1 (while below the plateau), meaning exactly one new JJ-multiplet is introduced at each step, until M=j1−j2M = j_1-j_2, where the degeneracy plateaus at 2j2+12j_2+1. At that point g(M)g(M) stops increasing — meaning the last new multiplet introduced is J=j1−j2=∣j1−j2∣J=j_1-j_2 = |j_1-j_2| (using j1≥j2j_1\ge j_2), and no further multiplets are introduced below the plateau, since all subsequent g(M)g(M) growth (there is none — it's flat then decreasing) is already saturated by the 2j2+12j_2+1 multiplets J=j1−j2,…,j1+j2J = j_1-j_2,\ldots,j_1+j_2 already accounted for.

Step 5 — Check the total. Summing dimensions of these multiplets:

∑J=∣j1−j2∣j1+j2(2J+1)=(2j2+1)(2j1+1)\sum_{J=|j_1-j_2|}^{j_1+j_2}(2J+1) = (2j_2+1)(2j_1+1)

which is an arithmetic series with 2j2+12j_2+1 terms, first term 2∣j1−j2∣+12|j_1-j_2|+1, last term 2(j1+j2)+12(j_1+j_2)+1; the sum of an arithmetic series equals (number of terms)×(average of first and last), giving exactly (2j2+1)×(2j1+1)(2j_2+1)\times(2j_1+1). This matches the uncoupled count exactly, confirming no multiplets are missing or double-counted. ■\blacksquare

Go deeper · Derivation B · 1 minBuilding coupled states with ladder operators — the j1=1,j2=12j_1=1,j_2=\tfrac12 example explicitly
Derivation BBuilding coupled states with ladder operators — the j1=1,j2=12j_1=1,j_2=\tfrac12 example explicitly

Goal. Construct all six coupled states ∣JM⟩|JM\rangle for j1=1, j2=12j_1=1,\ j_2=\tfrac12 explicitly as linear combinations of uncoupled kets ∣m1m2⟩|m_1 m_2\rangle, using only raising/lowering operators — the concrete machine behind the abstract triangle rule.

Step 1 — Identify the top state. The maximum MM is j1+j2=3/2j_1+j_2 = 3/2, achieved uniquely by m1=1,m2=+12m_1=1,m_2=+\tfrac12. By the argument of the general derivation, this state has definite J=3/2J=3/2:

∣J=32,M=32⟩=∣1,+1⟩ ∣12,+12⟩|J{=}\tfrac32,M{=}\tfrac32\rangle = |1,+1\rangle\,\big|\tfrac12,+\tfrac12\big\rangle

(using shorthand ∣m1⟩∣m2⟩|m_1\rangle|m_2\rangle).

Step 2 — Lower once to get J=3/2,M=1/2J=3/2,M=1/2. Apply the total lowering operator J−=J1−+J2−J_- = J_{1-}+J_{2-} to both sides. On the left, the single-particle ladder formula gives

J−∣J,M⟩=ℏJ(J+1)−M(M−1) ∣J,M−1⟩J_-|J,M\rangle = \hbar\sqrt{J(J+1)-M(M-1)}\,|J,M-1\rangle

with J=3/2,M=3/2J=3/2,M=3/2: the coefficient is ℏ32⋅52−32⋅12=ℏ154−34=ℏ3\hbar\sqrt{\tfrac32\cdot\tfrac52-\tfrac32\cdot\tfrac12}=\hbar\sqrt{\tfrac{15}{4}-\tfrac34}=\hbar\sqrt{3}. On the right, J1−J_{1-} acting on ∣1,+1⟩|1,+1\rangle gives ℏ1⋅2−1⋅0 ∣1,0⟩=ℏ2 ∣1,0⟩\hbar\sqrt{1\cdot2-1\cdot0}\,|1,0\rangle=\hbar\sqrt2\,|1,0\rangle (times ∣12,+12⟩|\tfrac12,+\tfrac12\rangle unchanged), and J2−J_{2-} acting on ∣12,+12⟩|\tfrac12,+\tfrac12\rangle gives ℏ12⋅32−12⋅(−12) ∣12,−12⟩=ℏ ∣12,−12⟩\hbar\sqrt{\tfrac12\cdot\tfrac32-\tfrac12\cdot(-\tfrac12)}\,|\tfrac12,-\tfrac12\rangle=\hbar\,|\tfrac12,-\tfrac12\rangle (times ∣1,+1⟩|1,+1\rangle unchanged). So:

ℏ3 ∣J=32,M=12⟩=ℏ2 ∣1,0⟩∣12,12⟩+ℏ ∣1,1⟩∣12,−12⟩\hbar\sqrt3\,|J{=}\tfrac32,M{=}\tfrac12\rangle = \hbar\sqrt2\,|1,0\rangle|\tfrac12,\tfrac12\rangle + \hbar\,|1,1\rangle|\tfrac12,-\tfrac12\rangle

Dividing through by ℏ3\hbar\sqrt3:

∣J=32,M=12⟩=23 ∣1,0⟩∣12,12⟩+13 ∣1,1⟩∣12,−12⟩|J{=}\tfrac32,M{=}\tfrac12\rangle = \sqrt{\tfrac23}\,|1,0\rangle\big|\tfrac12,\tfrac12\big\rangle + \sqrt{\tfrac13}\,|1,1\rangle\big|\tfrac12,-\tfrac12\big\rangle

Step 3 — Find J=1/2,M=1/2J=1/2,M=1/2 by orthogonality. At M=1/2M=1/2 there are exactly two uncoupled states, ∣1,0⟩∣12,12⟩|1,0\rangle|\tfrac12,\tfrac12\rangle and ∣1,1⟩∣12,−12⟩|1,1\rangle|\tfrac12,-\tfrac12\rangle, and exactly two coupled states must live in this two-dimensional subspace: the one just found (J=3/2J=3/2) and one orthogonal partner, which must be J=1/2,M=1/2J=1/2,M=1/2 (the only other allowed JJ by the triangle rule ∣1−12∣≤J≤1+12|1-\tfrac12|\le J\le 1+\tfrac12). Orthogonality (with a sign convention fixing the coefficient of the stretched configuration m1=j1m_1=j_1 to be positive) gives:

∣J=12,M=12⟩=13 ∣1,1⟩∣12,−12⟩−23 ∣1,0⟩∣12,12⟩|J{=}\tfrac12,M{=}\tfrac12\rangle = \sqrt{\tfrac13}\,|1,1\rangle\big|\tfrac12,-\tfrac12\big\rangle - \sqrt{\tfrac23}\,|1,0\rangle\big|\tfrac12,\tfrac12\big\rangle

Step 4 — Complete the J=3/2J=3/2 ladder by lowering again. Applying J−J_- once more to ∣J=32,M=12⟩|J{=}\tfrac32,M{=}\tfrac12\rangle (coefficient ℏ32⋅52−12(−12)=ℏ4=2ℏ\hbar\sqrt{\tfrac32\cdot\tfrac52-\tfrac12(-\tfrac12)}=\hbar\sqrt{4}=2\hbar on the left) and carrying through the same bookkeeping on the right yields, after normalizing:

∣J=32,M=−12⟩=13 ∣1,−1⟩∣12,12⟩+23 ∣1,0⟩∣12,−12⟩|J{=}\tfrac32,M{=}-\tfrac12\rangle = \sqrt{\tfrac13}\,|1,-1\rangle\big|\tfrac12,\tfrac12\big\rangle + \sqrt{\tfrac23}\,|1,0\rangle\big|\tfrac12,-\tfrac12\big\rangle
∣J=32,M=−32⟩=∣1,−1⟩∣12,−12⟩|J{=}\tfrac32,M{=}-\tfrac32\rangle = |1,-1\rangle\big|\tfrac12,-\tfrac12\big\rangle

and orthogonality at M=−1/2M=-1/2 gives the last state, ∣J=12,M=−12⟩=23∣1,−1⟩∣12,12⟩−13∣1,0⟩∣12,−12⟩|J{=}\tfrac12,M{=}-\tfrac12\rangle = \sqrt{\tfrac23}|1,-1\rangle|\tfrac12,\tfrac12\rangle - \sqrt{\tfrac13}|1,0\rangle|\tfrac12,-\tfrac12\rangle.

Step 5 — Verify closure. All six coupled states (44 from J=3/2J{=}3/2, 22 from J=1/2J{=}1/2) are now expressed purely in the six uncoupled kets, the transformation matrix is orthogonal (each row and column normalized, cross terms vanishing — checkable directly from the coefficients 2/3,1/3\sqrt{2/3},\sqrt{1/3} etc.), and dimension is conserved exactly as the figure above illustrated. ■\blacksquare

Worked examples

Example 1 — Foundational: coupling two spin-1/2 electrons. Two electrons each have j1=j2=12j_1=j_2=\tfrac12. The triangle rule gives

J=12+12, 12+12−1=1, 0J = \tfrac12+\tfrac12,\, \tfrac12+\tfrac12-1 = 1,\,0

so total spin S=1S=1 (triplet, 3 states) or S=0S=0 (singlet, 1 state), totaling 3+1=4=2×23+1=4=2\times2 states — matching (2⋅12+1)2=4(2\cdot\tfrac12+1)^2=4. Explicitly, using the ladder-operator method above (or simply pattern-matching the j1=j2=1/2j_1=j_2=1/2 Clebsch-Gordan table):

∣S=1,MS=1⟩=∣ ⁣↑↑⟩,∣S=1,MS=0⟩=12(∣ ⁣↑↓⟩+∣ ⁣↓↑⟩),∣S=1,MS=−1⟩=∣ ⁣↓↓⟩|S{=}1,M_S{=}1\rangle = |\!\uparrow\uparrow\rangle,\quad |S{=}1,M_S{=}0\rangle = \tfrac{1}{\sqrt2}\big(|\!\uparrow\downarrow\rangle+|\!\downarrow\uparrow\rangle\big),\quad |S{=}1,M_S{=}-1\rangle=|\!\downarrow\downarrow\rangle
∣S=0,MS=0⟩=12(∣ ⁣↑↓⟩−∣ ⁣↓↑⟩)|S{=}0,M_S{=}0\rangle = \tfrac{1}{\sqrt2}\big(|\!\uparrow\downarrow\rangle - |\!\downarrow\uparrow\rangle\big)

As a limiting-case check: setting one spin to "off" conceptually (i.e., taking j2→0j_2\to0 formally) would collapse the coupling to J=j1J=j_1 only, a single multiplet with the uncoupled and coupled bases coinciding — consistent with the general rule that coupling to j2=0j_2=0 never splits anything, since the triangle rule with j2=0j_2=0 forces J=j1J=j_1 exclusively. This singlet/triplet decomposition is precisely what underlies the Pauli-allowed spatial symmetries in the helium atom and in covalent bond formation between two electrons.

Example 2 — Applied: coupling LL and SS for a 3P^3P carbon-like term. For the ground configuration of carbon, 2p22p^2, Hund's rules and the Pauli principle restrict the allowed terms to 1S, 1D, 3P^1S, \,^1D,\,^3P (a standard exercise in term-symbol construction). Focus on the 3P^3P term: here L=1L=1 (P) and S=1S=1 (triplet, spin multiplicity 2S+1=32S+1=3). Coupling LL and SS via the same triangle rule (now applied to orbital and spin angular momentum rather than two orbital or two spin momenta) gives total JJ:

J=∣L−S∣,…,L+S=∣1−1∣,…,1+1=0, 1, 2J = |L-S|,\ldots, L+S = |1-1|,\ldots,1+1 = 0,\,1,\,2

So the 3P^3P term splits into three fine-structure levels: 3P0,3P1,3P2^3P_0, ^3P_1, ^3P_2, with degeneracies 2J+1=1,3,52J+1 = 1,3,5. As a numerical check on total degeneracy: the full 3P^3P term (before spin-orbit splitting) has degeneracy (2L+1)(2S+1)=3×3=9(2L+1)(2S+1)=3\times3=9, and indeed 1+3+5=91+3+5=9 — the L,SL,S multiplet redistributes exactly, level by level, into the coupled JJ basis, with no states lost or gained. (In reality, carbon's 3P0^3P_0 ground level lies about 16.4 cm−116.4\ \text{cm}^{-1} below 3P1^3P_1 and 43.4 cm−143.4\ \text{cm}^{-1} below 3P2^3P_2, an ordering and spacing set by the sign and magnitude of the spin-orbit constant — the subject that this coupling machinery feeds directly into.)

Example 3 — Exam/competition level: three coupled angular momenta and a nontrivial branching count. Consider coupling three angular momenta with j1=32, j2=1, j3=12j_1 = \tfrac32,\ j_2=1,\ j_3=\tfrac12 (a stylized model, e.g., relevant to hyperfine coupling of two electron-like momenta plus a nuclear spin). Coupling is done pairwise: first combine j1j_1 and j2j_2 into an intermediate j12j_{12}, then couple j12j_{12} with j3j_3 to get the final JJ. Step A. j1=32,j2=1j_1=\tfrac32, j_2=1: triangle rule gives j12=32+1,…,∣32−1∣=52, 32, 12j_{12} = \tfrac32+1,\ldots,|\tfrac32-1| = \tfrac52,\,\tfrac32,\,\tfrac12. Step B. For each j12j_{12}, couple with j3=12j_3=\tfrac12:

  • j12=52j_{12}=\tfrac52: J=3, 2J = 3,\,2
  • j12=32j_{12}=\tfrac32: J=2, 1J = 2,\,1
  • j12=12j_{12}=\tfrac12: J=1, 0J = 1,\,0

So the final list of total JJ values (with multiplicities, since J=2J=2 and J=1J=1 each arise from two different j12j_{12} "routes") is:

J=3 (×1),J=2 (×2),J=1 (×2),J=0 (×1)J=3\ (\times1),\quad J=2\ (\times2),\quad J=1\ (\times2),\quad J=0\ (\times1)

Dimension check. Total uncoupled dimension: (2⋅32+1)(2⋅1+1)(2⋅12+1)=4×3×2=24(2\cdot\tfrac32+1)(2\cdot1+1)(2\cdot\tfrac12+1) = 4\times3\times2=24. Coupled dimension:

1×(2⋅3+1)+2×(2⋅2+1)+2×(2⋅1+1)+1×(2⋅0+1)1\times(2\cdot3+1) + 2\times(2\cdot2+1) + 2\times(2\cdot1+1) + 1\times(2\cdot0+1)
=1×7+2×5+2×3+1×1=7+10+6+1=24= 1\times7 + 2\times5 + 2\times3 + 1\times1 = 7+10+6+1 = 24

which matches exactly, confirming the branching diagram is complete and correctly enumerated — and illustrating the important general point that when JJ is reachable by more than one intermediate-coupling route, its overall multiplicity in the final decomposition is genuinely greater than one, unlike the two-body case where every JJ appears exactly once.

Limiting cases, scope & common traps

Regime / limitWhat happensWhy
j2=0j_2 = 0Only J=j1J=j_1 survives, coupled = uncoupled basisTriangle rule collapses to a single value; coupling a "nothing" changes nothing
j1=j2=jj_1 = j_2 = jJJ ranges over 0,1,…,2j0,1,\ldots,2j, all appearing onceFull triangle rule range from ∣j−j∣=0|j-j|=0 to j+j=2jj+j=2j
j1≫j2j_1 \gg j_2 (classical limit)J≈j1J \approx j_1, small angular deviationsThe big vector is barely perturbed by the small one — recovers vector-addition intuition
Identical fermions, j1=j2=12j_1=j_2=\tfrac12Only symmetric S=1S=1 combines with antisymmetric spatial part, and vice versaPauli antisymmetry couples spin-coupling parity to spatial exchange symmetry
Three or more angular momentaOrder of pairwise coupling matters for intermediate labels (not final JJ set)Recoupling is described by 6-j symbols; final JJ spectrum is basis-independent but intermediate multiplet labels are not
Weak-field vs strong-field coupling (LS vs jj)LS coupling valid when residual electrostatic repulsion ≫ spin-orbit; jj coupling valid in the reverse (heavy atoms)Different terms in HH dominate, selecting a different "good" basis

The chemist behind it

The mathematics of angular momentum coupling was worked out by physicists rather than chemists, but its constant use in atomic spectroscopy makes it as much a tool of the physical chemist as of the physicist. The two names permanently attached to the coefficients are Alfred Clebsch and Paul Gordan, nineteenth-century mathematicians who developed the underlying theory of invariants of binary forms — group representation theory did not yet exist as such, and it was Eugene Wigner decades later who recognized that their algebraic formulas were exactly the coefficients needed to decompose products of rotation-group representations, cementing the "Clebsch-Gordan coefficient" as the standard tool of atomic and nuclear spectroscopy.

Origins · 1 minEugene Paul Wigner and how this was discovered
EP
Eugene Paul Wigner
1902–1995 · Princeton, mid-20th century group theory & quantum mechanics

Trained as a chemical engineer in Berlin before turning fully to theoretical physics, Wigner brought a chemist's comfort with symmetry and structure to the formal machinery of quantum mechanics. He recognized that the abstract classification of atomic and nuclear states — the term symbols, selection rules, and multiplet structures that spectroscopists had cataloged empirically — was in fact controlled entirely by the representation theory of the rotation group. His 1931 book Gruppentheorie und ihre Anwendung auf die Quantenmechanik der Atomspektren systematized angular momentum coupling, introduced the 3-j and 6-j symbols that generalize the Clebsch-Gordan coefficients to more complex coupling schemes, and laid the foundation for the Wigner-Eckart theorem governing matrix elements of tensor operators between coupled angular momentum states. He shared the 1963 Nobel Prize in Physics for these contributions to the theory of the atomic nucleus and elementary particles, in particular through the discovery and application of fundamental symmetry principles.

History
- 1868–1872 — Clebsch and Gordan develop invariant-theoretic identities for binary forms that later supply the coupling coefficients.
- 1925–1927 — Uhlenbeck and Goudsmit propose electron spin; Pauli formalizes spin operators, setting up two-angular-momentum problems in atomic structure.
- 1927 — Russell and Saunders classify atomic spectral terms via LS coupling, an empirical scheme angular momentum addition later justifies rigorously.
- 1931 — Wigner's Gruppentheorie book formalizes angular momentum coupling via representation theory, introducing 3-j and 6-j symbols.
- 1963 — Wigner shares the Nobel Prize in Physics for symmetry-principle contributions built substantially on this coupling formalism.
Why can't the total JJ just be any value between ∣j1−j2∣|j_1-j_2| and j1+j2j_1+j_2, continuously?

Quantum mechanically, J2J^2 is quantized: its eigenvalues must be J(J+1)ℏ2J(J+1)\hbar^2 for JJ a non-negative half-integer or integer, exactly as for any single angular momentum. This quantization is forced by the requirement that the ladder operators J±J_\pm terminate after a finite number of steps (otherwise you'd generate negative-norm states), the same argument used to derive the single-particle ∣j,m⟩|j,m\rangle spectrum from Ladder Operators and the Angular Momentum Algebra. Since J=J1+J2\mathbf J=\mathbf J_1+\mathbf J_2 obeys the identical commutation algebra as a single angular momentum, its eigenvalues must obey the identical discreteness — hence only integer-spaced JJ values within the classical bound ∣j1−j2∣≤J≤j1+j2|j_1-j_2|\le J\le j_1+j_2 are allowed, not a continuum.

Is the coupled basis always the "correct" one to use physically?

No — correctness is dictated entirely by the Hamiltonian. If the true HH commutes with J12,J22,J2,JzJ_1^2, J_2^2, J^2, J_z (for instance, isotropic exchange or spin-orbit coupling that depends on J1⋅J2\mathbf J_1\cdot\mathbf J_2), then the coupled basis diagonalizes HH and gives the true energy eigenstates. If instead HH commutes with J1zJ_{1z} and J2zJ_{2z} individually (e.g., a strong external magnetic field that couples separately to each moment, the Paschen-Back regime), the uncoupled basis is the physically appropriate one. Many real systems interpolate between these limits, requiring diagonalization of HH in either basis as an intermediate-coupling problem.

Do Clebsch-Gordan coefficients have a simple closed formula?

There is a general closed-form expression (due to Wigner) involving factorials and a sum over an auxiliary integer, but it is rarely used by hand for anything beyond the smallest cases. In practice, chemists and physicists either look up tabulated values for small j1,j2j_1,j_2 (widely available in standard references) or generate the coupled states directly by repeated application of ladder operators combined with orthogonality, exactly as shown in the worked derivation above — this is faster and more transparent than the general formula for most chemically relevant cases (spins of 12\tfrac12, 11; orbital LL up to a few units).

How does this relate to Russell-Saunders (LS) versus jj coupling in atomic spectroscopy?

LS coupling first combines all individual electron orbital angular momenta ℓi\ell_i into a total LL and all spins sis_i into a total SS (using the triangle rule repeatedly for multi-electron atoms), then combines LL and SS into JJ. This order is appropriate when electrostatic repulsion between electrons (which depends on LL and SS but not JJ directly) dominates over spin-orbit coupling — true for light and moderately heavy atoms. In heavy atoms, spin-orbit coupling on each electron individually can dominate, favoring instead coupling each electron's own ℓi\ell_i and sis_i into jij_i first, then combining the jij_i's into total JJ — this is jjjj coupling. Both schemes use identical Clebsch-Gordan mathematics; they differ only in which pairwise coupling order best matches the dominant physical interaction, a distinction elaborated fully in Spin-Orbit Coupling and Fine Structure.

What role do Clebsch-Gordan coefficients play in selection rules for spectroscopic transitions?

Via the Wigner-Eckart theorem, the matrix element of any spherical tensor operator (such as the electric dipole operator, which transforms as a rank-1 tensor) between two angular-momentum eigenstates factors into a Clebsch-Gordan coefficient (governing which J,MJ,M combinations are geometrically allowed) times a reduced matrix element independent of orientation quantum numbers. The Clebsch-Gordan coefficient vanishes identically unless the triangle rule is satisfied among the operator's rank and the initial/final JJ values, which is precisely the origin of angular momentum selection rules like ΔJ=0,±1\Delta J=0,\pm1 (with J=0→J=0J=0\to J=0 forbidden) for electric dipole transitions.

Can angular momentum coupling apply to orbital angular momentum of light, not just particles?

Yes. Photons carry both spin angular momentum (polarization, s=1s=1) and, for structured light beams (Laguerre-Gaussian modes), orbital angular momentum with integer values ℓ\ell. Coupling these — or coupling a photon's total angular momentum to an atom's electronic angular momentum during absorption/emission — follows the identical triangle-rule and Clebsch-Gordan formalism developed here, since photon angular momentum obeys the same rotation-group algebra as electron or nuclear angular momentum. This underlies selection rules for multipole radiation and modern optical-vortex and orbital-angular-momentum-multiplexing technologies.

Why does the total spin of two electrons matter so much for chemical bonding?

The Pauli exclusion principle requires the total electronic wavefunction (spatial × spin) to be antisymmetric under exchange of the two electrons. Coupling two spin-12\tfrac12 electrons gives a symmetric triplet (S=1S=1) or antisymmetric singlet (S=0S=0), as derived above. Antisymmetry of the total wavefunction then forces the spatial part to be antisymmetric when spins are triplet-coupled, and symmetric when singlet-coupled. A symmetric spatial wavefunction typically allows greater electron density between nuclei — this is exactly why the ground state of H₂ is a spin singlet with a symmetric (bonding) spatial orbital, while the corresponding triplet state is dissociative (antibonding character), tying the abstract angular momentum coupling directly to whether a covalent bond forms at all.

Before this

Where this leads