# Peeter Joot's (OLD) Blog.

• ## Archives

 Adam C Scott on avoiding gdb signal noise… Ken on Scotiabank iTrade RESP …… Alan Ball on Oops. Fixing a drill hole in P… Peeter Joot's B… on Stokes theorem in Geometric… Exploring Stokes The… on Stokes theorem in Geometric…

• 294,017

## Hydrogen like atom, and Laguerre polynomials.

Posted by peeterjoot on November 29, 2010

# Motivation.

For the hydrogen atom, after some variable substitutions the radial part of the Schr\”{o}dinger equation takes the form

\begin{aligned}\frac{d^2 R_l}{d\rho^2} + \frac{2}{\rho} \frac{d R_l}{d\rho} + \left( \frac{\lambda}{\rho} - \frac{l(l+1)}{\rho^2} - \frac{1}{{4}} \right) R_l = 0\end{aligned} \hspace{\stretch{1}}(1.1)

In [1] it is argued that the functions $R_l$ are of the form

\begin{aligned}R_l = \rho^s L(\rho) e^{-\rho/2}\end{aligned} \hspace{\stretch{1}}(1.2)

where $L$ is a polynomial in $\rho$, specifically Laguerre polynomials. Let’s look at some of those details a bit more closely.

# Guts

The first part of the argument comes from considering the $\rho \rightarrow \infty$ case, where Schr\”{o}dinger’s equation is approximately

\begin{aligned}\frac{d^2 R_l}{d\rho^2} - \frac{1}{{4}} R_l \approx 0.\end{aligned} \hspace{\stretch{1}}(2.3)

This large $\rho$ approximation has solutions $e^{\pm \rho/2}$, and we take the negative sign case as physically meaningful in order for the wave function to be normalizable.

Next it is argued that polynomial multiples of this will also be approximate solutions. Utilizing monomial multiple of the decreasing exponential as a trial solution, let’s compute how this fits into the radial Schr\”{o}dinger’s equation 1.1 above. Write

\begin{aligned}R_l = \rho^s e^{-\rho/2}\end{aligned} \hspace{\stretch{1}}(2.4)

The derivatives are

\begin{aligned}R_l' &= \rho^{s-1} \left( s -\frac{\rho}{2}\right) e^{-\rho/2} \\ R_l'' &=\rho^{s-2}\left( s (s-1) -s \rho +\frac{1}{4} \rho^2\right)e^{-\rho/2}\end{aligned}

and substitution yields

\begin{aligned}\rho^{s-2}e^{-\rho/2}\left((s - \rho) (s+1)+\lambda \rho- l(l+1)\right)\end{aligned} \hspace{\stretch{1}}(2.5)

There are two things that this can show. The first is that for $\rho \rightarrow \infty$ this produces a polynomial with degree $s-2$ and $s-1$ terms multiplied by the exponential, and we have approximately

\begin{aligned}\rho^{s-1}e^{-\rho/2}(\lambda - s - 1)\end{aligned} \hspace{\stretch{1}}(2.6)

The $s-1$ terms will dominate the polynomial, but the exponential dominate all, approaching zero for $\rho \rightarrow \infty$, just as the non-polynomial multiplied $e^{-\rho/2}$ approximate solution will. This confirms that in the limit this polynomial multiplied exponential still has the desired behavior in the large $\rho$ limit. Also observe that in the limit of small $\rho$ we have approximately

\begin{aligned}\rho^{s-2}e^{-\rho/2}\left(s (s+1) - l(l+1)\right)\end{aligned} \hspace{\stretch{1}}(2.7)

Since $\rho^{s-2} \rightarrow \infty$ as $\rho \rightarrow 0$, we require either a different trial solution, or $s=l$ to have a normalizable wavefunction.

Before settling on $s=l$ let’s compute the derivatives for a more general trial function, of the form 1.2, and substitute those. After a bit of computation we find

\begin{aligned}R_l' = \rho^{s-1} e^{-\rho/2} \left( \left( s - \frac{\rho}{2} \right) L + \rho L'\right)\end{aligned} \hspace{\stretch{1}}(2.8)

\begin{aligned}R_l'' = \rho^{s-2} e^{-\rho/2} \left(\left( s(s-1) - s \rho + \frac{\rho^2}{4} \right) L+\left( 2 s \rho -\rho^2 \right) L'+ \rho^2 L''\right)\end{aligned} \hspace{\stretch{1}}(2.9)

Putting these together and substitution back into 1.1 yields

\begin{aligned}0 = \rho^{s-2} e^{-\rho/2} \left(L \left( (s-\rho)(s+1) + \rho \lambda -l (l+1)\right)+\rho L' \left( 2 (s+1) -\rho \right)+ \rho^2 L''\right)\end{aligned} \hspace{\stretch{1}}(2.10)

In the $\rho \rightarrow 0$ limit where the $\rho^{s-2}$ terms dominate 2.11 becomes

\begin{aligned}0 \approx\rho^{s-2} L \left(s(s+1) - l(l+1)\right)\end{aligned} \hspace{\stretch{1}}(2.11)

Again, this provides the $s=l$ or $s = -(l+1)$ possibilities from the text, and we discard $s=-(l+1)$ due to non-normalizability. A side question. How does one solve integer equations like this?

## What remains?

With $s=l$ killing off the $\rho^{s-2}$ terms, what is our differential equation for $L$?

\begin{aligned}0 =\rho L''+L' \left( 2 (l+1) -\rho \right)+L \left( \lambda - (l+1) \right)\end{aligned} \hspace{\stretch{1}}(2.12)

Comparing this to [2] we have something pretty close to the stated differential equation for the Laguerre polynomial. Ours is of the form

\begin{aligned}0 =\rho L''+L' \left( m + 1 -\rho \right)+L n,\end{aligned} \hspace{\stretch{1}}(2.13)

where the differential equation in the wikipedia article has $m=0$. No change of variables involving a scalar multiplicative factor for $\rho$ appears to be able to get it into that form, and I am guessing this is the differential equation for the associated Laguerre polynomial (something not stated in the wikipedia article).

Let’s derive the recurrence relations for the coefficients, and work out the first few such polynomials to compare. Plugging in a polynomial of the form

\begin{aligned}L = \sum_{k=0}^r a_k \rho^{k},\end{aligned} \hspace{\stretch{1}}(2.14)

where $a_r$ is assumed to be non-zero. We also assume that this polynomial is not an infinite series (ruling out the infinite series with convergence arguments is covered nicely in the text).

we have for 2.13

\begin{aligned}0 &= \sum_{k=0}^r a_k\left(k (k-1) \rho^{k-1}+ k (m+1) \rho^{k-1}- k \rho^{k}+ n \rho^{k}\right) \\ &=\sum_{{k'}=1}^r\rho^{{k'}-1}a_{k'}{k'} \left({k'}-1 +(m +1)\right)+\sum_{k=0}^r\rho^{k}a_k\left(- k + n \right) \\ &=\sum_{{k}=0}^{r-1}\rho^{k}a_{k+1}(k+1) \left(k +(m+1)\right)+\sum_{k=0}^r\rho^{k}a_k\left(- k + n \right) \\ &=\sum_{{k}=0}^{r-1}\rho^{k}\Bigl(a_{k+1} (k+1) (k + m + 1)+a_k (n -k)\Bigr)+a_r (n-r) \rho^{r}\end{aligned}

Observe first that since we have assumed $a_r \ne 0$, we must have $r=n$. Requiring termwise equality with zero gives us the recurrance relation between the coefficents, for $k \in [0,n-1]$

\begin{aligned}a_{k+1} = a_k \frac{k - n}{ (k+1) (k + m + 1) }.\end{aligned} \hspace{\stretch{1}}(2.15)

Repeated application shows the pattern for these coefficients, and with $a_0=1$ we have

\begin{aligned}a_1 &= -\frac{n-0}{(1)(m+1)} \\ a_2 &= \frac{(n-1)(n-0)}{(2)(1)(m+2)(m+1)} \\ a_3 &= -\frac{(n-2)(n-1)(n-0)}{(3)(2)(1)(m+3)(m+2)(m+1)},\end{aligned}

With

\begin{aligned}a_k &= \frac{(-1)^k (n-(k-1))\cdots(n-1)(n-0)}{k!(m+k)\cdots(m+2)(m+1)} \\ &= \frac{(-1)^k n! m!}{k!(m+k)!(n-(k-1) -1)!},\end{aligned}

Or

\begin{aligned}a_k = \frac{(-1)^k n! m!}{k!(m+k)!(n-k)!}.\end{aligned} \hspace{\stretch{1}}(2.16)

Forming the complete series, we can get at the form of the associated Laguerre polynomials in the wikipedia article without too much trouble

\begin{aligned}L_n^m(\rho) &\propto 1 + \sum_{k=1}^n \frac{(-1)^k}{k!} \frac{n! m!}{(n-k)!(m+k)!} \rho^k \\ &\propto \frac{(n+m)!}{n!m!} + \sum_{k=1}^n \frac{(-1)^k}{k!} \frac{(n+m)!}{(n-k)!(m+k)!} \rho^k.\end{aligned}

Dropping the proportionality, this simplifies to just

\begin{aligned}L_n^m(\rho) = \sum_{k=0}^n \frac{(-1)^k}{k!} \binom{n+m}{m+k} \rho^k\end{aligned} \hspace{\stretch{1}}(2.17)

This isn’t neccessarily the form of the polynomials used in the text. To see if that is the case, we need to check the normalization.

According to the wikipedia article we have for the associated Laguerre polynomials as defined above

\begin{aligned}\int_0^{\infty}\rho^m e^{-\rho} L_n^{m}(\rho)L_{n'}^{m}(\rho)d\rho = \frac{(n+m)!}{n!}\delta_{n,{n'}}\end{aligned} \hspace{\stretch{1}}(2.18)

whereas in the text we have

\begin{aligned}\int_0^{\infty}\rho^{2l + 2} e^{-\rho} \left( L_{n+l}^{2l + 1}(\rho) \right)^2 d\rho = \frac{2n ((n+l)!)^3}{(n-l-1)!}.\end{aligned} \hspace{\stretch{1}}(2.19)

It seems clear that two different notations are being used. In this physical context of wave functions we want the normalization defined by

\begin{aligned}1 = \int_0^\infty \rho^2 R_l^2(\rho) d\rho = \int_0^\infty \rho^{2l + 2} e^{-\rho} L^2(\rho) d\rho\end{aligned} \hspace{\stretch{1}}(2.20)

Using the wikipedia notation, with

\begin{aligned}L(\rho) = A L_n^{2l+1},\end{aligned} \hspace{\stretch{1}}(2.21)

we want

\begin{aligned}1 &= \int \rho^{2l + 2} e^{-\rho} L^2(\rho) d\rho \\ &= A^2 \sum_{a,b=0}^n \frac{(-1)^{a+b}}{a!b!} \binom{n+2l+1}{2l+1+a}\binom{n+2l+1}{2l+1+b}\int_0^\infty d\rho \rho^{2l + 2 + a + b} e^{-\rho} \end{aligned}

Since $\int_0^\infty d\rho \rho^{a} e^{-\rho} = \Gamma(a+1) = a!$ we have

\begin{aligned}1 = A^2 \sum_{a,b=0}^n \frac{(-1)^{a+b}}{a!b!} \binom{n+m}{m+a}\binom{n+m}{m+b}(m + 1 + a + b)!\end{aligned} \hspace{\stretch{1}}(2.22)

It looks like there is probably some way to simplify this, and if so we’d be able to map the notation used (without definition) used in the text, to the notation used in the wikipedia article. If we don’t care about that, nor the specifics of the normalization constant then there’s not too much more to say.

This is an ugly kind of place to leave things, but that’s enough for today. It’s too bad that the text isn’t just more explicit, and it’s probably best to refer elsewhere for any more detail. With no specifics about the functions themselves in any form, one has to do that anyways.

# References

[1] BR Desai. Quantum mechanics with basic field theory. Cambridge University Press, 2009.

[2] Wikipedia. Laguerre polynomials — wikipedia, the free encyclopedia, 2010. [Online; accessed 29-November-2010]. http://en.wikipedia.org/w/index.php?title=Laguerre_polynomials&oldid=38%6787645.