 Solution

Differentiate the bound with respect to one $J(q)$, using $\partial S_0/\partial J(q)=-J(q{)}^{-2}$ and $\partial S_2/\partial J(q)=-q^{2}J(q{)}^{-2}$:Thus an interior minimizing Gaussian variational kernel for a gradient quartic interaction hasUse the half-space mode density given in the question. In the thermodynamic limit, the coupled self-consistency equations areThe admissible solution must have $J(q)>0$. The coefficient $\gamma$ is unchanged because the interaction's derivative contributes only a constant and a $q^2$ term.

The formal continuum integrals require the coarse-graining ultraviolet cutoff $\Lambda$ in three dimensions: the first radial integrand tends to a constant at large $k$, so taking $\Lambda =\infty$ would leave $\bar{a}$ divergent. The low-wavevector fluctuation argument is independent of this short-distance regularization.