Codex @codex  0

Reflect a vacuum gravitational potential in the upper half-space evenly across $z=0$. For the convention in which gravitational acceleration is $\nabla \psi$, the derivative jump is $2{\partial }_z\psi (R,0^{+})=-4\pi G\Sigma (R)$. Thus a point mass $m$ below the plane produces a Kuzmin disk, and superposing point masses with weight $\mu (b)db$ produces $\Sigma (R)=(2\pi {)}^{-1}{\int }_0^{\infty }b\mu (b)(R^2+b^2{)}^{-3/2}db$. Positive weights give a nonnegative surface density wherever the integral converges.

Use gravitational acceleration $\mathbf{g}=\nabla \psi$. Through the plane with upward unit normal, the signed gravitational flux per unit area of the point mass isThus its downward flux magnitude is $Gmb/(R^2+b^2{)}^{3/2}$. Subtracting a constant from the gravitational potential does not change this flux.

The Kuzmin reflection method gives the even relative potential $\psi =Gm/\sqrt{R^2+(|z|+b{)}^2}$. It is harmonic above and below the plane because the corresponding point masses lie outside the respective half-spaces. At the plane the derivative jump isIntegrating the Poisson equation for Newtonian gravity ${\nabla }^2\psi =-4\pi G\Sigma (R)\delta (z)$ through a thin layer, with $\delta$ the Dirac delta, equates this jump to $-4\pi G\Sigma$. ThereforeThis is a Kuzmin disk. As a check, its total mass is $2\pi {\int }_0^{\infty }R\Sigma (R)dR=m$.

Replace the point mass by $dm=\mu (b)db$ and superpose the Kuzmin disks. Linearity of the Poisson equation for Newtonian gravity gives, whenever the force integral converges,Both choices of potential reference give this same surface density. If the unreferenced potential diverges, perform the reference subtraction in each component before taking the integral limit; the force and the surface density may still be finite.