The Sobolev embedding theorem in dimension three, combined with the preceding estimate, gives
One direct justification of this embedding is a bounded Sobolev extension operator into , followed by Fourier inversion and the Cauchy-Schwarz inequality: bounds the inverse Fourier integral by the norm.
For radial forcing, the preceding part makes and . Set , and . The radial equation gives . Both zero boundary values imply
Moreover . Since , this proves the radial Poisson gradient estimate:
For a nonradial counterexample choose in the shell and a nonconstant . For sufficiently small , set and . These are smooth with support away from the boundary. Scaling gives
Thus the failure of L2-to-Linfinity Poisson gradient bounds in three dimensions is explicit, despite the valid and value-supremum estimates.