Laguerre and Associated Laguerre Polynomials

Pose a physics problem on the half-line $[0,\infty)$ — a radial coordinate, a reaction rate, a lifetime, an energy that runs from zero to infinity — and one weight function keeps reappearing: the exponential $e^{-x}$ . It is the only elementary weight that decays fast enough to tame every power of $x$ at infinity while leaving the origin untouched, and the orthogonal polynomials it selects are the Laguerre polynomials $L_n(x)$ . They are what the naive basis $1,x,x^2,\dots$ becomes after it is orthogonalized against $e^{-x}$ : a clean, complete family in $L^2([0,\infty),e^{-x}\,\mathrm{d}x)$ , each member the polynomial eigenfunction of a singular second-order operator whose spectrum is forced to be the non-negative integers. Insert one extra dial — replace the weight by $x^{\alpha}e^{-x}$ — and the family generalizes to the associated (or generalized) Laguerre polynomials $L_n^{\alpha}(x)$ , with $L_n^0=L_n$ . That single parameter $\alpha$ is not cosmetic: it is exactly the factor a radial quantum problem needs to swallow its centrifugal barrier. This is why $L_n^{\alpha}$ , not $L_n$ , sits inside the hydrogen atom, the three-dimensional isotropic oscillator, Morse-type potentials, and the mode functions of a laser cavity. Structurally the family is the textbook classical-orthogonal-polynomial system: a Rodrigues formula, a compact generating function, a three-term recurrence, mixed derivative relations, and a self-adjoint Sturm–Liouville ODE — five faces of one object, each derivable from any other. Underneath every normalization, moment, and special value sits the Gamma function, which supplies the norm $\Gamma(n+\alpha+1)/n!$ . The polynomials therefore stand at the crossroads of analysis, atomic physics, and numerical integration, where the zeros of $L_n$ become the optimal nodes of Gauss–Laguerre quadrature.

Key takeaways

The Laguerre polynomials $L_n(x)$ are the orthogonal polynomials on $[0,\infty)$ for the exponential weight $e^{-x}$ ; the associated family $L_n^{\alpha}(x)$ uses the weight $x^{\alpha}e^{-x}$ with $\alpha>-1$ .
One object, five faces: a Sturm–Liouville ODE, a Rodrigues formula, an explicit finite sum, a generating function, and a three-term recurrence — each recoverable from any other.
The parameter $\alpha$ is the dial that absorbs the centrifugal term, so $L_n^{\alpha}$ governs the hydrogen atom ( $\alpha=2\ell+1$ ), the 3D isotropic oscillator ( $\alpha=\ell+\tfrac12$ ), and every exponential-weight radial problem.
The degree is physical: $L_n^{\alpha}$ has exactly $n$ real positive zeros, and in hydrogen that count is the number of radial nodes $n_r=n-\ell-1$ , so $n=n_r+\ell+1$ .
The orthogonality norm $\int_0^\infty x^\alpha e^{-x}(L_n^\alpha)^2\,\mathrm{d}x=\Gamma(n+\alpha+1)/n!$ is a pure Gamma-function value, and the zeros of $L_n$ drive Gauss–Laguerre quadrature, exact to degree $2n-1$ .
Key idea: Laguerre polynomials are the orthogonal polynomials for the exponential weight on the half-line; the parameter $\alpha$ in $x^{\alpha}e^{-x}$ is precisely the dial that fits radial quantum problems, so $L_n^{\alpha}$ is the universal radial building block.

Where this lives today

The associated Laguerre polynomials are live in three active fields. In photonics, the Laguerre–Gaussian modes $u_{p\ell}(r,\phi)\propto (r/w)^{|\ell|}L_p^{|\ell|}(2r^2/w^2)\,e^{-r^2/w^2}e^{i\ell\phi}$ are the paraxial eigenmodes of a cylindrically symmetric cavity; the azimuthal index $\ell$ makes each photon carry orbital angular momentum $\ell\hbar$ , and the Laguerre index $p$ counts the radial rings. These optical vortices drive optical tweezers that spin microscopic particles and underpin high-dimensional quantum key distribution and OAM-multiplexed communication, where one beam encodes many bits by mode. In quantum optics, the Wigner quasiprobability of a Fock state $|n\rangle$ is $W_n(\alpha)=\frac{2}{\pi}(-1)^n e^{-2|\alpha|^2}L_n(4|\alpha|^2)$ : the polynomial sculpts the concentric rings and the negativities that are the fingerprint of nonclassical light, and cavity- and circuit-QED experiments reconstruct exactly these Laguerre-patterned Wigner functions by tomography. And in computational physics, Gauss–Laguerre quadrature remains the default for exponential-tail integrals — Fermi–Dirac and Bose–Einstein integrals in astrophysical equations of state, partition-function tails, Laplace-transform inversion — because its nodes deliver spectral accuracy for a handful of evaluations. An open frontier is the instability of very high-degree zeros, where asymptotic node generators built from the Mehler–Heine formula now replace naive root-finding.

The physical picture

Picture the weight $e^{-x}$ as a soft wall that leans everything toward the origin. On $[0,\infty)$ you cannot use the monomials $1,x,x^2,\dots$ as an orthogonal basis, because each of them grows without bound and their overlap integrals diverge. The weight fixes this: $\int_0^\infty x^m x^n e^{-x}\,\mathrm{d}x=(m+n)!$ is finite, so Gram–Schmidt applied to the monomials under this inner product terminates cleanly and produces a family of polynomials, one per degree. $L_n$ is simply what $x^n$ becomes after this orthogonalization, rescaled to the standard normalization $L_n(0)=1$ . Because each $L_n$ is genuinely orthogonal to every lower degree, it must change sign often — a degree- $n$ orthogonal polynomial has exactly $n$ simple real zeros, all inside $(0,\infty)$ , and that node-counting fact is the seed of everything physical that follows.

Figure

The first Laguerre polynomials Lₙ(x) on the half-line: L₁ crosses zero once, L₂ twice, L₃ three times — a degree-n orthogonal polynomial has exactly n real positive zeros, and those zeros become the radial nodes of the hydrogen wavefunction. L₂ vanishes at x = 2 − √2 ≈ 0.586 and x = 2 + √2 ≈ 3.414

The parameter $\alpha$ tunes how hard the wall pushes near the origin. The factor $x^{\alpha}$ in the weight $x^{\alpha}e^{-x}$ suppresses functions near $x=0$ more strongly as $\alpha$ grows — the peak of the weight sits at $x=\alpha$ — so raising $\alpha$ literally moves probability away from the origin. In radial quantum mechanics that is exactly the centrifugal barrier: a particle with angular momentum $\ell$ is kept away from the nucleus, and when the Coulomb radial equation is stripped to its bare form the leftover ODE is the associated Laguerre equation with $\alpha=2\ell+1$ . The degree of the surviving polynomial is the number of interior radial nodes, so the principal quantum number reads literally as "radial nodes plus angular index plus one". The zeros in the figure above are not an abstraction — they are the shells where the electron's radial probability density touches zero. The same zeros do a second job that has nothing to do with quantum mechanics. Because $L_n$ is orthogonal to all lower-degree polynomials against $e^{-x}$ , its $n$ zeros are the uniquely optimal sampling points for approximating $\int_0^\infty e^{-x}f(x)\,\mathrm{d}x$ : an $n$ -point rule placed there integrates every polynomial up to degree $2n-1$ with zero error. This is Gauss–Laguerre quadrature, and it is the workhorse for exponential-tail integrals across statistical mechanics and computational physics.

Figure

Gauss–Laguerre quadrature: the two zeros of L₂(x), at x₁ = 0.5858 and x₂ = 3.4142, are the optimal nodes for the weight e⁻ˣ on the half-line. With weights w₁ = 0.8536 and w₂ = 0.1464 the 2-point rule integrates every polynomial up to degree 2n − 1 = 3 exactly — for instance ∫₀^∞ x³e⁻ˣ dx = 3! = 6 is reproduced to the last digit

The core equations

Everything below is a different projection of the same eigenvalue problem. The natural starting point is the differential equation itself, because it is where the integer $n$ and the parameter $\alpha$ first appear and where the demand for a normalizable solution does the quantizing.

\mathcal{L}\,y\equiv x\,y''+(\alpha+1-x)\,y'+n\,y=0 \qquad (1)

where: $y$ — the solution $L_n^{\alpha}(x)$ ; $x$ — independent variable on $[0,\infty)$ ; $\alpha$ — weight parameter, $\alpha>-1$ ; $n$ — non-negative integer index, equal to the polynomial degree. Reading it: this is a linear second-order ODE with a regular singular point at $x=0$ and an irregular one at $x=\infty$ ; only for integer $n$ does one of its two solutions stay polynomial, and that polynomial is $L_n^{\alpha}$ . Multiplying (1) by the integrating factor $x^{\alpha}e^{-x}$ turns it into manifestly self-adjoint (Sturm–Liouville) form, $\frac{\mathrm{d}}{\mathrm{d}x}\!\big(x^{\alpha+1}e^{-x}y'\big)+n\,x^{\alpha}e^{-x}y=0$ , which exposes the weight $w(x)=x^{\alpha}e^{-x}$ and guarantees orthogonality. The next equation packages the polynomial in closed form: the Rodrigues formula builds $L_n^{\alpha}$ by differentiating the weight $n$ times.

L_n^{\alpha}(x)=\frac{x^{-\alpha}e^{x}}{n!}\,\frac{\mathrm{d}^n}{\mathrm{d}x^n}\!\left(x^{n+\alpha}e^{-x}\right) \qquad (2)

where: $n!$ — normalizer fixing the leading coefficient to $(-1)^n/n!$ ; $x^{-\alpha}e^{x}$ — prefactor that restores polynomial form after the differentiation. Reading it: the operator $\mathrm{d}^n/\mathrm{d}x^n$ acting on $x^{n+\alpha}e^{-x}$ lowers the power of $x$ and reproduces the weight, so dividing back out by $w$ leaves a degree- $n$ polynomial. Expanding the derivative with the Leibniz rule gives the explicit coefficients.

L_n^{\alpha}(x)=\sum_{k=0}^{n}(-1)^k\binom{n+\alpha}{n-k}\frac{x^k}{k!} \qquad (3)

where: $\binom{n+\alpha}{n-k}=\dfrac{\Gamma(n+\alpha+1)}{\Gamma(k+\alpha+1)\,(n-k)!}$ — the generalized binomial coefficient; $k$ — summation index and power of $x$ . Reading it: every coefficient is a ratio of Gamma functions, so each is an exactly computable rational (for integer $\alpha$ ) or Gamma value; the alternating signs are what make the polynomial oscillate and cross zero $n$ times. The whole family can be squeezed into a single function of two variables. Expanding it in powers of $t$ generates every $L_n^{\alpha}$ at once, and differentiating it is the fastest route to the recurrences.

G(x,t)=\frac{1}{(1-t)^{\alpha+1}}\exp\!\left(-\frac{xt}{1-t}\right)=\sum_{n=0}^{\infty}L_n^{\alpha}(x)\,t^n,\qquad |t|<1 \qquad (4)

where: $t$ — bookkeeping (expansion) variable, $|t|<1$ ; $G$ — generating function packaging all $L_n^{\alpha}$ . Reading it: the coefficient of $t^n$ is exactly $L_n^{\alpha}(x)$ , so any operation you can perform on $G$ (differentiate in $t$ , in $x$ , multiply by $x$ ) becomes an identity among the polynomials once you match powers of $t$ . Matching powers after $\partial_t G$ yields the three-term recurrence — the numerically stable engine for generating the polynomials to any degree.

(n+1)\,L_{n+1}^{\alpha}(x)=(2n+\alpha+1-x)\,L_n^{\alpha}(x)-(n+\alpha)\,L_{n-1}^{\alpha}(x) \qquad (5)

where: $L_{n\pm1}^{\alpha}$ — neighbouring degrees at fixed $\alpha$ . Reading it: seeded by $L_0^{\alpha}=1$ and $L_1^{\alpha}=\alpha+1-x$ , this builds the entire family with only multiplications and additions; the appearance of $x$ (multiplication by the coordinate) is the general Favard/three-term structure shared by all orthogonal polynomials. Two companion derivative relations, $x\frac{\mathrm{d}}{\mathrm{d}x}L_n^{\alpha}=n\,L_n^{\alpha}-(n+\alpha)L_{n-1}^{\alpha}$ and $\frac{\mathrm{d}}{\mathrm{d}x}L_n^{\alpha}=-L_{n-1}^{\alpha+1}$ , connect the degree and parameter ladders. Finally the property that justifies the name — orthogonality, with the Gamma-function norm that fixes every physical normalization.

\int_0^{\infty} x^{\alpha}e^{-x}\,L_m^{\alpha}(x)\,L_n^{\alpha}(x)\,\mathrm{d}x=\frac{\Gamma(n+\alpha+1)}{n!}\,\delta_{mn} \qquad (6)

where: $w(x)=x^{\alpha}e^{-x}$ — weight on $[0,\infty)$ , requiring $\alpha>-1$ for convergence; $\delta_{mn}$ — Kronecker delta; $\Gamma$ — Gamma function supplying the norm. Reading it: distinct degrees are orthogonal against this fixed weight, and the diagonal value $\Gamma(n+\alpha+1)/n!$ is the single number that turns any expansion coefficient or wavefunction normalization into a Gamma ratio.

Derivation AOrthogonality and the Gamma-function norm

The goal is to prove both halves of equation (6) from first principles: that polynomials of different degree are orthogonal against $w(x)=x^{\alpha}e^{-x}$ , and that the diagonal norm is exactly $\Gamma(n+\alpha+1)/n!$ . The off-diagonal statement follows from the self-adjointness of the Sturm–Liouville operator, while the diagonal value comes from feeding the Rodrigues formula into the integral and integrating by parts $n$ times until only a bare Gamma integral survives. Nothing here is assumed about the polynomials except that they solve equation (1) and are normalized so their leading coefficient is $(-1)^n/n!$ .

Step 1. Write equation (1) for degree $n$ and for degree $m$ in Sturm–Liouville form, $(x^{\alpha+1}e^{-x}L_n^{\alpha}{}')'=-n\,x^{\alpha}e^{-x}L_n^{\alpha}$ . Multiply the $n$ -equation by $L_m^{\alpha}$ , the $m$ -equation by $L_n^{\alpha}$ , subtract, and integrate over $[0,\infty)$ ; the left side collapses to a total derivative, giving $(n-m)\int_0^\infty w\,L_m^{\alpha}L_n^{\alpha}\,\mathrm{d}x=\big[x^{\alpha+1}e^{-x}(L_m^{\alpha}{}'L_n^{\alpha}-L_m^{\alpha}L_n^{\alpha}{}')\big]_0^\infty$ . Why: subtracting the two eigenvalue equations is the standard Sturm–Liouville trick that converts a product of solutions into a boundary term.

Step 2. Evaluate the boundary term. At infinity $e^{-x}$ kills it; at the origin $x^{\alpha+1}\to0$ because $\alpha>-1$ . Both ends vanish, so $(n-m)\int_0^\infty w\,L_m^{\alpha}L_n^{\alpha}\,\mathrm{d}x=0$ , and for $m\neq n$ the integral itself must be zero — that is orthogonality. Why: the condition $\alpha>-1$ is not decoration; it is exactly what makes the weight integrable at the origin and the boundary term die, which is why the whole family is only defined for $\alpha>-1$ .

Step 3. For the diagonal, put one factor of $L_n^{\alpha}$ in Rodrigues form (2). The prefactor $x^{-\alpha}e^{x}$ cancels the weight, leaving $\int_0^\infty w\,(L_n^{\alpha})^2\,\mathrm{d}x=\frac{1}{n!}\int_0^\infty L_n^{\alpha}\,\frac{\mathrm{d}^n}{\mathrm{d}x^n}\!\big(x^{n+\alpha}e^{-x}\big)\,\mathrm{d}x$ . Integrate by parts $n$ times; every boundary term carries a positive power of $x$ times $e^{-x}$ and vanishes, flipping the derivatives onto $L_n^{\alpha}$ with a sign $(-1)^n$ . Why: repeated integration by parts is the engine that trades the awkward $n$ -fold derivative of the weight for a single derivative of the polynomial, which we know explicitly.

Step 4. Only the $n$ -th derivative of $L_n^{\alpha}$ survives, and since its leading term is $(-1)^n x^n/n!$ that derivative is the constant $(-1)^n$ . The integral reduces to $\frac{(-1)^n}{n!}(-1)^n\int_0^\infty x^{n+\alpha}e^{-x}\,\mathrm{d}x=\frac{1}{n!}\Gamma(n+\alpha+1)$ , which is the boxed norm. Why: the last integral is the definition of the Gamma function, so the norm of a Laguerre polynomial is never anything more exotic than a Gamma value — this is the hook that ties the whole family to the Gamma-function article. *Verified: at $n=0$ the norm is $\Gamma(\alpha+1)/0!=\int_0^\infty x^\alpha e^{-x}\,\mathrm{d}x$ , consistent with $L_0^\alpha=1$ ; the code check reproduces $\Gamma(n+\alpha+1)/n!$ for several $(n,\alpha)$ .*

Go deeper · Derivation B · 2 minHydrogen radial wavefunction from the Laguerre equation

Derivation BHydrogen radial wavefunction from the Laguerre equation

The goal is to show that the bound-state radial equation of the Coulomb problem is the associated Laguerre equation in disguise, and to read the quantum numbers straight off the identification. The strategy is standard asymptotic stripping: the radial equation has a centrifugal power $r^\ell$ forced by the origin and an exponential decay $e^{-r/na_0}$ forced by infinity; pulling both factors out of $R(r)$ leaves a smooth function that can only be a polynomial, and matching that polynomial equation term by term to equation (1) pins $\alpha$ and the degree. This is the calculation that makes $L_n^{\alpha}$ the marquee radial building block of atomic physics.

Step 1. Start from the reduced radial Schrödinger equation for a Coulomb potential, $-\frac{\hbar^2}{2m}\frac{1}{r^2}(r^2R')'+\big[\frac{\hbar^2\ell(\ell+1)}{2mr^2}-\frac{e^2}{4\pi\varepsilon_0 r}\big]R=ER$ , with $E<0$ for a bound state and $m$ the reduced mass. Why: this is the exact radial equation after separating the angular part; the two bracketed terms are the centrifugal barrier and the Coulomb attraction that shape the polynomial.

Step 2. Introduce the dimensionless variable $\rho=2r/(na_0)$ with $a_0=4\pi\varepsilon_0\hbar^2/me^2$ and the energy quantization $E=-\frac{me^4}{2(4\pi\varepsilon_0)^2\hbar^2 n^2}$ , then substitute the ansatz $R(r)=\rho^{\ell}e^{-\rho/2}u(\rho)$ . Why: $\rho^\ell$ captures the small- $r$ behaviour set by the centrifugal term and $e^{-\rho/2}$ captures the large- $r$ decay set by the negative energy, so whatever is left in $u$ is neither singular nor exponentially large.

Step 3. Carrying the ansatz through and dividing out the common factors gives $\rho\,u''+(2\ell+2-\rho)\,u'+(n-\ell-1)\,u=0$ . Compare with equation (1): this is the associated Laguerre ODE with $\alpha=2\ell+1$ and integer index $n-\ell-1$ . Polynomial (hence normalizable) solutions exist only when $n-\ell-1\in\{0,1,2,\dots\}$ , so $n\ge\ell+1$ . Why: the requirement that the leftover series terminate is what quantizes $n$ into integers and forces $\ell<n$ — the same truncation argument that gave integer $n$ in the abstract ODE now delivers the entire hydrogen spectrum.

Step 4. Identify $u(\rho)=L_{n-\ell-1}^{2\ell+1}(\rho)$ , so the radial function is $R_{n\ell}(r)=N_{n\ell}\,\rho^{\ell}e^{-\rho/2}L_{n-\ell-1}^{2\ell+1}(\rho)$ with $\rho=2r/(na_0)$ . The degree $n-\ell-1$ equals the number of interior radial nodes $n_r$ , and $N_{n\ell}$ is fixed by the norm (6) together with the $r^2\,\mathrm{d}r$ measure, giving $N_{n\ell}\propto\sqrt{(n-\ell-1)!/[2n\,(n+\ell)!]}\,(2/na_0)^{3/2}$ . Why: the abstract orthogonality norm becomes the physical wavefunction normalization, and the polynomial degree becomes a directly observable node count — the bridge from pure special-function theory to spectroscopy. *Verified: the ground state $n=1,\ell=0$ gives $L_0^1=1$ and $R_{10}\propto e^{-r/a_0}$ , the standard nodeless 1s orbital; the 3s state uses $L_2^1(\rho)=3-3\rho+\rho^2/2$ , whose zeros $\rho=3\pm\sqrt3$ map to radial nodes $r=1.902\,a_0$ and $7.098\,a_0$ — exactly two, as $n_r=n-\ell-1=2$ demands.*

Go deeper · Derivation C · 2 minGauss–Laguerre quadrature from the zeros

Derivation CGauss–Laguerre quadrature from the zeros

The goal is to prove that placing $n$ sample points at the zeros of $L_n$ produces an integration rule for $\int_0^\infty e^{-x}f(x)\,\mathrm{d}x$ that is exact for every polynomial up to degree $2n-1$ — twice the degree an ordinary $n$ -point interpolatory rule would manage — and to exhibit the closed-form weights. The argument rests entirely on orthogonality: the trick is polynomial division, writing any high-degree integrand as a multiple of $L_n$ plus a remainder, and letting orthogonality annihilate the $L_n$ part. This is the theoretical content behind the second figure and the reason the exponential-tail integrals of statistical mechanics are so cheap to evaluate.

Step 1. Take the nodes to be the $n$ real, simple, positive zeros $x_1,\dots,x_n$ of $L_n(x)$ (their reality and simplicity follow from orthogonality on $[0,\infty)$ ). Define the weights by the interpolatory condition that the rule be exact for all polynomials of degree $<n$ , i.e. $w_i=\int_0^\infty e^{-x}\ell_i(x)\,\mathrm{d}x$ with $\ell_i$ the Lagrange basis polynomial for node $i$ . Why: fixing the weights this way guarantees at least degree- $(n-1)$ exactness by construction, so the only question left is how much further orthogonality pushes it.

Step 2. Let $p(x)$ be any polynomial with $\deg p\le 2n-1$ and divide by $L_n$ : $p=q\,L_n+r$ with $\deg q\le n-1$ and $\deg r\le n-1$ . Integrating against the weight, $\int_0^\infty e^{-x}p\,\mathrm{d}x=\int_0^\infty e^{-x}q\,L_n\,\mathrm{d}x+\int_0^\infty e^{-x}r\,\mathrm{d}x$ ; the first integral is zero because $L_n$ is orthogonal to every polynomial of degree $\le n-1$ , including $q$ . Why: this single application of orthogonality is the whole miracle — it silently removes the top half of the integrand's degree.

Step 3. At the nodes $L_n(x_i)=0$ , so $p(x_i)=r(x_i)$ , and the rule reproduces $\int_0^\infty e^{-x}r\,\mathrm{d}x=\sum_i w_i r(x_i)=\sum_i w_i p(x_i)$ exactly, since $r$ has degree $<n$ . Hence the rule integrates all $p$ with $\deg p\le 2n-1$ exactly. Using the derivative recurrence $x_iL_n'(x_i)=-(n+1)L_{n+1}(x_i)$ collapses the Christoffel weights to the closed form $w_i=x_i/[(n+1)^2 L_{n+1}(x_i)^2]>0$ . Why: positivity of every weight is what makes the rule numerically stable and the closed form is what makes it practical — no matrix inversions, just evaluations of the next polynomial at the nodes. *Verified: $n=1$ gives one node $x_1=1$ (root of $L_1=1-x$ ) with $w_1=1$ , reproducing $\int_0^\infty e^{-x}\,\mathrm{d}x=1$ ; $n=2$ gives $x_{1,2}=2\mp\sqrt2$ with $w_{1,2}=0.8536,\,0.1464$ , and $w_1x_1^3+w_2x_2^3=6.000=\int_0^\infty x^3e^{-x}\,\mathrm{d}x=3!$ , exact because $\deg 3\le 2n-1$ .*

Limiting cases

Special values of $n$ , $\alpha$ , and $x$ , together with the large- $n$ asymptotics, connect $L_n^{\alpha}$ to simpler objects and expose its physical content.

Case	Result	Meaning
$\alpha=0$	$L_n^{0}(x)=L_n(x)$	the associated family collapses to the ordinary Laguerre polynomials; the weight is the bare $e^{-x}$
$n=0$	$L_0^{\alpha}(x)=1$	the lowest member is constant for every $\alpha$ — the seed of the recurrence and the nodeless ground state
$x=0$	$L_n^{\alpha}(0)=\binom{n+\alpha}{n}=\dfrac{\Gamma(n+\alpha+1)}{n!\,\Gamma(\alpha+1)}$	the standard special value that pins normalization and matches the confluent-hypergeometric scale
large $n$ , fixed $x$	$L_n^{\alpha}(x)\sim \dfrac{n^{\alpha/2-1/4}}{\sqrt{\pi}}\,\dfrac{e^{x/2}}{x^{\alpha/2+1/4}}\cos\!\big(2\sqrt{nx}-\tfrac{\pi}{2}\alpha-\tfrac{\pi}{4}\big)$	Mehler–Heine asymptotics: oscillation with a $\sqrt{nx}$ phase — the bridge to Bessel behaviour near the origin
confluent hypergeometric	$L_n^{\alpha}(x)=\binom{n+\alpha}{n}\,{}_1F_1(-n;\alpha+1;x)$	$L_n^{\alpha}$ is a terminating Kummer function, embedding it in the broader hypergeometric hierarchy

Common trapTwo conventions and one index that trap everyone

There are two live normalization conventions, and mixing them silently corrupts every coefficient. The modern physics standard (Arfken–Weber, Stone–Goldbart, Abramowitz–Stegun) uses the polynomials with $L_n^{\alpha}(0)=\binom{n+\alpha}{n}$ and leading coefficient $(-1)^n/n!$ — the ones on this page. Some older texts use a factorial-rescaled version differing by a factor $(n+\alpha)!$ , so their $L_n^{\alpha}$ is $(n+\alpha)!$ times ours; a formula copied across the divide will be off by that factor. The second, deadlier trap is the hydrogen index. The radial function carries $L_{n-\ell-1}^{2\ell+1}$ — degree $n-\ell-1$ , upper parameter $2\ell+1$ — not $L_n^{2\ell+1}$ or $L_{n+\ell}^{2\ell+1}$ . Getting the pair wrong shifts the node count $n_r=n-\ell-1$ and destroys the normalization, and because both indices are integers the error is easy to make and hard to spot. Always check that your polynomial has exactly $n_r$ zeros.

Worked examples

1 · Building L₂ and L₃ and checking the recurrence

Construct $L_2(x)$ and $L_3(x)$ (ordinary case, $\alpha=0$ ) from the explicit sum (3), then verify $L_3$ against the three-term recurrence (5).

Step 1. For $L_2$ use $\binom{2}{2-k}$ with $\alpha=0$ : the $k=0,1,2$ terms are $\binom{2}{2}=1$ , $-\binom{2}{1}x=-2x$ , and $+\binom{2}{0}x^2/2!=x^2/2$ . So $L_2(x)=\mathbf{1-2x+\tfrac12 x^2}$ . Its zeros solve $x^2-4x+2=0$ , giving $x=2\pm\sqrt2\approx0.586,\,3.414$ — two real positive roots, as required. Spot-check: $L_2(0)=1$ , $L_2(2)=1-4+2=-1$ , matching the figure.

Step 2. For $L_3$ : the terms $k=0,1,2,3$ are $1$ , $-3x$ , $+\binom{3}{1}x^2/2=\tfrac32 x^2$ , and $-\binom{3}{0}x^3/6=-\tfrac16 x^3$ , so $L_3(x)=\mathbf{1-3x+\tfrac32 x^2-\tfrac16 x^3}$ .

Step 3. Verify with the recurrence at $n=2$ , $\alpha=0$ : $(n+1)L_{n+1}=(2n+1-x)L_n-n\,L_{n-1}$ becomes $3L_3=(5-x)L_2-2L_1$ with $L_1=1-x$ . Expand $(5-x)(1-2x+\tfrac12x^2)=5-11x+\tfrac92x^2-\tfrac12x^3$ and subtract $2L_1=2-2x$ to get $3L_3=3-9x+\tfrac92x^2-\tfrac12x^3$ . Dividing by 3 gives $L_3=\mathbf{1-3x+\tfrac32x^2-\tfrac16x^3}$ , matching Step 2 exactly. The explicit sum and the recurrence agree.

2 · Two-point Gauss–Laguerre quadrature, exact to degree three

Use the zeros of $L_2$ to integrate $\int_0^\infty x^3 e^{-x}\,\mathrm{d}x$ and confirm the exact value $3!=6$ .

Step 1. Nodes are the zeros of $L_2$ : $x_1=2-\sqrt2=\mathbf{0.585786}$ and $x_2=2+\sqrt2=\mathbf{3.414214}$ . The weights from $w_i=x_i/[(n+1)^2 L_{n+1}(x_i)^2]$ with $n=2$ and $L_3$ from Example 1 evaluate to $w_1=\mathbf{0.853553}$ and $w_2=\mathbf{0.146447}$ (and $w_1+w_2=1=\int_0^\infty e^{-x}\,\mathrm{d}x$ , a useful sanity check).

Step 2. The integrand is $f(x)=x^3$ , of degree $3=2n-1$ , so the rule is exact. Evaluate: $x_1^3=0.201009$ , $x_2^3=39.79899$ .

Step 3. Sum: $w_1x_1^3+w_2x_2^3=0.853553\times0.201009+0.146447\times39.79899=0.171573+5.828427=\mathbf{6.000000}$ . This matches $\int_0^\infty x^3e^{-x}\,\mathrm{d}x=\Gamma(4)=3!=6$ to every digit carried. A two-point rule reproduces a cubic-moment integral exactly, because the nodes are Laguerre zeros.

3 · Reading node counts off the hydrogen radial function

For the $n=3$ shell, write the Laguerre factor for the 3s, 3p, and 3d states and state how many radial nodes each has.

Step 1. The radial function carries $L_{n-\ell-1}^{2\ell+1}(\rho)$ with $\rho=2r/(3a_0)$ . For 3s ( $\ell=0$ ): $L_2^{1}(\rho)=\mathbf{3-3\rho+\tfrac12\rho^2}$ , degree 2. For 3p ( $\ell=1$ ): $L_1^{3}(\rho)=\mathbf{4-\rho}$ , degree 1. For 3d ( $\ell=2$ ): $L_0^{5}=\mathbf{1}$ , degree 0.

Step 2. Interior radial nodes equal the polynomial degree $n_r=n-\ell-1$ : 3s has $\mathbf{2}$ nodes, 3p has $\mathbf{1}$ , 3d has $\mathbf{0}$ . Check the sum rule $n_r+\ell+1=n$ : $2+0+1=3$ , $1+1+1=3$ , $0+2+1=3$ — all give the principal quantum number 3, as they must.

Step 3. Locate the 3s nodes: $L_2^1(\rho)=0\Rightarrow\rho^2-6\rho+6=0\Rightarrow\rho=3\pm\sqrt3$ , and $r=\tfrac{3a_0}{2}\rho$ gives $r=\mathbf{1.90\,a_0}$ and $\mathbf{7.10\,a_0}$ . The single 3p node sits at $\rho=4$ , i.e. $r=\mathbf{6\,a_0}$ . The degree of the Laguerre factor is literally the number of shells where the electron's radial density vanishes.

Physics of the portrait

The Laguerre polynomials are named for a geometer who met them almost in passing, while chasing an entirely different integral — and who would have been astonished to learn that his little family sits at the heart of the hydrogen atom, the laser, and the quantum theory of light.

Origins · 1 minEdmond Laguerre and how this was discovered

Edmond Laguerre

1834–1886 · orthogonal polynomials and the exponential-weight problem

Edmond Nicolas Laguerre was born in Bar-le-Duc in 1834 and trained at the École Polytechnique, where — after a spell as an artillery officer cut short by poor health — he returned as a tutor and later professor, taking a chair at the Collège de France in 1885, the year before his death. He was primarily a geometer and analyst, remembered for Laguerre geometry (oriented lines and circles), the Laguerre transformations of projective geometry, deep work on the theory of algebraic equations, and Laguerre's still-used iterative method for polynomial roots. The polynomials arrived through a side door: in 1879, studying the exponential integral $\int_x^\infty t^{-1}e^{-t}\,\mathrm{d}t$ and its continued-fraction expansion, he was led to the polynomials orthogonal against $e^{-x}$ on the half-line and worked out their properties. Special cases had been touched earlier — Abel examined the underlying integral, and Chebyshev's general theory already contained the seed — but Laguerre gave the exponential-weight family its systematic treatment. The generalization to the weight $x^{\alpha}e^{-x}$ is due to Nikolay Sonine in 1880 (hence the alternative name Sonine polynomials). The physics came half a century later: when Schrödinger solved the hydrogen atom in 1926, the associated Laguerre polynomials turned out to be exactly the radial factors of the wavefunctions — nineteenth-century analysis waiting for its twentieth-century application.

History

1826

— Abel examines the exponential-weight integral that seeds the family

1859

— Chebyshev develops the general theory of orthogonal polynomials on a weight

1879

— Laguerre introduces Lₙ(x) while studying the exponential integral ∫ₓ^∞ t⁻¹e⁻ᵗ dt

1880

— Sonine generalizes to the associated polynomials Lₙ^α, the Sonine polynomials

1926

— Schrödinger solves the hydrogen atom and Lₙ^α appears in the radial wavefunction

Going further

Texts. Arfken and Weber, Mathematical Methods for Physicists, ch. 13 (Laguerre functions) · Stone and Goldbart, Mathematics for Physics, Sturm–Liouville and special-functions chapters · Hassani, Mathematical Physics, chapters on classical orthogonal polynomials · Byron and Fuller, Mathematics of Classical and Quantum Physics, special-functions chapters · Riley, Hobson and Bence, Mathematical Methods for Physics and Engineering · Abramowitz and Stegun, Handbook of Mathematical Functions, ch. 22 · Sakurai and Napolitano, Modern Quantum Mechanics (hydrogen radial solution).

Provenance. All formulas are standard and cross-checked against Arfken–Weber ch. 13, Stone–Goldbart, and Abramowitz–Stegun ch. 22. A verifier should: (1) confirm the norm $\Gamma(n+\alpha+1)/n!$ by the Rodrigues integration-by-parts argument; (2) recompute $L_0,L_1,L_2,L_3$ from the explicit sum and check the three-term recurrence (Example 1); (3) expand the generating function to order $t^2$ and match $L_0,L_1,L_2$ ; (4) verify the hydrogen reduction gives $\alpha=2\ell+1$ , degree $n-\ell-1$ , and that $L_0^1=1$ reproduces the 1s orbital; (5) verify the 2-point Gauss–Laguerre nodes $x=2\mp\sqrt2$ , weights $0.8536,0.1464$ , and the exact result $\int_0^\infty x^3e^{-x}\,\mathrm{d}x=6$ . The Wigner-function value $W_n=\frac{2}{\pi}(-1)^n e^{-2|\alpha|^2}L_n(4|\alpha|^2)$ is the standard Cahill–Glauber convention.

Prerequisites. The Gamma Function · Orthogonal Polynomials: General Theory · Generating Functions for Special Functions · Recurrence Relations for Special Functions. Leads to. Associated Legendre Functions · Spherical Bessel Functions · The Gamma Function.

Ask the tutor: Ask me to derive any single formula on this page from the generating function, or to build the full 3s, 3p, or 3d hydrogen radial wavefunction and count its nodes with you step by step.

Why the exponential weight, and not some other decaying function?

Because $e^{-x}$ is the unique elementary weight that makes the half-line behave like a self-contained spectral problem. It decays fast enough that every moment $\int_0^\infty x^n e^{-x}\,\mathrm{d}x=n!$ is finite, so all the monomials live in the space and Gram–Schmidt closes on a clean polynomial family. Physically the exponential is not arbitrary either: it is the large-distance behaviour forced by any binding potential (a negative-energy bound state must decay exponentially), so the natural radial coordinate always arrives pre-dressed with an $e^{-r}$ factor. Pull that factor out and the leftover is a polynomial in the Laguerre family. Other weights give other families — $e^{-x^2}$ on the whole line gives Hermite, $(1-x^2)^\alpha$ on $[-1,1]$ gives Jacobi — but the half-line-with-exponential-decay geometry is Laguerre's alone.

What exactly does the parameter α do physically?

It tunes how strongly the weight suppresses the origin. Since the weight is $x^\alpha e^{-x}$ and peaks at $x=\alpha$ , raising $\alpha$ pushes the bulk of the distribution outward, away from $x=0$ . In a radial quantum problem that is precisely the centrifugal barrier: a particle with angular momentum $\ell$ feels a repulsive $\ell(\ell+1)/r^2$ term that keeps it away from the centre, and when you reduce the radial equation to Laguerre form that barrier reappears as $\alpha=2\ell+1$ in the hydrogen atom (or $\alpha=\ell+\tfrac12$ in the isotropic oscillator, working in the variable $r^2$ ). So $\alpha$ is not a free mathematical decoration — it is the angular momentum of the physical state, encoded in the weight.

How can one generating function contain the whole infinite family?

The generating function $G(x,t)=(1-t)^{-\alpha-1}\exp[-xt/(1-t)]$ is a single analytic function of two variables whose Taylor coefficients in $t$ are the polynomials: the coefficient of $t^n$ is $L_n^\alpha(x)$ . This is enormously economical, because any manipulation of $G$ becomes an identity among all the polynomials at once. Differentiate $G$ in $t$ and match powers, and you get the three-term recurrence; differentiate in $x$ , and you get the derivative relations; set $x=0$ , and you read off the special values $L_n^\alpha(0)$ . In practice a fluent physicist rarely memorizes the individual formulas — they remember $G$ and regenerate whatever identity they need by differentiating and comparing coefficients of $t^n$ .

Why do the polynomial's zeros matter so much?

Because the zeros do double duty. First, a degree- $n$ orthogonal polynomial has exactly $n$ simple zeros inside the interval, and in a radial wavefunction those zeros are the physical radial nodes — the shells where the electron's probability density vanishes, which is why $n_r=n-\ell-1$ is a directly countable feature of an orbital. Second, those same zeros are the mathematically optimal sampling points for numerical integration: a quadrature rule placed at the zeros of $L_n$ integrates polynomials up to degree $2n-1$ exactly, double what any other node placement achieves. So the zeros are simultaneously where the physics goes quiet and where the numerics goes maximally accurate — the same list of numbers wearing two hats.

Is the Laguerre index the same as the principal quantum number n?

No, and conflating them is the most common hydrogen error. The principal quantum number $n$ labels the energy level and satisfies $n\ge\ell+1$ . The Laguerre polynomial in the radial function is $L_{n-\ell-1}^{2\ell+1}$ : its degree is $n-\ell-1$ (the radial node count), and its upper parameter is $2\ell+1$ . So the polynomial's degree is $n$ minus $\ell$ minus one, not $n$ itself, and the two integers labelling the polynomial are both smaller than or different from the principal quantum number. If you write $L_n^{2\ell+1}$ by mistake you get a polynomial of the wrong degree, the wrong number of nodes, and a normalization that no longer matches the $r^2\,\mathrm{d}r$ measure — a silent, wavefunction-corrupting bug.

How are Laguerre polynomials related to the confluent hypergeometric function?

They are a terminating special case of it. The associated Laguerre polynomial equals $L_n^\alpha(x)=\binom{n+\alpha}{n}\,{}_1F_1(-n;\alpha+1;x)$ , where ${}_1F_1$ is Kummer's confluent hypergeometric function. The key is the first argument $-n$ : when it is a negative integer the hypergeometric series truncates after $n+1$ terms and becomes a polynomial. This embedding is more than bookkeeping — it places Laguerre inside the same family as the Coulomb wavefunctions, the Bessel functions (via limits), and the Hermite polynomials, and it means any identity proved for ${}_1F_1$ specializes instantly to a Laguerre identity. It is also the cleanest way to see the Mehler–Heine asymptotics in the limiting-cases table, since the large- $n$ limit of ${}_1F_1$ is a Bessel function.