The Radial Equation and Radial Wavefunctions

The explicit hydrogenic radial functions are far more than a textbook exercise; they are the working scaffold of modern atomic and molecular science. Every Slater-type orbital (STO) and every contracted Gaussian-type orbital (GTO) used in quantum chemistry software (Gaussian, ORCA, Psi4, NWChem) is designed to mimic the shape of $R_{nl}(r)$ — the correct nuclear cusp behavior for $s$ functions and the correct node structure for higher shells — because these exact hydrogenic shapes are the best available zeroth-order guess for how an electron distributes itself radially around a nucleus. In spectroscopy, the radial matrix elements built from $R_{nl}$ set the selection rules and intensities of atomic emission and absorption lines — the basis of atomic emission spectroscopy used in astrophysics to determine stellar composition, and in analytical chemistry (ICP-AES) to quantify trace metals. Rydberg atoms — atoms excited to very high $n$ — are giant, fragile systems whose radial wavefunctions can be tuned continuously; they are now engineered in cold-atom physics for quantum computing qubits and for ultra-sensitive electric-field sensors, precisely because the $n^2$ scaling of orbital size derived here makes their properties exquisitely tunable. In chemical bonding, the radial node count and orbital extent explain why 3d orbitals in transition metals are more contracted (relative to 4s) than a naive size argument would suggest — the reduced radial extent of $d$ orbitals compared to $s$ orbitals of similar energy underlies the entire pattern of transition metal chemistry, ligand field splitting, and catalytic reactivity. Finally, NMR and EPR spectroscopists rely on the nonvanishing $s$-orbital density at the nucleus (Fermi contact interaction) — a direct, quantitative consequence of $R_{n0}(0)\neq0$ derived straight from this equation — to interpret hyperfine coupling constants in radicals and to assign chemical shifts.

Ask the tutor: How would the radial equation and its quantization condition change if the electron also felt a small additional $1/r^2$ perturbing potential, and what would happen to the $n$-only energy degeneracy?

The hydrogen atom's Hamiltonian, written in spherical coordinates centered on the nucleus, separates cleanly because the Coulomb potential depends only on $r$, not on the angles $\theta,\phi$. Writing the full wavefunction as a product,

pushes all the angular momentum physics into the spherical harmonics $Y_l^m$, which are eigenfunctions of $\hat{L}^2$ and $\hat{L}_z$. What remains is a single ordinary differential equation in $r$ alone — the radial equation — that still "feels" the angular momentum through a centrifugal term. This is the payoff of separation of variables: a genuinely three-dimensional problem becomes, for each $l$, a one-dimensional problem no harder in structure than a particle in a box or a harmonic oscillator. Physically, an electron near a nucleus experiences two competing effects along the radial direction. The attractive Coulomb pull drags it inward; but unless $l=0$, angular momentum generates an outward centrifugal push, exactly as in classical orbital mechanics. These two effects combine into an effective one-dimensional potential energy curve, $V_{\text{eff}}(r)$, sketched below. Bound states are wavefunctions trapped in the well formed by this effective potential, oscillating or decaying according to whether they sit above or below the local potential energy — precisely the same logic used for any 1D bound-state problem, only now the "particle" is confined by $1/r$ attraction and $1/r^2$ repulsion rather than a well with hard walls. The mathematics of solving this equation mirrors the harmonic oscillator: strip off the correct asymptotic behavior (an exponential decay at large $r$, a power-law behavior at small $r$), leave a polynomial equation for what remains, and discover that a convergent (normalizable) polynomial solution exists only for special, integer-labeled energies. Those polynomials turn out to be the associated Laguerre polynomials, and the integer that labels them becomes the principal quantum number $n$. The chemistry payoff is a complete, exact catalogue of radial wavefunctions $R_{nl}(r)$: functions that predict orbital sizes, the number and location of radial nodes, and — combined with the angular part — the full three-dimensional electron density used to rationalize bonding, spectra, and periodic trends.

The full time-independent Schrödinger equation for an electron of mass $m_e$ bound to a nucleus of charge $+Ze$ is

where: $\nabla^2$ is the Laplacian in spherical coordinates, $Ze^2/4\pi\varepsilon_0 r$ is the Coulomb potential energy between electron and nucleus, and $E$ is the (as-yet unknown) total energy. Separating $\psi = R(r)Y_l^m(\theta,\phi)$ and using the known eigenvalue $\hbar^2 l(l+1)$ of $\hat L^2$ acting on $Y_l^m$ leaves the radial equation:

where: $R(r)$ is the radial wavefunction, $l(l+1)\hbar^2/2m_er^2$ is the centrifugal energy term inherited from separating out the angular part, and the bracket is the effective potential $V_{\text{eff}}(r)$ seen by the radial motion. A standard substitution, $u(r) = rR(r)$, converts this into a form identical to a 1D Schrödinger equation:

where: $u(r)$ behaves like a particle-in-a-box-style wavefunction confined by $V_{\text{eff}}(r)$; the boundary condition $u(0)=0$ is required for $R(r)=u(r)/r$ to remain finite at the origin. Solving this eigenvalue problem (worked out below) yields quantized bound-state energies

where: $n$ is the principal quantum number, and $Z$ is the nuclear charge (so $Z=1$ for hydrogen itself). This reproduces the Bohr formula exactly, but now derived from wave mechanics with $n$ appearing as a normalizability requirement rather than an assumed quantization rule. The normalized radial wavefunctions themselves are

where: $a_0 = 4\pi\varepsilon_0\hbar^2/(m_e e^2)$ is the Bohr radius, $\rho$ is a dimensionless scaled radius, and $L_{n+l}^{2l+1}(\rho)$ is an associated Laguerre polynomial of degree $n-l-1$ in $\rho$. (Sign and normalization conventions for the Laguerre polynomials vary between textbooks; the physically meaningful, convention-independent statement is that $R_{nl}^2$ integrates to 1 against $r^2\,dr$.)

Reading the equations: the energy depends only on $n$, not on $l$ or $m$ — an "accidental" degeneracy special to the pure $1/r$ Coulomb potential. The radial function's shape, however, depends on both $n$ and $l$: the power $\rho^l$ suppresses the wavefunction near the nucleus for higher angular momentum (centrifugal repulsion keeps $p$, $d$, $f$ electrons away from $r=0$), while the polynomial factor injects exactly $n-l-1$ radial nodes.

Example 1 — Foundational: writing down $R_{21}(r)$ for hydrogen and verifying its node count. For hydrogen, $Z=1$. Take $n=2, l=1$ (a 2p orbital). Since $n-l-1 = 0$, the associated Laguerre polynomial is a nonzero constant, so no radial node is expected. The standard result (which can be obtained from the general formula with $\rho = 2r/(2a_0) = r/a_0$) is

Check normalization numerically using $a_0 = 0.529\ \text{Å}$: substitute $x = r/a_0$,

Using $\int_0^\infty x^4e^{-x}dx = 4! = 24$, this gives $24/24 = 1$ exactly, confirming correct normalization. This function has one maximum (found by setting $dR_{21}/dr=0$): differentiating $r e^{-r/2a_0}$ gives a maximum at $r = 2a_0 = 1.058\ \text{Å}$, and no sign change anywhere for $r>0$ — zero radial nodes, exactly as $n-l-1=0$ predicts.

Example 2 — Applied: comparing orbital size for He⁺ versus H using the same $n,l$. He⁺ is a hydrogen-like ion with $Z=2$. Compare the most probable radius of the 1s electron in He⁺ to that in neutral H. For $l=0$, the most probable radius (found by maximizing $r^2R_{n0}^2$, the radial probability density) in the ground state satisfies

For H ($Z=1$): $r_{\text{mp}} = a_0 = 0.529\ \text{Å}$. For He⁺ ($Z=2$): $r_{\text{mp}} = a_0/2 = 0.2645\ \text{Å}$. The orbital shrinks by exactly a factor of 2, tracking the doubled nuclear charge — a direct, quantitative consequence of the $\rho = 2Zr/(na_0)$ scaling built into every $R_{nl}$. The binding energy meanwhile scales as $Z^2$: using $E_n = -13.6\ \text{eV}\cdot Z^2/n^2$,

He⁺'s 1s electron is bound four times more tightly, not twice — illustrating that orbital energy and orbital size respond to $Z$ with different powers, a distinction easy to miss when reasoning loosely about "bigger charge, smaller and lower."

Example 3 — Exam/competition level: locating the radial node of the 3s orbital of Li²⁺ and testing consistency. Li²⁺ is hydrogen-like with $Z=3$. For $n=3,l=0$, the number of radial nodes is $n-l-1 = 2$. The (unnormalized, but functionally exact) form is

(using $Z=3, n=3$, so $\rho = 2\cdot3\cdot r/(3a_0) = 2r/a_0$). Setting the quadratic bracket to zero to find the nodes:

Numerically, $\sqrt{108} = 10.392$, giving

Converting back to $r = \rho\,a_0/2$ (since $\rho = 2r/a_0$):

Two distinct positive nodes, exactly matching the prediction $n-l-1=2$. As a consistency check, both nodes must satisfy $0 < r < \infty$ and be non-degenerate — both conditions hold here, and if the discriminant $324-216=108$ had instead come out negative, that would signal an algebra error, since a valid $3s$ function must have precisely two real, positive radial nodes.

Solving the hydrogen atom's radial equation was the decisive early triumph of wave mechanics, carried out independently and almost simultaneously by Erwin Schrödinger in 1926, using his newly formulated wave equation together with the mathematical machinery of Laguerre polynomials that had already been studied by pure mathematicians for other purposes. Schrödinger's paper "Quantisierung als Eigenwertproblem" showed for the first time that Bohr's quantized energy levels emerge automatically from demanding that a wavefunction be single-valued, finite, and normalizable — no ad hoc quantization postulate required. This was a conceptual shift as large as the formula itself: energy quantization stopped being an assumption and became a mathematical consequence of confinement.