One dimension was training. Real atoms are three-dimensional, and the force an electron feels
from a nucleus depends only on distance — a central potential. That symmetry splits the
Schrödinger equation into a radial problem and a universal angular one, solved by the same ladder
trick that cracked the harmonic oscillator — now climbing angular momentum instead of energy. A
cat circling a warm laptop conserves L; so does an electron circling a proton — but
only the electron's comes in steps of ℏ.
Learning Objectives
After this lesson you will be able to:
Separate variables in the 3D Schrödinger equation for a central potential V(r), identifying the separation constant ℓ(ℓ+1).
Derive the commutator algebra[L^x,L^y]=iℏL^z from the canonical commutators and recognize its classical shadow, the Poisson bracket{Lx,Ly}=Lz.
Construct the ladder operatorsL^± and derive the eigenvalues of L^2 and L^z purely algebraically.
Explain why orbital angular momentum requires integer ℓ, while the algebra alone also permits half-integer values.
Use the spherical harmonicsYℓm — eigenvalue equations, orthonormality, the 2ℓ+1degeneracy — and critique the vector model via ⟨Lx2+Ly2⟩>0.
Intuition
Classically a central force exerts no torque, so L=r×p is conserved
and orbits flatten into planes (P.1.3 Hamiltonian Mechanics).
Quantum mechanics inherits the structure: rotational symmetry peels the angular dependence off the
wavefunction once and for all — the angular factors are the same for hydrogen, a 3D oscillator,
any V(r). But there is a twist: classically Lx,Ly,Lz are three numbers and you may know
them all; quantum mechanically they are noncommuting operators — you may know the length
∣L∣ and one component, and the other two are then irreducibly indeterminate. The
entire spectrum follows from the commutators alone, by the ladder logic of
P.5.3 The Harmonic Oscillator.
Theory
Separation of variables in a central potential
The 3D time-independent Schrödinger equation is −2mℏ2∇2ψ+Vψ=Eψ.
For V=V(r) use spherical coordinates, with the Laplacian (a standard chain-rule exercise, quoted):
Try ψ(r,θ,φ)=R(r)Y(θ,φ); substitute and multiply by −ℏ22mr2RY1:
depends only on r{R1drd(r2drdR)−ℏ22mr2[V(r)−E]}+depends only on θ,φY1{sinθ1∂θ∂(sinθ∂θ∂Y)+sin2θ1∂φ2∂2Y}=0.
A function of r plus a function of (θ,φ) vanishes identically only if each is
constant; write the constant — with foresight — as ℓ(ℓ+1). The radial piece keeps V(r)
and is deferred to P.6.2; the angular piece,
sinθ1∂θ∂(sinθ∂θ∂Y)+sin2θ1∂φ2∂2Y=−ℓ(ℓ+1)Y,
is universal — and, as we now show, it is the eigenvalue problem of angular momentum.
The angular momentum operators and their algebra
Quantizing L=r×p gives L^=r^×p^,
with components L^x=y^p^z−z^p^y, L^y=z^p^x−x^p^z,
L^z=x^p^y−y^p^x (no ordering ambiguity: each product pairs a coordinate
with a different momentum component). Using [x^i,p^j]=iℏδij, expand
[L^x,L^y] into four brackets; [y^p^z,x^p^z] and
[z^p^y,z^p^x] contain only mutually commuting factors and vanish, leaving
and cyclically [L^y,L^z]=iℏL^x, [L^z,L^x]=iℏL^y. This
mirrors the classical Poisson bracket {Lx,Ly}=Lz computed in
P.1.3 Hamiltonian Mechanics via
Dirac's rule {⋅,⋅}→[⋅,⋅]/iℏ — the algebra of rotations was there all along.
No two components commute, so no common eigenbasis exists: Lx and Lz cannot both be sharp
(the uncertainty relation behind this is developed in
1.3.2 Expectation & Uncertainty).
But L^2=L^x2+L^y2+L^z2 commutes with every component: using
[A^2,B^]=A^[A^,B^]+[A^,B^]A^,
so L^2 and L^zdo share simultaneous eigenfunctions — the maximal sharp angular
data. In spherical coordinates (a rotation about z changes only φ, so
∂φ=x∂y−y∂x):
Let ∣λ,μ⟩ satisfy L^2∣λ,μ⟩=λ∣λ,μ⟩
and L^z∣λ,μ⟩=μ∣λ,μ⟩. Then
L^z(L^±∣λ,μ⟩)=(L^±L^z±ℏL^±)∣λ,μ⟩=(μ±ℏ)L^±∣λ,μ⟩
while L^2(L^±∣λ,μ⟩)=λL^±∣λ,μ⟩:
a ladder in Lz, steps of ℏ, fixed total length. It is bounded, since
λ−μ2=⟨L^x2+L^y2⟩=∥L^x∣λ,μ⟩∥2+∥L^y∣λ,μ⟩∥2≥0.
So there is a top rung, L^+∣λ,μmax⟩=0 with μmax≡ℏℓ.
Expanding L^∓L^±=L^x2+L^y2±i[L^x,L^y]=L^2−L^z2∓ℏL^z
and applying L^−L^+ to the top rung: 0=λ−ℏ2ℓ2−ℏ2ℓ, i.e.
λ=ℏ2ℓ(ℓ+1). Likewise a bottom rung ℏℓˉ gives
λ=ℏ2ℓˉ(ℓˉ−1), so ℓˉ=−ℓ (the root ℓ+1 would put the
bottom above the top). Descending from ℏℓ to −ℏℓ in steps of ℏ requires
2ℓ to be a non-negative integer:
The commutators alone — no wavefunctions, no boundary conditions — fixed the entire spectrum,
and they permit both integer and half-integerℓ.
The differential route: why orbital ℓ is an integer
For angular momentum realized as −iℏr×∇ on wavefunctions, more is true.
Solving L^zΦ=ℏmΦ gives Φ=eimφ; but φ and
φ+2π label the same point of space, and single-valuedness
eim(φ+2π)=eimφ forces m∈Z — hence integer ℓ for
orbital angular momentum. The half-integer solutions are not garbage: they describe angular
momentum that is not motion through space — no φ, no constraint. That loophole stays
open until the Stern–Gerlach experiment forces us through it in
P.6.3 Magnetic Moments, Stern–Gerlach & Spin.
Spherical harmonics
The normalized simultaneous eigenfunctions of L^2 and L^z are the spherical
harmonics, built from the associated Legendre functions Pℓm:
They are orthonormal on the sphere — ∫02π∫0πYℓ′m′∗Yℓmsinθdθdφ=δℓℓ′δmm′
— and complete. For each ℓ there are 2ℓ+1 values of m with the same L^2
eigenvalue: a degeneracy every central potential inherits, since energy cannot depend on the
orientation of L when the potential has no preferred direction.
The vector model and its limits
A common mnemonic draws L as an arrow of length ℏℓ(ℓ+1) on a cone
about z, with fixed projection ℏm and smeared azimuth. The cone's opening is real physics:
⟨Lx2+Ly2⟩=⟨L^2−L^z2⟩=ℏ2[ℓ(ℓ+1)−m2]>0for ℓ>0,
even at m=ℓ (where it equals ℏ2ℓ): L never points exactly along z,
tilting at best to cosθmin=ℓ/ℓ(ℓ+1)<1. But don't take the cone
literally — nothing secretly precesses at a hidden azimuth.
Caution.∣L∣=ℏℓ(ℓ+1)>ℏℓ≥ℏm. The Bohr-era phrase
"the angular momentum is ℓℏ" is shorthand and strictly false — the length exceeds the
largest projection, which is exactly why L can never align with an axis. And a state
with definite Lz has indeterminateLx and Ly: not unknown values — no values.
Worked Examples
Example 1 — Sizing up ℓ=2
∣L∣=ℏ2⋅3=6ℏ≈2.449ℏ, with five projections
Lz=mℏ, m∈{−2,…,2}. Even at m=2 the vector tilts off-axis by
cosθ=2/6=0.8165, i.e. θ≈35.3∘, with transverse spread
⟨Lx2+Ly2⟩=ℏ2(6−4)=2ℏ2. The largest projection 2ℏ never
reaches the length 2.449ℏ.
Example 2 — Y10 through the machinery
Y10=3/4πcosθ has no φ-dependence, so L^zY10=0: m=0. ✓
For L^2 (the φ term drops):
so L^2Y10=2ℏ2Y10=ℏ2ℓ(ℓ+1)Y10 with ℓ=1. ✓ Lowering with
L^−=−ℏe−iφ(∂θ−icotθ∂φ):
L^−Y10=ℏ3/4πsinθe−iφ=ℏ2Y1−1, with
coefficient exactly ℏℓ(ℓ+1)−m(m−1)=ℏ2 — the ladder in action.
Hands-on (Python)
Visualize ∣Yℓm∣2 for ℓ=0,1,2 as polar surface plots (radius = magnitude):
import numpy as np
import matplotlib.pyplot as plt
from scipy.special import sph_harm
# NOTE scipy's argument order/naming: sph_harm(m, l, azimuthal, polar).
polar = np.linspace(0, np.pi, 101) # theta in physics texts
azim = np.linspace(0, 2*np.pi, 201) # varphi in physics texts
TH, PH = np.meshgrid(polar, azim)
fig = plt.figure(figsize=(13, 8))
for l inrange(3):
for m inrange(-l, l + 1):
rad = np.abs(sph_harm(m, l, PH, TH))**2# radius = |Y_l^m|^2
X, Y, Z = (rad*np.sin(TH)*np.cos(PH), rad*np.sin(TH)*np.sin(PH), rad*np.cos(TH))
ax = fig.add_subplot(3, 5, 5*l + m + 3, projection='3d') # center each row
ax.plot_surface(X, Y, Z, cmap='viridis', rstride=2, cstride=2)
ax.set_title(f"$\\ell={l},\\ m={m}$"); ax.set_axis_off()
plt.tight_layout(); plt.show()
# Expect: l=0 a sphere; l=1 dumbbells (m=0 along z, |m|=1 donuts); l=2 the# donut/dumbbell family. |Y|^2 is always symmetric about the z-axis.
Verify orthonormality ∫Yℓ′m′∗YℓmdΩ=δℓℓ′δmm′ by 2D quadrature:
defoverlap(l1, m1, l2, m2, n_th=400, n_ph=800):
th, ph = np.linspace(0, np.pi, n_th), np.linspace(0, 2*np.pi, n_ph)
TH, PH = np.meshgrid(th, ph, indexing='ij')
f = np.conj(sph_harm(m1, l1, PH, TH))*sph_harm(m2, l2, PH, TH)*np.sin(TH)
return np.trapezoid(np.trapezoid(f, ph, axis=1), th) # np.trapz on older NumPy
labels = [(l, m) for l inrange(3) for m inrange(-l, l + 1)] # 9 states
G = np.array([[overlap(*a, *b) for b in labels] for a in labels])
print(np.max(np.abs(G - np.eye(9)))) # ~1e-6: Gram matrix = identity
Exercises
E1 (easy). For ℓ=3: how many m values? Compute ∣L∣ and the smallest possible
angle between L and the z-axis.
Solution
2ℓ+1=7 values, m=−3,…,3; ∣L∣=ℏ12=23ℏ≈3.464ℏ.
At m=3: cosθ=3/12=3/2, so θ=30∘ exactly.
E2 (easy). Derive [L^z,L^x]=iℏL^y (a) directly from the canonical
commutators, (b) by cyclic relabeling.
Solution
(a) [x^p^y−y^p^x,y^p^z−z^p^y] has nonvanishing pieces
[x^p^y,y^p^z]=x^[p^y,y^]p^z=−iℏx^p^z and
[y^p^x,z^p^y]=z^p^x[y^,p^y]=iℏz^p^x; sum
=iℏ(z^p^x−x^p^z)=iℏL^y. (b) x→y→z→x preserves
[x^i,p^j]=iℏδij and maps L^x→L^y→L^z→L^x,
carrying [L^x,L^y]=iℏL^z into the other two identities.
E3 (medium). Prove L^∓L^±=L^2−L^z2∓ℏL^z and deduce
L^±∣ℓ,m⟩=ℏℓ(ℓ+1)−m(m±1)∣ℓ,m±1⟩.
Solution
$(\hat L_x \mp i\hat L_y)(\hat L_x \pm i\hat L_y) = \hat L_x^2 + \hat L_y^2 \pm i[\hat L_x,\hat L_y]
= \hat L^2 - \hat L_z^2 \mp \hbar\hat L_z$. Then
$\lVert\hat L_\pm\lvert \ell,m\rangle\rVert^2 = \langle \ell,m\rvert\hat L_\mp\hat L_\pm\lvert \ell,m\rangle
= \hbar^2[\ell(\ell+1) - m(m\pm1)]$; the standard positive phase gives the stated coefficient,
which vanishes at m=±ℓ: the ladder terminates itself.
E4 (medium). By explicit integration, verify that Y11 is normalized and orthogonal to Y10.
Solution
∫∣Y11∣2dΩ=8π3⋅2π∫0πsin3θdθ=8π3⋅2π⋅34=1 ✓.
Orthogonality: ∫Y10∗Y11dΩ∝∫02πeiφdφ=0 — the
azimuthal integral kills any pair with m=m′. ✓
E5 (hard). In ∣ℓ,m⟩, show ⟨Lx⟩=⟨Ly⟩=0 and
⟨Lx2⟩=⟨Ly2⟩=2ℏ2[ℓ(ℓ+1)−m2], then verify
ΔLxΔLy≥2ℏ∣⟨Lz⟩∣.
Solution
L^x=21(L^++L^−) maps ∣ℓ,m⟩ to orthogonal states, so
⟨Lx⟩=0 (same for Ly). In L^x2=41(L^+2+L^−2+L^+L^−+L^−L^+)
the L^±2 terms shift m by ±2 and drop; E3 gives
⟨L^+L^−+L^−L^+⟩=2ℏ2[ℓ(ℓ+1)−m2], so
⟨Lx2⟩=⟨Ly2⟩=2ℏ2[ℓ(ℓ+1)−m2]. Then
ΔLxΔLy≥2ℏ2∣m∣⟺ℓ(ℓ+1)≥∣m∣(∣m∣+1), true since
ℓ≥∣m∣ — with equality at m=±ℓ: the end rungs are minimum-uncertainty states.
Checkpoint
In ψ=RY, what is the separation constant, and what operator equation does Y satisfy?
State [L^x,L^y] and its classical counterpart. What rule connects them?
Why do L^2 and L^z share eigenfunctions while L^x and L^z do not?
Which step of which argument restricts orbital ℓ to integers — and why doesn't it bind spin?
Why can L never point exactly along the z-axis for ℓ>0?
Answers
ℓ(ℓ+1); the angular equation is L^2Y=ℏ2ℓ(ℓ+1)Y.
[L^x,L^y]=iℏL^z, mirroring {Lx,Ly}=Lz via Dirac's rule {⋅,⋅}→[⋅,⋅]/iℏ.
[L^2,L^z]=0 guarantees a common eigenbasis; [L^x,L^z]=−iℏL^y=0 forbids one.
Single-valuedness of eimφ under φ→φ+2π in the differential realization forces m∈Z. Spin is not a function of spatial angles, so the constraint never applies — the algebra's half-integers survive.
∣L∣=ℏℓ(ℓ+1) strictly exceeds the largest projection ℏℓ; equivalently ⟨Lx2+Ly2⟩=ℏ2[ℓ(ℓ+1)−m2]>0.
Further Reading
[Gri] Griffiths & Schroeter, §4.1 — separation of variables and the angular equation; this lesson's conventions.
[Gri] Griffiths & Schroeter, §4.3 — the algebraic theory of angular momentum; our ladder derivation follows it closely.
[Sha] Shankar, Ch. 12 — angular momentum from rotational symmetry; the algebra as the group theory of rotations.
[ER] Eisberg & Resnick, Ch. 7 — quantum numbers with the experimental context of one-electron atoms.