Gas In A Box Article Index for
Gas
Website Links For
Gas
 

Information About

Gas In A Box




equilibrium situation for a quantum ideal gas in a box which is a box
containing a large number particles which do not interact with each other except for
instantaneous thermalizing collisions. This simple model can be used to describe the
classical Ideal Gas as well as the various quantum ideal gases such as the ideal
massive Fermi Gas , the ideal massive Bose Gas as well as Black Body
radiation which may be treated as a massless Bose gas.

Using the results from either Maxwell-Boltzmann Statistics , Bose-Einstein Statistics
or Fermi-Dirac Statistics we use the Thomas-Fermi Approximation
and go to the limit of a very large box, and express the degeneracy of the energy
states as a differential, and summations over states as
integrals. We will then be in a position to calculate the thermodynamic properties of
the gas using the Partition Function or the Grand Partition Function . These
results will be carried out for both massive and massless particles. More complete
calculations will be left to separate articles, but some simple examples will be
given in this article.


THOMAS-FERMI APPROXIMATION FOR THE DEGENERACY OF STATES


For both massive and massless Particles In A Box , the states of a particle are
enumerated by a set of quantum numbers
{Link without Title} . The absolute value of the momentum is given by:

:p= rac{h}{2L}\sqrt{n_x^2+n_y^2+n_z^2}~~~~~~~~~n_i=1,2,3,\ldots

where ''h '' is Planck's Constant and ''L '' is the length of a side of the box.
We can think of each possible state of a particle as a point on a 3-dimensional
grid of positive integers. The distance from the origin to any point will be

:n=\sqrt{n_x^2+n_y^2+n_z^2}= rac{2Lp}{h}

Suppose each set of quantum numbers specify ''f '' states where ''f '' is
the number of internal degrees of freedom of the particle that can be altered by
collision. For example, a spin 1/2 particle would have ''f=2'', one for each spin
state. The Thomas-Fermi Approximation assumes that the quantum numbers are so large
that they may be considered to be a continuum. For large values of ''n '', we can
estimate the number of states with absolute value of momentum less than or equal to
''p '' from the above equation as

:
g=\left( rac{f}{8} ight) rac{4}{3}\pi n^3
= rac{4\pi f}{3} \left( rac{Lp}{h} ight)^3


which is just ''f '' times the volume of a sphere of radius ''n '' divided by eight
since we only consider the octant with positive ''n''. The number of
states with absolute value of momentum between ''p '' and ''p+dp '' is
therefore

:
dg= rac{\pi}{2}~f n^2\,dn = rac{4\pi fV}{h^3}~ p^2\,dp


where ''V=L3 '' is the volume of the box. Notice that in using this
continuum approximation, we have lost the ability to characterize the low-energy
states including the ground state where ''n=1''. For most cases this
will not be a problem, but when considering Bose-Einstein condensation, in which a
large portion of the gas is in or near the ground state, we will need to recover the
ability to deal with low energy states.

Without using the continuum approximation, the number of particles with
energy εi is given by

:
N_i = rac{g_i}{\Phi}


where

with β = ''1/kT '' with ''k '' being Boltzmann's Constant ,
''T '' being Temperature , and μ being the chemical potential. Using the
continuum approximation, the number of particles ''dN '' with energy between
''E '' and ''E+dE '' is now written:

:dN= rac{dg}{\Phi}


THE ENERGY DISTRIBUTION FUNCTION


We are now in a position to determine some distribution functions for the "gas in
a box". The distribution function for any variable A is PAdA and
is equal to the fraction of particles which have values for ''A '' between
''A '' and ''A+dA''

:P_A~dA = rac{dN}{N} = rac{dg}{N\Phi}

It follows that:

:\int_A P_A~dA = 1

The distribution function for the absolute value of the momentum is:

:P_p~dp = rac{Vf}{N}~ rac{4\pi}{h^3\Phi}~p^2dp

and the distribution function for the energy ''E '' is:

:
P_E~dE = P_p rac{dp}{dE}~dE


For a particle in a box (and for a free particle as well), the relationship between energy ''E '' and momentum ''p '' is different for massive and massless particles. For massive particles, we have

: E= rac{p^2}{2m}

while for massless particles:

:E=pc\,

where ''m '' is the mass of the particle and ''c '' is the speed of light.
Using these relationships we have:


  • For massive particles

    :
    dg = \left( rac{Vf}{\Lambda^3} ight)
    rac{2}{\sqrt{\pi}}~\beta^{3/2}E^{1/2}~dE


    :
    P_E~dE = rac{1}{N}\left( rac{Vf}{\Lambda^3} ight)
    rac{2}{\sqrt{\pi}}~ rac{\beta^{3/2}E^{1/2}}{\Phi}~dE


    where Λ is the Thermal Wavelength of the gas.

    :
    \Lambda =\sqrt{ rac{h^2 \beta }{2\pi m}}


    This is an important quantity, since when Λ is on the order of the
    interparticle distance ''(V/N)1/3 '', quantum effects begin to
    dominate and the gas can no longer be considered to be a Maxwell-Boltzmann gas.

  • For massless particles

    :dg = \left( rac{Vf}{\Lambda^3} ight) rac{1}{2}~\beta^3E^2~dE

    :
    P_E~dE = rac{1}{N}\left( rac{Vf}{\Lambda^3} ight)
    rac{1}{2}~ rac{\beta^3E^2}{\Phi}~dE


    where Λ is now the Thermal Wavelength for massless particles.

    :\Lambda = rac{ch\beta}{2\,\pi^{1/3}}



SPECIFIC EXAMPLES


The following sections give an example of results for some specific cases.


Massive Maxwell-Boltzmann particles


For this case:

:\Phi=e^{\beta(E-\mu)}\,

Integrating the energy distribution function and solving for ''N'' gives

:N = \left( rac{Vf}{\Lambda^3} ight)\,\,e^{\beta\mu}

Substituting into the original energy distribution function gives

:P_E~dE = 2 \sqrt{ rac{\beta^3 E}{\pi}}~e^{-\beta E}~dE

which are the same results obtained classically for the
Maxwell-Boltzmann Distribution . Further results can be found in the article on the
Classical Ideal Gas .


Massive Bose-Einstein particles


For this case:

:\Phi=e^{\beta \epsilon}/z-1\,

where ''z'' is defined as

:z=e^{\beta\mu}\,

Integrating the energy distribution function and solving for ''N'' gives
the particle number

:N = \left( rac{Vf}{\Lambda^3} ight) extrm{Li}_{3/2}(z)

Where Lis(z) is the Polylogarithm function and Λ is the
Thermal Wavelength . The polylogarithm term must always be positive
and real, which means its value will go from 0 to ζ(3/2) as ''z '' goes from
0 to 1. As the temperature drops towards zero, Λ will become larger and larger,
until finally Λ will reach a critical value Λc where ''z=1'' and

:N = \left( rac{Vf}{\Lambda_c^3} ight)\zeta(3/2)

The temperature at which Λ=Λc is the critical temperature. For
temperatures below this critical temperature, the above equation for the particle number
has no solution. The critical temperature is the temperature at which a Bose-Einstein
condensate begins to form. The problem is, as mentioned
above, that the ground state has been ignored in the continuum approximation. It turns
out, however, that the above equation for particle number expresses the number of bosons in excited states
rather well, and so we may write:

:
N= rac{g_0
z}{1-z}+\left( rac{Vf}{\Lambda^3} ight) extrm{Li}_{3/2}(z)


where the added term is the number of particles in the ground state. (The ground
state energy has been ignored.) This equation will hold down to zero temperature.
Further results can be found in the article on the ideal Bose Gas .


Massless Bose-Einstein particles (e.g. black body radiation)


The most common massless Bose gas is a gas of Photons in a Black Body . Taking the "box" to be a black body cavity, the photons are continually being absorbed and re-emitted by the walls. The number of photons is not a conserved quantity, since we are actually in a relativistic regime. In the derivation of Bose-Einstein Statistics , when the restraint on the number of particles is removed, this is effectively the same as setting the chemical potential (''μ'') to
zero. Therefore:

:
\Phi=e^{\beta E}-1\,


The spectral energy density (energy per unit volume per unit frequency) is

:U_
u~d
u = rac{NE}{V} P_E rac{dE}{d
u}~d
u
= rac{4\pi f h
u^3 }{c^3}~ rac{1}{e^{h
u/kT}-1}~d
u

where ''E''=''hν'' has been used, and since the radiation is the same in all directions, and propagates at the speed of light (''c''), the spectral intensity (energy/time/area/solid angle/frequency) is

:I_
u = rac{U_
u\,c}{4\pi}

which yields

:
I_
u~d
u = rac{f h
u^3 }{c^2}~ rac{1}{e^{h
u/kT}-1}~d
u


Since photons have two spin states, we have ''f=2'' which yields
Planck's Law Of Black Body Radiation . Note that if we had carried out this
procedure for massless Maxwell-Boltzmann particles, we would recover
Wien's Distribution which approximates a Planck's distribution for high temperature or low density.

Another massless Bose gas is given by the Debye Model for heat capacity. This considers a gas of Phonons in a box and differs from the development for photons in that the speed of the phonons is less than light speed, and there is a maximum allowed wavelength for each axis of the box. This means that the integration over phase space cannot be carried out to infinity, and instead of results being expressed in Polylogarithm s, they are expressed in the related Debye Function s.


Massive Fermi-Dirac particles (e.g. electrons in a metal)


For this case:

:\Phi=e^{\beta(E-\mu)}+1\,

Integrating the energy distribution function gives

:N=\left( rac{Vf}{\Lambda^3} ight)\left[- extrm{Li}_{3/2}(-z) ight]

Where again, Lis(z) is the Polylogarithm function and Λ is the
Thermal De Broglie Wavelength . Further results can be found in the article on the
ideal Fermi Gas .


REFERENCES


  • Huang, Kerson, "Statistical Mechanics", John Wiley & Sons, New York, 1967ě

  • A. Isihara, "Statistical Physics", Academic Press, New York, 1971

  • L. D. Landau and E. M. Lifshitz, "Statistical Physics, 3rd Edition Part 1", Butterworth-Heinemann, Oxford, 1996

  • Zijun Yan, "General thermal wavelength and its applications", Eur. J. Phys. 21 (2000) 625-631; http://www.iop.org/EJ/article/0143-0807/21/6/314/ej0614.pdf