The multipole expansion decomposes solutions to the diffusion equation in spherical coordinates into orthogonal series terms (spherical harmonics) characterized by degree $n$ and order $m$, representing physical configurations from monopoles ($n=0$) through dipoles, quadrupoles, and octopoles. These abstract mathematical components correspond directly to superpositions of point sources and sinks arranged with specific symmetries that eliminate lower-order moments (net source strength, dipole moment), resulting in radial decay rates proportional to $1/R^{n+1}$. This theory provides a formal mechanism for describing potential fields where boundary conditions require satisfying linear combinations of decaying harmonics without net energy flux from the origin.
The multipole expansion decomposes solutions to the diffusion equation in spherical coordinates into orthogonal series terms (spherical harmonics) characterized by degree $n$ and order $m$, represent…