Spin-weighted spherical harmonics

In special functions, a topic in mathematics, spin-weighted spherical harmonics are generalizations of the standard spherical harmonics and—like the usual spherical harmonics—are functions on the sphere. Unlike ordinary spherical harmonics, the spin-weighted harmonics are U(1) gauge fields rather than scalar fields: mathematically, they take values in a complex line bundle. The spin-weighted harmonics are organized by degree Template:Mvar, just like ordinary spherical harmonics, but have an additional spin weight Template:Mvar that reflects the additional U(1) symmetry. A special basis of harmonics can be derived from the Laplace spherical harmonics Y_lm, and are typically denoted by _sY_lm, where Template:Mvar and Template:Mvar are the usual parameters familiar from the standard Laplace spherical harmonics. In this special basis, the spin-weighted spherical harmonics appear as actual functions, because the choice of a polar axis fixes the U(1) gauge ambiguity. The spin-weighted spherical harmonics can be obtained from the standard spherical harmonics by application of spin raising and lowering operators. In particular, the spin-weighted spherical harmonics of spin weight s = 0 are simply the standard spherical harmonics:

${}_0{Y}_{lm}=Y_{lm}.$

Spaces of spin-weighted spherical harmonics were first identified in connection with the representation theory of the Lorentz group . They were subsequently and independently rediscovered by and applied to describe gravitational radiation, and again by as so-called "monopole harmonics" in the study of Dirac monopoles.

Spin-weighted functions

Regard the sphere S² as embedded into the three-dimensional Euclidean space R³. At a point x on the sphere, a positively oriented orthonormal basis of tangent vectors at x is a pair a, b of vectors such that

$\begin{array}{rl}\mathbf{𝐱}\cdot \mathbf{𝐚}=\mathbf{𝐱}\cdot \mathbf{𝐛}& =0\\ \mathbf{𝐚}\cdot \mathbf{𝐚}=\mathbf{𝐛}\cdot \mathbf{𝐛}& =1\\ \mathbf{𝐚}\cdot \mathbf{𝐛}& =0\\ \mathbf{𝐱}\cdot (\mathbf{𝐚}\times \mathbf{𝐛})& >0,\end{array}$

where the first pair of equations states that a and b are tangent at x, the second pair states that a and b are unit vectors, the penultimate equation that a and b are orthogonal, and the final equation that (x, a, b) is a right-handed basis of R³.

A spin-weight Template:Mvar function Template:Mvar is a function accepting as input a point x of S² and a positively oriented orthonormal basis of tangent vectors at x, such that

$f(\mathbf{𝐱},(\cos \theta )\mathbf{𝐚}-(\sin \theta )\mathbf{𝐛},(\sin \theta )\mathbf{𝐚}+(\cos \theta )\mathbf{𝐛})=e^{is\theta }f(\mathbf{𝐱},\mathbf{𝐚},\mathbf{𝐛})$

for every rotation angle Template:Mvar.

Following , denote the collection of all spin-weight Template:Mvar functions by B(s). Concretely, these are understood as functions Template:Mvar on C²\{0} satisfying the following homogeneity law under complex scaling

$f(\lambda z,\bar{\lambda }\overline{z})={\left(\frac{\bar{\lambda }}{\lambda }\right)}^{s}f(z,\overline{z}).$

This makes sense provided Template:Mvar is a half-integer.

Abstractly, B(s) is isomorphic to the smooth vector bundle underlying the antiholomorphic vector bundle O(2s) of the Serre twist on the complex projective line CP¹. A section of the latter bundle is a function Template:Mvar on C²\{0} satisfying

$g(\lambda z,\bar{\lambda }\overline{z})=\stackrel{2s}{\bar{\lambda }}g(z,\overline{z}).$

Given such a Template:Mvar, we may produce a spin-weight Template:Mvar function by multiplying by a suitable power of the hermitian form

$P(z,\overline{z})=z\cdot \overline{z}.$

Specifically, f = P^−sg is a spin-weight Template:Mvar function. The association of a spin-weighted function to an ordinary homogeneous function is an isomorphism.

The operator Template:Mvar

The spin weight bundles B(s) are equipped with a differential operator Template:Mvar (eth). This operator is essentially the Dolbeault operator, after suitable identifications have been made,

$\partial :\overline{\mathbf{𝐎}(2s)}\to {{ℰ}}^{1,0}\otimes \overline{\mathbf{𝐎}(2s)}\cong \overline{\mathbf{𝐎}(2s)}\otimes \mathbf{𝐎}(-2).$

Thus for f ∈ B(s),

$\mathrm{ð}f\stackrel{\text{def}}{=}{P}^{-s+1}\partial \left(P^{s}f\right)$

defines a function of spin-weight s + 1.

Spin-weighted harmonics

Just as conventional spherical harmonics are the eigenfunctions of the Laplace-Beltrami operator on the sphere, the spin-weight Template:Mvar harmonics are the eigensections for the Laplace-Beltrami operator acting on the bundles Template:Mathcal(s) of spin-weight Template:Mvar functions.

Representation as functions

The spin-weighted harmonics can be represented as functions on a sphere once a point on the sphere has been selected to serve as the North pole. By definition, a function Template:Mvar with spin weight Template:Mvar transforms under rotation about the pole via

$\eta \to e^{is\psi }\eta .$

Working in standard spherical coordinates, we can define a particular operator Template:Mvar acting on a function Template:Mvar as:

$\mathrm{ð}\eta =-{(\sin \theta )}^s\{\frac{\partial }{\partial \theta }+\frac{i}{\sin \theta }\frac{\partial }{\partial \varphi }\}\left[{(\sin \theta )}^{-s}\eta \right].$

This gives us another function of Template:Mvar and Template:Mvar. (The operator Template:Mvar is effectively a covariant derivative operator in the sphere.)

An important property of the new function Template:Mvar is that if Template:Mvar had spin weight Template:Mvar, Template:Mvar has spin weight s + 1. Thus, the operator raises the spin weight of a function by 1. Similarly, we can define an operator Template:Mvar which will lower the spin weight of a function by 1:

$\overline{\mathrm{ð}}\eta =-{(\sin \theta )}^{-s}\{\frac{\partial }{\partial \theta }-\frac{i}{\sin \theta }\frac{\partial }{\partial \varphi }\}\left[{(\sin \theta )}^s\eta \right].$

The spin-weighted spherical harmonics are then defined in terms of the usual spherical harmonics as:

${}_s{Y}_{lm}=\{\begin{array}{lll}\sqrt{\frac{(l-s)!}{(l+s)!}}{\mathrm{ð}}^s{Y}_{lm},& & 0\le s\le l;\\ \sqrt{\frac{(l+s)!}{(l-s)!}}{(-1)}^s{\overline{\mathrm{ð}}}^{-s}{Y}_{lm},& & -l\le s\le 0;\\ 0,& & l<|s|.\end{array}$

The functions _sY_lm then have the property of transforming with spin weight Template:Mvar.

Other important properties include the following:

$\begin{array}{rl}\mathrm{ð}\left({}_s{Y}_{lm}\right)& =+\sqrt{(l-s)(l+s+1)}{}_{s+1}{Y}_{lm};\\ \overline{\mathrm{ð}}\left({}_s{Y}_{lm}\right)& =-\sqrt{(l+s)(l-s+1)}{}_{s-1}{Y}_{lm};\end{array}$

Orthogonality and completeness

The harmonics are orthogonal over the entire sphere:

$\underset{S^2}{\int }{}_s{Y}_{lm}{}_s{\overline{Y}}_{l^{\prime }{m}^{\prime }}dS={\delta }_{l{l}^{\prime }}{\delta }_{m{m}^{\prime }},$

and satisfy the completeness relation

$\sum _{lm}{}_s{\overline{Y}}_{lm}({\theta }^{\prime },{\varphi }^{\prime }){}_s{Y}_{lm}(\theta ,\varphi )=\delta ({\varphi }^{\prime }-\varphi )\delta (\cos {\theta }^{\prime }-\cos \theta )$

Calculating

These harmonics can be explicitly calculated by several methods. The obvious recursion relation results from repeatedly applying the raising or lowering operators. Formulae for direct calculation were derived by . Note that their formulae use an old choice for the Condon–Shortley phase. The convention chosen below is in agreement with Mathematica, for instance.

The more useful of the Goldberg, et al., formulae is the following:

${}_s{Y}_{lm}(\theta ,\varphi )={(-1)}^{l+m-s}\sqrt{\frac{(l+m)!(l-m)!(2l+1)}{4\pi (l+s)!(l-s)!}}{\sin }^{2l}\left(\frac{\theta }{2}\right)e^{im\varphi }\times \sum _{r=0}^{l-s}{(-1)}^r(\genfrac{}{}{0ex}{}{l-s}{r})(\genfrac{}{}{0ex}{}{l+s}{r+s-m}){\cot }^{2r+s-m}\left(\frac{\theta }{2}\right).$

A Mathematica notebook using this formula to calculate arbitrary spin-weighted spherical harmonics can be found here.

With the phase convention here:

$\begin{array}{rl}{}_s{\overline{Y}}_{lm}& ={(-1)}^{s+m}{}_{-s}{Y}_{l(-m)}\\ {}_s{Y}_{lm}(\pi -\theta ,\varphi +\pi )& ={(-1)}^l{}_{-s}{Y}_{lm}(\theta ,\varphi ).\end{array}$

First few spin-weighted spherical harmonics

Analytic expressions for the first few orthonormalized spin-weighted spherical harmonics:

Spin-weight s = 1, degree l = 1

$\begin{array}{rl}{}_1{Y}_{10}(\theta ,\varphi )& =\sqrt{\frac{3}{8\pi }}\sin \theta \\ {}_1{Y}_{1\pm 1}(\theta ,\varphi )& =-\sqrt{\frac{3}{16\pi }}(1\mp \cos \theta )e^{\pm i\varphi }\end{array}$

Relation to Wigner rotation matrices

$D_{-ms}^l(\varphi ,\theta ,-\psi )={(-1)}^m\sqrt{\frac{4\pi }{2l+1}}{}_s{Y}_{lm}(\theta ,\varphi )e^{is\psi }$

This relation allows the spin harmonics to be calculated using recursion relations for the [[Wigner D-matrix|Template:Mvar-matrices]].

Triple integral

The triple integral in the case that s₁ + s₂ + s₃ = 0 is given in terms of the [[3-j symbol|3-Template:Mvar symbol]]:

$\underset{S^2}{\int }{}_{s_1}{Y}_{j_1{m}_1}{}_{s_2}{Y}_{j_2{m}_2}{}_{s_3}{Y}_{j_3{m}_3}=\sqrt{\frac{(2j_1+1)(2j_2+1)(2j_3+1)}{4\pi }}\left(\begin{array}{ccc}{j}_1& j_2& j_3\\ m_1& m_2& m_3\end{array}\right)\left(\begin{array}{ccc}{j}_1& j_2& j_3\\ -s_1& -s_2& -s_3\end{array}\right)$

See also

• Spherical basis

References

• .
• .
• ; (1963) Representations of the rotation and Lorentz groups and their applications (translation). Macmillan Publishers.
• (Note: As mentioned above, this paper uses a choice for the Condon-Shortley phase that is no longer standard.)
• .
• .