Hydrogenic Wavefunctions and Quantum Numbers

Solving the Coulomb radial equation hands you an energy ladder — but an energy is not a state. What does the electron actually look like? The answer is a complete catalogue: every bound state of every hydrogenic atom — hydrogen, He⁺, Li²⁺ — is one entry $\psi_{n\ell m} = R_{n\ell}(r)\,Y_\ell^m(\theta,\phi)$ , indexed by three integers that work like a postal address: $n$ picks the floor of the energy ladder, $\ell$ the angular shape, $m$ its orientation. Three integers, one state, no exceptions. The catalogue has a non-negotiable grammar. You will never meet a $1p$ or a $2d$ state — the rule $\ell \le n-1$ forbids them, enforced not by decree but by the same series-termination condition that quantized the energy. You can count a state's radial nodes without plotting anything: exactly $n-\ell-1$ . And the radial functions are built from a single polynomial family, the associated Laguerre polynomials, whose two-index bookkeeping is the most common source of wrong answers in the chapter — this article pins the conventions down and shows you how to audit any table entry yourself. The $(n, \ell, m)$ labels are the vocabulary of every atomic spectrum, every periodic-table argument, and every quantum-chemistry code: when a spectroscopist says "the $3d$ level" or a chemist draws a $p_x$ orbital, they are speaking the language assembled here.

Key takeaways

Every hydrogenic bound state factors as $\psi_{n\ell m}(r,\theta,\phi) = R_{n\ell}(r)\,Y_\ell^m(\theta,\phi)$ : a Coulomb radial part times a universal angular part.
Three integers label everything: $n = 1, 2, 3, \dots$ (energy), $\ell = 0, 1, \dots, n-1$ (angular momentum), $m = -\ell, \dots, +\ell$ (orientation). The bound $\ell \le n-1$ comes from series termination, not from an extra postulate.
The radial functions are $R_{n\ell} \propto e^{-Zr/na_0}\,(2Zr/na_0)^\ell\,L^{2\ell+1}_{n-\ell-1}(2Zr/na_0)$ — associated Laguerre polynomials dressed with the near-origin rise $r^\ell$ and the bound-state tail.
Radial nodes number exactly $n-\ell-1$ ; angular nodes number $\ell$ ; the total is always $n-1$ .
The energy $E_n = -13.6\ \mathrm{eV}\,Z^2/n^2$ ignores $\ell$ and $m$ entirely, so shell $n$ holds $n^2$ degenerate orbital states ( $2n^2$ with spin).
Laguerre conventions differ between textbooks by factorials and a sign — the bridge is derived here.
Key idea: one polynomial family plus one integer grammar generates every bound state of every one-electron atom.

The core insight

Separability turns the 3D hydrogenic problem into a universal angular part (spherical harmonics, labeled by $\ell, m$ ) times a Coulomb radial part (Laguerre polynomials, labeled by $n, \ell$ ); the polynomial-termination condition that makes $R_{n\ell}$ normalizable simultaneously quantizes the energy as $E_n \propto -Z^2/n^2$ and enforces $0 \le \ell \le n-1$ .

Where this lives today

The $(n, \ell, m)$ catalogue is the working coordinate system of computational atomic physics: quantum-chemistry packages expand molecular orbitals in hydrogen-like basis functions (Slater-type orbitals carry exactly the $r^\ell e^{-\zeta r}$ anatomy of equation (6)), and when sympy or Mathematica hands you `R_nl(3, 1, r)` it is evaluating (6) in the modern Laguerre convention — the verification habit taught here is a daily-practice defense, because old-convention tables still circulate. At the exotic end, circular Rydberg states ( $\ell = n-1$ , $n \approx 50$ ) are laboratory objects — micron-sized, node-free, millisecond-lived microwave-photon sensors and qubit candidates: the single-humped $R_{n,n-1}$ , peaked at $n^2a_0$ , is quantum mechanics' closest reconstruction of Bohr's orbit.

The physical picture

Figure

The quantum-number address book — every hydrogenic state on one grid

Read the grid row by row. Each shell $n$ opens exactly $n$ subshell doors, $\ell = 0$ through $n-1$ — the staircase widens by one door per floor, never overhanging — and each door hides $2\ell+1$ rooms, one per $m$ . All rooms on a floor share one energy: the row heights are the measured $-13.6\ \mathrm{eV}/n^2$ ladder, and adding across a row gives $1, 4, 9, 16$ — the $n^2$ count derived at The Hydrogen Atom: Coulomb Radial Equation, and explained at Accidental Degeneracy and the SO(4) Symmetry. Behind each address sits a concrete function: the angular factor comes from the fixed catalogue of Spherical Harmonics; the radial factor is where hydrogen's own physics lives:

Figure

Anatomy of a radial function — R₃₁ dissected into its three factors

Every $R_{n\ell}$ is the same three-part story. Near the origin it rises as $r^\ell$ — the centrifugal barrier's fingerprint, why only s-states touch the nucleus. Far away it dies as $e^{-Zr/na_0}$ , the universal bound-state tail whose reach $na_0/Z$ grows with the shell. In between sits a polynomial of degree $n-\ell-1$ with exactly that many positive roots — the radial nodes. Walk the s-ladder and watch the rule tick: $1s$ has none, $2s$ crosses zero once (at $r = 2a_0$ exactly), $3s$ twice ( $1.90\,a_0$ and $7.10\,a_0$ ) — each new floor buys one more sign change, exactly as the oscillation theory of Node Theorem and Energy Ordering predicts on the half-line. Keep the three-factor anatomy in mind and you can sketch any hydrogenic state — or debug any table — without solving anything.

Assembling the state

The separation of variables established at Separation of Variables in Spherical Coordinates dictates the product form; hydrogen's job is only to supply the radial factor. The complete bound state reads

\psi_{n\ell m}(r,\theta,\phi) = R_{n\ell}(r)\,Y_\ell^m(\theta,\phi) \qquad (1)

where: $R_{n\ell}$ — the Coulomb radial function carrying all the dynamics; $Y_\ell^m$ — the universal angular factor owned by Spherical Harmonics. Note the bookkeeping: $R$ never sees $m$ , $Y$ never sees $n$ — the factors communicate only through the shared $\ell$ , which is why the quantum numbers interlock at all. The three labels have sharply different jobs; their ranges interlock:

n = 1, 2, 3, \dots \qquad \ell = 0, 1, \dots, n-1 \qquad m = -\ell, \dots, +\ell \qquad (2)

where: $n$ — the principal quantum number, sole owner of the energy $E_n = -13.606\ \mathrm{eV}\,Z^2/n^2$ ; $\ell$ — the orbital quantum number, fixing $L^2 = \hbar^2\ell(\ell+1)$ ; $m$ — the magnetic quantum number, fixing $L_z = \hbar m$ . The ranges of $\ell$ and $m$ were settled upstream (integer character at Eigenfunctions of L_z and the Azimuthal Quantization and Associated Legendre Functions; the $m$ -range at Ladder Operators and the Angular Momentum Spectrum). What is new — special to the Coulomb problem — is the upper bound on $\ell$ . That bound is inherited, not imposed. The series at The Hydrogen Atom: Coulomb Radial Equation terminates when $n = j_{\max} + \ell + 1$ , with $j_{\max} \ge 0$ the surviving polynomial's degree. Read it backwards: for given $n$ ,

n_r \equiv n - \ell - 1 \ge 0 \quad \Longleftrightarrow \quad \ell \le n - 1 \qquad (3)

where: $n_r$ — the radial quantum number, the number of radial nodes. Equation (3) is the address book's entire grammar: $n = n_r + \ell + 1$ is a sum of radial oscillations plus angular momentum plus one. A $2d$ state would need $n_r = -1$ nodes — no such polynomial, no such state. Spectroscopy compresses $(n, \ell)$ into the notation you will use forever: the number is $n$ , the letter encodes $\ell$ — $s, p, d, f$ for $\ell = 0, 1, 2, 3$ (then $g, h, \dots$ alphabetically, skipping $j$ ). So $3d$ means $(3, 2)$ : zero radial nodes, five $m$ -rooms. The letters are fossils of pre-quantum line-series taxonomy (sharp, principal, diffuse, fundamental), but the code stuck.

The Laguerre machinery

The polynomial family that survives termination is not arbitrary: it is the associated Laguerre family, the classical orthogonal polynomials of the half-line with weight $x^\alpha e^{-x}$ . The modern definition (used by Griffiths 3rd edition, by every computer-algebra system, and by this article) is

L_q^{\alpha}(x) = \frac{x^{-\alpha}e^{x}}{q!}\,\frac{d^{q}}{dx^{q}}\!\left(e^{-x}x^{q+\alpha}\right) \qquad (4)

where: $q$ — the polynomial degree; $\alpha$ — the superscript tracking the weight. This is the Rodrigues form; the members you need most are $L_0^{\alpha} = 1$ and $L_1^{\alpha}(x) = 1+\alpha-x$ . Each $L_q^\alpha$ has degree $q$ with $q$ real positive roots, and the family's defining property is weight-orthogonality:

\int_0^\infty x^{\alpha}e^{-x}\,L_q^{\alpha}(x)\,L_{q'}^{\alpha}(x)\,dx = \frac{(q+\alpha)!}{q!}\,\delta_{qq'} \qquad (5)

where: the right-hand normalization $(q+\alpha)!/q!$ is what fixes every constant in the atom. For the hydrogen case $\alpha = 2\ell+1$ , $q = n-\ell-1$ : check the smallest instance yourself — $\int_0^\infty x^2 e^{-x}\,[L_2^2(x)]^2\,dx = 4!/2! = 12$ . Dressing the polynomial with the two asymptotic factors and normalizing gives the master formula of the article:

R_{n\ell}(r) = \sqrt{\left(\frac{2Z}{na_0}\right)^{3}\frac{(n-\ell-1)!}{2n\,(n+\ell)!}}\;e^{-Zr/na_0}\left(\frac{2Zr}{na_0}\right)^{\ell}L^{2\ell+1}_{n-\ell-1}\!\left(\frac{2Zr}{na_0}\right) \qquad (6)

where: $a_0$ — the Bohr radius (its anatomy at The Bohr Radius and Characteristic Atomic Scales); $Z$ — the nuclear charge, entering only through the combination $Zr/a_0$ (the one-substitution scaling story is Hydrogenic Ions and Scaling with Nuclear Charge's). With the convention (4), equation (6) is positive at small $r$ for every state — $R_{n\ell} \propto +r^\ell$ near the origin — which is the standard sign choice and the one every table below obeys. Be warned that (6) is convention-locked: change the Laguerre definition and both the constant and the sign change with it. That trap gets its own derivation card below.

Derivation ANormalizing the radial functions

The goal is the square root in equation (6): starting from $R_{n\ell} = N\,e^{-\rho/2}\rho^\ell L^{2\ell+1}_{n-\ell-1}(\rho)$ with $\rho = 2Zr/(na_0)$ , fix $N$ so that $\int_0^\infty R_{n\ell}^2\,r^2\,dr = 1$ — the normalization demanded by Boundary Conditions and Normalization of Radial Wavefunctions. The strategy is a change of variables into the Laguerre orthogonality integral (5). The one subtlety: the radial measure $r^2dr$ produces the weight $x^{\alpha+1}$ , one power higher than the orthogonality weight, so we need the first-moment cousin of (5) — that extra power is exactly where the $2n$ in the denominator comes from. Watch for it in Step 3.

Step 1. Change variables. With $\rho = 2Zr/(na_0)$ , $r = (na_0/2Z)\rho$ and $r^2\,dr = (na_0/2Z)^3\rho^2\,d\rho$ , so

\int_0^\infty R_{n\ell}^2\,r^2\,dr = N^2\left(\frac{na_0}{2Z}\right)^{3}\int_0^\infty e^{-\rho}\,\rho^{2\ell+2}\left[L^{2\ell+1}_{n-\ell-1}(\rho)\right]^2 d\rho.

Why: all the physics scales ( $Z$ , $n$ , $a_0$ ) now sit in one prefactor; what remains is a pure Laguerre integral with $\alpha = 2\ell+1$ but weight power $\alpha + 1$ .

Step 2. Reach for the right integral — the first-moment formula

\int_0^\infty e^{-x}x^{\alpha+1}\left[L_q^{\alpha}(x)\right]^2 dx = \frac{(q+\alpha)!}{q!}\,(2q+\alpha+1).

Why: it follows from (5) plus the three-term recurrence — multiplying by $x$ shifts $L_q^\alpha$ into a combination of $L_{q-1}, L_q, L_{q+1}$ , and orthogonality kills the off-diagonal pieces.

Step 3. Evaluate the bookkeeping. Here $q = n-\ell-1$ and $\alpha = 2\ell+1$ , so $q + \alpha = n+\ell$ and $2q + \alpha + 1 = 2(n-\ell-1) + 2\ell + 2 = 2n$ :

\int_0^\infty e^{-\rho}\rho^{2\ell+2}\left[L^{2\ell+1}_{n-\ell-1}\right]^2 d\rho = \frac{(n+\ell)!}{(n-\ell-1)!}\,2n.

Why: the promised $2n$ appears as $2q+\alpha+1$ — the "extra power of $\rho$ " premium. Note how tidy the Coulomb indices are: $q+\alpha$ collapses to $n+\ell$ with no leftover.

Step 4. Solve for $N$ :

N^2\left(\frac{na_0}{2Z}\right)^{3}\frac{(n+\ell)!\,2n}{(n-\ell-1)!} = 1 \;\;\Longrightarrow\;\; N = \sqrt{\left(\frac{2Z}{na_0}\right)^{3}\frac{(n-\ell-1)!}{2n\,(n+\ell)!}},

which is equation (6). The positive root is chosen so that $R_{n\ell}(r) \to +N\rho^\ell$ as $r \to 0$ — all table entries start positive. Verified: equation (6) reproduces sympy's `R_nl` symbolically for every $(n,\ell)$ with $n \le 4$ at $Z = 1$ and $Z = 3$ , and each tabulated $R_{n\ell}$ integrates to exactly 1 (verify script, sections 1 and 3).

The formula in hand, here is the working table — the six states you will use ninety percent of the time. Every entry has been machine-checked against equation (6), symbol by symbol:

state	$R_{n\ell}(r)$	nodes $n{-}\ell{-}1$
$1s$	$2\,a_0^{-3/2}\,e^{-r/a_0}$	0
$2s$	$\dfrac{1}{2\sqrt{2}}\,a_0^{-3/2}\left(2 - \dfrac{r}{a_0}\right)e^{-r/2a_0}$	1
$2p$	$\dfrac{1}{2\sqrt{6}}\,a_0^{-3/2}\,\dfrac{r}{a_0}\,e^{-r/2a_0}$	0
$3s$	$\dfrac{2}{81\sqrt{3}}\,a_0^{-3/2}\left(27 - 18\dfrac{r}{a_0} + 2\dfrac{r^2}{a_0^2}\right)e^{-r/3a_0}$	2
$3p$	$\dfrac{8}{27\sqrt{6}}\,a_0^{-3/2}\,\dfrac{r}{a_0}\left(1 - \dfrac{r}{6a_0}\right)e^{-r/3a_0}$	1
$3d$	$\dfrac{4}{81\sqrt{30}}\,a_0^{-3/2}\,\dfrac{r^2}{a_0^2}\,e^{-r/3a_0}$	0

For $Z \ne 1$ : replace $a_0 \to a_0/Z$ and multiply by $Z^{3/2}$ — the scaling story belongs to Hydrogenic Ions and Scaling with Nuclear Charge.

Go deeper · Derivation B · 2 minTwo Laguerre alphabets, one atom — and the sign hiding in the conventions

Derivation BTwo Laguerre alphabets, one atom — and the sign hiding in the conventions

The goal is to defuse the most common table-lookup error in the chapter: textbooks define the associated Laguerre polynomials in two incompatible ways, and formulas migrate between books with wrong constants — or a wrong overall sign — when the convention silently changes. We derive the bridge between the conventions and show what a correct radial formula looks like in each. This matters concretely: this article's own source printed equation (6) with a leading minus sign, an artifact of mixing the two alphabets (edition note in Going further). The strategy: define both conventions, prove the conversion identity, translate the formula both ways.

Step 1. State the two definitions. The modern convention is equation (4): $L_q^{\alpha}$ has degree $q$ , with $1/q!$ built in. Many older texts (Pauling–Wilson, Schiff, Messiah, Griffiths 2nd ed.) instead build from the unnormalized $\bar L_k(x) = e^x \frac{d^k}{dx^k}(e^{-x}x^k)$ and differentiate: $\bar L^{s}_{k}(x) = \frac{d^{s}}{dx^{s}}\bar L_k(x)$ , an object of degree $k-s$ whose subscript no longer equals its degree. Why: the old subscript tracks the parent polynomial, the modern one the degree — the first thing to check in any book.

Step 2. Derive the bridge. Both objects solve the same Laguerre ODE and have the same degree $k-s$ , so they are proportional; comparing leading coefficients gives

\bar L^{s}_{k}(x) = (-1)^{s}\,k!\;L^{s}_{k-s}(x).

The two alphabets differ by a factorial and by a sign $(-1)^s$ — and for hydrogen $s = 2\ell+1$ is always odd, so the sign is always $-1$ . Why: this is the entire mechanism by which a correct old-convention formula grows a minus sign that looks wrong in modern notation. For hydrogen the sign never hides: odd superscript, guaranteed flip.

Step 3. Translate the radial formula. In the old alphabet the hydrogen polynomial is $\bar L^{2\ell+1}_{n+\ell}$ — degree $n-\ell-1$ , as required — and the correct normalization is

R_{n\ell} = -\sqrt{\left(\frac{2Z}{na_0}\right)^{3}\frac{(n-\ell-1)!}{2n\,[(n+\ell)!]^{3}}}\;e^{-Zr/na_0}\,\rho^\ell\,\bar L^{2\ell+1}_{n+\ell}(\rho)

— note the cube on $(n+\ell)!$ and the leading minus, which together exactly cancel the $(-1)\,(n+\ell)!$ of the bridge and reproduce the positive modern formula (6). Why: the minus sign and the cubed factorial travel together. A formula with the minus but a first-power $(n+\ell)!$ — or with the cube but no minus — is a hybrid of the two conventions and is simply wrong by a sign or by factorials.

Step 4. Fix the sanity anchor: whatever the alphabet, $R_{10} = +2a_0^{-3/2}e^{-r/a_0}$ . One evaluation settles any table's convention in ten seconds. Verified: the bridge identity checked symbolically for $(k,s) = (2,1), (3,1), (4,2), (5,3)$ ; the old-convention formula of Step 3 shown equal to the modern (6) for all $n \le 3$ ; the hybrid (minus sign with first-power factorial) shown equal to $-R_{n\ell}$ , contradicting the positive $1s$ anchor (verify script, section 2).

Nodes, orthogonality, and the freedom inside a shell

How much structure can you read off an address without computing? Almost all of it. The polynomial $L^{2\ell+1}_{n-\ell-1}$ has degree $n-\ell-1$ , and all its roots are real, positive, and simple — each a radius where $R_{n\ell}$ changes sign:

\text{radial nodes} = n - \ell - 1, \qquad \text{angular nodes} = \ell, \qquad \text{total} = n - 1 \qquad (7)

where: angular nodes are the $\ell$ nodal cones and planes of $Y_\ell^m$ . The total $n-1$ depends on $n$ alone: energy buys nodes, and each shell spends its allowance differently — $3s$ puts both in the radius, $3d$ both in the angles, $3p$ splits them. The catalogue is also orthonormal: $\langle \psi_{n\ell m} | \psi_{n'\ell'm'} \rangle = \delta_{nn'}\delta_{\ell\ell'}\delta_{mm'}$ . The spherical harmonics kill any mismatch in $(\ell, m)$ ; what remains — same $\ell$ , different $n$ — is the radial statement

\int_0^\infty R_{n\ell}(r)\,R_{n'\ell}(r)\,r^2\,dr = \delta_{nn'} \qquad (8)

where: the weight $r^2$ is essential — it is the Sturm–Liouville weight of the radial problem, the same measure that Boundary Conditions and Normalization of Radial Wavefunctions built into the normalization doctrine (its Sturm–Liouville argument covers exactly this integral; we consume it rather than re-prove it). Note what (8) does not say: radial functions of different $\ell$ — say $R_{21}$ and $R_{31}$ versus $R_{20}$ — are orthogonal via the angular factor, not the radial integral; $\int R_{10}R_{21}\,r^2dr$ by itself vanishes for no reason and generally doesn't. One more freedom hides inside the shell: all $n^2$ states share one energy, so any normalized combination is an equally good stationary state — the eigenbasis within a degenerate level is a choice, the labeling freedom formalized by Simultaneous Diagonalization and Complete Sets of Commuting Observables. Chemistry exercises it constantly: the real orbitals $p_x = (Y_1^{-1} - Y_1^{1})/\sqrt 2$ and $p_y = i(Y_1^{-1} + Y_1^{1})/\sqrt 2$ re-basis the $m = \pm 1$ pair, trading definite $L_z$ for definite direction. Every chemistry dumbbell is a superposition of two perfectly good $m$ -eigenstates.

Limiting cases

Case	Result	Meaning
Ground state	$\psi_{100} = \pi^{-1/2}(Z/a_0)^{3/2}e^{-Zr/a_0}$	nodeless, spherical, positive — the anchor that fixes every sign convention
$r \to 0$	$R_{n\ell} \propto +r^{\ell}$	centrifugal suppression: only s-states reach the nucleus (contact terms, hyperfine)
$r \to \infty$	$R_{n\ell} \propto r^{n-1}e^{-Zr/na_0}$	universal exponential tail with reach $na_0/Z$ ; the polynomial's top power sets the prefactor
Circular states $\ell = n-1$	$R_{n,n-1} \propto r^{n-1}e^{-Zr/na_0}$ , no radial nodes	single-humped, peaked at $n^2a_0/Z$ — the most Bohr-like orbits
Nuclear charge $Z$	$\psi^{(Z)}_{n\ell m}(\mathbf r) = Z^{3/2}\psi^{(1)}_{n\ell m}(Z\mathbf r)$	one substitution rescales the whole catalogue — the story of Hydrogenic Ions and Scaling with Nuclear Charge
Large $n$	$\langle r \rangle \sim n^2 a_0$ , spacing $\to 0$	Rydberg states: micron-sized atoms, quasi-classical shells

Common trapn counts the shell, NOT the nodes — and the Laguerre subscript is neither

Three integers orbit every radial function and students routinely permute them. $n$ is the shell label: it sets the energy and equals nodes-plus- $\ell$ -plus-one. The radial node count is $n-\ell-1$ — a $4p$ state has 2 radial nodes, not 4. The polynomial degree (the modern Laguerre subscript) is that same $n-\ell-1$ , not $n$ : the $3d$ function contains $L_0^5$ , a constant — no radial wiggles at all. Worse, older tables write the subscript as $n+\ell$ (Derivation B's old convention), so the "same" symbol names different polynomials in different books. Before using any printed $R_{n\ell}$ : check the $1s$ sign anchor and confirm the polynomial degree equals $n-\ell-1$ . Ten seconds of bookkeeping saves an evening of sign-chasing.

Worked examples

1 · Build a state from the catalogue (confidence builder)

Assemble $\psi_{210}$ explicitly, verify its normalization over all of 3D space, and census its nodes.

Step 1. Look up the parts. The address $(n, \ell, m) = (2, 1, 0)$ needs $R_{21} = \frac{1}{2\sqrt 6}a_0^{-3/2}(r/a_0)e^{-r/2a_0}$ from the table and $Y_1^0 = \sqrt{3/4\pi}\cos\theta$ from the spherical-harmonic catalogue.

Step 2. Multiply and tidy. $\psi_{210} = R_{21}Y_1^0 = \frac{1}{2\sqrt 6}\sqrt{\frac{3}{4\pi}}\,a_0^{-3/2}\frac{r}{a_0}e^{-r/2a_0}\cos\theta$ . The constants combine: $\frac{1}{2\sqrt 6}\cdot\frac{\sqrt 3}{2\sqrt\pi} = \frac{1}{4\sqrt{2\pi}}$ , so

\psi_{210} = \frac{1}{4\sqrt{2\pi}}\,a_0^{-3/2}\,\frac{r}{a_0}\,e^{-r/2a_0}\cos\theta.

Step 3. Verify normalization over $d^3r$ , not just the radial integral: the angular part gives $2\pi\cdot\tfrac23 = \tfrac{4\pi}{3}$ , the radial part $\frac{1}{32\pi a_0^5}\int_0^\infty r^4e^{-r/a_0}dr = \frac{4!}{32\pi} = \frac{3}{4\pi}$ . Product: $\mathbf{1}$ . ✓

Step 4. Census the nodes without plotting: radial $n-\ell-1 = 0$ ( $L_0^3 = 1$ , a constant); angular $\ell = 1$ — the plane $\cos\theta = 0$ . Total $1$ . The $2p_z$ dumbbell: no radial zero, one nodal plane.

2 · Audit a table entry (exam level)

A handbook lists $R_{32}(r) = \frac{4}{81\sqrt{30}}a_0^{-3/2}(r/a_0)^2e^{-r/3a_0}$ . Audit it: confirm the structure against the master formula, then predict every feature of the $3d$ and $4d$ states a spectroscopist would ask for.

Step 1. Structure check against (6): for $(3,2)$ the polynomial index is $0$ , so $L_0^5 = 1$ — pure $r^2\times$ exponential, no wiggle ✓; rise $r^\ell = r^2$ ✓; tail $e^{-r/3a_0}$ matches $na_0$ ✓.

Step 2. Constant check: equation (6) gives $N = \sqrt{(2/3a_0)^3\cdot\frac{0!}{2\cdot3\cdot5!}}$ , and carrying the $(2/3)^2$ from $\rho^2 = (2r/3a_0)^2$ into the prefactor reproduces $\frac{4}{81\sqrt{30}} = 9.02\times10^{-3}$ exactly; the machine check confirms the entry equals (6) and integrates to 1. Entry verified.

Step 3. Spectroscopist's census. $3d$ : radial nodes $0$ , angular nodes $2$ , degeneracy $2\ell+1 = 5$ , energy $E_3 = -13.606/9 = \mathbf{-1.512\ eV}$ . It is a circular state ( $\ell = n-1$ ): single-humped, peaking at $n^2a_0 = 9a_0$ (set $\frac{d}{dr}[r\,R_{32}]^2 = 0$ — the maximum lands exactly at $9a_0$ ).

Step 4. Now $4d$ from pure grammar: radial nodes $4-2-1 = \mathbf{1}$ , five $m$ -states, $E_4 = -13.606/16 = \mathbf{-0.850\ eV}$ , total nodal surfaces $3$ . No integral computed — the address did all the work.

Origins · 1 minEdmond Laguerre and how this was discovered

Edmond Laguerre

1834–1886 · Paris

Edmond Laguerre never heard of quantum mechanics, atoms, or the electron — he died a year before J.J. Thomson's cathode-ray work began. A geometer at heart, he served as an artillery officer, then taught at the École Polytechnique, publishing over 140 papers while chronically ill. His polynomials emerged from an 1879 memoir on the integral $\int e^{-x}/x\,dx$ and continued fractions — approximation theory, as far from spectroscopy as mathematics gets; he proved the family orthogonal under the weight $e^{-x}$ on the half-line. Forty years after his death, Schrödinger's radial equation was found to terminate on exactly this family: the associated Laguerre polynomials became the radial alphabet of every atom in the universe, computed daily in every quantum-chemistry code. Laguerre joins Clebsch in the club of mathematicians drafted into physics decades after death — and even his polynomials' definition kept evolving without him: the convention chaos this article untangles is the residue of a family that outgrew its author's notation.

History

1879

— Laguerre studies the orthogonal polynomial family arising from continued-fraction expansions of the exponential integral.

1926

— Schrödinger solves the hydrogen atom; the radial solutions terminate on associated Laguerre polynomials.

1928

— Hartree begins systematic numerical tables of atomic radial functions, hydrogen's analytic forms as the benchmark.

1935

— Condon and Shortley's Theory of Atomic Spectra standardizes spectroscopic notation and phase conventions for atomic states.

1951

— Erdélyi's Higher Transcendental Functions codifies the modern Laguerre normalization, the convention computer algebra later adopted.

Going further

Texts. Griffiths & Schroeter, Introduction to Quantum Mechanics, 3rd ed., §4.2 · Sakurai & Napolitano, Modern Quantum Mechanics, 2nd ed., §3.7 · Shankar, Principles of Quantum Mechanics, 2nd ed., Ch. 13 · Cohen-Tannoudji, Diu & Laloë, Vol. 1, Ch. VII · Townsend, A Modern Approach to Quantum Mechanics, 2nd ed., Ch. 10.

Provenance. Equations cross-checked against Griffiths 3rd ed. §4.2 (radial wavefunctions Eq. 4.75, energy Eq. 4.70, expectation values Table 4.4), Sakurai & Napolitano 2nd ed. §3.7, Shankar 2nd ed. Ch. 13, and Cohen-Tannoudji Vol. 1 Ch. VII. Verifiers should recompute: (1) the radial recurrence and termination giving $n=j_{\max}+\ell+1$ ; (2) the degeneracy sum $\sum_{\ell=0}^{n-1}(2\ell+1)=n^2$ ; (3) standard constants $\mathrm{Ry}=13.606\,\mathrm{eV}$ and $a_0=0.529\times10^{-10}\,\mathrm{m}$ (CODATA); (4) $\langle r\rangle$ and $\langle 1/r\rangle$ formulas via the Kramers relation or direct Laguerre integrals. [verbatim from the AllQFT export] Edition note (2026-07-08). The AllQFT base's radial-wavefunction card printed equation (6) with a leading minus sign while using the modern factorial structure $(n-\ell-1)!/(n+\ell)!$ and the modern subscript $L^{2\ell+1}_{n-\ell-1}$ — an inconsistent hybrid of the two Laguerre conventions (the minus belongs to the old $[(n+\ell)!]^3$ / $\bar L^{2\ell+1}_{n+\ell}$ convention, as Derivation B shows), and one that contradicts the base's own positive $1s$ limit card. Machine verification: the base's card equals $-R_{n\ell}$ (sympy) for every $n \le 4$ ; the corrected positive form (6) is what this article ships. The termination derivation itself is owned by The Hydrogen Atom: Coulomb Radial Equation and consumed here; radial moments $\langle r^k \rangle$ and the base's expectation-value card are owned by Radial Probability Distribution and Expectation Values.

Prerequisites. Separation of variables in spherical coordinates · the radial equation · effective potential and the centrifugal barrier · the hydrogen Coulomb radial equation. Leads to. The Bohr energy spectrum and Rydberg constant · degeneracy and the role of the magnetic quantum number · accidental degeneracy and the SO(4) symmetry · radial probability distribution and expectation values · hydrogenic ions and scaling with nuclear charge.

Is n a physical observable or just a bookkeeping index?

Physical: $n$ labels the Hamiltonian eigenvalue — every spectral line is a difference of two $E_n$ . But its decomposition $n = n_r + \ell + 1$ is bookkeeping for how the state spends its energy: radial oscillations ( $n_r$ ) or angular momentum ( $\ell$ ). Two states of the same $n$ are identical in energy (until fine structure) yet geometrically opposite — compare the two-node $3s$ sphere with the node-free $3d$ doughnut.

Why does the radial function only rise as r^ℓ near the origin — where exactly does that come from?

From the centrifugal term $\hbar^2\ell(\ell+1)/2mr^2$ : near $r=0$ it dominates everything, and the regular solution behaves as $r^\ell$ (the indicial analysis at The Radial Equation). Angular momentum is an inertial wall — a state with $\ell \ne 0$ cannot sit on the axis it rotates about. The consequences are outsized: only s-electrons have contact density at the nucleus, which is why hyperfine structure, electron capture, and isotope shifts are s-state phenomena.

Do the R_nl for different n but the same ℓ overlap? They look similar.

Exactly orthogonal — but only under the $r^2$ weight, equation (8): the integral balances positive and negative regions to precisely zero, the nodes placed by the Sturm–Liouville structure to make that possible. That is why node counts strictly increase with $n$ — a function orthogonal to all its predecessors under a positive weight needs one more sign change.

Are the p_x and p_y orbitals of chemistry eigenstates of the hydrogen Hamiltonian?

Yes of energy, no of $L_z$ . They are exact eigenstates of $\hat H$ and $\hat L^2$ — degeneracy makes any shell-internal mix stationary — but carry no definite $m$ : measure $L_z$ on $p_x$ and you get $\pm\hbar$ with equal probability. Chemistry prefers them because bonds break rotational symmetry along definite directions; spectroscopy in a magnetic field prefers the $m$ -basis. Same shell, two bases, one physics.

The table's 2s function is negative for r > 2a₀. Is a negative wavefunction allowed?

Entirely — only $|\psi|^2$ is probability. But the sign structure is not decoration: it makes orthogonality (8) possible, and relative signs control interference when states superpose (bonding vs antibonding orbitals are a relative-sign story). What is conventional is each $R_{n\ell}$ 's overall sign; the standard choice makes $R_{n\ell}(0^+) > 0$ , and Derivation B's convention chaos is that choice made inconsistently.

Does the catalogue ψ_nlm cover ALL states of hydrogen?

No — complete for the bound spectrum only. Hydrogen also has a continuum of unbound states ( $E > 0$ ), and a complete basis needs both: expanding an arbitrary wavepacket takes the bound sum plus a continuum integral. The continuum's radial functions solve the same equation with oscillating asymptotics — free spherical waves distorted by the Coulomb tail. Forgetting the continuum is a classic error in sum-over-states quantities like polarizabilities.