Bose-Einstein distribution


The Bose-Einstein distribution is a Probability distribution that describes the behavior of a system of NN non-interacting bosons in thermal equilibrium. Its Probability density function is

ni=1eβ(εiμ)1=1z1eβεi1\langle n_{i} \rangle =\frac{1}{e^{\beta(\varepsilon_{i}-\mu)}-1}=\frac{1}{z^{-1}e^{\beta \varepsilon_{i}}-1}

ni\langle n_{i} \rangle is the average number of fermions in the ii-th single-particle state of energy εi\varepsilon_{i}, μ\mu is the system's chemical potential, β=1/kBT\beta=1/k_{B}T is the inverse temperature, with kBk_{B} the Boltzmann constant and TT the temperature, and zz is the fugacity. The normalization constant is

ini=N\sum_{i}n_{i}=N

Bosons are not subject to the Pauli exclusion principle, so ni\langle n_{i} \rangle can be any positive real.

Bose-Einstein condensation

Note that this function diverges if eβ(εiμ)=1e^{\beta(\varepsilon_{i}-\mu)}=1, which is to say μ=0\mu=0 and ε0\varepsilon\to 0. In this state, the ground state's occupation number blows up to infinity and all other states becomes unoccupied, which means that all bosons "collapsed" into the ground state. This is a phase transition called Bose-Einstein condensation and it only occurs in bosons due to the 1-1 term at the denominator. Neither fermions nor classical particles exhibit this behavior.