Spherical harmonics and integer l

In spherical coordinates the Laplacian splits cleanly: ∇² = (1/r²)∂ᵣ(r²∂ᵣ) − L̂²/(ħ²r²), where L̂² is exactly the angular operator −ħ²[(1/sinθ)∂_θ(sinθ∂_θ) + (1/sin²θ)∂²_φ]. Put ψ = R(r)Y(θ, φ) into −(ħ²/2m)∇²ψ + V(r)ψ = Eψ and the angular part detaches as L̂²Y = ħ²l(l+1)Y, with l(l+1) the separation constant. Notice what is missing from that equation. There is no V, no E, no particle mass — nothing that could carry an energy. The angular problem is pure geometry on the unit sphere, so its solutions are identical for hydrogen, for the isotropic oscillator, for a finite spherical well and for the shell-model nucleus. All the angular half ever exports to the radial equation is one number, entering as the centrifugal term ħ²l(l+1)/2mr². Solve it once and you have solved it for every central force there will ever be.

Write Y(θ, φ) = Θ(θ)Φ(φ). Because L̂_z = −iħ∂/∂φ commutes with L̂², we may demand L̂_zY = mħY as well, and that fixes Φ = e^(imφ)/√(2π) at once. What remains is the associated Legendre equation, which in x = cosθ reads d/dx[(1 − x²)dP/dx] + [l(l+1) − m²/(1 − x²)]P = 0. Pull out the factor (1 − x²)^(|m|/2) that the m² term demands and expand the rest as a power series; the recursion ratio a_(k+2)/aₖ tends to 1 for large k, and a series behaving that way diverges at x = ±1 — the two poles — unless it terminates. It terminates only when l = |m| + k for some non-negative integer k, so given integer m, integer l follows and l ≥ |m| comes free. The first few results are Y₀⁰ = 1/√(4π) = 0.2821, Y₁⁰ = √(3/4π)cosθ = 0.4886 cosθ, and Y₂⁰ = √(5/16π)(3cos²θ − 1).

Everything above hangs on m being an integer, and the reason usually given is that ψ must return to itself when φ advances by 2π. That is a postulate dressed as a derivation, and a suspect one. Half-integer m gives Φ(φ + 2π) = −Φ(φ), an overall sign, and an overall phase changes no probability, no expectation value and no interference pattern within a single tower. Worse for the argument, nature uses exactly that behaviour: a spin-½ spinor picks up a minus sign under a 2π rotation and needs 4π to come home, and that sign has been seen directly in neutron interferometry. Nor can the half-integer candidate be dismissed as too singular to normalise. The l = ½ top rung works out to f = (1/π)√(sinθ)e^(iφ/2), and ∫|f|²dΩ = (1/π²)(2π)∫₀^π sin²θ dθ = (1/π²)(2π)(π/2) = 1. Finite, smooth away from the poles, and perfectly respectable-looking. If half-integer l is to be excluded, something else must do the excluding.

Use L̂_± = ±ħe^(±iφ)(∂_θ ± i cotθ ∂_φ). Acting on Θ(θ)e^(imφ), L̂₊ returns ħ(Θ′ − m cotθ Θ)e^(i(m+1)φ) and L̂₋ returns −ħ(Θ′ + m cotθ Θ)e^(i(m−1)φ). Demand L̂₊f = 0 at the top of a j = ½ tower and you get Θ′/Θ = ½cotθ, so Θ ∝ √(sinθ) — the normalisable state just quoted. Lower it once: L̂₋f = −ħ(1/π)cosθ(sinθ)^(−1/2)e^(−iφ/2), whose norm is exactly ħ, matching the ħ√[j(j+1) − m(m−1)] = ħ√(¾ + ¼) = ħ that the algebra demands. Still consistent. Now the tower must stop, because m = −½ is the bottom rung and L̂₋ has to annihilate it. It does not: the bracket comes out −(sinθ)¹⁄² − cos²θ(sinθ)^(−3/2), which vanishes nowhere on 0 < θ < π. And that leftover is not even on the sphere — its norm integrand is sin²θ + 2cos²θ + cos⁴θ/sin²θ, whose last term integrates like ∫dθ/θ² at θ → 0. The ladder neither closes nor stays normalisable. That is the exclusion.

Pₗ^|m|(cosθ) factorises as (sinθ)^|m| times a polynomial of degree l − |m| in cosθ, so |Yₗ^m|² has l − |m| nodal cones at fixed polar angles. For Y₂⁰ they sit where 3cos²θ = 1, at θ = 54.74° and 125.26° — the magic angle that NMR spinning exploits. The azimuthal nodes need care. |Yₗ^m|² is independent of φ, so the complex harmonic has none at all; take the real combinations that chemists draw, carrying cos(mφ) or sin(mφ), and |m| nodal planes containing the z axis appear. In the real form the total comes to l angular nodal surfaces. Parity is a one-line calculation: r → −r means θ → π − θ and φ → φ + π, giving (−1)^(l+|m|) from the Legendre factor and (−1)^m from the exponential, hence (−1)^l overall. Because the electric-dipole operator is odd, that parity forbids Δl = 0 and every even Δl, and the dipole's own l = 1 angular character then caps |Δl| at 1, leaving Δl = ±1. And every l carries 2l + 1 states, an exact degeneracy for any central potential; hydrogen's extra n²-fold degeneracy is a peculiarity of the 1/r potential, not a general rule.

The argument excludes half-integer l for angular momentum built from position and momentum — anything, that is, whose eigenstates must be functions on the sphere. It says nothing about angular momentum that has no such representation. Spin obeys identical commutators, so the ladder derivation of j(j+1)ħ² and 2j+1 projections applies to it unchanged, but there is no Θ(θ) for the l = ½ calculation to break, and s = ½ survives untouched. That is the honest statement of where the sphere's authority ends. Two further limits are worth naming. The angular equation returns no energy at all, so nothing here predicts a spectrum: the 2l + 1 degeneracy within a level is exact, and everything else waits on V(r) and the radial equation. And that degeneracy is a statement about a rotationally symmetric Hamiltonian — switch on a magnetic field and it lifts, which is the entire content of the Zeeman effect coming next.