Particle on a Sphere and Spherical Harmonics

Introduction

A particle on a ring rotates in one plane. Let the particle roam over the whole surface of a sphere and rotation becomes three-dimensional. This is the model for a rotating diatomic molecule and, crucially, it supplies the angular part of every atomic orbital. The eigenfunctions, called spherical harmonics, are the mathematical origin of the shapes of s, p, d and f orbitals, so understanding them turns orbital pictures from memorised drawings into consequences of theory.

Core explanation

The model. A particle of mass m moves freely on a sphere of fixed radius r. The potential energy is zero and r never changes, so only angular motion contributes. The Hamiltonian contains only the angular part of the Laplacian, and since that part is L̂²/ħ² up to a sign, the Hamiltonian becomes simply

Ĥ = L̂²/(2I), with I = mr²

Energies. The wavefunctions must be eigenfunctions of L̂², which we know has eigenvalues l(l + 1)ħ². Therefore

E = l(l + 1)ħ²/(2I), l = 0, 1, 2, …

The levels lie at 0, 2, 6, 12, 20 … in units of ħ²/(2I). As on a ring, the lowest energy is zero.

Two boundary conditions, two quantum numbers. Motion around the z-axis (in φ) must satisfy the cyclic condition, giving the integer mₗ as for the ring. Motion from pole to pole (in θ) must give a wavefunction that stays finite at θ = 0 and θ = π. This second condition limits mₗ so that mₗ ≤ l. A three-dimensional rotation thus needs two quantum numbers, l and mₗ.

Spherical harmonics. The eigenfunctions are written Y(l, mₗ)(θ, φ) and factorise as Θ(θ) × e^(imₗφ). The first few, normalised over the surface of a sphere, are

Y(0, 0) = (1/4π)^½ Y(1, 0) = (3/4π)^½ cos θ Y(1, ±1) = ∓(3/8π)^½ sin θ e^(±iφ) Y(2, 0) = (5/16π)^½ (3cos²θ − 1)

Y(0, 0) is the same in every direction: this is the angular part of every s orbital. Y(1, 0) is proportional to cos θ, which equals z/r, so it has the shape of a p z orbital, with a nodal plane at θ = 90°. Y(2, 0) contains 3cos²θ − 1, the characteristic shape of the d z² orbital with two nodal cones.

Nodes. Each Y(l, mₗ) has exactly l angular nodes (planes or cones). This is why s orbitals have none, p orbitals one and d orbitals two.

Degeneracy. The energy depends only on l, not on mₗ, so each level is (2l + 1)-fold degenerate. Physically, the energy of rotation cannot depend on how the rotation axis is oriented in empty space; spherical symmetry guarantees the degeneracy.

Space quantisation. Because L z = mₗħ can take only 2l + 1 values, the angular momentum vector can adopt only 2l + 1 orientations relative to a chosen axis. This orientation restriction was demonstrated directly in the Stern–Gerlach experiment, where a beam of atoms split into discrete beams in an inhomogeneous magnetic field.

Formulae

E = l(l + 1)ħ²/(2I); degeneracy g = 2l + 1; L² = l(l + 1)ħ²; L z = mₗħ; mₗ ≤ l; number of angular nodes = l. For molecular rotation the letter J replaces l: E J = J(J + 1)ħ²/(2I).

Step-by-step reasoning

To analyse a state of a particle on a sphere:

1. Identify l; the energy is l(l + 1)ħ²/(2I). 2. List the 2l + 1 allowed values of mₗ. 3. Recognise the matching spherical harmonic and count its l angular nodes. 4. Relate the shape to the corresponding atomic orbital: l = 0 is s-like, 1 p-like, 2 d-like.

Visual explanation

Picture a globe. Y(0, 0) paints it a single uniform colour. Y(1, 0) paints the northern hemisphere positive and the southern negative, the equator being a node. Y(2, 0) paints both polar caps positive and a band around the equator negative, separated by two cones at about 54.7° from the poles.

Real-world analogy

The vibrations of a soap bubble or a ringing bell follow similar patterns: a whole-body pulsation, a sloshing from one side to the other, and more complex patterns with lines of stillness between regions moving in opposite directions. The lines of stillness are the counterparts of angular nodes.

Real-world example

Microwave spectroscopy of molecules such as carbon monoxide measures transitions between rotational levels E J = J(J + 1)ħ²/(2I). The lines are nearly equally spaced, about 3.84 cm⁻¹ apart for CO, and from that spacing the bond length of 113 pm is obtained. Radio astronomers use the same lines to map cold molecular clouds.

Why?

Why do angular shapes of orbitals appear in the hydrogen atom at all? Because the Coulomb potential depends only on r, the angular motion of the electron is exactly that of a particle on a sphere. The atom must therefore have the same spherical harmonics as its angular wavefunctions, whatever its radial behaviour.

Common misconception

"The p x, p y and p z orbitals correspond to mₗ = −1, 0 and +1." Only p z corresponds directly to mₗ = 0. The p x and p y orbitals are real combinations of the complex mₗ = +1 and −1 functions and do not have definite values of mₗ.

Worked example

Question: For a particle on a sphere, find the energy of the l = 2 level in units of ħ²/(2I), its degeneracy, and the energy gap to the l = 3 level.

Reasoning: E₂ = 2 × 3 = 6 units. Degeneracy = 2 × 2 + 1 = 5. E₃ = 3 × 4 = 12 units, so the gap is 12 − 6 = 6 units. In general the gap from l to l + 1 is 2(l + 1) units.

Answer: 6ħ²/(2I), five-fold degenerate, and a gap of 6ħ²/(2I), which equals 3ħ²/I.

Quick check

1. How many angular nodes does the spherical harmonic Y(3, 0) have, and which type of orbital shares it? Answer: It has three angular nodes, the number always equal to l, and it gives the angular shape of an f orbital.

Exam focus

Know E = l(l + 1)ħ²/(2I), the degeneracy 2l + 1 and the limit mₗ ≤ l. Be able to write Y(0, 0) and Y(1, 0) and link them to s and p z orbitals. Examiners also test the difference between one quantum number on a ring and two on a sphere.

Advanced insight

The real combinations of spherical harmonics used in chemistry, such as p x ∝ sin θ cos φ and d xy ∝ sin²θ sin 2φ, are exactly proportional to x/r, xy/r² and similar Cartesian expressions. This is why orbital labels mirror Cartesian functions. Any angular function on a sphere can be expanded as a sum of spherical harmonics, much as any periodic signal can be built from sines and cosines.

Summary

A particle on a sphere has Ĥ = L̂²/(2I), so its energies are l(l + 1)ħ²/(2I) with l = 0, 1, 2, … Cyclic motion in φ gives the integer mₗ; finiteness at the poles limits mₗ ≤ l. The eigenfunctions are spherical harmonics Y(l, mₗ), each with l angular nodes, and they provide the angular shapes of atomic orbitals. Each level is (2l + 1)-fold degenerate and angular momentum is space-quantised.

Practice questions

1. What is the energy of the l = 1 level of a particle on a sphere? Answer: E = 1 × 2 × ħ²/(2I) = ħ²/I. 2. Which condition limits mₗ to values between −l and +l? Answer: The requirement that the θ part of the wavefunction remains finite at the poles, θ = 0 and θ = π. 3. Show that Y(1, 0) has a nodal plane and identify it. Answer: Y(1, 0) is proportional to cos θ, which is zero at θ = 90°, so the xy-plane is a nodal plane. 4. Why is each level of the particle on a sphere degenerate? Answer: The energy depends only on the magnitude of the angular momentum, not its orientation, and in a spherically symmetric system all 2l + 1 orientations are equivalent.