4.4 Ideal quantum gases

Fundamentals of ideal quantum gases

Ideal quantum gases are systems of many non-interacting particles that obey quantum mechanics rather than classical mechanics. They serve as the starting models for understanding real quantum many-body systems, much like the classical ideal gas serves as a starting point in thermodynamics. The key departure from classical physics is that identical quantum particles are truly indistinguishable, and this single fact forces us into entirely different statistics.

Quantum statistical mechanics basics

Statistical mechanics for quantum systems replaces the classical phase space with discrete quantum states and energy levels. The Heisenberg uncertainty principle,

$\Delta x \, \Delta p \geq \frac{\hbar}{2}$

sets a fundamental limit on how finely you can specify a particle's position and momentum simultaneously. This means phase-space cells have a minimum volume of order $\hbar^3$ per degree of freedom.

When the thermal de Broglie wavelength of particles becomes comparable to the average interparticle spacing, the system enters the quantum degenerate regime. At this point, classical statistics breaks down and you must use quantum statistics. This typically happens at low temperatures or high densities.

Indistinguishability of particles

In quantum mechanics, two particles of the same species (e.g., two electrons) are fundamentally indistinguishable. Swapping them cannot produce a physically distinct state. This constraint forces the many-particle wavefunction to have definite symmetry under particle exchange:

• Symmetric wavefunctions → bosons
• Antisymmetric wavefunctions → fermions

This is not just a mathematical convenience. It has direct physical consequences: it determines which states are accessible and how particles distribute among energy levels. Classical particles, by contrast, are always treated as distinguishable, which is why classical counting (Maxwell-Boltzmann) overcounts states at low temperatures.

Bose-Einstein vs Fermi-Dirac statistics

The two types of quantum statistics follow directly from exchange symmetry:

• Bose-Einstein statistics apply to bosons (integer spin: 0, 1, 2, ...).
  • Any number of bosons can pile into the same quantum state.
  • Examples: photons, phonons, $^4\text{He}$ atoms, gluons.
• Fermi-Dirac statistics apply to fermions (half-integer spin: 1/2, 3/2, ...).
  • The Pauli exclusion principle limits each quantum state to at most one fermion.
  • Examples: electrons, protons, neutrons, $^3\text{He}$ atoms.

At high temperatures or low densities (where $\lambda_{dB} \ll$ interparticle spacing), both distributions converge to the classical Maxwell-Boltzmann distribution. The quantum corrections become negligible because the chance of two particles wanting the same state is vanishingly small.

Bose-Einstein condensation

Bose-Einstein condensation (BEC) occurs when a macroscopic fraction of bosons collapses into the single-particle ground state. It happens below a critical temperature and represents a genuine phase transition into a state of matter with macroscopic quantum coherence.

Bose-Einstein distribution function

The mean occupation number for bosons in a state with energy $\epsilon_i$ is:

$\langle n_i \rangle = \frac{1}{e^{(\epsilon_i - \mu)/k_B T} - 1}$

Here $\mu$ is the chemical potential, $k_B$ is Boltzmann's constant, and $T$ is temperature. Notice the $-1$ in the denominator. For this to stay non-negative, you need $\mu \leq \epsilon_0$ (the ground state energy). As $\mu \to \epsilon_0$, the ground state occupation diverges, which is exactly the onset of BEC.

Critical temperature

For a uniform 3D ideal Bose gas, the critical temperature is:

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

where $n$ is the particle number density, $m$ is the particle mass, and $\zeta(3/2) \approx 2.612$ is the Riemann zeta function evaluated at 3/2. Below $T_c$, the excited states can no longer accommodate all the particles, so the excess is forced into the ground state.

Condensate fraction

Below $T_c$, the fraction of particles in the ground state for a uniform 3D gas is:

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

This fraction grows smoothly from 0 at $T = T_c$ to 1 at $T = 0$. The derivative of the condensate fraction is discontinuous at $T_c$, which is the signature of a phase transition (specifically, a continuous transition with no latent heat).

Experimental realization

BEC was first achieved in 1995 by Eric Cornell and Carl Wieman (with $^{87}\text{Rb}$) and independently by Wolfgang Ketterle (with $^{23}\text{Na}$). Reaching the required nanokelvin temperatures demands a combination of laser cooling and evaporative cooling. The condensate is detected through time-of-flight imaging: when the trap is switched off, the condensate expands as a narrow peak in the momentum distribution, sharply distinct from the broad thermal cloud.

Fermi gases

Fermi gases are composed of fermions and governed by Fermi-Dirac statistics. The Pauli exclusion principle forces fermions to stack into progressively higher energy states, giving Fermi gases a nonzero energy and pressure even at absolute zero. This physics underlies the behavior of electrons in metals and the structure of compact stars.

Fermi-Dirac distribution function

The mean occupation number for fermions is:

$\langle n_i \rangle = \frac{1}{e^{(\epsilon_i - \mu)/k_B T} + 1}$

The only difference from the Bose-Einstein distribution is the $+1$ in the denominator, but the consequences are dramatic. This function is bounded between 0 and 1, enforcing the Pauli exclusion principle automatically. At $T = 0$, it becomes a sharp step function: states below $\mu$ are fully occupied, states above are empty.

Fermi energy and temperature

The Fermi energy $E_F$ is the energy of the highest occupied state at $T = 0$. For a uniform 3D gas:

$E_F = \frac{\hbar^2}{2m}(3\pi^2 n)^{2/3}$

The Fermi temperature is defined as $T_F = E_F / k_B$. It sets the scale for when quantum effects matter:

• $T \ll T_F$: the gas is degenerate and strongly quantum mechanical.
• $T \gg T_F$: the gas behaves classically.

For conduction electrons in a typical metal, $T_F \sim 10^4 - 10^5$ K, so electrons are deeply degenerate at room temperature.

Density of states

The density of states $g(\epsilon)$ counts how many single-particle states exist per unit energy. For a 3D uniform gas with spin degeneracy already included:

$g(\epsilon) = \frac{V}{2\pi^2}\left(\frac{2m}{\hbar^2}\right)^{3/2}\sqrt{\epsilon}$

This $\sqrt{\epsilon}$ dependence is specific to 3D free particles. In lower dimensions or with external potentials, the functional form changes. The density of states is essential for converting sums over states into integrals when computing thermodynamic quantities.

Degenerate vs non-degenerate regimes

• Degenerate regime ($T \ll T_F$): The Fermi-Dirac distribution approaches a step function. Nearly all states below $E_F$ are filled, and only a thin shell of states near $E_F$ (of width $\sim k_B T$) participates in thermal excitations. This is why the electronic specific heat of metals is so small at room temperature.
• Non-degenerate regime ($T \gg T_F$): Thermal energy smears out the distribution so broadly that the $+1$ in the denominator becomes irrelevant. The Fermi-Dirac distribution reduces to the Maxwell-Boltzmann distribution, and the gas behaves classically.

Quantum gases in harmonic traps

Most experiments on ultracold gases confine atoms in harmonic potentials (created by magnetic fields or focused laser beams) rather than in uniform boxes. The trapping potential modifies the density of states, which in turn changes the critical temperature, condensate fraction, and other thermodynamic quantities.

Density of states in traps

For a 3D isotropic harmonic trap with frequency $\omega$, the density of states is:

$g(\epsilon) = \frac{\epsilon^2}{2(\hbar\omega)^3}$

Compare this to the $\sqrt{\epsilon}$ dependence for a uniform gas. The $\epsilon^2$ dependence in a trap leads to different power-law scalings in nearly every thermodynamic quantity.

Bose-Einstein condensation in traps

The critical temperature for a trapped ideal Bose gas of $N$ particles is:

$T_c = \frac{\hbar\omega}{k_B}\left(\frac{N}{\zeta(3)}\right)^{1/3}$

where $\zeta(3) \approx 1.202$. Below $T_c$, the condensate fraction follows:

$\frac{N_0}{N} = 1 - \left(\frac{T}{T_c}\right)^3$

Note the exponent is 3 here (not 3/2 as in the uniform case). The condensate sits at the center of the trap in the harmonic oscillator ground state wavefunction, while the thermal (non-condensed) atoms form a broader cloud surrounding it.

Fermi gases in traps

For an isotropic harmonic trap, the Fermi energy is:

$E_F = \hbar\omega(6N)^{1/3}$

At $T = 0$, the density profile of a trapped Fermi gas follows an inverted parabola (the Thomas-Fermi profile), with a sharp edge at the classical turning point. Trapped degenerate Fermi gases are widely used to study strongly interacting fermionic systems, including the BEC-BCS crossover.

Thermodynamic properties

The thermodynamic behavior of quantum gases differs qualitatively from classical ideal gases, especially at low temperatures. The differences show up clearly in the temperature dependence of energy, heat capacity, and pressure.

Internal energy

For a Bose gas below $T_c$ (uniform 3D):

$U = \frac{3}{2}\frac{\zeta(5/2)}{\zeta(3/2)} N k_B T \left(\frac{T}{T_c}\right)^{3/2}$

For a degenerate Fermi gas ($T \ll T_F$, uniform 3D):

$U = \frac{3}{5}N E_F \left[1 + \frac{5\pi^2}{12}\left(\frac{T}{T_F}\right)^2 + \cdots \right]$

The leading term $\frac{3}{5}N E_F$ is the zero-point energy of the filled Fermi sea. The correction is quadratic in $T/T_F$, reflecting the fact that only fermions near the Fermi surface can be thermally excited.

Heat capacity

For a Bose gas below $T_c$:

$C_V \propto T^{3/2} \quad \text{(uniform)}, \qquad C_V \propto T^3 \quad \text{(harmonic trap)}$

For a degenerate Fermi gas:

$C_V = \frac{\pi^2}{2} N k_B \frac{T}{T_F}$

The linear-in-$T$ heat capacity of a Fermi gas is a hallmark of degeneracy. It arises because only a fraction $\sim T/T_F$ of the fermions near the Fermi surface can absorb thermal energy. This linear term is directly measured in the low-temperature specific heat of metals.

Pressure and equation of state

Both quantum gases satisfy $P = \frac{2U}{3V}$ (the same relation as the classical ideal gas, which follows from the $\epsilon \propto p^2$ dispersion). However, the internal energy $U$ itself differs from the classical value, so the pressure does too.

For a degenerate Fermi gas at $T = 0$:

$P = \frac{2}{5} n E_F$

This is nonzero even at absolute zero. There is no classical analog for this: a classical ideal gas has zero pressure at $T = 0$. This degeneracy pressure is what supports white dwarfs and neutron stars against gravitational collapse.

Applications and phenomena

Photon gas and blackbody radiation

Photons are massless bosons with zero chemical potential ($\mu = 0$, because photon number is not conserved). Applying Bose-Einstein statistics to photons in thermal equilibrium gives Planck's radiation law:

$u(\nu, T) = \frac{8\pi h \nu^3}{c^3} \frac{1}{e^{h\nu / k_B T} - 1}$

This spectral energy density explains the shape of the blackbody spectrum and historically motivated the development of quantum mechanics. The cosmic microwave background radiation, with its near-perfect blackbody spectrum at $T \approx 2.725$ K, is a direct application.

Electron gas in metals

Conduction electrons in metals form a degenerate Fermi gas with $T_F \sim 10^4$ K. This explains several otherwise puzzling observations:

• Electronic specific heat is much smaller than the classical prediction ($\frac{3}{2} N k_B$) because only electrons within $\sim k_B T$ of the Fermi energy participate.
• Pauli paramagnetism: the magnetic susceptibility is nearly temperature-independent, because only the thin shell of electrons near $E_F$ can flip their spins in response to a magnetic field.

Neutron stars and white dwarfs

These are extreme astrophysical examples of degenerate Fermi gases:

• White dwarfs are supported against gravity by electron degeneracy pressure. The maximum mass a white dwarf can have is the Chandrasekhar limit, $M_{\text{Ch}} \approx 1.44 \, M_\odot$, derived from the relativistic Fermi gas model.
• Neutron stars are supported by neutron degeneracy pressure (and nuclear forces). Their cores reach densities several times that of atomic nuclei.

Quantum degeneracy

A gas becomes quantum degenerate when the thermal de Broglie wavelength $\lambda_{dB} = \sqrt{2\pi\hbar^2 / (m k_B T)}$ becomes comparable to the mean interparticle spacing $n^{-1/3}$. At this point, the wave packets of neighboring particles overlap, and quantum statistics can no longer be ignored.

Quantum degeneracy pressure

In a degenerate Fermi gas, the Pauli exclusion principle forces particles into high-momentum states even at $T = 0$. The resulting pressure scales as:

$P \propto n^{5/3}$

for non-relativistic fermions. This pressure does not depend on temperature and cannot be removed by cooling. It is the mechanism that prevents white dwarfs and neutron stars from collapsing under their own gravity.

Pauli exclusion principle effects

The exclusion principle has consequences far beyond degeneracy pressure:

• At $T = 0$, fermions fill a Fermi sea in momentum space up to the Fermi momentum $p_F = \hbar(3\pi^2 n)^{1/3}$.
• The shell structure of atoms (and therefore all of chemistry) arises because electrons cannot share the same set of quantum numbers.
• In multi-electron systems, the antisymmetry requirement produces exchange interactions that have no classical counterpart.

Quantum vs classical behavior

The key contrasts between quantum and classical ideal gases:

• A classical gas has zero energy and zero pressure at $T = 0$. A Fermi gas does not.
• Bose-Einstein condensation is a purely quantum phase transition with no classical analog.
• Thermodynamic quantities (heat capacity, pressure, compressibility) follow different power laws in the quantum regime.
• At sufficiently high temperature, all quantum gases recover classical Maxwell-Boltzmann behavior. The crossover scale is set by $T_F$ for fermions and $T_c$ for bosons.

Experimental techniques

Laser cooling and trapping

Laser cooling exploits the momentum kick atoms receive when absorbing and re-emitting photons. The basic steps:

1. Doppler cooling: Six counter-propagating laser beams, red-detuned from an atomic resonance, preferentially slow atoms moving toward each beam. This creates an effective viscous force, cooling the gas to the Doppler limit (typically ~100 μK).
2. Magneto-optical trap (MOT): Adding a spatially varying magnetic field to the laser beams creates a position-dependent restoring force, trapping the cooled atoms.
3. Sub-Doppler cooling (e.g., Sisyphus cooling): Exploits polarization gradients to cool below the Doppler limit, reaching ~1 μK.

Evaporative cooling

Laser cooling alone cannot reach quantum degeneracy. Evaporative cooling bridges the gap:

1. Atoms are transferred to a conservative trap (magnetic or optical).
2. The trap depth is gradually lowered, allowing the most energetic atoms to escape.
3. The remaining atoms rethermalize through collisions at a lower temperature.
4. Repeating this process drives the temperature into the nanokelvin range.

This technique sacrifices atom number for lower temperature. Sympathetic cooling extends the method to species that don't collide efficiently by letting them thermalize with a co-trapped, efficiently cooled species.

Detection methods

• Time-of-flight imaging: The trap is switched off and the cloud expands freely. After a set time, an absorption image is taken. The resulting density profile maps the original momentum distribution. A BEC appears as a sharp central peak.
• In-situ absorption imaging: A resonant laser beam passes through the trapped cloud, and its shadow reveals the spatial density distribution.
• Bragg spectroscopy: Two laser beams with a controlled frequency difference probe the excitation spectrum and coherence properties of the gas.
• Noise correlation measurements: Shot-to-shot fluctuations in absorption images reveal quantum statistical correlations (bunching for bosons, antibunching for fermions).

Theoretical approaches

Grand canonical ensemble

The grand canonical ensemble is the natural framework for quantum gases because particle number fluctuates between quantum states. The grand partition function is:

$\ln \mathcal{Z} = -\sum_i \ln\left(1 \mp e^{-\beta(\epsilon_i - \mu)}\right)$

where the upper sign ($-$) applies to bosons and the lower sign ($+$) to fermions, and $\beta = 1/(k_B T)$. The chemical potential $\mu$ is fixed by requiring the average total particle number to equal $N$.

All thermodynamic quantities follow from derivatives of $\ln \mathcal{Z}$:

• Average particle number: $\langle N \rangle = \frac{1}{\beta}\frac{\partial \ln \mathcal{Z}}{\partial \mu}$
• Pressure: $P = \frac{1}{\beta V} \ln \mathcal{Z}$
• Internal energy: $U = -\frac{\partial \ln \mathcal{Z}}{\partial \beta}\bigg|_{\mu}$

Occupation numbers

The mean occupation of state $i$ is the central quantity in quantum statistical mechanics:

• Bosons: $\langle n_i \rangle = \frac{1}{e^{\beta(\epsilon_i - \mu)} - 1}$
• Fermions: $\langle n_i \rangle = \frac{1}{e^{\beta(\epsilon_i - \mu)} + 1}$

Once you know these, you compute any thermodynamic property by summing (or integrating, using the density of states) over all single-particle states. The entire structure of ideal quantum gas theory flows from these two formulas.

Partition function for quantum gases

Writing out the grand partition function explicitly:

• Bose gas: $\ln \mathcal{Z} = -\sum_i \ln\left(1 - e^{-\beta(\epsilon_i - \mu)}\right)$
• Fermi gas: $\ln \mathcal{Z} = \sum_i \ln\left(1 + e^{-\beta(\epsilon_i - \mu)}\right)$

In practice, the sum over states is converted to an integral using the density of states: $\sum_i \to \int_0^\infty g(\epsilon) \, d\epsilon$. This step is what connects the microscopic quantum statistics to macroscopic thermodynamic predictions. The choice of $g(\epsilon)$ (uniform box vs. harmonic trap vs. other geometry) determines the specific results for each physical system.