About the Project
14 Legendre and Related FunctionsApplications

§14.30 Spherical and Spheroidal Harmonics

Contents
  1. §14.30(i) Definitions
  2. §14.30(ii) Basic Properties
  3. §14.30(iii) Sums
  4. §14.30(iv) Applications

§14.30(i) Definitions

With l and m integers such that |m|\leq l, and \theta and \phi angles such that 0\leq\theta\leq\pi, 0\leq\phi\leq 2\pi,

14.30.2 Y_{l}^{m}\left(\theta,\phi\right)=\cos\left(m\phi\right)\mathsf{P}^{m}_{l}%
\left(\cos\theta\right)\text{ or }\sin\left(m\phi\right)\mathsf{P}^{m}_{l}%
\left(\cos\theta\right).

Y_{{l},{m}}\left(\theta,\phi\right) are known as spherical harmonics. Y_{l}^{m}\left(\theta,\phi\right) are known as surface harmonics of the first kind: tesseral for |m|<l and sectorial for |m|=l. Sometimes Y_{{l},{m}}\left(\theta,\phi\right) is denoted by i^{-l}\mathfrak{D}_{lm}(\theta,\phi); also the definition of Y_{{l},{m}}\left(\theta,\phi\right) can differ from (14.30.1), for example, by inclusion of a factor (-1)^{m}.

P^{m}_{n}\left(x\right) and Q^{m}_{n}\left(x\right) (x>1) are often referred to as the prolate spheroidal harmonics of the first and second kinds, respectively. P^{m}_{n}\left(ix\right) and Q^{m}_{n}\left(ix\right) (x>0) are known as oblate spheroidal harmonics of the first and second kinds, respectively. Segura and Gil (1999) introduced the scaled oblate spheroidal harmonics R_{n}^{m}(x)=e^{-i\pi n/2}P^{m}_{n}\left(ix\right) and T_{n}^{m}(x)=ie^{i\pi n/2}Q^{m}_{n}\left(ix\right) which are real when x>0 and n=0,1,2,\dots.

§14.30(ii) Basic Properties

Most mathematical properties of Y_{{l},{m}}\left(\theta,\phi\right) can be derived directly from (14.30.1) and the properties of the Ferrers function of the first kind given earlier in this chapter.

Special Values

14.30.4 Y_{{l},{m}}\left(0,\phi\right)=\begin{cases}\left(\dfrac{2l+1}{4\pi}\right)^{1%
/2},&m=0,\\
0,&|m|=1,2,3,\dots,\end{cases}

Symmetry

14.30.6 Y_{{l},{-m}}\left(\theta,\phi\right)=(-1)^{m}\overline{Y_{{l},{m}}\left(\theta%
,\phi\right)}.

Parity Operation

Orthogonality

14.30.8 \int_{0}^{2\pi}\!\!\int_{0}^{\pi}\overline{Y_{{l_{1}},{m_{1}}}\left(\theta,%
\phi\right)}Y_{{l_{2}},{m_{2}}}\left(\theta,\phi\right)\sin\theta\,\mathrm{d}%
\theta\,\mathrm{d}\phi=\delta_{l_{1},l_{2}}\delta_{m_{1},m_{2}}.

See also (34.3.22), and for further related integrals see Askey et al. (1986).

Herglotz generating function

The following is the Herglotz generating function

14.30.8_5 {\mathrm{e}}^{t{\mathbf{a}}\cdot{\mathbf{x}}}=\sqrt{4\pi}\sum_{n=0}^{\infty}%
\sum_{m=-n}^{n}\frac{t^{n}r^{n}\lambda^{m}Y_{{n},{m}}\left(\theta,\phi\right)}%
{\sqrt{(2n+1)(n+m)!(n-m)!}},

where \mathbf{a}=\left(\tfrac{1}{2\lambda}-\tfrac{\lambda}{2},-\tfrac{\mathrm{i}}{2%
\lambda}-\tfrac{\mathrm{i}\lambda}{2},1\right) and \mathbf{x}=\left(r\sin\theta\cos\phi,r\sin\theta\sin\phi,r\cos\theta\right).

§14.30(iii) Sums

Distributional Completeness

For a series representation of the product of two Dirac deltas in terms of products of spherical harmonics see §1.17(iii).

Addition Theorem

14.30.9 \mathsf{P}_{l}\left(\cos\theta_{1}\cos\theta_{2}+\sin\theta_{1}\sin\theta_{2}%
\cos\left(\phi_{1}-\phi_{2}\right)\right)=\frac{4\pi}{2l+1}\sum_{m=-l}^{l}%
\overline{Y_{{l},{m}}\left(\theta_{1},\phi_{1}\right)}Y_{{l},{m}}\left(\theta_%
{2},\phi_{2}\right).

See also (18.18.9) and (34.3.19).

§14.30(iv) Applications

In general, spherical harmonics are defined as the class of homogeneous harmonic polynomials. See Andrews et al. (1999, Chapter 9). The special class of spherical harmonics Y_{{l},{m}}\left(\theta,\phi\right), defined by (14.30.1), appear in many physical applications. As an example, Laplace’s equation \nabla^{2}W=0 in spherical coordinates (§1.5(ii)):

has solutions W(\rho,\theta,\phi)=\rho^{l}Y_{{l},{m}}\left(\theta,\phi\right), which are everywhere one-valued and continuous.

In the quantization of angular momentum the spherical harmonics Y_{{l},{m}}\left(\theta,\phi\right) are normalized solutions of the eigenvalue equations

14.30.11 \mathrm{L}^{2}Y_{{l},{m}}=\hbar^{2}l(l+1)Y_{{l},{m}},l=0,1,2,\dots,

and

14.30.11_5 \mathrm{L}_{z}Y_{{l},{m}}=\hbar mY_{{l},{m}},m=-l,-1+1,\dots,0,\dots,l-1,l,

where \hbar is the reduced Planck’s constant. Here, in spherical coordinates, \mathrm{L}^{2} is the squared angular momentum operator:

14.30.12 \mathrm{L}^{2}=-\hbar^{2}\left(\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}}\right),

and \mathrm{L}_{z} is the zcomponent of the angular momentum operator

14.30.13 \mathrm{L}_{z}=-\mathrm{i}\hbar\frac{\partial}{\partial\phi};

see Edmonds (1974, §2.5).

For applications in geophysics see Stacey (1977, §§4.2, 6.3, and 8.1).