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
where is an associated Laguerre polynomial and and are constants. This equation is not self-adjoint, but can be made so by multiplying by :
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
so the self-adjoint condition is satisfied.
We have seen that a self-adjoint operator operating on a function that satisfies certain boundary conditions is a Hermitian operator, and as such we can write the ODE in the form
where is the weighting function and is the eigenvalue. In this case
This means that the associated Laguerre polynomials for distinct values of must be orthogonal. That is if :
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):
Here the function is the coefficient of the second derivative term, so is
which is zero at both limits. Since the 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 . What we need to find is the normalization condition; that is, we need to find the integral of the square of a given polynomial:
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 so that
The variable here is just a dummy variable which allows the series expansion to be derived. The catch, of course, is how to obtain 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:
(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:
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 appears only in the exponential and the term. It can be evaluated by integrating by parts times, and noting that the integrated term vanishes at each stage due to the boundary conditions. We get (I cheated and used Maple):
where the final answer is obtained by putting everything over a common denominator and using .
Now what we need to do is expand this result as a series in powers of 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 (again, if you haven’t encountered binomial expansions, please accept this on faith):
In this case, the binomial coefficient must be generalized (since can be negative and non-integer):
In our case, , and so we get
so comparing coefficients of on both sides, the normalization integral comes out to
Multiplying this out we get
Plugging this back into 2 we get