Multipole Expansion of Radiation Fields

• Multipole expansion is a systematic framework that decomposes radiation fields into hierarchically ordered angular modes and source moments.
• It employs spherical harmonics and irreducible STF tensors to capture angular complexity, decay properties, and far-field behavior in electromagnetic, acoustic, and gravitational systems.
• This approach underpins conservation laws and gauge invariance, enabling precise extraction of observables such as radiated energy, momentum, and angular momentum.

Multipole expansions provide a systematic framework for representing the spatial structure of radiation fields—acoustic, electromagnetic, or gravitational—in terms of hierarchical source moments. By expanding solutions of the underlying wave equations (Helmholtz, Maxwell, or linearized Einstein equations) in spherical or tensor harmonics, the emitted fields at large distances can be precisely characterized according to their angular complexity and decay properties. This analytic machinery is central in both classical and quantum radiation theory, with deep connections to conservation laws, gauge invariance, selection rules, and the extraction of physical observables such as radiated energy, momentum, and angular momentum.

1. Mathematical Foundations of Multipole Expansion

The multipole expansion of a radiation field decomposes any solution of a homogeneous or inhomogeneous wave equation generated by a spatially localized source into a sum of spatially orthogonal angular modes. In the setting of electromagnetism, the vector potential in the far-zone is expressed as a series of electric and magnetic multipole fields (Alaee et al., 2017), while in linearized gravity, the metric perturbation is expanded in symmetric trace-free (STF) tensors representing mass-type and current-type moments (Zschocke, 2014).

For a time-harmonic source, the expansion of the field—for example, the electromagnetic vector potential—yields: $$\mathbf{A}(\mathbf{r}, \omega) \sim \sum_{l=1}^\infty \sum_{m=-l}^l \left[ a_{lm}^{(E)} \mathbf{N}_{lm}(k, \mathbf{r}) + a_{lm}^{(M)} \mathbf{M}_{lm}(k, \mathbf{r}) \right]$$ with spherical multipole coefficients $a_{lm}^{(E)}, a_{lm}^{(M)}$, and vector spherical harmonics $\mathbf{N}_{lm}, \mathbf{M}_{lm}$ (Riccardi et al., 2022). Each term corresponds to a definite angular momentum component, with the lowest orders (dipole and quadrupole) dominating the long-wavelength regime.

In Cartesian formalism, the expansion utilizes irreducible STF tensors. For scalar wave equations, the STF moments $Q_{L} = \int d^3x\, \rho(\mathbf{x}) x^{L}$ project out the radiative degrees of freedom directly (Ross, 2012). The connection between Cartesian and spherical expansions is algebraically exact in isotropic media but fails in general in anisotropic backgrounds (Boudec et al., 2024).

2. Source Moments and Their Physical Interpretation

The distinct multipole families—electric, magnetic, and in the electromagnetic context, toroidal and higher-order mean-square-radius (MSR) moments—are defined as weighted integrals over the charge and current density distributions. Explicitly, for the electromagnetic field (Alaee et al., 2017):

• Electric multipole moment:

$$Q_{lm}(\omega) = \int d^3r\, \rho(\mathbf{r}, \omega) j_l(kr) Y_{lm}^*(\hat{\mathbf{r}})$$

• Magnetic multipole moment:

$$M_{lm}(\omega) = \frac{1}{c} \int d^3r\, [\mathbf{r} \times \mathbf{J}(\mathbf{r}, \omega)]\, j_l(kr) Y_{lm}^*(\hat{\mathbf{r}})$$

with $j_l$ the spherical Bessel functions, and $Y_{lm}$ the scalar spherical harmonics.

Toroidal and MSR moments arise from subleading terms in the expansion of the Bessel functions and represent nontrivial current configurations with the same angular distribution as the parent electric/magnetic multipoles but distinct scaling in $k$ (wavenumber) (Li et al., 2018, Nemkov et al., 2018). The hierarchy:

• Electric multipole ($\sim (kr)^l$)
• Toroidal multipole ($\sim (kr)^{l+2}$)
• MSR corrections ($\sim (kr)^{l+2n}$, with $n \ge 1$)

Mass and current multipole moments in gravity involve source integrals over the stress-energy tensor, projected via STF tensors (Zschocke, 2014, Ross, 2012).

3. Far-Field Structure and Radiation Observables

In the radiation zone ($r \gg a$, with $a$ the source size), the field components decay as $1/r$ with leading angular and temporal structure set by the multipole order. For electromagnetic radiation, the electric and magnetic fields in the far field are (Alaee et al., 2017, Riccardi et al., 2022): $$\mathbf{E}(\mathbf{r}) \sim \frac{e^{ikr}}{r} \sum_{l,m} \left\{ A_{lm}(k) \mathbf{X}_{lm}(\hat{r}) + B_{lm}(k) \mathbf{Y}_{lm}(\hat{r}) \right\}$$ with explicit coefficients involving the electric, magnetic, and toroidal moments, and their MSR extensions (Nemkov et al., 2018).

The radiated power is quadratic in the highest time-derivative of each multipole moment: $$\frac{dE}{dt} = \sum_{l=1}^\infty \frac{(l+1)}{l\,l!\,(2l+1)!!} \left\langle \left| \frac{d^{l+1}}{dt^{l+1}} Q_L \right|^2 + \frac{l}{l+1} \left| \frac{d^{l+1}}{dt^{l+1}} J_L \right|^2 \right\rangle$$ for electromagnetism (Ross, 2012), and a structurally analogous expression for gravitation, with mass-/current-type tensor moments (Zschocke, 2014).

Radiation reaction and angular distributions follow straightforwardly from these harmonic decompositions. Non-radiating combinations ("anapoles"), constructed by destructive interference between ordinary and toroidal terms, exemplify configurations with vanishing far-zone fields (Li et al., 2018, Nemkov et al., 2018).

4. Rigorous Formulations and Gauge Structure

The physical content of the multipole expansion is ensured by decomposing the source tensors into irreducible STF components. This removes redundancy due to trace parts and aligns the expansion with the SO(3) representation theory: only the STF part of a rank-$l$ tensor contributes to radiation of multipole order $l$ (Zschocke, 2014, Ross, 2012, Riccardi et al., 2022).

Gauge invariance is systematically maintained. In electromagnetism, multipole charges can be interpreted as Noether charges associated to residual gauge symmetries surviving Lorenz gauge fixing, encoded in the conservation equations for multipole moments of the source and the field (Seraj, 2016). This structure generates an infinite hierarchy of constraints (Ward identities) that translate into exact relations among radiated fields (Compère et al., 31 Mar 2025).

In the quantum context, the multipole expansion of transition amplitudes employs irreducible tensors built from the current density, naturally coupling to photon spherical harmonics and enforcing angular-momentum selection rules. Both electric ($E_l$) and magnetic ($M_l$) multipole transitions arise, with well-defined parity and angular dependencies (Agre, 2014, Casini et al., 2024).

5. Applications and Extensions Across Physical Systems

Multipole expansion underlies key results across domains:

• Electromagnetic scattering and nano-optics: Accurate modeling of scattering, absorption, and radiation forces on particles (including in the Mie regime), requires exact multipole formulations, particularly spherical multipoles for isotropic spheres (Riccardi et al., 2022), and necessitates inclusion of toroidal and MSR moments in complex geometries (Li et al., 2018, Nemkov et al., 2018).
• Quantum radiation theory: Transition rates, photoemission angular distributions, and polarization phenomena use the multipole expansion of matrix elements connecting initial and final states (Agre, 2014, Casini et al., 2024).
• Gravitational radiation: The maintenance of gauge invariance and the uniqueness of radiative multipole moments (mass, current, Weyl-type in higher dimensions) is key for gravitational-wave modeling (Zschocke, 2014, Amalberti et al., 2023).
• Acoustics: In the problem of acoustic radiation force and torque on a particle, the entire formalism carries over: the (complex) beam-shape coefficients $a_{nm}$ governing incident field expansion in regular spherical waves, Mie-type scattering coefficients for the response, and force/torque expressed as convergent multipole sums (Gong et al., 2021).
• Anisotropic and metamaterial media: In backgrounds lacking full rotational symmetry, the equivalence between Cartesian and spherical multipole expansions breaks down, and completeness is only guaranteed for the Cartesian tensor expansion (Boudec et al., 2024).
• Quantum Hall systems and angular-momentum transfer: The far-field radiation of edge states in Dirac materials requires access to very high-order multipoles, determined by the electronic coherence and spatial structure of the edge (Gullans et al., 2017).

6. Formalism at the Level of the Action and Higher-Dimensional Generalizations

The effective field theory (EFT) approach formalizes the multipole expansion directly at the level of the action via a derivative expansion for the long-wavelength sector (Ross, 2012, Amalberti et al., 2023). The multipole moments enter as Wilson coefficients, encoding all source-structure dependence, and the action's structure guarantees manifest gauge invariance order by order. The extension to arbitrary spacetime dimensions introduces new multipolar families (e.g., Weyl-type moments in gravity for $d > 3$) and new coupling patterns (Amalberti et al., 2023).

7. Tabular Comparison: Multipole Types and Properties

Multipole Family	Source Structure	Field Angular Pattern	Dominance in Radiation
Electric ($Q_{lm}$)	Charge density, $\rho(\mathbf{r})$	$Y_{lm}$	Lowest allowed ($l$)
Magnetic ($M_{lm}$)	Current loops, $\mathbf{J}(\mathbf{r})$	VSH (axial parity)	Next to electric
Toroidal ($T_{lm}$)	Complex currents, MSRs	Same as electric, higher $k$-scaling	Subleading ($(kr)^{l+2}$)
MSR ($X_{lm}^{(n)}$)	Radial extension of above	Same as parent	Higher-order corrections
Mass ($I_L$), Current ($J_L$), Weyl ($K_L$)	Stress-energy tensor in gravity	STF tensors, VSH	Quadrupole and above (gravity)

This table summarizes the main multipole families, their physical source, and their scaling or angular dependence in radiating fields.

• (Alaee et al., 2017) "An electromagnetic multipole expansion beyond the long-wavelength approximation"
• (Li et al., 2018) "Origin of the anapole condition as revealed by a simple expansion beyond the toroidal multipole"
• (Nemkov et al., 2018) "Electromagnetic sources beyond common multipoles"
• (Riccardi et al., 2022) "Multipolar expansions for scattering and optical force calculations beyond the long wavelength approximation"
• (Ross, 2012) "Multipole expansion at the level of the action"
• (Zschocke, 2014) "A detailed proof of the fundamental theorem of STF multipole expansion in linearized gravity"
• (Gong et al., 2021) "Equivalence between angular spectrum-based and multipole expansion-based formulas of the acoustic radiation force and torque"
• (Agre, 2014) "Multipole expansions in quantum radiation theory"
• (Amalberti et al., 2023) "Multipole expansion at the level of the action in $d$-dimensions"
• (Seraj, 2016) "Multipole charge conservation and implications on electromagnetic radiation"
• (Compère et al., 31 Mar 2025) "Electromagnetic multipole expansions and the logarithmic soft photon theorem"
• (Gullans et al., 2017) "High-Order Multipole Radiation from Quantum Hall States in Dirac Materials"
• (Boudec et al., 2024) "Cartesian and spherical multipole expansions in anisotropic media"
• (Casini et al., 2024) "A unifying polarization formalism for electric- and magnetic-multipole interactions"
• (Wu et al., 2018) "Multipole analysis in the radiation field for linearized $f(R)$ gravity with irreducible Cartesian tensors"
• (Wu et al., 2018) "Multipole analysis for linearized $f(R)$ gravity with irreducible Cartesian tensors"

The multipole expansion thus unifies the treatment of the radiation field across physical contexts, enables systematic extraction of observable signatures, and provides a direct formal link between source structure and far-zone behavior. Its theoretical extensions underpin the effective field theory analysis of radiation in both electromagnetic and gravitational sectors.