# 6.7: Bose-Einstein Statistics

- Page ID
- 18918

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

Perhaps you thought the Fermi-Dirac results were strange: Non-interacting particles forming a collection as hard as steel. . . room temperature being effectively zero. Wait until you see the Bose-Einstein results.

## 6.7.1 Theory

For independent bosons, whether free or subject to an external potential, the mean occupation number function is

\[ b(\mathcal{E})=\frac{1}{e^{\beta(\mathcal{E}-\mu)}-1}.\]

To begin to understand this function, note that

\[ \begin{array}{cc}{\text { when: }} & {\text { we have: }} \\ \hline \mathcal{E}<\mu & {b(\mathcal{E}) \text { negative }} \\ {\mathcal{E}=\mu} & {b(\mathcal{E})=\infty} \\ {\mathcal{E}>\mu} & {b(\mathcal{E}) \text { positive }}\end{array}\]

Thus this function has the general character sketched below:

Although we have written \( b ( \mathcal{E})\) as a function of the continuous variable \( \mathcal{E}\), we will in fact have occasion to evaluate it only at the energy eigenvalues \(\epsilon_1,\epsilon_2, \epsilon_3,. . . , \epsilon_r,. . . \), eigenvalues which of course differ for different external potentials. It seems bizarre that b(E) can be negative, and indeed this is only a mathematical artifact: Recall that in our derivation of the Bose function we needed to assume that *µ* < \(\varepsilon_1\) in order to insure convergence (see equation 6.35). Evaluating \( b ( \mathcal{E})\) at any eigenvalue r will always result in a positive mean occupation number.

The character of the Bose function is dominated by the singularity at *\(\mathcal{E}\)* = µ, so in trying to understand the function and its physical implications one must first locate the chemical potential. This section will provide a tour of Bose-Einstein behavior with decreasing temperature, throughout which the chemical potential shifts to the right. (When we investigated Fermi-Dirac behavior, we started at *T* = 0 and toured with increasing temperature, so the chemical potential shifted left.) This rightward shift presents a potential problem, because as the temperature decreases µ might shift right all the way to \(\epsilon_1\), and we know that *µ* < 1. We will just have to go ahead and take the tour, being wary because we know that a *µ* = 1 road block might pop up right in front of us as we view the countryside. With any luck *µ* will not yet have reached \(\epsilon_1\) when our tour halts at *T* = 0.

For the case of free and independent bosons (subject to periodic boundary conditions), the ground level energy is \(\epsilon_1\) = 0. The natural first step is to find *µ*(*T, V, N*) by demanding that

\[ N=\int_{0}^{\infty} G(\mathcal{E}) b(\mathcal{E}) d \mathcal{E}.\]

Natural though this may be, caution is in order. Remember that the integral above is an approximation to the sum over discrete energy levels

\[ N=\sum_{r} b\left(\epsilon_{r}\right).\]

It is legitimate to replace the sum with the integral when the value of \( b ( \mathcal{E})\) changes little from one energy level to the next. We saw in section 6.5 that in the thermodynamic limit, the free-particle level spacings approach zero, so for the most part this approximation is excellent. . . even exact in the thermodynamic limit. But there is one exception: At \(\epsilon\) = *µ*, the Bose function \( b ( \mathcal{E})\) is infinite, so *b*(µ) = ∞ is very different from *b*(µ + δ), no matter how small the positive number δ is. Usually we can ignore this caution, because *µ* < 1 = 0. But if *µ* approaches 0 then we expect the approximation (6.71) to fail.

In summary, the integral (6.71) is a good approximation for the sum (6.72) as long as the integrand varies slowly. Now for any value of \(\mathcal{E}\) greater than µ, you can make b(\(\mathcal{E}\) + δ) very close to b(E) simply by choosing δ > 0 small enough. This is what happens in the thermodynamic limit. But for \(\mathcal{E}\) = µ, there is always a large difference between \( b ( \mathcal{E})\) = ∞ and *b*(\(\mathcal{E}\) + δ), which is finite. Thus the integral approximation will be a good one as long as we avoid \(\mathcal{E}\) = µ.

In situations where the sum can be legitimately replaced with the integral, we have

\[ N=\int_{0}^{\infty} G(\mathcal{E}) b(\mathcal{E}) d \mathcal{E}\]

\[ =V\left[\frac{\sqrt{2 m^{3}}}{2 \pi^{2} \hbar^{3}}\right] \int_{0}^{\infty} \sqrt{\mathcal{E}} \frac{1}{e^{\beta(\mathcal{E}-\mu)}-1} d \mathcal{E}.\]

Use of the obvious substitution \(x = \beta \mathcal{E}\) gives

\[ N=V\left[\frac{\sqrt{2 m^{3}}}{2 \pi^{2} \hbar^{3}}\right]\left(k_{B} T\right)^{3 / 2} \int_{0}^{\infty} \frac{x^{1 / 2}}{e^{x} e^{-\beta \mu}-1} d x,\]

and remembering the definition (5.4) of thermal de Broglie wavelength results in a more convenient expression

\[ N=\frac{2}{\sqrt{\pi}} \frac{V}{\lambda^{3}(T)} \int_{0}^{\infty} \frac{x^{1 / 2}}{e^{x} e^{-\beta \mu}-1} d x.\]

Note that the definite integral above is not “just a number”. . . it is a function of the product *βµ*. There is no closed-form expression for the integral, although it is readily evaluated numerically and can be found in tabulations of functions. However it is easy to find an analytic upper bound: Because *µ* < 0, we have e^{−βµ} > 1 whence

\[ \int_{0}^{\infty} \frac{x^{1 / 2}}{e^{x} e^{-\beta \mu}-1} d x<\int_{0}^{\infty} \frac{x^{1 / 2}}{e^{x}-1} d x.\]

The expression on the right *is* just a number, and a little bit of work (see problem 6.26) shows that it is the number

\[ \zeta\left(\frac{3}{2}\right) \frac{\sqrt{\pi}}{2},\]

where the Riemann zeta function is defined by

\[ \zeta(s)=\sum_{n=1}^{\infty} \frac{1}{n^{s}} \quad \text { and } \quad \zeta\left(\frac{3}{2}\right)=2.612375348 \dots\]

So, how does this upper bound help us? It shows that

\[ N=\frac{2}{\sqrt{\pi}} \frac{V}{\lambda^{3}(T)} \int_{0}^{\infty} \frac{x^{1 / 2}}{e^{x} e^{-\beta \mu}-1} d x<\frac{V}{\lambda^{3}(T)} \zeta\left(\frac{3}{2}\right),\]

but what is the significance of this result? Remember that the upper bound is just the value of the function when µ = 0, which is exactly where we expect a road block due to the breakdown of the integral approximation (6.71). Our hopes that we could avoid the issue have been dashed. The breakdown occurs at the temperature *T*_{0} that satisfies

\[ N=\frac{V}{\lambda^{3}\left(T_{0}\right)} \zeta\left(\frac{3}{2}\right) \quad \text { or } \quad \lambda^{3}\left(T_{0}\right)=\zeta\left(\frac{3}{2}\right) / \rho\]

or

\[ T_{0}(\rho)=\frac{h^{2}}{2 \pi m k_{B} \zeta^{2 / 3}\left(\frac{3}{2}\right)} \rho^{2 / 3}.\]

For temperatures above* T*_{0}, the chemical potential shifts right as the temperature falls, and the integral approximation (6.71) is legitimate (indeed, exact in the thermodynamic limit). But as the temperature decreases below *T*_{0}, the chemical potential sticks at µ = 0 and the integral approximation (6.71) must be modified.

What is the proper modification? The function \(b(\mathcal{E})\) is slowly varying for all values of \(\mathcal{E}\) except \(\mathcal{E}\) = µ, so the integral approximation is legitimate for all the energy levels except the ground level at 1 = 0 = µ. Only for the ground level is some other result needed, so we just add in the ground level occupancy by hand:

\[ N=\left\langle n_{1}\right\rangle+\int_{0}^{\infty} G(\mathcal{E}) b(\mathcal{E}) d \mathcal{E}.\]

Now, we have already seen that when µ = 0—the condition for validity of this equation—the integral can be evaluated exactly and we have

\[ N=\left\langle n_{1}\right\rangle+\frac{V}{\lambda^{3}(T)} \zeta\left(\frac{3}{2}\right).\]

Note that \( \langle n_1 \rangle\) here is not given by its traditional formula (6.45), because

\[ \frac{1}{e^{\beta(\mathcal{E}-\mu)}-1}=\frac{1}{e^{0}-1}=\infty.\]

Instead, equation (6.84) is the formula for \( \langle n_1 \rangle \) when *T* < *T*_{0}. The mean ground level occupancy \( \langle n_1 \rangle \) is an intensive quantity when *T* > *T*_{0} but an extensive quantity when *T* < *T*_{0}.

In summary, the correct normalization equation breaks into two parts, namely

\[ N=\left\{\begin{array}{ll}{\frac{2}{\sqrt{\pi}} \frac{V}{\lambda^{3}(T)} \int_{0}^{\infty} \frac{x^{1 / 2}}{e^{x} e^{-\beta \mu}-1} d x} & {\text { for } T>T_{0}(\rho)} \\ {\left\langle n_{1}\right\rangle+\frac{V}{\lambda^{3}(T)} \zeta\left(\frac{3}{2}\right)} & {\text { for } T<T_{0}(\rho)}\end{array}\right.\]

We should expect that each part will behave quite differently, i.e. we expect a sudden change of behavior as the temperature drops through *T*_{0}.

What does all this mean physically? A naive reading of equation (6.80) suggests an upper bound on the number of particles that can be placed into the volume *V*. This would be sensible if the particles were marbles with hard-core repulsions. But these are non-interacting particles! Surely we can add more particles just by throwing them into the container. Indeed we can do so, and the associated excess mean occupation number is due to the level that particles like best at low temperatures, namely the ground level. The ground level thus has a much higher mean occupancy than the first excited level, and this rapid variation of mean occupancy with energy renders the approximation of sum by integral invalid. The inequality (6.80) does not limit the number of particles in the system: instead it shows the domain within which it is legitimate to approximate the sum (6.72) by the integral (6.71).

The abrupt transition at *T*_{0}(ρ) is called Bose-Einstein condensation and the material at temperatures below *T*_{0}(ρ) is called the Bose condensate. These terms are unfortunate: they conjure images of a gas condensing into a liquid, in which circumstance the atoms separate into two different classes: those in the liquid and those remaining in the vapor. This suggests that in Bose-Einstein condensation too there are two classes of particles: those in the ground level and those in the excited levels. This picture is totally false. It is incorrect to say “one particle is in the ground level, another is in the fourth excited level”. In truth the individual particles are not in individual levels at all: instead the whole system is in a state produced by multiplying together the individual level wavefunctions (“building blocks”) and then symmetrizing them. The literature of Bose-Einstein statistics is full of statements like “at temperatures below *T*_{0}, any particle added to the system goes into the ground level.” Such statements are wrong. They should be replaced with statements like “at temperatures below *T*_{0}, any increase in particle number occurs through an increase in \( \langle n_1 \rangle\), the mean occupancy of the ground level.” Or alternatively, “at temperatures below *T*_{0}, it is very likely that many of the building blocks from which the system wavefunction is built are the ground level.” Or again, to be absolutely precise, “at temperatures below *T*_{0}, if the energy is measured then it is very likely that many of the building blocks from which the resulting energy eigenfunction is built are the ground level.” Read again the previous paragraph—the one that begins “What does all this mean physically?”. Notice that I never need to say that a particle “is in” or “goes into” a given level.

## 6.7.2 Experiment

References: M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman, E.A. Cornell, “Observation of Bose-Einstein condensation in a dilute atomic vapor”, *Science*, **269** (14 July 1995) 198–201; Randall G. Hullet at Rice; Malcolm W. Browne, “Two groups of physicists produce matter that Einstein postulated”, *New York Times*, 14 July 1995, page 1.

## 6.7.3 Problems

*6.23 Character of the Bose function *

What are the limits of the Bose function \(b( \mathcal{E})\) (equation 6.70) as \(\mathcal{E}\) → ±∞? Is the curvature of the function greater when the temperature is high or when it is low?

*6.24 Thermodynamics of the Bose condensate *

For temperatures less than the Bose condensation temperature *T*_{0}, find the energy, heat capacity, and entropy of an ideal gas of spin-zero bosons confined to a volume *V*. Write your answers in terms of the dimensionless integral

\[ I=\int_{0}^{\infty} \frac{x^{3 / 2}}{e^{x}-1} d x\]

but don't bother to evaluate it. Show that

\[ C_{V}=\frac{5}{2} \frac{E}{T} \quad \text { and } \quad S=\frac{5}{3} \frac{E}{T}.\]

*6.25 More thermodynamics of the Bose condensate *

For the system of the previous problem, show that

\[ F=-\frac{2}{3} E \quad \text { and } \quad p=\frac{2}{3} \frac{E}{V}.\]

From this show that at low temperatures, the pressure of a collection of free and independent bosons goes like *p* ∼ *T*^{5/2}. (This is always less than the classical pressure *p* ∼ *T*.)

*6.26 An integral important for Bose condensation *

Show that

\[ \int_{0}^{\infty} \frac{x^{1 / 2}}{e^{x}-1} d x=\zeta\left(\frac{3}{2}\right) \frac{\sqrt{\pi}}{2}\]

where

\[ \zeta(s)=\sum_{n=1}^{\infty} \frac{1}{n^{s}}.\]

Clue:

\[ \frac{1}{e^{x}-1}=\frac{1}{e^{x}\left(1-e^{-x}\right)}=e^{-x} \sum_{n=0}^{\infty} e^{-n x} \quad \text { for } \quad x>0.\]