Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

Posts Tagged ‘Bose condensation’

A final pre-exam update of my notes compilation for ‘PHY452H1S Basic Statistical Mechanics’, Taught by Prof. Arun Paramekanti

Posted by peeterjoot on April 22, 2013

Here’s my third update of my notes compilation for this course, including all of the following:

April 21, 2013 Fermi function expansion for thermodynamic quantities

April 20, 2013 Relativistic Fermi Gas

April 10, 2013 Non integral binomial coefficient

April 10, 2013 energy distribution around mean energy

April 09, 2013 Velocity volume element to momentum volume element

April 04, 2013 Phonon modes

April 03, 2013 BEC and phonons

April 03, 2013 Max entropy, fugacity, and Fermi gas

April 02, 2013 Bosons

April 02, 2013 Relativisitic density of states

March 28, 2013 Bosons

plus everything detailed in the description of my previous update and before.

Posted in Math and Physics Learning. | Tagged: , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , | 1 Comment »

PHY452H1S Basic Statistical Mechanics. Lecture 20: Bosons. Taught by Prof. Arun Paramekanti

Posted by peeterjoot on April 2, 2013

[Click here for a PDF of this post with nicer formatting (especially if my latex to wordpress script has left FORMULA DOES NOT PARSE errors.)]


Peeter’s lecture notes from class. May not be entirely coherent.


In order to maintain a conservation of particles in a Bose condensate as we decrease temperature, we are forced to change the chemical potential to compensate. This is illustrated in fig. 1.1.

Fig 1.1: Chemical potential in Bose condensation region


Bose condensatation occurs for T < T_{\mathrm{BEC}}. At this point our number density becomes (except at \mathbf{k} = 0)

\begin{aligned}n(\mathbf{k}) = \frac{1}{{e^{\beta \epsilon_\mathbf{k}} - 1}}.\end{aligned} \hspace{\stretch{1}}(1.2.1)

Except for \mathbf{k} = 0, n(\mathbf{k}) is well defined, and not described by this distribution. We are forced to say that

\begin{aligned}N = N_0 + \sum_{\mathbf{k} \ne 0} n(\mathbf{k}) = N_0 + V\int \frac{d^3 \mathbf{k}}{(2 \pi)^3} \frac{1}{{ e^{\beta \epsilon_\mathbf{k}} - 1 }}.\end{aligned} \hspace{\stretch{1}}(1.2.1)

Introducing the density of states, our density is

\begin{aligned}\rho = \rho_0 + \int_0^\infty d\epsilon \frac{N(\epsilon)}{e^{\beta \epsilon} - 1 },\end{aligned} \hspace{\stretch{1}}(1.2.3)


\begin{aligned}N(\epsilon) = \frac{1}{{4 \pi^2}} \left( \frac{2m}{\hbar} \right)^{3/2} \epsilon^{1/2}.\end{aligned} \hspace{\stretch{1}}(1.2.4)

We worked out last time that

\begin{aligned}\rho = \rho_0 + \rho \left( \frac{T}{T_{\mathrm{BEC}}} \right)^{3/2},\end{aligned} \hspace{\stretch{1}}(1.2.4)


\begin{aligned}\rho_0 = \rho \left( 1 - \left( \frac{T}{T_{\mathrm{BEC}}} \right) ^{3/2} \right).\end{aligned} \hspace{\stretch{1}}(1.2.6)

This is plotted in fig. 1.2.

Fig 1.2: Density variation with temperature for Bosons


\begin{aligned}\rho_0 = \frac{N_{\mathbf{k} = 0}}{V}.\end{aligned} \hspace{\stretch{1}}(1.7)

For T \ge T_{\mathrm{BEC}}, we have \rho_0 = 0. This condensation temperature is

\begin{aligned}T_{\mathrm{BEC}} \propto \rho^{2/3}.\end{aligned} \hspace{\stretch{1}}(1.8)

This is plotted in fig. 1.3.

Fig 1.3: Temperature vs pressure demarkation by T_BEC curve


There is a line for each density that marks the boundary temperature for which we have or do not have this condensation phenomina where \mathbf{k} = 0 states start filling up.

Specific heat: T < T_{\mathrm{BEC}}

\begin{aligned}\frac{E}{V} &= \int \frac{d^3 \mathbf{k}}{(2 \pi)^3} \frac{1}{{ e^{\beta \hbar^2 k^2/2m} - 1}}\frac{\hbar^2 k^2}{2m} \\ &= \int_0^\infty d\epsilon N(\epsilon) \frac{1}{{ e^{\beta \epsilon} - 1 }} \epsilon \\ &\propto \int_0^\infty d\epsilon \frac{\epsilon^{3/2}}{ e^{\beta \epsilon} - 1 } \\ &\propto \left( k_{\mathrm{B}} T \right)^{5/2},\end{aligned} \hspace{\stretch{1}}(1.9)

so that

\begin{aligned}\frac{C}{V} \propto \left( k_{\mathrm{B}} T \right)^{3/2}.\end{aligned} \hspace{\stretch{1}}(1.10)

Compare this to the classical and Fermionic specific heat as plotted in fig. 1.4.

Fig 1.4: Specific heat for Bosons, Fermions, and classical ideal gases


One can measure the specific heat in this Bose condensation phenomina for materials such as Helium-4 (spin 0). However, it turns out that Helium-4 is actually quite far from an ideal Bose gas.

Photon gas

A system that is much closer to an ideal Bose gas is that of a gas of photons. To a large extent, photons do not interact with each other. This allows us to calculate black body phenomina and the low temperature (cosmic) background radiation in the universe.

An important distinction between a photon sea and some of these other systems is that the photon number is actually not fixed.

Photon numbers are not “conserved”.

If a photon interacts with an atom, it can impart energy and disappear. An excited atom can emit a photon and change its energy level. In a thermodynamic system we can generally expect that introducing heat will generate more photons, whereas a cold sink will tend to generate fewer photons.

We have a few special details that distinguish photons that we’ll have to consider.

  1. spin 1.
  2. massless, moving at the speed of light.
  3. have two polarization states.

Because we do not have a constraint on the number of particles, we essentially have no chemical potential, even in the grand canonical scheme.


\begin{aligned}\lambda = \left\{\begin{array}{l l}+1 & \quad \mbox{Right circular polarization} \\ -1 & \quad \mbox{Left circular polarization}\end{array}\right.\end{aligned} \hspace{\stretch{1}}(1.11)

Our number density, since we have no chemical potential, is of the form

\begin{aligned}n_{\mathbf{k}, \lambda}= \frac{1}{{e^{\beta \epsilon_{\mathbf{k}, \lambda}} - 1 }},\end{aligned} \hspace{\stretch{1}}(1.12)

Observe that the average number of photons in this system is temperature dependent. Because this chemical potential is not there, it can be quite easy to work out a number of the thermodynamic results.

Photon average energy density

We’ll now calculate the average energy density of the photons. The energy of a single photon is

\begin{aligned}\epsilon_{\mathbf{k}, \lambda} = \hbar c k = \hbar \omega,\end{aligned} \hspace{\stretch{1}}(1.2.13)

so that the average energy density is

\begin{aligned}\frac{E}{V} &= \sum_{\mathbf{k}, \lambda} \frac{1}{{ e^{ \beta \epsilon_\mathbf{k}} - 1}} \epsilon_\mathbf{k}\rightarrow\underbrace{2}_{\text{number of polarizations}}\int \frac{d^3 \mathbf{k}}{(2 \pi)^3}\frac{ \hbar c k}{ e^{ \beta \epsilon_\mathbf{k}} - 1} \\ &= 2 \int_0^\infty d\epsilon \underbrace{\frac{1}{{(2 \pi)^3}} 4 \pi \frac{\epsilon^2}{(\hbar c)^3} }_{\text{Photon density of states}}\frac{\epsilon}{e^{\beta \epsilon} - 1} \\ &= \frac{1}{{\pi^2}} \frac{1}{{ (\hbar c)^3 }} \int_0^\infty d\epsilon \frac{\epsilon^3}{e^{\beta \epsilon} - 1}\end{aligned} \hspace{\stretch{1}}(1.2.13)

Mathematica tells us that this integral is

\begin{aligned}\int_0^\infty d\epsilon \frac{\epsilon^3}{e^{\beta \epsilon} - 1} =\frac{\pi ^4}{15 \beta ^4},\end{aligned} \hspace{\stretch{1}}(1.2.13)

for an end result of

\begin{aligned}\frac{E}{V} =\frac{\pi^2}{15} \frac{1}{{(\hbar c)^3}} \left( k_{\mathrm{B}} T \right)^4.\end{aligned} \hspace{\stretch{1}}(1.2.13)

Phonons and other systems

There is a very similar phenomina in matter. We can discuss lattice vibrations in a solid. These are called phonon modes, and will have the same distribution function where the only difference is that the speed of light is replaced by the speed of the sound wave in the solid. Once we understand the photon system, we are able to look at other Bose distributions such as these phonon systems. We’ll touch on this very briefly next time.

Posted in Math and Physics Learning. | Tagged: , , , , , , , , , , , | Leave a Comment »