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