Gaussian variational kernel for a gradient quartic interaction ($J(q)=\bar{a}+\bar{\kappa }{q}^2+\gamma q^4$)

For $H=\int [a{\varphi }^2/2+\kappa |\nabla \varphi {|}^2/2+\gamma ({\nabla }^2\varphi {)}^2/2+B{\varphi }^2|\nabla \varphi {|}^2]$ and an even positive trial kernel $J(q)$, define $S_0={\sum }_q^{+}1/J(q)$ and $S_2={\sum }_q^{+}{q}^2/J(q)$ using a positive-wavevector sum for a real field. Isserlis theorem gives $⟨H_4{⟩}_0=4B{S}_0{S}_2/V$. The Feynman-Bogoliubov inequality therefore has stationary kernelsThus both the mass and gradient coefficient shift. In the thermodynamic limit, their self-consistency equations areA coarse-graining ultraviolet cutoff $\Lambda$ is necessary for the first integral in three dimensions: its large-$k$ radial integrand tends to a constant. Physical solutions require $J(q)>0$.