next up previous contents
Next: The infinite potential well Up: Central potentials Previous: Introduction   Contents


Derivation of radial equation

Now, we have seen that the cartesian components of the momentum, ${\bf p}$, can be represented as (see Sect. 7.2)
\begin{displaymath}
p_i = -{\rm i}\,\hbar\,\frac{\partial}{\partial x_i}
\end{displaymath} (603)

for $i=1,2,3$, where $x_1\equiv x$, $x_2\equiv y$, $x_3\equiv z$, and ${\bf r}\equiv (x_1, x_2, x_3)$. Likewise, it is easily demonstrated, from the above expressions, and the basic definitions of the spherical polar coordinates [see Eqs. (527)-(529)], that the radial component of the momentum can be represented as
\begin{displaymath}
p_r \equiv \frac{{\bf p}\cdot{\bf r}}{r} = -{\rm i}\,\hbar\,\frac{\partial}{\partial r}.
\end{displaymath} (604)

Recall that the angular momentum vector, ${\bf L}$, is defined [see Eq. (508)]

\begin{displaymath}
{\bf L} = {\bf r}\times {\bf p}.
\end{displaymath} (605)

This expression can also be written in the following form:
\begin{displaymath}
L_i = \epsilon_{ijk}\,x_j\,p_k.
\end{displaymath} (606)

Here, the $\epsilon_{ijk}$ (where $i,j,k$ all run from 1 to 3) are elements of the so-called totally anti-symmetric tensor. The values of the various elements of this tensor are determined via a simple rule:
\begin{displaymath}
\epsilon_{ijk} = \left\{
\begin{array}{rcl}
0 &\mbox{\hspace...
...\ are anti-cyclic permutation of $1,2,3$}
\end{array}\right. .
\end{displaymath} (607)

Thus, $\epsilon_{123}=\epsilon_{231}=1$, $\epsilon_{321}=\epsilon_{132}=-1$, and $\epsilon_{112}=\epsilon_{131}=0$, etc. Equation (606) also makes use of the Einstein summation convention, according to which repeated indices are summed (from 1 to 3). For instance, $a_i\,b_i\equiv a_1\,b_1+a_2\,b_2+a_3\,b_3$. Making use of this convention, as well as Eq. (607), it is easily seen that Eqs. (605) and (606) are indeed equivalent.

Let us calculate the value of $L^2$ using Eq. (606). According to our new notation, $L^2$ is the same as $L_i\,L_i$. Thus, we obtain

\begin{displaymath}
L^2 = \epsilon_{ijk}\,x_j\,p_k\,\epsilon_{ilm}\,x_l\,p_m =
\epsilon_{ijk}\,\epsilon_{ilm}\,x_j\,p_k\,x_l\,p_m.
\end{displaymath} (608)

Note that we are able to shift the position of $\epsilon_{ilm}$ because its elements are just numbers, and, therefore, commute with all of the $x_i$ and the $p_i$. Now, it is easily demonstrated that
\begin{displaymath}
\epsilon_{ijk}\,\epsilon_{ilm}\equiv \delta_{jl}\,\delta_{km}-\delta_{jm}\,\delta_{kl}.
\end{displaymath} (609)

Here $\delta_{ij}$ is the usual Kronecker delta, whose elements are determined according to the rule
\begin{displaymath}
\delta_{ij} = \left\{
\begin{array}{rcl}
1 &\mbox{\hspace{1c...
...box{if $i$\ and $j$\ different}\\ [0.5ex]
\end{array}\right. .
\end{displaymath} (610)

It follows from Eqs. (608) and (609) that
\begin{displaymath}
L^2 = x_i\,p_j\,x_i\,p_j - x_i\,p_j\,x_j\,p_i.
\end{displaymath} (611)

Here, we have made use of the fairly self-evident result that $\delta_{ij}\,a_i\,b_j
\equiv a_i\,b_i$. We have also been careful to preserve the order of the various terms on the right-hand side of the above expression, since the $x_i$ and the $p_i$ do not necessarily commute with one another.

We now need to rearrange the order of the terms on the right-hand side of Eq. (611). We can achieve this by making use of the fundamental commutation relation for the $x_i$ and the $p_i$ [see Eq. (465)]:

\begin{displaymath}[x_i,p_j]= {\rm i}\,\hbar\,\delta_{ij}.
\end{displaymath} (612)

Thus,
$\displaystyle L^2$ $\textstyle =$ $\displaystyle x_i\left(x_i\,p_j - [x_i,p_j]\right) p_j
- x_i\,p_j\,\left(p_i\,x_j+[x_j,p_i]\right)$  
  $\textstyle =$ $\displaystyle x_i\,x_i\,p_j\,p_j - {\rm i}\,\hbar\,\delta_{ij}\,x_i\,p_j
-x_i\,p_j\,p_i\,x_j - {\rm i}\,\hbar\,\delta_{ij}\,x_i\,p_j$  
  $\textstyle =$ $\displaystyle x_i\,x_i\,p_j\,p_j -x_i\,p_i\,p_j\,x_j - 2\,{\rm i}\,\hbar\,x_i\,p_i.$ (613)

Here, we have made use of the fact that $p_j\,p_i=p_i\,p_j$, since the $p_i$ commute with one another [see Eq. (464)]. Next,
\begin{displaymath}
L^2 = x_i\,x_i\,p_j\,p_j - x_i\,p_i\left(x_j\,p_j - [x_j,p_j]\right) - 2\,{\rm i}\,\hbar\,x_i\,p_i.
\end{displaymath} (614)

Now, according to (612),
\begin{displaymath}[x_j,p_j]\equiv [x_1,p_1]+[x_2,p_2]+[x_3,p_3] = 3\,{\rm i}\,\hbar.
\end{displaymath} (615)

Hence, we obtain
\begin{displaymath}
L^2 = x_i\,x_i\,p_j\,p_j - x_i\,p_i\,x_j\,p_j + {\rm i}\,\hbar\,x_i\,p_i.
\end{displaymath} (616)

When expressed in more conventional vector notation, the above expression becomes
\begin{displaymath}
L^2 = r^2\,p^2 - ({\bf r}\cdot{\bf p})^2 + {\rm i}\,\hbar\,{\bf r}\cdot{\bf p}.
\end{displaymath} (617)

Note that if we had attempted to derive the above expression directly from Eq. (605), using standard vector identities, then we would have missed the final term on the right-hand side. This term originates from the lack of commutation between the $x_i$ and $p_i$ operators in quantum mechanics. Of course, standard vector analysis assumes that all terms commute with one another.

Equation (617) can be rearranged to give

\begin{displaymath}
p^2 = r^{-2}\left[({\bf r}\cdot{\bf p})^2- {\rm i}\,\hbar\,{\bf r}\cdot{\bf p}+L^2\right].
\end{displaymath} (618)

Now,
\begin{displaymath}
{\bf r}\cdot{\bf p} = r\,p_r = -{\rm i}\,\hbar\,r\,\frac{\partial}{\partial r},
\end{displaymath} (619)

where use has been made of Eq. (604). Hence, we obtain
\begin{displaymath}
p^2 = -\hbar^2\left[\frac{1}{r}\frac{\partial}{\partial r}\l...
...}\frac{\partial}{\partial r}- \frac{L^2}{\hbar^2\,r^2}\right].
\end{displaymath} (620)

Finally, the above equation can be combined with Eq. (602) to give the following expression for the Hamiltonian:
\begin{displaymath}
H = -\frac{\hbar^2}{2\,m}\left(\frac{\partial^2}{\partial r^...
...{\partial}{\partial r}- \frac{L^2}{\hbar^2\,r^2}\right)
+V(r).
\end{displaymath} (621)

Let us now consider whether the above Hamiltonian commutes with the angular momentum operators $L_z$ and $L^2$. Recall, from Sect. 8.3, that $L_z$ and $L^2$ are represented as differential operators which depend solely on the angular spherical polar coordinates, $\theta $ and $\phi$, and do not contain the radial polar coordinate, $r$. Thus, any function of $r$, or any differential operator involving $r$ (but not $\theta $ and $\phi$), will automatically commute with $L^2$ and $L_z$. Moreover, $L^2$ commutes both with itself, and with $L_z$ (see Sect. 8.2). It is, therefore, clear that the above Hamiltonian commutes with both $L_z$ and $L^2$.

Now, according to Sect. 4.10, if two operators commute with one another then they possess simultaneous eigenstates. We thus conclude that for a particle moving in a central potential, the eigenstates of the Hamiltonian are simultaneous eigenstates of $L_z$ and $L^2$. Now, we have already found the simultaneous eigenstates of $L_z$ and $L^2$--they are the spherical harmonics, $Y_{l,m}(\theta,\phi)$, discussed in Sect. 8.7. It follows that the spherical harmonics are also eigenstates of the Hamiltonian. This observation leads us to try the following separable form for the stationary wave-function:

\begin{displaymath}
\psi(r,\theta,\phi) = R(r)\,Y_{l,m}(\theta,\phi).
\end{displaymath} (622)

It immediately follows, from (535) and (536), and the fact that $L_z$ and $L^2$ both obviously commute with $R(r)$, that
$\displaystyle L_z\,\psi$ $\textstyle =$ $\displaystyle m\,\hbar\,\psi,$ (623)
$\displaystyle L^2\,\psi$ $\textstyle =$ $\displaystyle l\,(l+1)\,\hbar^2\,\psi.$ (624)

Recall that the quantum numbers $m$ and $l$ are restricted to take certain integer values, as explained in Sect. 8.6.

Finally, making use of Eqs. (601), (621), and (624), we obtain the following differential equation which determines the radial variation of the stationary wave-function:

\begin{displaymath}
-\frac{\hbar^2}{2\,m}\left(\frac{d^2}{d r^2}
+ \frac{2}{r}\f...
...- \frac{l\,(l+1)}{r^2}\right)R_{n,l}
+V\,R_{n,l} = E\,R_{n,l}.
\end{displaymath} (625)

Here, we have labeled the function $R(r)$ by two quantum numbers, $n$ and $l$. The second $l$ is, of course, related to the eigenvalue of $L^2$. [Note that the azimuthal quantum number, $m$, does not appear in the above equation, and, therefore, does not influence either the function $R(r)$ or the energy, $E$.] As we shall see, the first, $n$, is determined by the constraint that the radial wave-function be square-integrable.


next up previous contents
Next: The infinite potential well Up: Central potentials Previous: Introduction   Contents
Richard Fitzpatrick 2006-12-12