This is the Gibbs--Bogoliubov--Feynman inequality. In units with , is the exact dimensionless Hamiltonian, here the full interacting free-energy functional, and is a freely chosen trial Hamiltonian. Its partition function and free energy are
The notation
means an expectation in the trial Gibbs ensemble. The bound is therefore
Choose the translationally invariant Gaussian variational approximation
The superscript means that one member of each pair is included; this is necessary because a real field obeys . The constrained zero mode is omitted.
For each independent complex mode,
Consequently,
Substitution into the Gibbs--Bogoliubov--Feynman inequality gives
which is the required upper bound.