So I know this has been asked before, but I really do not get why the divergence of the gravitational field $\mathbf{g}$ that exists due to a point mass $M$ at the origin is zero. I mean, I get it mathematically, but not the intuition. So take a point mass $M$ at the origin. For any mass $m$ with position vector $\mathbf{r}$ we have $\mathbf{g}(\mathbf{r})=\frac{-GMm}{|\mathbf{r}|^3}\mathbf{r}$. Now if you compute $\nabla\cdot\mathbf{g}(\mathbf{r})$ then this can be easily shown to be zero. My question is what does this mean? If you 'draw' the vector field then you get a lot of arrows pointing at the origin and the magnitude of these vectors decrease the further away from the origin we get. Naively you might think we have a vector field with negative divergence, but this is offset by the magnitude of the gravitational field 'blowing up' as we approach the origin. So I get why there is a problem, but not the intuition of what the statement is saying.
The best I have come up with is this. If we think of the situation where we have a spherically symmetric mass distribution of radius $R$, then the above formula holds provided $|\mathbf{r}|>R$ and we can use the fact that this is a special case of Gauss' law of gravitation, $\nabla\cdot\mathbf{g}=-4\pi G\rho$ where $\rho$ is the mass density at that point where you compute $\nabla\cdot\mathbf{g}$. I have read that if we are 'outside' our spherical mass then $\rho=0$ (so we are in empty space) but why is the mass density zero. I thought when we derive $\mathbf{g}$ above we take a mass $m$ with position vector $\mathbf{r}$, so at that point we have a mass $m$. Therefore, surely the mass density isn't zero. Or do we mean the mass density of the thing that is causing the gravitational field? So in the point mass at the origin case, $\rho=0$ except at the origin (which is infinite as we have a point mass). In the spherically symmetric mass distribution case $\rho=0$ as soon as we go beyond the radius of the mass.
Edit: Upon further thought, could it also be linked with the following intuitive idea. If a mass $m$ is attracted to $M$ then it will pick up speed as it approaches $M$. As it won't suddenly come to a stop then it will continue travelling further on beyond $M$. So if you imagine a sphere encompassing our mass $M$, then for every arrow going into this surface, there is a corresponding arrow going out. So when we are outside this surface, the net divergence is zero? Or am I completely barking up the wrong tree here?


