I am attempting to look at some Liouville type theorems via a Fourier analysis approach and after looking at a baby problem I seem to be very confused. I assume this doesn't count as a research problem (but the original questions might be).
So here is the baby problem. Consider $ \Delta \phi(x,y)=0$ for $ (x,y) \in R^2_+=R \times (0,\infty)$ with $ \phi(x,0)=0$. I will impose the a gradient bound like $ | \nabla \phi(x,y)| \le C y^{-\alpha}$ where $ \alpha \in (0,1)$. So $ | \phi(x,y)| \le C_1 y^{1-\alpha}$.
So without the gradient bound one clearly has solutions like $ \phi(x,y)= C y$ and after imposing the bound one sees they must have $ C=0$. So with this explicit solution if we take a Fourier transform in $x$ we arrive at $ \hat{\phi}(\xi,y)= C y \delta_0(\xi)$ where $ \delta_0$ is the Dirac mass at $0$.
Okay we now move on to the equation and assume $ \phi$ is generic (so not the explicit solution above). Taking a Fourier transform in $x$ we arrive at $$ \partial_{yy} \hat{\phi}(\xi,y) = \xi^2 \hat{\phi}(\xi,y) $$ and solving this ode gives $$ \hat{\phi}(\xi,y)= A(\xi) e^{\xi y} + B(\xi) e^{-\xi y} $$ for $ \xi \neq0 $ and $ \hat{\phi}(0,y)=C_1 y + C_2$. (I am being quite careful here since I can't tell where the problem is coming from). So now if we impose the boundary condition we get $B(\xi)=- A(\xi)$ and $ C_2=0$. From this formula i don't see how one is getting the possible solution with the Dirac mass.