Skip to main content
Mathematics LibreTexts

6.5: Laplace’s Equation and Spherical Symmetry

  • 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}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    We have seen that Laplace's equation, \(\nabla^{2} u=0\), arises in electrostatics as an equation for electric potential outside a charge distribution and it occurs as the equation governing equilibrium temperature distributions. As we had seen in the last chapter, Laplace’s equation generally occurs in the study of potential theory, which also includes the study of gravitational and fluid potentials. The equation is named after Pierre-Simon Laplace (1749-1827) who had studied the properties of this equation. Solutions of Laplace’s equation are called harmonic functions.

    Example \(\PageIndex{1}\)

    Solve Laplace’s equation in spherical coordinates.


    We seek solutions of this equation inside a sphere of radius \(r\) subject to the boundary condition as shown in Figure \(\PageIndex{1}\). The problem is given by Laplace’s equation Laplace’s equation in spherical coordinates\(^{1}\) \[\frac{1}{\rho^{2}} \frac{\partial}{\partial \rho}\left(\rho^{2} \frac{\partial u}{\partial \rho}\right)+\frac{1}{\rho^{2} \sin \theta} \frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial u}{\partial \theta}\right)+\frac{1}{\rho^{2} \sin ^{2} \theta} \frac{\partial^{2} u}{\partial \phi^{2}}=0,\label{eq:1}\] where \(u=u(\rho, \theta, \phi)\).

    Figure \(\PageIndex{1}\): A sphere of radius \(r\) with the boundary condition \(u(r, \theta, \phi)=\) \(g(\theta, \phi)\).

    The boundary conditions are given by \[u(r, \theta, \phi)=g(\theta, \phi), \quad 0<\phi<2 \pi, \quad 0<\theta<\pi,\nonumber \] and the periodic boundary conditions \[u(\rho, \theta, 0)=u(\rho, \theta, 2 \pi), \quad u_{\phi}(\rho, \theta, 0)=u_{\phi}(\rho, \theta, 2 \pi),\nonumber \] where \(0<\rho<\infty\), and \(0<\theta<\pi\).


    The Laplacian in spherical coordinates is given in Problem ?? in Chapter 8.

    Figure \(\PageIndex{2}\): Definition of spherical coordinates \((\rho, \theta, \phi)\). Note that there are different conventions for labeling spherical coordinates. This labeling is used often in physics.

    As before, we perform a separation of variables by seeking product solutions of the form \(u(\rho, \theta, \phi)=R(\rho) \Theta(\theta) \Phi(\phi)\). Inserting this form into the Laplace equation, we obtain \[\frac{\Theta \Phi}{\rho^{2}} \frac{d}{d \rho}\left(\rho^{2} \frac{d R}{d \rho}\right)+\frac{R \Phi}{\rho^{2} \sin \theta} \frac{d}{d \theta}\left(\sin \theta \frac{d \Theta}{d \theta}\right)+\frac{R \Theta}{\rho^{2} \sin ^{2} \theta} \frac{d^{2} \Phi}{d \phi^{2}}=0 .\label{eq:2}\] Multiplying this equation by \(\rho^{2}\) and dividing by \(R \Theta \Phi\), yields \[\frac{1}{R} \frac{d}{d \rho}\left(\rho^{2} \frac{d R}{d \rho}\right)+\frac{1}{\sin \theta \Theta} \frac{d}{d \theta}\left(\sin \theta \frac{d \Theta}{d \theta}\right)+\frac{1}{\sin ^{2} \theta \Phi} \frac{d^{2} \Phi}{d \phi^{2}}=0 .\label{eq:3}\]

    Note that the first term is the only term depending upon \(\rho\). Thus, we can separate out the radial part. However, there is still more work to do on the other two terms, which give the angular dependence. Thus, we have \[-\frac{1}{R} \frac{d}{d \rho}\left(\rho^{2} \frac{d R}{d \rho}\right)=\frac{1}{\sin \theta \Theta} \frac{d}{d \theta}\left(\sin \theta \frac{d \Theta}{d \theta}\right)+\frac{1}{\sin ^{2} \theta \Phi} \frac{d^{2} \Phi}{d \phi^{2}}=-\lambda,\label{eq:4}\] where we have introduced the first separation constant. This leads to two equations: \[\frac{d}{d \rho}\left(\rho^{2} \frac{d R}{d \rho}\right)-\lambda R=0\label{eq:5}\] and \[\frac{1}{\sin \theta \Theta} \frac{d}{d \theta}\left(\sin \theta \frac{d \Theta}{d \theta}\right)+\frac{1}{\sin ^{2} \theta \Phi} \frac{d^{2} \Phi}{d \phi^{2}}=-\lambda .\label{eq:6}\]

    The final separation can be performed by multiplying the last equation by \(\sin ^{2} \theta\), rearranging the terms, and introducing a second separation constant: \[\frac{\sin \theta}{\Theta} \frac{d}{d \theta}\left(\sin \theta \frac{d \Theta}{d \theta}\right)+\lambda \sin ^{2} \theta=-\frac{1}{\Phi} \frac{d^{2} \Phi}{d \phi^{2}}=\mu .\label{eq:7}\] From this expression we can determine the differential equations satisfied by \(\Theta(\theta)\) and \(\Phi(\phi)\) : \[\sin \theta \frac{d}{d \theta}\left(\sin \theta \frac{d \Theta}{d \theta}\right)+\left(\lambda \sin ^{2} \theta-\mu\right) \Theta=0\label{eq:8}\] and \[\frac{d^{2} \Phi}{d \phi^{2}}+\mu \Phi=0 .\label{eq:9}\]


    Equation \(\eqref{eq:6}\) is a key equation which occurs when studying problems possessing spherical symmetry. It is an eigenvalue problem for \(Y(\theta, \phi)=\Theta(\theta) \Phi(\phi)\), \(L Y=-\lambda Y\), where \[L=\frac{1}{\sin \theta} \frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right)+\frac{1}{\sin ^{2} \theta} \frac{\partial^{2}}{\partial \phi^{2}} .\nonumber \] The eigenfunctions of this operator are referred to as spherical harmonics.

    We now have three ordinary differential equations to solve. These are the radial equation \(\eqref{eq:5}\) and the two angular equations \(\eqref{eq:8}\)-\(\eqref{eq:9}\). We note that all three are in Sturm-Liouville form. We will solve each eigenvalue problem subject to appropriate boundary conditions.

    The simplest of these differential equations is Equation \(\eqref{eq:9}\) for \(\Phi(\phi)\). We have seen equations of this form many times and the general solution is a linear combination of sines and cosines. Furthermore, in this problem \(u(\rho, \theta, \phi)\) is periodic in \(\phi\), \[u(\rho, \theta, 0)=u(\rho, \theta, 2 \pi), \quad u_{\phi}(\rho, \theta, 0)=u_{\phi}(\rho, \theta, 2 \pi) .\nonumber \] Since these conditions hold for all \(\rho\) and \(\theta\), we must require that \(\Phi(\phi)\) satisfy the periodic boundary conditions \[\Phi(0)=\Phi(2 \pi), \quad \Phi^{\prime}(0)=\Phi^{\prime}(2 \pi) .\nonumber \] The eigenfunctions and eigenvalues for Equation \(\eqref{eq:9}\) are then found as \[\Phi(\phi)=\{\cos m \phi, \sin m \phi\}, \quad \mu=m^{2}, \quad m=0,1, \ldots .\label{eq:10}\]

    Next we turn to solving equation, \(\eqref{eq:8}\). We first transform this equation in order to identify the solutions. Let \(x=\cos \theta\). Then the derivatives with respect to \(\theta\) transform as \[\frac{d}{d \theta}=\frac{d x}{d \theta} \frac{d}{d x}=-\sin \theta \frac{d}{d x} .\nonumber \] Letting \(y(x)=\Theta(\theta)\) and noting that \(\sin ^{2} \theta=1-x^{2}\), Equation \(\eqref{eq:8}\) becomes \[\frac{d}{d x}\left(\left(1-x^{2}\right) \frac{d y}{d x}\right)+\left(\lambda-\frac{m^{2}}{1-x^{2}}\right) y=0 .\label{eq:11}\] We further note that \(x \in[-1,1]\), as can be easily confirmed by the reader.

    This is a Sturm-Liouville eigenvalue problem. The solutions consist of a set of orthogonal eigenfunctions. For the special case that \(m=0\) Equation \(\eqref{eq:11}\) becomes \[\frac{d}{d x}\left(\left(1-x^{2}\right) \frac{d y}{d x}\right)+\lambda y=0 \text {. }\label{eq:12}\] In a course in differential equations one learns to seek solutions of this equation in the form \[y(x)=\sum_{n=0}^{\infty} a_{n} x^{n} .\nonumber \] This leads to the recursion relation \[a_{n+2}=\frac{n(n+1)-\lambda}{(n+2)(n+1)} a_{n} .\nonumber \] Setting \(n=0\) and seeking a series solution, one finds that the resulting series does not converge for \(x=\pm 1\). This is remedied by choosing \(\lambda=\ell(\ell+1)\) for \(\ell=0,1, \ldots\), leading to the differential equation \[\frac{d}{d x}\left(\left(1-x^{2}\right) \frac{d y}{d x}\right)+\ell(\ell+1) y=0 .\label{eq:13}\] We saw this equation in Chapter 5 in the form \[\left(1-x^{2}\right) y^{\prime \prime}-2 x y^{\prime}+\ell(\ell+1) y=0 .\nonumber \] The solutions of this differential equation are Legendre polynomials, denoted by \(P_{\ell}(x)\)

    For the more general case, \(m \neq 0\), the differential equation \(\eqref{eq:11}\) with \(\lambda=\ell(\ell+1)\) becomes \[\frac{d}{d x}\left(\left(1-x^{2}\right) \frac{d y}{d x}\right)+\left(\ell(\ell+1)-\frac{m^{2}}{1-x^{2}}\right) y=0 .\label{eq:14}\] The solutions of this equation are called the associated Legendre functions. The two linearly independent solutions are denoted by \(P_{\ell}^{m}(x)\) and \(Q_{\ell}^{m}(x)\). The latter functions are not well behaved at \(x=\pm 1\), corresponding to the north and south poles of the original problem. So, we can throw out these solutions in many physical cases, leaving \[\Theta(\theta)=P_{\ell}^{m}(\cos \theta)\nonumber \] as the needed solutions. In Table \(6.5\) we list a few of these.

    Table \(\PageIndex{1}\):Associated Legendre Functions, \(P_{n}^{m}(x)\).
    \(P_{n}^{m}(x)\) \(P_{n}^{m}(\cos \theta)\)
    \(P_{0}^{0}(x)\) 1 1
    \(P_{1}^{0}(x)\) \(x\) \(\cos \theta\)
    \(P_{1}^{1}(x)\) \(-\left(1-x^{2}\right)^{\frac{1}{2}}\) \(-\sin \theta\)
    \(P_{2}^{0}(x)\) \(\frac{1}{2}\left(3 x^{2}-1\right)\) \(\frac{1}{2}\left(3 \cos ^{2} \theta-1\right)\)
    \(P_{2}^{1}(x)\) \(-3 x\left(1-x^{2}\right)^{\frac{1}{2}}\) \(-3 \cos \theta \sin \theta\)
    \(P_{2}^{2}(x)\) \(3\left(1-x^{2}\right)\) \(3 \sin ^{2} \theta\)
    \(P_{3}^{0}(x)\) \(\frac{1}{2}\left(5 x^{3}-3 x\right)\) \(\frac{1}{2}\left(5 \cos ^{3} \theta-3 \cos \theta\right)\)
    \(P_{3}^{1}(x)\) \(-\frac{3}{2}\left(5 x^{2}-1\right)\left(1-x^{2}\right)^{\frac{1}{2}}\) \(-\frac{3}{2}\left(5 \cos ^{2} \theta-1\right) \sin \theta\)
    \(P_{3}^{2}(x)\) \(15 x\left(1-x^{2}\right)\) \(15 \cos \theta \sin ^{2} \theta\)
    \(P_{3}^{3}(x)\) \(-15\left(1-x^{2}\right)^{\frac{3}{2}}\) \(-15 \sin ^{3} \theta\)

    The associated Legendre functions are related to the Legendre polynomials by\(^{2}\) \[P_{\ell}^{m}(x)=(-1)^{m}\left(1-x^{2}\right)^{m / 2} \frac{d^{m}}{d x^{m}} P_{\ell}(x),\label{eq:15}\] for \(\ell=0,1,2, \ldots\) and \(m=0,1, \ldots, \ell\). We further note that \(P_{\ell}^{0}(x)=P_{\ell}(x)\), as one can see in the table. Since \(P_{\ell}(x)\) is a polynomial of degree \(\ell\), then for \(m>\ell, \frac{d^{m}}{d x^{m}} P_{\ell}(x)=0\) and \(P_{\ell}^{m}(x)=0\).


    The factor of \((-1)^{m}\) is known as the Condon-Shortley phase and is useful in quantum mechanics in the treatment of agular momentum. It is sometimes omitted by some

    Furthermore, since the differential equation only depends on \(m^{2}, P_{\ell}^{-m}(x)\) is proportional to \(P_{\ell}^{m}(x)\). One normalization is given by \[P_{\ell}^{-m}(x)=(-1)^{m} \frac{(\ell-m) !}{(\ell+m) !} P_{\ell}^{m}(x) .\nonumber \]

    The associated Legendre functions also satisfy the orthogonality condition \[\int_{-1}^{1} P_{\ell}^{m}(x) P_{\ell^{\prime}}^{m}(x) d x=\frac{2}{2 \ell+1} \frac{(\ell+m) !}{(\ell-m) !} \delta_{\ell \ell^{\prime}} .\label{eq:16}\]

    The last differential equation we need to solve is the radial equation. With \(\lambda=\ell(\ell+1), \ell=0,1,2, \ldots\), the radial equation \(\eqref{eq:5}\) can be written as \[\rho^{2} R^{\prime \prime}+2 \rho R^{\prime}-\ell(\ell+1) R=0 .\label{eq:17}\]

    The radial equation is a Cauchy-Euler type of equation. So, we can guess the form of the solution to be \(R(\rho)=\rho^{s}\), where \(s\) is a yet to be determined constant. Inserting this guess into the radial equation, we obtain the characteristic equation \[s(s+1)=\ell(\ell+1) .\nonumber \] Solving for \(s\), we have \[s=\ell,-(\ell+1) .\nonumber \] Thus, the general solution of the radial equation is \[R(\rho)=a \rho^{\ell}+b \rho^{-(\ell+1)} .\label{eq:18}\]

    We would normally apply boundary conditions at this point. The boundary condition \(u(r, \theta, \phi)=g(\theta, \phi)\) is not a homogeneous boundary condition, so we will need to hold off using it until we have the general solution to the three dimensional problem. However, we do have a hidden condition. Since we are interested in solutions inside the sphere, we need to consider what happens at \(\rho=0\). Note that \(\rho^{-(\ell+1)}\) is not defined at the origin. Since the solution is expected to be bounded at the origin, we can set \(b=0\). So, in the current problem we have established that \[R(\rho)=a \rho^{\ell} .\nonumber \]


    When seeking solutions outside the sphere, one considers the boundary condition \(R(\rho) \rightarrow 0\) as \(\rho \rightarrow \infty\). In this case, \(R(\rho)=\rho^{-(\ell+1)} .\)

    We have carried out the full separation of Laplace’s equation in spherical coordinates. The product solutions consist of the forms \[u(\rho, \theta, \phi)=\rho^{\ell} P_{\ell}^{m}(\cos \theta) \cos m \phi\nonumber \] and \[u(\rho, \theta, \phi)=\rho^{\ell} P_{\ell}^{m}(\cos \theta) \sin m \phi\nonumber \] for \(\ell=0,1,2, \ldots\) and \(m=0, \pm 1, \ldots, \pm \ell\). These solutions can be combined to give a complex representation of the product solutions as \[u(\rho, \theta, \phi)=\rho^{\ell} P_{\ell}^{m}(\cos \theta) e^{i m \phi} .\nonumber \] The general solution is then given as a linear combination of these product can be rewritten as solutions. As there are two indices, we have a double sum:\(^{3}\) \[u(\rho, \theta, \phi)=\sum_{\ell=0}^{\infty} \sum_{m=-\ell}^{\ell} a_{\ell m} \rho^{\ell} P_{\ell}^{m}(\cos \theta) e^{i m \phi} .\label{eq:19}\]


    While this appears to be a complex-valued solution, it can be rewritten as a sum over real functions. The inner sum contains terms for both \(m=k\) and to give a complex representation of the product solutions as \(m=-k\). Adding these contributions, we have that \[a_{\ell k} \rho^{\ell} P_{\ell}^{k}(\cos \theta) e^{i k \phi}+a_{\ell(-k)} \rho^{\ell} P_{\ell}^{-k}(\cos \theta) e^{-i k \phi}\nonumber\] can be rewritten as \[\left(A_{\ell k} \cos k \phi+B_{\ell k} \sin k \phi\right) \rho^{\ell} P_{\ell}^{k}(\cos \theta) .\nonumber \]

    Example \(\PageIndex{2}\): Laplace's Equation with Azimuthal Symmetry

    As a simple example we consider the solution of Laplace’s equation in which there is azimuthal symmetry. Let \[u(r, \theta, \phi)=g(\theta)=1-\cos 2 \theta \text {. }\nonumber \] This function is zero at the poles and has a maximum at the equator. So, this could be a crude model of the temperature distribution of the Earth with zero temperature at the poles and a maximum near the equator.

    Figure \(\PageIndex{3}\): A sphere of radius \(r\) with the boundary condition \[u(r,\theta ,\phi )=1-\cos 2\theta.\nonumber\]

    In problems in which there is no \(\phi\)-dependence, only the \(m=0\) terms of the general solution survives. Thus, we have that \[u(\rho, \theta, \phi)=\sum_{\ell=0}^{\infty} a_{\ell} \rho^{\ell} P_{\ell}(\cos \theta) .\label{eq:20}\] Here we have used the fact that \(P_{\ell}^{0}(x)=P_{\ell}(x)\). We just need to determine the unknown expansion coefficients, \(a_{\ell}\). Imposing the boundary condition at \(\rho=r\), we are lead to \[g(\theta)=\sum_{\ell=0}^{\infty} a_{\ell} r^{\ell} P_{\ell}(\cos \theta)\label{eq:21}\]

    This is a Fourier-Legendre series representation of \(g(\theta)\). Since the Legendre polynomials are an orthogonal set of eigenfunctions, we can extract the coefficients.

    In Chapter 5 we had proven that \[\int_{0}^{\pi} P_{n}(\cos \theta) P_{m}(\cos \theta) \sin \theta d \theta=\int_{-1}^{1} P_{n}(x) P_{m}(x) d x=\frac{2}{2 n+1} \delta_{n m} .\nonumber \] So, multiplying the expression for \(g(\theta)\) by \(P_{m}(\cos \theta) \sin \theta\) and integrating, we obtain the expansion coefficients: \[a_{\ell}=\frac{2 \ell+1}{2 r^{\ell}} \int_{0}^{\pi} g(\theta) P_{\ell}(\cos \theta) \sin \theta d \theta \text {. }\label{eq:22}\]

    Sometimes it is easier to rewrite \(g(\theta)\) as a polynomial in \(\cos \theta\) and avoid the integration. For this example we see that \[\begin{align} g(\theta) &=1-\cos 2 \theta\nonumber \\ &=2 \sin ^{2} \theta\nonumber \\ &=2-2 \cos ^{2} \theta .\label{eq:23} \end{align}\] Thus, setting \(x=\cos \theta\) and \(G(x)=g(\theta(x))\), we have \(G(x)=2-2 x^{2}\).

    We seek the form \[G(x)=c_{0} P_{0}(x)+c_{1} P_{1}(x)+c_{2} P_{2}(x),\nonumber \] where \(P_{0}(x)=1, P_{1}(x)=x\), and \(P_{2}(x)=\frac{1}{2}\left(3 x^{2}-1\right)\). Since \(G(x)=2-2 x^{2}\) does not have any \(x\) terms, we know that \(c_{1}=0\). So, \[2-2 x^{2}=c_{0}(1)+c_{2} \frac{1}{2}\left(3 x^{2}-1\right)=c_{0}-\frac{1}{2} c_{2}+\frac{3}{2} c_{2} x^{2} .\nonumber \]

    By observation we have \(c_{2}=-\frac{4}{3}\) and thus, \(c_{0}=2+\frac{1}{2} c_{2}=\frac{4}{3}\). Therefore, \(G(x)=\frac{4}{3} P_{0}(x)-\frac{4}{3} P_{2}(x) .\)

    We have found the expansion of \(g(\theta)\) in terms of Legendre polynomials, \[g(\theta)=\frac{4}{3} P_{0}(\cos \theta)-\frac{4}{3} P_{2}(\cos \theta) .\label{eq:24}\] Therefore, the nonzero coefficients in the general solution become \[a_{0}=\frac{4}{3}, \quad a_{2}=\frac{4}{3} \frac{1}{r^{2}},\nonumber \] and the rest of the coefficients are zero. Inserting these into the general solution, we have the final solution \[\begin{align} u(\rho, \theta, \phi) &=\frac{4}{3} P_{0}(\cos \theta)-\frac{4}{3}\left(\frac{\rho}{r}\right)^{2} P_{2}(\cos \theta)\nonumber \\ &=\frac{4}{3}-\frac{2}{3}\left(\frac{\rho}{r}\right)^{2}\left(3 \cos ^{2} \theta-1\right)\label{eq:25} \end{align}\]

    Spherical Harmonics

    The solutions of the angular parts of the problem are often combined into one function of two variables, as problems with spherical symmetry arise often, leaving the main differences between such problems confined to the radial equation. These functions are referred to as spherical harmonics, \(Y_{\ell m}(\theta, \phi)\), which are defined with a special normalization as \[Y_{\ell m}(\theta, \phi)=(-1)^{m} \sqrt{\frac{2 \ell+1}{4 \pi} \frac{(\ell-m) !}{(\ell+m) !}} P_{\ell}^{m}(\cos \theta) e^{i m \phi} .\nonumber \] These satisfy the simple orthogonality relation \[\int_{0}^{\pi} \int_{0}^{2 \pi} Y_{\ell m}(\theta, \phi) Y_{\ell^{\prime} m^{\prime}}^{*}(\theta, \phi) \sin \theta d \phi d \theta=\delta_{\ell \ell^{\prime}} \delta_{m m^{\prime}} .\nonumber\]


    \(Y_{\ell m}(\theta, \phi)\), are the spherical harmonics. Spherical harmonics are important in applications from atomic electron configurations to gravitational fields, planetary magnetic fields, and the cosmic microwave background radiation.

    As seen earlier in the chapter, the spherical harmonics are eigenfunctions of the eigenvalue problem \(L Y=-\lambda Y\), where \[L=\frac{1}{\sin \theta} \frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right)+\frac{1}{\sin ^{2} \theta} \frac{\partial^{2}}{\partial \phi^{2}} .\nonumber \] This operator appears in many problems in which there is spherical symmetry, such as obtaining the solution of Schrödinger’s equation for the hydrogen atom as we will see later. Therefore, it is customary to plot spherical harmonics. Because the \(Y_{\ell m}\) ’s are complex functions, one typically plots either the real part or the modulus squared. One rendition of \(\left|Y_{\ell m}(\theta, \phi)\right|^{2}\) is shown in Table \(\PageIndex{2}\) for \(\ell, m=0,1,2,3\).

    Table \(\PageIndex{2}\): The first few spherical harmonics, \(\left|Y_{\ell m}(\theta, \phi)\right|^{2}\)


    We could also look for the nodal curves of the spherical harmonics like we had for vibrating membranes. Such surface plots on a sphere are shown in Table \(\PageIndex{3}\). The colors provide for the amplitude of the \(|Y_{\ell m}(θ, φ)|^2\). We can match these with the shapes in Table \(\PageIndex{2}\) by coloring the plots with some of the same colors as shown in Table \(\PageIndex{3}\). However, by plotting just the sign of the spherical harmonics, as in Table \(\PageIndex{4}\), we can pick out the nodal curves much easier.

    Table \(\PageIndex{3}\): Spherical harmonic contours for \(\left|Y_{\ell m}(\theta, \phi)\right|^{2}\).


    Table \(\PageIndex{4}\): In these figures we show the nodal curves of \(\left|Y_{\ell m}(\theta, \phi)\right|^{2}\) Along the first column \((m=0)\) are the zonal harmonics seen as \(\ell\) horizontal circles. Along the top diagonal \((m=\ell)\) are the sectional harmonics. These look like orange sections formed from \(m\) vertical circles. The remaining harmonics are tesseral harmonics. They look like a checkerboard pattern formed from intersections of \(\ell-m\) horizontal circles and \(m\) vertical circles.


    Spherical, or surface, harmonics can be further grouped into zonal, sectoral, and tesseral harmonics. Zonal harmonics correspond to the \(m = 0\) modes. In this case, one seeks nodal curves for which \(P_{\ell} (\cos θ) = 0\). Solutions of this equation lead to constant \(θ\) values such that \(\cos θ\) is a zero of the Legendre polynomial, \(P_{\ell}(x)\). The zonal harmonics correspond to the first column in Table \(\PageIndex{4}\). Since \(P_{\ell}(x)\) is a polynomial of degree \(\ell\), the zonal harmonics consist of \(\ell\) latitudinal circles.

    Sectoral, or meridional, harmonics result for the case that \(m=\pm \ell\). For this case, we note that \(P_{\ell}^{\pm \ell}(x) \propto\left(1-x^{2}\right)^{m / 2}\). This function vanishes for \(x=\pm 1\), or \(\theta=0, \pi\). Therefore, the spherical harmonics can only produce nodal curves for \(e^{i m \phi}=0\). Thus, one obtains the meridians satisfying the condition \(A \cos m \phi+B \sin m \phi=0\). Solutions of this equation are of the form \(\phi=\) constant. These modes can be seen in Table \(\PageIndex{4}\) in the top diagonal and can be described as \(m\) circles passing through the poles, or longitudinal circles.

    Tesseral harmonics consist of the rest of the modes, which typically look like a checker board glued to the surface of a sphere. Examples can be seen in the pictures of nodal curves, such as Table \(\PageIndex{4}\). Looking in Table \(\PageIndex{4}\) along the diagonals going downward from left to right, one can see the same number of latitudinal circles. In fact, there are \(\ell-m\) latitudinal nodal curves in these figures.

    In summary, the spherical harmonics have several representations, as show in Tables \(\PageIndex{3}\)-\(\PageIndex{4}\). Note that there are \(\ell\) nodal lines, \(m\) meridional curves, and \(\ell − m\) horizontal curves in these figures. The plots in Table \(\PageIndex{2}\) are the typical plots shown in physics for discussion of the wavefunctions of the hydrogen atom. Those in \(\PageIndex{3}\) are useful for describing gravitational or electric potential functions, temperature distributions, or wave modes on a spherical surface. The relationships between these pictures and the nodal curves can be better understood by comparing respective plots. Several modes were separated out in Figures \(\PageIndex{4}\)-\(\PageIndex{9}\) to make this comparison easier.

    Figure \(\PageIndex{4}\): Zonal harmonics, \(\ell=1\), \(m=0\).
    Figure \(\PageIndex{5}\): Zonal harmonics, \(\ell=2\), \(m=0\).
    Figure \(\PageIndex{6}\): Sectoral harmonics, \(\ell=2\), \(m=2\).
    Figure \(\PageIndex{7}\): Tesseral harmonics, \(\ell=3\), \(m=1\).
    Figure \(\PageIndex{8}\): Sectoral harmonics, \(\ell=3\), \(m=3\).
    Figure \(\PageIndex{9}\): Tesseral harmonics, \(\ell=4\), \(m=3\).

    This page titled 6.5: Laplace’s Equation and Spherical Symmetry is shared under a CC BY-NC-SA 3.0 license and was authored, remixed, and/or curated by Russell Herman via source content that was edited to the style and standards of the LibreTexts platform.