time independent Schr\"odinger equation in spherical coordinatesTimeIndependentSchrodingerEquationInSphericalCoordinates2006-01-15 15:34:542008-03-26 09:15:47Definition% this is the default PlanetMath preamble. as your knowledge
% of TeX increases, you will probably want to edit this, but
% it should be fine as is for beginners.
% almost certainly you want these
\usepackage{amssymb}
\usepackage{amsmath}
\usepackage{amsfonts}
% used for TeXing text within eps files
%\usepackage{psfrag}
% need this for including graphics (\includegraphics)
%\usepackage{graphicx}
% for neatly defining theorems and propositions
%\usepackage{amsthm}
% making logically defined graphics
%\usepackage{xypic}
% there are many more packages, add them here as you need them
% define commands hereWhen writing the time independent Schr\"odinger equation in spherical coordinates, we need to plug the Laplacian in spherical coordinates into the time independent Schr\"odinger equation. The Laplacian was found to be
\[
\nabla _{sph}^{2} = \frac{1}{r^2} \frac{\partial}{\partial r}\left(r^2 \frac{\partial}{\partial r}\right) + \frac{1}{r^2 sin\theta} \frac{\partial}{\partial \theta} \left( sin \theta \frac{\partial}{\partial \theta}\right) + \frac{1}{r^2 sin^2 \theta} \frac{\partial^2}{\partial \phi^2}
\]
Using the three dimensional Schr\"odinger equation we then have
\[
\hat{H} \psi(r,\theta, \phi) = -\frac{\hbar^2}{2m} \left [ \frac{1}{r^2} \frac{\partial }{\partial r}\left(r^2 \frac{\partial \psi(r,\theta, \phi)}{\partial r}\right) \frac{1}{r^2 sin\theta} \frac{\partial}{\partial \theta} \left( sin \theta \frac{\partial \psi(r,\theta, \phi)}{\partial \theta}\right) + \frac{1}{r^2 sin^2 \theta} \frac{\partial^2 \psi(r,\theta, \phi)}{\partial \phi^2} \right ] + V(r,\theta, \phi) \psi(r,\theta, \phi) = E\psi(r, \theta, \phi)
\]
We can gain insight into this somewhat ugly equation by rewriting it using the square of the angular momentum operator in spherical polar coordinates:
\[
\hat{L}^2 = {1 \over \sin\theta}
{\partial\over\partial\theta}\left(\sin\theta{\partial\over\partial\theta}\right) +\frac{1}{\sin^2\theta}\frac{\partial^2}{\partial \phi^2}
\]
This leads to
\[
\left(-\frac{\hbar^2}{2m}\left({1 \over r^2} {\partial\over\partial r}\left(r^2{\partial\over\partial r}\right) \right)+\frac{1}{2m}\frac{\hat{L}^2}{r^2}+V(r, \theta, \phi)\right)\psi(r, \theta, \phi) = E\psi(r, \theta, \phi)
\]
\subsection{Spherically symmetric separable solution}
This equation is only exactly solvable if $V=V(r)$, a function without angular dependence. We then write $\psi(r, \theta, \phi) = R(r)Y(\theta, \phi)$ leading to the following equation:
\begin{align*}
\left(-\frac{\hbar^2}{2m}\left({1 \over r^2} {\partial\over\partial r}\left(r^2 \frac{\partial}{\partial r}\right)\right)+\frac{1}{2m}\frac{\hat{L}^2}{r^2}+V(r, \theta, \phi)\right)\psi(r, \theta, \phi)R(r) Y(\theta, \phi) &= E R(r) Y(\theta, \phi)\\
-\frac{\hbar^2}{2m} \left( Y(\theta, \phi) \left( \frac{1}{r^2}{\partial \over \partial r}\left(r^2 \frac{\partial}{\partial r}\right) \right)R(r) \right) + \frac{R(r)}{2m}\frac{\hat{L}^2 \,Y(\theta, \phi)}{r^2} + V(r)R(r)Y(\theta, \phi) &= E R(r)(Y(\theta, \phi)
\end{align*}
To solve this equation we need to remove the angular dependence. This is simply done by substituting the eigenfunctions of $\hat{L}^2$ into the equation. These are known to be the spherical harmonics, $Y^m_l(\theta, \phi)$. We also know that these have eigenvalues $\hbar^2 l(l+1)$, i.e.
\[
\hat{L}^2 \,Y_l^m(\theta, \phi) = \hbar^2 l(l+1) Y^m_l(\theta, \phi)
\]
We now substitute this result into the Schr\"odinger equation and divide through by a common factor of $Y_l^m(\theta, \phi)$
\[
\left( -\frac{\hbar^2}{2m} \left(\frac{1}{r^2}{\partial \over \partial r} \left(r^2 \frac{\partial}{\partial r}\right) +\frac{\hbar^2l(l+1)}{r^2}\right) + V(r) \right) R(r) = E R(r)
\]
This is the radial equation.\\
{\bf References}
[1] Griffiths, D. "Introduction to Quantum Mechanics" Prentice Hall, New Jersey, 1995.