2.4.18Thermodynamics & Statistical Mechanics (Advanced)

Bose-Einstein condensation — concept

1,876 words9 min readdifficulty · medium2 backlinks

WHY does BEC happen at all?

WHAT makes bosons special? They obey Bose–Einstein statistics: any number of bosons may occupy the same single-particle state. The mean occupation of a state with energy ε\varepsilon is

nˉ(ε)=1e(εμ)/kBT1.\bar n(\varepsilon) = \frac{1}{e^{(\varepsilon-\mu)/k_BT}-1}.

HOW do we constrain μ\mu? Occupations must be non-negative, so e(εμ)/kBT>1e^{(\varepsilon-\mu)/k_BT} > 1 for all states. The lowest energy is ε0=0\varepsilon_0 = 0, so we need

μ0(for bosons, with ground state at 0).\mu \le 0 \quad\text{(for bosons, with ground state at }0).

The ground-state occupation is N0=1eμ/kBT1.N_0 = \frac{1}{e^{-\mu/k_BT}-1}. As μ0\mu \to 0^-, N0N_0 \to \infty. This is the escape valve: when the excited states saturate, μ\mu pins itself just below 00 and the ground state quietly soaks up the rest.


Deriving the saturation number (from scratch)

Step 1 — Count states. Sum over excited states → integral with density of states. Why this step? For a large box the levels are dense; replacing the sum by an integral is the continuum approximation. The catch (remember it!): the integral misses the ε=0\varepsilon=0 state because the density of states g(ε)εg(\varepsilon)\propto\sqrt\varepsilon vanishes there — so we must track N0N_0 separately.

The density of states (derived from counting kk-shells, N(k)=V(2π)343πk3gN(k)=\frac{V}{(2\pi)^3}\cdot\frac{4}{3}\pi k^3 \cdot g): g(ε)dε=gV4π2(2m2)3/2εdε.g(\varepsilon)\,d\varepsilon = \frac{gV}{4\pi^2}\left(\frac{2m}{\hbar^2}\right)^{3/2}\sqrt{\varepsilon}\,d\varepsilon.

Step 2 — Maximum excited population occurs at μ=0\mu = 0: Nexcmax(T)=0g(ε)dεeε/kBT1.N_{\text{exc}}^{\max}(T)=\int_0^\infty \frac{g(\varepsilon)\,d\varepsilon}{e^{\varepsilon/k_BT}-1}.

Step 3 — Substitute x=ε/kBTx=\varepsilon/k_BT. Why? To pull all TT-dependence out front and expose a pure number. Nexcmax=gV4π2(2mkBT2)3/20xdxex1Γ(3/2)ζ(3/2).N_{\text{exc}}^{\max}=\frac{gV}{4\pi^2}\left(\frac{2mk_BT}{\hbar^2}\right)^{3/2}\underbrace{\int_0^\infty\frac{\sqrt{x}\,dx}{e^x-1}}_{\Gamma(3/2)\,\zeta(3/2)}.

Using Γ(3/2)=π2\Gamma(3/2)=\tfrac{\sqrt\pi}{2} and the thermal de Broglie wavelength λT=h2πmkBT\lambda_T=\dfrac{h}{\sqrt{2\pi m k_BT}}, this collapses beautifully to

WHAT it means: the excited states can hold at most gζ(3/2)\sim g\,\zeta(3/2) particles per thermal volume λT3\lambda_T^3. Push the real density higher (or cool TT so λT\lambda_T grows) and the rest condense.


The transition temperature

Set nexcmax(Tc)=nn_{\text{exc}}^{\max}(T_c)=n (the actual density) — the point where excited states are just full: n=gζ(3/2)λTc3    λTc3=gζ(3/2)n.n=g\frac{\zeta(3/2)}{\lambda_{T_c}^3}\;\Rightarrow\; \lambda_{T_c}^3=\frac{g\,\zeta(3/2)}{n}.

Solving for TT (using λTT1/2\lambda_T\propto T^{-1/2}):


Condensate fraction below TcT_c

For T<TcT<T_c, μ0\mu\approx 0, so Nexc(T)T3/2N_{\text{exc}}(T)\propto T^{3/2}. Compare to its value at TcT_c (where Nexc=NN_{\text{exc}}=N): Nexc(T)N=(TTc)3/2.\frac{N_{\text{exc}}(T)}{N}=\left(\frac{T}{T_c}\right)^{3/2}. Everything else is in the ground state:

Figure — Bose-Einstein condensation — concept

Worked examples


Common mistakes (Steel-manned)


Active recall

Recall Forecast-then-verify: predict before unfolding

Q: As you add atoms at fixed TT, what happens to μ\mu and where do extra atoms go past saturation? A: μ\mu rises toward 00 then sticks there; surplus atoms pour into the ε=0\varepsilon=0 ground state (the condensate). N0N_0 diverges as μ0\mu\to0^-, providing infinite capacity.

Recall Feynman: explain to a 12-year-old

Imagine a stadium where the cheap seats (excited states) can only hold a fixed crowd at a given "energy budget" (temperature). Identical, super-social fans (bosons) are happy to all sit on the same special VIP chair on the field (the ground state). When the cheap seats fill up and more fans arrive, they have only one place left: they ALL pile onto that one VIP chair. Suddenly a huge fraction of fans share a single seat — that's a Bose–Einstein condensate. It happens not because they shove each other, but because the rules let them share, and the cheap seats simply run out of room when it's cold.

What statistics do bosons obey for state occupation?
nˉ(ε)=1e(εμ)/kBT1\bar n(\varepsilon)=\frac{1}{e^{(\varepsilon-\mu)/k_BT}-1}
Why must μ0\mu\le 0 for an ideal Bose gas (ground state at 0)?
Else nˉ(ε0)\bar n(\varepsilon_0) would be negative; occupations must be non-negative.
What physically is the BEC condensate?
A macroscopic fraction of bosons occupying the single lowest-energy quantum state.
The phase-space density condition for BEC onset?
nλT3=gζ(3/2)2.612n\lambda_T^3 = g\,\zeta(3/2)\approx 2.612.
Formula for the critical temperature TcT_c?
kBTc=2π2m(ngζ(3/2))2/3k_BT_c=\frac{2\pi\hbar^2}{m}\left(\frac{n}{g\zeta(3/2)}\right)^{2/3}
Condensate fraction below TcT_c?
N0/N=1(T/Tc)3/2N_0/N = 1-(T/T_c)^{3/2}.
Why must the ground state be counted separately in the integral?
g(ε)ε0g(\varepsilon)\propto\sqrt\varepsilon\to0 at ε=0\varepsilon=0, so the integral misses the macroscopically occupied ground state.
Is BEC driven by interactions?
No — it is purely a quantum-statistical effect; ideal (non-interacting) bosons condense.
Why is there no BEC in a uniform 2D ideal Bose gas?
The saturation integral diverges, so excited states never saturate.
Value of ζ(3/2)\zeta(3/2)?
2.612\approx 2.612.
Condensate fraction at T=0.5TcT=0.5T_c?
1(0.5)3/20.6461-(0.5)^{3/2}\approx 0.646 (about 65%).

Connections

  • Bose-Einstein statistics — the occupation law that makes condensation possible.
  • Thermal de Broglie wavelengthλT=h/2πmkBT\lambda_T=h/\sqrt{2\pi mk_BT} sets the overlap criterion.
  • Density of states — its ε\sqrt\varepsilon form causes the ground-state subtlety.
  • Chemical potential — pinning at 00^- is the mechanism of condensation.
  • Riemann zeta functionζ(3/2)\zeta(3/2) appears from the saturation integral.
  • Phase transitions — BEC as a continuous transition with no latent heat from interactions.
  • Superfluidity — BEC underlies superfluid 4^4He and atomic condensates.
  • Fermi-Dirac statistics — contrast: Pauli exclusion forbids condensation for fermions.

Concept Map

allows

gives

non-negative requires

as mu to 0

excited states

misses eps=0 so

integrate at mu=0

sub x=eps/kT

with lambda_T

extras overflow into

escape valve for

holds

Bose-Einstein statistics

Any number share a state

Mean occupation n-bar

Chemical potential mu <= 0

Ground-state N0 diverges

Density of states sqrt eps

Track N0 separately

Max excited population

Yields zeta 3/2 factor

Saturation density

Bose-Einstein condensate

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, Bose–Einstein condensation ka core idea bahut simple hai. Boson particles (jaise photons, ya cold atoms like 87^{87}Rb) ka ek special property hai — woh ek hi quantum state mein jitne chaaho utne baith sakte hain, koi limit nahi (fermions ke ulta). Jab tum gas ko bahut zyada thanda karte ho, har atom ki "matter wave" (de Broglie wavelength λT\lambda_T) badi hoti jaati hai. Ek time aata hai jab ye waves padosi atoms ki waves se overlap karne lagti hain. Bas wahi moment BEC ka onset hai.

Ab WHY important hai: excited states (energy >0>0) ek fixed maximum number of particles hi rakh sakti hain at a given temperature — yeh hum derive karte hain aur answer milta hai nλT3=2.612n\lambda_T^3 = 2.612. Agar isse zyada atoms ghusao, toh extra atoms ke paas jaane ki ekmatra jagah bachti hai — ground state (ε=0\varepsilon=0). Sabhi surplus atoms ek hi lowest state mein dhad-dhad gir jaate hain. Yahi macroscopic pile-up hi condensate hai.

Critical temperature ka formula yaad rakho: kBTc=2π2m(n/gζ(3/2))2/3k_BT_c=\frac{2\pi\hbar^2}{m}(n/g\zeta(3/2))^{2/3}, aur condensate fraction N0/N=1(T/Tc)3/2N_0/N = 1-(T/T_c)^{3/2}. Real alkali gases mein density chhoti hoti hai isliye TcT_c sirf hundreds of nanokelvin hota hai — isliye Nobel-winning lab cooling chahiye. Ek important baat: BEC interactions ki wajah se nahi hota, yeh pure quantum statistics ka khel hai. Aur ek galti mat karna — integral mein ground state alag se count karna padta hai, kyunki density of states ε\sqrt\varepsilon ground state ko kha jaata hai.

Exam ke liye 80/20: phase-space density condition nλT32.612n\lambda_T^3\approx 2.612 yaad rakho, TcT_c formula, aur condensate fraction — bas yahi sab kuch cover kar deta hai.

Go deeper — visual, from zero

Test yourself — Thermodynamics & Statistical Mechanics (Advanced)

Connections