Consider the Laplacian $-\Delta$ on (in a suitable sense) twice differentiable functions subject to homogeneous Dirichlet boundary conditions $\mathscr{H}=\{f : f(0)=f(1)=0\}$. We can identify the Green's function of $-\Delta$ to be
$$\mathcal{G}(s, t) = \min(s, t) - st.$$
This Green's function will of course work for solving the Poisson equation on any linear subspace of $\mathscr{H}$, but it is the distributional form of the Green's function—we may be able to find a functional form of the Green's function that seem more appropriate for a given subspace.
For example, if we consider the subspace of functions $f\in\mathscr{H}$ such that $\int_0^1 f(t)\,\mathrm{d}t =0$, we can find a functional form of the Green's function on this subspace that obeys the integral condition:
$$\mathcal{G}(s, t) = \min(s, t) - st− 3(1 − s)s(1 − t)t.$$
Now for my question, if we consider the Green's function for the Laplacian on the functions $f\in\mathscr{H}$ which are also subject to homogeneous Neumann boundary conditions $f'(0)=f'(1)=0$, I cannot find a functional form of the Green's function that obeys both types of conditions. Does such a form exist?
Through some applications, I have an indication that the extra Neumann conditions are non-restricting in some sense, but it is all very unclear to me what is going on. Is this something that any of you have had experience with, or can any of you refer me to some relevant work?