# Laplace expansion (potential)

Jump to: navigation, search
See also Laplace expansion of determinant.

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 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

${\displaystyle {\frac {1}{|\mathbf {r} -\mathbf {r} '|}}=\sum _{\ell =0}^{\infty }{\frac {4\pi }{2\ell +1}}\sum _{m=-\ell }^{\ell }(-1)^{m}{\frac {r_{\scriptscriptstyle <}^{\ell }}{r_{\scriptscriptstyle >}^{\ell +1}}}Y_{\ell }^{-m}(\theta ,\varphi )Y_{\ell }^{m}(\theta ',\varphi ').}$

Here r has the spherical polar coordinates (r, θ, φ) and r' has ( r', θ', φ') with homogeneous polynomials of degree ℓ in (x,y,z). Further r< is min(r, r') and r> is max(r, r'). The function ${\displaystyle Y_{\ell }^{m}}$ is a normalized spherical harmonic function. The expansion takes a simpler form when written in terms of solid harmonics,

${\displaystyle {\frac {1}{|\mathbf {r} -\mathbf {r} '|}}=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }(-1)^{m}I_{\ell }^{-m}(\mathbf {r} )R_{\ell }^{m}(\mathbf {r} ')\quad {\hbox{with}}\quad |\mathbf {r} |>|\mathbf {r} '|.}$

## Derivation

One writes

${\displaystyle {\frac {1}{|\mathbf {r} -\mathbf {r} '|}}={\frac {1}{\sqrt {r^{2}+(r')^{2}-2rr'\cos \gamma }}}={\frac {1}{r{\sqrt {1+h^{2}-2h\cos \gamma }}}}\quad {\hbox{with}}\quad h\equiv {\frac {r'}{r}}.}$

We find here the generating function of the Legendre polynomials ${\displaystyle P_{\ell }(\cos \gamma )}$ :

${\displaystyle {\frac {1}{\sqrt {1+h^{2}-2h\cos \gamma }}}=\sum _{\ell =0}^{\infty }h^{\ell }P_{\ell }(\cos \gamma ).}$

Use of the spherical harmonic addition theorem

${\displaystyle P_{\ell }(\cos \gamma )={\frac {4\pi }{2\ell +1}}\sum _{m=-\ell }^{\ell }(-1)^{m}Y_{\ell }^{-m}(\theta ,\varphi )Y_{\ell }^{m}(\theta ',\varphi ')}$

gives the desired result.