Skip to content

The Weakly Interacting Bose Gas

A homogeneous three-dimensional Bose gas is weakly interacting when its scattering-length gas parameter na3na^3 is small. The Gross–Pitaevskii equation and the leading energy density are saddle-point results built from the matched coupling g=4πa/mg=4\pi a/m; depletion and Lee–Huang–Yang terms are fluctuation corrections, not part of that saddle.

Required background. Use ideal condensation, short-range scattering matching, and collective-field transformations.

For a repulsive gas at wavelengths much longer than the interaction range r0r_0, take

S=dtd3x[iψtψψ22mVextψ2g2ψ4],g=4πam.S=\int\mathrm dt\,\mathrm d^3x\left[ i\psi^*\partial_t\psi-\frac{|\boldsymbol\nabla\psi|^2}{2m} -V_{\mathrm{ext}}|\psi|^2-\frac{g}{2}|\psi|^4 \right], \qquad g=\frac{4\pi a}{m}.

The last equality is a low-energy matching condition. A regulated bare coupling depends on the cutoff; inserting 4πa/m4\pi a/m before the ultraviolet subtraction in a loop calculation double counts renormalization. At the saddle, ψ=eiμtΦ(x)\psi=e^{-i\mu t}\Phi(\mathbf x) gives

[22m+Vext+gΦ2]Φ=μΦ,d3xΦ2=N0.\left[-\frac{\nabla^2}{2m}+V_{\mathrm{ext}}+g|\Phi|^2\right]\Phi=\mu\Phi, \qquad \int \mathrm d^3x\,|\Phi|^2=N_0.

For a uniform condensate, n0=Φ2n_0=|\Phi|^2, μ=gn0\mu=gn_0, and

EMFV=gn022.\frac{E_{\mathrm{MF}}}{V}=\frac{g n_0^2}{2}.

The energy cost of a density variation balances kinetic and interaction terms at the healing length

ξ=12mgn0=18πan0.\xi=\frac{1}{\sqrt{2mgn_0}}=\frac{1}{\sqrt{8\pi a n_0}}.

Thus the contact description requires r0ξr_0\ll\xi and the dilute expansion requires na31na^3\ll1. The Gross–Pitaevskii construction and its normalization are derived in Pethick and Smith 2008, chs. 5–6.

When the condensate varies over a scale RξR\gg\xi, neglecting the kinetic term gives

n0(x)=μVext(x)gΘ ⁣(μVext(x)).n_0(\mathbf x)=\frac{\mu-V_{\mathrm{ext}}(\mathbf x)}{g} \,\Theta\!\left(\mu-V_{\mathrm{ext}}(\mathbf x)\right).

This Thomas–Fermi profile fails in a boundary layer of thickness order ξ\xi and near vortices, where gradients are essential. It also assumes a local three-dimensional coupling; tight confinement changes the matching.

The first omitted contribution to the homogeneous energy is relative order na3\sqrt{na^3}. Quantum depletion has the same characteristic small parameter Pitaevskii and Stringari 2016, chs. 4–5. Consequently a Gross–Pitaevskii solution may be spatially accurate while its energy is not accurate enough for a cancellation problem such as a droplet. The Bogoliubov page supplies quadratic fluctuations, and the dilute-expansion page performs the renormalized next order.

Restore \hbar and verify the dimensions of gg, μ\mu, and ξ\xi in three dimensions.

Solution

g=4π2a/mg=4\pi\hbar^2a/m has dimensions energy times volume. Hence gngn is an energy and equals the mean-field chemical potential. Balancing 2/(2mξ2)=gn\hbar^2/(2m\xi^2)=gn gives ξ=/2mgn\xi=\hbar/\sqrt{2mgn}, a length.