Fixed-mean-number Bogoliubov variational energy ($\mathcal{E}=\frac{gN}{2}+{\sum }_{k\ne 0}[({ϵ}_k+g){\sinh }^2{\theta }_k-g\sinh {\theta }_k\cosh {\theta }_k]$)

For a weakly depleted homogeneous Bose gas with repulsive contact interaction, set $g=uN/V$ and ${ϵ}_k={\hslash }^2{k}^2/(2m)$. A zero-mode coherent state times a two-mode squeezed vacuum for each opposite-momentum pair has mean particle number $N=N_0+{\sum }_{k\ne 0}{\sinh }^2{\theta }_k$. Substitution of $N_0=N-\sum {\sinh }^2{\theta }_k$ into the leading condensate interaction energy, followed by a consistent quadratic truncation, gives the displayed variational energy. Omitting this substitution loses the subtraction embodied in the number-conserving Bogoliubov quadratic Hamiltonian.

For each mode set $A_k={ϵ}_k+g$ and $E_k=\sqrt{A_k^2-g^2}$. If $A_k>|g|$, choose $\tanh (2{\theta }_k^{\ast })=g/A_k$. The mode contribution equals $\frac{1}{2}{E}_k\cosh [2({\theta }_k-{\theta }_k^{\ast })]-A_k/2$, proving a unique minimum $(E_k-A_k)/2$. At fixed ultraviolet cutoff this lowers the energy relative to the uncorrelated condensate. The approximation requires small quantum depletion of a Bose condensate; an unregulated three-dimensional ground-state energy also requires renormalization of a three-dimensional dilute-gas contact coupling.