Fact-checked by Grok 2 weeks ago

Multipole expansion

In physics, the multipole expansion is a mathematical series representation of the scalar or vector potential produced by a localized charge or current distribution, approximating the far-field behavior by decomposing it into successive terms of increasing angular order and decreasing radial dependence, starting from the monopole (ℓ=0), dipole (ℓ=1), quadrupole (ℓ=2), and higher multipoles.^[1] This expansion arises from the Taylor series development of the inverse distance in the potential integral, typically using Legendre polynomials or spherical harmonics to handle the angular dependence, allowing the potential V(R) at a distant point R from a compact source of size much smaller than R to be written as V(R) = Σ_{ℓ=0}^∞ [1/(4πε₀ R^{ℓ+1})] ∫ ρ(r') r'^ℓ P_ℓ(cosα) d³r', where ρ is the charge density and α is the angle between R and r'.^[2] The monopole term captures the net charge Q_net as V_monopole = Q_net / (4πε₀ R), vanishing for neutral systems, while the dipole term involves the dipole moment p = ∫ r' ρ(r') d³r' with V_dipole = (p · R̂) / (4πε₀ R²), and the quadrupole term uses the traceless tensor Q_{ij} = ∫ (3x_i' x_j' - r'^2 δ_{ij}) ρ(r') d³r' for V_quadrupole = [\hat{R} · Q · \hat{R}] / [2 (4πε₀ R³)].^[3] The technique is essential in electrostatics for simplifying calculations of potentials and fields from complex distributions, such as molecules or nuclei, where direct integration is impractical, and it enables far-field approximations by truncating at the lowest non-vanishing term for desired accuracy.^[4] Beyond electrostatics, multipole expansions extend to magnetostatics via magnetic multipole moments, gravitational potentials in general relativity, and quantum mechanics for atomic and nuclear interactions, with the expansion parameter typically being (source size / observation distance), ensuring convergence outside the source.^[2] Higher-order terms like octupoles (ℓ=3) describe finer asymmetries in the source, but their contributions diminish rapidly with distance, making the expansion a powerful tool for both analytical and numerical modeling in diverse physical contexts.^[1]

Mathematical Foundations

General Definition

The multipole expansion provides a series representation of scalar or vector potentials generated by localized source distributions at distant observation points, analogous to a Taylor series but featuring inverse powers of the distance r from the source origin to the field point. This method decomposes the potential into contributions from successively higher-order moments of the source, offering approximations that improve with increasing distance from the source relative to its spatial extent. It is fundamental in fields like electromagnetism and gravity for analyzing far-field behaviors without solving the full source distribution explicitly.^[2] In its general one-dimensional form, applicable along a line or axis of symmetry, the scalar potential \phi(\mathbf{r}) due to a source density \rho(\mathbf{r}') (with r > r') is expressed as

\phi(\mathbf{r}) = \sum_{n=0}^{\infty} \frac{1}{r^{n+1}} \int (r')^n P_n(\cos \alpha) \rho(\mathbf{r}') \, dV',

where P_n denotes the Legendre polynomials and \alpha is the angle between \mathbf{r} and \mathbf{r}'.^[2] For full three-dimensional cases with arbitrary angular dependence, this extends to a sum over spherical harmonics Y_l^m(\theta, \phi), replacing the axial symmetry with a complete angular basis to capture all directional variations.^[5] The successive terms represent distinct multipole orders: the n=0 (monopole) term captures the net source strength, such as total charge; the n=1 (dipole) term accounts for the first-order asymmetry or moment of the distribution; the n=2 (quadrupole) and higher terms describe more refined spatial variations, serving as perturbative corrections that diminish rapidly for large r.^[2] This hierarchical structure enables truncation at low orders for practical computations while retaining accuracy far from the source. The multipole expansion originated with George Green's 1828 essay on the mathematical analysis of electricity and magnetism, where it was introduced for electrostatic potentials.^[6]

Convergence Criteria

The multipole expansion of the electrostatic potential converges absolutely when the distance r from the expansion center to the observation point exceeds the maximum distance r'_{\max} of any source charge from the center, ensuring the observation point lies outside the smallest sphere enclosing all sources. This condition, r > r'_{\max}, guarantees the validity of the underlying Taylor series expansion of the reciprocal distance $1/|\mathbf{r} - \mathbf{r}'| in powers of r'/r. For a localized charge distribution confined within a sphere of radius R, the expansion holds for all r > R.^[7]^[2] The series converges, exhibiting rapid convergence as r increases, with successive terms diminishing in magnitude for large observation distances. However, inside the source region (r < r'_{\max}), the expansion generally diverges due to the proximity of singularities from the source charges. To quantify truncation errors, the remainder after the term of order \ell follows from the Taylor theorem applied to the generating function, bounding the error by the magnitude of the next term. For practical estimates, if the expansion is truncated after the $2^\ell-pole (order \ell), the relative error scales as O\left( (R/r)^{\ell+1} \right), where R is the source radius, highlighting the expansion's utility for far-field approximations.^[7]^[2]^[8] In two dimensions, the multipole expansion for the logarithmic electrostatic potential corresponds to a Laurent series in the complex plane, with convergence governed by the annulus between the expansion center and the nearest source singularity, analyzable via analytic continuation principles from complex analysis. This formulation underscores the series' geometric convergence radius and its limitations near sources.^[7]^[9]

Electrostatic Potential Expansion

Cartesian Form

The electrostatic potential \phi(\mathbf{r}) due to a localized charge distribution \rho(\mathbf{r}') at a point \mathbf{r} far from the origin is given by Coulomb's law in integral form:

\phi(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \int \frac{\rho(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|} \, dV'.

This expression holds for any charge distribution confined within a region much smaller than r = |\mathbf{r}|.^[10] To obtain the multipole expansion in Cartesian coordinates, expand the denominator $1/|\mathbf{r} - \mathbf{r}'| using a multivariate Taylor series around \mathbf{r}' = 0, treating \mathbf{r} as fixed and \mathbf{r}' as the small displacement. The general term in the expansion is

\frac{1}{|\mathbf{r} - \mathbf{r}'|} = \sum_{n=0}^{\infty} \frac{(-1)^n}{n!} r'_i r'_j \cdots r'_k \frac{\partial^n}{\partial r_i \partial r_j \cdots \partial r_k} \left( \frac{1}{r} \right),

where summation over repeated indices i, j, \ldots from 1 to 3 (corresponding to x, y, z) is implied, and the n-th order term involves n factors of r'_l and n partial derivatives. This series converges for r' < r. Substituting into the potential integral yields the multipole series by interchanging the sum and integral, assuming the charge distribution is localized.^[11] The lowest-order term (n=0) is the monopole contribution, where Q = \int \rho(\mathbf{r}') \, dV' is the total charge. The first-order (n=1) term gives the dipole, with components p_i = \int r'_i \rho(\mathbf{r}') \, dV'. For the second-order (n=2) quadrupole term, the Cartesian tensor is defined in traceless form as

Q_{ij} = \int (3 r'_i r'_j - \delta_{ij} r'^2) \rho(\mathbf{r}') \, dV',

which is symmetric (Q_{ij} = Q_{ji}) and traceless (Q_{ii} = 0), capturing the deviation from spherical symmetry with five independent components. Higher-order terms follow similarly but grow more complex.^[12] The resulting multipole expansion of the potential, truncated at quadrupole order, is

\phi(\mathbf{r}) \approx \frac{1}{4\pi\epsilon_0} \left[ \frac{Q}{r} + \frac{\mathbf{p} \cdot \mathbf{r}}{r^3} + \frac{1}{6} Q_{ij} \frac{\partial^2}{\partial r_i \partial r_j} \left( \frac{1}{r} \right) + \cdots \right],

where the derivatives act on $1/r, and \partial_i \partial_j (1/r) = (3 r_i r_j - \delta_{ij} r^2)/r^5. Each term falls off as $1/r^{n+1} for the n-th multipole, allowing truncation based on distance and charge asymmetry.^[10] This Cartesian form is intuitive for systems with rectangular or cubic symmetry, as the tensor components align directly with coordinate axes, facilitating computation for anisotropic distributions. However, it becomes cumbersome for higher multipoles due to the increasing number of tensor components (e.g., 15 for octupole before symmetrization), often requiring irreducible representations for efficiency.^[10]

Spherical Harmonic Form

The spherical harmonic form of the multipole expansion leverages the rotational invariance of the electrostatic potential by separating the radial and angular dependencies using spherical coordinates. This approach is particularly suited for systems with spherical symmetry, providing a compact representation of the angular variations through the orthogonal basis of spherical harmonics Y_{lm}(\theta, \phi). The expansion begins with the separation of variables for the reciprocal distance in the potential integral, which forms the foundation for expressing the potential due to a localized charge distribution. The key identity for the expansion of $1/|\mathbf{r} - \mathbf{r}'| assumes r > r' (i.e., the observation point is outside the charge distribution) and is given by

\frac{1}{|\mathbf{r} - \mathbf{r}'|} = \sum_{l=0}^{\infty} \sum_{m=-l}^{l} \frac{4\pi}{2l+1} \frac{r_<^l}{r_>^{l+1}} Y_{lm}^*(\theta', \phi') Y_{lm}(\theta, \phi),

where r_< = \min(r, r'), r_> = \max(r, r'), and the spherical harmonics satisfy the orthogonality relation \int Y_{lm}^*(\theta', \phi') Y_{l'm'}(\theta', \phi') \, d\Omega' = \delta_{ll'} \delta_{mm'}.^[13] This expansion exploits the addition theorem for spherical harmonics, reducing the angular dependence to Legendre polynomials when m=0, but generalizing to full azimuthal dependence for arbitrary m. For the exterior electrostatic potential \phi(\mathbf{r}) generated by a charge density \rho(\mathbf{r}') confined within a region of radius a < r, the multipole series becomes

\phi(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \sum_{l=0}^{\infty} \sum_{m=-l}^{l} \frac{4\pi}{2l+1} \frac{q_{lm}}{r^{l+1}} Y_{lm}(\theta, \phi),

where the multipole moments are defined as

q_{lm} = \int r'^l Y_{lm}^*(\theta', \phi') \rho(\mathbf{r}') \, dV'.

In this convention, the factor \frac{4\pi}{2l+1} is included explicitly for direct correspondence with the harmonics' normalization.^[13] The l=0 term represents the monopole (total charge), l=1 the dipole, and higher l higher-order multipoles, with each l contributing $2l+1 independent moments due to the azimuthal quantum number m. This spherical form relates to the Cartesian tensor expansion through the theory of irreducible representations of the rotation group SO(3). The multipole moments q_{lm} transform as irreducible spherical tensors of rank l, corresponding to the $2^l-pole in the traditional nomenclature (e.g., l=2 for quadrupole). Unlike the Cartesian form, which uses traceless symmetric tensors of rank l with (2l+1) independent components after removing lower-order traces, the spherical harmonics provide a manifestly irreducible basis, facilitating rotationally invariant computations. The orthogonality of spherical harmonics underpins the computational efficiency of this expansion, enabling fast evaluation via series truncation and analytic translations between multipole and local expansions. This property forms the basis for fast multipole methods (FMM), which achieve O(N) complexity for N-body interactions in three dimensions by hierarchically grouping charges and using spherical harmonic rotations for far-field approximations.

Multipole Moments

Charge Moments

The monopole moment of a charge distribution \rho(\mathbf{r}) is defined as the total charge

Q = \int \rho(\mathbf{r}) \, dV,

which is independent of the origin chosen for the coordinates.^[14] This scalar quantity represents the net charge of the distribution and, when nonzero, provides the leading term in the far-field electrostatic potential, equivalent to that of a point charge at the origin.^[14] The dipole moment is a vector defined by

\mathbf{p} = \int \mathbf{r} \, \rho(\mathbf{r}) \, dV,

which quantifies the effective separation of positive and negative charges in the distribution.^[14] Physically, it arises from the first moment of the charge density and becomes the dominant contribution to the potential when the monopole vanishes. However, \mathbf{p} depends on the choice of origin unless Q = 0, as shifting the origin by \mathbf{a} transforms it to \mathbf{p}' = \mathbf{p} + Q \mathbf{a}.^[15] Higher-order multipole moments describe more detailed aspects of the charge distribution's asymmetry and are represented as symmetric traceless tensors of increasing rank. The quadrupole moment, for instance, is the second-order tensor

Q_{ij} = \int (3 x_i x_j - \delta_{ij} r^2) \, \rho(\mathbf{r}) \, dV,

which captures deviations from spherical symmetry in the charge arrangement.^[1] Similarly, the octupole and subsequent moments follow analogous tensor forms for third and higher orders. These moments generally transform under translation of the origin; for example, the dipole moment vanishes at the center of charge \mathbf{r}_c = \frac{1}{Q} \int \mathbf{r} \, \rho(\mathbf{r}) \, dV when Q \neq 0, and higher moments obey relations derived from Steiner's theorem, akin to the parallel axis theorem for inertia tensors.^[16] In terms of units and scaling, the l-th multipole moment has dimensions of charge times length to the power l, reflecting the integral's dependence on the spatial extent of the distribution: for a system of characteristic size L, the moment scales as Q L^l./Physical_Properties_of_Matter/Atomic_and_Molecular_Properties/Intermolecular_Forces/Multipole_Expansion) This scaling ensures that higher moments diminish in influence at large distances, justifying their use in asymptotic expansions of the potential.

Molecular Conventions

In quantum chemistry, multipole moments for molecules are derived from the total charge density \rho(\mathbf{r}), which combines nuclear and electronic contributions: \rho(\mathbf{r}) = \sum_A Z_A \delta(\mathbf{r} - \mathbf{R}_A) - \sum_i |\psi_i(\mathbf{r})|^2, where Z_A and \mathbf{R}_A are the charge and position of nucleus A, and the electronic density is obtained from the squared modulus of the molecular wavefunction \psi_i for occupied orbitals (with the negative sign accounting for electron charge).^[17] This density is integrated to compute moments, with the choice of origin typically at the center of mass for neutral molecules, where the dipole moment is origin-independent (since Q=0), though higher moments generally remain origin-dependent unless all lower-order moments vanish at that point, or at a specific nuclear coordinate for ionic or asymmetric systems.^[18] Conventions for the sign and direction of multipole moments vary between traditions, notably the Buckingham convention prevalent in molecular physics and the IUPAC standard in chemical spectroscopy. Under the Buckingham convention, the dipole moment is defined as \boldsymbol{\mu}_\alpha = \sum_A Z_A R_{A\alpha} - \int r_\alpha \rho_\mathrm{el}(\mathbf{r}) \, d\mathbf{r}, where the nuclear term is positive and the electronic integral carries a negative sign, yielding a vector pointing from negative to positive charge regions.^[17] In contrast, the IUPAC convention aligns the dipole vector direction from the negative to the positive charge center, consistent with spectroscopic measurements, though the arrow symbol in chemical diagrams often points oppositely for illustrative purposes.^[19] For higher-order moments, molecular physics employs traceless Cartesian tensors to eliminate origin dependence and simplify anisotropy descriptions. The quadrupole moment tensor Q_{\alpha\beta} is rendered traceless such that Q_{xx} + Q_{yy} + Q_{zz} = 0, with the anisotropy often quantified as \Delta Q = Q_{zz} - (Q_{xx} + Q_{yy})/2, particularly useful for linear molecules where off-diagonal elements vanish.^[20] Computing these moments via ab initio methods introduces challenges due to basis set dependence, where incomplete basis sets lead to systematic errors in the electronic density and thus in dipole and quadrupole values; convergence requires augmented, diffuse functions for accurate long-range properties.^[21] Experimentally, moments are determined through spectroscopic techniques like Stark effect in microwave spectroscopy for dipoles or collision-induced absorption for quadrupoles, providing benchmarks against theoretical predictions.^[22] These conventions were standardized in the 1970s through foundational work on spherical tensor formulations and molecular electrostatics, notably by Buckingham, enabling consistent comparisons across quantum chemical calculations and experiments.^[23]

Interactions Between Distributions

Interaction Energy

The electrostatic interaction energy U between two non-overlapping charge distributions \rho_A(\mathbf{r}) and \rho_B(\mathbf{r}') is given by the double integral

U = \frac{1}{4\pi\epsilon_0} \iint \frac{\rho_A(\mathbf{r}) \rho_B(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|} \, dV \, dV',

where \epsilon_0 is the vacuum permittivity. This expression arises from integrating Coulomb's law over the continuous charge densities of the distributions. For well-separated distributions, where the distance R between their centers greatly exceeds the sizes of the distributions, the denominator |\mathbf{r} - \mathbf{r}'| can be expanded using the multipole series. The resulting interaction energy approximates to

U \approx \sum_{l_a=0}^\infty \sum_{l_b=0}^\infty \frac{1}{R^{l_a + l_b + 1}} T_{l_a l_b} q_{l_a}^A q_{l_b}^B,

where q_{l_a}^A and q_{l_b}^B are the multipole moments of orders l_a and l_b for distributions A and B, respectively, and T_{l_a l_b} is the interaction tensor that encodes the angular dependence and incorporates the factor $1/(4\pi\epsilon_0). This double sum captures contributions from all pairwise interactions between multipoles of the two distributions, converging rapidly for large R. The leading term in this expansion is the monopole-monopole interaction, U_{\text{mono-mono}} = \frac{Q_A Q_B}{4\pi\epsilon_0 R}, where Q_A and Q_B are the total charges of the distributions; this dominates for charged systems at long range. The next term, the monopole-dipole interaction, is U_{\text{mono-dip}} = \frac{Q_A}{4\pi\epsilon_0} \mathbf{p}_B \cdot \nabla \left( \frac{1}{R} \right), where \mathbf{p}_B is the dipole moment of distribution B; it arises when one distribution is charged and the other has a permanent dipole. The dipole-dipole term follows as U_{\text{dip-dip}} = \frac{1}{4\pi\epsilon_0} \mathbf{p}_A \cdot \nabla \nabla \left( \frac{1}{R} \right) \cdot \mathbf{p}_B, which is crucial for neutral polar molecules and scales as $1/R^3. Higher-order terms, such as quadrupole-quadrupole, contribute at shorter ranges or for systems lacking lower moments. The multipole expansion exhibits symmetry properties related to the parity of the interacting moments: interactions between multipoles of even total order (e.g., monopole-monopole or dipole-dipole) are even under inversion and dominate electrostatic repulsion or attraction between like-charged distributions, while odd total order terms (e.g., monopole-dipole) are odd and describe induction effects in perturbation theory. Van der Waals dispersion forces, arising from second-order perturbation theory involving correlated fluctuations, correspond to even-order terms beyond the first-order electrostatic expansion. A general framework for this interaction expansion in spherical tensors was developed by Stone in 1978, providing compact expressions for the electrostatic energy using irreducible spherical tensor representations of the multipole moments and interaction operators, which facilitate efficient computation of orientation-dependent forces and torques between molecules.^[24]

Expansion of Potential

The multipole expansion of the electrostatic potential due to a localized charge distribution A, evaluated at points within the region of another localized charge distribution B, provides an approximate representation of the potential field when the separation between the distributions greatly exceeds their individual extents. This expansion is particularly useful for distant systems, where the potential varies slowly across B relative to the separation vector \mathbf{R} from the origin of A to the origin of B. Consider a point in the region of B at position \mathbf{x} = \mathbf{R} + \mathbf{s}, where \mathbf{s} is the position vector relative to B's origin and |\mathbf{s}| \ll |\mathbf{R}|. The potential \phi_A(\mathbf{x}) generated by distribution A can then be expressed via the multivariable Taylor series expansion around \mathbf{R}:

\phi_A(\mathbf{R} + \mathbf{s}) = \sum_{n=0}^{\infty} \frac{1}{n!} \left( \mathbf{s} \cdot \nabla \right)^n \phi_A(\mathbf{R}),

with the gradient operator \nabla acting on the coordinates of the evaluation point and evaluated at \mathbf{R}. Here, \phi_A(\mathbf{R}) itself admits a multipole expansion in Cartesian form, expressing the potential at B's origin as a series in inverse powers of R:

\phi_A(\mathbf{R}) = \frac{1}{4\pi \epsilon_0} \left[ \frac{Q_A}{R} + \mathbf{p}_A \cdot \nabla \left( \frac{1}{R} \right) + \frac{1}{2} Q_{ij}^A \partial_i \partial_j \left( \frac{1}{R} \right) + \cdots \right],

where Q_A is the monopole moment (total charge), \mathbf{p}_A the dipole moment, and Q_{ij}^A the traceless quadrupole tensor of A, with higher-order terms following analogously. Substituting this multipole series into the Taylor expansion yields \phi_A(\mathbf{R} + \mathbf{s}) as a double expansion: the leading multipoles of A contracted with tensorial derivatives of $1/R at the origin of B, further Taylor-expanded in powers of \mathbf{s} to approximate the potential across B. This structure allows truncation at low orders for practical computations when higher multipoles contribute negligibly.^[25] For scenarios involving close or partially overlapping but non-penetrating distributions—such as atomic or molecular systems where standard spherical expansions diverge—a multipole expansion in bipolar coordinates proves effective. Bipolar coordinates (\tau, \sigma, \phi) are defined with foci at the origins of A and B, separated by distance d < R, transforming the interaction via elliptic coordinates along the axis joining the centers; the potential is then expanded as a series of bipolar harmonics, ensuring convergence in the interpenetrating region without assuming large separation. This potential expansion underpins perturbation theory in quantum mechanics, where the interaction between distant subsystems is treated as a perturbative correction to unperturbed Hamiltonians, with multipole terms providing the leading-order long-range couplings.

Applications and Examples

Physical Systems

In electrostatics, the multipole expansion provides a powerful framework for approximating the potential due to localized charge distributions at large distances, enabling the analysis of interactions between systems such as ions and atoms. For instance, the interaction energy between a point-charge ion and a neutral atom can be expanded in terms of the atom's multipole moments, starting with the monopole term for the ion's charge and including the atom's induced dipole and higher-order responses for short-range corrections.^[26] This approach is particularly useful in classical models of scattering and van der Waals forces, where the leading terms capture the dominant physics while higher multipoles account for angular dependencies and asymmetries in the charge distribution.^[27] The Earth's global electric field, arising from atmospheric charge separations, can also be approximated via multipole expansions to model regional variations, though the dominant fair-weather field is often treated as a near-uniform downward-directed component with small dipole perturbations from ionospheric sources.^[28] In gravitational physics, the Newtonian potential for planetary fields is expanded in spherical harmonics to account for deviations from spherical symmetry, with the leading correction being the quadrupole term characterized by the zonal harmonic coefficient J_2. For Earth, this term arises primarily from rotational oblateness, flattening the planet at the poles and bulging at the equator, yielding J_2 \approx 1.0827 \times 10^{-3}, which influences satellite orbits and tidal dynamics.^[29] The full expansion takes the form

\Phi(r, \theta) = -\frac{GM}{r} \left[ 1 - \sum_{l=2}^{\infty} J_l \left( \frac{R}{r} \right)^l P_l(\cos \theta) \right],

where M is the planet's mass, R its reference radius, and P_l are Legendre polynomials; higher terms like J_4 are much smaller, on the order of $10^{-5}.^[30] In magnetostatics, multipole expansions describe the field from current distributions, with the magnetic dipole moment defined as

\mathbf{m} = \frac{1}{2} \int \mathbf{r} \times \mathbf{J}(\mathbf{r}) \, dV,

where \mathbf{J} is the current density; this moment quantifies the effective loop area and strength for localized sources like atomic currents or planetary dynamos.^[31] Far from the source, the magnetic field approximates the dipole form

\mathbf{B}(\mathbf{r}) \approx \frac{\mu_0}{4\pi} \frac{3(\mathbf{m} \cdot \hat{\mathbf{r}}) \hat{\mathbf{r}} - \mathbf{m}}{r^3},

which captures the field's $1/r^3 decay and angular variation, essential for modeling Earth's geomagnetic field dominated by its internal dipole. A classic example is the electrostatic potential of a uniformly charged sphere of total charge Q and radius R. Outside the sphere (r > R), the potential is purely monopolar, equivalent to a point charge at the center:

\phi(r) = \frac{1}{4\pi \epsilon_0} \frac{Q}{r},

with all higher multipole moments vanishing due to spherical symmetry. Inside the sphere (r < R), the potential includes higher-order terms when expanded around the center, reflecting the quadratic variation:

\phi(r) = \frac{1}{4\pi \epsilon_0} \frac{Q}{2R} \left( 3 - \left( \frac{r}{R} \right)^2 \right),

which can be expressed as a series involving even Legendre polynomials starting from the monopole.^[10] Experimental validation of gravitational multipole expansions has been advanced since the 1960s through satellite measurements, with the GRACE mission (2002–2017) and its successor GRACE Follow-On (GRACE-FO, launched 2018 and operational as of 2025) providing high-resolution maps of Earth's gravity field variations, enabling precise determination of low-degree multipoles like J_2 and their temporal changes due to mass redistributions such as ice melt and groundwater depletion.^[32]^[33] These data confirm the expansion's accuracy to scales of hundreds of kilometers, with J_2 variations tracked at the $10^{-11} level relative to its mean value.^[34]

Chemical and Biological Contexts

In quantum chemistry, multipole expansions play a crucial role in computing intermolecular forces through methods like symmetry-adapted perturbation theory (SAPT), where electrostatic interactions are derived from multipole moments obtained via Hartree-Fock (HF) or density functional theory (DFT) wavefunctions.^[35] These moments enable the decomposition of interaction energies into physically interpretable components, such as electrostatics, induction, and dispersion, providing insights into molecular associations at the quantum level.^[36] For instance, dipole-dipole terms within the multipole expansion dominate the electrostatic contribution to hydrogen bonding, accurately capturing the attraction between electronegative atoms and hydrogen in systems like water dimers.^[37] In drug design, multipole expansions enhance the modeling of protein-ligand binding energies by representing electrostatic interactions more precisely than simple point charges. The Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field employs permanent atomic multipoles up to quadrupole order, combined with inducible dipoles, to simulate polarization effects in biomolecular complexes.^[38] This approach has demonstrated superior accuracy in predicting binding affinities for protein-ligand systems, outperforming fixed-charge models in scenarios involving charged or polar ligands.^[39] Biological applications leverage multipole expansions to elucidate electrostatic mechanisms in key processes. In ion channels, such as potassium channels, multipole moments derived from partial charge distributions inform selectivity by quantifying ion-protein electrostatic interactions at the selectivity filter, reproducing observed preferences for K⁺ over Na⁺.^[40] Similarly, for DNA base pairing, distributed multipole analyses predict stable Watson-Crick geometries through accurate electrostatic potentials between bases, highlighting the role of higher-order moments in stabilizing hydrogen-bonded pairs.^[41] Advanced methods like distributed multipole analysis (DMA), developed by Stone and Alderton in the 1980s, distribute multipoles across atomic sites to overcome limitations of centered expansions, yielding more accurate intermolecular potentials for molecular simulations.^[42] However, standard HF or DFT computations of these moments often underestimate electron correlation effects, which influence higher-order multipoles and require post-HF methods, such as coupled-cluster theory, for improved accuracy in correlated systems.^[43]

General Properties

Uniqueness and Completeness

The uniqueness of the multipole expansion for the electrostatic potential outside a charge distribution follows from the uniqueness theorem for solutions to Laplace's equation in the exterior domain. For a localized charge distribution with compact support, the potential satisfies ∇²Φ = 0 for r greater than the radius enclosing all charges, and specifying the behavior at infinity (Φ → 0 as r → ∞) uniquely determines the solution among all harmonic functions vanishing at infinity.^[44] This uniqueness extends to the multipole representation, where the coefficients (multipole moments) are uniquely fixed by the boundary values on a sphere surrounding the sources. The Kelvin inversion theorem maps interior solutions to exterior ones while preserving harmonicity.^[45] The completeness of the multipole expansion arises because the spherical harmonics {Y_{ℓm}(θ, φ)} form a complete orthonormal basis for the Hilbert space L² of square-integrable functions on the unit sphere. Any harmonic function in the exterior region can thus be uniquely expanded as a series in spherical harmonics times radial factors (1/r^{ℓ+1}), ensuring that the multipole series fully represents the potential without omission of components. The multipole expansion provides a far-field representation of the potential, valid for observation points distant from the sources, while its dual—the local expansion—offers a near-field description centered away from the sources, using positive powers of r. These dual representations are interconvertible via translation operators in the fast multipole method, allowing efficient evaluation in both regimes. Multipole moments depend on the choice of expansion origin, introducing a gauge-like ambiguity under translations, but physical observables such as interaction energies remain invariant due to compensatory changes in higher-order moments.^[44] For a charge density ρ with compact support within a sphere of radius a, the infinite multipole series exactly equals the true potential for all r > a, converging term-by-term to the solution of Poisson's equation outside the sources.^[44]

Relation to Other Expansions

The multipole expansion bears a close mathematical resemblance to Laurent series expansions, particularly in the context of far-field representations. In two dimensions, the potential generated by a localized charge distribution admits a Laurent series expansion in powers of $1/z, where z is the complex position coordinate; the terms with negative powers encapsulate the multipolar contributions, providing an asymptotic description for large |z|. In three dimensions, this concept generalizes to an expansion in inverse powers of r modulated by spherical harmonics, serving as the spherical analog of the 2D Laurent series for capturing angular dependencies in the far field. In contrast to Fourier series or transforms, which decompose fields from periodic or extended sources into plane-wave components, multipole expansions are specifically formulated for compact, localized distributions, enabling hierarchical approximations that decay rapidly with distance. This distinction is leveraged in computational methods, such as the fast multipole method (FMM) developed by Greengard and Rokhlin, which accelerates the summation of interactions from N sources and targets from O(N^2) to O(N) complexity by translating and combining multipole expansions across octree levels. Multipole expansions extend beyond scalar electrostatic potentials to vector and tensor fields in broader physical theories. In electromagnetism, the vector potential \mathbf{A} is expanded in a multipolar form, often involving the curl \nabla \times \mathbf{A} to separate transverse components into electric and magnetic multipoles of order l, conventionally termed $2^l-poles (e.g., dipole for l=1, quadrupole for l=2). In general relativity, similar expansions apply to the metric tensor using tensor spherical harmonics, decomposing the gravitational field of isolated sources into mass (electric-type) and current (magnetic-type) multipoles that characterize spacetime curvature at large distances. In quantum mechanics, multipole expansions describe radiative transitions in atomic and molecular spectra, where the interaction Hamiltonian couples to electromagnetic fields via multipolar operators; selection rules, derived from angular momentum conservation and parity considerations, allow transitions where the change in total angular momentum satisfies |ΔJ| ≤ l ≤ J_i + J_f for 2^l-pole radiation (with J_i + J_f ≥ 1), with parity change for electric multipoles if l is odd (no change if even) and the opposite for magnetic multipoles.^[46] A fundamental distinction from plane-wave expansions lies in the asymptotic behavior: multipole terms decay as $1/r^{l+1} in the far field, allowing truncation for distant observations, whereas plane-wave representations maintain oscillatory propagation without inherent decay, suiting scattering or uniform wave problems.

References

[1]
[PDF] Today in Physics 217: multipole expansion
Oct 21, 2002 · Dipole moment is defined the same way in cgs and MKS. Expressions for potential and field still need a factor of to convert from cgs to MKS. -q.
[2]
[PDF] MULTIPOLE EXPANSION - UT Physics
Hence, let's define the spherical harmonics of multipoles according to. Mℓ,m def. = r. 4π. 2ℓ + 1 ZZZ. rℓ × Y. ∗. ℓ,m(θ, φ) × ρ(r, θ, φ) d. 3. Vol(r, θ, φ). (68).
[3]
[PDF] Multipole Expansion of the Electrostatic Potential—CE
define the components of a 3 × 3 matrix or second-rank tensor called the quadrupole moment t. Q. This quantity behaves like the moment of inertia. I r rr dm.
[4]
[PDF] Outlines - Physics - University of Central Florida
Oct 13, 2016 · We will also explain (c) what is multipole expansion and why we are interested in multipole expansion. Electric dipole. A charge q is located ...
[5]
Multipole Expansion (and Legendre Polynomials) | PHY309
They arise naturally when you separate variables in spherical coordinates · They arise naturally when you use Coulomb's equation for potential, and consider it ...
[6]
[PDF] Green 1828 - The Application of Mathematical Analysis
MATHEMATICAL ANALYSIS TO THE THEORIES OF ELECTRICITY. AND MAGNETISM. INTRODUCTORY OBSERVATIONS. A HE object of this Essay is to submit to Mathematical Analysis ...
[7]
The definition of multipole moments for extended bodies
Multipole Expansion: Unifying Formalism for Earth and Planetary Gravitational Dynamics ... Henri Poincaré,11, 221–237. MathSciNet Google Scholar. Beiglböck ...<|control11|><|separator|>
[8]
[PDF] Math of Multipole Expansion - UT Physics
any r<R. If we analytically continue it to the complex r, it would converge ... the angular dependence of the multipole expansion using the following.
[9]
[PDF] Classical Electrodynamics - Duke Physics
antenna is on the order of λ, the multipole expansion is an extremely accurate and rapidly converging approximation. That is, after all, why we use it so ...
[10]
[PDF] Exponential Convergence for Multipole and Local Expansions and ...
This paper derives multipole and local expansions for 2D layered media, proves their exponential convergence, and shows convergence depends on polarization ...Missing: criteria | Show results with:criteria
[11]
Multipole Expansion - Richard Fitzpatrick
Multipole expansion is a type of expansion for the electrostatic potential of a charge distribution, using multipole moments. The first three are monopole, ...
[12]
[PDF] B. Svetitsky, October 2021 TAYLOR EXPANSION OF THE ...
The Cartesian multipole expansion, valid at large distances r, can be derived by expanding. 1/|r − r/| in powers of r/. Treating r as a constant, ...
[13]
[PDF] Multipole expansion of Electrostatic Potential - Squarespace
Multipole expansion in Cartesian Coordinates. One can expand the potential expression (1) in Cartesian coordinates. In this method, we expand the function. 1. | ...<|control11|><|separator|>
[14]
[PDF] Lecture 8 Notes, Electromagnetic Theory I
Multipole expansion uses spherical harmonics to approximate potential far from a localized charge. The l=0 term is the monopole, l=1 is dipole, etc.
[15]
26. Electric Multipoles - Galileo and Einstein
In fact, charge distributions are often given in Cartesian coordinates, so it's useful to express these (θ,ϕ) multipole moments in Cartesian coordinates:.
[16]
[PDF] Lecture Notes 8 - High Energy Physics
Linear superposition of multipole electric fields!!! Thus, we see that, for observation / field point distances far away from the (arbitrary) localized electric ...
[17]
[PDF] arXiv:1105.3117v1 [gr-qc] 16 May 2011
May 16, 2011 · Then, by means of parallel axis theorem (Huygens-Steiner theorem), they apply general expressions for the quadrupole moment and their time ...
[18]
[PDF] Molecular Electric Moments Calculated By Using Natural Orbital ...
Theoretical calculations are therefore essential but challenging for quantum chemistry methods. ... convention of Buckingham [3, 4] is expressed by the formula.
[19]
Ab initio calculations on the multipole moments and dipole ...
There have been only three known experimental studies of the spectroscopy of the PO+ molecular ion. Dressler [3] measured eight bands in the A–X emission ...
[20]
[PDF] "Conventions, Symbols, Quantities, Units and Constants for High ...
This document summarizes conventions, symbols, quantities, units, and constants for high-resolution molecular spectroscopy, including quantity calculus, SI ...
[21]
Multipole-moment effects in ion–molecule reactions at low ... - NIH
The electric quadrupole moment of N2 can be described by a traceless, diagonal rank-2 tensor with elements Qzz = −2Qxx = −2Qyy, and a negative value of Qzz (= − ...
[22]
Performance and basis set dependence of density functional theory ...
A study is presented concerning the performance and basis set dependence of density functional methods in the calculation of dipole and quadrupole moments.
[23]
An accurate ab initio study of the multipole moments and ...
The accuracy of determining multipole moments and polarizabilities from collision induced spectra is discussed.
[24]
Spherical tensor theory of molecular multipole moments and ...
A spherical tensor theory of molecular multiple moments and molecular polarizabilities is developed. It is shown that the spherical tensor form has several ...
[25]
[PDF] Multipole Expansion of the Electrostatic Potential - UT Physics
Multipole expansion decomposes the electrostatic potential into a series of multipole potentials, such as monopole, dipole, and quadrupole, based on charge ...Missing: criteria | Show results with:criteria
[26]
[PDF] ATOMS - NASA Technical Reports Server (NTRS)
The object of this calculation is to examine the contribution to the interaction energy between two neutral hydrogenic atoms in their ground state because ...
[27]
[PDF] Multipoles; Macroscopic Media; Dielectrics
Dec 23, 2000 · The monopole moment is just the total charge of the distribution. Thus the monopole term gives the potential due to the charge as a whole. Since ...
[28]
[PDF] The Multipole Structure of Earth's STEP Signal
Mar 25, 2004 · The rotational flattening of Earth and geographical irregularities of our planet's crust then produces a multipole structure for the Equivalence.
[29]
[PDF] 3.02 Potential Theory and Static Gravity Field of the Earth - Elsevier
yields a multipole expansion (so called from electro- statics) of the potential. The spherical harmonic. (Stokes) coefficients are multipoles of the density.
[30]
[PDF] gravity.pdf - atmo.arizona.edu
The Earth gravitational acceleration we feel includes the effects of the Earth's spin. ... Plugging in values of J2 = 1.08270x10-3 and a3 ω. 2/GM = 3.46775x10-3 ...
[31]
[PDF] 3. Magnetostatics - DAMTP
The first equation of magnetostatics, r ⇥ B = µ0J. (3.3) is known as ... where m is the magnetic dipole moment of the current distribution. Suddenly ...
[32]
[PDF] GGM02 – An Improved Earth Gravity Field Model from GRACE
GGM02 is still relatively short, the J2 harmonic was constrained to its long-term (multi-decadal) mean value from EGM96 (using the original EGM96 sigma for J2),.
[33]
The gravity recovery and climate experiment: Mission overview and ...
May 8, 2004 · The GRACE mission is designed to track changes in the Earth's gravity field for a period of five years. Launched in March 2002, ...Missing: multipoles | Show results with:multipoles
[34]
SAPT codes for calculations of intermolecular interaction energies
May 13, 2020 · Symmetry-adapted perturbation theory (SAPT) is a method for calculations of intermolecular (noncovalent) interaction energies.
[35]
A Quantum Chemical Topology Picture of Intermolecular ...
Jul 19, 2021 · Based on the Interacting Quantum Atoms approach, we present herein a conceptual and theoretical framework of short-range electrostatic interactions.
[36]
Multipole Expansion of Atomic Electron Density Fluctuation ... - NIH
Nov 14, 2024 · In this study, we present a multipole-extended second-order DFTB (mDFTB2) method that takes into account atomic dipole and quadrupole interactions.
[37]
Current Status of the AMOEBA Polarizable Force Field
AMOEBA is shown to be especially successful on protein−ligand binding and computational X-ray crystallography where polarization and accurate electrostatics are ...The AMOEBA Force Field · Protein−Ligand Binding · Software Infrastructure for...
[38]
The Polarizable Atomic Multipole-based AMOEBA Force Field for ...
AMOEBA utilizes an interactive atomic dipole induction scheme where the field produced by permanent multipoles and induced dipoles induces a dipole at each ...
[39]
Incidence of partial charges on ion selectivity in potassium channels
Jan 23, 2006 · The charges derived from the potential generally reproduce the first multipole moments well and optimally also intermolecular interactions.
[40]
What base pairings can occur in DNA? A distributed multipole study ...
Ab initio distributed multipole electrostatic calculations are used to predict likely nucleic acid base pair structures for both the gas phase and within a ...Missing: expansion | Show results with:expansion
[41]
Distributed multipole analysis: Methods and applications
Abstract. The conventional multipole expansion gives a description of electrostatic interactions which is only useful at long distances.
[42]
Accurate electric multipole moment, static polarizability and ...
Electron correlation effects on D1⌽ and D2⌽ are less im- portant than those observed for the quadrupole moment. It is also surprising that for such a high ...<|control11|><|separator|>
[43]
Classical Electrodynamics, 3rd Edition
**Summary of Multipole Expansion Uniqueness, Completeness, Origin Dependence from Jackson's Book Description:**
[44]
(PDF) Kelvin transformation and inverse multipoles in electrostatics
Jan 6, 2025 · The inversion in the sphere or Kelvin transformation, which exchanges the radial coordinate for its inverse, is used as a guide to relate ...Missing: uniqueness | Show results with:uniqueness