research28 min studyLesson 67 in Physics

Bose condensation and Bogoliubov theory

Derive the excitation scale of a dilute repulsive Bose gas and separate condensation, coherence, and superfluid response.

On this page

What you’ll understand

  • State the dilute-gas assumptions behind a contact interaction.
  • Recover sound-like and particle-like limits of the Bogoliubov spectrum.
  • Identify what condensate occupation alone does not establish.

When many bosons occupy the same one-particle mode, a collective field becomes a useful organizing variable. Interactions make its fluctuations more interesting than those of an ideal gas: long-wavelength excitations can carry sound rather than follow a free-particle parabola. We will analyze a homogeneous three-dimensional dilute gas at zero temperature, where weak depletion gives a controlled starting point. This precise setting matters because low dimensions, strong coupling, and thermal fluctuations can change the answer.

Condensation is a statement about a density matrix

The one-body density matrix records correlations between removing a boson at one point and adding it elsewhere. Bose condensation means one eigenvalue of that matrix grows proportionally to particle number in the thermodynamic limit. For a uniform gas that mode is the zero-momentum state. This definition works even in a number-conserving description where the expectation value of the field operator vanishes. Choosing a condensate phase is a convenient symmetry-breaking description, not a claim that exact particle conservation has failed.

For short-range interactions at low energy, replace microscopic details by a coupling g=4πℏ2as/mg=4\pi\hbar^2a_s/m fixed by the scattering length as>0a_s>0. This replacement requires momenta below the inverse interaction range, and weak depletion in three dimensions is governed by nas3≪1na_s^3\ll1. A contact Hamiltonian needs an appropriate ultraviolet prescription beyond the leading formulas used here. The physical scattering length, rather than an arbitrary bare delta-function strength, is the measurable input.

ρ1(r,r′)=⟨ψ†(r′)ψ(r)⟩,g=4πℏ2asm,nas3≪1\rho_1(\mathbf r,\mathbf r')=\langle\psi^\dagger(\mathbf r')\psi(\mathbf r)\rangle,\qquad g=\frac{4\pi\hbar^2a_s}{m},\qquad na_s^3\ll1
The density n is number per volume. The coupling g has units energy times volume. The dilute, repulsive, homogeneous three-dimensional gas is the model under discussion.↔ Scroll wide equations horizontally.

Pairing fluctuations produces collective sound

Write the field as a macroscopic condensate contribution plus a fluctuation, and retain the Hamiltonian through quadratic order in nonzero-momentum fluctuations. Interactions couple creation of momentum k\mathbf k to creation of −k-\mathbf k, so ordinary single-mode occupation operators no longer diagonalize the quadratic expression. A bosonic Bogoliubov transformation mixes creation and annihilation operators while preserving their commutator. Keeping that algebra is essential: a fermionic rotation would give the wrong normalization.

Diagonalization yields a positive excitation energy for repulsive coupling. At small momentum, the free kinetic energy is much smaller than the interaction energy, and the spectrum becomes linear with sound speed c=gn0/mc=\sqrt{gn_0/m}. At large momentum within the effective model, it approaches a particle-like branch with an interaction shift. The crossover is smooth. Quadratic theory neglects interactions among these quasiparticles, so it does not by itself predict their finite lifetime or the whole finite-temperature phase diagram.

ϵk=ℏ2k22m,Ek=ϵk(ϵk+2gn0),Ek∼k→0ℏkgn0m\epsilon_k=\frac{\hbar^2k^2}{2m},\qquad E_k=\sqrt{\epsilon_k(\epsilon_k+2gn_0)},\qquad E_k\underset{k\to0}{\sim}\hbar k\sqrt{\frac{gn_0}{m}}
n₀ is condensate density. At leading order in small depletion n₀ can be replaced by n, but they are conceptually different quantities.↔ Scroll wide equations horizontally.

Coherence is related to, but different from, flow

Condensate fraction measures the macroscopic eigenvalue of a one-body density matrix. Superfluid density measures response to a slow phase twist or moving boundary. They are not interchangeable definitions. Interactions can deplete the condensate while leaving a homogeneous Galilean-invariant ground state fully superfluid at zero temperature. An ideal Bose condensate, conversely, does not acquire a finite Landau critical velocity from a quadratic free-particle dispersion. Response requires its own calculation.

Dimension and order of limits also matter. A homogeneous interacting one-dimensional gas generally lacks conventional long-range condensate order in its thermodynamic ground state, and a uniform two-dimensional gas at positive temperature requires a different treatment of long-wavelength phase fluctuations. Finite traps and finite observation lengths modify these statements. Before using a condensate ansatz, specify geometry, temperature, interaction scale, and the observable being predicted. A useful mean field begins a controlled explanation; it does not erase fluctuations.

Locate the excitation crossover

Use dimensionless energy units with gn₀ = 1. Compare modes whose free kinetic energies are ε = 0.02 and ε = 8.

  1. For the first mode, E=0.02(0.02+2)=0.0404≈0.2010E=\sqrt{0.02(0.02+2)}=\sqrt{0.0404}\approx0.2010.

  2. The sound approximation gives E≈2ϵ=0.2000E\approx\sqrt{2\epsilon}=0.2000, about half a percent below the quadratic-theory result.

  3. For the second mode, E=8(8+2)=80≈8.9443E=\sqrt{8(8+2)}=\sqrt{80}\approx8.9443. The high-energy approximation ϵ+gn0=9\epsilon+gn_0=9 is close; the sound formula gives 4 and is inappropriate.

The result

The same dispersion contains two regimes. Its approximations must be selected by ε/(gn₀), while the validity of Bogoliubov theory itself additionally requires weak depletion.

Explore the model

Interactive lab

Interactive lab · waves

Inspect phase coherence in a simpler wave model

Vary the phase difference and amplitude ratio of sin(2x−t)+a sin(2x−t+phase). The exact sum displays coherent addition at a common wave number. These are classical waves; the lab does not calculate condensate depletion, a Bogoliubov spectrum, or superfluid density.

Inspect phase coherence in a simpler wave modelTwo waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Current readings: Sum at x = 0: 0.866; Resultant amplitude: 1.732; Phase difference: 60°. All parameters are adjustable with the labeled controls below.−3.14−1.5701.573.14−3.4−1.701.73.4x · radiansAmplitude
SumFirst componentSecond component
Sum at x = 0
0.866
Resultant amplitude
1.732
Phase difference
60°
t = 0 s

Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forRelative phase affects coherent sums. Establishing a quantum condensate or superfluid response needs the density-matrix and response calculations described above.

Check your understanding

For a dilute repulsive gas, what does the small-k Bogoliubov spectrum imply?

Try it yourself

At fixed mass and interaction coupling, increase condensate density by a factor of four. How does the leading sound speed change? State one condition you must recheck.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Use the square root in c and the independent diluteness condition.

Compare with the solution

The sound speed doubles because c is proportional to √n₀. The gas parameter naₛ³ increases with density, so weak depletion and the dilute-gas approximation must be checked again. A larger sound speed is not permission to extrapolate the leading formula into a strongly interacting regime.

Take it with you

  • Condensation is macroscopic occupation of a one-body mode.
  • Weak interactions turn long-wavelength fluctuations into sound-like quasiparticles.
  • Condensate fraction and superfluid response answer different questions.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.