Sub­sec­tions


A.6 Con­stant spher­i­cal po­ten­tials

This ad­den­dum de­scribes the so­lu­tions of the Hamil­ton­ian eigen­value prob­lem in spher­i­cal co­or­di­nates if the po­ten­tial is con­stant.

These so­lu­tions are im­por­tant for many rea­sons. For ex­am­ple, you might want to cre­ate a sim­pli­fied model for the hy­dro­gen atom that way. To do so, you could, for ex­am­ple, as­sume that the po­ten­tial en­ergy has a con­stant neg­a­tive value up to say the Bohr ra­dius and is zero be­yond it. That is not re­ally a very good model for the hy­dro­gen atom. How­ever, it works much bet­ter for nu­cle­ons in atomic nu­clei, chap­ter 14.12.

The so­lu­tions in this note are also im­por­tant for de­scrib­ing ex­per­i­ments in which par­ti­cles are scat­tered from some tar­get, {A.30}. And more fun­da­men­tally, they give the en­ergy states of def­i­nite an­gu­lar mo­men­tum for par­ti­cles in empty space.


A.6.1 The eigen­value prob­lem

The Hamil­ton­ian eigen­value prob­lem is

\begin{displaymath}
- \frac{\hbar^2}{2m} \nabla^2 \psi + V \psi = E \psi
\end{displaymath}

In this note it is as­sumed that the po­ten­tial $V$ is a con­stant in the ra­dial re­gion of in­ter­est.

To clean the prob­lem up a bit, take the po­ten­tial en­ergy term to the other side, and also the co­ef­fi­cient of the Lapla­cian. That pro­duces

\begin{displaymath}
- \nabla^2 \psi = \frac{p_{\rm {c}}^2}{\hbar^2} \psi
\end{displaymath}

where

\begin{displaymath}
p_{\rm {c}} \equiv \sqrt{2m(E-V)}
\end{displaymath}

The con­stant $p_{\rm {c}}$ is what clas­si­cal physics would take to be the lin­ear mo­men­tum of a par­ti­cle with to­tal en­ergy $E$ and po­ten­tial en­ergy $V$.


A.6.2 The eigen­func­tions

Be­cause the po­ten­tial is spher­i­cally sym­met­ric like for the hy­dro­gen atom, the eigen­func­tions are of sim­i­lar form:

\begin{displaymath}
\fbox{$\displaystyle
\psi_{Elm}(r,\theta,\phi) = R_{El}(r) Y_l^m(\theta,\phi)
$} %
\end{displaymath} (A.17)

Here the $Y_l^m$ are again the spher­i­cal har­mon­ics. These eigen­func­tions have def­i­nite square an­gu­lar mo­men­tum $l(l+1)\hbar^2$ where $l$ is an non­neg­a­tive in­te­ger. They also have def­i­nite an­gu­lar mo­men­tum in the cho­sen $z$-​di­rec­tion equal to $m\hbar$, where $m$ is an in­te­ger that sat­is­fies $\vert m\vert$ $\raisebox{-.3pt}{$\leqslant$}$ $l$.

The ra­dial func­tions $R_{El}$ in the eigen­func­tions $\psi_{Elm}$ are dif­fer­ent from those of the hy­dro­gen atom. De­pend­ing on what­ever is eas­i­est in a given ap­pli­ca­tion, they can be writ­ten in two ways, {D.16}.

The first way is as

\begin{displaymath}
\fbox{$\displaystyle
R_{El} = c_{\rm{s}} j_l(p_{\rm{c}}r/\...
...\rm{c}}r/\hbar)
\qquad p_{\rm{c}} \equiv \sqrt{2m(E-V)}
$} %
\end{displaymath} (A.18)

Here $c_{\rm {s}}$ and $c_{\rm {c}}$ are ar­bi­trary con­stants. The func­tions $j_l$ and $n_l$ are called the “spher­i­cal Bessel func­tions” of the first and sec­ond kinds. The $n_l$ are also called the “Neu­mann func­tions” and might in­stead be in­di­cated by $y_l$ or $\eta_l$.

Ex­pres­sions for these Bessel func­tions can be found in ad­vanced math­e­mat­i­cal hand­books, [1]:

\begin{displaymath}
\fbox{$\displaystyle
j_l(x)
= (-x)^l \left(\frac{1}{x} \frac{{\rm d}}{{\rm d}x}\right)^l \frac{\sin x}{x}
$} %
\end{displaymath} (A.19)


\begin{displaymath}
\fbox{$\displaystyle
n_l(x)
= - (-x)^l \left(\frac{1}{x} \frac{{\rm d}}{{\rm d}x}\right)^l\frac{\cos x}{x}
$} %
\end{displaymath} (A.20)

The spher­i­cal Bessel func­tions are of­ten con­ve­nient in a re­gion of con­stant po­ten­tial that in­cludes the ori­gin, be­cause the Bessel func­tions of the first kind $j_l$ give the so­lu­tions that are fi­nite at the ori­gin. (To see that, note that the Tay­lor se­ries of $\sin{x}$ di­vided by $x$ is a power se­ries in $x^2$, and that $x{\rm d}{x}$ $\vphantom0\raisebox{1.5pt}{$=$}$ $\frac12{\rm d}{x}^2$.) In par­tic­u­lar for small $x$:

\begin{displaymath}
\fbox{$\displaystyle
j_l(x)
= \frac{2^l l!}{(2l+1)!}x^l + O(x^{l+2})
\equiv \frac{x^l}{(2l+1)!!} + O(x^{l+2})
$} %
\end{displaymath} (A.21)

Here !! is one of these un­for­tu­nate no­ta­tions. The sec­ond ! means that all even fac­tors are dropped from the fac­to­r­ial.

The Bessel func­tions of the sec­ond kind are sin­gu­lar at the ori­gin and nor­mally do not ap­pear if the ori­gin is part of the con­sid­ered re­gion.

Also, the spher­i­cal Bessel func­tions are real for real $x$. How­ever, in a re­gion where the po­ten­tial $V$ is larger than the en­ergy $E$ of the par­ti­cles, the ar­gu­ment of the Bessel func­tions in (A.18) will be imag­i­nary.

The other way to write the ra­dial func­tions is as

\begin{displaymath}
\fbox{$\displaystyle
R_{El} = c_{\rm{f}} h_l^{(1)}(p_{\rm{...
...\rm{c}}r/\hbar)
\qquad p_{\rm{c}} \equiv \sqrt{2m(E-V)}
$} %
\end{displaymath}

where $h_l^{(1)}$ and $h_l^{(2)}$ are called the “spher­i­cal Han­kel func­tions.”

The spher­i­cal Han­kel func­tions can again be found in ad­vanced math­e­mat­i­cal hand­books, [1]:

\begin{displaymath}
\fbox{$\displaystyle
h_l^{(1)}(x) = - {\rm i}(-x)^l\left(\...
...\right)^l
\frac{e^{{\rm i}x}}{x} = j_l(x)+{\rm i}n_l(x)
$} %
\end{displaymath} (A.22)


\begin{displaymath}
\fbox{$\displaystyle
h_l^{(2)}(x) = {\rm i}(-x)^l \left(\f...
...right)^l
\frac{e^{-{\rm i}x}}{x} = j_l(x)-{\rm i}n_l(x)
$} %
\end{displaymath} (A.23)

The given ex­pres­sions in terms of the spher­i­cal Bessel func­tions are read­ily in­verted to give the Bessel func­tions in terms of the Han­kel func­tions,
\begin{displaymath}
j_l(x) = \frac{h_l^{(1)}(x)+h_l^{(2)}(x)}{2}
\qquad
n_l(x) = \frac{h_l^{(1)}(x)-h_l^{(2)}(x)}{2{\rm i}} %
\end{displaymath} (A.24)

For large $x$ the spher­i­cal Han­kel func­tions can be ap­prox­i­mated as

\begin{displaymath}
\fbox{$\displaystyle
h_l^{(1)}(x) \sim (-i)^{l+1} \frac{e^...
...qquad
h_l^{(2)}(x) \sim i^{l+1} \frac{e^{-{\rm i}x}}{x}
$} %
\end{displaymath} (A.25)

This as­ymp­totic be­hav­ior tends to make the Han­kel func­tions more con­ve­nient far from the ori­gin. Ex­po­nen­tials are math­e­mat­i­cally sim­pler and more fun­da­men­tal than the sines and cosines in the as­ymp­totic be­hav­ior of the Bessel func­tions.


A.6.3 About free space so­lu­tions

The most im­por­tant case for which the en­ergy eigen­func­tions of the pre­vi­ous sub­sec­tion ap­ply is for par­ti­cles in empty space. They de­scribe en­ergy states with def­i­nite square and $z$ an­gu­lar mo­men­tum. How­ever, some­times par­ti­cles in empty space are bet­ter de­scribed by states of def­i­nite lin­ear mo­men­tum. And in some cases, like in scat­ter­ing prob­lems, you need both types of so­lu­tion. Then you also need to con­vert be­tween them.

The en­ergy states in empty space with def­i­nite square and $z$ an­gu­lar mo­men­tum are

\begin{displaymath}
\fbox{$\displaystyle
\psi_{Elm} = j_l(p_{\rm{c}}r/\hbar) Y_l^m(\theta,\phi)
\qquad p_{\rm{c}} \equiv \sqrt{2mE}
$} %
\end{displaymath} (A.26)

These states have square an­gu­lar mo­men­tum $l(l+1)\hbar^2$ and an­gu­lar mo­men­tum in the cho­sen $z$-​di­rec­tion $m\hbar$. They are non­sin­gu­lar at the ori­gin.

A state that has def­i­nite lin­ear mo­men­tum $p_{\rm {c}}$ purely in the $z$-​di­rec­tion has an en­ergy eigen­func­tion

\begin{displaymath}
\fbox{$\displaystyle
\psi_{{\vec k}} = e^{{\rm i}p_{\rm{c}} z/\hbar}
\qquad p_{\rm{c}} \equiv \sqrt{2mE}
$} %
\end{displaymath} (A.27)

This eigen­func­tion is not nor­mal­ized, and can­not be nor­mal­ized. How­ever, nei­ther can the eigen­func­tion $\psi_{Elm}$ above be. It is the curse of eigen­func­tions in in­fi­nite empty space. An in­tro­duc­tion to the ad­just­ments that must be made to deal with this is given in chap­ter 7.9.

It is some­times nec­es­sary to write a lin­ear mo­men­tum eigen­func­tion of the form (A.27) in terms of an­gu­lar mo­men­tum ones of the form (A.26). Rayleigh worked out the cor­rect com­bi­na­tion a very long time ago, {D.16}:

\begin{displaymath}
\fbox{$\displaystyle
e^{{\rm i}p_{\rm{c}} z/\hbar}
= \sum...
...a)
\qquad
c_{{\rm{w}},l} = {\rm i}^l \sqrt{4\pi(2l+1)}
$} %
\end{displaymath} (A.28)

Note that only eigen­func­tions with $m$ $\vphantom0\raisebox{1.5pt}{$=$}$ 0 are needed.