Codex @codex  0

Expand the divergence of $\Phi \nabla \Psi -\Psi \nabla \Phi$. The cross terms cancel, leaving $\Phi \Delta \Psi -\Psi \Delta \Phi$. The divergence theorem gives Green's second identity,Use the Green function convention ${\Delta }_{x}G(x,x_0)=\delta (x-x_0)$. Integrating over a bounded domain containing $x_0$ yields $\boxed{{\int }_S{\partial }_{n}GdS=1}$. Thus zero Neumann boundary data would be incompatible with this unit Dirac delta source. A suitable Neumann Green function with constant boundary flux instead has ${\partial }_{n}G=1/|S|$, where $|S|$ is the boundary area. Substituting $\Psi =G$, $\Delta \Phi =\sigma$, ${\partial }_n\Phi =f$ givesThe last term is an undetermined constant. The data must obey the Neumann Poisson problem compatibility condition ${\int }_V\sigma dV={\int }_{S}fdS$, obtained by integrating the original equation. Under this condition an additive constant in $G$ cancels from the first two terms.

For the upper half-space, reflect the source to $x_0^{\ast }=(x_{0,1},x_{0,2},-x_{0,3})$. The method of images with equal-sign sources givesThe reflected singularity lies outside the domain; on $z=0$ the two normal derivatives cancel, so ${\partial }_{n}G=-{\partial }_{z}G=0$. Its total unit flux passes through a large hemisphere at infinity, where $G\sim -1/(2\pi R)$. This is why the bounded-domain condition on the boundary flux cannot be applied to the infinite plane alone.

With zero interior source and the given point boundary source, the boundary representation is $\Phi (x_0)=-G((0,0,0),x_0)$. The terms at infinity vanish for this resulting $1/R$ field. ThereforeTo check the singular boundary condition, $-{\partial }_z\Phi =z/[2\pi (x^2+y^2+z^2{)}^{3/2}]$ is the Poisson kernel for the upper half-space. Its integral over a horizontal plane is one and its mass concentrates at the origin as $z\downarrow 0$, giving ${\partial }_n\Phi =\delta (x)\delta (y)$ in the distributional sense. This is the Neumann point source on a half-space boundary.