Laplace expansion (potential)

In physics, the Laplace expansion of a 1/r - type potential is applied to expand Newton's gravitational potential or Coulomb's electrostatic potential. In quantum mechanical calculations on atoms the expansion is used in the evaluation of integrals of the interelectronic repulsion.

The expansion
The Laplace expansion is in fact the expansion of the inverse distance between two points. Let the points have position vectors r and r', then the Laplace expansion is
 * $$ \frac{1}{|\mathbf{r}-\mathbf{r}'|} = \sum_{\ell=0}^\infty \frac{4\pi}{2\ell+1}  \sum_{m=-\ell}^{\ell} (-1)^m \frac{r_^\ell }{r_^{\ell+1} } Y^{-m}_\ell(\theta, \varphi) Y^{m}_\ell(\theta', \varphi'). $$

Here   r has the spherical polar coordinates (r, &theta;, &phi;) and r' has ( r', &theta;', &phi;'). Further r&lt; is min(r, r') and r&gt; is max(r, r'). The function $$Y^m_{\ell}$$ is a normalized spherical harmonic function. The expansion takes a simpler form when written in terms of solid harmonics,
 * $$ \frac{1}{|\mathbf{r}-\mathbf{r}'|} = \sum_{\ell=0}^\infty  \sum_{m=-\ell}^{\ell} (-1)^m  I^{-m}_\ell(\mathbf{r}) R^{m}_\ell(\mathbf{r}')\quad\hbox{with}\quad |\mathbf{r}| > |\mathbf{r}'|, $$

where $$R^{m}_\ell$$ is a regular solid harmonic:
 * $$ R^m_{\ell}(\mathbf{r}) \equiv \sqrt{\frac{4\pi}{2\ell+1}}\; r^\ell Y^m_{\ell}(\theta,\varphi),  $$

and $$I^{m}_\ell$$ is an irregular solid harmonic:
 * $$ I^m_{\ell}(\mathbf{r}) \equiv \sqrt{\frac{4\pi}{2\ell+1}} \; \frac{ Y^m_{\ell}(\theta,\varphi)}{r^{\ell+1}} . $$

Derivation
The derivation of this expansion is simple. One writes
 * $$ \frac{1}{|\mathbf{r}-\mathbf{r}'|} = \frac{1}{\sqrt{r^2 + (r')^2 - 2 r r' \cos\gamma}} =   \frac{1}{r_ \sqrt{1 + h^2 - 2 h \cos\gamma}} \quad\hbox{with}\quad h \equiv \frac{r_}{r_} .   $$

We find here the generating function of the Legendre polynomials $$P_\ell(\cos\gamma)$$ :
 * $$ \frac{1}{\sqrt{1 + h^2 - 2 h \cos\gamma}} = \sum_{\ell=0}^\infty h^\ell P_\ell(\cos\gamma). $$

Use of the spherical harmonic addition theorem
 * $$ P_{\ell}(\cos \gamma) = \frac{4\pi}{2\ell + 1} \sum_{m=-\ell}^{\ell} (-1)^m Y^{-m}_{\ell}(\theta, \varphi) Y^m_{\ell}(\theta', \varphi') $$

gives the desired result.