 Solution

The differential Gauss law is $\nabla \cdot \mathbf{E}=\rho /{ϵ}_0$. On a large ball $B_R$, the product rule and $\nabla \varphi =-\mathbf{E}$ giveThe divergence theorem consequently expresses the electrostatic energy asFor a localized charge distribution with the electric potential chosen to vanish at infinity, $\varphi =O(R^{-1})$ and $\mathbf{E}=O(R^{-2})$, so the surface term is $O(R^{-1})$ and tends to zero. More generally, it suffices to assume that this boundary integral vanishes. Passing to all space givesThe physical decay and potential reference are relevant: $\mathbf{E}\to 0$ by itself does not specify the additive constant in $\varphi$. Adding a constant $C$ changes the charge-based expression by $CQ/2$ when the total charge is $Q$, while leaving the field unchanged. The zero-at-infinity potential is the usual reference implicit in this all-space electrostatic energy formula.