TSC said:
How do you prove that equation?
One may do so using the test function method.
We observe that the Laplacian of 1/r is zero every except possibly at the origin. Then we assume a test function f with nice properties. We then integrate, splitting the integral up into two regions - one a small sphere around the origin with radius a, and another everything else (which integrates to 0).
\int_0^{2\pi} d\phi \int_0^\pi d\theta \int_0^a \mathop{dr} r^2 \sin \theta f(\mathbf{r})\boldsymbol{\nabla}^2 \frac{1}{r}
We shrink the radius of the sphere down to zero. We argue as the radius becomes arbitrarily small, value of \mathbf{0}, f(\mathbf{r}) is essentially constant within the sphere so we can take it outside the integral. So we are left with
f(\mathbf{0})\int_0^{2\pi} d\phi \int_0^\pi d\theta \int_0^a \mathop{dr} r^2 \sin \theta\, \boldsymbol{\nabla} \cdot \left(\boldsymbol{\nabla} \frac{1}{r}\right).
Next we use the divergence theorem, noting \boldsymbol{\nabla} \frac{1}{r} = -\frac{\hat{\mathbf{r}}}{r^2}. The dot product with the unit normal to the sphere is the same everywhere, -\left.\frac{\hat{\mathbf{r}} \cdot \mathbf{n}}{r^2} \right|_{r =a} = -\frac{1}{a^2}.
f(\mathbf{0})\int \int \left(\boldsymbol{\nabla} \frac{1}{r}\right)\cdot \mathbf{n} \mathop{dS} =f(\mathbf{0})\int \int \left( -\frac{1}{a^2} \right)\mathop{dS}
Since it's the same everywhere on the sphere to obtain the answer we just multiply it by the surface area of the sphere, to obtain the final result 4\pi a^2\left(-\frac{1}{a^2} \right)f(\mathbf{0}) = -4\pi f(\mathbf{0}).
So we determined
\int \int \int dV f(\mathbf{r})\boldsymbol{\nabla}^2 \frac{1}{r} = -4\pi f(\mathbf{0}).