The hydrogen atom is the only multi-particle system with an exact analytical solution to the Schrödinger equation. The wavefunctions ψ_{nlm} are products of radial functions R_{nl}(r) and spherical harmonics Y_l^m(θ,φ), each labeled by three quantum numbers: principal (n), angular momentum (l), and magnetic (m). Energy levels depend only on n and go as E_n = −13.6/n² eV. The radial probability distribution P(r) = r²|R_{nl}|² reveals where electrons are most likely to be found, directly explaining orbital shapes and the concept of shells.
Plot radial probability distributions for s, p, and d orbitals and count nodes — n−l−1 radial nodes and l angular nodes. Connect each quantum number to a physical property: n → energy and size, l → shape, m → orientation.
The hydrogen atom holds a unique place in quantum chemistry: it is the only atom for which the Schrödinger equation can be solved exactly, producing closed-form wavefunctions. Everything you know about atomic orbitals — their shapes, their quantum numbers, their energies — derives directly from this solution.
The wavefunction ψ_{nlm}(r,θ,φ) factors into two independent pieces: a radial part R_{nl}(r) that depends only on distance from the nucleus, and an angular part Y_l^m(θ,φ) — a spherical harmonic — that describes the directional shape. The three quantum numbers encode distinct physical information: n determines the energy (E_n = −13.6/n² eV) and the overall size of the orbital; l determines the shape (l = 0 is spherical, l = 1 has a dumbbell shape, l = 2 is cloverleaf); and m determines orientation in space. Notice that for hydrogen, only n matters for energy — all the l and m sub-levels with the same n are exactly degenerate, a special symmetry of the 1/r Coulomb potential that disappears in multi-electron atoms.
Nodes are the zeros of the wavefunction — surfaces where the electron has exactly zero probability density. A radial node is a sphere where R_{nl} = 0; there are n−l−1 of them. An angular node is a plane or cone where Y_l^m = 0; there are l of them. Total nodes = n−1. The 2p orbital (n=2, l=1) has zero radial nodes and one angular node (the nodal plane). The 3d orbital (n=3, l=2) has zero radial nodes and two angular nodes. Counting nodes is a powerful consistency check.
A critical conceptual distinction: the wavefunction ψ can take negative values, but this does not mean the electron is "excluded" from those regions. Probability density is |ψ|², which is always non-negative. A negative lobe of ψ is just as accessible to the electron as a positive lobe of the same magnitude. The sign of ψ carries phase information that only becomes observable when two orbitals interact — it determines whether overlap leads to bonding (same sign, constructive) or antibonding (opposite sign, destructive) combinations.
Finally, to correctly describe where the electron actually lives, use the radial probability distribution P(r) = r²|R_{nl}|² rather than |ψ|² alone. The r² factor accounts for the fact that a thin spherical shell of thickness dr has a volume 4πr² dr that grows with radius. For the 1s orbital, the wavefunction amplitude is largest at r = 0, but the most probable radius (the peak of P(r)) is the Bohr radius a₀ = 0.529 Å — because near the nucleus, the small shell volume makes electron detection unlikely despite high amplitude. This is the quantitative picture behind the familiar statement that "electrons occupy orbitals" rather than "electrons sit at the nucleus."