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 [itex]f[/itex] 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).
[itex]\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}[/itex]
We shrink the radius of the sphere down to zero. We argue as the radius becomes arbitrarily small, value of [itex]\mathbf{0}[/itex], [itex]f(\mathbf{r})[/itex] is essentially constant within the sphere so we can take it outside the integral. So we are left with
[itex]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)[/itex].
Next we use the divergence theorem, noting [itex]\boldsymbol{\nabla} \frac{1}{r} = -\frac{\hat{\mathbf{r}}}{r^2}[/itex]. The dot product with the unit normal to the sphere is the same everywhere, [itex]-\left.\frac{\hat{\mathbf{r}} \cdot \mathbf{n}}{r^2} \right|_{r =a} = -\frac{1}{a^2}[/itex].
[itex]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}[/itex]
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 [itex]4\pi a^2\left(-\frac{1}{a^2} \right)f(\mathbf{0}) = -4\pi f(\mathbf{0})[/itex].
So we determined
[itex]\int \int \int dV f(\mathbf{r})\boldsymbol{\nabla}^2 \frac{1}{r} = -4\pi f(\mathbf{0})[/itex].