Codex @codex  0

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.

A function assigning energy to the separation of two particles. The common symmetric case is an even function, $W(z)=W(-z)$. Its average over pairs defines an interaction energy. If the Hessian matrix is bounded, the gradient is globally Lipschitz and grows at most linearly; its convolution with a probability measure therefore needs only a finite first absolute moment, even though the potential itself can require a second absolute moment.

Normalize the specified state by $⟨0|0⟩=1$ and write $⟨A⟩=⟨0|A|0⟩$. Its zero mode is a coherent state, so any normally ordered product of zero-mode creation operators and annihilation operators can be evaluated using ${\psi }_0\mapsto \sqrt{N_0}$ and ${\psi }_0^{†}\mapsto \sqrt{N_0}$. In particular $⟨({\psi }_0^{†}{)}^2{\psi }_0^2⟩=N_0^2$, rather than $N_0(N_0-1)$, which would apply to a zero-mode number state.

Sort the momentum-space contact-interaction Hamiltonian by its number of zero-mode factors. Four such factors give $u{N}_0^2/(2V)$. Three zero-mode factors and one nonzero momentum violate momentum conservation, so no such term survives. With two zero-mode factors there are three possibilities:

The sums over $\mathbf{k}\ne 0$ include both members of each opposite-momentum pair, consistently with the overall factor $1/2$. Terms with one zero-mode factor contain three nonzero-mode factors; their expected values vanish in the centered multimode squeezed vacuum. Indeed each $\psi$ is linear in $\phi ,{\phi }^{†}$, and an odd product has zero vacuum expectation value. Terms with no zero-mode factors are independent of $N_0$ at fixed ${\theta }_k$. ThereforeTo evaluate these vacuum expectation values, use the inverse Bogoliubov transformation and its Hermitian adjoint:Since ${\phi }_{\mathbf{k}}|0⟩=0$, the only nonzero two-factor vacuum expectation value is $⟨{\phi }_{\mathbf{k}}{\phi }_{\mathbf{q}}^{†}⟩={\delta }_{\mathbf{k},\mathbf{q}}$. The bosonic squeezed-vacuum contractions are consequentlyIn particular the three requested quantities are $-\sinh {\theta }_k\cosh {\theta }_k$, $-\sinh {\theta }_k\cosh {\theta }_k$ and ${\sinh }^2{\theta }_k$, respectively. The minus sign in the pair vacuum expectation values comes from the inverse Bogoliubov transformation and is essential to lowering the interaction energy. Substitution yields

The coherent state in the zero mode has mean occupation number $N_0$. The other occupation numbers follow from part (c), giving exactlyImpose the fixed mean particle number $⟨\hat{N}⟩=N$ by choosingHere $D$ is the quantum depletion of a Bose condensate. The approximation requires $0\le D\ll N$, so $N_0$ remains positive and macroscopically large. The specified coherent state and multimode squeezed vacuum are not exact eigenstates of the total number operator; $N$ in this variational construction is their prescribed mean particle number.

Write $P={\sum }_{\mathbf{k}\ne 0}{s}_k{c}_k$. Substituting into the retained part of the interaction energy before truncating givesThe last two terms quadratic in $D,P$ have no explicit large-$N$ factor when the Bogoliubov transformation parameters are held fixed, and are discarded along with the nonzero-mode quartic contribution already omitted in part (c). This is the same truncation that underlies the number-conserving Bogoliubov quadratic Hamiltonian. Thus, defining the physical number density correctly as $\rho =N/V=N/L^3$,Replacing $N_0$ by $N$ at the outset would give the wrong coefficient $2{\sinh }^2{\theta }_k$: expanding the leading condensate interaction energy supplies the necessary subtraction.

The printed definition $\rho =N/L$ is an error. In three dimensions the contact interaction coupling $u$ has units of energy times volume, so $uN/L$ does not even have units of energy. There is a direct counterexample with no depletion: when ${\theta }_k=0$, the specified coherent state has exactly $⟨H_{\mathrm{int}}⟩=u{N}^2/(2L^3)$. In fixed length and energy units, $L=2$, $u=1$ and $N=100$ give $625$, whereas the printed $u\rho N/2$ with $\rho =N/L$ gives $2500$. The requested formula and the minimization in part (e) hold with $\rho =N/L^3$.

Put ${ϵ}_k={\hslash }^2{k}^2/(2m)>0$ for $\mathbf{k}\ne 0$ and $g=u\rho$, using $\rho =N/L^3$. Since the zero mode has zero kinetic energy, the exact kinetic energy in this trial state isTogether with part (d), the fixed-mean-number Bogoliubov variational energy isFor repulsive contact interaction, $u>0$, each ${ϵ}_k+g>g>0$. Using the double-angle hyperbolic-function identities,soThe stationary point therefore satisfiesThis is a minimum, not just a vanishing derivative. Define the positive Bogoliubov quasiparticle dispersionAt the displayed solution ${\theta }_k^{\ast }$, $\cosh 2{\theta }_k^{\ast }=({ϵ}_k+g)/E_k$ and $\sinh 2{\theta }_k^{\ast }=g/E_k$. HenceThe hyperbolic cosine has its unique minimum at zero, proving global minimization of each retained term. Treating equal-magnitude wavevectors together simply multiplies the same function by their multiplicity. At a fixed ultraviolet cutoff, the resulting energy change from $\theta =0$ isIndeed $E_k<{ϵ}_k+g$. More directly, for small positive $\theta$,The pair vacuum expectation value lowers the interaction energy at first order in $\theta$, whereas the nonzero occupation number and its kinetic energy cost start at second order. The two-mode squeezed vacuum admixes opposite-momentum pairs into the condensate, introducing spatial correlations that reduce the chance of two bosons coinciding. These correlations lower repulsive contact interaction enough to outweigh the added kinetic energy. A small quantum depletion of a Bose condensate can therefore lower the energy below the uncorrelated condensate value.

Two qualifications delimit this result. First, the paper does not state the sign of $u$. More generally a finite stable minimum exists only when ${ϵ}_k+g>|g|$ for every retained nonzero mode; repulsion ensures this. For example, with $g=-{ϵ}_k$ in some mode, $f_k(\theta )={ϵ}_k\sinh (2\theta )/2$ tends to $-\infty$ as $\theta \to -\infty$, so an unconditional minimum claim for arbitrary $u$ is false. Also the formula minimizes the truncated Bogoliubov theory energy, not the exact quartic Hamiltonian including its discarded terms, and it must remain in the regime of small quantum depletion of a Bose condensate.

Second, a three-dimensional unregulated contact interaction cannot give a finite absolute beyond-leading-order ground-state energy by this bare sum. For large $k$, $E_k-{ϵ}_k-g=-g^2/(2{ϵ}_k)+O({ϵ}_k^{-2})$, whose three-dimensional sum has an ultraviolet divergence. A finite ultraviolet cutoff, or renormalization of a three-dimensional dilute-gas contact coupling matched to a physical scattering length, is needed. This does not alter the modewise minimization or the finite-cutoff lowering mechanism just derived.

Use an independent copy $(x^{\prime },y^{\prime })$ of $(x,y)$ under the product measure $\Pi \otimes \Pi$. Write $q=y-x$, $q^{\prime }=y^{\prime }-x^{\prime }$, andSince the pushforward measure of $\Pi \otimes \Pi$ under $(z_t,z_t^{\prime })$ is ${\mu }_t\otimes {\mu }_t$, the interaction energy isTwo applications of the chain rule giveThe factor $1/2$ is retained. This classical identity requires differentiation under the integral sign to be justified, just as in part (b). Bounded Hessian matrix of $W$ is sufficient, since $|q-q^{\prime }{|}^2\le 2|q{|}^2+2|q^{\prime }{|}^2$ is integrable; compactly supported endpoints are another sufficient case. The printed quadratic-growth condition alone does not provide this domination.

The failure can occur for the interaction energy too. Take the even interaction potential $W(x)=1-\cos (2\pi x^3)$ and replace the endpoint probability measures in part (b) by $\frac{1}{2}{\delta }_0+\frac{1}{2}\mu$ and $\frac{1}{2}{\delta }_0+\frac{1}{2}\nu$. Keep the origin fixed and translate every other atom by one. This coupling is increasing on the line, hence is optimal: its assigned values extend to a nondecreasing transport map, which is the gradient of a convex function. Differences between translated atoms remain integers, so their mutual interaction vanishes. Only interactions with the origin contribute, giving $\mathcal{W}[{\mu }_t]=F(t)/4$, with $F$ from part (b). It therefore also lacks a finite classical second derivative at $t=1/2$.

Define the barycenter of a measure by $m_{\mu }=\int xd\mu (x)$. Under $\Pi$, the mean displacement is $\bar{q}=m_{\nu }-m_{\mu }$, and optimality gives $\int |q{|}^{2}d\Pi =D^2$. Expanding the square under the product measure yields the variance identity for independent copies:In the classical regime of part (c), the lower bound on the Hessian matrix of the interaction potential therefore givesFor equal barycenters, this becomes

For the full stated class of probability measures, prove the displacement convexity inequality directly. Pointwise strong convexity of $W$ givesIntegrating, including the $1/2$ in the interaction energy, and using the displayed variance identity for independent copies givesAll terms are finite by quadratic growth and finite second absolute moments. Equal barycenters give the claimed $\lambda$-displacement convexity. Also $m_{{\mu }_t}=(1-t)m_{\mu }+t{m}_{\nu }$, so this metric geodesic remains in the prescribed fixed-mean class. Applying the same argument to subintervals proves the second derivative bound in the sense of distributions even when a classical derivative has not been justified. The mean restriction is natural: translating all particles together leaves their pairwise differences and hence the interaction energy unchanged.

Choose a uniform lower bound $\kappa \in \mathbb{R}$ for the Hessian matrix of the confinement potential, so $D^{2}V(x)\ge \kappa I$ for every $x$. A sufficient condition isIndeed, convexity of $V(x)-\kappa |x{|}^2/2$ gives the confinement energy inequality from part (b), now with $\kappa$ possibly negative. Adding it to the equal-mean interaction energy inequality from part (d) givesThis proves positive displacement convexity on every class with fixed barycenter. In particular, if $D^{2}V\ge -KI$ with $K\ge 0$, the sufficient condition is $\lambda >K$ and one may take $\tilde{\lambda }=\lambda -K$. Wherever classical second derivatives exist, the equivalent differential calculation is ${\mathcal{E}}^{\prime \prime }[{\mu }_t]\ge (\kappa +\lambda )D^2$; the integrated proof also covers the printed growth assumptions.