Ideal Bose Gas Equation of State

1. The Grand Canonical Ensemble for Bosons

Just as with fermions, the grand canonical ensemble provides the most convenient framework for analyzing a gas of indistinguishable bosons. We consider a system of non-interacting bosons in thermal and chemical equilibrium with a reservoir at temperature $T$ and chemical potential $\mu$.

The fundamental quantity is the Bose-Einstein distribution, which gives the average number of particles $⟨n_{\epsilon }⟩$ in a single-particle state of energy $\epsilon$:

where $\beta =1/(k_{B}T)$. For the occupation number to be a positive real number, the denominator must be positive. This requires $e^{\beta (\epsilon -\mu )}>1$, which implies $\epsilon -\mu >0$ for all available energy states $\epsilon$. If we set the ground-state energy to be ${\epsilon }_0=0$, this imposes a critical constraint on the chemical potential of a Bose gas:

The fugacity, $z=e^{\beta \mu }$, is therefore restricted to the range $0<z<1$.

The following video lecture develops these initial concepts and proceeds to the crucial step of calculating the total number of particles in the system.

As the video highlights, the continuous density of states for a 3D gas, $g(\epsilon )\propto \sqrt{\epsilon }$, vanishes at $\epsilon =0$. A naive conversion of the sum over all states to an integral would incorrectly assign zero weight to the ground state. However, as $\mu \to 0^{-}$ (i.e., $z\to 1$), the occupation of the ground state, $n_0$, can become macroscopically large. This is the central feature of a Bose gas and necessitates splitting the total particle number $N$ into two parts: particles in the ground state ($N_0$) and particles in all excited states ($N_{ex}$).

For spin-$s$ bosons, this gives:

This separation is the key to correctly describing the system at all temperatures.

2. The Full Equation of State

With the formalism for the particle number established, we can now proceed to the complete derivation of the equation of state. This involves finding expressions for the particle density $n=N/V$ and the pressure $P$, both as functions of temperature $T$ and fugacity $z$. The following resource provides a rigorous and complete derivation, which aligns well with your mathematical background.

By following the derivation, you arrive at the parametric equation of state for an ideal Bose gas. For spin-$s$ bosons, it is given by the following pair of equations:

1. Particle Density: $n=\frac{N}{V}=\frac{2s+1}{{\lambda }^3}{g}_{3/2}(z)+n_0$
2. Pressure: $\frac{P}{k_{B}T}=\frac{2s+1}{{\lambda }^3}{g}_{5/2}(z)$

Here, $\lambda =\sqrt{2\pi {\hslash }^2/(m{k}_{B}T)}$ is the thermal de Broglie wavelength, and $n_0=(2s+1)/V(z^{-1}-1)$ is the density of particles in the ground state. These two equations, parameterized by the fugacity $z$, form the complete equation of state, valid at any temperature. To find the pressure of a gas with a given density $n$ and temperature $T$, one must first solve the particle density equation for $z$ and then substitute that value into the pressure equation.

3. Energy, Pressure, and the High-Temperature Limit

The derivation also reveals a simple relationship between the internal energy $U$ and the pressure $P$.

The result, $U=\frac{3}{2}PV$, is identical to that for a classical monatomic ideal gas and a non-relativistic Fermi gas. It is a general consequence of the quadratic dispersion relation $\epsilon \propto k^2$ for non-relativistic particles in three dimensions.

To gain physical intuition for the parametric equation of state, it's instructive to examine its behavior in the high-temperature (or low-density) limit. In this regime, quantum effects are expected to be small, and we should recover the classical ideal gas law. This corresponds to the limit where the fugacity $z\ll 1$.

The expansion yields:

This is a significant result. The first quantum correction due to Bose-Einstein statistics is negative, meaning the pressure of a Bose gas is slightly lower than that of a classical gas at the same density and temperature. This can be interpreted as an effective "attraction" between bosons arising from their statistical tendency to cluster together. This stands in direct contrast to the case of fermions, where the Pauli exclusion principle leads to a positive correction, an effective "repulsion" that increases the pressure.

Conclusion

In this lesson, we have constructed the equation of state for a gas of non-interacting bosons, valid across all temperatures. We navigated the subtleties of the grand canonical ensemble for bosons, particularly the need to treat the ground state with special care.

Key Takeaways:

• The chemical potential for a Bose gas must be negative ($\mu <0$), which constrains the fugacity to $0<z<1$.
• The macroscopic occupation of the ground state is possible as $\mu \to 0$. Therefore, the total particle number must be separated into ground-state and excited-state contributions: $N=N_0+N_{ex}$.
• The equation of state for an ideal Bose gas is given by a pair of parametric equations relating pressure $P$ and number density $n$ through the fugacity $z$:
  $\frac{P}{k_{B}T}=\frac{2s+1}{{\lambda }^3}{g}_{5/2}(z)\text{and}n=\frac{2s+1}{{\lambda }^3}{g}_{3/2}(z)+n_0$
• In the high-temperature limit, this equation reduces to the classical ideal gas law with a negative first-order quantum correction, reflecting an effective statistical attraction between bosons.

The machinery we have assembled today is precisely what we need to investigate the most striking prediction for ideal Bose gases. In our next lesson, we will explore what happens in the low-temperature, high-density regime. We will find that there is a critical temperature below which the excited states become "full," forcing a macroscopic number of particles into the ground state—a remarkable phase transition known as Bose-Einstein Condensation.