I am reading Galactic Dynamics by Binney and Termaine, and I'm struggling with the following passage.
According to Poisson's equation the resulting change in potential $\delta\Phi(\mathbf x)$ satisfies $\nabla^2(\delta\Phi)=4\pi G(\delta\rho)$, so $$\delta W=\frac{1}{4\pi G}\int\mathrm d^3\mathbf x\ \Phi\nabla^2(\delta\Phi) \tag{2.14}$$ Using the divergence theorem in the form $(\mathrm B.45)$, we may write this as $$\delta W= \frac{1}{4\pi G}\int\Phi\nabla(\delta\Phi) \cdot\mathrm d^2\mathbf S -\frac{1}{4\pi G}\int\mathrm d^3\mathbf x\ \nabla\Phi\cdot\nabla(\delta\Phi), \tag{2.15}$$ where the surface integral vanishes because $\Phi\propto r^{-1}$, $|\nabla\delta\Phi|\propto r^{-2}$ as $r\to\infty$, so the integrand $\propto r^{-3}$ while the total surface area $\propto r^2$. But $\nabla\Phi\cdot\nabla(\delta\phi)=\tfrac12\delta(\nabla\Phi\cdot\nabla\Phi)=\tfrac12\delta|(\nabla\Phi)|^2$. Hence $$\delta W= -\frac{1}{8\pi G}\delta\left(\int\mathrm d^3\mathbf x\ |\nabla\Phi|^2\right). \tag{2.16}$$
I am trying to understand the integral being carried out there. I understand why surface integral is vanishing as $r$ tends to infinity. But like same reasoning why volume integral is not vanishing though it has inverse proportionality with $r$.
Similarly there are two instances in the book where he reasons alike.
What is being compromised when volume integral is considered? Any help would be very great.