Multipole expansion

Last updated

A multipole expansion is a mathematical series representing a function that depends on angles—usually the two angles used in the spherical coordinate system (the polar and azimuthal angles) for three-dimensional Euclidean space, ${\displaystyle \mathbb {R} ^{3}}$. Similarly to Taylor series, multipole expansions are useful because oftentimes only the first few terms are needed to provide a good approximation of the original function. The function being expanded may be real or complex-valued and is defined either on ${\displaystyle \mathbb {R} ^{3}}$, or less often on ${\displaystyle \mathbb {R} ^{n}}$ for some other ${\displaystyle n}$.

Contents

Multipole expansions are used frequently in the study of electromagnetic and gravitational fields, where the fields at distant points are given in terms of sources in a small region. The multipole expansion with angles is often combined with an expansion in radius. Such a combination gives an expansion describing a function throughout three-dimensional space. [1]

The multipole expansion is expressed as a sum of terms with progressively finer angular features (moments). The first (the zeroth-order) term is called the monopole moment, the second (the first-order) term is called the dipole moment, the third (the second-order) the quadrupole moment, the fourth (third-order) term is called the octupole moment, and so on. Given the limitation of Greek numeral prefixes, terms of higher order are conventionally named by adding "-pole" to the number of poles—e.g., 32-pole (rarely dotriacontapole or triacontadipole) and 64-pole (rarely tetrahexacontapole or hexacontatetrapole). [2] [3] [4] A multipole moment usually involves powers (or inverse powers) of the distance to the origin, as well as some angular dependence.

In principle, a multipole expansion provides an exact description of the potential, and generally converges under two conditions: (1) if the sources (e.g. charges) are localized close to the origin and the point at which the potential is observed is far from the origin; or (2) the reverse, i.e., if the sources are located far from the origin and the potential is observed close to the origin. In the first (more common) case, the coefficients of the series expansion are called exterior multipole moments or simply multipole moments whereas, in the second case, they are called interior multipole moments.

Expansion in spherical harmonics

Most commonly, the series is written as a sum of spherical harmonics. Thus, we might write a function ${\displaystyle f(\theta ,\varphi )}$ as the sum

${\displaystyle f(\theta ,\varphi )=\sum _{\ell =0}^{\infty }\,\sum _{m=-\ell }^{\ell }\,C_{\ell }^{m}\,Y_{\ell }^{m}(\theta ,\varphi )}$

where ${\displaystyle Y_{\ell }^{m}(\theta ,\varphi )}$ are the standard spherical harmonics, and ${\displaystyle C_{\ell }^{m}}$ are constant coefficients which depend on the function. The term ${\displaystyle C_{0}^{0}}$ represents the monopole; ${\displaystyle C_{1}^{-1},C_{1}^{0},C_{1}^{1}}$ represent the dipole; and so on. Equivalently, the series is also frequently written [5] as

${\displaystyle f(\theta ,\varphi )=C+C_{i}n^{i}+C_{ij}n^{i}n^{j}+C_{ijk}n^{i}n^{j}n^{k}+C_{ijk\ell }n^{i}n^{j}n^{k}n^{\ell }+\dots }$

where the ${\displaystyle n^{i}}$ represent the components of a unit vector in the direction given by the angles ${\displaystyle \theta }$ and ${\displaystyle \varphi }$, and indices are implicitly summed. Here, the term ${\displaystyle C}$ is the monopole; ${\displaystyle C_{i}}$ is a set of three numbers representing the dipole; and so on.

In the above expansions, the coefficients may be real or complex. If the function being expressed as a multipole expansion is real, however, the coefficients must satisfy certain properties. In the spherical harmonic expansion, we must have

${\displaystyle C_{\ell }^{-m}=(-1)^{m}C_{\ell }^{m\ast }\,.}$

In the multi-vector expansion, each coefficient must be real:

${\displaystyle C=C^{\ast };\ C_{i}=C_{i}^{\ast };\ C_{ij}=C_{ij}^{\ast };\ C_{ijk}=C_{ijk}^{\ast };\ \ldots }$

While expansions of scalar functions are by far the most common application of multipole expansions, they may also be generalized to describe tensors of arbitrary rank. [6] This finds use in multipole expansions of the vector potential in electromagnetism, or the metric perturbation in the description of gravitational waves.

For describing functions of three dimensions, away from the coordinate origin, the coefficients of the multipole expansion can be written as functions of the distance to the origin, ${\displaystyle r}$—most frequently, as a Laurent series in powers of ${\displaystyle r}$. For example, to describe the electromagnetic potential, ${\displaystyle V}$, from a source in a small region near the origin, the coefficients may be written as:

${\displaystyle V(r,\theta ,\varphi )=\sum _{\ell =0}^{\infty }\,\sum _{m=-l}^{\ell }C_{\ell }^{m}(r)\,Y_{\ell }^{m}(\theta ,\varphi )=\sum _{j=1}^{\infty }\,\sum _{\ell =0}^{\infty }\,\sum _{m=-l}^{\ell }{\frac {D_{\ell ,j}^{m}}{r^{j}}}\,Y_{\ell }^{m}(\theta ,\varphi ).}$

Applications

Multipole expansions are widely used in problems involving gravitational fields of systems of masses, electric and magnetic fields of charge and current distributions, and the propagation of electromagnetic waves. A classic example is the calculation of the exterior multipole moments of atomic nuclei from their interaction energies with the interior multipoles of the electronic orbitals. The multipole moments of the nuclei report on the distribution of charges within the nucleus and, thus, on the shape of the nucleus. Truncation of the multipole expansion to its first non-zero term is often useful for theoretical calculations.

Multipole expansions are also useful in numerical simulations, and form the basis of the fast multipole method of Greengard and Rokhlin, a general technique for efficient computation of energies and forces in systems of interacting particles. The basic idea is to decompose the particles into groups; particles within a group interact normally (i.e., by the full potential), whereas the energies and forces between groups of particles are calculated from their multipole moments. The efficiency of the fast multipole method is generally similar to that of Ewald summation, but is superior if the particles are clustered, i.e. the system has large density fluctuations.

Multipole expansion of a potential outside an electrostatic charge distribution

Consider a discrete charge distribution consisting of N point charges qi with position vectors ri. We assume the charges to be clustered around the origin, so that for all i: ri<rmax, where rmax has some finite value. The potential V(R), due to the charge distribution, at a point R outside the charge distribution, i.e., |R|>rmax, can be expanded in powers of 1/R. Two ways of making this expansion can be found in the literature: The first is a Taylor series in the Cartesian coordinates x, y, and z, while the second is in terms of spherical harmonics which depend on spherical polar coordinates. The Cartesian approach has the advantage that no prior knowledge of Legendre functions, spherical harmonics, etc., is required. Its disadvantage is that the derivations are fairly cumbersome (in fact a large part of it is the implicit rederivation of the Legendre expansion of 1/|rR|, which was done once and for all by Legendre in the 1780s). Also it is difficult to give a closed expression for a general term of the multipole expansionusually only the first few terms are given followed by an ellipsis.

Expansion in Cartesian coordinates

Let ${\displaystyle v}$ satisfy ${\displaystyle v(x)=v(-x)}$. Then the Taylor expansion of v(rR) around the origin r = 0 can be written

${\displaystyle v(\mathbf {r} -\mathbf {R} )=v(\mathbf {R} )-\sum _{\alpha =x,y,z}r_{\alpha }v_{\alpha }(\mathbf {R} )+{\frac {1}{2}}\sum _{\alpha =x,y,z}\sum _{\beta =x,y,z}r_{\alpha }r_{\beta }v_{\alpha \beta }(\mathbf {R} )-\ldots +\dots }$

with

${\displaystyle v_{\alpha }(\mathbf {R} )\equiv \left({\frac {\partial v(\mathbf {r} -\mathbf {R} )}{\partial r_{\alpha }}}\right)_{\mathbf {r} =\mathbf {0} }\quad {\hbox{and}}\quad v_{\alpha \beta }(\mathbf {R} )\equiv \left({\frac {\partial ^{2}v(\mathbf {r} -\mathbf {R} )}{\partial r_{\alpha }\partial r_{\beta }}}\right)_{\mathbf {r} =\mathbf {0} }.}$

If v(rR) satisfies the Laplace equation

${\displaystyle \left(\nabla ^{2}v(\mathbf {r} -\mathbf {R} )\right)_{\mathbf {r} =\mathbf {0} }=\sum _{\alpha =x,y,z}v_{\alpha \alpha }(\mathbf {R} )=0}$

then the expansion can be rewritten in terms of the components of a traceless Cartesian second rank tensor:

${\displaystyle \sum _{\alpha =x,y,z}\sum _{\beta =x,y,z}r_{\alpha }r_{\beta }v_{\alpha \beta }(\mathbf {R} )={\frac {1}{3}}\sum _{\alpha =x,y,z}\sum _{\beta =x,y,z}(3r_{\alpha }r_{\beta }-\delta _{\alpha \beta }r^{2})v_{\alpha \beta }(\mathbf {R} ),}$

where δαβ is the Kronecker delta and r2|r|2. Removing the trace is common, because it takes the rotationally invariant r2 out of the second rank tensor.

Example

Consider now the following form of v(rR):

${\displaystyle v(\mathbf {r} -\mathbf {R} )\equiv {\frac {1}{|\mathbf {r} -\mathbf {R} |}}.}$

Then by direct differentiation it follows that

${\displaystyle v(\mathbf {R} )={\frac {1}{R}},\quad v_{\alpha }(\mathbf {R} )=-{\frac {R_{\alpha }}{R^{3}}},\quad {\hbox{and}}\quad v_{\alpha \beta }(\mathbf {R} )={\frac {3R_{\alpha }R_{\beta }-\delta _{\alpha \beta }R^{2}}{R^{5}}}.}$

Define a monopole, dipole, and (traceless) quadrupole by, respectively,

${\displaystyle q_{\mathrm {tot} }\equiv \sum _{i=1}^{N}q_{i},\quad P_{\alpha }\equiv \sum _{i=1}^{N}q_{i}r_{i\alpha },\quad {\hbox{and}}\quad Q_{\alpha \beta }\equiv \sum _{i=1}^{N}q_{i}(3r_{i\alpha }r_{i\beta }-\delta _{\alpha \beta }r_{i}^{2}),}$

and we obtain finally the first few terms of the multipole expansion of the total potential, which is the sum of the Coulomb potentials of the separate charges: [7] :137138

{\displaystyle {\begin{aligned}4\pi \varepsilon _{0}V(\mathbf {R} )&\equiv \sum _{i=1}^{N}q_{i}v(\mathbf {r} _{i}-\mathbf {R} )\\&={\frac {q_{\mathrm {tot} }}{R}}+{\frac {1}{R^{3}}}\sum _{\alpha =x,y,z}P_{\alpha }R_{\alpha }+{\frac {1}{2R^{5}}}\sum _{\alpha ,\beta =x,y,z}Q_{\alpha \beta }R_{\alpha }R_{\beta }+\dots \end{aligned}}}

This expansion of the potential of a discrete charge distribution is very similar to the one in real solid harmonics given below. The main difference is that the present one is in terms of linearly dependent quantities, for

${\displaystyle \sum _{\alpha }v_{\alpha \alpha }=0\quad {\hbox{and}}\quad \sum _{\alpha }Q_{\alpha \alpha }=0.}$

NOTE: If the charge distribution consists of two charges of opposite sign which are an infinitesimal distance d apart, so that d/R ≫ (d/R)2, it is easily shown that the only non-vanishing term in the expansion is

${\displaystyle V(\mathbf {R} )={\frac {1}{4\pi \varepsilon _{0}R^{3}}}(\mathbf {P} \cdot \mathbf {R} ),}$

the electric dipolar potential field.

Spherical form

The potential V(R) at a point R outside the charge distribution, i.e. |R| > rmax, can be expanded by the Laplace expansion:

${\displaystyle V(\mathbf {R} )\equiv \sum _{i=1}^{N}{\frac {q_{i}}{4\pi \varepsilon _{0}|\mathbf {r} _{i}-\mathbf {R} |}}={\frac {1}{4\pi \varepsilon _{0}}}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }(-1)^{m}I_{\ell }^{-m}(\mathbf {R} )\sum _{i=1}^{N}q_{i}R_{\ell }^{m}(\mathbf {r} _{i}),}$

where ${\displaystyle I_{\ell }^{-m}(\mathbf {R} )}$ is an irregular solid harmonic (defined below as a spherical harmonic function divided by ${\displaystyle R^{\ell +1}}$) and ${\displaystyle R_{\ell }^{m}(\mathbf {r} )}$ is a regular solid harmonic (a spherical harmonic times r). We define the spherical multipole moment of the charge distribution as follows

${\displaystyle Q_{\ell }^{m}\equiv \sum _{i=1}^{N}q_{i}R_{\ell }^{m}(\mathbf {r} _{i}),\quad \ -\ell \leq m\leq \ell .}$

Note that a multipole moment is solely determined by the charge distribution (the positions and magnitudes of the N charges).

A spherical harmonic depends on the unit vector ${\displaystyle {\hat {R}}}$. (A unit vector is determined by two spherical polar angles.) Thus, by definition, the irregular solid harmonics can be written as

${\displaystyle I_{\ell }^{m}(\mathbf {R} )\equiv {\sqrt {\frac {4\pi }{2\ell +1}}}{\frac {Y_{\ell }^{m}({\hat {R}})}{R^{\ell +1}}}}$

so that the multipole expansion of the field V(R) at the point R outside the charge distribution is given by

{\displaystyle {\begin{aligned}V(\mathbf {R} )&={\frac {1}{4\pi \varepsilon _{0}}}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }(-1)^{m}I_{\ell }^{-m}(\mathbf {R} )Q_{\ell }^{m}\\&={\frac {1}{4\pi \varepsilon _{0}}}\sum _{\ell =0}^{\infty }\left[{\frac {4\pi }{2\ell +1}}\right]^{1/2}\;{\frac {1}{R^{\ell +1}}}\;\sum _{m=-\ell }^{\ell }(-1)^{m}Y_{\ell }^{-m}({\hat {R}})Q_{\ell }^{m},\qquad R>r_{\mathrm {max} }.\end{aligned}}}

This expansion is completely general in that it gives a closed form for all terms, not just for the first few. It shows that the spherical multipole moments appear as coefficients in the 1/R expansion of the potential.

It is of interest to consider the first few terms in real form, which are the only terms commonly found in undergraduate textbooks. Since the summand of the m summation is invariant under a unitary transformation of both factors simultaneously and since transformation of complex spherical harmonics to real form is by a unitary transformation, we can simply substitute real irregular solid harmonics and real multipole moments. The  =0 term becomes

${\displaystyle V_{\ell =0}(\mathbf {R} )={\frac {q_{\mathrm {tot} }}{4\pi \varepsilon _{0}R}}\quad {\hbox{with}}\quad q_{\mathrm {tot} }\equiv \sum _{i=1}^{N}q_{i}.}$

This is in fact Coulomb's law again. For the  =1 term we introduce

${\displaystyle \mathbf {R} =(R_{x},R_{y},R_{z}),\quad \mathbf {P} =(P_{x},P_{y},P_{z})\quad {\hbox{with}}\quad P_{\alpha }\equiv \sum _{i=1}^{N}q_{i}r_{i\alpha },\quad \alpha =x,y,z.}$

Then

${\displaystyle V_{\ell =1}(\mathbf {R} )={\frac {1}{4\pi \varepsilon _{0}R^{3}}}(R_{x}P_{x}+R_{y}P_{y}+R_{z}P_{z})={\frac {\mathbf {R} \cdot \mathbf {P} }{4\pi \varepsilon _{0}R^{3}}}={\frac {{\hat {R}}\cdot \mathbf {P} }{4\pi \varepsilon _{0}R^{2}}}.}$

This term is identical to the one found in Cartesian form.

In order to write the  =2 term, we have to introduce shorthand notations for the five real components of the quadrupole moment and the real spherical harmonics. Notations of the type

${\displaystyle Q_{z^{2}}\equiv \sum _{i=1}^{N}q_{i}\;{\frac {1}{2}}(3z_{i}^{2}-r_{i}^{2}),}$

can be found in the literature. Clearly the real notation becomes awkward very soon, exhibiting the usefulness of the complex notation.

Interaction of two non-overlapping charge distributions

Consider two sets of point charges, one set {qi} clustered around a point A and one set {qj} clustered around a point B. Think for example of two molecules, and recall that a molecule by definition consists of electrons (negative point charges) and nuclei (positive point charges). The total electrostatic interaction energy UAB between the two distributions is

${\displaystyle U_{AB}=\sum _{i\in A}\sum _{j\in B}{\frac {q_{i}q_{j}}{4\pi \varepsilon _{0}r_{ij}}}.}$

This energy can be expanded in a power series in the inverse distance of A and B. This expansion is known as the multipole expansion of UAB.

In order to derive this multipole expansion, we write rXY = rYrX, which is a vector pointing from X towards Y. Note that

${\displaystyle \mathbf {R} _{AB}+\mathbf {r} _{Bj}+\mathbf {r} _{ji}+\mathbf {r} _{iA}=0\quad \iff \quad \mathbf {r} _{ij}=\mathbf {R} _{AB}-\mathbf {r} _{Ai}+\mathbf {r} _{Bj}.}$

We assume that the two distributions do not overlap:

${\displaystyle |\mathbf {R} _{AB}|>|\mathbf {r} _{Bj}-\mathbf {r} _{Ai}|{\text{ for all }}i,j.}$

Under this condition we may apply the Laplace expansion in the following form

${\displaystyle {\frac {1}{|\mathbf {r} _{j}-\mathbf {r} _{i}|}}={\frac {1}{|\mathbf {R} _{AB}-(\mathbf {r} _{Ai}-\mathbf {r} _{Bj})|}}=\sum _{L=0}^{\infty }\sum _{M=-L}^{L}\,(-1)^{M}I_{L}^{-M}(\mathbf {R} _{AB})\;R_{L}^{M}(\mathbf {r} _{Ai}-\mathbf {r} _{Bj}),}$

where ${\displaystyle I_{L}^{M}}$ and ${\displaystyle R_{L}^{M}}$ are irregular and regular solid harmonics, respectively. The translation of the regular solid harmonic gives a finite expansion,

${\displaystyle R_{L}^{M}(\mathbf {r} _{Ai}-\mathbf {r} _{Bj})=\sum _{\ell _{A}=0}^{L}(-1)^{L-\ell _{A}}{\binom {2L}{2\ell _{A}}}^{1/2}}$
${\displaystyle \times \sum _{m_{A}=-\ell _{A}}^{\ell _{A}}R_{\ell _{A}}^{m_{A}}(\mathbf {r} _{Ai})R_{L-\ell _{A}}^{M-m_{A}}(\mathbf {r} _{Bj})\;\langle \ell _{A},m_{A};L-\ell _{A},M-m_{A}\mid LM\rangle ,}$

where the quantity between pointed brackets is a Clebsch–Gordan coefficient. Further we used

${\displaystyle R_{\ell }^{m}(-\mathbf {r} )=(-1)^{\ell }R_{\ell }^{m}(\mathbf {r} ).}$

Use of the definition of spherical multipoles Qm
and covering of the summation ranges in a somewhat different order (which is only allowed for an infinite range of L) gives finally

{\displaystyle {\begin{aligned}U_{AB}={}&{\frac {1}{4\pi \varepsilon _{0}}}\sum _{\ell _{A}=0}^{\infty }\sum _{\ell _{B}=0}^{\infty }(-1)^{\ell _{B}}{\binom {2\ell _{A}+2\ell _{B}}{2\ell _{A}}}^{1/2}\\[5pt]&\times \sum _{m_{A}=-\ell _{A}}^{\ell _{A}}\sum _{m_{B}=-\ell _{B}}^{\ell _{B}}(-1)^{m_{A}+m_{B}}I_{\ell _{A}+\ell _{B}}^{-m_{A}-m_{B}}(\mathbf {R} _{AB})\;Q_{\ell _{A}}^{m_{A}}Q_{\ell _{B}}^{m_{B}}\;\langle \ell _{A},m_{A};\ell _{B},m_{B}\mid \ell _{A}+\ell _{B},m_{A}+m_{B}\rangle .\end{aligned}}}

This is the multipole expansion of the interaction energy of two non-overlapping charge distributions which are a distance RAB apart. Since

${\displaystyle I_{\ell _{A}+\ell _{B}}^{-(m_{A}+m_{B})}(\mathbf {R} _{AB})\equiv \left[{\frac {4\pi }{2\ell _{A}+2\ell _{B}+1}}\right]^{1/2}\;{\frac {Y_{\ell _{A}+\ell _{B}}^{-(m_{A}+m_{B})}\left({\widehat {\mathbf {R} }}_{AB}\right)}{R_{AB}^{\ell _{A}+\ell _{B}+1}}},}$

this expansion is manifestly in powers of 1/RAB. The function Yml is a normalized spherical harmonic.

Molecular moments

All atoms and molecules (except S-state atoms) have one or more non-vanishing permanent multipole moments. Different definitions can be found in the literature, but the following definition in spherical form has the advantage that it is contained in one general equation. Because it is in complex form it has as the further advantage that it is easier to manipulate in calculations than its real counterpart.

We consider a molecule consisting of N particles (electrons and nuclei) with charges eZi. (Electrons have a Z-value of −1, while for nuclei it is the atomic number). Particle i has spherical polar coordinates ri, θi, and φi and Cartesian coordinates xi, yi, and zi. The (complex) electrostatic multipole operator is

${\displaystyle Q_{\ell }^{m}\equiv \sum _{i=1}^{N}eZ_{i}\;R_{\ell }^{m}(\mathbf {r} _{i}),}$

where ${\displaystyle R_{\ell }^{m}(\mathbf {r} _{i})}$ is a regular solid harmonic function in Racah's normalization (also known as Schmidt's semi-normalization). If the molecule has total normalized wave function Ψ (depending on the coordinates of electrons and nuclei), then the multipole moment of order ${\displaystyle \ell }$ of the molecule is given by the expectation (expected) value:

${\displaystyle M_{\ell }^{m}\equiv \langle \Psi \mid Q_{\ell }^{m}\mid \Psi \rangle .}$

If the molecule has certain point group symmetry, then this is reflected in the wave function: Ψ transforms according to a certain irreducible representation λ of the group ("Ψ has symmetry type λ"). This has the consequence that selection rules hold for the expectation value of the multipole operator, or in other words, that the expectation value may vanish because of symmetry. A well-known example of this is the fact that molecules with an inversion center do not carry a dipole (the expectation values of ${\displaystyle Q_{1}^{m}}$ vanish for m = −1, 0, 1). For a molecule without symmetry, no selection rules are operative and such a molecule will have non-vanishing multipoles of any order (it will carry a dipole and simultaneously a quadrupole, octupole, hexadecapole, etc.).

The lowest explicit forms of the regular solid harmonics (with the Condon-Shortley phase) give:

${\displaystyle M_{0}^{0}=\sum _{i=1}^{N}eZ_{i},}$

(the total charge of the molecule). The (complex) dipole components are:

${\displaystyle M_{1}^{1}=-{\sqrt {\tfrac {1}{2}}}\sum _{i=1}^{N}eZ_{i}\langle \Psi |x_{i}+iy_{i}|\Psi \rangle \quad {\hbox{and}}\quad M_{1}^{-1}={\sqrt {\tfrac {1}{2}}}\sum _{i=1}^{N}eZ_{i}\langle \Psi |x_{i}-iy_{i}|\Psi \rangle .}$
${\displaystyle M_{1}^{0}=\sum _{i=1}^{N}eZ_{i}\langle \Psi |z_{i}|\Psi \rangle .}$

Note that by a simple linear combination one can transform the complex multipole operators to real ones. The real multipole operators are of cosine type ${\displaystyle C_{\ell }^{m}}$ or sine type ${\displaystyle S_{\ell }^{m}}$. A few of the lowest ones are:

{\displaystyle {\begin{aligned}C_{1}^{0}&=\sum _{i=1}^{N}eZ_{i}\;z_{i}\\C_{1}^{1}&=\sum _{i=1}^{N}eZ_{i}\;x_{i}\\S_{1}^{1}&=\sum _{i=1}^{N}eZ_{i}\;y_{i}\\C_{2}^{0}&={\frac {1}{2}}\sum _{i=1}^{N}eZ_{i}\;(3z_{i}^{2}-r_{i}^{2})\\C_{2}^{1}&={\sqrt {3}}\sum _{i=1}^{N}eZ_{i}\;z_{i}x_{i}\\C_{2}^{2}&={\frac {1}{3}}{\sqrt {3}}\sum _{i=1}^{N}eZ_{i}\;(x_{i}^{2}-y_{i}^{2})\\S_{2}^{1}&={\sqrt {3}}\sum _{i=1}^{N}eZ_{i}\;z_{i}y_{i}\\S_{2}^{2}&={\frac {2}{3}}{\sqrt {3}}\sum _{i=1}^{N}eZ_{i}\;x_{i}y_{i}\end{aligned}}}

Note on conventions

The definition of the complex molecular multipole moment given above is the complex conjugate of the definition given in this article, which follows the definition of the standard textbook on classical electrodynamics by Jackson, [7] :137 except for the normalization. Moreover, in the classical definition of Jackson the equivalent of the N-particle quantum mechanical expectation value is an integral over a one-particle charge distribution. Remember that in the case of a one-particle quantum mechanical system the expectation value is nothing but an integral over the charge distribution (modulus of wavefunction squared), so that the definition of this article is a quantum mechanical N-particle generalization of Jackson's definition.

The definition in this article agrees with, among others, the one of Fano and Racah [8] and Brink and Satchler. [9]

Examples

There are many types of multipole moments, since there are many types of potentials and many ways of approximating a potential by a series expansion, depending on the coordinates and the symmetry of the charge distribution. The most common expansions include:

Examples of 1/R potentials include the electric potential, the magnetic potential and the gravitational potential of point sources. An example of a lnR potential is the electric potential of an infinite line charge.

General mathematical properties

Multipole moments in mathematics and mathematical physics form an orthogonal basis for the decomposition of a function, based on the response of a field to point sources that are brought infinitely close to each other. These can be thought of as arranged in various geometrical shapes, or, in the sense of distribution theory, as directional derivatives.

Multipole expansions are related to the underlying rotational symmetry of the physical laws and their associated differential equations. Even though the source terms (such as the masses, charges, or currents) may not be symmetrical, one can expand them in terms of irreducible representations of the rotational symmetry group, which leads to spherical harmonics and related sets of orthogonal functions. One uses the technique of separation of variables to extract the corresponding solutions for the radial dependencies.

In practice, many fields can be well approximated with a finite number of multipole moments (although an infinite number may be required to reconstruct a field exactly). A typical application is to approximate the field of a localized charge distribution by its monopole and dipole terms. Problems solved once for a given order of multipole moment may be linearly combined to create a final approximate solution for a given source.

Related Research Articles

In mathematics and physics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace who first studied its properties. This is often written as

The moment of inertia, otherwise known as the mass moment of inertia, angular mass, second moment of mass, or most accurately, rotational inertia, of a rigid body is a quantity that determines the torque needed for a desired angular acceleration about a rotational axis, akin to how mass determines the force needed for a desired acceleration. It depends on the body's mass distribution and the axis chosen, with larger moments requiring more torque to change the body's rate of rotation.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space under the operation of composition. By definition, a rotation about the origin is a transformation that preserves the origin, Euclidean distance, and orientation. Every non-trivial rotation is determined by its axis of rotation and its angle of rotation. Composing two rotations results in another rotation; every rotation has a unique inverse rotation; and the identity map satisfies the definition of a rotation. Owing to the above properties, the set of all rotations is a group under composition. Rotations are not commutative, making it a nonabelian group. Moreover, the rotation group has a natural structure as a manifold for which the group operations are smoothly differentiable; so it is in fact a Lie group. It is compact and has dimension 3.

In mathematics and physical science, spherical harmonics are special functions defined on the surface of a sphere. They are often employed in solving partial differential equations in many scientific fields.

In physics, a moment is an expression involving the product of a distance and physical quantity, and in this way it accounts for how the physical quantity is located or arranged.

Slater-type orbitals (STOs) are functions used as atomic orbitals in the linear combination of atomic orbitals molecular orbital method. They are named after the physicist John C. Slater, who introduced them in 1930.

In atomic physics, hyperfine structure is defined by small shifts in otherwise degenerate energy levels and the resulting splittings in those energy levels of atoms, molecules, and ions, due to interaction between the nucleus and electron clouds.

The rigid rotor is a mechanical model of rotating systems. An arbitrary rigid rotor is a 3-dimensional rigid object, such as a top. To orient such an object in space requires three angles, known as Euler angles. A special rigid rotor is the linear rotor requiring only two angles to describe, for example of a diatomic molecule. More general molecules are 3-dimensional, such as water, ammonia, or methane. The rigid-rotor Schroedinger equation is discussed in Section 11.2 on pages 240-253 of the textbook by Bunker and Jensen.

A quadrupole or quadrapole is one of a sequence of configurations of things like electric charge or current, or gravitational mass that can exist in ideal form, but it is usually just part of a multipole expansion of a more complex structure reflecting various orders of complexity.

In mathematics, Gegenbauer polynomials or ultraspherical polynomialsC(α)
n
(x) are orthogonal polynomials on the interval [−1,1] with respect to the weight function (1 − x2)α–1/2. They generalize Legendre polynomials and Chebyshev polynomials, and are special cases of Jacobi polynomials. They are named after Leopold Gegenbauer.

The electromagnetic wave equation is a second-order partial differential equation that describes the propagation of electromagnetic waves through a medium or in a vacuum. It is a three-dimensional form of the wave equation. The homogeneous form of the equation, written in terms of either the electric field E or the magnetic field B, takes the form:

Spherical multipole moments are the coefficients in a series expansion of a potential that varies inversely with the distance R to a source, i.e., as 1/R. Examples of such potentials are the electric potential, the magnetic potential and the gravitational potential.

In quantum mechanics, Landau quantization refers to the quantization of the cyclotron orbits of charged particles in a uniform magnetic field. As a result, the charged particles can only occupy orbits with discrete, equidistant energy values, called Landau levels. These levels are degenerate, with the number of electrons per level directly proportional to the strength of the applied magnetic field. It is named after the Soviet physicist Lev Landau.

In physics, the Laplace expansion of potentials that are directly proportional to the inverse of the distance, such as Newton's gravitational potential or Coulomb's electrostatic potential, expresses them in terms of the spherical Legendre polynomials. In quantum mechanical calculations on atoms the expansion is used in the evaluation of integrals of the inter-electronic repulsion.

In physics and mathematics, the solid harmonics are solutions of the Laplace equation in spherical polar coordinates, assumed to be (smooth) functions . There are two kinds: the regular solid harmonics, which vanish at the origin and the irregular solid harmonics, which are singular at the origin. Both sets of functions play an important role in potential theory, and are obtained by rescaling spherical harmonics appropriately:

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.

In the mathematical study of rotational symmetry, the zonal spherical harmonics are special spherical harmonics that are invariant under the rotation through a particular fixed axis. The zonal spherical functions are a broad extension of the notion of zonal spherical harmonics to allow for a more general symmetry group.

In linear algebra, a raising or lowering operator is an operator that increases or decreases the eigenvalue of another operator. In quantum mechanics, the raising operator is sometimes called the creation operator, and the lowering operator the annihilation operator. Well-known applications of ladder operators in quantum mechanics are in the formalisms of the quantum harmonic oscillator and angular momentum.

In pure and applied mathematics, quantum mechanics and computer graphics, a tensor operator generalizes the notion of operators which are scalars and vectors. A special class of these are spherical tensor operators which apply the notion of the spherical basis and spherical harmonics. The spherical basis closely relates to the description of angular momentum in quantum mechanics and spherical harmonic functions. The coordinate-free generalization of a tensor operator is known as a representation operator.

References

1. Edmonds, A. R. (1960). . Princeton University Press.
2. Auzinsh, Marcis; Budker, Dmitry; Rochester, Simon (2010). Optically polarized atoms : understanding light-atom interactions. Oxford: New York. p. 100. ISBN   9780199565122.
3. Okumura, Mitchio; Chan, Man-Chor; Oka, Takeshi (2 January 1989). "High-resolution infrared spectroscopy of solid hydrogen: The tetrahexacontapole-induced transitions" (PDF). Physical Review Letters. 62 (1): 32–35. Bibcode:1989PhRvL..62...32O. doi:10.1103/PhysRevLett.62.32. PMID   10039541.
4. Ikeda, Hiroaki; Suzuki, Michi-To; Arita, Ryotaro; Takimoto, Tetsuya; Shibauchi, Takasada; Matsuda, Yuji (3 June 2012). "Emergent rank-5 nematic order in URu2Si2". Nature Physics. 8 (7): 528–533. arXiv:. Bibcode:2012NatPh...8..528I. doi:10.1038/nphys2330.
5. Thompson, William J. Angular Momentum. John Wiley & Sons, Inc.
6. Thorne, Kip S. (April 1980). "Multipole Expansions of Gravitational Radiation" (PDF). Reviews of Modern Physics. 52 (2): 299–339. Bibcode:1980RvMP...52..299T. doi:10.1103/RevModPhys.52.299.
7. Jackson, John David (1975). (2d ed.). New York: Wiley. ISBN   047143132X.
8. U. Fano and G. Racah, Irreducible Tensorial Sets, Academic Press, New York (1959). p. 31
9. D. M. Brink and G. R. Satchler, Angular Momentum, 2nd edition, Clarendon Press, Oxford, UK (1968). p. 64. See also footnote on p. 90.