Laguerre polynomials
{{Short description|Sequence of differential equation solutions}}
In mathematics, the Laguerre polynomials, named after Edmond Laguerre (1834–1886), are nontrivial solutions of Laguerre's differential equation:
y = y(x)
which is a second-order linear differential equation. This equation has nonsingular solutions only if {{mvar|n}} is a non-negative integer.
Sometimes the name Laguerre polynomials is used for solutions of
where {{mvar|n}} is still a non-negative integer.
Then they are also named generalized Laguerre polynomials, as will be done here (alternatively associated Laguerre polynomials or, rarely, Sonine polynomials, after their inventor{{cite journal|title=Recherches sur les fonctions cylindriques et le développement des fonctions continues en séries|author=N. Sonine|journal=Math. Ann.|date=1880|volume=16| issue=1|pages=1–80|doi=10.1007/BF01459227|s2cid=121602983|url=http://www.digizeitschriften.de/dms/img/?PPN=PPN235181684_0016&DMDID=dmdlog8}} Nikolay Yakovlevich Sonin).
More generally, a Laguerre function is a solution when {{mvar|n}} is not necessarily a non-negative integer.
The Laguerre polynomials are also used for Gauss–Laguerre quadrature to numerically compute integrals of the form
These polynomials, usually denoted {{math|L0}}, {{math|L1}}, ..., are a polynomial sequence which may be defined by the Rodrigues formula,
reducing to the closed form of a following section.
They are orthogonal polynomials with respect to an inner product
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 that is larger by a factor of n
Recursive definition, closed form, and generating function
One can also define the Laguerre polynomials recursively, defining the first two polynomials as
and then using the following recurrence relation for any {{math|k ≥ 1}}:
Furthermore,
In solution of some boundary value problems, the characteristic values can be useful:
The closed form is
The generating function for them likewise follows,
The operator form is
Polynomials of negative index can be expressed using the ones with positive index:
class="wikitable" style="margin:0.5em auto"
|+A table of the Laguerre polynomials | |
width="20%" | n
! | |
---|---|
align="center" | 0 | |
align="center" | 1 | |
align="center" | 2
| | |
align="center" | 3
| | |
align="center" | 4
| | |
align="center" | 5
| | |
align="center" | 6
| | |
align="center" | 7
| | |
align="center" | 8
| | |
align="center" | 9
| | |
align="center" | 10
| | |
align="center" | n
| |
Generalized Laguerre polynomials
For arbitrary real α the polynomial solutions of the differential equationA&S p. 781
are called generalized Laguerre polynomials, or associated Laguerre polynomials.
One can also define the generalized Laguerre polynomials recursively, defining the first two polynomials as
and then using the following recurrence relation for any {{math|k ≥ 1}}:
The simple Laguerre polynomials are the special case {{math|1=α = 0}} of the generalized Laguerre polynomials:
The Rodrigues formula for them is
= \frac{x^{-\alpha}}{n!}\left( \frac{d}{dx}-1\right)^nx^{n+\alpha}.
The generating function for them is
= Properties =
- Laguerre functions are defined by confluent hypergeometric functions and Kummer's transformation asA&S p. 509 where is a generalized binomial coefficient. When {{mvar|n}} is an integer the function reduces to a polynomial of degree {{mvar|n}}. It has the alternative expressionA&S p. 510 in terms of Kummer's function of the second kind.
- The closed form for these generalized Laguerre polynomials of degree {{mvar|n}} isA&S p. 775 derived by applying Leibniz's theorem for differentiation of a product to Rodrigues' formula.
- Laguerre polynomials have a differential operator representation, much like the closely related Hermite polynomials. Namely, let and consider the differential operator . Then .{{citation needed|date=April 2023}}
- The first few generalized Laguerre polynomials are:
class="wikitable" style="margin:0.5em auto" |
width="20%"| n
! |
---|
align="center" | 0
| |
align="center" | 1
| |
align="center" | 2
| |
align="center" | 3
| |
align="center" | 4
| |
align="center" | 5
| |
align="center" | 6
| |
align="center" | 7
| |
align="center" | 8
| |
align="center" | 9
| |
align="center" | 10
| |
- The coefficient of the leading term is {{math|(−1)n/n
! }}; - The constant term, which is the value at 0, is
- If {{math|α}} is non-negative, then Ln(α) has n real, strictly positive roots (notice that is a Sturm chain), which are all in the interval {{citation needed|date=September 2011}}
- The polynomials' asymptotic behaviour for large {{mvar|n}}, but fixed {{mvar|α}} and {{math|x > 0}}, is given bySzegő, p. 198.D. Borwein, J. M. Borwein, R. E. Crandall, "Effective Laguerre asymptotics", SIAM J. Numer. Anal., vol. 46 (2008), no. 6, pp. 3285–3312 {{doi|10.1137/07068031X}}
\begin{align}
& 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}}} \sin\left(2 \sqrt{nx}- \frac{\pi}{2}\left(\alpha-\frac{1}{2} \right) \right)+O\left(n^{\frac{\alpha}{2}-\frac{3}{4}}\right), \\[6pt]
& L_n^{(\alpha)}(-x) = \frac{(n+1)^{\frac{\alpha}{2}-\frac{1}{4}}}{2\sqrt{\pi}} \frac{e^{-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),
\end{align}
and summarizing by where 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
where the contour circles the origin once in a counterclockwise direction without enclosing the essential singularity at 1
= Recurrence relations =
The addition formula for Laguerre polynomials:{{Cite web |title=DLMF: §18.18 Sums ‣ Classical Orthogonal Polynomials ‣ Chapter 18 Orthogonal Polynomials |url=https://dlmf.nist.gov/18.18 |access-date=2025-03-18 |website=dlmf.nist.gov}}
Laguerre's polynomials satisfy the recurrence relations
in particular
and
or
moreover
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{align}
They can be used to derive the four 3-point-rules
L_n^{(\alpha)}(x) &= L_n^{(\alpha+1)}(x) - L_{n-1}^{(\alpha+1)}(x) = \sum_{j=0}^k {k \choose j}(-1)^j L_{n-j}^{(\alpha+k)}(x), \\[10pt]
n L_n^{(\alpha)}(x) &= (n + \alpha )L_{n-1}^{(\alpha)}(x) - x L_{n-1}^{(\alpha+1)}(x), \\[10pt]
& \text{or } \\
\frac{x^k}{k!}L_n^{(\alpha)}(x) &= \sum_{i=0}^k (-1)^i {n+i \choose i} {n+\alpha \choose k-i} L_{n+i}^{(\alpha-k)}(x), \\[10pt]
n L_n^{(\alpha+1)}(x) &= (n-x) L_{n-1}^{(\alpha+1)}(x) + (n+\alpha)L_{n-1}^{(\alpha)}(x) \\[10pt]
x L_n^{(\alpha+1)}(x) &= (n+\alpha)L_{n-1}^{(\alpha)}(x)-(n-x)L_n^{(\alpha)}(x);
\end{align}
combined they give this additional, useful recurrence relations
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{align}
Since is a monic polynomial of degree in ,
there is the partial fraction decomposition
\frac{n!\,L_n^{(\alpha)}(x)}{(\alpha+1)_n}
&= 1- \sum_{j=1}^n (-1)^j \frac{j}{\alpha + j} {n \choose j}L_n^{(-j)}(x) \\
&= 1- \sum_{j=1}^n \frac{x^j}{\alpha + j}\,\,\frac{L_{n-j}^{(j)}(x)}{(j-1)!} \\
&= 1-x \sum_{i=1}^n \frac{L_{n-i}^{(-\alpha)}(x) L_{i-1}^{(\alpha+1)}(-x)}{\alpha +i}.
\end{align}
The second equality follows by the following identity, valid for integer i and {{mvar|n}} and immediate from the expression of in terms of Charlier polynomials:
For the third equality apply the fourth and fifth identities of this section.
= Derivatives =
Differentiating the power series representation of a generalized Laguerre polynomial {{mvar|k}} times leads to
(-1)^k L_{n-k}^{(\alpha+k)}(x) & \text{if } k\le n, \\
0 & \text{otherwise.}
\end{cases}
This points to a special case ({{math|1=α = 0}}) of the formula above: for integer {{math|1=α = k}} the generalized polynomial may be written
the shift by {{mvar|k}} sometimes causing confusion with the usual parenthesis notation for a derivative.
Moreover, the following equation holds:
which generalizes with Cauchy's formula to
The derivative with respect to the second variable {{mvar|α}} has the form,{{Cite journal | doi=10.1080/10652469708819127 | title = Identities for families of orthogonal polynomials and special functions| journal=Integral Transforms and Special Functions | volume=5| issue=1–2| pages=69–102|year = 1997|last1 = Koepf|first1 = Wolfram| citeseerx=10.1.1.298.7657}}
The generalized Laguerre polynomials obey the differential equation
which may be compared with the equation obeyed by the kth derivative of the ordinary Laguerre polynomial,
where for this equation only.
In Sturm–Liouville form the differential equation is
which shows that {{math|L{{su|b=n|p=(α)}}}} is an eigenvector for the eigenvalue {{mvar|n}}.
= Orthogonality =
The generalized Laguerre polynomials are orthogonal over {{closed-open|0, ∞}} with respect to the measure with weighting function {{math|xα e−x}}:{{Cite web | url=http://mathworld.wolfram.com/AssociatedLaguerrePolynomial.html | title=Associated Laguerre Polynomial}}
which follows from
If denotes the gamma distribution then the orthogonality relation can be written as
The associated, symmetric kernel polynomial has the representations (Christoffel–Darboux formula){{citation needed|date=October 2011}}
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}}\\[4pt]
& =\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}} \\[4pt]
&= \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{align}
recursively
Moreover,{{clarify|post-text=Limit as n goes to infinity?|date=January 2016}}
Turán's inequalities can be derived here, which is
The following integral is needed in the quantum mechanical treatment of the hydrogen atom,
= Series expansions =
Let a function have the (formal) series expansion
Then
The series converges in the associated Hilbert space {{math|L2[0, ∞)}} if and only if
== Further examples of expansions==
Monomials are represented as
while binomials have the parametrization
This leads directly to
for the exponential function. The incomplete gamma function has the representation
In quantum mechanics
In quantum mechanics the Schrödinger equation for the hydrogen-like atom is exactly solvable by separation of variables in spherical coordinates. The radial part of the wave function is a (generalized) Laguerre polynomial.{{Cite book|title=Quantum Mechanics in Chemistry|last=Ratner, Schatz|first=Mark A., George C.|publisher=Prentice Hall|year=2001|location=0-13-895491-7| pages=90–91}}
Vibronic transitions in the Franck-Condon approximation can also be described using Laguerre polynomials.{{Cite journal|last1=Jong|first1=Mathijs de|last2=Seijo|first2=Luis|last3=Meijerink|first3=Andries| last4=Rabouw |first4=Freddy T.| date=2015-06-24|title=Resolving the ambiguity in the relation between Stokes shift and Huang–Rhys parameter |url=https://pubs.rsc.org/en/content/articlelanding/2015/cp/c5cp02093j|journal=Physical Chemistry Chemical Physics|language=en| volume=17 |issue=26|pages=16959–16969|doi=10.1039/C5CP02093J|pmid=26062123|bibcode=2015PCCP...1716959D|hdl=1874/321453|s2cid=34490576 | issn=1463-9084|hdl-access=free}}
Multiplication theorems
Erdélyi gives the following two multiplication theorems C. Truesdell, "[http://www.pnas.org/cgi/reprint/36/12/752.pdf On the Addition and Multiplication Theorems for the Special Functions]", Proceedings of the National Academy of Sciences, Mathematics, (1950) pp. 752–757.
& t^{n+1+\alpha} e^{(1-t) z} L_n^{(\alpha)}(z t)=\sum_{k=n}^\infty {k \choose n}\left(1-\frac 1 t\right)^{k-n} L_k^{(\alpha)}(z), \\[6pt]
& e^{(1-t)z} L_n^{(\alpha)}(z t)=\sum_{k=0}^\infty \frac{(1-t)^k z^k}{k!}L_n^{(\alpha+k)}(z).
\end{align}
Relation to Hermite polynomials
The generalized Laguerre polynomials are related to the Hermite polynomials:
H_{2n}(x) &= (-1)^n 2^{2n} n! L_n^{(-1/2)} (x^2) \\[4pt]
H_{2n+1}(x) &= (-1)^n 2^{2n+1} n! x L_n^{(1/2)} (x^2)
\end{align}
where the {{math|Hn(x)}} are the Hermite polynomials based on the weighting function {{math|exp(−x2)}}, the so-called "physicist's version."
Because of this, the generalized Laguerre polynomials arise in the treatment of the quantum harmonic oscillator.
Applying the addition formula,
=\sum_{m_1+\cdots+m_r=n} \prod_{i=1}^r H_{2m_i}(z_i).
Relation to hypergeometric functions
The Laguerre polynomials may be defined in terms of hypergeometric functions, specifically the confluent hypergeometric functions, as
where is the Pochhammer symbol (which in this case represents the rising factorial).
Hardy–Hille formula
The generalized Laguerre polynomials satisfy the Hardy–Hille formulaSzegő, p. 102.W. A. Al-Salam (1964), [https://projecteuclid.org/euclid.dmj/1077375084 "Operational representations for Laguerre and other polynomials"], Duke Math J. 31 (1): 127–142.
where the series on the left converges for and . Using the identity
(see generalized hypergeometric function), this can also be written as
where denotes the modified Bessel function of the first kind, defined as
I_\alpha(z) = \sum_{k=0}^\infty \frac{1}{k!\, \Gamma(k+\alpha+1)} \left(\frac{z}{2}\right)^{2k+\alpha}
This formula is a generalization of the Mehler kernel for Hermite polynomials, which can be recovered from it by setting the Hermite polynomials as a special case of the associated Laguerre polynomials.
Substitute and take the limit, we obtain Szegő, page 102, Equation (5.1.16)
\sum_{n=0}^\infty
\frac{t^n}{\Gamma(n+1+\alpha)} L_n^{(\alpha)}(x)
=
\frac{e^t}{(-xt)^{\alpha/2}}I_{\alpha}(2\sqrt{-xt}).
Physics convention
The generalized Laguerre polynomials are used to describe the quantum wavefunction for hydrogen atom orbitals.{{cite book |last1=Griffiths |first1=David J. |title=Introduction to quantum mechanics |date=2005 |publisher=Pearson Prentice Hall |location=Upper Saddle River, NJ |isbn=0131118927 |edition=2nd}}{{cite book |last1=Sakurai |first1=J. J. |title=Modern quantum mechanics |date=2011 |publisher=Addison-Wesley |location=Boston |isbn=978-0805382914 |edition=2nd}}{{cite book |last1=Merzbacher |first1=Eugen |title=Quantum mechanics |date=1998 |publisher=Wiley |location=New York |isbn=0471887021 |edition=3rd}} The convention used throughout this article expresses the generalized Laguerre polynomials as {{cite book |last1=Abramowitz |first1=Milton |title=Handbook of mathematical functions, with formulas, graphs, and mathematical tables |date=1965 |publisher=Dover Publications |location=New York |isbn=978-0-486-61272-0}}
where is the confluent hypergeometric function.
In the physics literature, the generalized Laguerre polynomials are instead defined as
The physics version is related to the standard version by
There is yet another, albeit less frequently used, convention in the physics literature {{cite book |last1=Schiff |first1=Leonard I. |title=Quantum mechanics |date=1968 |publisher=McGraw-Hill |location=New York |isbn=0070856435 |edition=3d}}{{cite book |last1=Messiah |first1=Albert |title=Quantum Mechanics. |date=2014 |publisher=Dover Publications |isbn=9780486784557}}{{cite book |last1=Boas |first1=Mary L. |title=Mathematical methods in the physical sciences |date=2006 |publisher=Wiley |location=Hoboken, NJ |isbn=9780471198260 |edition=3rd}}
Umbral calculus convention
Generalized Laguerre polynomials are linked to Umbral calculus by being Sheffer sequences for when multiplied by . In Umbral Calculus convention,{{Cite journal |last1=Rota |first1=Gian-Carlo |last2=Kahaner |first2=D |last3=Odlyzko |first3=A |date=1973-06-01 |title=On the foundations of combinatorial theory. VIII. Finite operator calculus |journal=Journal of Mathematical Analysis and Applications |language=en |volume=42 |issue=3 |pages=684–760 |doi=10.1016/0022-247X(73)90172-8 |issn=0022-247X|doi-access=free }} the default Laguerre polynomials are defined to bewhere are the signless Lah numbers. is a sequence of polynomials of binomial type, ie they satisfy
See also
- Orthogonal polynomials
- Rodrigues' formula
- Angelescu polynomials
- Bessel polynomials
- Denisyuk polynomials
- Transverse mode, an important application of Laguerre polynomials to describe the field intensity within a waveguide or laser beam profile.
Notes
{{Reflist|35em}}
References
- {{Abramowitz_Stegun_ref|22|773}}
- G. Szegő, Orthogonal polynomials, 4th edition, Amer. Math. Soc. Colloq. Publ., vol. 23, Amer. Math. Soc., Providence, RI, 1975.
- {{dlmf|id=18|title=Orthogonal Polynomials|first=Tom H.|last=Koornwinder|first2=Roderick S. C.|last2= Wong|first3=Roelof | last3=Koekoek ||first4=René F. |last4=Swarttouw}}
- B. Spain, M.G. Smith, Functions of mathematical physics, Van Nostrand Reinhold Company, London, 1970. Chapter 10 deals with Laguerre polynomials.
- {{springer|title=Laguerre polynomials|id=p/l057310}}
- Eric W. Weisstein, "[http://mathworld.wolfram.com/LaguerrePolynomial.html Laguerre Polynomial]", From MathWorld—A Wolfram Web Resource.
- {{cite book|author=George Arfken and Hans Weber| title= Mathematical Methods for Physicists| publisher=Academic Press| year=2000| isbn = 978-0-12-059825-0 }}
External links
- {{cite web|author=Timothy Jones|url=http://www.physics.drexel.edu/~tim/open/hydrofin | title=The Legendre and Laguerre Polynomials and the elementary quantum mechanical model of the Hydrogen Atom}}
- {{MathWorld|title=Laguerre polynomial|id=LaguerrePolynomial}}
{{Authority control}}