Spherical Harmonics

Open any chemistry book and you meet them: the round s orbital, the dumbbell-shaped p orbitals, the four-lobed d orbitals. Those iconic shapes are not artists' impressions — they are the graphs of a single family of mathematical functions, the spherical harmonics Y_\ell^m(\theta,\varphi). They are what you get when you ask: what are the natural vibration patterns of a sphere?

A circle's natural patterns are \cos(m\varphi) and \sin(m\varphi) — the Fourier modes. A sphere's natural patterns are the spherical harmonics. They are the angular "shape functions" that fall out whenever the Laplacian \nabla^2 is separated in spherical coordinates, so they turn up everywhere spheres do: the electron clouds of atoms, the ripples of the cosmic microwave background, the lumps in the Earth's gravity field, the tone of a ringing bell. This page is about that one family.

The angular part of the Laplacian

Write \nabla^2 in spherical coordinates (r,\theta,\varphi) and it splits into a radial piece and an angular piece. The angular piece — call it \nabla^2_{\theta\varphi} — is an operator on the surface of the sphere, and the spherical harmonics are exactly its eigenfunctions:

\nabla^2_{\theta\varphi}\,Y_\ell^m(\theta,\varphi) = -\ell(\ell+1)\,Y_\ell^m(\theta,\varphi).

That eigenvalue -\ell(\ell+1) is the whole reason \ell is quantised to whole numbers, and in quantum mechanics it is precisely the squared angular momentum: \hat{L}^2\,Y_\ell^m = \hbar^2\,\ell(\ell+1)\,Y_\ell^m. Separating the two angles further, Y_\ell^m = \Theta(\theta)\,\Phi(\varphi), gives a beautifully structured formula:

The \varphi part is just a Fourier mode e^{im\varphi} — the circle's harmonics wrapped around the equator. The \theta part is the associated Legendre function. In fact, when m=0 the harmonic loses its \varphi dependence and reduces to an ordinary Legendre polynomial: Y_\ell^0 \propto P_\ell(\cos\theta). The special functions all connect.

The lowest harmonics — and the orbital shapes

The first few, written out, are strikingly simple:

Y_0^0 = \frac{1}{\sqrt{4\pi}}, \qquad Y_1^0 = \sqrt{\frac{3}{4\pi}}\,\cos\theta, \qquad Y_2^0 = \sqrt{\frac{5}{16\pi}}\,\big(3\cos^2\theta - 1\big).

Y_0^0 is a constant — perfectly round, no angular variation. That is the s orbital. Y_1^0 \propto \cos\theta is positive on top and negative below, pinched to zero at the equator — the dumbbell of a p orbital. Y_2^0 \propto 3\cos^2\theta-1 adds a waist, giving the peanut-and-ring of a d orbital. The figure plots the polar cross-sections r = |Y_\ell^0(\theta)| for \ell = 0, 1, 2 — the distance from the origin is the size of the harmonic in that direction. Reveal them one at a time.

Count the places where each curve pinches back to the origin: that number — the angular nodes — is exactly \ell. More \ell means a wigglier, more structured pattern on the sphere, just as a higher harmonic means a wigglier standing wave on a string.

Orthonormal on the sphere

Like sines and cosines on an interval and Legendre polynomials on [-1,1], the spherical harmonics are an orthonormal set — but now the integration runs over the whole surface of the sphere, with the area element d\Omega = \sin\theta\,d\theta\,d\varphi:

\int_0^{2\pi}\!\!\int_0^{\pi} Y_\ell^m(\theta,\varphi)\,\overline{Y_{\ell'}^{m'}(\theta,\varphi)}\;\sin\theta\,d\theta\,d\varphi = \delta_{\ell\ell'}\,\delta_{mm'}.

(The bar is a complex conjugate, because the e^{im\varphi} factor makes the harmonics complex.) This is the master property. It means any function on the sphere — a temperature map of the Earth, the brightness of the microwave sky, an electron's angular wavefunction — can be expanded as a sum f(\theta,\varphi) = \sum_{\ell,m} c_{\ell m}\,Y_\ell^m(\theta,\varphi), and each coefficient is extracted by a single integral, exactly as with Fourier and Legendre series. Cosmologists do this to the entire universe: the famous "power spectrum" of the cosmic microwave background is just the list of |c_{\ell m}|^2.

Worked examples

Example 1 — counting harmonics. How many spherical harmonics share the degree \ell = 3? The order runs m = -3,-2,-1,0,1,2,3, so there are 2\ell+1 = 7 of them. (These are the seven f orbitals.)

Example 2 — the eigenvalue. What squared angular momentum does Y_2^m carry? Use \hat{L}^2 Y = \hbar^2\ell(\ell+1)Y with \ell=2: \hbar^2\cdot 2\cdot 3 = 6\hbar^2, so |\mathbf{L}| = \sqrt{6}\,\hbar, independent of m.

Example 3 — the link to Legendre. Write Y_2^0 in terms of a Legendre polynomial. Since P_2(x)=\tfrac12(3x^2-1) and x=\cos\theta, we have Y_2^0 = \sqrt{\tfrac{5}{4\pi}}\,P_2(\cos\theta) — the m=0 harmonic is a scaled Legendre polynomial.

Example 4 — angular nodes. How many angular nodes does Y_4^0 have? A degree-\ell harmonic has \ell angular nodes, so Y_4^0 has 4 — four latitude circles where it vanishes (a g orbital).

It is 2\ell+1, not 2\ell. The order m runs over every integer from -\ell to +\ell including zero: -\ell, \dots, -1, 0, 1, \dots, +\ell. Forgetting the m=0 term (or the "+1") is the classic miscount — it is why there are 3 p orbitals, 5 d orbitals and 7 f orbitals, not 2, 4 and 6.

Two more traps. First, spherical harmonics are functions of angle only (\theta,\varphi) — they are the angular part of a solution; the radial dependence is a separate factor. Second, because of the e^{im\varphi} the harmonics are complex for m\ne 0. The real, lobed pictures in chemistry books (p_x, p_y) are combinations Y_\ell^{m} \pm Y_\ell^{-m} — the same information, rearranged into real form.

Because the hydrogen atom's Schrödinger equation separates just like Laplace's equation. Its solution factorises as \psi(r,\theta,\varphi) = R_{n\ell}(r)\,Y_\ell^m(\theta,\varphi) — a radial part times a spherical harmonic. The probability of finding the electron in a given direction is |\psi|^2, whose angular shape is |Y_\ell^m|^2. So the dumbbells and cloverleaves you were shown in chemistry are literally plots of |Y_\ell^m(\theta,\varphi)|^2.

The quantum numbers you memorised are the harmonic's labels in disguise: \ell is the orbital angular-momentum quantum number (s, p, d, f = \ell=0,1,2,3) and m is the magnetic quantum number that splits under a field. The whole architecture of the periodic table — why shells hold 2, 6, 10, 14 electrons — is the counting 2(2\ell+1) of spherical harmonics, doubled for spin.