Bose-Einstein Statistics and Bose Gas

In the landscape of quantum statistical mechanics, the behavior of many-body systems is fundamentally dictated by the intrinsic nature of the particles involved. Based on the quantum mechanical rules they obey, particles are categorized into two distinct classes: Fermions and Bosons. While Fermions are governed by the Pauli Exclusion Principle—preventing any two particles from occupying the same quantum state—Bosons behave in a strikingly different manner. They are "social" particles, capable of occupying the same quantum state simultaneously. This unique property gives rise to Bose-Einstein Statistics and leads to extraordinary macroscopic phenomena, most notably Bose-Einstein Condensation (BEC).
The defining characteristic of a Boson lies in the symmetry of its wave function. For a system composed of $N$ identical particles, the total wave function $\Psi$ must remain unchanged when the coordinates of any two particles are interchanged:

$$\Psi(x_1, x_2, \dots, x_i, \dots, x_j, \dots, x_N) = \Psi(x_1, x_2, \dots, x_j, \dots, x_i, \dots, x_N)$$

This exchange symmetry dictates that Bosons must possess integer spin (e.g., $0, 1, 2, \dots$).

Crucially, because Bosons are truly identical and indistinguishable, we cannot apply the counting methods used in classical Maxwell-Boltzmann statistics. In classical mechanics, we might distinguish between "particle A in state 1 and particle B in state 2" versus "particle B in state 1 and particle A in state 2." In Bose-Einstein statistics, such distinctions are meaningless. We are concerned only with the occupancy number—how many particles reside in a specific energy level—rather than which specific particle is where.

Derivation of the Bose-Einstein Distribution

To determine how particles are distributed across various energy levels, we utilize the Grand Canonical Ensemble. This framework allows the number of particles in a system to fluctuate, which is mathematically convenient for treating individual energy levels as subsystems.

Consider a system with multiple energy levels $\epsilon_i$. Each level can accommodate an arbitrary number of particles, $n_i$. For a specific level $\epsilon_i$, the grand canonical partition function $z_i$ is given by the sum over all possible occupancy numbers:

$$z_i = \sum_{n_i=0}^{\infty} e^{-\beta (\epsilon_i - \mu) n_i}$$

Here, $\beta = 1/k_B T$ (where $k_B$ is the Boltzmann constant and $T$ is temperature) and $\mu$ is the chemical potential. This expression is a geometric series. For the sum to converge, we must satisfy the condition $\epsilon_i - \mu > 0$ for all $i$. This implies that the chemical potential $\mu$ must always be strictly less than the lowest energy level of the system ($\epsilon_0$).

The sum of the geometric series yields:

$$z_i = \frac{1}{1 - e^{-\beta (\epsilon_i - \mu)}}$$

From this, we can derive the average occupancy $\langle n_i \rangle$ of the energy level $\epsilon_i$ using the standard thermodynamic relation:

$$\langle n_i \rangle = -\frac{1}{\beta} \frac{\partial \ln z_i}{\partial \epsilon_i} = \frac{1}{e^{\beta(\epsilon_i - \mu)} - 1}$$

This is the Bose-Einstein distribution function. The presence of the "$-1$" in the denominator is the mathematical engine behind the unique behavior of Bosons. As the energy $\epsilon_i$ approaches the chemical potential $\mu$, the denominator approaches zero, causing the occupancy number to grow rapidly. This explains why Bosons tend to cluster in low-energy states.

The Ideal Bose Gas Model

An ideal Bose gas is a theoretical model consisting of non-interacting Bosons. To describe the macroscopic properties of such a gas in three-dimensional space, we transition from discrete energy levels to a continuous density of states, $g(\epsilon)$.

For non-relativistic particles of mass $m$ in a volume $V$, the density of states is typically expressed as:

$$g(\epsilon) d\epsilon = C V \epsilon^{1/2} d\epsilon$$

where $C$ is a constant involving $\hbar$ and $m$. The total number of particles $N$ and the total internal energy $U$ are then calculated by integrating the distribution function over all energies:

$$N = \int_{0}^{\infty} \frac{g(\epsilon)}{e^{\beta(\epsilon - \mu)} - 1} d\epsilon$$
$$U = \int_{0}^{\infty} \frac{\epsilon \cdot g(\epsilon)}{e^{\beta(\epsilon - \mu)} - 1} d\epsilon$$

As the temperature $T$ decreases, the chemical potential $\mu$ must increase to keep the total particle number $N$ constant. Eventually, $\mu$ approaches the ground state energy $\epsilon_0$ (often set to zero).

Bose-Einstein Condensation (BEC)

Bose-Einstein Condensation is a remarkable quantum phase transition that occurs when the temperature of a Bose gas drops below a certain critical threshold.

The Physics of Condensation

When the temperature falls below a critical value $T_c$, the excited states (higher energy levels) can no longer accommodate the total number of particles in the system. The "excess" particles are forced to collapse into the ground state $\epsilon_0$.

At this point, a macroscopic fraction of the total particles occupies a single quantum state. This results in a single, coherent quantum wave function that describes a large portion of the system, effectively allowing quantum mechanical effects to be observed on a macroscopic scale.

Mathematically, for $T < T_c$, the integral approximation for $N$ fails because the ground state contribution becomes significant. We must split the total particle count into two parts:

$$N = N_0 + N_{ex}$$

where $N_0$ is the number of particles in the ground state (the condensate) and $N_{ex}$ is the number of particles in the excited states.

Determining the Critical Temperature $T_c$

By evaluating the limit where $\mu \to 0$ for the excited state integral, we can derive the expression for the critical temperature:

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

In this formula, $n = N/V$ is the particle density and $\zeta(3/2) \approx 2.612$ is the Riemann Zeta function. Two key insights emerge:

• Density Dependence: Higher particle densities lead to higher critical temperatures.
• Mass Dependence: Lighter particles undergo condensation more easily (at higher temperatures).

Characteristics of the Condensate

Once the system enters the BEC state, it exhibits exotic properties such as superfluidity, where the fluid flows without viscosity. The system essentially becomes a "dual" entity: a collection of "classical-like" particles in excited states and a "quantum condensate" in the ground state. This phenomenon has been experimentally validated using ultra-cold atomic gases, such as Rubidium-87.

Summary and Comparison

To clarify the distinction between the two fundamental types of quantum statistics, the following table provides a direct comparison:

Bose-Einstein statistics serve as the theoretical cornerstone for understanding not only the behavior of matter but also the nature of light (photons being massless Bosons with $\mu=0$). It remains a vital tool in modern condensed matter physics, driving research into superconductivity, superfluidity, and the frontiers of ultra-cold atomic physics.