Skip to main content
Physics LibreTexts

4.4: Orbital Eigenfunctions in 3-D

  • Page ID
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    The Angular Momentum Operators in Spherical Polar Coordinates

    The angular momentum operator


    In spherical polar coordinates,

    \[ \begin{matrix} x=r\sin\theta\cos\phi \\ y=r\sin\theta\sin\phi \\ z=r\cos\theta \\ ds^2=dr^2+r^2d\theta^2+r^2\sin^2\theta d\phi^2 \end{matrix}\label{4.4.1}\]

    the gradient operator is

    \[ \vec{\nabla}=\hat{r}\dfrac{\partial}{\partial r}+\hat{\theta}\dfrac{1}{r}\dfrac{\partial}{\partial \theta}+\hat{\phi}\dfrac{1}{r\sin\theta}\dfrac{\partial}{\partial \phi} \label{4.4.2}\]

    where now the little hats denote unit vectors: \(\hat{r}\) is radially outwards, \(\hat{\theta}\) points along a line of longitude away from the north pole (and therefore in the direction of increasing \(\theta\) ) and \(\hat{\phi}\) points along a line of latitude in an anticlockwise direction as seen looking down on the north pole (that is, in the direction of increasing \(\phi\)).


    Here \(\hat{r},\hat{\theta},\hat{\phi}\) form an orthonormal local basis, and

    \[ \hat{r}\times\hat{\theta}=\hat{\phi}, \;\; \hat{r}\times\hat{\phi}=-\hat{\theta} \label{4.4.3}\]

    as should be clear from the diagram.


    \[ \vec{r}\times\vec{\nabla}=\hat{\phi}\dfrac{\partial}{\partial \theta}-\hat{\theta}\dfrac{1}{\sin\theta}\dfrac{\partial}{\partial \phi} \label{4.4.4}\]

    (Explicitly, \(\hat{\phi}=(-\sin\phi,\cos\phi,0)\) and \(\hat{\theta}=(\cos\theta\cos\phi, \cos\theta\sin\phi,-\sin\theta)\).)

    The vector \(\hat{\phi}\) has zero component in the z-direction, the vector \(\hat{\theta}\) has component \(-\sin\theta\) in the z-direction, so we can immediately conclude that

    \[ L_z=(\vec{r}\times\vec{p})_z=-(i\hbar\vec{r}\times\vec{\nabla})_z=-i\hbar\dfrac{\partial}{\partial \phi} \label{4.4.5}\]

    just as in the two-dimensional case.

    The operator \[ L^2=-\hbar^2\left( \hat{\phi}\dfrac{\partial}{\partial \theta}-\hat{\theta}\dfrac{1}{\sin\theta}\dfrac{\partial}{\partial \phi} \right) \cdot \left( \hat{\phi}\dfrac{\partial}{\partial \theta}-\hat{\theta}\dfrac{1}{\sin\theta}\dfrac{\partial}{\partial \phi} \right). \label{4.4.6}\]

    To evaluate this expression, we use \(\hat{\phi^2}=1,\; \hat{\theta^2}=1,\; \hat{\phi}\cdot\hat{\theta}=0\) but we must also check the effects of the differential operators in the first expression on the variables in the second, including the unit vectors.

    From the explicit coordinate expressions for the unit vectors, or by staring at the diagram, you should be able to establish the following: \(\partial\hat{\phi}/\partial\theta=0, \; \partial\hat{\theta}/\partial\theta\) is in the r-direction, \(\partial\hat{\phi}/\partial\phi\) is a horizontal unit vector pointing inwards perpendicular to \(\hat{\phi}\), and having component \(-\cos\theta\) in the \(\hat{\theta}\)-direction, \(\partial\hat{\theta}/\partial\phi=\hat{\phi}\cos\theta\).

    Therefore, the only "differentiation of a unit vector" term that contributes to \(L^2\) is

    \[ \hbar^2\hat{\theta}\dfrac{1}{\sin\theta}\cdot\dfrac{\partial\hat{\phi}}{\partial \phi}\dfrac{\partial}{\partial \theta}=-\hbar^2\cot\theta\dfrac{\partial}{\partial \theta}. \label{4.4.7}\]

    The \(\hat{\phi}\dfrac{\partial}{\partial \theta}\) acting on the \(\sin\theta\) in \(-\hat{\theta}\dfrac{1}{\sin\theta}\dfrac{\partial}{\partial\phi}\) contributes nothing because \(\hat{\phi}\cdot\hat{\theta}=0\).


    \[ \begin{matrix} L^2=-\hbar^2\left( \dfrac{\partial^2}{\partial\theta^2}+\cot\theta\dfrac{\partial}{\partial\theta}+\dfrac{1}{\sin^2\theta}\dfrac{\partial^2}{\partial\phi^2} \right) \\ =-\hbar^2\left( \dfrac{1}{\sin\theta}\dfrac{\partial}{\partial\theta}\sin\theta\dfrac{\partial}{\partial\theta}+\dfrac{1}{\sin^2\theta}\dfrac{\partial^2}{\partial\phi^2} \right) \end{matrix} \label{4.4.8}\]

    Now, we know that \(L^2\) and \(L_z\) have a common set of eigenkets (since they commute) and we've already established that those of \(L_z\) are \(\Phi_m(\phi)=e^{im\phi}/\sqrt{2\pi}\), with \(m\) an integer, so the eigenkets of \(L^2\) must have this same \(f\) dependence, so they must be of the form \(\Theta^m_l(\theta)\Phi(\phi)\), where \(\Theta^m_l(\theta)\) is a (suitably normalized) solution of the equation

    \[ \dfrac{1}{\sin\theta}\dfrac{d}{d\theta}\sin\theta\dfrac{d\Theta^m_l(\theta)}{d\theta}-\dfrac{m^2}{\sin^2\theta}\Theta^m_l(\theta)=-l(l+1)\Theta^m_l(\theta) \label{4.4.9}\]

    more conveniently written \[ \sin\theta\dfrac{d}{d\theta}\sin\theta\dfrac{d\Theta^m_l(\theta)}{d\theta}+(l(l+1)\sin^2\theta-m^2)\Theta^m_l(\theta)=0. \label{4.4.10}\]

    To summarize: the solutions to this differential equation, with integer \(l,\; m,\; |m|\le l,\) will (together with \(\Phi_m(\phi)\)) give the complete set of eigenstates of \(L^2\), \(L_z\) in the coordinate representation.

    Finding the m = l Eigenket of \(L^2\), \(L_z\)

    Recall now that for the simple harmonic oscillator, the easiest wave function to find was that of the ground state, the solution of the simple linear equation \(\hat{a}\Psi_0=0\) (as well as being a solution of the quadratic Schrödinger equation, of course). The other state wave functions could then be found by applying the creation operator in differential form the necessary number of times.

    A similar strategy works here: we can easily find the highest state on the \(l\) ladder, \(m = l\), the state \(|l,l\rangle\), since it satisfies the linear equation \(L_+|l,l\rangle=0\), where \(L_+=L_x+iL_y\). We just need to cast this equation in coordinate form. In Cartesian coordinates, \(L_+=-i\hbar(\hat{r}\times\vec{\nabla})_+\), and we've already shown that

    \[ \vec{r}\times\vec{\nabla}=\hat{\phi}\dfrac{\partial}{\partial \theta}-\hat{\theta}\dfrac{1}{\sin\theta}\dfrac{\partial}{\partial \phi}. \label{4.4.10B}\]

    Therefore \[ (\vec{r}\times\vec{\nabla})_+=\hat{\phi_+}\dfrac{\partial}{\partial \theta}-\hat{\theta_+}\dfrac{1}{\sin\theta}\dfrac{\partial}{\partial \phi}, \label{4.4.11}\]

    and using \(\hat{\phi}=(-\sin\phi,\cos\phi,0)\), \(\hat{\theta}=(\cos\theta\cos\phi, \cos\theta\sin\phi,-\sin\theta\) we see that \(\hat{\theta_+}\), the component of \(\hat{\theta}\) in the + direction, is \(\phi_+=\phi_x+i\phi_y=ie^{i\phi}\), and similarly \(\hat{\theta_+}=\cos\theta e^{i\phi}\).

    So \[ \begin{matrix} L_+=\hbar e^{i\phi}\left( \dfrac{\partial}{\partial\theta}+i\cot\theta\dfrac{\partial}{\partial\phi}\right) \\ L_-=-\hbar e^{-i\phi}\left( \dfrac{\partial}{\partial\theta}-i\cot\theta\dfrac{\partial}{\partial\phi}\right) \end{matrix}\label{4.4.12}\]

    and \(L_+|l,l\rangle=0\) becomes \[ \left( \dfrac{\partial}{\partial\theta}+i\cot\theta\dfrac{\partial}{\partial\phi}\right) \Theta^l_l(\theta)e^{il\phi}=0 \label{4.4.13}\]

    That is, \[ \left( \dfrac{d}{d\theta}-l\cot\theta\right) \Theta^l_l(\theta)=0 \label{4.4.14}\]

    The solution to this equation is \[ \Theta^l_l(\theta)=N(\sin\theta)^l \label{4.4.15}\]

    where \(N\) is the normalization constant. The \(m\neq l\) wave functions are generated by applying the lowering operator \(L_-\).

    Normalizing the m = l Eigenket

    The standard notation for the normalized eigenkets \(|l,m\rangle\) is \(Y^m_l(\theta,\phi)=\Theta^m_l(\theta)\Phi_m(\phi)\). These functions, being eigenkets of Hermitian operators with different eigenvalues, must satisfy \[ \int^{\pi}_{\theta=0}\int^{2\pi}_{\phi=0} Y^{m'^{*}}_{l'}(\theta,\phi)Y^m_l(\theta,\phi)\sin\theta d\theta d\phi=\int Y^{m'^*}_{l'}(\theta,\phi)Y^m_l(\theta,\phi)d\Omega=\delta_{l'l}\delta_{m'm} \; . \label{4.4.16}\]

    So, our first job is to normalize \(\Theta^l_l(\theta)=N(\sin\theta)^l\) (taking \(\Phi_l(\phi)=e^{il\phi}/\sqrt{2\pi}\) already normalized) \[ |N|^2\int^{\pi}_0 (\sin\theta)^{2l+1}d\theta=1 \label{4.4.17}\]

    The integral can be evaluated using the substitution \(\mu=\cos\theta\) to give \(\int^1_{-1}(1-\mu^2)^l d\mu\), then making the further substitution \(u=\dfrac{1}{2}(1-\mu)\) to give \(2^{2l+1}\int^1_0 u^l(1-u)^l du\), which can be integrated by parts to give \[ |N|^2 2^{2l+1}(l!)^2/(2l+1)!=1. \label{4.4.18}\]

    Therefore \[ Y^l_l(\theta,\phi)=(-1)^l \left( \dfrac{(2l+1)!}{4\pi}\right)^{1/2}\dfrac{1}{2^l l!}(\sin\theta)^l e^{il\phi}=c_l(\sin\theta)^l e^{il\phi} \label{4.4.19}\]

    where we have fixed the sign in accord with the standard convention, and we will denote the rather cumbersome normalization constant by \(c_l\).

    Notice that for large values of \(l\), this function is heavily weighted around the equator, as we would expect -- for a given total angular momentum one gets a maximum component in the z-direction when the motion is concentrated in the x, y plane. This looks like a Bohr orbit.

    Finding the Rest of the Eigenkets: the Details

    Now that \(|l,l\rangle\) is normalized, we can automatically produce correctly normalized \(|l,m\rangle\)'s, since we know the matrix element of the lowering operator between normalized states. We don't have to do any more integrals.

    For example, \(L_-|l,l\rangle=\hbar\sqrt{2l}|l,l-1\rangle\), equivalently (the \(\hbar\)'s of course cancel) \[ Y^{l-1}_l(\theta,\phi)=\dfrac{(-1)}{\sqrt{2l}}e^{-i\phi}\left( \dfrac{\partial}{\partial\theta}-i\cot\theta\dfrac{\partial}{\partial\phi}\right) Y^l_l . \label{4.4.20}\]

    That is, \[ \begin{matrix} Y^{l-1}_l(\theta,\phi)=c_l e^{-i\phi}\left( \dfrac{\partial}{\partial\theta}-i\cot\theta\dfrac{\partial}{\partial\phi}\right) \sin^l\theta\cdot e^{il\phi} \\ =-c_l e^{i(l-1)\phi}\sqrt{2l}\sin^{l-1}\theta\cos\theta \end{matrix}\label{4.4.21}\]

    (both terms giving equal contributions).

    Note that this function is actually zero on the equator, but for large \(l\) it peaks close to the equator (on both sides).

    In principle, we can reapply this differential operator over and over to generate all the \(|l,m\rangle\) states, but this gets very messy. However, there is a neat theorem concerning the lowering operator that makes it all straightforward: \[ L_-e^{im\phi}f(\theta)=e^{i(m-1)\phi}\left( \sin^{1-m}\theta\dfrac{d}{d(\cos\theta)}\sin^m\theta \right) f(\theta) \label{4.4.22}\]

    Exercise: prove this.

    So \[ L_-e^{il\phi}\sin^l\theta=e^{i(l-1)\phi}\left( \sin^{1-l}\theta\dfrac{d}{d(\cos\theta)}\sin^l\theta \right) \sin^l(\theta) \label{4.4.23}\]

    and applying the operator again, \[ \begin{matrix} (L_-)^2e^{il\phi}\sin^l\theta=L_-e^{i(l-1)\phi}\left( \sin^{1-l}\theta\dfrac{d}{d(\cos\theta)}\sin^l\theta \right) \sin^l(\theta) \\ =e^{i(l-2)\phi}\left( \sin^{2-l}\theta\dfrac{d}{d(\cos\theta)}\sin^{l-1}\theta \right) \cdot\left( \sin^{1-l}\theta\dfrac{d}{d(\cos\theta)}\sin^l\theta \right) \sin^l(\theta) \\ =e^{i(l-2)\phi}\left( \sin^{2-l}\theta\dfrac{d^2}{d^2(\cos\theta)}\sin^l\theta \right) \sin^l(\theta). \end{matrix} \label{4.4.23}\]

    So the point of introducing this odd-looking representation of the lowering operator is that the \(\sin^{1-l}\theta\) term in the middle is exactly canceled when the operator is applies twice, and similar cancellations occur on repeating the operation, giving the (relatively) simple representation: \[ Y^m_l(\theta,\phi)=c_l\sqrt{\dfrac{(l+m)!}{(2l)!(l-m)!}}e^{im\phi}\sin^{-m}\theta\dfrac{d^{l-m}}{d(\cos\theta)^{l-m}}\sin^{2l}\theta \label{4.4.24}\]

    (Where did all those factorials come from? They're the product of all the inverse square root factors in \[ |l,m-1\rangle=\dfrac{1}{\sqrt{(l+m)(l-m+1)}}L_-|l,m\rangle \label{4.4.25}\]

    for the number of lowerings necessary.)

    Note that for \(m = 0\) the function is \[ Y^0_l(\theta,\phi)=c_l\sqrt{\dfrac{1}{(2l)!}}\dfrac{d^l}{d(\cos\theta)^l}\sin^{2l}\theta, \label{4.4.26}\]

    and in fact not a function of \(\phi\) at all. This isn't surprising, since it has zero angular momentum about the z-direction, the appropriate \(\Phi(\phi)\) is just constant.

    For \(m=-l\), the differentiation becomes trivial, because, writing \(\cos\theta=\mu\), the differentiation becomes \(\dfrac{d^{2l}}{d\mu^{2l}}(1-\mu^2)^l\) and only the \(\mu^{2l}\) term survives, giving \[ Y^{-l}_l(\theta,\phi)=(-1)^l c_l e^{-il\phi} \sin^l\theta . \label{4.4.27}\]

    Of course, this could also have been found from the linear equation \(L_-|l,-l\rangle=0\), and we could have equally generated all the states by applying \(L_+\) to this state. In fact, this gives a different -- but of course equivalent -- expression for the \(Y^m_l(\theta,\phi)\): \[ Y^m_l(\theta,\phi)=(-1)^m c_l \sqrt{\dfrac{(l-m)!}{(2l)!(l+m)!}}e^{im\phi}\sin^{m}\theta\dfrac{d^{l+m}}{d(\cos\theta)^{l+m}}\sin^{2l}\theta \label{4.4.28}\]

    (from Messiah, page 522).

    Relating the Yl m's to the Legendre Functions

    The Legendre polynomials \(P_n(\cos\theta)\) are defined by: \[ \begin{matrix} P_n(\cos\theta)=\dfrac{1}{2^n n!}\dfrac{d^n}{d(\cos\theta)^n}\sin^{2n}\theta, \; or\; \\ P_n(\mu)=\dfrac{1}{2^n n!}\dfrac{d^n}{d\mu^n}(1-\mu^2)^n \end{matrix}\label{4.4.29}\]

    where \(\mu=\cos\theta\), so \(d\mu=-\sin\theta d\theta\). From this form, it is easy to show that \(P_n(1)=1\) (all \(n\) differentiations must take out a \((1-\mu^2)\) factor to give a nonzero contribution), and \(P_n(\mu)\) must have \(n\) zeros in the interval (-1, 1). \(P_n(\mu)\) alternates between an even function and an odd function.

    The normalization of the \(P_n(\mu)\)'s is \[ \begin{matrix} \int^1_{-1} (P_n(\mu))^2 d\mu=\left( \dfrac{1}{2^n n!}\right)^2 \int^1_{-1} \dfrac{d^n}{d\mu^n}(\mu^2-1)^n\dfrac{d^n}{d\mu^n}(\mu^2-1)^n d\mu \\ =(-1)^n\left( \dfrac{1}{2^n n!}\right)^2 \int^1_{-1} (\mu^2-1)^n\dfrac{d^{2n}}{d\mu^{2n}}(\mu^2-1)^n d\mu \\ =(2n)! \left( \dfrac{1}{2^n n!}\right)^2 \int^1_{-1} (\mu^2-1)^n d\mu \\ =\dfrac{2}{2n+1} \end{matrix} \label{4.4.30}\]

    where in that last line we used the result for the integral obtained earlier in this lecture for normalizing \(Y^l_l\).

    Doing the same repeated integration by parts for two different Legendre polynomials proves they are orthogonal, \[ \int^1_{-1} P_m(\mu)P_n(\mu)d\mu=0, \; m\neq n \label{4.4.31}\]

    The associated Legendre functions are defined (for n and m zero or positive integers, \(n\ge m\)image090.png) by: \[ \begin{matrix} P^m_n(\mu)=(1-\mu^2)^{m/2}\dfrac{d^m}{d\mu^m}P_n(\mu) \\ =(-1)^n \dfrac{(1-\mu^2)^{m/2}}{2^n n!}\dfrac{d^{n+m}}{d\mu^{n+m}}(1-\mu^2)^n . \end{matrix} \label{4.4.32}\]

    Following Messiah in requiring \(Y^0_l(0,0)\) be real and positive, we find \[ Y^0_l(\theta,\phi)=\sqrt{\dfrac{2l+1}{4\pi}}P_l(\cos\theta) \label{4.4.33}\]

    where the coefficient just reflects the differing normalization conventions. Similarly, the spherical harmonics with nonzero m are proportional to the associated Legendre functions (the odd ones are not polynomials in \(\cos\theta\), despite Shankar p. 337, since they include odd powers of \(\sin\theta\)), \[ Y^m_l(\theta,\phi)=\sqrt{\dfrac{(2l+1)(l-m)!}{4\pi(l+m)!}}(-1)^m e^{im\phi} P^m_l(\cos\theta). \label{4.4.34}\]

    The Spherical Harmonics as a Basis

    We have found explicit expressions for the spherical harmonics: an orthonormal set of eigenfunctions of \(L^2\) and \(L_z\) defined on the surface of a sphere, \[ \int^{\pi}_{\theta=0}\int^{2\pi}_{\phi=0} Y^{m'^{*}}_{l'}(\theta,\phi)Y^m_l(\theta,\phi)\sin\theta d\theta d\phi=\int Y^{m'^*}_{l'}(\theta,\phi)Y^m_l(\theta,\phi)d\Omega=\delta_{l'l}\delta_{m'm} \; . \label{4.4.16}\]

    They form a complete set: \[ \sum^{\infty}_{l=0}\sum^{l}_{m=-l} |l,m\rangle \langle l,m|=I \label{4.4.35}\]

    or \[ \sum^{\infty}_{l=0}\sum^{l}_{m=-l} Y^{m^*}_l(\theta,\phi)Y^{m}_l(\theta',\phi')=\delta(\cos\theta-\cos\theta')\delta(\phi-\phi')=\delta(\Omega-\Omega') \label{4.4.36}\]

    in the notation of Messiah, where W refers to a point on the spherical surface.

    (Formal proof of the completeness is given in Byron and Fuller, Mathematics of Classical and Quantum Physics.)

    The above equation could also be written \[ \sum^{\infty}_{l=0}\sum^{l}_{m=-l} \langle\theta,\phi|l,m\rangle\langle l,m|\theta',\phi'\rangle=\langle\theta,\phi|\theta',\phi'\rangle=\delta(\cos\theta-\cos\theta')\delta(\phi-\phi') \label{4.4.37}\]

    where the ket \(|\theta',\phi'\rangle\) is to be understood as a localized ket, the spherical-surface version of \(|x\rangle\), normalized by its delta function inner product with the bra \(\langle\theta,\phi|\), exactly analogous to \(\langle x|x'\rangle=\delta(x-x')\), bearing in mind that the infinitesimal area element is \(-d(\cos\theta)d\phi\), (a positive quantity in the relevant interval, 0 to \(\pi\)).

    This completeness means that any reasonable function on the surface of the sphere can be expressed as a sum over spherical harmonics with appropriate coefficients, in other words the spherical generalization of a Fourier series.

    In fact, \(L^2\) is equivalent to \(\nabla^2\) on the spherical surface, so the \(Y^m_l\) are the eigenfunctions of the operator \(\nabla^2\). Just as in one dimension the eigenfunctions of \(d^2/dx^2\) have the spatial dependence of the eigenmodes of a vibrating string, the spherical harmonics have the spatial dependence of the eigenmodes of a vibrating spherical balloon. Of course, to describe the displacement of the balloon skin (which must be real!) with these eigenfunctions, we can no longer use the eigenfunctions of the z-component of angular momentum, since they are complex except in the trivial zero case. We must rearrange the eigenfunctions of \(L^2\), for example replacing the pair \(e^{i\phi},\; e^{-i\phi}\) with \(\cos\phi,\sin\phi\). These real solutions, essentially \[ \dfrac{1}{\sqrt{2}}(|l,l\rangle\pm|l,-l\rangle) \label{4.4.38}\]

    , have \(l\) nodal lines (zeroes) of longitude. Moving down one notch in \(|m|\), the (real) state with \(|m|=l-1\) has \(l-1\) longitudinal nodes, but has added a latitudinal node: the equator. Then \(|m|=l-2\) has \(l-2\) longitudinal nodes, 2 latitudinal nodal lines -- there are always \(l\) nodal lines total.

    Some of these modes of vibration have been observed in the sun after a sunspot storm. The spherical harmonics are also used in analyzing the cosmic background radiation.

    Some Low Order Spherical Harmonics

    Let's look in more detail at the lowest order spherical harmonics. For the first few, the normalization of the highest state \(|l,l\rangle\) is pretty easy to do from scratch: factoring out the \(\phi\) dependence as before, \(Y^m_l(\theta,\phi)=\Theta^m_l(\theta)\Phi_m(\phi)\), and taking the normalized \(\Phi_m(\phi)=e^{im\phi}/\sqrt{2\pi}\), the \(\theta\) normalization for \(|l,l\rangle\) is just

    \[ |N|^2\int^{\pi}_0 (\sin\theta)^{2l+1}d\theta=1 \label{4.4.17}\]

    This is easily accomplished for \(l=0,1,2\).

    All we then need is \[ L_{\pm}=\pm\hbar e^{\pm i\phi}\left( \dfrac{\partial}{\partial\theta}\pm i\cot\theta\dfrac{\partial}{\partial\phi}\right) \label{4.4.39}\]

    , \[ L_-|l,m\rangle=\hbar\sqrt{l(l+1)-m(m-1)}|l,m-1\rangle \label{4.4.40}\]

    , and finally the sign convention that \(Y^0_l(0,0)\) be real and positive.

    With a few elementary steps, it can be established that: \[ \begin{matrix} Y^0_0=\dfrac{1}{\sqrt{4\pi}} \\ Y^1_1=-\sqrt{\dfrac{3}{8\pi}}\sin\theta e^{i\phi} \\ Y^0_1=\sqrt{\dfrac{3}{4\pi}}\cos\theta \\ Y^{-1}_1=\sqrt{\dfrac{3}{8\pi}}\sin\theta e^{-i\phi} \end{matrix} \label{4.4.41}\]

    \[ \begin{matrix} Y^2_2=\sqrt{\dfrac{15}{32\pi}}\sin^2\theta e^{2i\phi},\; Y^1_2=-\sqrt{\dfrac{15}{8\pi}}\sin\theta\cos\theta e^{i\phi},\; Y^0_2=\sqrt{\dfrac{5}{16\pi}}(3\cos^2\theta-1) \\ Y^{-2}_2=\sqrt{\dfrac{15}{32\pi}}\sin^2\theta e^{-2i\phi},\; Y^{-1}_2=\sqrt{\dfrac{15}{8\pi}}\sin\theta\cos\theta e^{-i\phi},\; \end{matrix} \label{4.4.42}\]


    It is often useful to write the \(Y^m_l\) in terms of Cartesian coordinates, \[ (x,y,z)=(r\sin\theta\cos\phi,\; r\sin\theta\sin\phi,\; r\cos\theta) \label{4.4.43}\]

    so \[ Y^1_1(x,y,z)=-\sqrt{\dfrac{3}{8\pi}}\cdot \dfrac{x+iy}{r},\; Y^0_1(x,y,z)=\sqrt{\dfrac{3}{4\pi}}\cdot \dfrac{z}{r},\; Y^{-1}_1(x,y,z)=\sqrt{\dfrac{3}{8\pi}}\cdot \dfrac{x-iy}{r} \label{4.4.44}\]

    and \[ Y^2_2=\sqrt{\dfrac{15}{32\pi}} \dfrac{(x+iy)^2}{r^2},\; Y^1_2=-\sqrt{\dfrac{15}{8\pi}}\dfrac{(x+iy)z}{r^2},\; Y^0_2=\sqrt{\dfrac{5}{16\pi}}\dfrac{(3z^2-1)}{r^2},\; etc. \label{4.4.45}\]

    The Y1m as a Basis of the l = 1 Subspace

    The \(Y^m_1\) are the \(l = 1\) eigenstates of \(L^2\) and \(L_z\). But what if we'd chosen to look for the common eigenstates of \(L^2\) and \(L_x\) instead? What \(l = 1\) state has zero angular momentum component in the direction of the x-axis? Clearly it will be \[ \sqrt{\dfrac{3}{4\pi}}\cdot\dfrac{x}{r} \label{4.4.46}\]

    , in other words the previous \(Y^0_1\) with z replaced by x, because after all, our labeling of axes was arbitrary.

    Now, \[ \sqrt{\dfrac{3}{4\pi}}\cdot\dfrac{x}{r} \; is \; just\; (1/\sqrt{2})(-Y^1_1+Y^{-1}_1). \label{4.4.47}\]

    In fact, any \(l = 1\) state, with a specified component in any direction, can be written as \[ \alpha_1|1,1\rangle+\alpha_0|1,0\rangle+\alpha_{-1}|1,-1\rangle=\sum \alpha_m|1,m\rangle. \label{4.4.48}\]

    This can be seen as follows: an \(l = 1\) state has to be linear in \(x/r,\; y/r,\; z/r\) (any quadratic term would give rise to \(e^{2i\phi}\) about an appropriate axis, call that the z-axis, so \(m = 2\) and \(l\) must be 2 or greater), and any such state can be written as a linear combination of \[ (x+iy)/\sqrt{2}r,\; (x-iy)/\sqrt{2}r,\; z/r \label{4.4.49}\]

    The bottom line, then, is that the \(Y^m_1\) do indeed provide a complete basis for the \(l = 1\) space of eigenstates of \(L^2\).

    Representing the Rotation Operator Within the l = 1 Subspace

    Recall that we originally introduced the angular momentum operator \(\vec{J}\) by defining it as the generator of infinitesimal rotations when acting on any wave function, including multicomponent wave functions. We found, using the commutativity properties of ordinary rotations, that the vector components of \(\vec{J}\) had to satisfy \([J_x,J_y]=i\hbar J_z\), etc., and from that we deduced the possible sets of eigenvalues of the commuting pair of operators \(\vec{J}^2\), \(J_z\) were \(j(j+1)\hbar^2\) for \(\vec{J}^2\), with \(j\) an integer of half an odd integer, and for each such \(j\) the allowed eigenvalues of \(J_z\) were

    \[ m\hbar,\; m=-j,-j+1,\dots,+j. \label{4.4.50}\]

    Back to the \(l = 1\) angular wave functions: we have established that any such function can be written

    \[ \alpha_1|1,1\rangle+\alpha_0|1,0\rangle+\alpha_{-1}|1,-1\rangle=\sum \alpha_m|1,m\rangle. \label{4.4.48b}\]

    and so is a vector in a three-dimensional space spanned by the set \(|1,m\rangle,\; m=1,0,-1\). In other words, the wave function is a three-component object. The angular momentum operator must therefore be a matrix operator in this three-dimensional space, such that, by definition, the effect of an infinitesimal rotation on the multicomponent wave function is:

    \[ R(\delta\vec{\theta})\Psi_{l=1}(\theta,\phi)=e^{-\dfrac{i}{\hbar}\delta\vec{\theta}\cdot\hat{\vec{J}}} \begin{pmatrix} \alpha_1 \\ \alpha_0 \\ \alpha_{-1} \end{pmatrix}=\begin{pmatrix} \alpha_1' \\ \alpha_0' \\ \alpha_{-1}' \end{pmatrix} \label{4.4.51}\]

    The unitary rotation operator acting in the \(l = 1\) subspace,

    \[ U(R(\vec{\theta}))=e^{-\dfrac{i\vec{\theta}\cdot\vec{J}}{\hbar}} \label{4.4.52}\]

    , has to be a \(3\times3\) matrix. The standard notation for its matrix elements is:

    \[ D^{(1)}_{m'm}(R(\vec{\theta}))=\langle 1,m'|e^{-\dfrac{i\vec{\theta}\cdot\vec{J}}{\hbar}}|1,m\rangle \label{4.4.53}\]

    so the rotated ket is

    \[ \alpha_{m'}'=\sum_{m',m}D^{(1)}_{m'm}\alpha_m,\;\; or\; \alpha'=D\alpha \label{4.4.54}\]

    To evaluate this matrix explicitly, we must expand the exponential and we need the matrix elements of \(J_z\), \(J_+\), \(J_-\) between the states \(|1,m\rangle\) -- which we already know.

    Now, the basis of the three-dimensional space is just the common eigenkets of \(\vec{J}^2\), \(J_z\), in this case identical to \(\vec{L}^2\), \(L_z\). We know the matrix elements of \(J_z\), \(J_+\), \(J_-\) between states \(|j,m\rangle\) from the earlier lecture, so it is simple to find the matrix representations of the components of \(J\) in this space:

    \[ J^{(1)}_x=\dfrac{\hbar}{\sqrt{2}}\begin{pmatrix} 0&1&0 \\ 1&0&1 \\ 0&1&0 \end{pmatrix},\;\; J^{(1)}_y=\dfrac{i\hbar}{\sqrt{2}}\begin{pmatrix} 0&-1&0 \\ 1&0&-1 \\ 0&1&0 \end{pmatrix},\;\; J^{(1)}_z=\hbar\begin{pmatrix} 1&0&0 \\ 0&0&0 \\ 0&0&-1 \end{pmatrix}. \label{4.4.55}\]

    We have added the superscript (1) because this representation of the infinitesimal rotation operators is specific to \(j = 1\) (representations for general values of \(j\) are as \((2j+1)\times(2j+1)\) matrices, reflecting the dimensionality of the space spanned by the \(2j + 1\) distinct \(m\) values).

    Expanding the exponential is not difficult, because by inspection \((J^{(1)}_z/\hbar)^3=(J^{(1)}_z/\hbar)\), so from spherical symmetry \((\hat{\vec{n}}\cdot\vec{J}^{(1)}/\hbar)^3=(\hat{\vec{n}}\vec{J}^{(1)}/\hbar)\) for a unit vector in any direction. The result is:

    \[ D^{(1)}(R(\vec{\theta}))=e^{-\dfrac{i\theta\hat{\vec{n}}\cdot\vec{J}}{\hbar}}=I+(\cos\theta-1)\left( \dfrac{\hat{\vec{n}}\cdot\vec{J}}{\hbar}\right)^2-i\sin\theta\left( \dfrac{\hat{\vec{n}}\cdot\vec{J}}{\hbar}\right). \label{4.4.56}\]

    One other point we should note: at the end of the linear algebra lecture, we discussed rotations about the z-axis in ordinary (x, y, z) space. Obviously, if we label a point in the (x, y) plane using the complex number \(x + iy\), a rotation by an angle \(\theta\) about the z-axis will move the point in such a way that the new label is \(e^{i\theta}(x+iy)\). The angle in this case has the opposite sign to that given by the operator above: the reason is that when we write the eigenstate as \[ -\sqrt{\dfrac{3}{8\pi}}\cdot \dfrac{x+iy}{r} \label{4.4.57}\]

    this is a function of position in the plane, not a point in the plane, so for the reasons discussed at the beginning of the first Angular Momentum lecture, the sign is opposite.

    This page titled 4.4: Orbital Eigenfunctions in 3-D is shared under a not declared license and was authored, remixed, and/or curated by Michael Fowler via source content that was edited to the style and standards of the LibreTexts platform; a detailed edit history is available upon request.