Skip to content

Bose-Einstein Condensation

For bosons, the average occupation of a single-particle state of energy ε\varepsilon is

nε=1eβ(εμ)1\langle n_\varepsilon \rangle = \frac{1}{e^{\beta(\varepsilon - \mu)} - 1}

The chemical potential must satisfy με0\mu \leq \varepsilon_0 (the lowest single-particle energy) to prevent negative occupation numbers.

5.2 Density of States and Critical Temperature

Section titled “5.2 Density of States and Critical Temperature”

For a 3D free Bose gas with ε=2k2/(2m)\varepsilon = \hbar^2 k^2 / (2m)The density of states is g(ε)=(V/4π2)(2m/2)3/2εg(\varepsilon) = (V/4\pi^2)(2m/\hbar^2)^{3/2}\sqrt{\varepsilon}. The number of particles in excited states (ε>0\varepsilon > 0) is

Nex=0g(ε)dεeβε1=V(mkBT2π2)3/2ζ ⁣(32)N_{\mathrm{ex} = \int_0^\infty \frac{g(\varepsilon)\, d\varepsilon}{e^{\beta \varepsilon} - 1} = V\left(\frac{mk_BT}{2\pi\hbar^2}\right)^{3/2}\,\zeta\!\left(\frac{3}{2}\right)}

Where ζ(3/2)2.612\zeta(3/2) \approx 2.612 is the Riemann zeta function.

Theorem 5.1 (BEC critical temperature). The maximum number of particles that can be accommodated in excited states is achieved at μ=0\mu = 0. When NN exceeds this maximum, the excess condenses into the ground state. The critical temperature is

Tc=2π2mkB(nζ(3/2))2/3T_c = \frac{2\pi\hbar^2}{mk_B}\left(\frac{n}{\zeta(3/2)}\right)^{2/3}

Where n=N/Vn = N/V.

Proof. Setting N=NexmaxN = N_{\mathrm{ex}^{\max}} at μ=0\mu = 0 and solving for TT:

n=(mkBTc2π2)3/2ζ ⁣(32)n = \left(\frac{mk_B T_c}{2\pi\hbar^2}\right)^{3/2}\,\zeta\!\left(\frac{3}{2}\right)

Tc=2π2mkB(nζ(3/2))2/3T_c = \frac{2\pi\hbar^2}{mk_B}\left(\frac{n}{\zeta(3/2)}\right)^{2/3} \qquad \blacksquare

Below TcT_c, μ0\mu \approx 0 and the condensate fraction is

N0N=1(TTc)3/2\frac{N_0}{N} = 1 - \left(\frac{T}{T_c}\right)^{3/2}

This follows from N0=NNexN_0 = N - N_{\mathrm{ex}} with μ=0\mu = 0:

Nex=N(TTc)3/2N_{\mathrm{ex} = N\left(\frac{T}{T_c}\right)^{3/2}}

5.4 Thermodynamic Properties below TcT_c

Section titled “5.4 Thermodynamic Properties below TcT_cTc​”

The energy below TcT_c:

U=0εg(ε)dεeβε1=V(mkBT2π2)3/2(kBT)32ζ ⁣(52)Γ ⁣(52)U = \int_0^\infty \frac{\varepsilon\, g(\varepsilon)\, d\varepsilon}{e^{\beta\varepsilon} - 1} = V\left(\frac{mk_BT}{2\pi\hbar^2}\right)^{3/2}\,(k_BT)\,\frac{3}{2}\,\zeta\!\left(\frac{5}{2}\right) \cdot \Gamma\!\left(\frac{5}{2}\right)

=32NkBTcζ ⁣(52)/ζ ⁣(32)(TTc)5/2= \frac{3}{2}\,Nk_BT_c\,\zeta\!\left(\frac{5}{2}\right)\Big/\zeta\!\left(\frac{3}{2}\right)\,\left(\frac{T}{T_c}\right)^{5/2}

The heat capacity:

CV=154NkBζ ⁣(52)/ζ ⁣(32)(TTc)3/2T3/2C_V = \frac{15}{4}\,Nk_B\,\zeta\!\left(\frac{5}{2}\right)\Big/\zeta\!\left(\frac{3}{2}\right)\,\left(\frac{T}{T_c}\right)^{3/2} \propto T^{3/2}

This contrasts with the constant CV=32NkBC_V = \frac{3}{2}Nk_B above TcT_c (equipartition). There is a cusp (discontinuity in the derivative) at TcT_cCharacteristic of a phase transition.

Problem. Estimate TcT_c for a gas of N=104N = 10^4 rubidium-87 atoms confined in a harmonic trap with frequency ωho=2π×100\omega_{\mathrm{ho} = 2\pi \times 100} Hz.

Solution

For a harmonic trap, the effective density of states is g(ε)=ε2/(23ωho3)g(\varepsilon) = \varepsilon^2/(2\hbar^3\omega_{\mathrm{ho}^3)}. The critical temperature in a harmonic trap is:

kBTc=ωho(Nζ(3))1/3k_BT_c = \hbar\omega_{\mathrm{ho}\left(\frac{N}{\zeta(3)}\right)^{1/3}}

kBTc=(1.055×1034)(2π×100)(1041.202)1/3k_BT_c = (1.055 \times 10^{-34})(2\pi \times 100)\left(\frac{10^4}{1.202}\right)^{1/3}

=(6.63×1032)(20.1)=1.33×1030J= (6.63 \times 10^{-32})(20.1) = 1.33 \times 10^{-30}\,\mathrm{J}

Tc=1.33×10301.381×10239.6×108K96nKT_c = \frac{1.33 \times 10^{-30}}{1.381 \times 10^{-23}} \approx 9.6 \times 10^{-8}\,\mathrm{K} \approx 96\,\mathrm{nK}

This is consistent with the 1995 BEC experiments by Cornell and Wieman (JILA) and Ketterle (MIT), who achieved BEC at temperatures of a few hundred nanokelvin. \blacksquare


QuantityExpressionPhysical Meaning
Critical temperatureTc=2π2mkB(nζ(3/2))2/3T_c = \frac{2\pi\hbar^2}{mk_B}\left(\frac{n}{\zeta(3/2)}\right)^{2/3}Onset of macroscopic occupation
Condensate fractionN0N=1(T/Tc)3/2\frac{N_0}{N} = 1 - (T/T_c)^{3/2}Order parameter below TcT_c
Energy below TcT_cUNkBTc(T/Tc)5/2U \propto Nk_B T_c (T/T_c)^{5/2}Deviates from equipartition
Heat capacityCVT3/2C_V \propto T^{3/2} below TcT_cSignature of BEC phase
de Broglie wavelengthλth=h/2πmkBT\lambda_{\text{th}} = h/\sqrt{2\pi m k_B T}BEC occurs when nλth32.612n\lambda_{\text{th}}^3 \approx 2.612
  1. BEC is not a classical condensation: BEC is a purely quantum phenomenon driven by Bose statistics, not by interparticle interactions. An ideal Bose gas condenses, whereas a classical gas would not.
  2. Finite-size effects: The critical temperature derived assumes the thermodynamic limit (NN \to \infty, VV \to \infty, nn fixed). For finite traps with N104N \sim 10^4, there are corrections of order N1/3N^{-1/3}.
  3. Dimensionality matters: In 2D, the density of states is constant and the integral for NexN_{\text{ex}} diverges at μ=0\mu = 0 only logarithmically. Strict BEC does not occur in 2D uniform gases (Mermin—Wagner—Hohenberg theorem).
  4. Interactions modify TcT_c: Repulsive interactions slightly suppress TcT_c relative to the ideal gas prediction. The shift is ΔTc/Tc(n1/3as)\Delta T_c/T_c \propto (n^{1/3}a_s), where asa_s is the scattering length.
  • Atom lasers: A BEC releases coherent matter waves, analogous to an optical laser. Coherence lengths exceeding 1 mm have been demonstrated.
  • Precision measurement: BEC interferometry measures gravitational acceleration, rotations, and fundamental constants with extreme sensitivity.
  • Superfluid helium: Liquid 4^4He below 2.17 K exhibits superfluidity, with approximately 10% of atoms in the condensate (strongly interacting, unlike the ideal gas model).
  • Quantum simulation: Optical lattices loaded with BEC simulate the Hubbard model, enabling studies of quantum phase transitions.
  • Slow light: Electromagnetically induced transparency in BEC reduces light speed to metres per second.
  • Superconductivity: The BCS ground state is a condensate of Cooper pairs (composite bosons). The BCS—BEC crossover connects fermionic pairing to molecular BEC.
  • Quantum field theory: BEC is an example of spontaneous symmetry breaking — the U(1)U(1) phase symmetry of the matter field is broken, giving rise to a Goldstone mode (Bogoliubov phonon).
  • Statistical mechanics: The BEC transition is a textbook example of a phase transition driven purely by statistics, requiring no interactions.

Summary Table: Ideal Bose Gas vs Ideal Fermi Gas

Section titled “Summary Table: Ideal Bose Gas vs Ideal Fermi Gas”
PropertyBose GasFermi Gas
Statisticsni=(eβ(ϵiμ)1)1\langle n_i \rangle = (e^{\beta(\epsilon_i-\mu)} - 1)^{-1}ni=(eβ(ϵiμ)+1)1\langle n_i \rangle = (e^{\beta(\epsilon_i-\mu)} + 1)^{-1}
μ\mu constraintμϵ0\mu \leq \epsilon_0μ\mu unrestricted (can be positive at T=0T=0)
T=0T=0 stateAll particles in ground stateFilled up to ϵF\epsilon_F
Low-TT heat capacityCVT3/2C_V \propto T^{3/2}CVTC_V \propto T
Phase transitionBEC at TcT_cNo phase transition
High-TT limitMaxwell—BoltzmannMaxwell—Boltzmann