Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

How can one prove the following identity: $$ V_k(r_1, r_2) = {2k+1\over 2 r_1 r_2}\int_{|r_1 - r_2|}^{r_1+r_2} e^{-{r\over D}} P_k\left(r_1^2 - r^2 + r_2^2 \over 2 r_1 r_2 \right) d r= $$ $$ =(2k+1) {I_{k+{1\over2}}\left(r_<\over D\right) \over \sqrt{r_<}} {K_{k+{1\over2}}\left(r_>\over D\right) \over \sqrt{r_>}} $$ ? Here $r_<=\min(r_1, r_2)$, $r_>=\max(r_1, r_2)$, $I_\nu$ and $K_\nu$ are the modified Bessel functions of the first and second kind, $P_k(x)$ are Legendre polynomials and $k=0, 1, 2, 3, \dots$.

I verified the identity numerically, so I know that it works, but I didn't figure out how one could prove it. I know it has something to do with Gegenbauer's addition theorem and the equation 10.23.8 at [1]. The $V_k(r_1, r_2)$ is just the coefficient in the Legendre expansion of: $$ {e^{{|{\bf r}_1 - {\bf r}_2|}\over D}\over |{\bf r}_1 - {\bf r}_2|} =\sum_{k=0}^\infty V_k(r_1, r_2) P_k(\cos\theta) $$ Note: when $D\to\infty$, then $V_k$ becomes ${r_<^k \over r_>^{k+1}}$, which follows from the properties of the Bessel functions and this is just the well known multipole expansion [2].



share|cite|improve this question
Cross-listed to – Qmechanic Oct 18 '12 at 20:55

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.