Spherical Harmonics

Definition — angular Laplacian eigenfunctions

Write the 3D Laplacian in spherical coordinates (r, θ, φ) — radial coordinate r, polar angle θ ∈ [0, π], azimuth φ ∈ [0, 2π):

∇² = (1/r²) ∂/∂r (r² ∂/∂r) + (1/r²) ∇²_S²

∇²_S² = (1/sin θ) ∂/∂θ (sin θ ∂/∂θ) + (1/sin² θ) ∂²/∂φ²

∇²_S² is the angular Laplacian — the part that acts only on the directions, ignoring radius. Spherical harmonics are its eigenfunctions:

−∇²_S² Yₗᵐ(θ, φ) = ℓ(ℓ + 1) Yₗᵐ(θ, φ)

where ℓ is a non-negative integer. For each ℓ, the eigenvalue ℓ(ℓ+1) is the same for 2ℓ + 1 different eigenfunctions, labelled by m ∈ {−ℓ, −ℓ+1, …, ℓ−1, +ℓ}. ℓ counts the total number of angular wiggles; m counts how those wiggles split between latitude (θ-direction) and longitude (φ-direction).

Explicit formula

Separating variables Yₗᵐ(θ, φ) = Θ(θ) Φ(φ) reduces the eigenvalue equation to two ODEs. The φ equation gives Φ(φ) = e^{i m φ} (periodicity forces m to be an integer). The θ equation is the associated Legendre equation, whose regular solutions are the associated Legendre functions Pₗᵐ(cos θ). Putting it together:

Yₗᵐ(θ, φ) = √[((2ℓ+1)/(4π)) · ((ℓ−m)!/(ℓ+m)!)] · Pₗᵐ(cos θ) · e^{i m φ}

The first few:

Y₀⁰     = 1 / (2 √π)
Y₁⁰     = √(3/(4π)) cos θ
Y₁^{±1} = ∓ √(3/(8π)) sin θ e^{±i φ}
Y₂⁰     = √(5/(16π)) (3 cos²θ − 1)
Y₂^{±1} = ∓ √(15/(8π)) sin θ cos θ e^{±i φ}
Y₂^{±2} =   √(15/(32π)) sin²θ e^{±2 i φ}

Orthonormality and completeness

Spherical harmonics form a complete orthonormal basis for L²(S²):

∫_{S²} Yₗᵐ*(θ, φ) Yₗ′ᵐ′(θ, φ) dΩ = δ_{ℓℓ′} δ_{mm′}

f(θ, φ) = ∑_{ℓ=0}^∞ ∑_{m=−ℓ}^{+ℓ} a_{ℓm} Yₗᵐ(θ, φ)

a_{ℓm} = ∫_{S²} Yₗᵐ*(θ, φ) f(θ, φ) dΩ

Any square-integrable function on the sphere expands uniquely as a sum of Yₗᵐ — Fourier series, but for the sphere instead of the circle. The coefficients a_{ℓm} are the "spherical harmonic transform" of f.

Why 2ℓ + 1 functions per ℓ — group theory

The sphere has rotational symmetry SO(3). The Yₗᵐ for fixed ℓ span an irreducible representation of SO(3) of dimension 2ℓ + 1. Rotating the sphere mixes the m's among themselves but doesn't change ℓ — the eigenvalue ℓ(ℓ+1) is rotation-invariant. The number 2ℓ + 1 is forced by the irreducible-representation theory of SO(3).

This is why m loses physical meaning under rotation: it labels how the function looks in a particular coordinate system, not an invariant feature. In contrast, the "total power" Σ_m |a_{ℓm}|² is rotation-invariant — the basis for the CMB power spectrum C_ℓ.

Hydrogen orbitals

The time-independent Schrödinger equation for hydrogen, in spherical coordinates, separates:

ψ_{nℓm}(r, θ, φ) = R_{nℓ}(r) Yₗᵐ(θ, φ)

The radial wavefunction R_{nℓ}(r) involves Laguerre polynomials and depends on principal quantum number n and ℓ. The angular wavefunction is exactly a spherical harmonic. The chemistry shorthand:

• ℓ = 0 (s-orbital): 1 function — Y₀⁰ — spherically symmetric.
• ℓ = 1 (p-orbital): 3 functions — p_x, p_y, p_z (real combinations of Y₁ᵐ) — three "dumbbells".
• ℓ = 2 (d-orbital): 5 functions — d_{xy}, d_{xz}, d_{yz}, d_{x²−y²}, d_{z²} — cloverleaves and donuts.
• ℓ = 3 (f-orbital): 7 functions — more elaborate shapes.

The textbook orbital lobes you see are real linear combinations of Yₗᵐ (e.g., p_x = (Y₁¹ + Y₁⁻¹)/√2 up to sign). The degeneracy of each principal level n is n² states (from ℓ = 0 to n−1, each ℓ contributing 2ℓ + 1), giving 1, 4, 9, 16 … states in the n = 1, 2, 3, 4 shells — exactly the periodic-table block widths.

Cosmic Microwave Background

The CMB is a near-uniform 2.725 K blackbody filling the sky, but with tiny temperature anisotropies — ΔT/T ≈ 10⁻⁵. Decompose the temperature pattern T(θ, φ) − T̄ into spherical harmonics:

ΔT(θ, φ) = ∑_{ℓ,m} a_{ℓm} Yₗᵐ(θ, φ)

C_ℓ = (1/(2ℓ+1)) ∑_{m=−ℓ}^{+ℓ} |a_{ℓm}|²

C_ℓ is the angular power spectrum — the variance contributed by multipole ℓ. ℓ ≈ 180° / (angular scale), so ℓ = 2 is whole-sky, ℓ = 200 is ~1° (the first acoustic peak), ℓ = 2500 is ~0.07° (Silk damping scale).

Planck (2013–2018) measured C_ℓ from ℓ = 2 to ℓ ≈ 2500 with cosmic-variance precision over most of that range. The position, height, and ratios of the acoustic peaks determine the cosmological parameters: Ω_b ≈ 0.049, Ω_c ≈ 0.265, Ω_Λ ≈ 0.685, H₀ ≈ 67.4 km/s/Mpc, n_s ≈ 0.965. All from a spherical-harmonic decomposition.

Planetary gravity multipoles

Outside an isolated mass distribution, Laplace's equation ∇²V = 0 has separable solutions r^ℓ Yₗᵐ (interior) and r^{−ℓ−1} Yₗᵐ (exterior). The exterior gravity potential expands as

V(r, θ, φ) = − (GM/r) [ 1 + ∑_{ℓ≥1, m} (R/r)^ℓ (C_{ℓm} cos m φ + S_{ℓm} sin m φ) Pₗᵐ(cos θ) ]

The C_{ℓm}, S_{ℓm} (Stokes coefficients) describe how the body deviates from a perfect point mass at scale R^ℓ / r^ℓ. Earth's most important coefficient is J₂ = −C₂₀ ≈ 1.08 × 10⁻³ — the equatorial bulge from rotation. Satellites GRACE and GOCE measured Earth's gravity field to ℓ ≈ 360 — a global 0.5° gravity map, revealing ice loss, sea-level rise, and water-table changes from orbit.

The addition theorem

One of the most useful identities:

∑_{m=−ℓ}^{+ℓ} Yₗᵐ*(θ₁, φ₁) Yₗᵐ(θ₂, φ₂) = ((2ℓ + 1)/(4π)) Pₗ(cos γ)

where γ is the angle between the two directions and Pₗ is a Legendre polynomial. Two consequences:

1. The CMB angular power spectrum C_ℓ is rotation-invariant — the m-sum on the left has no preferred axis, only γ.
2. The Coulomb kernel expands as 1/|x − x′| = Σ_ℓ (r<^ℓ / r>^{ℓ+1}) Pₗ(cos γ) — the multipole expansion at the heart of classical electrostatics.

Spherical harmonics vs related bases

	Spherical harmonics Yₗᵐ	Fourier series	Hermite functions	Bessel functions
Domain	2-sphere S²	Circle / [0, 2π)	Real line ℝ	Disk / cylinder
Operator	−∇²_S²	−d²/dθ²	−d²/dx² + x²	Radial part of ∇²
Eigenvalues	ℓ(ℓ+1)	n²	2n + 1	Zeros of J_ν
Symmetry group	SO(3)	SO(2) ≃ U(1)	Heisenberg group	SO(2) on θ
Number per "level"	2ℓ + 1	1	1	1
Used in	Atoms, CMB, geodesy	Periodic signals	Quantum harmonic oscillator	Drum vibration, fiber modes
Coordinate	(θ, φ)	θ	x	r

Angular momentum eigenstates

In quantum mechanics, the orbital angular-momentum operators are L_x, L_y, L_z (vector cross-product of position and momentum). The total angular momentum squared L² = L_x² + L_y² + L_z² is proportional to −ℏ² ∇²_S². Spherical harmonics are simultaneous eigenfunctions:

L² Yₗᵐ = ℏ² ℓ(ℓ + 1) Yₗᵐ
L_z Yₗᵐ = ℏ m Yₗᵐ

ℓ is total angular-momentum quantum number; m is the z-projection. Selection rules for atomic spectral transitions (Δℓ = ±1, Δm = 0, ±1) come from matrix elements of the position operator between Yₗᵐ states — specifically, the integral of three spherical harmonics, expressible via Clebsch–Gordan or 3j symbols.

Where Yₗᵐ appears

• Atomic and molecular physics. Orbitals are R_{nℓ}(r) Yₗᵐ(θ, φ). All bonding (sp³, sp², sp hybridization) is real linear-combination arithmetic in Y-space.
• Cosmology. CMB analysis (Planck, WMAP, COBE) expands the sky map in Yₗᵐ. The power spectrum encodes the universe's content, geometry, and history.
• Geodesy. Earth's gravity, magnetic field, and topography all expand in spherical harmonics. EGM2008 (Earth Gravitational Model) goes to ℓ = 2190 — multi-cm resolution from satellite data.
• Seismology. Free oscillations of the Earth are modes labelled by (n, ℓ, m). After major earthquakes, the planet rings in spherical-harmonic modes for weeks.
• Computer graphics. BRDFs (bidirectional reflectance distributions) on a hemisphere expand in Yₗᵐ; "precomputed radiance transfer" projects environment lighting onto a low-ℓ spherical-harmonic basis.
• Quantum chemistry. Multi-electron wavefunctions are antisymmetric sums of Slater determinants over orbitals; the angular parts always involve Yₗᵐ.
• Crystal-field theory. Splitting of d-orbital energies in a crystal lattice is a tabulated decomposition of Y₂ᵐ under the lattice's symmetry group.

Common mistakes

• Conflating physics and mathematics conventions. Physicists usually take θ as polar (from z) and φ as azimuth; mathematicians often swap these. The Condon–Shortley phase factor (−1)^m is included in some conventions, not others. Signs of Yₗ^{±m} differ between Griffiths, Jackson, Sakurai, and SciPy.
• Forgetting the (sin θ) factor in dΩ. The integration measure on S² is dΩ = sin θ dθ dφ, not dθ dφ. Skip the sin θ and your orthogonality integrals come out wrong.
• Confusing m with the magnetic quantum number outside QM. In CMB analysis, m is just an azimuthal index, not anything magnetic. The notation is shared but the meaning depends on context.
• Mistaking real and complex Yₗᵐ. Complex Yₗᵐ are eigenfunctions of L_z (used in QM). The real combinations Y_x, Y_y, Y_z (used in chemistry orbital plots) are real but no longer L_z-eigenstates.
• Trying to interpret m geometrically as a "z-component count" outside QM. ℓ counts nodal-line total; m counts how many nodes are longitudinal vs latitudinal. Without QM context, m has no separate physical meaning.