Past exam of the advanced quantum theory: path integrals and many-particle physics course of the University of Oxford / 2024 / 2

Let $V=L^3$ be the box volume. The number operator is the sum of the occupation numbers of all momentum eigenstates, including the zero mode:In the Fourier series for the contact interaction, the two creation operators supply phase $e^{-i({\mathbf{k}}_1+{\mathbf{k}}_2)\cdot \mathbf{r}}$ and the two annihilation operators supply $e^{i({\mathbf{k}}_3+{\mathbf{k}}_4)\cdot \mathbf{r}}$. Fourier orthogonality with periodic boundary conditions givesThe four normalization factors together contribute $V^{-2}$, leaving the momentum-space contact-interaction HamiltonianThe Kronecker delta enforces momentum conservation. An equivalent expression that resolves the momentum transfer isBoth are already in normal ordering; no extra commutator constant is introduced.

Taking the Hermitian adjoint gives ${\phi }_{\mathbf{k}}=c_k{\psi }_{\mathbf{k}}+s_k{\psi }_{-\mathbf{k}}^{†}$. For nonzero wavevectors, the canonical commutation relations implyWhen the Kronecker delta is nonzero, $k=q$, so the hyperbolic Pythagorean identity $c_k^2-s_k^2=1$ givesFor the two annihilation operators, the only potentially nonzero terms areHere $\mathbf{q}=-\mathbf{k}$ implies equal magnitudes and hence equal $c,s$, because ${\theta }_k$ depends only on $|\mathbf{k}|$. ThusThe corresponding commutators of two creation operators vanish by the Hermitian adjoint. This verifies that the Bogoliubov transformation preserves the full bosonic canonical commutation relations.

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.