Chapter 12: Angular Momentum Algebra: Raising, Lowering, and the General Theory

"The algebraic approach to angular momentum is one of the most beautiful things in all of physics. You start with three simple commutation relations and end up with the entire spectrum of angular momentum — including half-integer values that you never would have guessed from orbital motion." — J. J. Sakurai, paraphrased

"In quantum mechanics, algebra is not the servant of physics — it is physics." — P. A. M. Dirac

This chapter is the crown jewel of angular momentum theory. We are about to do something extraordinary: starting from nothing more than three commutation relations — abstract algebraic rules that define what we mean by "angular momentum" in quantum mechanics — we will derive the complete eigenvalue spectrum, construct every matrix representation, and build the rotation operators that describe how quantum states transform under physical rotations.

The result is one of the most powerful demonstrations of a theme that runs throughout this textbook: mathematics and physics are inseparable in quantum mechanics. The formalism is not merely a tool for calculating physical predictions. The formalism is the physics. The algebra itself tells us that angular momentum can take half-integer values — a prediction that has no classical analogue and that experiment spectacularly confirms.

If you have worked through Chapter 8 (Dirac notation and matrix representations), you have everything you need. We will use Dirac notation exclusively. If you recall the orbital angular momentum from Chapter 5, that will provide valuable physical intuition, but our derivation will not depend on it. Indeed, the entire point is that we can derive more than Chapter 5 ever could, precisely because we abandon the specific representation in favor of pure algebra.

🏃 Fast Track: If you are confident with the commutation relations and want to get to the heart of the derivation, skip to Section 12.4 (the raising and lowering operators) and Section 12.5 (the eigenvalue spectrum). These two sections contain the essential algebraic argument. Return to Sections 12.1–12.3 as needed for background.

Let us lay out the roadmap for this chapter. We begin in Sections 12.1–12.2 by establishing the commutation relations and the Casimir operator $\hat{J}^2$. Section 12.3 sets up the eigenvalue problem. The central derivation occupies Sections 12.4–12.5: the construction of ladder operators and the derivation of the eigenvalue spectrum. Sections 12.6–12.7 develop the practical tools — matrix representations and rotation matrices. Sections 12.8–12.9 connect the general theory back to the specific cases of orbital angular momentum and spin. Section 12.10 summarizes and introduces the toolkit module.

12.1 Angular Momentum from Commutation Relations Alone

The Classical Prelude

In classical mechanics, angular momentum is defined as $\mathbf{L} = \mathbf{r} \times \mathbf{p}$, a quantity associated with rotational motion. Its components are:

$$L_x = yp_z - zp_y, \quad L_y = zp_x - xp_z, \quad L_z = xp_y - yp_x$$

This definition is tied to the concept of a particle orbiting in space — it requires position and momentum coordinates to even write down. In quantum mechanics, we can do something far more general.

🔄 Retrieval Practice (Ch 5): Before reading further, try to recall from Chapter 5: What are the eigenvalues of $\hat{L}^2$ and $\hat{L}_z$ for orbital angular momentum? What restriction applies to the quantum number $l$? Why does it have to be a non-negative integer?

The answer — $l = 0, 1, 2, \ldots$ with $m = -l, -l+1, \ldots, l$ — came from solving a differential equation (the angular part of the Schrödinger equation in spherical coordinates). That approach is powerful but limited: it can only find angular momenta that correspond to functions of the angles $\theta$ and $\phi$, which forces $l$ to be an integer.

Nature, however, is not so restrictive. Electrons carry spin angular momentum with $s = 1/2$, and no amount of solving differential equations in position space will ever find it. To discover spin, we need a more general approach — one that extracts everything about angular momentum from the algebraic structure alone, without ever mentioning wavefunctions, coordinates, or differential equations.

💡 Key Insight: The algebraic approach to angular momentum is not merely an alternative to the differential-equation approach. It is more general. It finds solutions that the differential-equation approach misses — namely, the half-integer angular momenta. This is why Sakurai and other modern treatments begin here rather than with spherical harmonics.

The Fundamental Postulate

We define a generalized angular momentum $\hat{\mathbf{J}} = (\hat{J}_x, \hat{J}_y, \hat{J}_z)$ as any set of three Hermitian operators satisfying the commutation relations:

$$[\hat{J}_i, \hat{J}_j] = i\hbar \sum_k \epsilon_{ijk} \hat{J}_k$$

Here $\epsilon_{ijk}$ is the Levi-Civita symbol ($+1$ for even permutations of $xyz$, $-1$ for odd, $0$ if any two indices are equal). Written explicitly, the three independent relations are:

$$[\hat{J}_x, \hat{J}_y] = i\hbar \hat{J}_z$$

$$[\hat{J}_y, \hat{J}_z] = i\hbar \hat{J}_x$$

$$[\hat{J}_z, \hat{J}_x] = i\hbar \hat{J}_y$$

These are the angular momentum commutation relations. Any triple of Hermitian operators satisfying them is, by definition, an angular momentum. Orbital angular momentum $\hat{\mathbf{L}}$ satisfies them. Spin angular momentum $\hat{\mathbf{S}}$ satisfies them. Total angular momentum $\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}$ satisfies them. Isospin in nuclear physics satisfies them. The formalism we develop applies to all of these simultaneously.

⚠️ Common Misconception: The commutation relations are not derived from $\hat{\mathbf{L}} = \hat{\mathbf{r}} \times \hat{\mathbf{p}}$ and then "generalized." The logic runs the other way: the commutation relations are the definition of angular momentum in quantum mechanics. That $\hat{\mathbf{L}} = \hat{\mathbf{r}} \times \hat{\mathbf{p}}$ happens to satisfy them is a consequence, not a premise.

A Concrete Verification: Orbital Angular Momentum

Before we explore the consequences of the commutation relations, let us verify that the familiar orbital angular momentum $\hat{\mathbf{L}} = \hat{\mathbf{r}} \times \hat{\mathbf{p}}$ satisfies them. Using $[\hat{x}_i, \hat{p}_j] = i\hbar\delta_{ij}$ (the canonical commutation relation from Chapter 6):

$$[\hat{L}_x, \hat{L}_y] = [\hat{y}\hat{p}_z - \hat{z}\hat{p}_y, \hat{z}\hat{p}_x - \hat{x}\hat{p}_z]$$

Expanding this product requires careful use of the canonical commutation relations. Only terms involving the same position-momentum pair contribute:

$$= \hat{y}\hat{p}_x[\hat{p}_z, \hat{z}] + \hat{x}\hat{p}_y[\hat{z}, \hat{p}_z] = \hat{y}\hat{p}_x(-i\hbar) + \hat{x}\hat{p}_y(i\hbar) = i\hbar(\hat{x}\hat{p}_y - \hat{y}\hat{p}_x) = i\hbar\hat{L}_z$$

The calculation is straightforward but instructive — the canonical commutation relations of position and momentum give rise to the angular momentum commutation relations. This is not an accident: it reflects the fact that orbital angular momentum generates rotations, and rotations of position and momentum are interconnected.

🔴 Warning: The derivation above works only for operators that are functions of $\hat{\mathbf{r}}$ and $\hat{\mathbf{p}}$. Spin angular momentum $\hat{\mathbf{S}}$ is not a function of $\hat{\mathbf{r}}$ and $\hat{\mathbf{p}}$ — it is an independent degree of freedom. Yet it satisfies the same commutation relations. The commutation relations are the deep fact; the specific realization in terms of position and momentum is a special case.

Why These Commutation Relations?

The angular momentum commutation relations are not arbitrary. They are the Lie algebra of the rotation group $SO(3)$ (or more precisely, its double cover $SU(2)$). This means they encode the mathematical structure of rotations in three-dimensional space.

🔄 Retrieval Practice (Ch 10): Recall from Chapter 10 that continuous symmetries are associated with conserved quantities via Noether's theorem, and that the generators of a symmetry transformation satisfy specific algebraic relations. Angular momentum is the generator of rotations, and the commutation relations above are exactly the statement that rotations about different axes do not commute — the order in which you rotate matters.

Consider two successive infinitesimal rotations: first by angle $d\phi_1$ about the $z$-axis, then by $d\phi_2$ about the $x$-axis. If you instead reverse the order, the final state differs by a rotation about the $y$-axis of order $d\phi_1 d\phi_2$. The commutation relation $[\hat{J}_x, \hat{J}_z] = -i\hbar \hat{J}_y$ (equivalently $[\hat{J}_z, \hat{J}_x] = i\hbar \hat{J}_y$) is the infinitesimal version of this geometric fact.

🔵 Historical Note: The algebraic approach to angular momentum was pioneered independently by Werner Heisenberg and Wolfgang Pauli in 1925-1926, in the context of matrix mechanics. The elegant derivation we present here was crystallized by Dirac in his Principles of Quantum Mechanics (1930) and refined by Schwinger and Sakurai. The recognition that the commutation relations are the Lie algebra of $SU(2)$ came from the mathematical physics community, particularly Weyl and Wigner.

12.2 The Angular Momentum Commutation Relations: $[\hat{J}_x, \hat{J}_y] = i\hbar\hat{J}_z$

Cyclic Structure

The three commutation relations have a beautiful cyclic structure. If we label the components $(1, 2, 3) \leftrightarrow (x, y, z)$, then:

$$[\hat{J}_1, \hat{J}_2] = i\hbar \hat{J}_3$$

and every relation is obtained by cyclic permutation: $1 \to 2 \to 3 \to 1$. This means once you know one relation, you know all three — the algebra is completely specified by its structure constants, which are $i\hbar \epsilon_{ijk}$.

It is worth writing out all three explicitly, because they are used so frequently that you should know them by heart:

$$[\hat{J}_x, \hat{J}_y] = i\hbar \hat{J}_z \qquad [\hat{J}_y, \hat{J}_z] = i\hbar \hat{J}_x \qquad [\hat{J}_z, \hat{J}_x] = i\hbar \hat{J}_y$$

Note the cyclic pattern: in each relation, the operator on the right is the "next" component in the cycle $x \to y \to z \to x$. If you reverse the order of any commutator, the sign flips: $[\hat{J}_y, \hat{J}_x] = -i\hbar \hat{J}_z$. This antisymmetry of the commutator is what makes the Levi-Civita symbol $\epsilon_{ijk}$ the natural language for expressing these relations.

📊 By the Numbers: The angular momentum commutation relations appear in virtually every branch of quantum physics. A survey of published physics papers would find them cited more frequently than almost any other algebraic relation in quantum mechanics, rivaled perhaps only by the canonical commutation relation $[\hat{x}, \hat{p}] = i\hbar$. Both encode fundamental features of the quantum world: the latter encodes the uncertainty principle for position and momentum, the former encodes the non-commutativity of rotations.

The Casimir Operator $\hat{J}^2$

Define the total angular momentum squared:

$$\hat{J}^2 \equiv \hat{J}_x^2 + \hat{J}_y^2 + \hat{J}_z^2$$

This is the Casimir operator of the angular momentum algebra. It has a remarkable property: it commutes with every component of $\hat{\mathbf{J}}$.

Proof that $[\hat{J}^2, \hat{J}_z] = 0$:

$$[\hat{J}^2, \hat{J}_z] = [\hat{J}_x^2, \hat{J}_z] + [\hat{J}_y^2, \hat{J}_z] + [\hat{J}_z^2, \hat{J}_z]$$

The last term vanishes (any operator commutes with itself). For the first term, we use the identity $[\hat{A}^2, \hat{B}] = \hat{A}[\hat{A}, \hat{B}] + [\hat{A}, \hat{B}]\hat{A}$:

$$[\hat{J}_x^2, \hat{J}_z] = \hat{J}_x[\hat{J}_x, \hat{J}_z] + [\hat{J}_x, \hat{J}_z]\hat{J}_x = \hat{J}_x(-i\hbar \hat{J}_y) + (-i\hbar \hat{J}_y)\hat{J}_x = -i\hbar(\hat{J}_x\hat{J}_y + \hat{J}_y\hat{J}_x)$$

Similarly:

$$[\hat{J}_y^2, \hat{J}_z] = \hat{J}_y[\hat{J}_y, \hat{J}_z] + [\hat{J}_y, \hat{J}_z]\hat{J}_y = \hat{J}_y(i\hbar \hat{J}_x) + (i\hbar \hat{J}_x)\hat{J}_y = i\hbar(\hat{J}_y\hat{J}_x + \hat{J}_x\hat{J}_y)$$

Adding these:

$$[\hat{J}^2, \hat{J}_z] = -i\hbar(\hat{J}_x\hat{J}_y + \hat{J}_y\hat{J}_x) + i\hbar(\hat{J}_y\hat{J}_x + \hat{J}_x\hat{J}_y) = 0 \quad \checkmark$$

By the cyclic symmetry, $[\hat{J}^2, \hat{J}_x] = 0$ and $[\hat{J}^2, \hat{J}_y] = 0$ as well. Therefore:

$$\boxed{[\hat{J}^2, \hat{J}_i] = 0 \quad \text{for all } i = x, y, z}$$

💡 Key Insight: $\hat{J}^2$ commutes with every component of $\hat{\mathbf{J}}$, but the components do not commute with each other. This means we can find simultaneous eigenstates of $\hat{J}^2$ and one component (conventionally $\hat{J}_z$), but not of two different components. The choice of $\hat{J}_z$ is conventional — any component would work equally well, and the physics cannot depend on which axis we call $z$.

Incompatibility of Different Components

Since $[\hat{J}_x, \hat{J}_y] = i\hbar \hat{J}_z \neq 0$ (in general), $\hat{J}_x$ and $\hat{J}_y$ cannot be simultaneously diagonalized. This is the angular momentum uncertainty relation:

$$\Delta J_x \cdot \Delta J_y \geq \frac{\hbar}{2} |\langle \hat{J}_z \rangle|$$

If a state has definite $\hat{J}_z$ with $m \neq 0$, then $\hat{J}_x$ and $\hat{J}_y$ must be uncertain. This is not a limitation of our measurement apparatus — it is a fundamental property of the algebra.

✅ Checkpoint: Verify for yourself that $[\hat{J}^2, \hat{J}_x] = 0$ by repeating the calculation above with $\hat{J}_z$ replaced by $\hat{J}_x$. The cyclic symmetry should make the algebra straightforward.

The Maximum Set of Commuting Observables

The operators $\hat{J}^2$ and $\hat{J}_z$ form a complete set of commuting observables (CSCO) for the angular momentum algebra. This means:

1. They commute with each other: $[\hat{J}^2, \hat{J}_z] = 0$.
2. Their simultaneous eigenvalues uniquely label every state within a given angular momentum subspace.
3. Adding any other component of $\hat{\mathbf{J}}$ to the CSCO is impossible, because $\hat{J}_x$ and $\hat{J}_y$ do not commute with $\hat{J}_z$.

The choice of $\hat{J}_z$ is conventional — we could equally well use $\hat{J}_x$ or $\hat{J}_y$ (or any component along an arbitrary axis $\hat{n}$). However, once we choose a quantization axis, the other components become genuinely uncertain. This is one of the most important differences between classical and quantum angular momentum: classically, all three components of $\mathbf{L}$ are simultaneously well-defined; quantum mechanically, only one can be sharp at a time.

The reason $\hat{J}_z$ is the conventional choice is historical and practical. In atomic physics, the quantization axis is typically defined by an external magnetic field (which conventionally points along $z$). The Zeeman effect — the splitting of spectral lines in a magnetic field — directly measures $m$-values, the eigenvalues of $\hat{J}_z$. But the physics is rotationally invariant until an external field breaks the symmetry.

🔵 Historical Note: The labels "magnetic quantum number" for $m$ and "azimuthal quantum number" for $l$ (or $j$) come from the Zeeman effect. Pieter Zeeman discovered the effect in 1896, and Hendrik Lorentz explained it classically. The quantum explanation required the full angular momentum algebra — and the anomalous Zeeman effect, which splits lines into more components than the classical theory predicts, provided crucial evidence for electron spin.

12.3 $\hat{J}^2$ and $\hat{J}_z$: Simultaneous Eigenstates $|j,m\rangle$

Setting Up the Eigenvalue Problem

Since $\hat{J}^2$ and $\hat{J}_z$ commute, they can be simultaneously diagonalized. We write their simultaneous eigenstates as $|j, m\rangle$ (we will determine the allowed values of $j$ and $m$ shortly):

$$\hat{J}^2 |j, m\rangle = \lambda_{j} |j, m\rangle$$

$$\hat{J}_z |j, m\rangle = m\hbar |j, m\rangle$$

We have factored out one power of $\hbar$ in the $\hat{J}_z$ eigenvalue for dimensional reasons ($\hat{J}_z$ has dimensions of angular momentum = $\hbar$, so $m$ is dimensionless). For $\hat{J}^2$, the eigenvalue $\lambda_j$ has dimensions of $\hbar^2$; we will shortly show it takes the form $j(j+1)\hbar^2$, but we do not assume this yet.

We require these states to be orthonormal:

$$\langle j', m' | j, m \rangle = \delta_{j'j} \delta_{m'm}$$

Hermiticity Constraints

Since all three $\hat{J}_i$ are Hermitian, $\hat{J}^2$ is Hermitian as well (sum of squares of Hermitian operators). Its eigenvalues $\lambda_j$ must be real and non-negative:

$$\lambda_j = \langle j, m | \hat{J}^2 | j, m \rangle = \langle j, m | (\hat{J}_x^2 + \hat{J}_y^2 + \hat{J}_z^2) | j, m \rangle$$

Each term $\langle j, m | \hat{J}_i^2 | j, m \rangle = \| \hat{J}_i |j, m\rangle \|^2 \geq 0$, so:

$$\lambda_j \geq 0$$

Furthermore, since $\hat{J}^2 = \hat{J}_x^2 + \hat{J}_y^2 + \hat{J}_z^2 \geq \hat{J}_z^2$, we have:

$$\lambda_j \geq m^2 \hbar^2$$

This critical inequality tells us that for a given $\lambda_j$, the quantum number $m$ is bounded. It cannot grow without limit. This boundedness is the seed from which the entire spectrum grows.

🔗 Connection (Ch 8): The argument that eigenvalues of $\hat{A}^\dagger \hat{A}$ are non-negative was established in Chapter 8. Here we use the same logic: $\hat{J}_i^2 = \hat{J}_i^\dagger \hat{J}_i$ because $\hat{J}_i$ is Hermitian, so its expectation value is the squared norm of $\hat{J}_i |j,m\rangle$.

A Note on Natural Units

For the derivations in this chapter, it is convenient to work in natural units where $\hbar = 1$. In these units, angular momentum is dimensionless, and the eigenvalue equations become:

$$\hat{J}^2 |j, m\rangle = \lambda_j |j, m\rangle, \quad \hat{J}_z |j, m\rangle = m |j, m\rangle$$

We can always restore factors of $\hbar$ at the end by dimensional analysis: wherever $m$ appears as an eigenvalue, multiply by $\hbar$; wherever $\lambda_j$ appears, multiply by $\hbar^2$. We will switch back and forth as needed, always stating clearly which convention is in use.

12.4 Raising and Lowering Operators: $\hat{J}_{\pm} = \hat{J}_x \pm i\hat{J}_y$

Construction

Define the raising and lowering operators (collectively, ladder operators):

$$\boxed{\hat{J}_+ \equiv \hat{J}_x + i\hat{J}_y, \qquad \hat{J}_- \equiv \hat{J}_x - i\hat{J}_y}$$

These are not Hermitian — they are adjoints of each other:

$$\hat{J}_+^\dagger = \hat{J}_-$$

The names "raising" and "lowering" will be justified momentarily.

🔄 Retrieval Practice (Ch 8): Recall the ladder operators $\hat{a}$ and $\hat{a}^\dagger$ for the harmonic oscillator from Chapter 8. The construction here is directly analogous: we combine two non-commuting Hermitian operators into a pair of non-Hermitian operators that act as "steps" on a ladder of eigenstates. The harmonic oscillator ladder had no upper bound; the angular momentum ladder, as we shall see, has both an upper and a lower bound.

Key Commutation Relations

Let us compute the commutators that make $\hat{J}_\pm$ so powerful.

$\hat{J}_\pm$ with $\hat{J}_z$:

$$[\hat{J}_z, \hat{J}_+] = [\hat{J}_z, \hat{J}_x + i\hat{J}_y] = [\hat{J}_z, \hat{J}_x] + i[\hat{J}_z, \hat{J}_y]$$

Using $[\hat{J}_z, \hat{J}_x] = i\hbar \hat{J}_y$ and $[\hat{J}_z, \hat{J}_y] = -i\hbar \hat{J}_x$:

$$[\hat{J}_z, \hat{J}_+] = i\hbar \hat{J}_y + i(-i\hbar \hat{J}_x) = i\hbar \hat{J}_y + \hbar \hat{J}_x = \hbar(\hat{J}_x + i\hat{J}_y) = \hbar \hat{J}_+$$

Similarly:

$$[\hat{J}_z, \hat{J}_-] = -\hbar \hat{J}_-$$

These two results, boxed for emphasis, are:

$$\boxed{[\hat{J}_z, \hat{J}_+] = +\hbar \hat{J}_+, \qquad [\hat{J}_z, \hat{J}_-] = -\hbar \hat{J}_-}$$

$\hat{J}_\pm$ with $\hat{J}^2$:

Since $[\hat{J}^2, \hat{J}_x] = 0$ and $[\hat{J}^2, \hat{J}_y] = 0$:

$$[\hat{J}^2, \hat{J}_\pm] = 0$$

$\hat{J}_+$ with $\hat{J}_-$:

$$[\hat{J}_+, \hat{J}_-] = [\hat{J}_x + i\hat{J}_y, \hat{J}_x - i\hat{J}_y] = -i[\hat{J}_x, \hat{J}_y] + i[\hat{J}_y, \hat{J}_x] = -2i[\hat{J}_x, \hat{J}_y] = -2i(i\hbar \hat{J}_z) = 2\hbar \hat{J}_z$$

$$\boxed{[\hat{J}_+, \hat{J}_-] = 2\hbar \hat{J}_z}$$

The Ladder Property

Now we prove that $\hat{J}_+$ raises and $\hat{J}_-$ lowers the eigenvalue $m$.

Theorem: If $|j, m\rangle$ is an eigenstate of $\hat{J}_z$ with eigenvalue $m\hbar$, then $\hat{J}_+ |j, m\rangle$ is an eigenstate of $\hat{J}_z$ with eigenvalue $(m+1)\hbar$ (or it is the zero vector).

Proof: We want to compute $\hat{J}_z (\hat{J}_+ |j, m\rangle)$. Using $[\hat{J}_z, \hat{J}_+] = \hbar \hat{J}_+$, which means $\hat{J}_z \hat{J}_+ = \hat{J}_+ \hat{J}_z + \hbar \hat{J}_+$:

$$\hat{J}_z (\hat{J}_+ |j, m\rangle) = (\hat{J}_+ \hat{J}_z + \hbar \hat{J}_+) |j, m\rangle = \hat{J}_+ (m\hbar |j, m\rangle) + \hbar \hat{J}_+ |j, m\rangle = (m+1)\hbar (\hat{J}_+ |j, m\rangle)$$

Therefore $\hat{J}_+ |j, m\rangle$ has $\hat{J}_z$ eigenvalue $(m+1)\hbar$. Similarly, since $[\hat{J}^2, \hat{J}_+] = 0$:

$$\hat{J}^2 (\hat{J}_+ |j, m\rangle) = \hat{J}_+ (\hat{J}^2 |j, m\rangle) = \lambda_j (\hat{J}_+ |j, m\rangle)$$

So $\hat{J}_+ |j, m\rangle$ has the same $\hat{J}^2$ eigenvalue. It moves us up the $m$-ladder while keeping us within the same $j$-multiplet:

$$\hat{J}_+ |j, m\rangle = c_+^{j,m} |j, m+1\rangle$$

where $c_+^{j,m}$ is a normalization constant (possibly zero). Analogously:

$$\hat{J}_- |j, m\rangle = c_-^{j,m} |j, m-1\rangle$$

💡 Key Insight: The ladder operators $\hat{J}_\pm$ change the $\hat{J}_z$ eigenvalue by exactly $\pm \hbar$ without changing the $\hat{J}^2$ eigenvalue. They navigate within a fixed $j$-multiplet, stepping up or down the rungs of the $m$-ladder one unit at a time.

Useful Identities

Before proceeding, let us establish identities that will be needed for the spectrum derivation. Express $\hat{J}^2$ in terms of ladder operators:

$$\hat{J}_+ \hat{J}_- = (\hat{J}_x + i\hat{J}_y)(\hat{J}_x - i\hat{J}_y) = \hat{J}_x^2 + \hat{J}_y^2 - i[\hat{J}_x, \hat{J}_y] = \hat{J}_x^2 + \hat{J}_y^2 + \hbar \hat{J}_z$$

Therefore:

$$\hat{J}^2 = \hat{J}_+ \hat{J}_- + \hat{J}_z^2 - \hbar \hat{J}_z$$

Similarly:

$$\hat{J}^2 = \hat{J}_- \hat{J}_+ + \hat{J}_z^2 + \hbar \hat{J}_z$$

These two forms will both be used momentarily.

✅ Checkpoint: Verify the identity $\hat{J}_- \hat{J}_+ = \hat{J}^2 - \hat{J}_z^2 - \hbar \hat{J}_z$ by expanding $\hat{J}_- \hat{J}_+ = (\hat{J}_x - i\hat{J}_y)(\hat{J}_x + i\hat{J}_y)$ directly.

12.5 The Eigenvalue Spectrum: $j = 0, \frac{1}{2}, 1, \frac{3}{2}, \ldots$ and $m = -j$ to $+j$

This is the crown jewel of the chapter — arguably the single most beautiful algebraic derivation in quantum mechanics. We will prove that the commutation relations alone determine the complete eigenvalue spectrum.

The Argument

Step 1: The ladder must terminate.

For a given $\hat{J}^2$ eigenvalue $\lambda_j$, we established that $\lambda_j \geq m^2 \hbar^2$. This means $m$ is bounded: $|m| \leq \sqrt{\lambda_j}/\hbar$. Since $\hat{J}_+$ raises $m$ by 1 and $\hat{J}_-$ lowers $m$ by 1, the ladder of $m$-values must terminate in both directions. There exists a maximum value $m_{\max}$ and a minimum value $m_{\min}$.

At the top of the ladder:

$$\hat{J}_+ |j, m_{\max}\rangle = 0$$

At the bottom:

$$\hat{J}_- |j, m_{\min}\rangle = 0$$

Step 2: Determine the eigenvalue from the top.

Apply $\hat{J}_- \hat{J}_+$ to $|j, m_{\max}\rangle$:

$$\hat{J}_- \hat{J}_+ |j, m_{\max}\rangle = 0$$

Using the identity $\hat{J}_- \hat{J}_+ = \hat{J}^2 - \hat{J}_z^2 - \hbar \hat{J}_z$:

$$(\hat{J}^2 - \hat{J}_z^2 - \hbar \hat{J}_z) |j, m_{\max}\rangle = 0$$

$$(\lambda_j - m_{\max}^2 \hbar^2 - m_{\max} \hbar^2) |j, m_{\max}\rangle = 0$$

Since $|j, m_{\max}\rangle \neq 0$:

$$\lambda_j = \hbar^2 m_{\max}(m_{\max} + 1)$$

Step 3: Determine the eigenvalue from the bottom.

Apply $\hat{J}_+ \hat{J}_-$ to $|j, m_{\min}\rangle$:

$$\hat{J}_+ \hat{J}_- |j, m_{\min}\rangle = 0$$

Using $\hat{J}_+ \hat{J}_- = \hat{J}^2 - \hat{J}_z^2 + \hbar \hat{J}_z$:

$$(\lambda_j - m_{\min}^2 \hbar^2 + m_{\min} \hbar^2) |j, m_{\min}\rangle = 0$$

$$\lambda_j = \hbar^2 m_{\min}(m_{\min} - 1)$$

Step 4: Equate the two expressions.

$$m_{\max}(m_{\max} + 1) = m_{\min}(m_{\min} - 1)$$

Rearranging:

$$m_{\max}^2 + m_{\max} = m_{\min}^2 - m_{\min}$$

$$m_{\max}^2 - m_{\min}^2 + m_{\max} + m_{\min} = 0$$

$$(m_{\max} - m_{\min})(m_{\max} + m_{\min}) + (m_{\max} + m_{\min}) = 0$$

$$(m_{\max} + m_{\min})(m_{\max} - m_{\min} + 1) = 0$$

The second factor $m_{\max} - m_{\min} + 1 > 0$ (since $m_{\max} \geq m_{\min}$), so:

$$\boxed{m_{\min} = -m_{\max}}$$

The ladder is symmetric about zero!

Step 5: Identify the quantum numbers.

Define $j \equiv m_{\max}$. Then $m_{\min} = -j$, the eigenvalue is:

$$\boxed{\lambda_j = j(j+1)\hbar^2}$$

and $m$ ranges from $-j$ to $+j$ in integer steps (since each application of $\hat{J}_+$ increases $m$ by 1). The number of steps from $m_{\min} = -j$ to $m_{\max} = j$ is $2j$, which must be a non-negative integer. Therefore:

$$\boxed{j = 0, \frac{1}{2}, 1, \frac{3}{2}, 2, \frac{5}{2}, \ldots}$$

For each $j$, there are $2j + 1$ values of $m$:

$$m = -j, -j+1, -j+2, \ldots, j-1, j$$

This is the complete eigenvalue spectrum of angular momentum. Let us write it in full glory:

$$\hat{J}^2 |j, m\rangle = j(j+1)\hbar^2 |j, m\rangle$$

$$\hat{J}_z |j, m\rangle = m\hbar |j, m\rangle$$

$$j = 0, \tfrac{1}{2}, 1, \tfrac{3}{2}, 2, \ldots; \qquad m = -j, -j+1, \ldots, j-1, j$$

💡 Key Insight: The half-integer values $j = 1/2, 3/2, 5/2, \ldots$ emerge automatically from the algebra. We did not put them in by hand. The commutation relations demanded them. This is remarkable: pure algebra predicts the existence of particles with spin-1/2 (electrons, protons, neutrons, quarks) and spin-3/2 (the $\Delta$ baryons). The mathematics knew about these particles before physicists did.

⚠️ Common Misconception: "The eigenvalue of $\hat{J}^2$ is $j^2 \hbar^2$." No! It is $j(j+1)\hbar^2$. This is always strictly greater than $(m\hbar)^2_{\max} = j^2 \hbar^2$. The angular momentum vector can never be perfectly aligned along any axis — there is always some transverse component. For $j = 1/2$: $\hat{J}^2$ has eigenvalue $\frac{3}{4}\hbar^2$, while the maximum $\hat{J}_z$ is only $\frac{1}{2}\hbar$. The "length" of the angular momentum vector is $\sqrt{3/4}\,\hbar \approx 0.87\hbar$, but its maximum projection is $0.5\hbar$.

Normalization Constants for the Ladder Operators

We need the constants $c_\pm^{j,m}$ in $\hat{J}_\pm |j, m\rangle = c_\pm^{j,m} |j, m \pm 1\rangle$.

Compute $\| \hat{J}_+ |j, m\rangle \|^2$:

$$\langle j, m | \hat{J}_-\hat{J}_+ | j, m \rangle = \langle j, m | (\hat{J}^2 - \hat{J}_z^2 - \hbar \hat{J}_z) | j, m \rangle = [j(j+1) - m^2 - m]\hbar^2 = [j(j+1) - m(m+1)]\hbar^2$$

This must equal $|c_+^{j,m}|^2$, so choosing the conventional positive real phase:

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

Similarly:

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

Sanity checks:

• $\hat{J}_+ |j, j\rangle = \hbar\sqrt{(j-j)(j+j+1)} |j, j+1\rangle = 0$ ✓ (top of ladder)
• $\hat{J}_- |j, -j\rangle = \hbar\sqrt{(j-j)(-j+j+1)} |j, -j-1\rangle = 0$ ✓ (bottom of ladder)
• $\hat{J}_+ |1/2, -1/2\rangle = \hbar\sqrt{(1)(1)} |1/2, 1/2\rangle = \hbar |1/2, 1/2\rangle$ ✓

✅ Checkpoint: Verify that $\hat{J}_- |1, 1\rangle = \hbar\sqrt{2} |1, 0\rangle$ and $\hat{J}_- |1, 0\rangle = \hbar\sqrt{2} |1, -1\rangle$ using the formula above.

Worked Example: Building the $j = 1$ Multiplet from the Top

Let us construct all three states $|1, 1\rangle$, $|1, 0\rangle$, $|1, -1\rangle$ starting from the top of the ladder, to illustrate the procedure concretely.

Starting point: $|1, 1\rangle$ is given (the top state). We verify: $\hat{J}_+ |1, 1\rangle = \hbar\sqrt{(1-1)(1+1+1)} |1, 2\rangle = 0$. Good — we are at the top.

Step down to $|1, 0\rangle$:

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

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

Step down to $|1, -1\rangle$:

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

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

Verify the bottom: $\hat{J}_- |1, -1\rangle = \hbar\sqrt{(1-1)(1+1+1)} |1, -2\rangle = 0$. We have reached the bottom of the ladder.

Consistency check: The number of states is $2(1) + 1 = 3$. The $m$-values are $\{1, 0, -1\}$, as expected.

This procedure generalizes: for any $j$, start from $|j, j\rangle$ and apply $\hat{J}_-$ repeatedly, collecting the normalization constants. After exactly $2j$ steps, you reach $|j, -j\rangle$, and one more application of $\hat{J}_-$ gives zero.

Expectation Values in the State $|j, m\rangle$

Several expectation values in the eigenstates $|j, m\rangle$ follow immediately from the algebra and will be used repeatedly in subsequent chapters:

$$\langle j, m | \hat{J}_z | j, m \rangle = m\hbar$$

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

$$\langle j, m | \hat{J}_x | j, m \rangle = 0, \qquad \langle j, m | \hat{J}_y | j, m \rangle = 0$$

The last two follow because $\hat{J}_x = (\hat{J}_+ + \hat{J}_-)/2$ and $\hat{J}_y = (\hat{J}_+ - \hat{J}_-)/2i$, and $\hat{J}_\pm$ connects $|j, m\rangle$ only to $|j, m \pm 1\rangle$, which is orthogonal to $|j, m\rangle$. So the diagonal matrix elements of $\hat{J}_x$ and $\hat{J}_y$ vanish.

For the second moments:

$$\langle j, m | \hat{J}_x^2 | j, m \rangle = \langle j, m | \hat{J}_y^2 | j, m \rangle = \frac{1}{2}[j(j+1) - m^2]\hbar^2$$

This follows from $\hat{J}_x^2 + \hat{J}_y^2 = \hat{J}^2 - \hat{J}_z^2$ and the symmetry $\langle \hat{J}_x^2 \rangle = \langle \hat{J}_y^2 \rangle$ (which holds because $|j, m\rangle$ is an eigenstate of $\hat{J}_z$, giving no preferred direction in the $xy$-plane).

The uncertainties are therefore:

$$\Delta J_x = \Delta J_y = \hbar\sqrt{\frac{j(j+1) - m^2}{2}}$$

For the state $|j, j\rangle$: $\Delta J_x = \Delta J_y = \hbar\sqrt{j/2}$. The angular momentum "cone" has a nonzero opening angle that decreases as $j$ increases — this is the angular momentum version of the correspondence principle.

12.6 Matrix Representations of Angular Momentum

Now we convert the abstract algebra into explicit matrices. For a given $j$, the Hilbert space spanned by $\{|j, m\rangle : m = -j, \ldots, j\}$ has dimension $2j+1$. In this basis, every angular momentum operator becomes a $(2j+1) \times (2j+1)$ matrix.

🔄 Retrieval Practice (Ch 8): Recall from Chapter 8 that the matrix element of an operator $\hat{A}$ in a basis $\{|n\rangle\}$ is $A_{mn} = \langle m | \hat{A} | n \rangle$. We will use the basis ordering $|j, j\rangle, |j, j-1\rangle, \ldots, |j, -j\rangle$ (descending $m$, the standard convention in most textbooks).

The $j = 1/2$ Representation

For $j = 1/2$, the basis is $\{|1/2, 1/2\rangle, |1/2, -1/2\rangle\}$, which we abbreviate as $\{|+\rangle, |-\rangle\}$ or $\{|\uparrow\rangle, |\downarrow\rangle\}$. The matrices are $2 \times 2$.

$\hat{J}_z$:

$$(\hat{J}_z)_{m'm} = \langle j, m' | \hat{J}_z | j, m \rangle = m\hbar \,\delta_{m'm}$$

$$\hat{J}_z = \frac{\hbar}{2} \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}$$

$\hat{J}_+$:

The only nonzero matrix element is $\langle +|\hat{J}_+|-\rangle = \hbar\sqrt{(1/2+1/2)(1/2-1/2+1)} = \hbar$:

$$\hat{J}_+ = \hbar \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}$$

$\hat{J}_-$:

$$\hat{J}_- = \hbar \begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix}$$

$\hat{J}_x$ and $\hat{J}_y$ from $\hat{J}_x = (\hat{J}_+ + \hat{J}_-)/2$ and $\hat{J}_y = (\hat{J}_+ - \hat{J}_-)/2i$:

$$\hat{J}_x = \frac{\hbar}{2} \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}, \qquad \hat{J}_y = \frac{\hbar}{2} \begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}$$

These are the Pauli matrices (up to the factor $\hbar/2$): $\hat{J}_i = (\hbar/2)\sigma_i$, where:

$$\sigma_x = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}, \quad \sigma_y = \begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}, \quad \sigma_z = \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}$$

🔗 Connection (Ch 13 preview): These Pauli matrices will be the central objects of Chapter 13, where we study spin-1/2 in full detail. Here we note simply that they emerged inevitably from the angular momentum algebra applied to the case $j = 1/2$.

$\hat{J}^2$:

$$\hat{J}^2 = j(j+1)\hbar^2 \hat{I} = \frac{3}{4}\hbar^2 \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}$$

$\hat{J}^2$ is proportional to the identity matrix within a fixed-$j$ subspace — this is a general feature, not special to $j = 1/2$.

The $j = 1$ Representation

For $j = 1$, the basis is $\{|1, 1\rangle, |1, 0\rangle, |1, -1\rangle\}$. The matrices are $3 \times 3$.

$\hat{J}_z$:

$$\hat{J}_z = \hbar \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & -1 \end{pmatrix}$$

$\hat{J}_+$: Using $\hat{J}_+ |1, m\rangle = \hbar\sqrt{(1-m)(2+m)} |1, m+1\rangle$:

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

$$\hat{J}_+ = \hbar\sqrt{2} \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix}$$

$\hat{J}_-$:

$$\hat{J}_- = \hbar\sqrt{2} \begin{pmatrix} 0 & 0 & 0 \\ 1 & 0 & 0 \\ 0 & 1 & 0 \end{pmatrix}$$

$\hat{J}_x$ and $\hat{J}_y$:

$$\hat{J}_x = \frac{\hbar}{\sqrt{2}} \begin{pmatrix} 0 & 1 & 0 \\ 1 & 0 & 1 \\ 0 & 1 & 0 \end{pmatrix}, \qquad \hat{J}_y = \frac{\hbar}{\sqrt{2}} \begin{pmatrix} 0 & -i & 0 \\ i & 0 & -i \\ 0 & i & 0 \end{pmatrix}$$

$\hat{J}^2$:

$$\hat{J}^2 = 1(1+1)\hbar^2 \hat{I} = 2\hbar^2 \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix}$$

✅ Checkpoint: Verify that $[\hat{J}_x, \hat{J}_y] = i\hbar \hat{J}_z$ using the $j = 1$ matrices above. This is a matrix multiplication exercise that confirms the entire algebraic structure is consistently represented.

Worked Example: The $j = 3/2$ Representation

For $j = 3/2$, the basis is $\{|3/2, 3/2\rangle, |3/2, 1/2\rangle, |3/2, -1/2\rangle, |3/2, -3/2\rangle\}$, giving $4 \times 4$ matrices. Let us construct $\hat{J}_+$.

We need $\hat{J}_+ |3/2, m\rangle = \hbar\sqrt{(3/2 - m)(3/2 + m + 1)} |3/2, m+1\rangle$:

• $m = 1/2$: $\hbar\sqrt{(1)(3)} = \hbar\sqrt{3}$
• $m = -1/2$: $\hbar\sqrt{(2)(2)} = 2\hbar$
• $m = -3/2$: $\hbar\sqrt{(3)(1)} = \hbar\sqrt{3}$
• $m = 3/2$: $0$

$$\hat{J}_+ = \hbar \begin{pmatrix} 0 & \sqrt{3} & 0 & 0 \\ 0 & 0 & 2 & 0 \\ 0 & 0 & 0 & \sqrt{3} \\ 0 & 0 & 0 & 0 \end{pmatrix}$$

$$\hat{J}_z = \frac{\hbar}{2} \begin{pmatrix} 3 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & -1 & 0 \\ 0 & 0 & 0 & -3 \end{pmatrix}$$

$$\hat{J}_x = \frac{\hbar}{2} \begin{pmatrix} 0 & \sqrt{3} & 0 & 0 \\ \sqrt{3} & 0 & 2 & 0 \\ 0 & 2 & 0 & \sqrt{3} \\ 0 & 0 & \sqrt{3} & 0 \end{pmatrix}$$

Note the pattern: the matrix elements of $\hat{J}_+$ appear on the superdiagonal, those of $\hat{J}_-$ on the subdiagonal, and $\hat{J}_z$ is always diagonal. This structure persists for all $j$.

General Pattern

For arbitrary $j$, the matrices are $(2j+1) \times (2j+1)$ with:

• $(\hat{J}_z)_{m'm} = m\hbar \, \delta_{m'm}$ — diagonal
• $(\hat{J}_+)_{m',m} = \hbar\sqrt{(j-m)(j+m+1)} \, \delta_{m',m+1}$ — superdiagonal
• $(\hat{J}_-)_{m',m} = \hbar\sqrt{(j+m)(j-m+1)} \, \delta_{m',m-1}$ — subdiagonal

These formulas, combined with $\hat{J}_x = (\hat{J}_+ + \hat{J}_-)/2$ and $\hat{J}_y = (\hat{J}_+ - \hat{J}_-)/2i$, give explicit matrix representations for any angular momentum.

📊 By the Numbers: The dimension of the representation grows linearly with $j$: dimension $= 2j + 1$. For $j = 50$ (relevant in molecular spectroscopy), the matrices are $101 \times 101$. For $j = 100$ (relevant in nuclear physics of high-spin states), they are $201 \times 201$. Even though the formulas are simple, practical calculations with large $j$ are best done computationally — which is why the Python toolkit (Section 12.10) is essential.

12.7 Rotation Matrices and the Wigner $D$-Functions

Angular momentum operators are generators of rotations. Let us now construct the rotation operators and their matrix representations.

Rotation Operators

A rotation by angle $\phi$ about the axis $\hat{n}$ is represented by the unitary operator:

$$\hat{R}(\hat{n}, \phi) = \exp\left(-\frac{i\phi}{\hbar} \hat{n} \cdot \hat{\mathbf{J}}\right) = \exp\left(-\frac{i\phi}{\hbar}(n_x \hat{J}_x + n_y \hat{J}_y + n_z \hat{J}_z)\right)$$

This is a unitary operator because $\hat{\mathbf{J}}$ is Hermitian. It transforms kets as $|j, m\rangle \to \hat{R} |j, m\rangle$ and operators as $\hat{A} \to \hat{R}\hat{A}\hat{R}^\dagger$.

🔄 Retrieval Practice (Ch 10): In Chapter 10, we established that the generator of spatial translations is momentum $\hat{p}$, with the translation operator $\hat{T}(a) = e^{-ia\hat{p}/\hbar}$. The rotation operator has exactly the same structure, with angular momentum replacing linear momentum and angle replacing displacement. Symmetry generators and their associated unitary transformations always have this exponential form.

Euler Angle Parameterization

Any rotation in 3D can be decomposed into three successive rotations using Euler angles $(\alpha, \beta, \gamma)$:

1. Rotate by $\alpha$ about the $z$-axis
2. Rotate by $\beta$ about the (new) $y$-axis
3. Rotate by $\gamma$ about the (new) $z$-axis

The corresponding operator is:

$$\hat{R}(\alpha, \beta, \gamma) = e^{-i\alpha \hat{J}_z/\hbar} \, e^{-i\beta \hat{J}_y/\hbar} \, e^{-i\gamma \hat{J}_z/\hbar}$$

The ranges are: $0 \leq \alpha < 2\pi$, $0 \leq \beta \leq \pi$, $0 \leq \gamma < 2\pi$.

The Wigner $D$-Matrix

The matrix elements of the rotation operator in the $|j, m\rangle$ basis are the Wigner $D$-functions:

$$D^{(j)}_{m'm}(\alpha, \beta, \gamma) \equiv \langle j, m' | \hat{R}(\alpha, \beta, \gamma) | j, m \rangle$$

Since $\hat{J}_z |j, m\rangle = m\hbar |j, m\rangle$, the $z$-rotations contribute only phases:

$$D^{(j)}_{m'm}(\alpha, \beta, \gamma) = e^{-im'\alpha} \, d^{(j)}_{m'm}(\beta) \, e^{-im\gamma}$$

where the reduced rotation matrix (or small $d$-matrix) is:

$$d^{(j)}_{m'm}(\beta) = \langle j, m' | e^{-i\beta \hat{J}_y/\hbar} | j, m \rangle$$

The full content of the rotation is in $d^{(j)}_{m'm}(\beta)$, which depends only on the polar angle $\beta$.

Explicit Formulas

For $j = 1/2$:

Using $\hat{J}_y = (\hbar/2)\sigma_y$ and the identity $e^{-i\theta \sigma_y/2} = \cos(\theta/2)\hat{I} - i\sin(\theta/2)\sigma_y$:

$$d^{(1/2)}(\beta) = \begin{pmatrix} \cos(\beta/2) & -\sin(\beta/2) \\ \sin(\beta/2) & \cos(\beta/2) \end{pmatrix}$$

Note the half-angle dependence: a rotation by $\beta = 2\pi$ gives $d^{(1/2)}(2\pi) = -\hat{I}$, not $+\hat{I}$. The state picks up a sign: $|j, m\rangle \to -|j, m\rangle$. A full $4\pi$ rotation is needed to return to the original state. This is the famous spinor sign change under $2\pi$ rotation, which distinguishes half-integer from integer angular momentum and is the hallmark of fermionic behavior.

🧪 Experiment: The $4\pi$ periodicity of spinors was confirmed experimentally by Werner, Colella, Overhauser, and Staudenmann (1975) using neutron interferometry. They rotated neutron spins by a controllable angle using a magnetic field and observed the phase shift. A $2\pi$ rotation shifted the interference pattern by exactly half a fringe, confirming the sign change. This is not merely a mathematical curiosity — it is a measurable physical effect.

For $j = 1$:

$$d^{(1)}(\beta) = \begin{pmatrix} \frac{1+\cos\beta}{2} & -\frac{\sin\beta}{\sqrt{2}} & \frac{1-\cos\beta}{2} \\[6pt] \frac{\sin\beta}{\sqrt{2}} & \cos\beta & -\frac{\sin\beta}{\sqrt{2}} \\[6pt] \frac{1-\cos\beta}{2} & \frac{\sin\beta}{\sqrt{2}} & \frac{1+\cos\beta}{2} \end{pmatrix}$$

Verify: $d^{(1)}(0) = \hat{I}_3$ ✓. $d^{(1)}(2\pi) = +\hat{I}_3$ ✓ (integer $j$ returns to itself under $2\pi$ rotation).

Worked Example: Rotating a Spin-1/2 State

To make the rotation matrices concrete, let us rotate the state $|\uparrow\rangle = |1/2, 1/2\rangle$ by angle $\beta = \pi/2$ about the $y$-axis. This is a common operation in NMR and quantum computing.

The rotation matrix is:

$$d^{(1/2)}(\pi/2) = \begin{pmatrix} \cos(\pi/4) & -\sin(\pi/4) \\ \sin(\pi/4) & \cos(\pi/4) \end{pmatrix} = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 & -1 \\ 1 & 1 \end{pmatrix}$$

Acting on $|\uparrow\rangle = \binom{1}{0}$:

$$d^{(1/2)}(\pi/2) \binom{1}{0} = \frac{1}{\sqrt{2}} \binom{1}{1}$$

The rotated state is $\frac{1}{\sqrt{2}}(|\uparrow\rangle + |\downarrow\rangle)$, which is an eigenstate of $\hat{J}_x$ with eigenvalue $+\hbar/2$. This makes physical sense: a $90°$ rotation about $y$ takes the $z$-axis to the $x$-axis, so a state pointing along $z$ becomes a state pointing along $x$.

The probability of measuring spin-up along $z$ after the rotation is:

$$P(\uparrow) = |d^{(1/2)}_{1/2, 1/2}(\pi/2)|^2 = \cos^2(\pi/4) = 1/2$$

And the probability of spin-down is also $1/2$. The rotation has created an equal superposition.

🧪 Experiment: This rotation is precisely the operation performed by a "$\pi/2$ pulse" in NMR spectroscopy. A short burst of radiofrequency radiation, tuned to the Larmor frequency and applied along the $y$-axis in the rotating frame, rotates the nuclear spin from alignment with the static field ($z$) to the transverse plane. The resulting precession generates the NMR signal detected by the receiver coil.

The General Wigner Formula

For arbitrary $j$, the reduced rotation matrix elements are given by the Wigner formula:

$$d^{(j)}_{m'm}(\beta) = \sum_s \frac{(-1)^{m'-m+s} \sqrt{(j+m')!(j-m')!(j+m)!(j-m)!}}{(j+m-s)! \, s! \, (m'-m+s)! \, (j-m'-s)!} \left(\cos\frac{\beta}{2}\right)^{2j-m'+m-2s} \left(\sin\frac{\beta}{2}\right)^{m'-m+2s}$$

where the sum runs over all integers $s$ for which the factorial arguments are non-negative. This formula is unwieldy for hand calculation but straightforward to implement computationally (see Section 12.10).

Properties of the $D$-Matrices

The Wigner $D$-matrices satisfy several important properties:

1. Unitarity: $\sum_{m''} D^{(j)*}_{m''m'} D^{(j)}_{m''m} = \delta_{m'm}$
2. Group property: $D^{(j)}(R_1) D^{(j)}(R_2) = D^{(j)}(R_1 R_2)$
3. Orthogonality: $\int d\Omega \, D^{(j_1)*}_{m_1' m_1} D^{(j_2)}_{m_2' m_2} = \frac{8\pi^2}{2j_1+1} \delta_{j_1 j_2} \delta_{m_1' m_2'} \delta_{m_1 m_2}$
4. Symmetry relations: $d^{(j)}_{m'm}(\beta) = (-1)^{m'-m} d^{(j)}_{mm'}(\beta) = d^{(j)}_{-m,-m'}(\beta)$

💡 Key Insight: The Wigner $D$-matrices are the irreducible representations of the rotation group. Every possible way a quantum state can transform under rotations is captured by some $D^{(j)}$ matrix. The integer-$j$ representations are "single-valued" ($2\pi$ rotation = identity), while the half-integer-$j$ representations are "double-valued" ($2\pi$ rotation = minus identity; $4\pi$ needed for identity). Both types are genuine, physically realized representations.

Worked Example: Rotation of a Spin-1 State

Consider a spin-1 particle in the state $|1, 1\rangle$ — maximum angular momentum along $z$. Rotate it by $\beta = \pi/2$ about the $y$-axis. What is the probability of each outcome if we then measure $\hat{J}_z$?

The rotated state is:

$$|\psi'\rangle = d^{(1)}(\pi/2) |1, 1\rangle = d^{(1)}_{11}(\pi/2)|1,1\rangle + d^{(1)}_{01}(\pi/2)|1,0\rangle + d^{(1)}_{-1,1}(\pi/2)|1,-1\rangle$$

Using the $d^{(1)}(\pi/2)$ matrix:

$$d^{(1)}_{11}(\pi/2) = \frac{1 + \cos(\pi/2)}{2} = \frac{1}{2}$$

$$d^{(1)}_{01}(\pi/2) = -\frac{\sin(\pi/2)}{\sqrt{2}} = -\frac{1}{\sqrt{2}}$$

$$d^{(1)}_{-1,1}(\pi/2) = \frac{1 - \cos(\pi/2)}{2} = \frac{1}{2}$$

Therefore:

$$|\psi'\rangle = \frac{1}{2}|1, 1\rangle - \frac{1}{\sqrt{2}}|1, 0\rangle + \frac{1}{2}|1, -1\rangle$$

The measurement probabilities are:

$$P(m = 1) = 1/4, \quad P(m = 0) = 1/2, \quad P(m = -1) = 1/4$$

Note that the most probable outcome is $m = 0$, not $m = 1$. This makes physical sense: rotating the angular momentum from $z$ to $x$ (a $90°$ rotation about $y$) means the angular momentum is now along $x$, which is equally likely to be "up" or "down" along $z$, with a peak at zero.

The expectation value is:

$$\langle \hat{J}_z \rangle = (1/4)(1) + (1/2)(0) + (1/4)(-1) = 0$$

as expected — the angular momentum is now perpendicular to $z$.

✅ Checkpoint: Verify that $P(m=1) + P(m=0) + P(m=-1) = 1$. Also verify that $\langle \hat{J}_z^2 \rangle = (1/4)(1) + (1/2)(0) + (1/4)(1) = 1/2$, giving $\Delta J_z = \hbar/\sqrt{2}$.

12.8 Orbital Angular Momentum as a Special Case (Integer $j$ Only)

Recovering Orbital Angular Momentum

The general theory allows $j = 0, 1/2, 1, 3/2, \ldots$. But when we studied orbital angular momentum $\hat{\mathbf{L}} = \hat{\mathbf{r}} \times \hat{\mathbf{p}}$ in Chapter 5, we found only integer values $l = 0, 1, 2, \ldots$. Why?

The reason is that orbital angular momentum is realized on a specific Hilbert space — the space of square-integrable functions of position, $L^2(\mathbb{R}^3)$. In this realization, the angular momentum operators are differential operators:

$$\hat{L}_z = -i\hbar \frac{\partial}{\partial \phi}$$

The eigenvalue equation $\hat{L}_z f = m\hbar f$ has solutions $f(\phi) = e^{im\phi}$. The requirement that $f$ be single-valued under $\phi \to \phi + 2\pi$ demands $e^{2\pi i m} = 1$, which forces $m$ to be an integer. Since $m$ ranges from $-l$ to $l$, $l$ must also be an integer.

⚠️ Common Misconception: "Half-integer angular momentum is impossible because wavefunctions must be single-valued." This is wrong. The correct statement is: orbital angular momentum is integer because the position-space representation requires single-valued functions. But spin angular momentum is not associated with spatial wavefunctions at all — it lives in a separate, internal Hilbert space where the single-valuedness constraint does not apply.

The Connection to Spherical Harmonics

The simultaneous eigenstates of $\hat{L}^2$ and $\hat{L}_z$ in position space are the spherical harmonics $Y_l^m(\theta, \phi)$:

$$\langle \theta, \phi | l, m \rangle = Y_l^m(\theta, \phi)$$

The ladder operators in position space are:

$$\hat{L}_\pm = \hbar e^{\pm i\phi} \left(\pm \frac{\partial}{\partial \theta} + i\cot\theta \frac{\partial}{\partial \phi}\right)$$

The algebraic relation $\hat{L}_+ |l, m\rangle = \hbar\sqrt{(l-m)(l+m+1)} |l, m+1\rangle$ becomes:

$$\hat{L}_+ Y_l^m(\theta, \phi) = \hbar\sqrt{(l-m)(l+m+1)} \, Y_l^{m+1}(\theta, \phi)$$

This provides a practical method for constructing all spherical harmonics: start from $Y_l^l(\theta, \phi)$ (which can be found by solving $\hat{L}_+ Y_l^l = 0$) and apply $\hat{L}_-$ repeatedly to generate $Y_l^{l-1}, Y_l^{l-2}, \ldots, Y_l^{-l}$.

🔗 Connection (Ch 5): In Chapter 5, we derived spherical harmonics by solving the angular part of the Schrödinger equation — a second-order differential equation with boundary conditions. The algebraic approach gives the same results more elegantly: solve one first-order equation ($\hat{L}_+ Y_l^l = 0$), then apply $\hat{L}_-$ repeatedly.

Worked Example: Constructing $Y_1^m$ from the Algebra

Let us derive the $l = 1$ spherical harmonics using the algebraic method, as a concrete illustration of the connection between the abstract algebra and the position-space functions.

Step 1: Find $Y_1^1$ by solving $\hat{L}_+ Y_1^1 = 0$. In position space:

$$\hat{L}_+ = \hbar e^{i\phi}\left(\frac{\partial}{\partial\theta} + i\cot\theta\frac{\partial}{\partial\phi}\right)$$

Writing $Y_1^1 = f(\theta) e^{i\phi}$ (since $\hat{L}_z Y_1^1 = \hbar Y_1^1$ forces the $\phi$-dependence to be $e^{i\phi}$), the equation $\hat{L}_+ Y_1^1 = 0$ becomes:

$$\frac{df}{d\theta} - \cot\theta \cdot f = 0$$

This first-order ODE has the solution $f(\theta) = C \sin\theta$, so:

$$Y_1^1(\theta, \phi) = C \sin\theta \, e^{i\phi}$$

Normalizing over the sphere gives $C = -\sqrt{3/(8\pi)}$ (the conventional Condon-Shortley phase), yielding:

$$Y_1^1(\theta, \phi) = -\sqrt{\frac{3}{8\pi}} \sin\theta \, e^{i\phi}$$

Step 2: Apply $\hat{L}_-$ to get $Y_1^0$:

$$\hat{L}_- Y_1^1 = \hbar\sqrt{(1+1)(1-1+1)} Y_1^0 = \hbar\sqrt{2} Y_1^0$$

Computing $\hat{L}_- Y_1^1$ in position space and dividing by $\hbar\sqrt{2}$ gives:

$$Y_1^0(\theta, \phi) = \sqrt{\frac{3}{4\pi}} \cos\theta$$

Step 3: Apply $\hat{L}_-$ once more to get $Y_1^{-1}$, and verify $\hat{L}_- Y_1^{-1} = 0$.

This procedure is faster than solving the associated Legendre equation directly, because it replaces a second-order ODE with boundary conditions by a first-order ODE plus algebraic stepping.

Comparison Table: Orbital vs. General Angular Momentum

✅ Checkpoint: For orbital angular momentum with $l = 2$, list all five $m$-values. Write the action of $\hat{L}_+$ on each basis state, including the normalization constant. Verify that $\hat{L}_+ |2, 2\rangle = 0$.

12.9 Preview: Spin as the Half-Integer Case

What the Algebra Tells Us

The general theory allows $j = 1/2, 3/2, 5/2, \ldots$, and these cannot be orbital angular momentum. What physical systems carry these quantum numbers?

The answer is spin — intrinsic angular momentum that is built into the particle itself, not associated with motion through space. We will develop spin fully in Chapter 13; here we preview the key features that follow directly from the algebra.

For $j = s = 1/2$ (electrons, protons, neutrons, quarks):

• The Hilbert space is two-dimensional: $\{|\uparrow\rangle, |\downarrow\rangle\} = \{|1/2, +1/2\rangle, |1/2, -1/2\rangle\}$
• The spin operators are $\hat{S}_i = (\hbar/2)\sigma_i$ (Pauli matrices)
• A $2\pi$ rotation produces a sign flip: $\hat{R}(2\pi) = -\hat{I}$
• The magnetic moment is $\boldsymbol{\mu} = -g_s \mu_B \hat{\mathbf{S}}/\hbar$ with $g_s \approx 2$

For $j = s = 1$ (photons, W/Z bosons, $\rho$ mesons):

• The Hilbert space is three-dimensional
• The spin operators are the $3 \times 3$ matrices from Section 12.6
• A $2\pi$ rotation returns the state to itself

For $j = s = 3/2$ ($\Delta$ baryons, $\Omega^-$ hyperon):

• The Hilbert space is four-dimensional
• A $2\pi$ rotation produces a sign flip (half-integer)

The Spin-Statistics Connection

One of the deepest results in theoretical physics — proved rigorously by Pauli in 1940 using relativistic quantum field theory — is that integer-spin particles obey Bose-Einstein statistics (bosons), while half-integer-spin particles obey Fermi-Dirac statistics (fermions). This spin-statistics theorem connects the algebraic property $(-1)^{2j}$ (the sign under $2\pi$ rotation) to the behavior of identical particles (Chapter 15).

💡 Key Insight: The algebraic fact that $j$ can be half-integer is not a mathematical curiosity — it is the foundation of all of chemistry (the Pauli exclusion principle, which depends on electrons being spin-1/2 fermions) and all of nuclear physics (where both half-integer and integer spin particles play essential roles). Remove half-integer angular momentum from the theory, and atoms cannot exist.

The Physical Evidence for Half-Integer Angular Momentum

The algebraic theory predicts half-integer $j$. What is the experimental evidence?

Stern-Gerlach experiment (1922): Silver atoms passed through an inhomogeneous magnetic field split into exactly two beams. For orbital angular momentum $l$, the number of beams would be $2l + 1$, which is always odd. Two beams require $2j + 1 = 2$, i.e., $j = 1/2$. This was the first observation of spin (though it was not understood as such until 1925).

Anomalous Zeeman effect: The splitting of spectral lines in a magnetic field showed patterns that could not be explained by integer angular momentum alone. The inclusion of electron spin ($s = 1/2$) and the Lande $g$-factor resolved all anomalies.

Fine structure of hydrogen: The hydrogen spectrum shows doublet structure (two closely spaced lines where one is expected) that is explained by spin-orbit coupling between the electron's orbital angular momentum ($l$) and its spin ($s = 1/2$). This is discussed in detail in Chapter 18.

The periodic table: The Pauli exclusion principle requires that no two electrons occupy the same quantum state. For a given $(n, l, m_l)$, there are exactly two allowed values of $m_s = \pm 1/2$, permitting exactly two electrons per orbital. This explains the structure of the periodic table — why the $s$-subshell holds 2 electrons, the $p$-subshell holds 6, and so on.

Particle physics: The Standard Model classifies all known elementary particles by spin: electrons, quarks, and neutrinos have spin $1/2$; photons, $W$ and $Z$ bosons, and gluons have spin $1$; the Higgs boson has spin $0$. The $\Delta^{++}$ baryon has spin $3/2$. All of these are representations of the angular momentum algebra we have developed.

What Cannot Be Determined Algebraically

The commutation relations tell us the possible values of $j$ and $m$, and how states transform under rotations. They do not tell us:

• Which value of $j$ a given particle carries (this is determined by experiment, or by deeper theories like the Standard Model)
• The magnetic moment — the relationship between angular momentum and magnetic moment requires additional physical input ($g$-factors)
• The dynamics — how angular momentum evolves in time depends on the Hamiltonian, not just the algebra
• The number of particles with each spin — the algebra tells us the representation theory, but not the particle content of the universe

12.10 Summary and Project Checkpoint

Summary

We have derived the complete algebraic theory of angular momentum. The logical chain is:

1. Starting point: Three commutation relations $[\hat{J}_i, \hat{J}_j] = i\hbar\epsilon_{ijk}\hat{J}_k$
2. Casimir operator: $\hat{J}^2$ commutes with all $\hat{J}_i$
3. Simultaneous eigenstates: $|j, m\rangle$ with $\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$
4. Ladder operators: $\hat{J}_\pm$ raise/lower $m$ by 1 within a $j$-multiplet
5. Spectrum: $j = 0, 1/2, 1, 3/2, \ldots$ and $m = -j, -j+1, \ldots, j$
6. Matrix representations: Explicit $(2j+1) \times (2j+1)$ matrices for any $j$
7. Rotation matrices: $D^{(j)}_{m'm}(\alpha, \beta, \gamma)$ from the exponential of angular momentum operators
8. Orbital angular momentum is the integer-$j$ special case; spin provides the half-integer cases

⚖️ Interpretation: The derivation in this chapter is one of the strongest arguments for the view that mathematical structure has physical content. We did not put half-integer angular momentum into the theory — we derived it from three algebraic relations that encode the geometry of rotations. Either the mathematics "knows" about physical reality in a deep way, or physical reality is constrained to follow the mathematics. Either way, the relationship between formalism and physics in quantum mechanics is far more intimate than in classical mechanics.

Connections to What Follows

The Generality of the Algebraic Method

Let us pause to appreciate the remarkable generality of what we have done. The same algebraic machinery applies to:

• Orbital angular momentum of electrons in atoms ($l = 0, 1, 2, \ldots$)
• Spin angular momentum of elementary particles ($s = 0, 1/2, 1, 3/2, \ldots$)
• Total angular momentum $\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}$ of an electron in an atom
• Nuclear spin $\hat{\mathbf{I}}$ of atomic nuclei (which can be integer or half-integer depending on the nucleus)
• Isospin in nuclear and particle physics (a mathematical quantity that obeys the same algebra, even though it has nothing to do with spatial rotations)
• Color charge in quantum chromodynamics (the gauge group $SU(3)$ generalizes the $SU(2)$ algebra we have studied)

In every case, the eigenvalue spectrum is $j(j+1)\hbar^2$ for $\hat{J}^2$ and $m\hbar$ for $\hat{J}_z$, with $j = 0, 1/2, 1, 3/2, \ldots$ and $m = -j, \ldots, j$. The ladder operators work identically, the matrix representations are constructed by the same algorithm, and the rotation matrices are computed by the same exponential formula. This universality is why angular momentum algebra is the most widely used mathematical tool in all of quantum physics.

Progressive Project Checkpoint: Angular Momentum Module (Toolkit v1.2)

Add the following to your Quantum Simulation Toolkit:

Module: angular_momentum.py

Functions to implement:

1. j_matrices(j) — Return $\hat{J}_x, \hat{J}_y, \hat{J}_z, \hat{J}_+, \hat{J}_-, \hat{J}^2$ as $(2j+1) \times (2j+1)$ NumPy arrays for any half-integer or integer $j$.

2. rotation_matrix(j, alpha, beta, gamma) — Return the $(2j+1) \times (2j+1)$ Wigner $D$-matrix for Euler angles $(\alpha, \beta, \gamma)$.

3. wigner_d(j, beta) — Return the reduced rotation matrix $d^{(j)}(\beta)$.

4. check_commutation(j) — Verify $[\hat{J}_x, \hat{J}_y] = i\hbar\hat{J}_z$ (and cyclic permutations) numerically for a given $j$. Return True if all relations hold to machine precision.

5. eigenvalue_test(j) — Verify $\hat{J}^2 = j(j+1)\hbar^2 \hat{I}$ and that $\hat{J}_z$ eigenvalues are $m\hbar$ for $m = -j, \ldots, j$.

Tests:

• All commutation relations hold for $j = 1/2, 1, 3/2, 2, 5/2$
• Rotation matrices are unitary
• $d^{(1/2)}(2\pi) = -\hat{I}_2$ and $d^{(1)}(2\pi) = +\hat{I}_3$
• Matrix elements match the analytical formulas for $j = 1/2$ and $j = 1$

See code/project-checkpoint.py for the complete implementation.

Key Formulas at a Glance

End-of-Chapter Reflection

Before moving on, take a moment to appreciate what we have accomplished. Starting from three lines of algebra — the commutation relations — we derived:

• The complete eigenvalue spectrum (including the existence of half-integer angular momentum)
• The action of ladder operators with exact normalization constants
• Explicit matrix representations for any $j$
• The rotation matrices that describe how quantum states transform under rotations

No other branch of physics yields so much from so little. This is the power of algebraic methods in quantum mechanics, and it is a pattern we will see again: in the addition of angular momenta (Chapter 14), in the classification of elementary particles (Chapter 32), and in the structure of quantum field theories (Chapter 35).

The theme of this chapter — that mathematics and physics are inseparable — will echo through every subsequent chapter. In quantum mechanics, the formalism is not a tool for computing predictions about a separately existing physical world. The formalism is the physics. Learn the algebra, and you learn the universe.

🔵 Historical Note: When Pauli used the algebraic method to derive the hydrogen spectrum in 1926 (before Schrödinger's wave equation!), he demonstrated that the abstract algebraic approach was not merely an alternative — it was, in some cases, the superior method. Heisenberg, who founded matrix mechanics, was so impressed by the power of the algebraic approach that he remarked, "The beauty of the mathematical structure is the surest guide to truth." Whether or not you agree with the metaphysics, the pragmatic lesson is clear: master the algebra.