Skip to content

Bogoliubov Theory and Bose Quasiparticles

Quadratic fluctuations about a dilute condensate mix particle creation with hole annihilation. A canonical Bogoliubov transformation diagonalizes that quadratic Hamiltonian and gives Ek=ϵk(ϵk+2gn0)E_k=\sqrt{\epsilon_k(\epsilon_k+2gn_0)}: phonons at kξ1k\xi\ll1, free particles at kξ1k\xi\gg1, and a gapless Goldstone mode enforced by symmetry.

Required background. Use the matched weak-gas saddle, the Goldstone pole argument, and Dyson/self-energy conventions.

Write a0N0a_{\mathbf0}\simeq\sqrt{N_0} in the symmetry-breaking representation and retain terms quadratic in ak0a_{\mathbf k\ne0}. With ϵk=k2/(2m)\epsilon_k=k^2/(2m) and μ=gn0\mu=gn_0,

K2=k0[(ϵk+gn0)akak+gn02(akak+akak)].K_2=\sum_{\mathbf k\ne0}\left[ (\epsilon_k+gn_0)a_{\mathbf k}^\dagger a_{\mathbf k} +\frac{gn_0}{2}\left(a_{\mathbf k}^\dagger a_{-\mathbf k}^\dagger+a_{\mathbf k}a_{-\mathbf k}\right) \right].

Set ak=ukbkvkbka_{\mathbf k}=u_k b_{\mathbf k}-v_k b^\dagger_{-\mathbf k} with uk2vk2=1u_k^2-v_k^2=1. Cancellation of anomalous terms yields

Ek=ϵk(ϵk+2gn0),uk2=12(ϵk+gn0Ek+1),vk2=12(ϵk+gn0Ek1).E_k=\sqrt{\epsilon_k(\epsilon_k+2gn_0)}, \quad u_k^2=\frac12\left(\frac{\epsilon_k+gn_0}{E_k}+1\right), \quad v_k^2=\frac12\left(\frac{\epsilon_k+gn_0}{E_k}-1\right).

The commutator condition is the normalization check; the positive root is selected by stability. At small momentum, EkckE_k\simeq ck with c=gn0/mc=\sqrt{gn_0/m}. At large momentum, Ek=ϵk+gn0+O(k2)E_k=\epsilon_k+gn_0+O(k^{-2}). These limits independently test both the factor of two and the chemical-potential subtraction Altland and Simons 2023, § 5.2, pp. 242–257.

The quasiparticle vacuum contains bare particles:

ndep=kvk2=83πn0n0a3n_{\mathrm{dep}}=\int_{\mathbf k}v_k^2 =\frac{8}{3\sqrt\pi}\,n_0\sqrt{n_0a^3}

in three dimensions at zero temperature. The small ratio ndep/n0n_{\mathrm{dep}}/n_0 is a self-consistency check, not a definition of superfluid density.

The density operator couples with amplitude ukvku_k-v_k, giving the leading static structure factor

S(k)=ϵkEk.S(k)=\frac{\epsilon_k}{E_k}.

Thus S(k)k/(2mc)S(k)\sim k/(2mc) in the phonon regime and tends to one in the particle regime. The compressibility sum rule reproduces c2=gn/mc^2=gn/m.

In Nambu language, exact gaplessness requires the Hugenholtz–Pines relation μ=Σ11(0)Σ12(0)\mu=\Sigma_{11}(0)-\Sigma_{12}(0), not the separate vanishing of either self-energy. The original theorem and its assumptions are given in Hugenholtz and Pines 1959, pp. 489–506.

This is the leading quadratic theory for na31na^3\ll1 and momenta below the interaction-range scale. It is not exact at strong depletion, in one dimension, or near a critical point. Using nn or n0n_0 in a subleading expression without matching the perturbative order can move terms between orders. The renormalized correction is derived next.

Verify from the formulas for uku_k and vkv_k that uk2vk2=1u_k^2-v_k^2=1 and that S(k)=ϵk/EkS(k)=\epsilon_k/E_k has the correct high-momentum limit.

Solution

Subtracting the two coherence factors cancels the common ratio and leaves one. Since Ek/ϵk=1+2gn0/ϵk1E_k/\epsilon_k=\sqrt{1+2gn_0/\epsilon_k}\to1 for kk\to\infty, S(k)1S(k)\to1, as required for resolving individual particles at short wavelength.

  • Alexander Altland and Ben Simons, Condensed Matter Field Theory, 3rd ed., Cambridge University Press (2023), § 5.2, pp. 242–257, doi:10.1017/9781108781244.
  • N. M. Hugenholtz and David Pines, “Ground-State Energy and Excitation Spectrum of a System of Interacting Bosons,” Physical Review 116 (1959) 489–506, doi:10.1103/PhysRev.116.489.