The Multipole Expansion
( \newcommand{\kernel}{\mathrm{null}\,}\)
A multipole expansion is a mathematical series representing a function that depends on angles—usually the two angles on a sphere. These series are useful because they can often be truncated, meaning that only the first few terms need to be retained for a good approximation to the original function. Multipole expansions are very frequently used 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.
The multipole expansion is expressed as a sum of terms with progressively finer angular features. For example, the initial term—called the zeroth, or monopole, moment—is a constant, independent of angle. The following term—the first, or dipole, moment—varies once from positive to negative around the sphere. Higher-order terms (like the quadrupole and octupole) vary more quickly with angles. A multipole moment usually involves powers (or inverse powers) of the distance to the origin, as well as some angular dependence.
Setting up the System
Consider an arbitrary charge distribution ρ(r′). We wish to find the electrostatic potential due to this charge distribution at a given point r. We assume that this point is at a large distance from the charge distribution, that is if r′ varies over the charge distribution, then r≫r′.
Now, the Coulomb potential for an arbitrary charge distribution is given by
V(r)=14πϵ0∫V′ρ(r′)|r−r′|dV′
Here,
|r−r′|=√|r2−2r⋅r′+r′2|=r√|1−2ˆr⋅r′r+(r′r)2|
where ˆr=r/r
Thus, using the fact that r is much larger than r′, we can write
1|r−r′|=1r1√|1−2ˆr⋅r′r+(r′r)2|
and using the binomial expansion (see below),
1√|1−2ˆr⋅r′r+(r′r)2|=1−ˆr⋅r′r+12r2(r′2−3(ˆr⋅r′)2)+O(r′r)3
(we neglect the third and higher order terms).
Binomial Theorem
The binomial theory can be used to expand specific functions into an infinite series:
(1+x)s=∞∑n=0s!n!(s−n)!xn=1+s1!x+s(s−1)2!x2+s(s−1)(s−2)3!x3+…
Equation ??? can be rewritten as
1|r−r′|=1r1√1+ϵ
where
ϵ=−2ˆr⋅r′r+(r′r)2
Applying the Binomial Theorem to Equation ??? (with s=−1/2) results in
1|r−r′|=1r(1−12ϵ+38ϵ2−516ϵ3+…)
Equation ??? originates from substituting Equation ??? into Equation ???.
The Expansion
Inserting Equation ??? into Equation ??? shows that the potential can be written as
V(r)=14πϵ0r∫V′ρ(r′)(1−ˆr⋅r′r+12r2(3(ˆr⋅r′)2−r′2)+O(r′r)3)dV′
We write this as
V(r)=Vmon(r)+Vdip(r)+Vquad(r)+…
The first (the zeroth-order) term in the expansion is called the monopole moment, the second (the first-order) term is called the dipole moment, the third (the second-order) is called the quadrupole moment, the fourth (third-order) term is called the octupole moment, and the fifth (fourth-order) term is called the hexadecapole moment. 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 (i.e., dotriacontapole) and 64-pole (hexacontatetrapole).
These moments can be expanded thusly
Vmon(r)=14πϵ0r∫V′ρ(r′)dV′Vdip(r)=−14πϵ0r2∫V′ρ(r′)(ˆr⋅r′)dV′Vquad(r)=18πϵ0r3∫V′ρ(r′)(3(ˆr⋅r′)2−r′2)dV′
and so on.
In principle, a multipole expansion provides an exact description of the potential and generally converges under two conditions:
- 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
- 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.
The Monopole Scalar
Observe that
Vmon(r)=14πϵ0r∫V′ρ(r′)dV′=q4πϵ0r
is a scalar, (actually the total charge in the distribution) and is called the electric monopole. This term indicates point charge electrical potential with charge q.
The Dipole Vector
If a charge distribution has a net total charge, it will tend to look like a monopole (point charge) from large distances. We can write
Vdip(r)=−ˆr4πϵ0r2⋅∫V′ρ(r′)r′dV′
The vector
p=∫V′ρ(r′)r′dV′
is called the electric dipole. And its magnitude is called the dipole moment of the charge distribution. This terms indicates the linear charge distribution geometry of a dipole electrical potential.
The Quadrupole Tensor
Let ˆr and r′ be expressed in Cartesian coordinates as (r1,r2,r3) and (x1,x2,x3). Then, (ˆr⋅r′)2=(rixi)2=rirjxixj.
We define a dyad to be the tensor ˆrˆr given by
(ˆrˆr)ij=rirj
Define the Quadrupole tensor as
T=∫V′ρ(r′)(3(r′r′)−Ir′2)dV′
Then, we can write Vqua as the tensor contraction
Vqua(r)=−ˆrˆr4πϵ0r3::T
this term indicates the three dimensional distribution of a quadruple electrical potential.