4.2 Bose-Einstein statistics

Bose-Einstein statistics describe how particles with integer spin behave in quantum systems. This framework explains phenomena like condensation and superfluidity, where many particles occupy the same quantum state at low temperatures. It underpins key concepts in condensed matter physics and quantum optics, revealing how bosons like photons and certain atoms behave collectively.

Fundamentals of Bose-Einstein statistics

Bose-Einstein statistics govern the thermal equilibrium behavior of bosons, particles that don't obey the Pauli exclusion principle. Unlike classical Maxwell-Boltzmann statistics, which treat particles as distinguishable, this framework accounts for the fundamental indistinguishability of quantum particles and predicts striking collective phenomena at low temperatures.

Bosons and indistinguishability

Bosons are particles with integer spin (0, 1, 2, ...) whose wavefunctions are symmetric under particle exchange. This symmetry has a dramatic consequence: any number of bosons can pile into the same quantum state simultaneously. Common examples include photons (spin 1), gluons (spin 1), W/Z bosons (spin 1), and composite particles like $^4\text{He}$ atoms, where the total spin adds up to an integer.

This unlimited state occupancy is what distinguishes bosons from fermions and drives all the collective behavior covered in this topic.

Quantum statistics principles

At high temperatures and low densities, bosons behave much like classical particles, and Maxwell-Boltzmann statistics work fine. The quantum nature becomes important when the thermal de Broglie wavelength of particles becomes comparable to the inter-particle spacing. At that point:

• Particle indistinguishability can no longer be ignored
• The wave-like nature of particles leads to overlapping wavefunctions
• The system begins to "feel" quantum statistics, and the occupation of low-energy states grows dramatically
• At sufficiently low temperatures, a macroscopic fraction of particles condenses into the lowest energy state

Bose-Einstein distribution function

The average number of particles occupying a single-particle energy state $E_i$ is given by the Bose-Einstein distribution:

$n_i = \frac{1}{e^{(E_i - \mu)/k_BT} - 1}$

where:

• $n_i$ is the mean occupation number of state $i$
• $E_i$ is the energy of state $i$
• $\mu$ is the chemical potential
• $k_B$ is the Boltzmann constant
• $T$ is the temperature

Notice the $-1$ in the denominator (contrast this with $+1$ for Fermi-Dirac). This minus sign allows occupation numbers to exceed 1 and, in fact, to diverge. As $\mu$ approaches the ground state energy from below, the occupation of that state grows without bound, signaling the onset of Bose-Einstein condensation.

Bose-Einstein condensation

Bose-Einstein condensation (BEC) occurs when a macroscopic number of bosons collapse into the single-particle ground state. Below a critical temperature, the excited states simply can't accommodate all the particles, so the excess "condenses" into the lowest energy level. The result is a coherent matter wave with long-range order.

Critical temperature

The critical temperature $T_c$ marks the onset of condensation. For an ideal (non-interacting) Bose gas in three dimensions:

$T_c = \frac{2\pi\hbar^2}{mk_B}\left(\frac{n}{\zeta(3/2)}\right)^{2/3}$

• $n$ is the particle number density
• $m$ is the particle mass
• $\zeta(3/2) \approx 2.612$ is the Riemann zeta function evaluated at 3/2

Higher density and lower mass both push $T_c$ upward. Above $T_c$, the system behaves as a normal (though quantum-corrected) gas. Below $T_c$, a growing fraction of particles occupies the ground state.

Condensate fraction

Below $T_c$, the fraction of particles in the ground state is:

$\frac{N_0}{N} = 1 - \left(\frac{T}{T_c}\right)^{3/2} \quad \text{for } T < T_c$

where $N_0$ is the number of condensed particles and $N$ is the total particle number. At $T = 0$, the condensate fraction reaches unity for an ideal Bose gas (all particles in the ground state). In real interacting systems, quantum depletion reduces this even at $T = 0$.

Superfluid properties

Bose-Einstein condensates exhibit superfluidity, characterized by:

• Zero viscosity: the fluid flows without dissipation
• Infinite thermal conductivity
• Quantized vortices: circulation is restricted to integer multiples of $h/m$
• Persistent currents that flow indefinitely

The two-fluid model describes the system as a mixture of a normal component (thermally excited particles) and a superfluid component (the condensate). This model explains phenomena like the fountain effect and second sound (a temperature wave) observed in liquid helium-4.

Applications of Bose-Einstein statistics

Photon gas

Photons are massless spin-1 bosons, and blackbody radiation is a photon gas in thermal equilibrium. Because photons can be freely created and destroyed, their chemical potential is zero ($\mu = 0$). Plugging this into the Bose-Einstein distribution gives Planck's radiation law:

$\langle n(\nu) \rangle = \frac{1}{e^{h\nu/k_BT} - 1}$

This distribution explains the spectrum of blackbody radiation, the physics of stimulated emission in lasers, and the cosmic microwave background radiation (a near-perfect blackbody at 2.725 K).

Phonons in solids

Phonons are quantized lattice vibrations in crystalline solids. Like photons, they are bosons with $\mu = 0$ (their number isn't conserved). Applying Bose-Einstein statistics to phonons yields the Debye model of specific heat, which correctly predicts:

• The $T^3$ dependence of specific heat at low temperatures
• Deviations from the classical Dulong-Petit law ($C_V = 3Nk_B$) that appears at high temperatures
• Thermal conductivity and sound propagation in materials

Phonons also play a central role in conventional superconductivity through electron-phonon coupling (Cooper pairing).

Liquid helium-4

$^4\text{He}$ undergoes a phase transition at 2.17 K, known as the lambda point (named for the lambda-shaped specific heat curve). Below this temperature, it becomes a superfluid (He-II) with zero viscosity and the ability to support quantized vortices and second sound. While the lambda transition is related to BEC, strong inter-particle interactions in liquid helium mean it's not a simple ideal-gas condensation. Only about 8% of atoms occupy the ground state even at $T = 0$.

Thermodynamic properties

Specific heat

The specific heat $C_V = \left(\frac{\partial U}{\partial T}\right)_V$ of a Bose gas shows distinctive behavior near the condensation transition. For an ideal Bose gas in 3D, $C_V$ has a cusp (discontinuity in its derivative) at $T_c$. In real systems like liquid helium-4, this becomes the famous lambda anomaly, a sharp, divergent peak at the superfluid transition. At low temperatures in solids, Bose-Einstein statistics for phonons explain why $C_V$ drops well below the classical Dulong-Petit value.

Entropy and free energy

The Helmholtz free energy is $F = U - TS$, and the equilibrium state minimizes $F$ at constant $T$ and $V$. As temperature drops toward absolute zero, the entropy of a Bose gas decreases rapidly, consistent with the third law of thermodynamics. (Classical statistics would incorrectly predict finite entropy at $T = 0$.) Near phase transitions, the free energy and its derivatives can exhibit discontinuities or singularities that characterize the order of the transition.

Chemical potential

The chemical potential $\mu$ controls the average particle number through:

$N = \sum_i \frac{1}{e^{(E_i - \mu)/k_BT} - 1}$

For a Bose gas, $\mu$ must always be less than or equal to the ground state energy $E_0$ (otherwise occupation numbers would go negative, which is unphysical). As $T$ decreases toward $T_c$, $\mu$ rises toward $E_0$ from below. At and below $T_c$, $\mu$ pins to $E_0$, and the ground state acquires macroscopic occupation. This behavior of $\mu$ is the mathematical signature of condensation.

Quantum gases

Ideal Bose gas

The ideal Bose gas is the simplest model: non-interacting bosons in a box (or harmonic trap). Despite its simplicity, it captures the essential physics of BEC, including the critical temperature, condensate fraction scaling, and the specific heat anomaly. It serves as the starting point for perturbative treatments of interactions.

Weakly interacting Bose gas

Real atomic gases have interactions, even if weak. In the mean-field approximation, the condensate wavefunction $\Psi(\mathbf{r}, t)$ obeys the Gross-Pitaevskii equation:

$i\hbar \frac{\partial \Psi}{\partial t} = \left(-\frac{\hbar^2}{2m}\nabla^2 + V_{\text{ext}} + g|\Psi|^2\right)\Psi$

where $g$ characterizes the interaction strength. Interactions modify the condensation temperature, change the excitation spectrum (giving rise to Bogoliubov quasiparticles with a linear, phonon-like dispersion at low momenta), and are essential for understanding superfluidity in dilute gases.

Bose-Einstein vs Fermi-Dirac statistics

Both reduce to Maxwell-Boltzmann statistics in the high-temperature, low-density (classical) limit.

Experimental observations

Cold atom experiments

The first gaseous BECs were created in 1995 by Eric Cornell and Carl Wieman (with $^{87}\text{Rb}$) and independently by Wolfgang Ketterle (with $^{23}\text{Na}$), earning them the 2001 Nobel Prize. These experiments use laser cooling followed by evaporative cooling in magnetic or optical traps to reach temperatures in the nanokelvin range.

Key capabilities of cold atom experiments:

• Tunable inter-particle interactions via Feshbach resonances (using external magnetic fields)
• Creation of optical lattices for simulating condensed matter systems
• Direct imaging of condensate wavefunctions and momentum distributions
• Serving as quantum simulators for complex many-body Hamiltonians

Laser cooling techniques

Reaching BEC requires cooling atoms through several stages:

1. Doppler cooling: Photon momentum kicks from red-detuned laser beams slow atoms to millikelvin temperatures
2. Sub-Doppler cooling (e.g., Sisyphus cooling): Exploits optical pumping in polarization gradients to reach microkelvin temperatures
3. Evaporative cooling: Selectively removes the most energetic atoms from a trap, allowing the remaining gas to rethermalize at a lower temperature. This final stage typically reaches the nanokelvin regime needed for condensation

Detection methods

• Time-of-flight imaging: Turn off the trap and let the cloud expand freely. The resulting density profile maps the original momentum distribution. A BEC appears as a sharp peak at zero momentum.
• In-situ imaging: Measures the real-space density profile while the atoms are still trapped
• Bragg spectroscopy: Uses two-photon transitions to probe the excitation spectrum and coherence properties
• Interference experiments: Overlapping two condensates produces interference fringes, demonstrating long-range phase coherence
• Noise correlation measurements: Reveal quantum statistical effects (bunching for bosons) in the density fluctuations

Mathematical formalism

Partition function

The grand canonical partition function for a Bose gas factorizes over single-particle states:

$\ln \Xi = -\sum_i \ln\left(1 - e^{-(E_i - \mu)/k_BT}\right)$

All thermodynamic quantities follow from this. The free energy is $F = -k_BT \ln \Xi$, and derivatives give the entropy, pressure, and particle number. The key step in deriving the Bose-Einstein distribution is maximizing the entropy subject to constraints on total energy and particle number.

Grand canonical ensemble

The grand canonical ensemble is the natural framework for BEC because particle number fluctuates between the condensate and excited states. The grand partition function is:

$\Xi = \sum_N Z_N \, e^{\mu N / k_BT}$

where $Z_N$ is the canonical partition function for $N$ particles. The chemical potential $\mu$ acts as a Lagrange multiplier controlling the average particle number. This formalism treats BEC as a genuine phase transition and simplifies calculations considerably for the ideal Bose gas.

Occupation number statistics

For a single energy level $E_i$, the probability of finding $n$ bosons in that state is:

$P_i(n) = \left(1 - e^{-(E_i - \mu)/k_BT}\right) e^{-n(E_i - \mu)/k_BT}$

This is a geometric distribution, meaning the most probable occupation is always $n = 0$, but the tail extends to arbitrarily large $n$. The variance in occupation number satisfies $\langle (\Delta n_i)^2 \rangle = \langle n_i \rangle (\langle n_i \rangle + 1)$, which exceeds the Poisson value $\langle n_i \rangle$. This super-Poissonian fluctuation is the statistical origin of photon bunching (the Hanbury Brown-Twiss effect) and is a direct signature of Bose-Einstein statistics.

Bose-Einstein in low dimensions

2D and 1D systems

Dimensionality profoundly affects condensation. The Mermin-Wagner theorem states that true long-range order (and therefore true BEC) cannot exist in 2D at any finite temperature for systems with continuous symmetry. Instead, 2D Bose systems undergo a Berezinskii-Kosterlitz-Thouless (BKT) transition, where the system develops quasi-long-range order through the binding and unbinding of vortex-antivortex pairs.

In 1D, quantum fluctuations are even more dominant. The Tonks-Girardeau gas is a strongly interacting 1D Bose gas that, remarkably, maps onto a non-interacting Fermi gas in its density correlations (though not in momentum space). These reduced-dimensional systems are realized experimentally using tightly confining optical traps and atom chips.

Quantum phase transitions

Quantum phase transitions occur at $T = 0$ and are driven by varying a system parameter (like interaction strength or lattice depth) rather than temperature. The most celebrated example in bosonic systems is the superfluid-to-Mott-insulator transition in an optical lattice:

• At shallow lattice depth, tunneling dominates and the system is a superfluid with delocalized atoms
• At deep lattice depth, interactions dominate and atoms localize with integer filling per site (Mott insulator)
• Near the critical point, the system exhibits universal scaling behavior

These transitions are governed by quantum fluctuations and provide a testing ground for quantum many-body theory.

Topological effects

Low-dimensional Bose systems with specific interactions or geometries can host topological phases:

• Bosonic analogues of the quantum Hall effect
• Exotic quasiparticles like anyons with fractional statistics (neither bosonic nor fermionic)
• Topologically protected states that are robust against local perturbations, making them candidates for quantum computing

These phases are realized in synthetic quantum systems, such as cold atoms subjected to artificial gauge fields created by laser-assisted tunneling or rotation.

Advanced topics

Non-equilibrium dynamics

Much of BEC physics assumes thermal equilibrium, but some of the richest phenomena occur out of equilibrium:

• Quantum quenches: Suddenly changing a system parameter and watching the subsequent dynamics
• Prethermalization: The system reaches a quasi-steady state before eventually thermalizing (if it thermalizes at all)
• Quantum turbulence: Tangled networks of quantized vortices, the quantum analogue of classical turbulence
• Dynamical phase transitions: Non-analytic behavior in time-evolved observables

These studies address fundamental questions about how isolated quantum systems approach thermal equilibrium.

Bose-Einstein in optical lattices

Optical lattices are created by interfering laser beams to form a periodic potential for cold atoms. Loading a BEC into an optical lattice realizes the Bose-Hubbard model:

$H = -J\sum_{\langle i,j \rangle} a_i^\dagger a_j + \frac{U}{2}\sum_i n_i(n_i - 1)$

where $J$ is the tunneling amplitude and $U$ is the on-site interaction energy. By tuning the lattice depth, you control $J/U$ and can drive the superfluid-to-Mott-insulator transition. Optical lattices also enable the study of many-body localization, disorder effects, and frustrated magnetism.

Quantum simulation applications

BECs in engineered potentials serve as quantum simulators, physical systems that mimic the behavior of other quantum systems too complex for classical computation. Active areas include:

• Simulating condensed matter models (Hubbard models, spin chains, gauge theories)
• Analogue gravity: BEC phonons propagating in a flowing condensate experience an effective spacetime metric, enabling laboratory studies of Hawking radiation analogues
• Exploring many-body entanglement and quantum information protocols

These platforms bridge atomic physics, condensed matter, and high-energy physics, offering experimental access to problems that remain theoretically intractable.