Deriving Equation (5.61) via Laguerre Polynomials

Equation (5.61) (in the numbering of Griffiths' textbook, and analogous numberings elsewhere) gives the completely normalized radial wavefunction $R_{n\ell }(r)$ of the hydrogen atom. After the radial Schrödinger equation is solved by a power-series/asymptotic-factorization method, the polynomial part of the solution is identified as an associated Laguerre polynomial $L_{n-\ell -1}^{2\ell +1}$. The remaining task — fixing the overall multiplicative constant so that ${\int }_0^{\infty }|R_{n\ell }{|}^2{r}^{2}dr=1$ — requires a specific weighted integral of Laguerre polynomials. The surprising fact is that this integral is not simply the orthogonality relation; it involves an extra factor of the argument, producing the characteristic $2n$ and $(n+\ell )!$ structure. Mastering this derivation is the canonical bridge between special-function theory and physical normalization.

The Core Insight

The radial wavefunction is a Laguerre polynomial times an exponential and a power. Its squared norm reduces to the integral

Formal Statement

Associated Laguerre Polynomials

The associated (generalized) Laguerre polynomials are defined by Rodrigues' formula

equivalently by the finite sum

Here $q=n-\ell -1\ge 0$ is the degree and $p=2\ell +1$ is the order for the hydrogen problem.

(Caution: a competing convention writes $L_n^k$ with different factorial normalization. Griffiths uses [1]–[2]; physicists must always check which convention multiplies their normalization integral. See Laguerre Polynomial Conventions.)

The Hydrogen Radial Equation in Dimensionless Form

With $a=4\pi {ϵ}_0{\hslash }^2/m{e}^2$ the Bohr radius, the bound-state energies are $E_n=-[m{e}^4/(32{\pi }^2{ϵ}_0^2{\hslash }^2)]/n^2$. Defining $\kappa =\sqrt{-2m{E}_n}/\hslash =1/(na)$ and $\rho =2\kappa r=2r/(na)$, the radial function written as $R=u/r$ satisfies a reduced equation whose normalizable solution is

Equation (5.61): The Normalized Radial Wavefunction

The normalized hydrogen radial wavefunctions are

normalized so that

The leading minus sign is a phase convention (Griffiths' choice) and carries no physical content.

Radial probability density $r^2 R_{n\ell}^2$ for the hydrogen ground state $n{=}1,\ell{=}0$, showing the Bohr-radius peak

Why It Holds / Derivation

From Normalization Integral to Equation (5.61)

Reinstating ${\rho }^{\ell }=(2r/na{)}^{\ell }$ and $e^{-\rho /2}=e^{-r/na}$, and choosing the conventional overall minus sign, yields exactly Equation (5.61), i.e. [4]. $■$

The Integral is Not the Orthogonality Relation

A frequent fatal error is to use the orthogonality integral

directly in Step 2. But the radial measure $r^{2}dr$ produces $x^{p+1}$, one power higher than the orthogonality weight $x^p$. That extra $x$ is precisely what generates the factor $2n=2q+p+1$. Forgetting it loses the $2n$ and gives a wrong normalization.

Edmond Laguerre | 1834–1886

Edmond Nicolas Laguerre, a French artillery officer turned mathematician at the École Polytechnique, introduced the polynomials now bearing his name in his 1879 study of the integral ${\int }_x^{\infty }{e}^{-t}/tdt$. He was a master of analytic geometry and continued fractions. The associated (generalized) Laguerre polynomials $L_q^p$ were systematized later (and by Sonine/Sonin in Russia, hence "Sonine polynomials"); their starring role in quantum mechanics came only in 1926 when Schrödinger's wave equation for hydrogen demanded exactly these terminating polynomial solutions. Laguerre never saw the atomic application, dying four decades before the quantum theory that immortalized his name in every chemistry classroom.

Building Intuition

Think of the normalization constant as the answer to a probability bookkeeping question: "How much amplitude must I dial in so that the electron is somewhere with total probability one?" Three competing scales fight inside Eq. (5.61):

• The $(2/na{)}^3$ factor is dimensional — it turns the dimensionless $\rho$-integral back into a physical $r$-integral with units of (length)${}^{-3}$, so that $|R{|}^2{r}^{2}dr$ is dimensionless.
• The $(n-\ell -1)!/(n+\ell )!$ ratio is the "shape penalty": higher-order Laguerre polynomials oscillate more and spread the amplitude over many radial lobes, so the per-lobe height shrinks.
• The $1/2n$ factor reflects how spread out high-$n$ orbits are. The electron in a Rydberg state ($n$ large) lives far from the nucleus over a huge volume, so the wavefunction must be small everywhere — and $2n$ in the denominator under the square root enforces that.

The deep point: all the "physics" of normalization is encoded in the single number $2q+p+1=2n$ from the recurrence relation. The factorials are just the orthogonality baseline.

Worked Examples

Ground State $n=1,\ell=0$

State $n=2,\ell=1$

Common Misconceptions

Convention-Dependence of the Polynomial

The numerical coefficients in Eq. (5.61) are valid only for the Rodrigues convention [1]. Some texts (and Mathematica's LaguerreL[q,p,x]) use a different normalization in which $L_q^p$ differs by a factor of $(q+p)!$ or $q!$. If you import a polynomial from one convention into Griffiths' normalization constant, your $|R{|}^2$ will not integrate to 1. Always self-check by verifying [5] numerically for $n=1$.

Using $r\,dr$ Instead of $r^2\,dr$

Suppose a student mistakenly normalizes with the two-dimensional measure ${\int }_0^{\infty }|R{|}^{2}rdr$ (perhaps confusing it with a 2D radial problem). For the ground state this gives

Connections

Equation (5.61) sits at a crossroads of several theories:

• Associated Laguerre Polynomials supply the polynomial factor and the master integral [7]; the same polynomials appear in the 3D isotropic harmonic oscillator.
• Hydrogen Atom Schrödinger Equation produces the radial ODE [3]; the energy quantization $E_n\propto -1/n^2$ emerges from the termination condition that makes $v(\rho )$ a polynomial.
• Spherical Harmonics carry the angular dependence; the full eigenfunction is ${\psi }_{n\ell m}=R_{n\ell }(r)Y_{\ell }^m(\theta ,\varphi )$, with each factor independently normalized.
• Rodrigues' Formula is the unifying tool for evaluating both the orthogonality [11] and the master integral [7] via integration by parts.
• SO(4) Symmetry of Hydrogen explains the otherwise-accidental $\ell$-degeneracy that makes the energy depend only on $n$, the same $n$ appearing as $2n$ in the normalization.

Real-World Applications

• Atomic spectroscopy: Transition probabilities (oscillator strengths, Einstein $A$ coefficients) require dipole matrix elements $\int R_{n^{\prime }{\ell }^{\prime }}r{R}_{n\ell }{r}^{2}dr$, which are tractable only because the $R_{n\ell }$ are correctly normalized Laguerre expressions.
• Quantum chemistry basis sets: Slater-type and hydrogenic orbitals built from Eq. (5.61) seed self-consistent-field (Hartree–Fock) calculations of molecular structure.
• Photoionization cross-sections: Bound-to-continuum integrals in X-ray absorption rely on the normalized bound radial functions.
• Rydberg-atom physics: High-$n$ states used in quantum computing and precision metrology have radial extent scaling as $n^{2}a$, directly read off the $(2/na{)}^3$ scale factor.
• Astrophysical line formation: Hydrogen Balmer/Lyman series intensities in stellar atmospheres depend on level populations and transition strengths computed from these wavefunctions.
• Antihydrogen and exotic atoms: Muonic and pionic atoms reuse the identical algebra with rescaled Bohr radius $a\to a(m_e/\mu )$.
• Semiconductor donor states: Shallow hydrogenic impurity levels in silicon are modeled with effective-mass hydrogenic wavefunctions of the same form.

Historical Development

• 1879 — Edmond Laguerre introduces the basic Laguerre polynomials while studying the exponential-integral function.
• 1880s — Nikolai Sonine generalizes them to the associated/generalized polynomials $L_q^p$ (hence "Sonine polynomials" in older Russian literature).
• 1885 — Rodrigues-type and orthogonality properties consolidated within the emerging theory of classical orthogonal polynomials.
• 1926 (Jan) — Erwin Schrödinger publishes the wave equation and solves hydrogen, recognizing the radial solutions as associated Laguerre polynomials in his first Annalen der Physik communication.
• 1926 — Wolfgang Pauli derives the hydrogen spectrum independently via matrix mechanics and the Runge–Lenz vector, foreshadowing the SO(4) explanation of degeneracy.
• 1929 — Condon and Shortley begin systematizing phase conventions (the source of the leading minus sign debate) in their atomic-spectra work, culminating in their 1935 treatise.
• 1935 — The Theory of Atomic Spectra (Condon & Shortley) fixes standard normalization and phase conventions used by physicists thereafter.
• 1957 — Bethe and Salpeter's Quantum Mechanics of One- and Two-Electron Atoms tabulates the radial functions and matrix-element integrals in the modern reference form.
• 1995/2005 — Griffiths' Introduction to Quantum Mechanics assigns the now-canonical "Equation (5.61)" numbering used in countless undergraduate courses.

Subtleties & Limits

Several finer points reward reflection:

• Convention sensitivity means there is no single "right" Eq. (5.61) across textbooks — the factor structure depends on the chosen Laguerre normalization, the phase, and even whether reduced mass is used. A definitive cross-text reference must always state its convention explicitly.
• The continuum ($E>0$) states are not covered by [4]; for scattering states the Laguerre polynomials become non-terminating confluent hypergeometric functions and require delta-function ("energy") normalization instead of [5]. The smooth matching of bound to continuum normalization remains pedagogically delicate.
• Relativistic corrections (Dirac equation) replace the simple Laguerre structure with a coupled pair of radial functions whose normalization involves Laguerre polynomials of non-integer order — a genuine generalization where [7] no longer applies verbatim.
• Finite nuclear size and QED self-energy shift the very levels whose wavefunctions [4] describes; whether the nonrelativistic normalization remains a useful zeroth-order basis for high-precision (Lamb-shift) calculations is a practical question of perturbation-theory convergence.

See Also