Laguerre polynomials – normalization

Required math: calculus

Required physics: none

Reference: Arfken, George B. & Weber, Hans J. (2005), Mathematical Methods for Physicists, 6th Edition, Academic Press – Sec 13.2.

The associated Laguerre ODE turns up in physics in several places, most notably in the solution of the Schrödinger equation for the hydrogen atom. The ODE is

\displaystyle  x\frac{d^{2}L_{n}^{k}(x)}{dx^{2}}+(k+1-x)\frac{dL_{n}^{k}(x)}{dx}+nL_{n}^{k}(x)=0
where {L_{n}^{k}(x)} is an associated Laguerre polynomial and {k} and {n} are constants. This equation is not self-adjoint, but can be made so by multiplying by {e^{-x}x^{k}}:

\displaystyle  e^{-x}x^{k+1}\frac{d^{2}L_{n}^{k}(x)}{dx^{2}}+(k+1-x)e^{-x}x^{k}\frac{dL_{n}^{k}(x)}{dx}+ne^{-x}x^{k}L_{n}^{k}(x)=0

The self-adjoint condition is that the derivative of the function multiplying the second derivative term equals the function multiplying the first derivative term. In this case

\displaystyle   \frac{d}{dx}\left(e^{-x}x^{k+1}\right) \displaystyle  = \displaystyle  -e^{-x}x^{k+1}+(k+1)e^{-x}x^{k}
\displaystyle  \displaystyle  = \displaystyle  (k+1-x)e^{-x}x^{k}

so the self-adjoint condition is satisfied.

We have seen that a self-adjoint operator {D} operating on a function that satisfies certain boundary conditions is a Hermitian operator, and as such we can write the ODE in the form

\displaystyle  DL_{n}^{k}(x)+\lambda w(x)L_{n}^{k}(x)=0
where {w(x)} is the weighting function and {\lambda} is the eigenvalue. In this case

\displaystyle   w(x) \displaystyle  = \displaystyle  e^{-x}x^{k}
\displaystyle  \lambda \displaystyle  = \displaystyle  n

This means that the associated Laguerre polynomials for distinct values of {n} must be orthogonal. That is if {m\ne n}:

\displaystyle  \int_{0}^{\infty}e^{-x}x^{k}L_{n}^{k}(x)L_{m}^{k}(x)\: dx=0\ \ \ \ \ (1)

Note that this condition relies on choosing the limits of integration so that the boundary conditions are satisfied, and in this case these boundary conditions are (see post on self-adjoint operators):

\displaystyle  L_{n}^{k}(x)p_{0}L_{m}^{k}(x)'-L_{m}^{k}(x)p_{0}L_{n}^{k}(x)'\Big|_{0}^{\infty}=0

the function {p_{0}(x)}
is the coefficient of the second derivative term, so is

\displaystyle  p_{0}(x)=e^{-x}x^{k+1} which is zero at both limits. Since the {L_{n}^{k}(x)} are polynomials, the exponential dominates the behaviour at infinity, so the boundary conditions are satisfied.

This shows that the Laguerre polynomials are orthogonal, provided they are integrated with the appropriate weighting function {w(x)}. What we need to find is the normalization condition; that is, we need to find the integral of the square of a given polynomial:

\displaystyle  \int_{0}^{\infty}e^{-x}x^{k}\left(L_{n}^{k}(x)\right)^{2}dx

To do this, we can make use of the generating function for the Laguerre polynomials. Generating functions are somewhat magical, and deriving them is often trickier than using them. In the case of the Laguerre polynomials, the generating function is obtained by the use of Cauchy’s theorem for integrals in the complex plane, which would take us too far afield for a single post. So I will ask the reader to accept the formula and leave the derivation for another post.

But first, what exactly is the generating function? It is a compact formula which, when expanded in a series, allows individual Laguerre polynomials to be picked off as the coefficients of individual powers in the series. That is, we define the generating function {g(x,z)} so that

\displaystyle  g(x,z)=\sum_{n=0}^{\infty}L_{n}^{k}(x)z^{n}

The variable {z} here is just a dummy variable which allows the series expansion to be derived. The catch, of course, is how to obtain {g(x,z)} so it gives just the right series. It is here that I’m asking the reader to accept the function on faith at this stage. The generating function for associated Laguerre polynomials is:

\displaystyle  \frac{e^{-xz/(1-z)}}{(1-z)^{k+1}}=\sum_{n=0}^{\infty}L_{n}^{k}(x)z^{n}
(As I said, it’s not something you can pull out of a hat.)

Now that we have the generating function, we can use it to derive the normalization condition. We square the above equation, multiply it by the weighting function and then integrate the result:

\displaystyle   \int_{0}^{\infty}e^{-x}x^{k}\frac{e^{-2xz/(1-z)}}{(1-z)^{2k+2}}dx \displaystyle  = \displaystyle  \int_{0}^{\infty}e^{-x}x^{k}\left(\sum_{n=0}^{\infty}L_{n}^{k}(x)z^{n}\right)^{2}dx
\displaystyle  \displaystyle  = \displaystyle  \sum_{n=0}^{\infty}z^{2n}\int_{0}^{\infty}e^{-x}x^{k}\left(L_{n}^{k}(x)\right)^{2}dx

The second line is obtained by using the orthogonality condition 1, since this causes all the cross terms to vanish out of the sum once they are integrated.

The integral on the left is not that hard to do, since the integration variable {x} appears only in the exponential and the {x^{k}} term. It can be evaluated by integrating by parts {k} times, and noting that the integrated term vanishes at each stage due to the boundary conditions. We get (I cheated and used Maple):

\displaystyle   \int_{0}^{\infty}e^{-x}x^{k}\frac{e^{-2xz/(1-z)}}{(1-z)^{2k+2}}dx \displaystyle  = \displaystyle  \frac{k!}{[z(1-z)]^{k+1}}\left(2+\frac{1-z}{z}\right)^{-(k+1)}
\displaystyle  \displaystyle  = \displaystyle  \frac{k!}{(1-z^{2})^{k+1}}

where the final answer is obtained by putting everything over a common denominator and using {(1-z)(1+z)=1-z^{2}}.

Now what we need to do is expand this result as a series in powers of {z} and match up corresponding terms with the desired integral in the series above. The series corresponding to this term is a negative binomial expansion, so we can use the formula, valid for any (even complex) number {r} (again, if you haven’t encountered binomial expansions, please accept this on faith):

\displaystyle  (t+y)^{r}=\sum_{n=0}^{\infty}{r \choose n}t^{r-n}y^{n}

In this case, the binomial coefficient must be generalized (since {r} can be negative and non-integer):

\displaystyle  {r \choose n}=\frac{r(r-1)\ldots (r-n+1)}{n!}

In our case, {r=-k-1}, {t=1} and {y=-z^{2}} so we get

\displaystyle   \frac{k!}{(1-z^{2})^{k+1}} \displaystyle  = \displaystyle  k!\sum_{n=0}^{\infty}{-k-1 \choose n}(-z^{2})^{n}
\displaystyle  \displaystyle  = \displaystyle  k!\sum_{n=0}^{\infty}{-k-1 \choose n}(-1)^{n}z^{2n}

so comparing coefficients of {z^{2n}} on both sides, the normalization integral comes out to

\displaystyle  \int_{0}^{\infty}e^{-x}x^{k}\left(L_{n}^{k}(x)\right)^{2}dx=k!{-k-1 \choose n}(-1)^{n}\ \ \ \ \ (2)

Multiplying this out we get

\displaystyle   {-k-1 \choose n} \displaystyle  = \displaystyle  \frac{(-k-1)(-k-2)\ldots (-k-1-n+1)}{n!}
\displaystyle  \displaystyle  = \displaystyle  (-1)^{n}\frac{(k+1)(k+2)\ldots (k+n)}{n!}
\displaystyle  \displaystyle  = \displaystyle  (-1)^{n}\frac{(n+k)!}{n!k!}

Plugging this back into 2 we get

\displaystyle   \int_{0}^{\infty}e^{-x}x^{k}\left(L_{n}^{k}(x)\right)^{2}dx \displaystyle  = \displaystyle  (-1)^{2n}\frac{(n+k)!}{n!}
\displaystyle  \displaystyle  = \displaystyle  \frac{(n+k)!}{n!}

3 thoughts on “Laguerre polynomials – normalization

  1. Pingback: Hermite polynomials – recursion relations « Physics tutorials

  2. Pingback: Hydrogen atom – series solution « Physics tutorials

  3. maestro13

    The norm of Laguerre polynomials can be determined far far easier – no generating series needed and also no joggling around with equations as done above.
    It involves only a couple of things: (1) the expression relating a Laguerre polynomial to the nth derivative of e^-x . x^(n+k), (2) partial integration, (3) the expression dL(n,k)/dx = L(n-1, k+1), (4) induction.
    Try it yourself!


Leave a Reply

Your email address will not be published. Required fields are marked *

You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <s> <strike> <strong>