Skip to content

Spherical harmonics are the angular wavefunctions that diagonalize L2L^2 and LzL_z. They are the orbital angular momentum basis functions on the sphere.

With the Condon–Shortley phase convention, they are written

Yℓm(θ,ϕ),ℓ=0,1,2,…,m=−ℓ,…,ℓ.Y_\ell^m(\theta,\phi), \qquad \ell=0,1,2,\ldots, \qquad m=-\ell,\ldots,\ell.

Spherical harmonics satisfy

L2Yℓm=ℏ2ℓ(ℓ+1)Yℓm,L^2Y_\ell^m = \hbar^2\ell(\ell+1)Y_\ell^m,

and

LzYℓm=ℏmYℓm.L_zY_\ell^m = \hbar mY_\ell^m.

Thus the labels ℓ\ell and mm are angular momentum quantum numbers:

j→ℓfor orbital angular momentum.j\to\ell \quad \text{for orbital angular momentum}.

Orbital ℓ\ell is integer because ordinary scalar wavefunctions on the sphere must be single-valued.

The inner product on the unit sphere uses the measure sin⁡θ dθ dϕ\sin\theta\,d\theta\,d\phi:

∫02π∫0πYℓm(θ,ϕ)∗Yℓ′m′(θ,ϕ)sin⁡θ dθ dϕ=δℓℓ′δmm′.\int_0^{2\pi}\int_0^\pi Y_\ell^m(\theta,\phi)^* Y_{\ell'}^{m'}(\theta,\phi) \sin\theta\,d\theta\,d\phi = \delta_{\ell\ell'}\delta_{mm'}.

Forgetting the measure is one of the most common errors.

Square-integrable angular functions can be expanded in spherical harmonics:

f(θ,ϕ)=∑ℓ=0∞∑m=−ℓℓcℓmYℓm(θ,ϕ).f(\theta,\phi) = \sum_{\ell=0}^{\infty} \sum_{m=-\ell}^{\ell} c_{\ell m}Y_\ell^m(\theta,\phi).

For a general square-integrable function, convergence here means convergence in the L2L^2 norm; pointwise convergence needs additional hypotheses. This is the spherical analogue of Fourier expansion. The quantum meaning is that the angular part of a wavefunction can be decomposed into orbital angular momentum eigenstates.

To remove an otherwise common sign ambiguity, define the phase-free associated Legendre function by

Pℓm(x)=(1−x2)m/2dmPℓ(x)dxm,0≤m≤ℓ.P_\ell^m(x)=(1-x^2)^{m/2}\frac{d^mP_\ell(x)}{dx^m},\qquad 0\leq m\leq\ell.

Then the normalized Condon–Shortley harmonics are

Yℓm(θ,ϕ)=(−1)m2ℓ+14π(ℓ−m)!(ℓ+m)!Pℓm(cos⁡θ)eimϕ,Y_\ell^m(\theta,\phi) =(-1)^m\sqrt{\frac{2\ell+1}{4\pi}\frac{(\ell-m)!}{(\ell+m)!}} P_\ell^m(\cos\theta)e^{im\phi},

with negative orders defined by Yℓ−m=(−1)m(Yℓm)∗Y_\ell^{-m}=(-1)^m(Y_\ell^m)^*. For example,

Y00=14π,Y10=34πcos⁡θ,Y_0^0=\frac{1}{\sqrt{4\pi}},\qquad Y_1^0=\sqrt{\frac{3}{4\pi}}\cos\theta, Y1±1=∓38πsin⁡θ e±iϕ.Y_1^{\pm1}=\mp\sqrt{\frac{3}{8\pi}}\sin\theta\,e^{\pm i\phi}.

This agrees with the associated-Legendre convention used here. Other references, including DLMF §14.30, absorb (−1)m(-1)^m into the polar function instead. Do not insert the phase twice: compare an explicit low-order function before combining tables of harmonics and Clebsch–Gordan coefficients.

Spherical harmonics have definite parity:

Yℓm(π−θ,ϕ+π)=(−1)ℓYℓm(θ,ϕ).Y_\ell^m(\pi-\theta,\phi+\pi) = (-1)^\ell Y_\ell^m(\theta,\phi).

This is why orbital angular momentum ℓ\ell controls parity for scalar central-potential wavefunctions.

The spherical harmonic addition theorem is

∑m=−ℓℓYℓm(Ω)∗Yℓm(Ω′)=2ℓ+14πPℓ(cos⁡γ),\sum_{m=-\ell}^{\ell} Y_\ell^m(\Omega)^* Y_\ell^m(\Omega') = \frac{2\ell+1}{4\pi} P_\ell(\cos\gamma),

where γ\gamma is the angle between the two directions Ω\Omega and Ω′\Omega'. This formula is central in partial waves, multipole expansions, and rotationally invariant kernels.

The polynomial PℓP_\ell is the Legendre Polynomial of degree ℓ\ell.

For a central potential, rotational symmetry lets wavefunctions separate as

ψ(r,θ,ϕ)=Rnℓ(r)Yℓm(θ,ϕ)\psi(r,\theta,\phi) = R_{n\ell}(r)Y_\ell^m(\theta,\phi)

when the radial problem has discrete labels. The spherical harmonic carries the angular momentum information; the radial function carries the dynamics of the specific potential.

In the Coulomb problem, the allowed angular multiplets inside each principal shell are summarized in Hydrogen Atom Angular Structure.

  • Forgetting the sin⁡θ\sin\theta measure.
  • Confusing ℓ\ell and mm.
  • Mixing phase conventions.
  • Treating visual orbital lobes as literal particle trajectories.
  • Assuming spherical harmonics are only hydrogen-atom objects; they are angular momentum basis functions.
  1. Verify the normalization of Y10Y_1^0. Then sum ∣Y1m∣2|Y_1^m|^2 over m=−1,0,1m=-1,0,1.
Solution

The normalization integral is

34π(2π)∫0πcos⁡2θsin⁡θ dθ=32∫−11u2 du=1.\frac{3}{4\pi}(2\pi)\int_0^\pi\cos^2\theta\sin\theta\,d\theta =\frac32\int_{-1}^{1}u^2\,du=1.

Using the explicit functions above gives

∑m=−11∣Y1m∣2=34π(cos⁡2θ+sin⁡2θ)=34π.\sum_{m=-1}^{1}|Y_1^m|^2 =\frac{3}{4\pi}(\cos^2\theta+\sin^2\theta) =\frac{3}{4\pi}.

An equal incoherent mixture of the three normalized states therefore has angular probability density 1/(4π)1/(4\pi). Individual members need not be isotropic; the complete equally weighted multiplet is.

  1. Use the parity rule to determine the parity of Y3mY_3^m.
Solution

The parity factor is (−1)ℓ(-1)^\ell. For ℓ=3\ell=3, this is −1-1, so every Y3mY_3^m has odd parity.

  • G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physicists, 7th ed., Academic Press, 2012.
  • A. R. Edmonds, Angular Momentum in Quantum Mechanics, Princeton University Press, 1957.
  • National Institute of Standards and Technology, Digital Library of Mathematical Functions, §14.30, Spherical and Spheroidal Harmonics, accessed 2026-10-13.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • D. A. Varshalovich, A. N. Moskalev, and V. K. Khersonskii, Quantum Theory of Angular Momentum, World Scientific, 1988.