Jump to content

Vector spherical harmonics

From Wikipedia, the free encyclopedia

In mathematics, vector spherical harmonics (VSH) are an extension of the scalar spherical harmonics for use with vector fields. The components of the VSH are complex-valued functions expressed in the spherical coordinate basis vectors.

Definition

[edit]

Several conventions have been used to define the VSH.[1][2][3][4][5] We follow that of Barrera et al.. Given a scalar spherical harmonic Yℓm(θ, φ), we define three VSH:

  • 𝐘ℓm=Yℓm𝐫̂,
  • 𝜳ℓm=r∇Yℓm,
  • 𝜱ℓm=𝐫×∇Yℓm,

with 𝐫̂ being the unit vector along the radial direction in spherical coordinates and 𝐫 the vector along the radial direction with the same norm as the radius, i.e., 𝐫=r𝐫̂. The radial factors are included to guarantee that the dimensions of the VSH are the same as those of the ordinary spherical harmonics and that the VSH do not depend on the radial spherical coordinate.

The interest of these new vector fields is to separate the radial dependence from the angular one when using spherical coordinates, so that a vector field admits a multipole expansion

𝐄=∑ℓ=0∞∑m=−ℓℓ(Eℓmr(r)𝐘ℓm+Eℓm(1)(r)𝜳ℓm+Eℓm(2)(r)𝜱ℓm).

The labels on the components reflect that Eℓmr is the radial component of the vector field, while Eℓm(1) and Eℓm(2) are transverse components (with respect to the radius vector 𝐫).

In physics

[edit]

In physics, the vector spherical harmonics 𝐘j,ℓ,smjare defined as spin s=1 eigenfunctions of the angular momentum operators J2,Jz,L2, and S2, where 𝐉=𝐋+𝐒 is the total angular momentum.[6] They are written as𝐘j,ℓ,1mj(𝐤)=∑mℓ=−ℓ+ℓ ∑ms=−1+1⟨j mj|ℓ 1 mℓ ms⟩Yℓmℓ(𝐤)𝐞̂ms,which are linear combinations of the scalar spherical harmonics Yℓmℓ with the vector angular momentum basis𝐞̂±1=∓𝐱̂±i𝐲̂2,𝐞̂0=𝐳̂.using the Clebsch-Gordan coefficients ⟨j mj|ℓ 1 mℓ ms⟩.

Because vector bosons such as the photon are spin-one, the vector spherical harmonics are commonly used in physics to describe vector and pseudovector interactions, such as electromagnetic transitions, in atomic and nuclear systems. They are a special (s=1) case of the spin spherical harmonics.

To derive these relations, one begins with the plane-wave expansion for plane waves with vector polarization.

Main properties

[edit]

Symmetry

[edit]

Like the scalar spherical harmonics, the VSH satisfy

𝐘ℓ,−m=(−1)m𝐘ℓm∗,𝜳ℓ,−m=(−1)m𝜳ℓm∗,𝜱ℓ,−m=(−1)m𝜱ℓm∗,

which cuts the number of independent functions roughly in half. The star indicates complex conjugation.

Orthogonality

[edit]

The VSH are orthogonal in the usual three-dimensional way at each point 𝐫:

𝐘ℓm(𝐫)⋅𝜳ℓm(𝐫)=0,𝐘ℓm(𝐫)⋅𝜱ℓm(𝐫)=0,𝜳ℓm(𝐫)⋅𝜱ℓm(𝐫)=0.

They are also orthogonal in Hilbert space:

∫𝐘ℓm⋅𝐘ℓ′m′∗dΩ=δℓℓ′δmm′,∫𝜳ℓm⋅𝜳ℓ′m′∗dΩ=ℓ(ℓ+1)δℓℓ′δmm′,∫𝜱ℓm⋅𝜱ℓ′m′∗dΩ=ℓ(ℓ+1)δℓℓ′δmm′,∫𝐘ℓm⋅𝜳ℓ′m′∗dΩ=0,∫𝐘ℓm⋅𝜱ℓ′m′∗dΩ=0,∫𝜳ℓm⋅𝜱ℓ′m′∗dΩ=0.

An additional result at a single point 𝐫 (not reported in Barrera et al., 1985) is, for all ℓ,m,ℓ′,m′,

𝐘ℓm(𝐫)⋅𝜳ℓ′m′(𝐫)=0,𝐘ℓm(𝐫)⋅𝜱ℓ′m′(𝐫)=0.

Vector multipole moments

[edit]

The orthogonality relations allow one to compute the spherical multipole moments of a vector field as

Eℓmr=∫𝐄⋅𝐘ℓm∗dΩ,Eℓm(1)=1ℓ(ℓ+1)∫𝐄⋅𝜳ℓm∗dΩ,Eℓm(2)=1ℓ(ℓ+1)∫𝐄⋅𝜱ℓm∗dΩ.

The gradient of a scalar field

[edit]

Given the multipole expansion of a scalar field

ϕ=∑ℓ=0∞∑m=−ℓℓϕℓm(r)Yℓm(θ,ϕ),

we can express its gradient in terms of the VSH as

∇ϕ=∑ℓ=0∞∑m=−ℓℓ(dϕℓmdr𝐘ℓm+ϕℓmr𝜳ℓm).

Divergence

[edit]

For any multipole field we have

∇⋅(f(r)𝐘ℓm)=(dfdr+2rf)Yℓm,∇⋅(f(r)𝜳ℓm)=−ℓ(ℓ+1)rfYℓm,∇⋅(f(r)𝜱ℓm)=0.

By superposition we obtain the divergence of any vector field:

∇⋅𝐄=∑ℓ=0∞∑m=−ℓℓ(dEℓmrdr+2rEℓmr−ℓ(ℓ+1)rEℓm(1))Yℓm.

We see that the component on Φℓm is always solenoidal.

Curl

[edit]

For any multipole field we have

∇×(f(r)𝐘ℓm)=−1rf𝜱ℓm,∇×(f(r)𝜳ℓm)=(dfdr+1rf)𝜱ℓm,∇×(f(r)𝜱ℓm)=−ℓ(ℓ+1)rf𝐘ℓm−(dfdr+1rf)𝜳ℓm.

By superposition we obtain the curl of any vector field:

∇×𝐄=∑ℓ=0∞∑m=−ℓℓ(−ℓ(ℓ+1)rEℓm(2)𝐘ℓm−(dEℓm(2)dr+1rEℓm(2))𝜳ℓm+(−1rEℓmr+dEℓm(1)dr+1rEℓm(1))𝜱ℓm).

Laplacian

[edit]

The action of the Laplace operator Δ=∇⋅∇ separates as follows:

Δ(f(r)𝐙ℓm)=(1r2∂∂rr2∂f∂r)𝐙ℓm+f(r)Δ𝐙ℓm, where 𝐙ℓm=𝐘ℓm,𝜳ℓm,𝜱ℓm and

Δ𝐘ℓm=−1r2(2+ℓ(ℓ+1))𝐘ℓm+2r2𝜳ℓm,Δ𝜳ℓm=2ℓ(ℓ+1)r2𝐘ℓm−1r2ℓ(ℓ+1)𝜳ℓm,Δ𝜱ℓm=−1r2ℓ(ℓ+1)𝜱ℓm.

Also note that this action becomes symmetric, i.e. the off-diagonal coefficients are equal to 2r2ℓ(ℓ+1), for properly normalized VSH.

Examples

[edit]
𝜳1m
𝜳2m
𝜳3m
𝜱1m
𝜱2m
𝜱3m
Visualizations of the real parts of ℓ=1,2,3 VSHs. Click to expand.

First vector spherical harmonics

[edit]
  • ℓ=0. 𝐘00=14π𝐫̂,𝜳00=𝟎,𝜱00=𝟎.
  • ℓ=1. 𝐘10=34πcos⁡θ𝐫̂,𝐘11=−38πeiφsin⁡θ𝐫̂, 𝜳10=−34πsin⁡θ𝜽̂,𝜳11=−38πeiφ(cos⁡θ𝜽̂+i𝝋̂), 𝜱10=−34πsin⁡θ𝝋̂,𝜱11=38πeiφ(i𝜽̂−cos⁡θ𝝋̂).
  • ℓ=2. 𝐘20=145π(3cos2θ−1)𝐫̂,𝐘21=−158πsin⁡θcos⁡θeiφ𝐫̂,𝐘22=14152πsin2θe2iφ𝐫̂. 𝜳20=−325πsin⁡θcos⁡θ𝜽̂,𝜳21=−158πeiφ(cos⁡2θ𝜽̂+icos⁡θ𝝋̂),𝜳22=158πsin⁡θe2iφ(cos⁡θ𝜽̂+i𝝋̂). 𝜱20=−325πsin⁡θcos⁡θ𝝋̂,𝜱21=158πeiφ(icos⁡θ𝜽̂−cos⁡2θ𝝋̂),𝜱22=158πsin⁡θe2iφ(−i𝜽̂+cos⁡θ𝝋̂).

Expressions for negative values of m are obtained by applying the symmetry relations.

Applications

[edit]

Electrodynamics

[edit]

The VSH are especially useful in the study of multipole radiation fields. For instance, a magnetic multipole is due to an oscillating current with angular frequency ω and complex amplitude

𝐉̂=J(r)𝜱ℓm,

and the corresponding electric and magnetic fields, can be written as

𝐄̂=E(r)𝜱ℓm,𝐁̂=Br(r)𝐘ℓm+B(1)(r)𝜳ℓm.

Substituting into Maxwell equations, Gauss's law is automatically satisfied

∇⋅𝐄̂=0,

while Faraday's law decouples as

∇×𝐄̂=−iω𝐁̂⇒{ℓ(ℓ+1)rE=iωBr,dEdr+Er=iωB(1).

Gauss' law for the magnetic field implies

∇⋅𝐁̂=0⇒dBrdr+2rBr−ℓ(ℓ+1)rB(1)=0,

and Ampère–Maxwell's equation gives

∇×𝐁̂=μ0𝐉̂+iμ0ε0ω𝐄̂⇒−Brr+dB(1)dr+B(1)r=μ0J+iωμ0ε0E.

In this way, the partial differential equations have been transformed into a set of ordinary differential equations.

Alternative definition

[edit]
Angular part of magnetic and electric vector spherical harmonics. Red and green arrows show the direction of the field. Generating scalar functions are also presented, only the first three orders are shown (dipoles, quadrupoles, octupoles).

In many applications, vector spherical harmonics are defined as fundamental set of the solutions of vector Helmholtz equation in spherical coordinates.[7][8]

In this case, vector spherical harmonics are generated by scalar functions, which are solutions of scalar Helmholtz equation with the wavevector 𝐤. ψemn=cos⁡mφPnm(cos⁡ϑ)zn(kr)ψomn=sin⁡mφPnm(cos⁡ϑ)zn(kr) here Pnm(cos⁡θ) are the associated Legendre polynomials, and zn(kr) are any of the spherical Bessel functions.

Vector spherical harmonics are defined as:

longitudinal harmonics
𝐋eomn=𝜵ψeomn
magnetic harmonics
𝐌oemn=∇×(𝐫ψoemn)
electric harmonics
𝐍oemn=∇×𝐌oemnk

Here we use harmonics real-valued angular part, where m≥0, but complex functions can be introduced in the same way.

Let us introduce the notation ρ=kr. In the component form vector spherical harmonics are written as: 𝐌emn(k,𝐫)=−msin⁡(θ)sin⁡(mφ)Pnm(cos⁡(θ))zn(ρ)𝐞θ−cos⁡(mφ)dPnm(cos⁡(θ))dθzn(ρ)𝐞φ 𝐌omn(k,𝐫)=msin⁡(θ)cos⁡(mφ)Pnm(cos⁡(θ))zn(ρ)𝐞θ−sin⁡(mφ)dPnm(cos⁡(θ))dθzn(ρ)𝐞φ

𝐍emn(k,𝐫)=zn(ρ)ρcos⁡(mφ)n(n+1)Pnm(cos⁡(θ))𝐞𝐫+cos⁡(mφ)dPnm(cos⁡(θ))dθ1ρddρ[ρzn(ρ)]𝐞θ−msin⁡(mφ)Pnm(cos⁡(θ))sin⁡(θ)1ρddρ[ρzn(ρ)]𝐞φ

𝐍omn(k,𝐫)=zn(ρ)ρsin⁡(mφ)n(n+1)Pnm(cos⁡(θ))𝐞𝐫+sin⁡(mφ)dPnm(cos⁡(θ))dθ1ρddρ[ρzn(ρ)]𝐞θ+mcos⁡(mφ)Pnm(cos⁡(θ))sin⁡(θ)1ρddρ[ρzn(ρ)]𝐞φ There is no radial part for magnetic harmonics. For electric harmonics, the radial part decreases faster than angular, and for big ρ can be neglected. We can also see that for electric and magnetic harmonics angular parts are the same up to permutation of the polar and azimuthal unit vectors, so for big ρ electric and magnetic harmonics vectors are equal in value and perpendicular to each other.

Longitudinal harmonics: 𝐋eomn(k,𝐫)=∂∂rzn(kr)Pnm(cos⁡θ)cossinmφ𝐞r+1rzn(kr)∂∂θPnm(cos⁡θ)cossinmφ𝐞θ∓mrsin⁡θzn(kr)Pnm(cos⁡θ)sincosmφ𝐞φ

Orthogonality

[edit]

The solutions of the Helmholtz vector equation obey the following orthogonality relations:[8] ∫02π∫0π𝐋oemn⋅𝐋oemnsin⁡ϑdϑdφ=(1+δm,0)2π(2n+1)2(n+m)!(n−m)!k2{n[zn−1(kr)]2+(n+1)[zn+1(kr)]2}∫02π∫0π𝐌oemn⋅𝐌oemnsin⁡ϑdϑdφ=(1+δm,0)2π2n+1(n+m)!(n−m)!n(n+1)[zn(kr)]2∫02π∫0π𝐍oemn⋅𝐍oemnsin⁡ϑdϑdφ=(1+δm,0)2π(2n+1)2(n+m)!(n−m)!n(n+1){(n+1)[zn−1(kr)]2+n[zn+1(kr)]2}∫0π∫02π𝐋oemn⋅𝐍oemnsin⁡ϑdϑdφ=(1+δm,0)2π(2n+1)2(n+m)!(n−m)!n(n+1)k{[zn−1(kr)]2−[zn+1(kr)]2}

All other integrals over the angles between different functions or functions with different indices are equal to zero.

Rotation and inversion

[edit]
Illustration of the transformation of vector spherical harmonics under rotations. One can see that they are transformed in the same way as the corresponding scalar functions.

Under rotation, vector spherical harmonics are transformed through each other in the same way as the corresponding scalar spherical functions, which are generating for a specific type of vector harmonics. For example, if the generating functions are the usual spherical harmonics, then the vector harmonics will also be transformed through the Wigner D-matrices[9][10][11] D̂(α,β,γ)𝐘JM(s)(θ,φ)=∑M′=−JJ[DMM′(J)(α,β,γ)]∗𝐘JM′(s)(θ,φ), The behavior under rotations is the same for electrical, magnetic and longitudinal harmonics.

Under inversion, electric and longitudinal spherical harmonics behave in the same way as scalar spherical functions, i.e. Î𝐍JM(θ,φ)=(−1)J𝐍JM(θ,φ), and magnetic ones have the opposite parity: Î𝐌JM(θ,φ)=(−1)J+1𝐌JM(θ,φ),

Fluid dynamics

[edit]

In the calculation of the Stokes' law for the drag that a viscous fluid exerts on a small spherical particle, the velocity distribution obeys Navier–Stokes equations neglecting inertia, i.e.,

0=∇⋅𝐯,𝟎=−∇p+η∇2𝐯,

with the boundary conditions

𝐯={𝟎r=a,−𝐔0r→∞.

where U is the relative velocity of the particle to the fluid far from the particle. In spherical coordinates this velocity at infinity can be written as

𝐔0=U0(cos⁡θ𝐫̂−sin⁡θ𝜽̂)=U0(𝐘10+𝜳10).

The last expression suggests an expansion in spherical harmonics for the liquid velocity and the pressure

p=p(r)Y10,𝐯=vr(r)𝐘10+v(1)(r)𝜳10.

Substitution in the Navier–Stokes equations produces a set of ordinary differential equations for the coefficients.

Integral relations

[edit]

Here the following definitions are used:

Yemn=cos⁡mφPnm(cos⁡θ)Yomn=sin⁡mφPnm(cos⁡θ)

𝐗eomn(𝐤k)=∇×(𝐤Yoemn(𝐤k))

𝐙oemn(𝐤k)=i𝐤k×𝐗eomn(𝐤k) In case, when instead of zn are spherical Bessel functions, with help of plane wave expansion one can obtain the following integral relations:[12]

𝐍pmn(1)(k,𝐫)=i−n4π∫𝐙pmn(𝐤k)ei𝐤⋅𝐫dΩk

𝐌pmn(1)(k,𝐫)=i−n4π∫𝐗pmn(𝐤k)ei𝐤⋅𝐫dΩk

In case, when zn are spherical Hankel functions, one should use the different formulae.[13][12] For vector spherical harmonics the following relations are obtained:

𝐌pmn(3)(k,𝐫)=i−n2πk∬−∞∞dk‖ei(kxx+kyy±kzz)kz𝐗pmn(𝐤k)

𝐍pmn(3)(k,𝐫)=i−n2πk∬−∞∞dk‖ei(kxx+kyy±kzz)kz𝐙pmn(𝐤k) where kz=k2−kx2−ky2, index (3) means, that spherical Hankel functions are used.

See also

[edit]

References

[edit]
  1. ↑ Barrera, R G; Estevez, G A; Giraldo, J (1985-10-01). "Vector spherical harmonics and their application to magnetostatics". European Journal of Physics. 6 (4). IOP Publishing: 287–294. Bibcode:1985EJPh....6..287B. doi:10.1088/0143-0807/6/4/014. ISSN 0143-0807. S2CID 250894245.
  2. ↑ Carrascal, B; Estevez, G A; Lee, Peilian; Lorenzo, V (1991-07-01). "Vector spherical harmonics and their application to classical electrodynamics". European Journal of Physics. 12 (4). IOP Publishing: 184–191. Bibcode:1991EJPh...12..184C. doi:10.1088/0143-0807/12/4/007. ISSN 0143-0807. S2CID 250886412.
  3. ↑ Hill, E. L. (1954). "The Theory of Vector Spherical Harmonics" (PDF). American Journal of Physics. 22 (4). American Association of Physics Teachers (AAPT): 211–214. Bibcode:1954AmJPh..22..211H. doi:10.1119/1.1933682. ISSN 0002-9505. S2CID 124182424. Archived from the original (PDF) on 2020-04-12.
  4. ↑ Weinberg, Erick J. (1994-01-15). "Monopole vector spherical harmonics". Physical Review D. 49 (2). American Physical Society (APS): 1086–1092. arXiv:hep-th/9308054. Bibcode:1994PhRvD..49.1086W. doi:10.1103/physrevd.49.1086. ISSN 0556-2821. PMID 10017069. S2CID 6429605.
  5. ↑ P.M. Morse and H. Feshbach, Methods of Theoretical Physics, Part II, New York: McGraw-Hill, 1898-1901 (1953)
  6. ↑ Cohen-Tannoudji, Claude; Diu, Bernard; Laloë, Franck (2020). "Complement B-XIX: Angular momentum of radiation". Quantum mechanics. Vol. III (2nd ed.). Weinheim, Germany: Wiley-VCH Verlag GmbH & Co. p. 2053. ISBN 978-3-527-34553-3.
  7. ↑ Bohren, Craig F. and Donald R. Huffman, Absorption and scattering of light by small particles, New York : Wiley, 1998, 530 p., ISBN 0-471-29340-7, ISBN 978-0-471-29340-8 (second edition)
  8. 1 2 Stratton, J. A. (1941). Electromagnetic Theory. New York: McGraw-Hill.
  9. ↑ D. A. Varhalovich, A. N. Moskalev, and V. K. Khersonskii, Quantum Theory of Angular Momentum [in Russian], Nauka, Leningrad (1975)
  10. ↑ Zhang, Huayong; Han, Yiping (2008). "Addition theorem for the spherical vector wave functions and its application to the beam shape coefficients". J. Opt. Soc. Am. B. 25 (2): 255–260. Bibcode:2008JOSAB..25..255Z. doi:10.1364/JOSAB.25.000255.
  11. ↑ Stein, Seymour (1961). "Addition theorems for spherical wave functions". Quarterly of Applied Mathematics. 19 (1): 15–24. doi:10.1090/qam/120407.
  12. 1 2 Stout, B. (2012). Popov, E (ed.). "Spherical harmonic lattice sums for gratings" (PDF). Institut Fresnel, Universite d'Aix-Marseille 6. Gratings: theory and numeric applications.
  13. ↑ Wittmann, R. C. (1988). "Spherical wave operators and the translation formulas". IEEE Transactions on Antennas and Propagation. 36 (8): 1078–1087. Bibcode:1988ITAP...36.1078W. doi:10.1109/8.7220.
[edit]