Ideal Bose Gas

A Bose–Einstein condensate (BEC) is a state of matter of a dilute gas of bosons cooled to temperatures very close to 0K (usually ~100nK in experiments). Under such conditions, a large fraction of bosons occupies the lowest quantum state, at which point macroscopic quantum phenomena become apparent.

Overview

In this section, we summarize some basic and useful thermodynamic results for Bose-Einstein condensation in a uniform, non-interacting gas of bosons. Most of the discussion here will be limited to 3D case. Physics for lower dimensions will be mentioned in the end.

Thermodynamics of a Bose Gas

Phase Space Density

The fundamental difference between a BEC and a classical gas is the occupancy of a single-particle state. In a classical gas, the mean occupation number for a single quantum state satisfies the Boltzmann distribution $\langle n_{\nu }\rangle =e^{-(\epsilon _{\nu }-\mu )/kT}$ which is much less than unity. This feature is qualitatively captured by the ${\textit {Phase-spaceDensity}}$ defined as (3D, homogeneous gas)

$\rho _{D}=n\lambda _{T}^{3}\,,$

where $\lambda _{T}=(2\pi \hbar ^{2}/mkT)^{1/2}$ is the thermal de Broglie wavelength.

Some typical parameters for

• Classical thermal gas
  • Atom density $n\sim 10^{25}{\text{m}}^{-3}$
  • Interatomic distance $(1/n)^{1/3}\sim 3{\text{nm}}$
  • Thermal de Broglie wavelength Failed to parse (syntax error): {\displaystyle \lambda_T \sim 10^{-2} \text{nm}（T = 300K) }
  • $\rho _{D}\sim 10^{-8}$
• BEC in dilute gas
  • Atom density $n\sim 10^{20}{\text{m}}^{-3}$
  • Interatomic distance $(1/n)^{1/3}\sim 10^{2}{\text{nm}}$
  • Thermal de Broglie wavelength Failed to parse (syntax error): {\displaystyle \lambda_T \sim 10^{3} \text{nm}（T = 100nK) }
  • $\rho _{D}\sim 10^{2}$

The Bose-Einstein Distribution

For non-interacting bosons in thermodynamic equilibrium, the mean occupation number of the single-particle state $\nu$ is

$f(\epsilon ,\mu ,T):=\langle n_{\nu }\rangle ={\frac {1}{e^{(\epsilon _{\nu }-\mu )/kT}-1}}\,.$

At high temperature, the chemical potential lies below $\epsilon _{min}$. As temperature is lower, the chemical potential rises until it reaches $\epsilon _{min}$ and the mean occupation numbers increase.

Transition Temperature - Semi-classical Picture

When $\nu =\epsilon _{min}$ , the occupation number on the ground state can be arbitrarily large, indicating the emrgence of a condensate. The corresponding temperature is the transition temperature $T_{c}$. $T_{c}$ can be calculated with the critieria that the maximum number of particles can be held in the excited states is equal to the total particle number $N$. In the semi-classical limit where the sum over all states is replaced by an integral and simple assumption that $\epsilon _{min}=0$ we have

$N_{ex}=\int \limits _{0}^{\infty }f(\epsilon ,\mu =0,T_{c})g(\epsilon )d\epsilon \,.$

Here $g(\epsilon )$ is the density of states. The number of atoms in the ground state is $N_{0}=N-N_{ex}$. The form of the transition temperature $T_{c}$ and therefore the condensate atom number depends strongly on the form of $g(\epsilon )$ which is affected by the dimension, trapping potential and the dispersion of the system. Under the most general assumption that $g(\epsilon )\propto \epsilon ^{\alpha -1}$, we reach

$N\propto C_{\alpha }(kT_{c})^{\alpha }\int \limits _{0}^{\infty }dx{\frac {x^{\alpha -1}}{e^{x}-1}}\,.$

where $x=\epsilon /(kT_{c})$. Straightforwardly we have a simple scaling function

$kT_{c}\propto N^{\frac {1}{\alpha }}\,.$

Some common cases are summarized below

	cases:	3D box	2D box	3D Harmonic	2D Harmonic
	$\alpha$	$3/2$	$1,diverge$	$3$	$2$

We see that the semi-classical picture is already good enough to capture some basic condensate physics.

Thermodynamic Properties

The thermodynamic properties $H$can be readily calculated from the Bose distributions and sum over all the states.

$H=\int \limits _{0}^{\infty }H(\epsilon )f(\epsilon ,\mu ,T)g(\epsilon )d\epsilon \,.$

For example, the total energy

$E=\int \limits _{0}^{\infty }\epsilon f(\epsilon ,\mu ,T)g(\epsilon )d\epsilon \propto (kT)^{\alpha }\int \limits _{0}^{\infty }{\frac {x^{\alpha }}{e^{x}-1}}dx\,.$

We therefore can obtain a scaling law for all the importnat thermodynamic quantities as listed below assuming that $g(\epsilon )\propto \epsilon ^{\alpha -1}$ and $H\propto T^{\delta }$

Thermodynamic Property:	E	$C_{V}=\partial E/\partial T$	Entropy S	2D Harmonic
$\delta$	$\alpha +1$	$\alpha$	$\alpha$	$2$