Codex @codex  0

Write $\mathcal{V}=L^3$ and $N_0=N-N_{\mathrm{ex}}$, where $N_{\mathrm{ex}}={\sum }_{\mathbf{k}\ne 0}{b}_{\mathbf{k}}^{†}{b}_{\mathbf{k}}$. Macroscopic occupation permits $b_0,b_0^{†}\simeq \sqrt{N_0}$, with relative commutator corrections of order $1/N_0$. Keep the fixed-total-number correction before replacing $n_0$ by $n$.

The four-condensate term is $V_0{N}_0(N_0-1)/(2\mathcal{V})$. Its leading expansion gives $V_0{N}^2/(2\mathcal{V})-n{V}_0{N}_{\mathrm{ex}}$. Terms with two condensate operators give normal scattering $n(V_0+V_{\mathbf{k}})b_{\mathbf{k}}^{†}{b}_{\mathbf{k}}$ and conversion of two condensate particles into opposite-momentum excitations, with coefficient $n{V}_{\mathbf{k}}/2$. The $n{V}_0$ normal contribution cancels against the change of condensate energy at fixed $N$. Hence the number-conserving Bogoliubov quadratic Hamiltonian isUse a real even interaction transform, $V_{\mathbf{k}}=V_{-\mathbf{k}}$, as for a real inversion-symmetric pair potential. The constant is the condensate interaction energy; the diagonal part combines particle kinetic energy with condensate scattering; the anomalous pair terms create and annihilate correlated opposite momenta while transferring particles to or from the condensate.

At fixed volume, counting replaced condensate operators gives quadratic vertices proportional to $N_0$, omitted cubic vertices proportional to $\sqrt{N_0}$ and quartic vertices with no condensate factor. There are also the finite-$N_0$ correction in $N_0(N_0-1)$ and corrections from the condensate-number expansion. Their physical control parameter is small depletion, $N_{\mathrm{ex}}/N\ll 1$, together with weak interaction; simply dropping them is not justified by low temperature in a strongly interacting gas. This is the weakly interacting Bose gas approximation.

For a dilute three-dimensional weakly interacting Bose gas with repulsive short-range scattering, use a contact coupling $g=4\pi {\hslash }^2{a}_s/m>0$ and $n{a}_s^3\ll 1$. In imaginary time ${\tau }_E$, the grand-canonical action isThe uniform saddle has $n_0=\mu /g$. Its real-time equation is the Gross–Pitaevskii equation; quantum depletion means condensate and total density need not coincide exactly.

Writing $\psi =\sqrt{n_0}+\delta \psi$ and retaining quadratic fluctuations gives a Nambu kernel with diagonal ${ϵ}_k+g{n}_0$, anomalous coupling $g{n}_0$ and ${ϵ}_k={\hslash }^2{k}^2/(2m)$. The Bogoliubov quasiparticle dispersion isPhase and density fluctuations form the gapless sound mode associated with broken global $U(1)$ symmetry. At large momentum the spectrum crosses towards particle behaviour; the healing length is $\xi =\hslash /\sqrt{2mg{n}_0}$. The depletion integral $\int d^{3}k{v}_k^2/(2\pi {)}^3$ is finite in three dimensions and gives $(n-n_0)/n\simeq 8\sqrt{n{a}_s^3}/(3\sqrt{\pi })$ to leading order. The contact coupling denotes a scattering-length-renormalized interaction; naive ultraviolet zero-point integrals require renormalization.

The condensate phase gives superfluid velocity ${\mathbf{v}}_s=(\hslash /m)\nabla \theta$ and quantized circulation $\oint {\mathbf{v}}_s\cdot d\mathbf{l}=2\pi \hslash \ell /m$. The Landau criterion, $v_c={\min }_k{E}_k/(\hslash k)$, gives $c_s$ for this homogeneous Bogoliubov spectrum. Boundaries and vortex nucleation can cause dissipation below that value, so it is not a universal measured critical speed. Condensation and superfluid stiffness are related but distinct; dimensional fluctuations can remove true condensate order while retaining superfluid responses.

Use the contact-interaction Bose-gas field integral for a homogeneous three-dimensional weakly interacting Bose gas with repulsive short-range interaction as the standard model. Let particle mass be $m_b$, scattering length $a_s>0$, physical coupling $g=4\pi {\hslash }^2{a}_s/m_b$ and gas parameter $n{a}_s^3\ll 1$. These are symbolic model definitions, not supplied numerical parameters. The grand-canonical coherent-state integral has periodic fields and actionThe contact term is the normal symbol of $g{\hat{\psi }}^{†2}{\hat{\psi }}^2/2$. In three dimensions a contact interaction needs a regulator and coupling matched to the measured scattering length; the displayed $g$ is the low-energy physical coupling. This field integral is useful because an ordered saddle and its quantum/thermal fluctuations follow from the same action.

Stationarity with respect to $\bar{\psi }$ gives the imaginary-time Gross–Pitaevskii equation. A uniform saddle exists at $\mu >0$ with ${\psi }_0=\sqrt{n_0}{e}^{i{\theta }_0}$ and condensate density $n_0=\mu /g$. Its arbitrary phase reflects global $U(1)$ symmetry; a selected condensate fixes a phase in the thermodynamic symmetry-breaking description. Write $\psi =\sqrt{n_0}+\eta$ after choosing that phase. Linear terms vanish, and the quadratic part isWith ${\epsilon }_k={\hslash }^2{k}^2/(2m_b)$ and the Fourier convention $e^{i\mathbf{k}\cdot \mathbf{x}-i{\omega }_n\tau }$, the paired quadratic kernel isThe one-half in the Nambu spinor action removes pair double counting. Analytic continuation of the inverse kernel gives the Bogoliubov dispersionAt high momentum it approaches the free-particle kinetic energy plus its mean-field shift. The crossover scale is the healing length $\xi =\hslash /\sqrt{2m_{b}g{n}_0}$. Gaussian integration adds a fluctuation determinant to the saddle free energy. Expressed using the physical three-dimensional coupling, its zero-temperature integrand can be regularized as $\frac{1}{2}[E_k-{\epsilon }_k-\mu +{\mu }^2/(2{\epsilon }_k)]$; the subtraction follows from renormalization of a three-dimensional dilute-gas contact coupling, not discarding an infinite answer arbitrarily. The thermal contribution is $k_{B}T\int d^{3}k/(2\pi {)}^3\log (1-e^{-\beta E_k})$. A calculation beyond mean field must keep that matching and its thermodynamic derivatives consistent.

The density-phase action of a Bose gas makes the physical content clear, $\psi =\sqrt{n_0+\delta n}{e}^{i\theta }$. Up to a periodic density total derivative, the long-wavelength quadratic action containsnumber-phase conjugacy pairs the density fluctuation and phase in the quadratic density-phase action. Dropping the density-gradient term only in the long-wave limit and completing the density Gaussian square gives the Bose-gas phase-only actionThis has the same sound pole. A term $i\hslash n_0{\partial }_{\tau }\theta$ is also present before restricting to the smooth zero-winding sector; although locally a total derivative, it cannot in general be ignored in topologically nontrivial sectors. The density and phase variables describe the same coupled Bogoliubov branch, not a second independent massive particle simply because the density has a static restoring cost.^[ref]

A phase gradient carries superflow ${\mathbf{v}}_s=(\hslash /m_b)\nabla \theta$. Its energy cost defines the superfluid number density. At leading mean-field order the stiffness uses $n_0$, but condensate density and superfluid number density are different observables: for a Galilean-invariant continuum gas at zero temperature the latter is the full density even when interactions deplete the condensate. Bogoliubov coefficients give $v_k^2=[({\epsilon }_k+g{n}_0)/E_k-1]/2$ and a small zero-temperature depletion of relative order $\sqrt{n{a}_s^3}$ in three dimensions. Large depletion or strong long-wavelength phase fluctuations would invalidate the small-fluctuation saddle approximation.

The Landau criterion for creating excitations is $v_L={\inf }_{k>0}{E}_k/(\hslash k)=c_s$ for this homogeneous dispersion. Thus repulsion produces a linear sound mode and a positive critical velocity from the Landau criterion, whereas an ideal quadratic gas has zero value of this criterion. Vortices, obstacles and nucleation can reduce an actual flow threshold below $c_s$; the criterion alone is not a complete theory of metastable flow. At finite temperature thermal quasiparticles reduce the superfluid number density. Dimensionality and finite size also matter: the uniform three-dimensional weak-depletion condensate expansion must not be transferred unchanged to a one- or two-dimensional thermodynamic system. The field integral provides systematic corrections only within its dilute, stable and sufficiently weak-fluctuation regime.

The Hartree product state underlying the Gross–Pitaevskii equation has the factorized one-body coherence of a Hartree condensate. Choosing its nodeless ground state orbital real givesIn particular, the Thomas–Fermi approximation for a condensate from part (d) predictsThe curve starts at $\mu /g$, decreases with zero slope at the origin, and reaches zero with an infinite slope at $R$. Restoring the kinetic energy in the full Gross–Pitaevskii equation rounds that edge and produces a nonzero tail. This corrects the Thomas–Fermi approximation for a condensate, but the resulting one-body density matrix still factorizes.

The additional missing physics is quantum depletion. Interactions scatter pairs of bosons into orbitals outside the Bose-Einstein condensate even at zero temperature. In terms of the true Bose-Einstein condensate orbital ${\phi }_0$ the exact one-body density matrix decomposes asThe single-orbital Hartree product state sets $N_0=N$ and omits ${\rho }_{\mathrm{ex}}$ altogether. In a weakly interacting Bose gas, Bogoliubov theory adds a comparatively short-range noncondensed contribution near the origin, together with reduced weight in the broad Bose-Einstein condensate component. The central value remains the total number density, not just the condensate density. Its numerical shift relative to the original Gross–Pitaevskii equation solution also depends on the interaction-induced redistribution of number density; it is not fixed by the condensate fraction alone.

A positive density curvature $g$ makes uniform static density fluctuations costly. It does not imply an additional independent gapped excitation: the first-order term makes density and phase conjugate. For a simple nonrelativistic weakly interacting Bose gas, their coupled quadratic dynamics yields one positive Bogoliubov spectrum, which is gapless at zero momentum. The word massive for the density field describes its static curvature in this setting.