Laguerre polynomials

In mathematics, the Laguerre polynomials, named after Edmond Laguerre (1834 - 1886), are solutions of Laguerre's equation:

${\displaystyle xy''+(1-x)y'+ny=0}$

which is a second-order linear differential equation. This equation has nonsingular solutions only if n is a non-negative integer.

The associated Laguerre polynomials (alternatively, but rarely, named Sonin polynomials, after their inventor[1] Nikolay Yakovlevich Sonin) are solutions of

${\displaystyle xy''+(\alpha +1-x)y'+ny=0~.}$

The Laguerre polynomials are also used for Gaussian quadrature to numerically compute integrals of the form

${\displaystyle \int _{0}^{\infty }f(x)e^{-x}\,dx.}$

These polynomials, usually denoted L0L1, ..., are a polynomial sequence which may be defined by the Rodrigues formula,

${\displaystyle L_{n}(x)={\frac {e^{x}}{n!}}{\frac {d^{n}}{dx^{n}}}\left(e^{-x}x^{n}\right)={\frac {1}{n!}}\left({\frac {d}{dx}}-1\right)^{n}x^{n},}$

reducing to the closed form of a following section.

They are orthogonal polynomials with respect to an inner product

${\displaystyle \langle f,g\rangle =\int _{0}^{\infty }f(x)g(x)e^{-x}\,dx.}$

The sequence of Laguerre polynomials n! Ln is a Sheffer sequence,

${\displaystyle {\frac {d}{dx}}L_{n}=\left({\frac {d}{dx}}-1\right)L_{n-1}.}$

The Rook polynomials in combinatorics are more or less the same as Laguerre polynomials, up to elementary changes of variables. Further see the Tricomi–Carlitz polynomials.

The Laguerre polynomials arise in quantum mechanics, in the radial part of the solution of the Schrödinger equation for a one-electron atom. They also describe the static Wigner functions of oscillator systems in quantum mechanics in phase space. They further enter in the quantum mechanics of the Morse potential and of the 3D isotropic harmonic oscillator.

Physicists sometimes use a definition for the Laguerre polynomials which is larger by a factor of n! than the definition used here. (Likewise, some physicist may use somewhat different definitions of the so-called associated Laguerre polynomials.)

The first few polynomials

These are the first few Laguerre polynomials:

The first six Laguerre polynomials.

Recursive definition, closed form, and generating function

One can also define the Laguerre polynomials recursively, defining the first two polynomials as

${\displaystyle L_{0}(x)=1}$
${\displaystyle L_{1}(x)=1-x}$

and then using the following recurrence relation for any k ≥ 1:

${\displaystyle L_{k+1}(x)={\frac {(2k+1-x)L_{k}(x)-kL_{k-1}(x)}{k+1}}.}$

The closed form is

${\displaystyle L_{n}(x)=\sum _{k=0}^{n}{\binom {n}{k}}{\frac {(-1)^{k}}{k!}}x^{k}.}$

The generating function for them likewise follows,

${\displaystyle \sum _{n}^{\infty }t^{n}L_{n}(x)={\frac {1}{1-t}}e^{-{\frac {tx}{1-t}}}.}$

Generalized Laguerre polynomials

For arbitrary real α the polynomial solutions of the differential equation [2]

${\displaystyle x\,y''+(\alpha +1-x)\,y'+n\,y=0}$

are called generalized Laguerre polynomials, or associated Laguerre polynomials.

One can also define the generalized Laguerre polynomials recursively, defining the first two polynomials as

${\displaystyle L_{0}^{(\alpha )}(x)=1}$
${\displaystyle L_{1}^{(\alpha )}(x)=1+\alpha -x}$

and then using the following recurrence relation for any k ≥ 1:

${\displaystyle L_{k+1}^{(\alpha )}(x)={\frac {(2k+1+\alpha -x)L_{k}^{(\alpha )}(x)-(k+\alpha )L_{k-1}^{(\alpha )}(x)}{k+1}}.}$

The simple Laguerre polynomials are included in the associated polynomials, through α = 0,

${\displaystyle L_{n}^{(0)}(x)=L_{n}(x).}$

The Rodrigues formula for them is

{\displaystyle {\begin{aligned}L_{n}^{(\alpha )}(x)&={x^{-\alpha }e^{x} \over n!}{d^{n} \over dx^{n}}\left(e^{-x}x^{n+\alpha }\right)\\&=x^{-\alpha }{\frac {({\frac {d}{dx}}-1)^{n}}{n!}}x^{n+\alpha }.\end{aligned}}}

The generating function for them is

${\displaystyle \sum _{n}^{\infty }t^{n}L_{n}^{(\alpha )}(x)={\frac {1}{(1-t)^{\alpha +1}}}e^{-{\frac {tx}{1-t}}}.}$
The first few associated Laguerre polynomials

Explicit examples and properties of the associated Laguerre polynomials

${\displaystyle L_{n}^{(\alpha )}(x):={n+\alpha \choose n}M(-n,\alpha +1,x).}$
When Template:Mvar is an integer the function reduces to a polynomial of degree Template:Mvar. It has the alternative expression[4]
${\displaystyle L_{n}^{(\alpha )}(x)={\frac {(-1)^{n}}{n!}}U(-n,\alpha +1,x)}$
in terms of Kummer's function of the second kind.
• The closed form for these associated Laguerre polynomials of degree Template:Mvar is[5]
${\displaystyle L_{n}^{(\alpha )}(x)=\sum _{i=0}^{n}(-1)^{i}{n+\alpha \choose n-i}{\frac {x^{i}}{i!}}}$
derived by applying Leibniz's theorem for differentiation of a product to Rodrigues' formula.
• The first few generalized Laguerre polynomials are:
{\displaystyle {\begin{aligned}L_{0}^{(\alpha )}(x)&=1\\L_{1}^{(\alpha )}(x)&=-x+\alpha +1\\L_{2}^{(\alpha )}(x)&={\frac {x^{2}}{2}}-(\alpha +2)x+{\frac {(\alpha +2)(\alpha +1)}{2}}\\L_{3}^{(\alpha )}(x)&={\frac {-x^{3}}{6}}+{\frac {(\alpha +3)x^{2}}{2}}-{\frac {(\alpha +2)(\alpha +3)x}{2}}+{\frac {(\alpha +1)(\alpha +2)(\alpha +3)}{6}}\end{aligned}}}
${\displaystyle L_{n}^{(\alpha )}(0)={n+\alpha \choose n}\approx {\frac {n^{\alpha }}{\Gamma (\alpha +1)}};}$

{{#invoke:Category handler|main}}{{#invoke:Category handler|main}}[citation needed] }}

${\displaystyle L_{n}^{(\alpha )}(x)={\frac {n^{{\frac {\alpha }{2}}-{\frac {1}{4}}}}{\sqrt {\pi }}}{\frac {e^{\frac {x}{2}}}{x^{{\frac {\alpha }{2}}+{\frac {1}{4}}}}}\cos \left(2{\sqrt {nx}}-{\frac {\pi }{2}}\left(\alpha +{\frac {1}{2}}\right)\right)+O\left(n^{{\frac {\alpha }{2}}-{\frac {3}{4}}}\right),}$
${\displaystyle L_{n}^{(\alpha )}(-x)={\frac {(n+1)^{{\frac {\alpha }{2}}-{\frac {1}{4}}}}{2{\sqrt {\pi }}}}{\frac {e^{-{\frac {x}{2}}}}{x^{{\frac {\alpha }{2}}+{\frac {1}{4}}}}}e^{2{\sqrt {x(n+1)}}}\cdot \left(1+O\left({\frac {1}{\sqrt {n+1}}}\right)\right),}$
and summarizing by
${\displaystyle {\frac {L_{n}^{(\alpha )}\left({\frac {x}{n}}\right)}{n^{\alpha }}}\approx e^{\frac {x}{2n}}\cdot {\frac {J_{\alpha }\left(2{\sqrt {x}}\right)}{{\sqrt {x}}^{\alpha }}},}$
where ${\displaystyle J_{\alpha }}$ is the Bessel function.

As a contour integral

Given the generating function specified above, the polynomials may be expressed in terms of a contour integral

${\displaystyle L_{n}^{(\alpha )}(x)={\frac {1}{2\pi i}}\oint {\frac {e^{-{\frac {xt}{1-t}}}}{(1-t)^{\alpha +1}\,t^{n+1}}}\;dt,}$

where the contour circles the origin once in a counterclockwise direction.

Recurrence relations

The addition formula for Laguerre polynomials:[8]

${\displaystyle L_{n}^{(\alpha +\beta +1)}(x+y)=\sum _{i=0}^{n}L_{i}^{(\alpha )}(x)L_{n-i}^{(\beta )}(y)}$.

Laguerre's polynomials satisfy the recurrence relations

${\displaystyle L_{n}^{(\alpha )}(x)=\sum _{i=0}^{n}L_{n-i}^{(\alpha +i)}(y){\frac {(y-x)^{i}}{i!}},}$

in particular

${\displaystyle L_{n}^{(\alpha +1)}(x)=\sum _{i=0}^{n}L_{i}^{(\alpha )}(x)}$

and

${\displaystyle L_{n}^{(\alpha )}(x)=\sum _{i=0}^{n}{\alpha -\beta +n-i-1 \choose n-i}L_{i}^{(\beta )}(x),}$

or

${\displaystyle L_{n}^{(\alpha )}(x)=\sum _{i=0}^{n}{\alpha -\beta +n \choose n-i}L_{i}^{(\beta -i)}(x);}$

moreover

{\displaystyle {\begin{aligned}L_{n}^{(\alpha )}(x)-\sum _{j=0}^{\Delta -1}{n+\alpha \choose n-j}(-1)^{j}{\frac {x^{j}}{j!}}&=(-1)^{\Delta }{\frac {x^{\Delta }}{(\Delta -1)!}}\sum _{i=0}^{n-\Delta }{\frac {n+\alpha \choose n-\Delta -i}{(n-i){n \choose i}}}L_{i}^{(\alpha +\Delta )}(x)\\[6pt]&=(-1)^{\Delta }{\frac {x^{\Delta }}{(\Delta -1)!}}\sum _{i=0}^{n-\Delta }{\frac {n+\alpha -i-1 \choose n-\Delta -i}{(n-i){n \choose i}}}L_{i}^{(n+\alpha +\Delta -i)}(x)\end{aligned}}}

They can be used to derive the four 3-point-rules

combined they give this additional, useful recurrence relations

{\displaystyle {\begin{aligned}L_{n}^{(\alpha )}(x)&=\left(2+{\frac {\alpha -1-x}{n}}\right)L_{n-1}^{(\alpha )}(x)-\left(1+{\frac {\alpha -1}{n}}\right)L_{n-2}^{(\alpha )}(x)\\[10pt]&={\frac {\alpha +1-x}{n}}L_{n-1}^{(\alpha +1)}(x)-{\frac {x}{n}}L_{n-2}^{(\alpha +2)}(x)\end{aligned}}}

A somewhat curious identity, valid for integer i and Template:Mvar, is

${\displaystyle {\frac {(-x)^{i}}{i!}}L_{n}^{(i-n)}(x)={\frac {(-x)^{n}}{n!}}L_{i}^{(n-i)}(x);}$

it may be used to derive the partial fraction decomposition

{\displaystyle {\begin{aligned}{\frac {L_{n}^{(\alpha )}(x)}{n+\alpha \choose n}}&=1-\sum _{j=1}^{n}{\frac {x^{j}}{\alpha +j}}{\frac {L_{n-j}^{(j)}(x)}{(j-1)!}}\\&=1-\sum _{j=1}^{n}(-1)^{j}{\frac {j}{\alpha +j}}{n \choose j}L_{n}^{(-j)}(x)\\&=1-x\sum _{i=1}^{n}{\frac {L_{n-i}^{(-\alpha )}(x)L_{i-1}^{(\alpha +1)}(-x)}{\alpha +i}}\end{aligned}}}

Derivatives of generalized Laguerre polynomials

Differentiating the power series representation of a generalized Laguerre polynomial k times leads to

${\displaystyle {\frac {d^{k}}{dx^{k}}}L_{n}^{(\alpha )}(x)=(-1)^{k}L_{n-k}^{(\alpha +k)}(x).}$

This points to a special case (α = 0) of the formula above: for integer α = k the generalized polynomial may be written

${\displaystyle L_{n}^{(k)}(x)=(-1)^{k}{\frac {d^{k}L_{n+k}(x)}{dx^{k}}},}$

the shift by k sometimes causing confusion with the usual parenthesis notation for a derivative.

Moreover, this following equation holds

${\displaystyle {\frac {1}{k!}}{\frac {d^{k}}{dx^{k}}}x^{\alpha }L_{n}^{(\alpha )}(x)={n+\alpha \choose k}x^{\alpha -k}L_{n}^{(\alpha -k)}(x),}$

which generalizes with Cauchy's formula to

${\displaystyle L_{n}^{(\alpha ')}(x)=(\alpha '-\alpha ){\alpha '+n \choose \alpha '-\alpha }\int _{0}^{x}{\frac {t^{\alpha }(x-t)^{\alpha '-\alpha -1}}{x^{\alpha '}}}L_{n}^{(\alpha )}(t)\,dt.}$

The derivative with respect to the second variable Template:Mvar has the form,[9]

${\displaystyle {\frac {d}{d\alpha }}L_{n}^{(\alpha )}(x)=\sum _{i=0}^{n-1}{\frac {L_{i}^{(\alpha )}(x)}{n-i}}.}$

This is evident from the contour integral representation below.

The generalized associated Laguerre polynomials obey the differential equation

${\displaystyle xL_{n}^{(\alpha )\prime \prime }(x)+(\alpha +1-x)L_{n}^{(\alpha )\prime }(x)+nL_{n}^{(\alpha )}(x)=0,}$

which may be compared with the equation obeyed by the kth derivative of the ordinary Laguerre polynomial,

${\displaystyle xL_{n}^{(k)\prime \prime }(x)+(k+1-x)L_{n}^{(k)\prime }(x)+(n-k)L_{n}^{(k)}(x)=0,}$

where ${\displaystyle L_{n}^{(k)}(x)\equiv {\frac {d^{k}L_{n}(x)}{dx^{k}}}}$ for this equation only.

In Sturm–Liouville form the differential equation is

${\displaystyle -\left(x^{\alpha +1}e^{-x}\cdot L_{n}^{(\alpha )}(x)^{\prime }\right)^{\prime }=n\cdot x^{\alpha }e^{-x}\cdot L_{n}^{(\alpha )}(x),}$

which shows that is an eigenvector for the eigenvalue Template:Mvar.

Orthogonality

The associated Laguerre polynomials are orthogonal over [0, ∞) with respect to the measure with weighting function xα ex:[10]

${\displaystyle \int _{0}^{\infty }x^{\alpha }e^{-x}L_{n}^{(\alpha )}(x)L_{m}^{(\alpha )}(x)dx={\frac {\Gamma (n+\alpha +1)}{n!}}\delta _{n,m},}$

which follows from

${\displaystyle \int _{0}^{\infty }x^{\alpha '-1}e^{-x}L_{n}^{(\alpha )}(x)dx={\alpha -\alpha '+n \choose n}\Gamma (\alpha ').}$

If ${\displaystyle \Gamma (x,\alpha +1,1)}$ denoted the Gamma distribution then the orthogonality relation can be written as

${\displaystyle \int _{0}^{\infty }L_{n}^{(\alpha )}(x)L_{m}^{(\alpha )}(x)\Gamma (x,\alpha +1,1)dx={n+\alpha \choose n}\delta _{n,m},}$

The associated, symmetric kernel polynomial has the representations (Christoffel–Darboux formula){{ safesubst:#invoke:Unsubst||date=__DATE__ |\$B= {{#invoke:Category handler|main}}{{#invoke:Category handler|main}}[citation needed] }}

{\displaystyle {\begin{aligned}K_{n}^{(\alpha )}(x,y)&:={\frac {1}{\Gamma (\alpha +1)}}\sum _{i=0}^{n}{\frac {L_{i}^{(\alpha )}(x)L_{i}^{(\alpha )}(y)}{\alpha +i \choose i}}\\&{=}{\frac {1}{\Gamma (\alpha +1)}}{\frac {L_{n}^{(\alpha )}(x)L_{n+1}^{(\alpha )}(y)-L_{n+1}^{(\alpha )}(x)L_{n}^{(\alpha )}(y)}{{\frac {x-y}{n+1}}{n+\alpha \choose n}}}\\&{=}{\frac {1}{\Gamma (\alpha +1)}}\sum _{i=0}^{n}{\frac {x^{i}}{i!}}{\frac {L_{n-i}^{(\alpha +i)}(x)L_{n-i}^{(\alpha +i+1)}(y)}{{\alpha +n \choose n}{n \choose i}}};\end{aligned}}}

recursively

${\displaystyle K_{n}^{(\alpha )}(x,y)={\frac {y}{\alpha +1}}K_{n-1}^{(\alpha +1)}(x,y)+{\frac {1}{\Gamma (\alpha +1)}}{\frac {L_{n}^{(\alpha +1)}(x)L_{n}^{(\alpha )}(y)}{\alpha +n \choose n}}.}$

Moreover,

${\displaystyle y^{\alpha }e^{-y}K_{n}^{(\alpha )}(\cdot ,y)\to \delta (y-\cdot ),}$

in the associated L2[0, ∞)-space.

Turán's inequalities can be derived here, which is

${\displaystyle L_{n}^{(\alpha )}(x)^{2}-L_{n-1}^{(\alpha )}(x)L_{n+1}^{(\alpha )}(x)=\sum _{k=0}^{n-1}{\frac {\alpha +n-1 \choose n-k}{n{n \choose k}}}L_{k}^{(\alpha -1)}(x)^{2}>0.}$

The following integral is needed in the quantum mechanical treatment of the hydrogen atom,

${\displaystyle \int _{0}^{\infty }x^{\alpha +1}e^{-x}\left[L_{n}^{(\alpha )}(x)\right]^{2}dx={\frac {(n+\alpha )!}{n!}}(2n+\alpha +1).}$

Series expansions

Let a function have the (formal) series expansion

${\displaystyle f(x)=\sum _{i=0}^{\infty }f_{i}^{(\alpha )}L_{i}^{(\alpha )}(x).}$

Then

${\displaystyle f_{i}^{(\alpha )}=\int _{0}^{\infty }{\frac {L_{i}^{(\alpha )}(x)}{i+\alpha \choose i}}\cdot {\frac {x^{\alpha }e^{-x}}{\Gamma (\alpha +1)}}\cdot f(x)\,dx.}$

The series converges in the associated Hilbert space L2[0, ∞) if and only if

${\displaystyle \|f\|_{L^{2}}^{2}:=\int _{0}^{\infty }{\frac {x^{\alpha }e^{-x}}{\Gamma (\alpha +1)}}|f(x)|^{2}dx=\sum _{i=0}^{\infty }{i+\alpha \choose i}|f_{i}^{(\alpha )}|^{2}<\infty .}$

Further examples of expansions

Monomials are represented as

${\displaystyle {\frac {x^{n}}{n!}}=\sum _{i=0}^{n}(-1)^{i}{n+\alpha \choose n-i}L_{i}^{(\alpha )}(x),}$

while binomials have the parametrization

${\displaystyle {n+x \choose n}=\sum _{i=0}^{n}{\frac {\alpha ^{i}}{i!}}L_{n-i}^{(x+i)}(\alpha ).}$

${\displaystyle e^{-\gamma x}=\sum _{i=0}^{\infty }{\frac {\gamma ^{i}}{(1+\gamma )^{i+\alpha +1}}}L_{i}^{(\alpha )}(x)\qquad {\text{convergent iff }}\Re (\gamma )>-{\tfrac {1}{2}}}$

for the exponential function. The incomplete gamma function has the representation

${\displaystyle \Gamma (\alpha ,x)=x^{\alpha }e^{-x}\sum _{i=0}^{\infty }{\frac {L_{i}^{(\alpha )}(x)}{1+i}}\qquad \left(\Re (\alpha )>-1,x>0\right).}$

Multiplication theorems

Erdélyi gives the following two multiplication theorems [11]

${\displaystyle t^{n+1+\alpha }e^{(1-t)z}L_{n}^{(\alpha )}(zt)=\sum _{k=n}{k \choose n}\left(1-{\frac {1}{t}}\right)^{k-n}L_{k}^{(\alpha )}(z),}$
${\displaystyle e^{(1-t)z}L_{n}^{(\alpha )}(zt)=\sum _{k=0}{\frac {(1-t)^{k}z^{k}}{k!}}L_{n}^{(\alpha +k)}(z).}$

Relation to Hermite polynomials

The generalized Laguerre polynomials are related to the Hermite polynomials:

{\displaystyle {\begin{aligned}H_{2n}(x)&=(-1)^{n}2^{2n}n!L_{n}^{(-1/2)}(x^{2})\\H_{2n+1}(x)&=(-1)^{n}2^{2n+1}n!xL_{n}^{(1/2)}(x^{2})\end{aligned}}}

where the Hn(x) are the Hermite polynomials based on the weighting function exp(−x2), the so-called "physicist's version."

Because of this, the generalized Laguerre polynomials arise in the treatment of the quantum harmonic oscillator.

Relation to hypergeometric functions

The Laguerre polynomials may be defined in terms of hypergeometric functions, specifically the confluent hypergeometric functions, as

${\displaystyle L_{n}^{(\alpha )}(x)={n+\alpha \choose n}M(-n,\alpha +1,x)={\frac {(\alpha +1)_{n}}{n!}}\,_{1}F_{1}(-n,\alpha +1,x)}$

where ${\displaystyle (a)_{n}}$ is the Pochhammer symbol (which in this case represents the rising factorial).

Poisson Kernel

${\displaystyle \sum _{n=0}^{\infty }{\frac {n!L_{n}^{(\alpha )}(x)L_{n}^{(\alpha )}(y)r^{n}}{\Gamma \left(1+\alpha +n\right)}}={\frac {\exp \left(-{\frac {\left(x+y\right)r}{1-r}}\right)I_{\alpha }\left({\frac {2{\sqrt {xyr}}}{1-r}}\right)}{\left(xyr\right)^{\frac {\alpha }{2}}\left(1-r\right)}},\qquad \alpha >-1,\quad |r|<1.}$

Notes

1. {{#invoke:Citation/CS1|citation |CitationClass=journal }}
2. A&S p. 781
3. A&S p.509
4. A&S p.510
5. A&S p. 775
6. G. Szegő, "Orthogonal polynomials", 4th edition, Amer. Math. Soc. Colloq. Publ., vol. 23, Amer. Math. Soc., Providence, RI, 1975, p. 198.
7. D. Borwein, J. M. Borwein, R. E. Crandall, "Effective Laguerre asymptotics", SIAM J. Numer. Anal., vol. 46 (2008), no. 6, pp. 3285-3312, http://dx.doi.org/10.1137/07068031X
8. A&S equation (22.12.6), p. 785
9. W. Koepf, "Identities for families of orthogonal polynomials and special functions.", Integral Transforms and Special Functions 5, (1997) pp.69-102. (Theorem 10)
10. A&S p. 774
11. C. Truesdell, "On the Addition and Multiplication Theorems for the Special Functions", Proceedings of the National Academy of Sciences, Mathematics, (1950) pp.752-757.

References

• Template:Abramowitz Stegun ref
• Template:Dlmf
• B. Spain, M.G. Smith, Functions of mathematical physics, Van Nostrand Reinhold Company, London, 1970. Chapter 10 deals with Laguerre polynomials.
• {{#invoke:citation/CS1|citation

|CitationClass=citation }}

|CitationClass=book }}

• S. S. Bayin (2006), Mathematical Methods in Science and Engineering, Wiley, Chapter 3.