Radial Distribution Functions

Radial distribution functions remain the standard first picture chemists reach for whenever atomic size, shielding, or orbital energy ordering needs explaining. In inorganic and physical chemistry courses and research alike, plots of $P(r)$ for valence ns, np, and (n−1)d orbitals are used to rationalize the anomalous stability of half-filled and filled subshells, the lanthanide contraction (where poor shielding by 4f electrons, visible as a compact, deeply penetrating radial distribution, pulls in the whole valence shell), and the relativistic contraction of s orbitals in heavy elements like gold and mercury, where relativistic effects further sharpen the inner hump of the 6s radial distribution. In spectroscopy and photoionization physics, the radial matrix elements that govern transition intensities and photoionization cross-sections are integrals over exactly this kind of radial density, weighted by additional powers of $r$; the shape of $P(r)$ — how much density sits at small versus large $r$ — directly controls how strongly an electron couples to incoming or outgoing radiation. In computational chemistry, radial distribution functions (and their multi-electron analogues obtained from self-consistent-field or density-functional calculations) are routinely plotted to check whether a basis set correctly reproduces core penetration, since basis-set errors near the nucleus can distort computed hyperfine couplings, NMR shielding constants, and core-ionization energies. Materials scientists use closely related radial concepts — atom-centered pair distribution functions from X-ray or neutron scattering — as the experimental, condensed-phase analogue of the same "probability by shell" logic introduced here for a single atom.

Ask the tutor: Can you show me how the radial distribution function for 2s develops its inner hump and node numerically, step by step from the wavefunction?

Start from a wavefunction you already trust: $\psi_{n\ell m}(r,\theta,\phi) = R_{n\ell}(r)\,Y_\ell^m(\theta,\phi)$. The full probability density is $|\psi|^2 = R_{n\ell}^2(r)\,|Y_\ell^m|^2$, and because the spherical harmonics are themselves normalized over the angular variables, the probability of finding the electron in a differential volume element $r^2\sin\theta\,dr\,d\theta\,d\phi$ is $|\psi|^2$ times that volume element. If you only care about distance from the nucleus — not direction — you integrate away the angles. Since $\int |Y_\ell^m|^2\,d\Omega = 1$, what survives is

This is the probability of finding the electron anywhere in the thin spherical shell between $r$ and $r+dr$, regardless of angle. It is a genuinely different object from $R^2(r)$: the radial density, which tells you the probability per unit volume at a point, versus $P(r)$, the probability per unit radius summed over the whole shell. For every hydrogen 1s-type orbital, $R^2(r)$ is largest at $r=0$ — the density is sharply peaked at the nucleus — yet $P(r)$ vanishes at $r=0$ because the shell volume $4\pi r^2$ vanishes there too. The competition between a decaying exponential ($R^2 \sim e^{-2Zr/na_0}$) and a growing polynomial ($r^2$, or higher powers for larger $\ell$) is exactly what sculpts the characteristic single- or multi-humped shape of $P(r)$. Two structural facts follow directly from the mathematical form of $R_{n\ell}(r)$, which contains an associated Laguerre polynomial of degree $n-\ell-1$ multiplied by $r^\ell$ and $e^{-Zr/na_0}$. First, the number of interior zeros of $R_{n\ell}$ — and hence of $P(r)$ — equals $n-\ell-1$; these are the radial nodes, spherical shells where the probability is exactly zero. Second, near the origin $R_{n\ell}(r) \sim r^\ell$, so $P(r) \sim r^{2\ell+2}$ for small $r$: higher angular momentum pushes the density away from the nucleus, while $s$ orbitals ($\ell=0$) are the only ones with any density directly at the nucleus itself (though even there $P(0)=0$ because of the $r^2$ factor — it is $R^2(0)$ that is nonzero for s states).

The 2s orbital is the simplest case that shows a radial node, and it is worth visualizing before any algebra, because the node is what distinguishes "an electron that is sometimes far away" from "an electron with genuinely zero probability at one particular radius, and two separated humps of probability on either side."

The small inner hump matters chemically far more than its modest area suggests: it represents the ability of a 2s electron to penetrate underneath the 1s shell and feel a larger effective nuclear charge, part of why 2s lies lower in energy than 2p in multi-electron atoms — the physical seed of Effective Nuclear Charge and Shielding.

The governing definitions and results are compact; each is worth committing to memory alongside its meaning.

where: $R_{n\ell}(r)$ is the (real, normalized) radial wavefunction for quantum numbers $n,\ell$; $P(r)$ has units of inverse length so that $P(r)\,dr$ is a dimensionless probability. Reading: $P(r)$ answers "if I could measure the electron's distance from the nucleus many times, what histogram would I get?" — it is the physically measurable radial profile.

where: the integral runs over all radii from the nucleus to infinity. Reading: this is just the normalization of $R_{n\ell}$ itself, restated with the $r^2$ weight made explicit — it guarantees $P(r)$ is a true probability density in $r$.

where: $n$ is the principal quantum number, $\ell$ the azimuthal (angular momentum) quantum number. Reading: radial nodes are spheres of zero probability; this count follows purely from the degree of the associated Laguerre polynomial inside $R_{n\ell}$, and it is exactly what remains once you subtract the $\ell$ angular nodes (nodal planes/cones) from the total $n-1$ nodes of the full orbital.

where: $r_{mp}$ is the most-probable radius, the location of a maximum (not necessarily the global one for multi-humped functions) in $P(r)$. Reading: standard calculus — differentiate $P(r)$, set it to zero, discard trivial or minimum roots, and check the sign of the second derivative or simply compare peak heights when several maxima exist.

where: $\langle r \rangle$ is the mean radius, an expectation value — generally not equal to $r_{mp}$, because $P(r)$ is not symmetric about its peak. Reading: chemists often conflate "most likely distance" and "average distance"; for hydrogen 1s these differ by 50%, a gap you can verify directly.

Example 1 — Foundational: locating the most-probable radius and node count for hydrogen 3p. The normalized 3p radial function (with $Z=1$, hydrogen, in units where $\rho = r/a_0$) is

First, count radial nodes using $n-\ell-1 = 3-1-1=1$: there should be exactly one radial node, and indeed the factor $(6-\rho)$ vanishes once, at $\rho = 6$, i.e. $r = 6a_0$. Everything between $r=0$ and $r=6a_0$ is one sign of $R_{3p}$; everything beyond is the opposite sign (irrelevant once squared, but it confirms the node is a genuine sign change, not just a flattening). To build $P(r)$, square and multiply by $r^2 = a_0^2\rho^2$:

This function is zero at $\rho=0$ (fourth-order contact — very strongly suppressed near the nucleus, consistent with $P(r)\sim r^{2\ell+2}=r^4$ for $\ell=1$), zero again at $\rho=6$ (the radial node), and must therefore have (at least) two humps: a smaller inner one between $0$ and $6$, and a larger outer one beyond $6$. Differentiating $f(\rho) = \rho^4(6-\rho)^2 e^{-2\rho/3}$ and setting $f'(\rho)=0$ leads to a quartic in $\rho$ after factoring out $e^{-2\rho/3}$ and $\rho^3(6-\rho)$ (common factors giving trivial roots at $\rho=0,6$):

Expanding: $24 - 4\rho - 2\rho - \frac{2}{3}(6\rho - \rho^2) = 0 \Rightarrow 24 - 6\rho - 4\rho + \frac{2}{3}\rho^2 = 0$, i.e.

Solving the quadratic:

giving $\rho = 3$ or $\rho = 12$. The root $\rho=3$ lies inside the node ($0<3<6$) and is the smaller, inner-hump maximum; $\rho=12$ lies beyond the node and is the dominant outer-hump maximum, the conventionally quoted "most-probable radius" of 3p, $r_{mp} \approx 12\,a_0$. (Comparing $f(3)$ and $f(12)$ numerically confirms the outer peak is the taller of the two, as expected since the outer lobe of an $n=3$ orbital carries most of the probability.)

Example 2 — Applied: penetration and the 2s/2p energy split in lithium. In hydrogen, all orbitals with the same $n$ are degenerate — 2s and 2p have identical energy. In real multi-electron atoms like lithium (Z=3, configuration $1s^2 2s^1$), they are not, and radial distribution functions explain why. Compare the radial distribution functions (using hydrogen-like shapes as a first approximation, with an effective charge) for 2s and 2p:

The 2s function has a node at $\rho=2$ and a small inner hump for $0<\rho<2$ — probability density that sits inside the 1s core, where the two 1s electrons only partially cancel the nuclear charge. That inner hump means a 2s electron spends part of its time much closer to the +3 nucleus than a 2p electron ever does, since $P_{2p}\propto \rho^4$ vanishes to fourth order at the origin and has no separate inner hump at all — it must build up smoothly from zero with no chance to poke through the core. Because the 2s electron partially penetrates the 1s shielding cloud, it feels a larger effective nuclear charge on average and is therefore more tightly bound (lower in energy) than the 2p electron, even though both have $n=2$. This is precisely why the lithium ground-state configuration is $1s^2\,2s^1$ rather than $1s^2\,2p^1$, and it generalizes into the familiar $ns < np < nd < nf$ energy ordering within a shell for many-electron atoms — the qualitative content of Effective Nuclear Charge and Shielding, derived here purely from the shapes of $P(r)$.

Example 3 — Competition-level: mean radius of a general hydrogen-like orbital, and a numerical check for $n=4,\ell=2$ (4d). A standard closed-form result (derivable via the same Gamma-function machinery as above, using the full associated-Laguerre normalization) gives the general expectation value

Apply this to $n=4,\ell=2$ (a 4d orbital) in a hydrogen-like ion with $Z=3$ (doubly-ionized... actually take $\text{Li}^{2+}$, a genuine one-electron ion with $Z=3$, so the formula applies exactly, no shielding approximations needed):

Now find the number of radial nodes and estimate whether the outermost peak of $P(r)$ lies near this mean. Radial nodes: $n-\ell-1 = 4-2-1=1$, so 4d has exactly one radial node — one inner hump, one dominant outer hump. Because there is only a single subordinate inner hump (not two, as there would be for 4s with $n-\ell-1=3$ nodes), the mean $\langle r\rangle = 7a_0$ should lie reasonably close to, but somewhat larger than, the outer peak $r_{mp}$ — exactly the same qualitative skew seen for 1s (where $\langle r\rangle/r_{mp}=1.5$), just less pronounced because higher-$\ell$ orbitals for a given $n$ have radial distributions concentrated in a narrower band (less inner-hump probability to skew the mean). As a limiting-case sanity check, set $\ell = n-1$ (the maximum possible angular momentum for given $n$ — no radial nodes at all, a single clean hump, e.g. 1s, 2p, 3d, 4f). Then

For $n=1,\ell=0$ this gives $\langle r\rangle = \dfrac{a_0}{Z}\cdot\dfrac{1\cdot 3}{2} = \dfrac{3}{2}\dfrac{a_0}{Z}$ — exactly the 1s result derived above by direct integration, confirming the general formula reduces correctly to the special case already checked independently.

The radial distribution function as a chemically usable object emerged from the marriage of Erwin Schrödinger's wave mechanics with the practical need, pressed hardest by chemists and spectroscopists, to translate abstract $\psi$ functions into pictures of "where the electron is." Schrödinger himself, in his 1926 papers solving the hydrogen atom, obtained the radial functions $R_{n\ell}(r)$ directly; it was the broader community of atomic physicists and chemists — including Linus Pauling and E. Bright Wilson, whose 1935 textbook Introduction to Quantum Mechanics with Applications to Chemistry popularized plots of $r^2R^2$ for exactly the pedagogical reasons developed here — who established the radial distribution plot as the standard chemist's window into atomic structure.