Spherical harmonics: Definition from Answers.com

In mathematics, the spherical harmonics are the angular portion of a set of solutions to Laplace's equation. Represented in a system of spherical coordinates, Laplace's spherical harmonics $Y_\ell^m$ are a specific set of spherical harmonics that forms an orthogonal system, first introduced by Pierre Simon de Laplace.^[1] Spherical harmonics are important in many theoretical and practical applications, particularly in the computation of atomic orbital electron configurations, representation of gravitational fields, geoids, and the magnetic fields of planetary bodies and stars, and characterization of the cosmic microwave background radiation. In 3D computer graphics, spherical harmonics play a special role in a wide variety of topics including indirect lighting (ambient occlusion, global illumination, precomputed radiance transfer, etc.) and in recognition of 3D shapes.

History

Spherical harmonics were first investigated in connection with the Newtonian potential of Newton's law of universal gravitation in three-dimensions. In 1782, Pierre-Simon de Laplace had, in his Mécanique Céleste, determined that the gravitational potential associated to a set of point-masses m_i located at points x_i was given by

$P(\mathbf{x}) = \sum_i \frac{m_i}{|\mathbf{x}_i - \mathbf{x}|}.$

Each term in the above summation is an individual Newtonian potential for a point mass. Just prior to that time, Adrien-Marie Legendre had investigated the expansion of the Newtonian potential in powers of r = |x| and r₁ = |x₁|. He discovered that

$\frac{1}{|\mathbf{x}_1 - \mathbf{x}|} = P_0(\cos\gamma)\frac{1}{r_1} + P_1(\cos\gamma)\frac{r}{r_1^2} + P_2(\cos\gamma)\frac{r^2}{r_1^3}+\cdots$

where γ is the angle between the vectors x and x₁. The coefficients P_i are the Legendre polynomials, and they are a special case of spherical harmonics. Subsequently, in his 1782 memoire, Laplace investigated these coefficients using spherical coordinates to represent the angle γ between x₁ and x.

In 1867, William Thomson (Lord Kelvin) and Peter Guthrie Tait introduced the solid spherical harmonics in their Treatise on Natural Philosophy, and also first introduced the name of "spherical harmonics" for these functions. The solid harmonics were solutions of Laplace's equation

$\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} = 0$

homogeneous of degree ℓ. By examining Laplace's equation in spherical coordinates, Thomson and Tait recovered Laplace's spherical harmonics. The term "Laplace's coefficients" was employed by William Whewell to describe the particular system of solutions introduced along these lines, whereas others reserved this designation for the zonal spherical harmonics that had properly been introduced by Laplace and Legendre.

The 19th century development of Fourier series made possible the solution of a wide variety of physical problems in rectangular domains, such as the solution of the heat equation and wave equation. This could be achieved by expansion of functions in series of trigonometric functions. Whereas the trigonometric functions in a Fourier series represent the fundamental modes of vibration in a string, the spherical harmonics represent the fundamental modes of vibration of a sphere in much the same way. Many aspects of the theory of Fourier series could be generalized by taking expansions in spherical harmonics rather than trigonometric functions. This was a boon for problems possessing spherical symmetry, such as those of celestial mechanics originally studied by Laplace and Legendre.

The prevalence of spherical harmonics already in physics set the stage for their later importance in the 20th century birth of quantum mechanics. The spherical harmonics are eigenfunctions of the square of the orbital angular momentum operator

$-i\hbar\mathbf{r}\times\nabla,$

and therefore they represent the different quantized configurations of atomic orbitals.

Laplace's spherical harmonics

Real (Laplace) spherical harmonics $Y_{\ell m}$ , for $\ell=0$ to 4 (top to bottom) and

m = 0

to 4 (left to right). The negative order harmonics $Y_{\ell,-m}$ are rotated about the

z

axis by $90^\circ/m$ with respect to the positive order ones.

Laplace's equation in spherical coordinates is:

$\nabla^2 f = {1 \over r^2}{\partial \over \partial r}\left(r^2 {\partial f \over \partial r}\right) + {1 \over r^2\sin\theta}{\partial \over \partial \theta}\left(\sin\theta {\partial f \over \partial \theta}\right) + {1 \over r^2\sin^2\theta}{\partial^2 f \over \partial \varphi^2} = 0$

(see also del in cylindrical and spherical coordinates). Consider the problem of finding solutions of the form ƒ(r,θ,φ) = r^ℓΘ(θ)Φ(φ), where ℓ is a non-negative integer.^[2] For ƒ of this special form, Laplace's equation reduces to

${\Phi(\varphi) \over \sin\theta}{d \over d\theta}\left(\sin\theta {d\Theta \over d\theta} \right) + {\Theta(\theta) \over \sin^2 \theta}{d^2\Phi \over d\varphi^2} + \ell(\ell+1)\Theta(\theta)\Phi(\varphi) = 0.$

Using the technique of separation of variables, two differential equations result:

$\frac{1}{\Phi(\varphi)} \frac{d^2 \Phi(\varphi)}{d\varphi^2} = -m^2$

$\ell(\ell+1)\sin ^2(\theta) + \frac{\sin(\theta)}{\Theta(\theta)} \frac{d}{d\theta} \left [ \sin(\theta) \frac{d\Theta}{d\theta} \right ] = m^2$

for some m. Because Φ must be a periodic function whose period evenly divides 2π, m is necessarily an integer. Hence, the angular solutions can be shown to be a product of trigonometric functions, here represented as a complex exponential, and associated Legendre functions:

$Y_\ell^m (\theta, \varphi ) = N \, e^{i m \varphi } \, P_\ell^m (\cos{\theta} ),$

which fulfill

$\nabla^2 Y_\ell^m (\theta, \varphi ) = -\frac{\ell (\ell + 1 )}{r^2} Y_\ell^m (\theta, \varphi ).$

Here $Y_\ell^m$ is called a spherical harmonic function of degree ℓ and order m, $P_\ell^m$ is an associated Legendre function, N is a normalization constant, and θ and φ represent colatitude and longitude, respectively. In particular, the colatitude θ, or polar angle, ranges from 0 at the North Pole to π at the South Pole, assuming the value of π/2 at the Equator, and the longitude φ, or azimuth, may assume all values with 0 ≤ φ < 2π.

When Laplace's equation is solved on the surface of the sphere, the periodic boundary conditions in φ, as well as regularity conditions at both the north and south poles, ensure that the degree ℓ and order m are integers that satisfy ℓ ≥ 0 and |m| ≤ ℓ. In contrast, if the function ƒ were only to have been defined for θ ≤ θ₀, then the resulting spherical cap harmonics would have been defined for integer order, but non-integer degree. The general solution to Laplace's equation in a ball centered at the origin is a linear combination of the spherical harmonic functions multiplied by the appropriate scale factor r^ℓ,

$f(r, \theta, \varphi) = \sum_{\ell=0}^\infty \sum_{m=-\ell}^\ell f_\ell^m \, r^\ell \, Y_\ell^m (\theta, \varphi ),$

where the $f_\ell^m$ are constants and the factors $r^\ell \, Y_\ell^m$ are known as solid harmonics. Such an expansion is valid in the ball

$r < R = 1/\limsup_{\ell\to\infty} |f_\ell^m|^{1/\ell}.$

Conventions

Orthogonality and normalization

Several different normalizations are in common use for the Laplace spherical harmonic functions. In physics and seismology, these functions are generally defined as

$Y_\ell^m( \theta , \varphi ) = \sqrt{{(2\ell+1)\over 4\pi}{(\ell-m)!\over (\ell+m)!}} \, P_\ell^m ( \cos{\theta} ) \, e^{i m \varphi }$

which are orthonormal

$\int_{\theta=0}^\pi\int_{\varphi=0}^{2\pi}Y_\ell^m \, Y_{\ell'}^{m'*} \, d\Omega=\delta_{\ell\ell'}\, \delta_{mm'},$

where δ_aa = 1, δ_ab = 0 if a ≠ b, (see Kronecker delta) and dΩ = sinθ dφ dθ. The disciplines of geodesy and spectral analysis use

$Y_\ell^m( \theta , \varphi ) = \sqrt{{(2\ell+1) }{(\ell-m)!\over (\ell+m)!}} \, P_\ell^m ( \cos{\theta} )\, e^{i m \varphi }$

which possess unit power

${1 \over 4 \pi} \int_{\theta=0}^\pi\int_{\varphi=0}^{2\pi}Y_\ell^m \, Y_{\ell'}^{m'*} d\Omega=\delta_{\ell\ell'}\, \delta_{mm'}.$

The magnetics community, in contrast, uses Schmidt semi-normalized harmonics

$Y_\ell^m( \theta , \varphi ) = \sqrt{{(\ell-m)!\over (\ell+m)!}} \, P_\ell^m ( \cos{\theta} ) \, e^{i m \varphi }$

which have the normalization

$\int_{\theta=0}^\pi\int_{\varphi=0}^{2\pi}Y_\ell^m \, Y_{\ell'}^{m'*}d\Omega={4 \pi \over (2 \ell + 1)}\delta_{\ell\ell'}\, \delta_{mm'}.$

In quantum mechanics this normalization is often used, too, and is there named Racah's normalization after Giulio Racah.

Using the identity (see associated Legendre functions)

$P_\ell ^{-m} = (-1)^m \frac{(\ell-m)!}{(\ell+m)!} P_\ell ^{m}$

it can be shown that all of the above normalized spherical harmonic functions satisfy

$Y_\ell^{m*} (\theta, \varphi) = (-1)^m Y_\ell^{-m} (\theta, \varphi),$

where the superscript * denotes complex conjugation. Alternatively, this equation follows from the relation of the spherical harmonic functions with the Wigner D-matrix.

Condon-Shortley phase

One source of confusion with the definition of the spherical harmonic functions concerns a phase factor of (-1)^m, commonly referred to as the Condon-Shortley phase in the quantum mechanical literature. In the quantum mechanics community, it is common practice to either include this phase factor in the definition of the associated Legendre functions, or to append it to the definition of the spherical harmonic functions. There is no requirement to use the Condon-Shortley phase in the definition of the spherical harmonic functions, but including it can simplify some quantum mechanical operations, especially the application of raising and lowering operators. The geodesy and magnetics communities never include the Condon-Shortley phase factor in their definitions of the spherical harmonic functions.

Real form

A real basis of spherical harmonics is given by setting instead

$Y_{\ell m} = \begin{cases} Y_\ell^0 & \mbox{if } m=0\\ {1\over\sqrt2}\left(Y_\ell^m+(-1)^m \, Y_\ell^{-m}\right) = \sqrt{2} N_{(\ell,m)} P_\ell^m(\cos \theta) \cos m\varphi & \mbox{if } m>0 \\ {1\over i\sqrt2}\left(Y_\ell^{-m}-(-1)^{m}\, Y_\ell^{m}\right) = \sqrt{2} N_{(\ell,m)} P_\ell^{-m}(\cos \theta) \sin m\varphi &\mbox{if } m<0. \end{cases}$

where $N_{(\ell,m)}$ denotes the normalization constant as a function of ℓ and $m$ . The real form requires only associated Legendre functions $P_\ell^{|m|}$ of non-negative |m|. The harmonics with m > 0 are said to be of cosine type, and those with m < 0 of sine type. These real spherical harmonics are sometimes known as tesseral spherical harmonics.^[3] These functions have the same normalization properties as the complex ones above. See here for a list of real spherical harmonics up to and including $\ell = 5$ . Note, however, that the listed functions differ by the phase (-1)^m from the phase given in this article.

Spherical harmonics expansion

The Laplace spherical harmonics form a complete set of orthonormal functions and thus form an orthonormal basis of the Hilbert space of square-integrable functions. On the unit sphere, any square-integrable function can thus be expanded as a linear combination of these:

$f(\theta,\varphi)=\sum_{\ell=0}^{\infty} \sum_{m=-\ell}^\ell f_\ell^m \, Y_\ell^m(\theta,\varphi).$

This expansion holds in the sense of mean-square convergence — convergence in L² of the sphere — which is to say that

$\lim_{N\to\infty} \int_0^{2\pi}\int_0^\pi \left|f(\theta,\varphi)-\sum_{\ell=0}^N f_\ell^m Y_\ell^m(\theta,\varphi)\right|^2\sin\theta\, d\theta d\phi = 0.$

The expansion coefficients are the analogs of Fourier coefficients, and can be obtained by multiplying the above equation by the complex conjugate of a spherical harmonic, integrating over the solid angle $\Omega\!\,$ , and utilizing the above orthogonality relationships. This is justified rigorously by basic Hilbert space theory. For the case of orthonormalized harmonics, this gives:

$f_\ell^m=\int_{\Omega} f(\theta,\varphi)\, Y_\ell^{m*}(\theta,\varphi)d\Omega = \int_0^{2\pi}d\varphi\int_0^{\pi}d\theta\sin\theta f(\theta,\varphi)Y_\ell^{m*} (\theta,\varphi).$

If the coefficients decay in ℓ sufficiently rapidly — for instance, exponentially — then the series also converges uniformly to ƒ.

A real square-integrable function ƒ can be expanded in terms of the real harmonics Y_ℓm above as a sum

$f(\theta, \varphi) = \sum_{\ell=0}^\infty \sum_{m=-\ell}^\ell f_{\ell m} \, Y_{\ell m}(\theta, \varphi).$

Convergence of the series holds again in the same sense.

Spectrum analysis

Power spectrum in signal processing

The total power of a function ƒ is defined in the signal processing literature as the integral of the function squared, divided by the area of its domain. Using the orthonormality properties of the real unit-power spherical harmonic functions, it is straightforward to verify that the total power of a function defined on the unit sphere is related to its spectral coefficients by a generalization of Parseval's theorem:

$\frac{1}{4 \, \pi} \int_\Omega f(\Omega)^2\, d\Omega = \sum_{\ell=0}^\infty S_{f\!f}(\ell),$

where

$S_{f\!f}(\ell) = \sum_{m=-\ell}^\ell f_{\ell m}^2$

is defined as the angular power spectrum. In a similar manner, one can define the cross-power of two functions as

$\frac{1}{4 \, \pi} \int_\Omega f(\Omega) \, g(\Omega) \, d\Omega = \sum_{\ell=0}^\infty S_{fg}(\ell),$

where

$S_{fg}(\ell) = \sum_{m=-\ell}^\ell f_{\ell m} g_{\ell m}$

is defined as the cross-power spectrum. If the functions ƒ and g have a zero mean (i.e., the spectral coefficients ƒ₀₀ and g₀₀ are zero), then S_ƒƒ(ℓ) and S_ƒg(ℓ) represent the contributions to the function's variance and covariance for degree ℓ, respectively. It is common that the (cross-)power spectrum is well approximated by a power law of the form

$S_{f\!f}(\ell) = C \, \ell^{\beta}.$

When β = 0, the spectrum is "white" as each degree possesses equal power. When β < 0, the spectrum is termed "red" as there is more power at the low degrees with long wavelengths than higher degrees. Finally, when β > 0, the spectrum is termed "blue". The condition on the order of growth of S_ƒƒ(ℓ) is related to the order of differentiability of ƒ in the next section.

Differentiability properties

One can also understand the differentiability properties of the original function ƒ in terms of the asymptotics of S_ƒƒ(ℓ). In particular, if S_ƒƒ(ℓ) decays faster than any rational function of ℓ as ℓ → ∞, then ƒ is infinitely differentiable. If, furthermore, S_ƒƒ(ℓ) decays exponentially, then ƒ is actually real analytic on the sphere.

The general technique is to use the theory of Sobolev spaces. Statements relating the growth of the S_ƒƒ(ℓ) to differentiability are then similar to analogous results on the growth of the coefficients of Fourier series. Specifically, if

$\sum_{\ell=0}^\infty (1+\ell^2)^s S_{ff}(\ell) < \infty,$

then ƒ is in the Sobolev space H^s(S²). In particular, the Sobolev embedding theorem implies that ƒ is infinitely differentiable provided that

$S_{ff}(\ell) = O(\ell^{-s})\quad\rm{as\ }\ell\to\infty$

for all s.

Algebraic properties

Addition theorem

A mathematical result of considerable interest and use is called the addition theorem for spherical harmonics. This is a generalization of the trigonometric identity

cos(θ' - θ) = cosθ'cosθ + sinθsinθ'

in which the role of the trigonometric functions appearing on the right-hand side is played by the spherical harmonics and that of the left-hand side is played by the Legendre polynomials.

Consider two unit vectors x and y, having spherical coordinates (θ,φ) and (θ′,φ′), respectively. The addition theorem states

$P_\ell( \mathbf{x}\cdot\mathbf{y} ) = \frac{4\pi}{2\ell+1}\sum_{m=-\ell}^\ell Y_{\ell m}^*(\theta',\varphi') \, Y_{\ell m}(\theta,\varphi).$

(1)

where P_ℓ is the Legendre polynomial of order ℓ. This expression is valid for both real and complex harmonics.^[4] The result can be proven analytically, using the properties of the Poisson kernel in the unit ball, or geometrically by applying a rotation to the vector y so that it points along the z-axis, and then directly calculating the right-hand side.^[5]

In particular, when x = y, this gives Unsöld's theorem^[6]

$\sum_{m=-\ell}^\ell Y_{\ell m}^*(\theta,\varphi) \, Y_{\ell m}(\theta,\varphi) = \frac{2\ell + 1}{4\pi}$

which generalizes the identity cos²θ + sin²θ = 1 to two dimensions.

In the expansion (1), the left-hand side P_ℓ(x·y) is a constant multiple of the degree ℓ zonal spherical harmonic. From this perspective, one has the following generalization to higher dimensions. Let Y_j be an arbitrary orthonormal basis of the space H_ℓ of degree ℓ spherical harmonics on the n-sphere. Then $Z^{(\ell)}_{\mathbf{x}}$ , the degree ℓ zonal harmonic corresponding to the unit vector x, decomposes as^[7]

$Z^{(\ell)}_{\mathbf{x}}({\mathbf{y}}) = \sum_{j=1}^{\dim(\mathbf{H}_\ell)}\overline{Y_j({\mathbf{x}})}\,Y_j({\mathbf{y}})$

(2)

Furthermore, the zonal harmonic $Z^{(\ell)}_{\mathbf{x}}({\mathbf{y}})$ is given as a constant multiple of the appropriate Gegenbauer polynomial:

$Z^{(\ell)}_{\mathbf{x}}({\mathbf{y}}) = C_\ell^{((n-1)/2)}({\mathbf{x}}\cdot {\mathbf{y}})$

(3)

Combining (2) and (3) gives (1) in dimension n = 2 when x and y are represented in spherical coordinates. Finally, evaluating at x = y gives the functional identity

$\frac{\dim \mathbf{H}_\ell}{\omega_{n-1}} = \sum_{j=1}^{\dim(\mathbf{H}_\ell)}|Y_j({\mathbf{x}})|^2$

where ω_n−1 is the volume of the (n−1)-sphere.

Clebsch-Gordan coefficients

Main article: Clebsch-Gordan coefficients

The Clebsch-Gordan coefficients are the coefficients appearing in the expansion of the product of two spherical harmonics in terms of spherical harmonics itself. A variety of techniques are available for doing essentially the same calculation, including the Wigner 3-jm symbol, the Racah coefficients, and the Slater integrals. Abstractly, the Clebsch-Gordan coefficients express the tensor product of two irreducible representations of the rotation group as a sum of irreducible representations: suitably normalized, the coefficients are then the multiplicities.

Parity

Main article: Parity (physics)

The spherical harmonics have well defined parity in the sense that they are either even or odd with respect to reflection about the origin. Reflection about the origin is represented by the operator $P\Psi(\vec{r}) = \Psi(-\vec{r})$ . For the spherical angles, ${θ,φ}$ this corresponds to the replacement ${π - θ,π + φ}$ . The associated Legendre functions gives (−1)^ℓ-m and from the exponential function we have (−1)^m, giving together for the spherical harmonics a parity of (-1)^ℓ:

$Y_\ell^m(\theta,\phi) \rightarrow Y_\ell^m(\pi-\theta,\pi+\phi) = (-1)^\ell Y_\ell^m(\theta,\phi)$

This remains true for spherical harmonics in higher dimensions: applying a point reflection to a spherical harmonic of degree ℓ changes the sign by a factor of (−1)^ℓ.

Visualization of the spherical harmonics

Schematic representation of $Y_{\ell m}$ on the unit sphere and its nodal lines. $Y_{\ell m}$ is equal to 0 along

m

great circles passing through the poles, and along $\ell-m$ circles of equal latitude. The function changes sign each time it crosses one of these lines.

3D color plot of the spherical harmonics of degree

n = 5

The Laplace spherical harmonics $Y_\ell^m$ can be visualized by considering their "nodal lines", that is, the set of points on the sphere where $Y_\ell^m$ vanishes. Nodal lines of $Y_\ell^m$ are comprised of circles: some are latitudes and others are longitudes. One can determine the number of nodal lines of each type by counting the number of zeros of $Y_\ell^m$ in the latitudinal and longitudinal directions independently. For the latitudinal direction, the associated Legendre functions possess ℓ−|m| zeros, whereas for the longitudinal direction, the trigonometric sin and cos functions possess 2|m| zeros.

When the spherical harmonic order m is zero (upper-left in the figure), the spherical harmonic functions do not depend upon longitude, and are referred to as zonal. Such spherical harmonics are a special case of zonal spherical functions. When ℓ = |m| (bottom-right in the figure), there are no zero crossings in latitude, and the functions are referred to as sectoral. For the other cases, the functions checker the sphere, and they are referred to as tesseral.

More general spherical harmonics of degree ℓ are not necessarily those of the Laplace basis $Y_\ell^m$ , and their nodal sets can be of a fairly general kind.^[8]

List of spherical harmonics

Main article: Table of spherical harmonics

Analytic expressions for the first few orthonormalized Laplace spherical harmonics that use the Condon-Shortley phase convention:

$Y_{0}^{0}(\theta,\varphi)={1\over 2}\sqrt{1\over \pi}$

$Y_{1}^{-1}(\theta,\varphi)={1\over 2}\sqrt{3\over 2\pi} \, \sin\theta \, e^{-i\varphi}$

$Y_{1}^{0}(\theta,\varphi)={1\over 2}\sqrt{3\over \pi}\, \cos\theta$

$Y_{1}^{1}(\theta,\varphi)={-1\over 2}\sqrt{3\over 2\pi}\, \sin\theta\, e^{i\varphi}$

$Y_{2}^{-2}(\theta,\varphi)={1\over 4}\sqrt{15\over 2\pi} \, \sin^{2}\theta \, e^{-2i\varphi}$

$Y_{2}^{-1}(\theta,\varphi)={1\over 2}\sqrt{15\over 2\pi}\, \sin\theta\, \cos\theta\, e^{-i\varphi}$

$Y_{2}^{0}(\theta,\varphi)={1\over 4}\sqrt{5\over \pi}\, (3\cos^{2}\theta-1)$

$Y_{2}^{1}(\theta,\varphi)={-1\over 2}\sqrt{15\over 2\pi}\, \sin\theta\,\cos\theta\, e^{i\varphi}$

$Y_{2}^{2}(\theta,\varphi)={1\over 4}\sqrt{15\over 2\pi}\, \sin^{2}\theta \, e^{2i\varphi}$

$Y_{3}^{0}(\theta,\varphi)={1\over 4}\sqrt{7\over \pi}\, (5\cos^{3}\theta-3\cos\theta)$

Higher dimensions

The classical spherical harmonics are defined as functions on the unit sphere S² inside three-dimensional Euclidean space. Spherical harmonics can be generalized to higher dimensional Euclidean space Rⁿ as follows.^[9] Let P_ℓ denote the space of homogeneous polynomials of degree ℓ in n variables. That is, a polynomial P is in P_ℓ provided that

$P(\lambda \mathbf{x}) = \lambda^\ell P(\mathbf{x}).$

Let A_ℓ denote the subspace of P_ℓ consisting of all harmonic polynomials; these are the solid spherical harmonics. Let H_ℓ denote the space of functions on the unit

$S^{n-1} = \{\mathbf{x}\in\mathbb{R}^n\,\mid\, |x|=1\}$

obtained by restriction from A_ℓ.

The following properties hold:

The spaces H_ℓ are dense in the set of continuous functions on Sⁿ⁻¹ with respect to the uniform topology, by the Stone-Weierstrass theorem. As a result, they are also dense in the space L²(Sⁿ⁻¹) of square-integrable functions on the sphere.

For all ƒ ∈ H_ℓ, one has

$\Delta_{S^{n-1}}f = -\ell(\ell+n-2)f.$

where Δ_Sⁿ⁻¹ is the Laplace-Beltrami operator on Sⁿ⁻¹. This operator is the analog of the angular part of the Laplacian in three dimensions; to wit, the Laplacian in n dimensions decomposes as

$\nabla^2 = r^{1-n}\frac{\partial}{\partial r}r^{n-1}\frac{\partial}{\partial r} + r^{-2}\Delta_{S^{n-1}}.$

It follows from the Stokes theorem and the preceding property that the spaces H_ℓ are orthogonal with respect to the inner product from L²(Sⁿ⁻¹). That is to say,

$\int_{S^{n-1}} f\bar{g}\,d\Omega = 0$

for ƒ ∈ H_ℓ and g ∈ H_k for k ≠ ℓ.

Conversely, the spaces H_ℓ are precisely the eigenspaces of Δ_Sⁿ⁻¹. In particular, an application of the spectral theorem to the Riesz potential $\Delta_{S^{n-1}}^{-1}$ gives another proof that the spaces H_ℓ are pairwise orthogonal and complete in L²(Sⁿ⁻¹).

Every homogeneous polynomial P ∈ P_ℓ can be uniquely written in the form

$P(x) = P_\ell(x) + |x|^2P_{\ell-2} + \cdots + \begin{cases} |x|^\ell P_0 & \ell \rm{\ even}\\ |x|^{\ell-1} P_1(x) & \ell\rm{\ odd} \end{cases}$

where P_j ∈ A_j. In particular,

$\dim \mathbf{H}_\ell = \binom{n+\ell-1}{n-1}-\binom{n+\ell-3}{n- 1}.$

An orthogonal basis of spherical harmonics in higher dimensions can be constructed inductively by the method of separation of variables, by solving the Sturm-Liouville problem for the spherical Laplacian

$\Delta_{S^{n-1}} = \sin^{2-n}\phi\frac{\partial}{\partial\phi}\sin^{n-2}\phi\frac{\partial}{\partial\phi} + \sin^{-2}\phi \Delta_{S^{n-2}}$

where φ is the axial coordinate in a spherical coordinate system on Sⁿ⁻¹.

Connection with representation theory

The space H_ℓ of spherical harmonics of degree ℓ is a representation of the symmetry group of rotations around a point (SO(3)) and its double-cover SU(2). Indeed, rotations act on the two-dimensional sphere, and thus also on H_ℓ by function composition

$\psi \mapsto \psi\circ\rho$

for ψ a spherical harmonic and ρ a rotation. The representation H_ℓ is an irreducible representation of SO(3).

The elements of H_ℓ arise as the restrictions to the sphere of elements of A_ℓ: harmonic polynomials homogeneous of degree ℓ on three-dimensional Euclidean space R³. By polarization of ψ ∈ A_ℓ, there are coefficients $\psi_{i_1\dots i_\ell}$ symmetric on the indices, uniquely determined by the requirement

$\psi(x_1,\dots,x_n) = \sum_{i_1\dots i_\ell}\psi_{i_1\dots i_\ell}x_{i_1}\cdots x_{i_\ell}.$

The condition that ψ be harmonic is equivalent to the assertion that the tensor $\psi_{i_1\dots i_\ell}$ must be trace free on every pair of indices. Thus as an irreducible representation of SO(3), H_ℓ is isomorphic to the space of traceless symmetric tensors of degree ℓ.

More generally, the analogous statements hold in higher dimensions: the space H_ℓ of spherical harmonics on the n-sphere is the irreducible representation of SO(n+1) corresponding to the traceless symmetric ℓ-tensors. However, whereas every irreducible tensor representation of SO(2) and SO(3) is of this kind, the special orthogonal groups in higher dimensions have additional irreducible representations that do not arise in this manner.

The special orthogonal groups have additional spin representations that are not tensor representations, and are typically not spherical harmonics. An exception are the spin representation of SO(3): strictly speaking these are representations of the double cover SU(2) of SO(3). In turn, SU(2) is identified with the group of unit quaternions, and so coincides with the 3-sphere. The spaces of spherical harmonics on the 3-sphere are certain spin representations of SO(3), with respect to the action by quaternionic multiplication.

Generalizations

The angle-preserving symmetries of the two-sphere are described by the group of Möbius transformations PSL(2,C). With respect to this group, the sphere is equivalent to the usual Riemann sphere. The group PSL(2,C) is isomorphic to the (proper) Lorentz group, and its action on the two-sphere agrees with the action of the Lorentz group on the celestial sphere in Minkowski space. The analog of the spherical harmonics for the Lorentz group is given by the hypergeometric series; furthermore, the spherical harmonics can be re-expressed in terms of the hypergeometric series, as SO(3) = PSU(2) is a subgroup of PSL(2,C).

More generally, hypergeometric series can be generalized to describe the symmetries of any symmetric space; in particular, hypergeometric series can be developed for any Lie group.^[10]^[11]^[12]^[13]

Notes

^ A historical account of various approaches to spherical harmonics in three-dimensions can be found in Chapter IV of MacRobert 1967. The term "Laplace spherical harmonics" is in common use; see Courant & Hilbert 1962 and Meijer & Bauer 2004.
^ More generally, one can consider solutions of the form R(r)Θ(θ)Φ(φ) and derive the form of R from that of Θ by solving a Sturm-Liouville equation with appropriate boundary conditions afterwards, but assuming that R has the required form in advance simplifies the exposition. For additional details, see Courant & Hilbert 1967
^ Watson & Whittaker 1927, p. 392.
^ This is valid for any orthonormal basis of spherical harmonics of degree ℓ. For unit power harmonics it is necessary to remove the factor of $4π$ .
^ Watson & Whittaker 1927, p. 395
^ Unsöld 1927
^ Stein & Weiss 1971, §IV.2
^ Eremenko, Jakobson & Nadirashvili 2007
^ Solomentsev 2001; Stein & Weiss 1971, §Iv.2
^ N. Vilenkin, Special Functions and the Theory of Group Representations, Am. Math. Soc. Transl., vol. 22, (1968).
^ J. D. Talman, Special Functions, A Group Theoretic Approach, (based on lectures by E.P. Wigner), W. A. Benjamin, New York (1968).
^ W. Miller, Symmetry and Separation of Variables, Addison-Wesley, Reading (1977).
^ A. Wawrzyńczyk, Group Representations and Special Functions, Polish Scientific Publishers. Warszawa (1984).

References

Cited references

Courant, Richard; Hilbert, David (1962), Methods of Mathematical Physics, Volume I, Wiley-Interscience .
Eremenko, Alexandre; Jakobson, Dmitry; Nadirashvili, Nikolai (2007), "On nodal sets and nodal domains on $S 2$ and $\mathbb{R}^2$ ", Université de Grenoble. Annales de l'Institut Fourier 57 (7): 2345–2360, MR 2394544, ISSN 0373-0956
MacRobert, T.M. (1967), Spherical harmonics: An elementary treatise on harmonic functions, with applications, Pergamon Press .
Meijer, Paul Herman Ernst; Bauer, Edmond (2004), Group theory: The application to quantum mechanics, Dover, ISBN 9780486437989 .
Solomentsev, E.D. (2001), "Spherical harmonics", in Hazewinkel, Michiel, Encyclopaedia of Mathematics, Kluwer Academic Publishers, ISBN 978-1556080104 .
Stein, Elias; Weiss, Guido (1971), Introduction to Fourier Analysis on Euclidean Spaces, Princeton, N.J.: Princeton University Press, ISBN 978-0-691-08078-9 .
Unsöld, Albrecht (1927), "Beiträge zur Quantenmechanik der Atome", Annalen der Physik 387 (3): 355–393, doi:10.1002/andp.19273870304 .
Watson, G. N.; Whittaker, E. T. (1927), A Course of Modern Analysis, Cambridge University Press .

General references

E.W. Hobson, The Theory of Spherical and Ellipsoidal Harmonics, (1955) Chelsea Pub. Co., ISBN 978-0828401043.
C. Müller, Spherical Harmonics, (1966) Springer, Lecture Notes in Mathematics, Vol. 17, ISBN 978-3-540-03600-5.
A.R. Edmonds, Angular Momentum in Quantum Mechanics, (1957) Princeton University Press, ISBN 0-691-07912-9.
E. U. Condon and G. H. Shortley, The Theory of Atomic Spectra, (1970) Cambridge at the University Press, ISBN 0-521-09209-4, See chapter 3.
J.D. Jackson, Classical Electrodynamics, ISBN 0-471-30932-X
Albert Messiah, Quantum Mechanics, volume II. (2000) Dover. ISBN 0-486-40924-4.
D. A. Varshalovich, A. N. Moskalev, V. K. Khersonskii Quantum Theory of Angular Momentum,(1988) World Scientific Publishing Co., Singapore, ISBN 9971-5-0107-4
Weisstein, Eric W., "Spherical harmonics" from MathWorld.

External links

This article's external links may not follow Wikipedia's content policies or guidelines. Please improve this article by removing excessive or inappropriate external links.

Interactive calculator of spherical harmonics on Tal Carmon's Research Homepage
Spherical harmonics applied to Acoustic Field analysis on Trinnov Audio's research page
Spherical Harmonics by Stephen Wolfram and Nodal Domains of Spherical Harmonics by Michael Trott, the Wolfram Demonstrations Project
An accessible introduction to spherical harmonics (by J. B. Calvert)
Citizendium:Spherical harmonics
OpenGL Spherical harmonics demo
Allen McNamara's spherical harmonics animations
Thorsten Becker's spherical harmonics animations

Software

This entry is from Wikipedia, the leading user-contributed encyclopedia. It may not have been reviewed by professional editors (see full disclaimer)

Donate to Wikimedia

Spherical harmonics

Contents