0

I want to calculate

$\nabla\cdot(\frac{\hat{r}}{r^2})$

Naively applying the expression for divergence in spherical coordinates gives

$\nabla\cdot(\frac{\hat{r}}{r^2})=\frac{1}{r^2}\frac{\partial}{\partial r}(r^2 *\frac{1}{r^2})=\frac{1}{r^2}\frac{\partial}{\partial r}(1)=0$

However, based on the context in which this problem arose, I believe there should be a delta function. You can get this by noticing:

$\nabla^2\frac{1}{r}=\nabla\cdot(\nabla\frac{1}{r})=-\nabla\cdot (\frac{\hat{r}}{r^2})$

Thus

$\nabla\cdot(\frac{\hat{r}}{r^2})=-\nabla^2\frac{1}{r}=4\pi\delta(\vec{r})$

Can anyone point out where the error is? I suspect there is an issue cancelling out the $r^2$ in the first method due to that covering up the zero-in-denominator issue, but I would like to understand precisely what goes wrong.

mkn
  • 59

1 Answers1

0

The divergence is indeed $0$ on $\Bbb{R}^3\setminus\{0\}$, so there's no error there. The only reason the dirac delta appears is because of the desire to "force" the divergence theorem to hold true even when it classically doesn't apply anymore. In other words, we WANT the equation $\int_{\Omega}\text{div}(F)\,dV = \int_{\partial \Omega}F\cdot n\,dA$ to hold true. The surface integral for this radial vector field is ALWAYS $4\pi$ for any surface which "encloses" the origin. Therefore, the only way this can hold true is if we define $\text{div}(F)=4\pi \delta$

peek-a-boo
  • 55,725
  • 2
  • 45
  • 89
  • Why the factor of $4 \pi$? It should be just $\text{div}(F)= \delta$ since u get the $4 \pi$ from the $sin(\theta)d \theta d \phi$ – Sagigever Jun 05 '21 at 11:46
  • @Sagigever If $\text{div}(F)=\delta$ then for any $\Omega$ containing the origin in its interior, $\int_{\Omega}\text{div}(F),dV = 1$, but we want it to equal $4\pi$ since $\int_{\partial \Omega}F\cdot n,dA=4\pi$. Hence, we need $\text{div}(F)=4\pi \delta$. – peek-a-boo Jun 05 '21 at 11:58
  • But $\delta(r)$ it doesnt depend on the integral of $\theta$ and $\phi$ according to my understanding. those integrals will give u $4\pi$ (assuming u take a sphere around the origin) – Sagigever Jun 05 '21 at 12:40
  • @Sagigever I'm talking about the thee dimensional dirac $\delta$, i.e such that if you integrate over any region containing the origin, you get $1$.. – peek-a-boo Jun 05 '21 at 14:16