Degenerate Gases

What if we want to incorporate rotation and vibration?In this case the energy of the system becomes eeneint where en is the energy of the orbital and eint is the energy associated with the rotation and vibration of the particle.If we require that the average occupancy of an orbital is small, n << 1, then the grand sum for the orbital is just

formula

or

formula

But the summation is just the partition function of the internal states, Zint.Thus the grand sum becomes


 

formula
(14.1)

We can think of Zint = Zrot + Zvib + Zelectronic.

What is the distribution function?

formula

but <N(en)> << 1 in the classical regime, so


 

formula
(14.2)

What about the quantum regime, formula?We can reach the quantum regime by increasing the concentration or lowering the temperature.The temperature for the quantum regime is found from nq to be any temperature below


 

formula
(14.3)

A gas in the quantum regime with a temperature t << t0 is often called a degenerate gas.

One of the early indications that something was wrong with classical physics came from the entropy for a classical gas.Nerst found that as formulas diverges as ln t.With the advent of quantum mechanics, this problem was resolved since now the entropy is calculated using either the Fermi-Dirac or Bose-Einstein distribution functions, both of which go to a single ground state as formula.Physically we say that the entropy was squeezed out on cooling the quantum gas.

Fermionic Degenerate Gases

Lets consider an electron gas.In this case it follows Fermi-Dirac statistics, and the ground state has a distribution function of 

formula

where I have used the fact that meF at t = 0.We would like to calculate eF.To do this, recall that if we consider waves we found that the energy could be written in terms of the wave number kformula, where k = (nxi + nyj + nzk) with nx, ny, nz = ±1, ±2, ±3, .…Thus, in k-space, each state is a point.We can apply this to our Fermi gas since we know from quantum mechanics that we can associate a wave with a particle.The density of states was calculated earlier and was found to be

D = 2V/(2p)3

where the factor of two comes from the fact that each fermion has two possible spin states associated with an energy state.Now let the temperature go to zero.At t = 0, all of the particles have fallen into the lowest energy available to them within the constraints of the Pauli exculsion priciple.This means that in k-space the sphere of radius kF is totally occupied with particles.Then the total number of particles is given by

formula

This can be inverted to find the radius kF


 

KF = (3p2n)1/3
(14.4)

Substituting this into ek yields


 

formula
(14.5)

We can use ek to determine the total energy in the ground state


 

formula
(14.6)

Thus we see that the average kinetic energy per particle is 0.6 eF.At constant N, as the energy increases the volume decreases, so that the Fermi energy gives a repulsive contribution to the binding of any material; that is, the Fermi energy tends to increase the volume.In most metals, and in white dwarf and neutron stars, it is the most important repulsive interaction.It is balanced in metals by the Coulomb attraction between electrons and ions and in stars by gravitational attraction.

How did we get D(e)de?

formula

where N(k) is the number of states with wave vector less than or equal to k.

formula

So

formula

Now recall that

formula

so that

formula

and thus


 

formula
(14.7)

Notice that, since formula, and recalling that formula, we get that the average occupancy looks like

graph

To derive the heat capacity of a degenerate fermionic gas, we need to consider the effect of raising the temperature.As the temperature goes from 0 to t, the energy in the system is raised from U0 to U(t) .The difference in energy can be calculated from


 

formula
(14.8)

Using the fact that

formula

(14.8) can be rewritten as

formula

or


 

formula
(14.9)

The first integral on the right hand side gives the energy needed to take electrons from eF to the orbitals of energy eeF and the second integral gives the energy needed to bring the electrons to eF from orbitals below eF.

Once we have DU, we can find the heat capacity by differentiating DU with respect to t.The only temperature dependent terms are the f(e), and so we get

formula

It is a good approximation to evaluate D(e) = D(eF), allowing it to be removed from the integral.Finally, we can also assume that the chemical potential varies slowly with t, allowing us to replace m with eF in f(e).Substituting these two approximations into the integral and evaluating, we get


 

formula
(14.10)

 Another way of obtaining this result is to notice that we can associate t with the magnitude of the energy change, and D(eF)t with the number of electrons whose energy changes.Then

formula

and so

formula

Returning to (14.10), if we recall that formula, then we see that

formula

whertFeF.

Bosonic Degenerate Gases

Now lets look at bosonic gases.Bosons differ from fermions in that there is no limit on the number of particles that can occupy the ground state.Thus, we get a remarkable effect with bosonic gases, namely that there is a certain transition temperature below which a substantial fraction of particles will occupy the ground orbital.Any other orbital will be occupied by a relatively negligible number of particles.This is known as Einstein condensation.To understand Einstein condensation, look at the behavior of the chemical potential at low temperatures.
Let the energy of the ground orbital be e0 = 0.The distribution function of the ground state is

formula

Whenformula, the occupancy of the ground orbital becomes equal to the total number of particles in the system, so that


 

formula
(14.11)

where we have used the series expansion formula.This can be rearranged to yield


 

mt/N
(14.12)

asformula.Note that the chemical potential is negative.This is because the chemical potential in a boson system must always be lower in energy than the ground state orbital in order that the occupancy of every orbital be non-negative.

Example:

Consider an atom in a cube of side l. The energy states of the atom are given by

formula

where nx, ny and nz are positive integers.In units of formula, we see that the energy of the ground state is e(111) = 3, while the set of next lowest orbitals all have energy e(211) = 6.Thus the lowest excitation energy of the atom is

formula

If we let the atom be a normal helium atom with a mass of m(4He) = 6.6 x l0-24 g and the size of the cube be l = 1 cm, then De = 2.48 x l0-30 erg.This corresponds to a temperature of 1.8 x l0-14 K.This seems to be an insignificant temperature, but for N = 1022 atoms, at t = 1 mK, m = -1.4 x 10-41 erg, where e(111) is taken as the zero energy orbital.Thus l is much closer to e(111) than to any other state, and e[e(111) - m]/t is much closer to 1 than e[e(211) - m]/t, so e(111) dominates the distribution function.The fractional occupancy of the first excited orbital is given by

formula

Thus we see that the Bose-Einstein distribution favors the situation where the greatest part of the population is left in the ground orbital at sufficiently low temperatures.The particles left in the ground orbital, as long as their numbers are >> 1, are called the Bose-Einstein condensate.

Temperature and Orbital Occupancy

What is the relationship between temperature and the occupancy of an orbital?Recall that

formula

Then the total number of 4He atoms in the ground and excited orbitals is given by


 

formula
(14.13)

where N0(t) is the number of atoms in the ground orbital at temperature t.We call N0 the number of atoms in the condensed phase and formulathe number of atoms in the normal phase.Using the Bose-Einstein distribution function and setting e = 0, we see that


 

N0(t) = l/(1 - l)
(14.14)

Similarly, Ne is found to be

formula

We can evaluate this integral by noting that whenever N0 >> 1, i.e. at sufficiently low temperatures, the formula.Setting l = 1, we get 


 

formula
(14.15)

We can write N0 and Ne as fractional occupancies by dividing each by the total number of atoms.This leads to


 

formula
(14.16)

and


 

Ne = 2.612 nq/n
(14.17)

Einstein Condensation Temperature

We define the Einstein condensation temperature,tE, as the temperature for which the number of atoms in the excited states is equal to the total number of atoms, that is Ne = N.From (14.15) we see that tE is given by

 
formula
(14.18)

This allows (14.17) to be rewritten as

formula

The number of atoms in the ground orbital is found from


 

formula
(14.19)

A graph of the fractional occupancy of the two phases yields

graph

Example:

The condensation temperature for liquid helium is found from

formula

where VM is the molar volume in cm3/mole and M is the molecular weight.For liquid helium VM = 27.6 cm3/mole and M = 4, so TE = 3.1 K.In reality, a transition is observed in liquid 4He at 2.17 K.A graph of the heat capacity as a function of temperature shows this transition clearly

graph

It is believed that at 2.17 K there is a condensation of liquid4He into the ground orbital of the system.This is different from the condensation in coordinate space that occurs in the condensation of a gas into a liquid.The interatomic forces that lead to liquefaction of 4He at 4.2 K does not destroy the effects of boson condensation, so in this respect the liquid behaves as a gas.The condensation is normally not permitted for fermions, but pairs of fermions may act as bosons, as in the superconductivity of electron pairs (Cooper pairs) in metals.In liquid 4He, we call the state above 2.17 K He I, and the state below 2.17 K He II.

The new state of liquid 4He below 2.17 K has some interesting properties.The viscosity drops to essentially zero and the thermal conductivity becomes very high.Thus, below 2.17 K,we say that 4He acts as a superfluid.If we look at a phase diagram of 4He, we see that it has two triple points

graph

The upper triple point is the transition of the solid, which does not exist below a pressure of 25 atm, into either He I or He II.The triple point is 1.74 K and about 30 atm.The second triple point is where the vapor and liquid phases meet.It is interesting to note that at normal pressure the solid phase does not ever occur, even at absolute zero.The two triple points are connected by a line which separates the regions of He I and He II.

How does liquid 4He move without viscosity?We can investigate this by realizing that the superfluid component of 4He behaves as if it were in a vacuum.The N0 atoms are condensed into the ground orbital and have no excitation energy.The superfluid has an energy only when the center of mass of the superfluid is given a velocity relative to the laboratory reference frame.The condensed component of the superfluid will flow with zero viscosity so long as the flow does not create excitations in the superfluid, i.e. as long as there are no transitions from the ground orbital to the excited orbitals.Thus, the criterion for superfluidity involves the energy and momentum relationship of the excitations in He Il.If the excited orbitals were like the orbitals of free atoms, then we could show that superfluidity would not be expected.However, the low energy excitations do not resemble free particle excitations, but are longitudinal sound waves (longitudinal phonons).We call the low lying excited states elementary excitations, and in their particle aspects they are called quasiparticles.

Elementary Excitations

Consider a body of mass M0 falling with velocity V down a column of liquid helium at rest at absolute zero.In order to generate an elementary excitation of energy ek and momentum formula, we must satisfy

 
 
½M0V2 = ½M0V'2ek
(14.20)

where V', is the velocity after the excitation, and


 
 

M0V = M0V' + formula
(14.21)

In order to find the minimum value of V which can generate an excitation, lets solve(14.21) for V' and use it to eliminate the kinetic energy in (14.20).Solving (14.21) forV' we get

formula

or

formula

Substituting this into (14.20) yields


 

formula
(14.22)

Whenformula, we can get the minimum of V as

formula

This can be simplified if we assume that M0 is large compared to the phonon momentum, so that


 

Vminek/formula
(14.23)

A body moving slower than this will not be able to create excitations in the liquid, so the motion will be frictionless.We can see this result geometrically if we plot the energy of excitation as a function of momentum.If we construct a straight line which just meets the curve from below, the slope of that line is the minimum velocity.

graph

If formula, as for the excitation of a free atom, the straight line has zero slope and the critical velocity is zero.The energy of a low energy phonon in liquid He II is formula, where vs is the velocity of sound in liquid He II.The critical velocity is equal to sound if Vcrit = vs is valid for all wavevectors, which it isn't for He II.The observed critical flow velocities are considerably lower than the velocity of sound, presumably because the plot of ekversesformulamay turn downward at very high formula.The actual spectrum of elementary excitations in liquid He II has been determined by observations on the inelastic scattering of slow neutrons.For the wavevectors covered by these neutron experiments, the critical velocity is VC = 5 x 103 cm/s.