1 Lecture 7 Gaseous systems composed of molecules with internal motion. • Monatomic molecules. • Diatomic molecules. • Fermi gas. • Electron gas. • Heat capacity of electron gas.
1
Lecture 7Lecture 7
Gaseous systems composed of molecules with internal motion.
• Monatomic molecules.
• Diatomic molecules.
• Fermi gas.
• Electron gas.
• Heat capacity of electron gas.
2
Gaseous systems composed of molecules with internal motion
Gaseous systems composed of molecules with internal motionIn most of our studies so far we have consider only the translation part of the molecular motion.
Though this aspect of motion is invariably present in a gaseous system, other aspects, which are essentially concerned with the internal motioninternal motion of the molecules, also exist. It is only natural that in the calculation of the physical properties of such a system, contributions arising from these motions are also taken into account.In doing so, we shall assume here that
a) the effects of the intermolecular interactions are negligible and
3
n=N/VnMkT
nh here ,1
) 2(3
2/3
3
(7.1)
is fulfilled; effectively, this makes our system an ideal, Boltzmannian gas.
b) the nondegeneracy criterion
Under these assumptions, which hold sufficiently well in a large number of practical applications, the partition functionpartition function of the system is given by
NZN
Z 1!
1 (7.2)
where
)(31 Tj
VZ
(7.3)
21
2/ MkTh
4
The factor in brackets is the transitional partition function of a molecule, while the factor j(T) is supposed to be the partition function corresponding to the internal motions. The latter may be written as
kT
ii
iegTj /)( (7.4)
where i is the molecular energy associated with an
internal state of motion (which characterized by the quantum numbers i), while gi represents the
degeneracy of that state.
The contributions made by the internal motions of the molecules to the various thermodynamic quantities of the system follow straightforwardly from the function j(T). We obtain
5
(7.5)
(7.6)
(7.7)
(7.8)
Fint= - N kT lnj
int= - kT lnj
j
TTjNk ln
lnint
Eint=NkT2 jT
ln
j
TT
TNkCV ln
2
int
(7.9)
6
Thus the central problem in this study consists of deriving an explicit expression for the function j(T) from a knowledge of the internal states of the molecules. For this purpose, we note that the internal state of a molecule is determined by:
• electronic state,• state of nuclei, • vibrational state and• rotational state.
Rigorously speaking, these four modes of excitation mutually interact; in many cases, however, they can be treated independently of one another. We then write
j(T)=jelec(T) jnuc(T) jvib(T) jrot(T) (7.10)
with the result that the net contribution made by the internal motions to the various thermodynamic quantities of the system is given by a simple sum of the four respective contributions.
7
Monatomic moleculesMonatomic molecules
At the very outset we should note that we cannot consider a monatomic gas except at temperatures such that the thermal energy kTkT is small in comparison with the ionization energy EEionion; for
different atoms, this amounts to the condition: T<<ET<<Eionion/k/k101044-10-1055 ooK.K. At these temperatures the number of ionized atoms in the gas would be quite insignificantquite insignificant. The same would be true for atoms in excited states, for the reason that separation of any of the excited states from the ground state of the atom is generally comparable to the ionization energy itself. Thus, we may regard all the atoms of the gas to be in their (electronic) ground state.
8
Now, there is a well-known class of atoms, namely He, Ne, A,..., which, in their ground state, possess neither orbital angular momentum nor spin (L=S=0L=S=0).
Their (electronic) ground state is clearly a singlet: ggee=1=1. The nucleus, however, possesses a degeneracy, which arises from the possibility of different orientations of the nuclear spin. (As is well known, the presence of the nuclear spin gives rise to the so-called hyperfine structure in the electronic statehyperfine structure in the electronic state.
However the intervals of this structure are such that for practically all temperatures of interest they are small in comparison with kTkT.) If the value of this spin is SSnn, the corresponding degeneracy factor
ggnn=2S=2Snn+1+1. Moreover, a monatomic molecule is
incapable of having any vibrational or rotational vibrational or rotational statesstates
9
The internal partition function (7.10) of such a molecule is therefore given by
j(T)=ggr.st.=ge gn=2Sn+1 (7.11)
Equations (7.4-7.9) then tell us that the internal the internal
motionsmotions in this case contribute only towards properties such as the chemical potential and the chemical potential and the
entropy of the gasentropy of the gas; they do not make contribution towards properties such as the internal energy and the specific heat.In other cases, the ground state of the atom may possess both orbital angular momentum and spin (LL0,S0,S00-- as, for example, in the case of alkali atoms), the ground state would then possess a definite fine structure.
Fint= - N kT lnj
int= - kT lnj
j
TTjNk ln
lnint
Eint=NkT2 jT
ln
j
TT
TNkCV ln
2
int
10
J
kTelect
JeJTj /)12()( (7.12)
The intervals of this structure are in general, comparable with kT; hence, in the evaluation of the partition function, the energies of the various components of the fine structure must be taken into account. Since these components differ from one another in the value of the total angular momentum J, the relevant partition function may be written as
The forgoing expressions simplifies considerably in the following limiting cases:
kT >> all J ; then
J
elect SLJTj )12)(12()12()( (7.13)
11
kT<< all J ; thenkT
elect eJTj /0
0)12()( (7.14)
where J0 is the total angular momentum, and 0 the
energy of the atom in the lowest state.
In their case, the electronic motion makes no
contribution towards the specific heat of the gas.
And, in view of the fact that both at high
temperatures the specific heat tends to be equal to
the translational value 3/2 Nk, it must be passing
through a maximum at a temperature comparable
to the separation of the fine levels.
Needless to say, the multiplicity (2Sn+1) introduced
by the nuclear spin must be taken into account in
each case.
12
Diatomic moleculesDiatomic molecules
Now, just as we could not consider a monatomic gas except at temperatures for which kTkT is small compared with the energy of ionization, for similar reasons one may not consider a diatomic gas except at temperatures for which kTkT is small compared with the energy of dissociation; for different molecules this amounts once again to the
condition: T<<ET<<Edissdiss/k/k101044-10-1055 ooKK. At this temperatures the number of dissociated molecules in the gas would be quite insignificant.
13
Cv=(Cv)elec+(Cv)vib+(Cv)rot (7.15)
The heat capacitance of the diatomic gas is consist from three parts
At the same time, in most cases, there would be practically no molecules in the excited states as well, for the separation of any of these states from the ground state of the molecule is in general comparable to the dissociation energy itself.
Let us consider them consequently.
kTelec eggTj /
10)( (7.16)
In the case of electron contribution the electronic partition function can be written as follows
14
where gg00 and gg11 are degeneracy factors of the two
components while is their separation energy.
])/(1][)/(1[
)/()(
/01
/10
2
kTkTelecv eggegg
kTNkC
(7.17
)
We note that this contribution vanishes both for T<<T<</k/k and for T>>T>>/k/k and has a maximum value for a certain temperature /k/k; cf. the corresponding situation in the case of monatomic atom.
The contribution made by (7.16) towards the various thermodynamic properties of the gas can be readily calculated with the help of the formula (7.4-7.9). In particular we obtain for the contribution towards specific heat
j
TT
TNkCv ln
2
int
15
Let us now consider the effect of vibrational states of the molecules on the thermodynamic properties of the gas. To have an idea of the temperature range, over which this effect would be significant, we note that the magnitude of the corresponding quantum of energy, namely , for different diatomic gases is of order of 101033 ooKK. Thus we would obtain full contributions (consistent with the dictates of the equipartition theorem) at temperatures of the order of 101044 ooKK or more, and practically no contribution at temperatures of the order of 101022 ooKK or less. We assume, however, that the temperature is not high enough to excite vibrational states of large energy; the oscillations of the nuclei are then small in amplitude and hence harmonic.
16
ke
eT
NkC T
Tv
vibV v
v
v2/
/2 ;
)1()(
(7.18)
We note that for T>>T>>vv the vibrational specific heat is very nearly equal to the equipatition value NkNk; otherwise, it is always less than Nk. In particular, for T<<T<<vv , the specific heat tends to zero (see Figure 7.1); the vibrational degrees of freedom are then said to be "frozenfrozen".
The energy levels for a mode of frequency are then given by the well-known expression (n+1/2) (n+1/2)
hh/2/2. The evaluation of the vibrational partition
function jjvibvib(T)(T) is quite elementary. In view of the
rapid convergence of the series involved, the summation may formally be extended to n=n=. The corresponding contributions towards the various thermodynamic properties of the system are given by eqn.(4.64 -4.69). In particular, we have
17
Figure 7.1 The vibrational specific heat of a gas of diatomic molecules. At T=v the specific heat is already about 93 % of the equipartition value.
0 0.5 1.0 1.5 2.0
0.5
1.0
T/V
Cvib/Nk
18
Finally, we consider the effect of • the states of the nuclei and • the rotational states of the molecule: wherever necessary, we shall take into account the mutual interaction of these modes.
This interaction is on no relevance in the case of the heternuclearheternuclear molecules, such as AB; it is, however, important in the case of homonuclearhomonuclear molecules, such as AA.In the case of heternuclear molecules the states of the nuclei may be treated separately from the rotational states of the molecule. Proceeding in the same manner as for the monatomic molecules we conclude that the effect of the nuclear states is adequately taken care of through degeneracy factor gn. Denoting the spins of the two nuclei by SA and SB,
this factor is given by
19
gn= (2SA+1)(2SB+1) (7.19)
As before, we obtain a finite contribution towards the chemical potential and the entropy of the gas but none towards the internal energy and specific heat.Now, the rotational levels of a linear "rigid" with two degrees of freedom (for the axis of rotation) and the principle moments of inertia (I, I, 0), are given by
Illrot 2/)1( 2 (7.20)
here I=M(r0)2 , where M=m1m2/(m1+m2) is the reduced mass of the nuclei and r0 the equilibrium distance between them. The rotational partition function of the molecule is then given by
20
IkTlll
IkTlllTj
rr
l
lrot
2 ;)1(exp)12(=
2)1(exp)12()(
2
0
2
0
(7.21)
For T>>r the spectrum of the rotational states may
be approximated by a continuum.
The summation (7.21) is the replaced by integration:
r
rrot
Tdl
TlllTj
0
)1(exp)12()( (7.22)
21
The rotational specific heat is the given by
(CV)rot=Nk (7.23)
which is indeed consistent with equipartition theorem.
A better evaluation of the sum (7.21) can be made with the help of the Euler-Maclaurin formula
....)0()0(''')0(')0()()( 302401
7201
121
0 0
21
V
n
ffffdxxfnf
Putting Txxxxf r /)1(exp)12()(
one can obtain
....315
4
15
1
3
1)(
2
TT
TTj rr
rrot
(7.24)
22
which is the so-called Mulholland's formulaMulholland's formula; as expected, the main term of this formula is identical with the classical partition function (7.22). The corresponding result for the specific heat is
....1)(
3
94516
2
451
TTNkC rr
rotV (7.25)
which shows that at high temperatures the
rotational specific heat decreases with
temperatures and ultimately tends to the classical
value NkNk.
23
Thus, at high (but finite) temperatures the rotational specific heat of diatomic gas is greater than the classical value. On the other hand, it must got to zero as T 0. We, therefore, conclude that it must pass through at least one maximum. (See Figure 7.2)
Fig.7.2. The rotational specific heat of a gas of heteronuclear diatomic molecules.
0 0.5 1.0 1.5 2.0
0.5
1.0
T/V
Cvib/Nk
24
In the opposite limiting case, namely for T<<T<<rr ,
one may retain only the first few terms of the sum (7.21); then
....531)( /6/2 TTrot
rr eeTj (7.26)
whence one obtains, in the lowest approximation
TrrotV
reT
NkC /22
12)(
(7.27)
Thus, as TT 0 0, the specific heat drops exponentially to zero (Fig. 7.2). Now we can conclude that at low temperatures the rotational degrees of freedom of the molecules are also "frozen"."frozen".
25
Fermi gasFermi gasFermi gasFermi gas
Let us consider the perfect gas composed of fermions. Let us consider in this case the behavior of the Fermi functionFermi function given by equation (5.46)
1
1)(
)(
ef
for the case when the assembly is at the absolute zero of temperature and the Fermi energy is FF(0)(0)..
for > >FF(0),(0), {{--FF(0)}/kT=(0)}/kT= while for < <FF(0), {(0), {--FF(0)}/kT=-(0)}/kT=-.
When T=0 T=0 the quantity {{--FF(0)}/kT(0)}/kT has two possible values:
26
for <F(0),
for >F(0),
11
1)(
e
f
01
1)(
e
f
(7.28)
There are therefore two possible values of the Fermi function:
Classical tail
kT
T=0
0
1
0.5
0
f()
Figure 7.3. Fermi-Dirac distribution function plotted at absolute zero and at a low temperature kT<<. The Fermi level o at T=0 is shown.
27
Equation (7.28) implies that, at the absolute zero of temperature, the probability that a state with energy <<FF(0)(0) is occupied is unity, i.e such states are all occupied. Conversely, all states with energies >>FF(0)(0) will be empty. The form of f(f()) at T=0T=0 is shown as a function of energy in Figure 7.3.
Thus with only one fermion one fermion allowed per stateallowed per state, all the lowest states will be occupied until the fermions are all accommodated. The Fermi Fermi levellevel, in this case, is simply the highest occupied state and above this energy level the states are unoccupied.
Classical tail
kT
T=0
0
1
0.5
0
f()
This behavior may be explained as following. At the absolute zero of temperature, the fermionsfermions will necessarily occupy the lowest available energy occupy the lowest available energy states. states.
28
Classical tail
kT
T=0
0
1
0.5
0
f()
Figure 7.3. Fermi-Dirac distribution function plotted at absolute zero and at a low temperature kT<<. The Fermi level o at T=0 is shown.
For the temperatures T<<TT<<TFF/k/kFF/k/k the Fermi-Dirac Fermi-Dirac
distributiondistribution behavior is shown in the Fig. 7.3 by bold line. The Fermi temperature TTFF and the Fermi energy
FF are defined by the indicated identities. The Fermi
energy FF is defined as the value of the chemical
potential at the absolute temperature FF (0(0).
29
We note that f=1/2f=1/2 when ==. The distribution for T=0T=0 cuts off abruptly at ==, but at a finite temperature the distribution fuzzes out over a width of the order of several kTkT. At high energies -->>kT>>kT the distribution has a classical form.
The value of the chemical potential is a function of temperature, although at low temperature for an ideal Fermi gas the temperature dependence of may often be neglected.
The determination of (()) is often the most tedious stage of a statistical problem, particularly in ionization problems. We note that is essentially a normalization parameter and that the value must be chosen to make the total number of particles come out properly.
30
An important analytic property of ff at low temperatures is that -df/d-df/d is approximately a delta function. We recall the central property of the Dirac delta function (x-a):(x-a):
)()()( aFaxxF (7.29)
Now consider the integral
d
fF
0
)(
At low temperatures -df/d-df/d is very large for and is small elsewhere. Unless F(F()) is rapidly varying in this neighborhood we may replace it by F(F()) and the integral becomes
31
)0()()]()[(
)( 0
0
fFfFdf
F
(7.30)
)(
)(0
Fd
fF
But at low temperatures f(0)f(0)11, so that
(7.31)
a result similar to (7.29).
32
Electron gasElectron gas
The conducting electrons in a metal may be considered as nearly free, moving in a constant potential field like the particles of an ideal gas. Electrons have half-integral spinhalf-integral spin, and hence the Fermi-Dirac statisticsFermi-Dirac statistics are applicable to an ideal gas of electrons.
fk=f(k-) (7.32)
1
1)(
/
Tef (7.33)
We use kk to specify the state of the electron, and kk its energy. Let be the chemical potential of the electron. Each state can accommodate at most one electron. Let ffkk be the FDFD average population of state kk.
33
The energy distribution of the states is an important property. Let
)]([1
)( 0 kkV
g (7.34)
(T=0) (7.35)
The function g(g()) is the energy distribution of the states per unit volume, which we simply call density of states, and is the energy with respect
to 00. The calculation of g(g()) gives
2/102
2
03
3)](2[)
2(
)2(2)(
m
m
m
ppdg (7.36)
34
The thermodynamic properties of this model can be largely expressed in terms of f and g, e.g. the density N/V of the electrons and the energy density E/V are
)()(/ dfgVN (7.37)
))(()(/ 0 dfgVE (7.38)
If T=0, then all the low energy states are filled up to the Fermi surface. Above this surface all the states are empty. The energy at the Fermi surface is 0, i.e the chemical potential when T=0, and is always denoted by F:
F(T=0)=0(7.39)
The Fermi surface can be thought of as spherical surface in the momentum space of the electrons. The radius of the sphere pF is called Fermi momentum
FFm
p 2
2(7.40)
35
There are N N states with energy less than F:
NVpF 3
3
)2(
2
3
4
(7.41)
i.e (volume of sphere in momentum spacevolume of sphere in momentum space) (volumevolume) (spin state (=2spin state (=2)) (h)(h)33, with h/2h/2 =1 =1. Hence
pF=(32n)1/3 (7.42)
where nnN/VN/V. Let a=h/(2a=h/(2 p pFF)) is approximately the average distance between the electrons, hence:
2
2~
maF
(7.43)
is approximately the zero-point energy of each electron. This zero-point energy is a result of the wave nature of the electron or a necessary result of the uncertainty principle. To fix an electron to within a space of size aa, its momentum would have to be of order h/ (2h/ (2 a). a).
36
In most metals, the distance between electrons is about 1010-8-8 cmcm and FF ~1 eV ~10 ~1 eV ~104 o4 oKK (See Table 7.1). Therefore, at ordinary temperatures, T<<T<<FF, i.e. the temperature is very low, only electrons very close to the Fermi surface can be excited and most of the electrons remain inside the sphere experiencing no changes. 3/22
2)/3(
2VN
mF
F
Fp
VNk 3/12 )/3(
F
Fp
v
kT F
F
(7.44)
(7.45)
(7.46)
(7.47)
37
Table 7.1 Properties of the electron gas at the Fermi surface.
N/V cm-3 kF cm-1 vF cm/sec F eV TF=F/k
Li 4.61022 1.1108 1.3108 4.7 5.5104
Na 2.5 0.9 1.1 3.1 3.7
K 1.34 0.73 0.85 2.1 2.4
Rb 1.08 0.68 0.79 1.8 2.1
Cs 0.86 0.63 0.73 1.5 1.8
Cu 8.50 1.35 1.56 7.0 8.2
Ag 5.76 1.19 1.38 5.5 6.4
Au 5.90 1.2 1.39 5.5 6.4
38
Heat CapacityHeat Capacity
Only a very small portion of the electrons is influenced by temperature. Hence the concept of holes appears naturally. The states below the Fermi surface are nearly filled, and empty states are rare. We shall call an empty an empty state a holestate a hole. Now this model becomes a new mixed gas of holes new mixed gas of holes together with electrons above Fermi surfacetogether with electrons above Fermi surface. (We shall call these outer electrons.)
The momentum of a hole is less than ppFF, while that
of the outer electrons is larger than ppFF. The lower
the temperature TT , the more dilute the gas is. At T=0T=0, this gas disappears.
39
At T=0T=0 the total energy is zero, i.e. there are no holes or outer electrons. The holes are also fermions because each state has at most one hole. Hence, a state of energy --'' can produce a hole of energy ''.. The average population of the hole is (for each state):
1-f(-'-)= f('+) (7.48)
Now the origin of the energy is shifted to the energy at the Fermi surface, i.e. =0=0. The energy of a hole is the energy required taking an electron from inside the Fermi surface to the outside.
As T<<T<<FF, the energy of the holes or the outer
electrons cannot exceed TT by too much. In this interval of energy, g(g()) is essentially unchanged, i.e.
40
g()g(0)=mpF/2 (7.49)
Hence the energy distribution is the same for the holes or the outer electrons. Therefore, 00. All the calculations can now be considerably simplified, e.g. the total energy is
0
)()0(2)]0()([1 dfgETEV
(7.50)
where 2g(0) is the density of states of the holes plus the outer electrons.
The energy of a hole cannot exceed FF, but FF>>T>>T and
so the upper limit of the integral in (7.50) can be taken to be . This integration is easy:
41
0
22
2
0121
)( Tdxe
xTdf
x (7.51)
Substituting in (7.50), and differentiating once, we get heat capacity
TgTV
CV
)0(3
E 11 2
(7.52)
This result is completely different from that of the
ideal gas in which CCvv=3/2 N=3/2 N.
In that case each gas molecule contribute a heat capacity of 3/23/2. Now only a small portion of the electrons is involved in motion and the number of
active electrons is about NT/NT/FF
42
F
NTdfVg
~)()0(
0
(7.53)
Hence CCvv~N(T/~N(T/FF)).. From (7.52) we get
FFv T
TkN
TNC
22
22
(7.54)
Each active electron contribute approximately 11 to heat capacity CC.
43
Figure 7.4 Electron (hole) energy versus momentum.
F
-F
energy of a hole
energy of anouter elctron
pF