Spherical harmonics

From Wikipedia, the free encyclopedia

In mathematics, the spherical harmonics are the angular portion of an orthogonal set of solutions to Laplace's equation represented in a system of spherical coordinates. Spherical harmonics are important in many theoretical and practical applications, particularly in the computation of atomic electron configurations, the representation of the gravitational field, geoid, and magnetic field of planetary bodies, as well as characterization of the cosmic microwave background radiation.

[edit] Introduction

Laplace's equation in spherical coordinates is:

<math> \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</math>

(see also nabla in cylindrical and spherical coordinates). For <math>f(r,\theta,\phi)=R(r)\Theta(\theta)\Phi(\phi)</math>, the angular portion of Laplace's equation satisfies

<math>{\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} + l(l+1)\Theta(\theta)\Phi(\varphi) = 0</math>.

Using the technique of separation of variables, the angular solutions can be shown to be a products of trigonometric functions and associated Legendre functions:

<math> Y_\ell^m (\theta, \varphi ) = N \, e^{i m \varphi } \, P_\ell^m (\cos{\theta} )</math>,

where <math>Y_\ell^m</math> is a called a spherical harmonic function of degree <math>\ell</math> and order <math>m</math>, <math>P_\ell^m</math> is an associated Legendre function, <math>N</math> is a normalization constant, and <math>\theta</math> and <math>\varphi</math> represent colatitude and longitude, respectively. The spherical coordinates used in this article are consistent those used by physicists, but differ from those employed by mathematicians (see spherical coordinates). In particular, the colatitude <math>\theta</math>, or polar angle, ranges from <math>0\leq\theta\leq\pi</math> and the longitude <math>\varphi</math>, or azimuth, ranges from <math>0\leq\varphi<2\pi</math>.

When Laplace's equation is solved on the surface of the sphere, the periodic boundary conditions in <math>\varphi</math>, as well as regularity conditions at both the north and south poles, ensure that the degree <math>\ell</math> and order <math>m</math> are integers that satisfy <math>\ell \ge 0</math> and <math>|m| \le \ell</math>. In contrast, if the function <math>f</math> were only to have been defined for <math>\theta \le \theta_0</math>, then the resulting spherical cap harmonics would have been defined for integer order, but non-integer degree. The general solution to Laplace's equation is a linear combination of the spherical harmonic functions multplied by the solutions of <math>R(r)</math>:

<math> f(r, \theta, \varphi) = \sum_{\ell=0}^\infty \sum_{m=-\ell}^\ell r^{-1-\ell} \, f_\ell^m \, Y_\ell^m (\theta, \varphi ) +

\sum_{\ell=0}^\infty \sum_{m=-\ell}^\ell  r^\ell  \, f_\ell^{m'} \, Y_\ell^m (\theta, \varphi ) </math>,

where <math>f_\ell^m</math> and <math>f_\ell^{m'}</math> are constants. The terms in the first summation approach zero as <math>r</math> goes to infinity, whereas the terms in the second summation approach zero at the origin.

[edit] Normalizations

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

<math> 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 } </math>

which are orthonormalized

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

where δ_aa = 1, δ_ab = 0 if a ≠ b, and <math>d\Omega=\sin\theta\,d\varphi\,d\theta</math>. The disciplines of geodesy and spectral analysis use

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

which possess unit power

<math>{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'}</math>.

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

<math> Y_\ell^m( \theta , \varphi ) = \sqrtTemplate:(\ell-m)!\over (\ell+m)! \, P_\ell^m ( \cos{\theta} ) \, e^{i m \varphi }</math>

which have the normalization

<math> \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'}</math>.

Using the identity (see associated Legendre functions)

<math>P_\ell ^{-m} = (-1)^m \frac{(\ell-m)!}{(\ell+m)!} P_\ell ^{m}</math>

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

<math>Y_\ell^{m*} (\theta, \varphi) = (-1)^m Y_\ell^{-m} (\theta, \varphi)</math>,

where the superscript * denotes complex conjugation.

[edit] Condon-Shortley Phase

One source of confusion with the definition of the spherical harmonic functions concerns a phase factor of <math>(-1)^m\,</math>, 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. It should be noted that 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 mechnical operations. It should be noted that the geodesy and magnetics communities never include the Condon-Shortley phase factor in their definitions of the spherical harmonic functions.

[edit] Spherical harmonics expansion

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

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

This expansion is exact as long as <math>\ell</math> goes to infinity. Truncation errors will arise when limiting the sum over <math>\ell</math> to a finite bandwidth <math>L</math>. The expansion coefficients can be obtained by multiplying the above equation by the complex conjugate of a spherical harmonic, integrating over the solid angle <math>\Omega\!\,</math>, and utilizing the above orthogonality relationships. For the case of orthonormalized harmonics, this gives

<math>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)</math>.

An alternative set of spherical harmonics for real functions may be obtained by taking the set

<math>Y_{\ell m} = \begin{cases}

Y_\ell^0 \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\quad\quad\mbox{ if } m=0\\ {1\over\sqrt2}\left(Y_\ell^m+(-1)^m \, Y_\ell^{-m}\right) = \sqrt{2} N P_\ell^m(\theta) \cos m\varphi \qquad\quad\quad \mbox{if } m>0 \\ {1\over i\sqrt2}\left(Y_\ell^{|m|}-(-1)^{|m|}\, Y_\ell^{-|m|}\right) = \sqrt{2} N P_\ell^{|m|}(\theta) \sin |m|\varphi \quad\mbox{ if } m<0. \end{cases} </math>

These functions have the same normalization properties as the complex ones above. In this notation, a real square-integrable function can be expressed as an infinite sum of real spherical harmonics as

<math> f(\theta, \varphi) = \sum_{\ell=0}^\infty \sum_{m=-\ell}^\ell f_{lm} \, Y_{lm}(\theta, \varphi) </math>.

[edit] Spectrum Analysis

The total power of a function <math>f</math> is defined in the signal processing literature as the integral of the function squared, divided by the area it spans. 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:

<math>\frac{1}{4 \, \pi} \int_\Omega f(\Omega)^2\, d\Omega = \sum_{l=0}^\infty S_{f\!f}(l)</math>,

where

<math>S_{f\!f}(l) = \sum_{m=-l}^l f_{lm}^2 </math>

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

<math>\frac{1}{4 \, \pi} \int_\Omega f(\Omega) \, g(\Omega) \, d\Omega = \sum_{l=0}^\infty S_{fg}(l)</math>,

where

<math>S_{fg}(l) = \sum_{m=-l}^l f_{lm} g_{lm} </math>

is defined as the cross-power spectrum. If the functions <math>f</math> and <math>g</math> have a zero mean (i.e., the spectral coefficients <math>f_{00}</math> and <math>g_{00}</math> are zero), then <math>S_{f\!f}(l)</math> and <math>S_{fg}(l)</math> represent the contributions to the function's variance and covariance for degree <math>\ell</math>, respectively. It is common that the (cross-)power spectrum is well approximated by a power law of the form

<math>S_{f\!f}(l) = C \, \ell^{\beta}</math>.

When <math>\beta=0</math>, the spectrum is "white" as each degree possesses equal power. When <math>\beta<0</math>, the spectrum is termed "red" as there is more power at the low degrees with long wavelengths than higher degrees. Finally, when <math>\beta>0</math>, the spectrum is termed "blue".

[edit] Addition theorem

A mathematical result of considerable interest and use is called the addition theorem for spherical harmonics. Two vectors r and r', with spherical coordinates <math>(r,\theta,\varphi)</math> and <math>(r ',\theta ',\varphi ')</math>,respectively, have an angle <math>\gamma</math> between them given by

<math>\cos\gamma=\cos\theta\cos\theta'+\sin\theta\sin\theta'\cos(\varphi-\varphi')</math>.

The addition theorem expresses a Legendre polynomial of order <math>l</math> in the angle <math>\gamma</math> in terms of products of two spherical harmonics with angular coordinates <math>(\theta,\varphi)</math> and <math>(\theta',\varphi')</math>:

<math>P_l( \cos \gamma ) = \frac{4\pi}{2l+1}\sum_{m=-l}^l Y_{lm}^*(\theta',\phi') \, Y_{lm}(\theta,\varphi) </math>.

This expression is valid for both real and complex harmonics. However, it should be emphasized that the quoted form above is valid only for the orthonormalized spherical harmonics. For unit power harmonics it is only necessary to remove the factor of <math>4 \pi</math>.

[edit] Visualization of the spherical harmonics

The spherical harmonics are easily visualized by counting the number of zero crossings they possess in both the latitudinal and longitudinal directions. For the latitudinal direction, the associated Legendre functions possess <math>l-|m|</math> zeros, whereas for the longitudinal direction, the trigonomentric <math>\sin</math> and <math>\cos</math> functions possess <math>2|m|</math> zeros.

When the spherical harmonic order <math>m</math> is zero, the spherical harmonic functions do not depend upon longitude, and are referred to as zonal. When <math>l=|m|</math>, 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.

[edit] First few spherical harmonics

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

<math>Y_{0}^{0}(\theta,\varphi)={1\over 2}\sqrt{1\over \pi}</math>

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

<math>Y_{1}^{0}(\theta,\varphi)={1\over 2}\sqrt{3\over \pi}\, \cos\theta</math>

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

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

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

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

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

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

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

More spherical harmonics up to Y₁₀

[edit] Generalizations

The spherical harmonics map can be seen as representations of the symmetry group of rotations around a point (SO(3)) and its double-cover SU(2). As such they capture the symmetry of the two-dimensional sphere (or two-sphere). Each set of spherical harmonics with a given value for the l-parameter map onto a different irreducible representation of SO(3).

In addition, the two-sphere is equivalent to the Riemann sphere. The complete set of symmetries of the Riemann sphere are described by the Mobius transformation group SL(2,C), of which the Lorentz group is a representation. The analog of the spherical harmonics for the Lorentz group are given by the hypergeometric series; indeed, the spherical harmonics can be re-expressed in terms of the hypergeometric series, as SO(3) is a subgroup of SL(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<ref> N. Vilenkin, Special Functions and the Theory of Group Representations, Am. Math. Soc. Transl., vol. 22, (1968). </ref><ref> J. D. Talman, Special Functions, A Group Theoretic Approach, (based on lectures by E.P. Wigner), W. A. Benjamin, New York (1968).</ref><ref> W. Miller, Symmetry and Separation of Variables, Addison-Wesley, Reading (1977). </ref><ref> A. Wawrzyńczyk, Group Representations and Special Functions, Polish Scientific Publishers. Warszawa (1984). </ref>

[edit] See also

• Clebsch-Gordan coefficients
• Harmonic function
• Rotation group
• Sturm-Liouville theory

[edit] References

[edit] Software

• SHTOOLS: Fortran 95 software archive
• HEALPIX: Fortran 90 and C++ software archive
• SpherePack: Fortran 77 software archivecs:Sférické harmonické funkce

de:Kugelflächenfunktionen fr:Harmonique sphérique it:Armoniche sferiche ja:球面調和関数