# Multipole expansion (interaction)

In physics, two static, non-polarizable, electric charge distributions interact via Coulomb's law. When the charge distributions do not overlap the interaction can expanded in a convergent series of powers of 1/*R*. Here *R* is the distance between the charge distributions, for instance, the distance between the centers of charge of the two distributions. Because the factors of 1/*R* in this expansion are products of multipoles of the two charge distributions, this series is known as the **multipole expansion** of the interaction energy.

## Mathematical formulation

A most important application of the multipole series is in the field of intermolecular forces, and because molecules are distributions consisting of discrete charges (electrons and nuclei), we consider the interaction between two discrete charge contributions.
So, let us introduce two sets of point charges, one set {*q*_{ i} } clustered around a point *A* and one set {*q*_{ j} } clustered around a point *B*.
By Coulomb's law the total electrostatic interaction energy *U*_{AB} between the two distributions is

where ε_{0} is the permittivity of the vacuum and *r*_{ ij} is the distance between particle *i* and particle *j*.
This energy can be expanded in a power series— the **multipole expansion**—in the inverse distance of *A* and *B*.

In order to derive this multipole expansion, we write **r**_{XY} ≡ **r**_{Y} − **r**_{X}, which is a vector pointing from *X* towards *Y*. Note that

We assume that the two distributions do not overlap:

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

where and are irregular and regular solid harmonics, respectively. The translation of the regular solid harmonic gives a finite expansion,

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

Use of the definition of spherical multipoles *Q*^{m}_{l} and covering of the summation ranges in a somewhat different order (which is only allowed for an infinite range of *L*) gives finally

This is the **multipole expansion** of the interaction energy of two non-overlapping charge distributions which are a distance *R*_{AB} apart.
Here the notation

is used for the Clebsch-Gordan series. Since

this expansion is manifestly in powers of 1/*R _{AB}*. The function Y

^{m}

_{l}is a normalized spherical harmonic, depending on the unit vector along the vector

**R**

_{AB}. Recall in this connection that a unit vector is in one-to-one correspondence with two spherical polar angles.

Finally, note that the existence of non-vanishing multipole moments is a matter of symmetry. For instance, a neutral, spherically symmetric charge distribution does not have any non-vanishing multipole moment. The interaction between two such non-overlapping systems is zero. An example of this is the interaction between two noble gas atoms that are so far apart that the overlap between their electronic wave functions is negligible.