Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Смирнова Н.А. Методы статистической термодинамики в физической химии учеб. пособие

.pdf
Скачиваний:
190
Добавлен:
27.10.2023
Размер:
20.12 Mб
Скачать

отсчитываемая от этого нулевого значения, может быть представлена как сумма энергии поступательного движения молекулы и энергии внутренних движений (вращения молекулы как целого, колебаний ядер, возбужденных электронных состояний, — подробнее см. гл. I X , а также § 5 настоящей главы) Изучать электронные состояния мож­ но только на базе квантовой теории (см. гл. IX , § 4). При невысоких температурах вклад возбужденных электронных состояний в сред­ нее значение энергии обычно незначителен. В настоящей главе, ко­ торая посвящена классической теории идеального газа, мы ограни­ чимся рассмотрением поступательного движения молекул, вращения

молекул как целого и

внутримолекулярных

колебаний и запишем:

 

s = £пост + еар + екол-

( I V . 15)

В случае одноатомного

газа примем, что

 

 

г = Е п о с т .

(IV . 16)

Если газ находится во внешнем силовом поле, то следует к правой части выражения (IV. 15) или (IV. 16) добавить слагаемое, определяю­ щее потенциальную энергию молекулы в этом поле.

Кинетическая энергия

 

 

 

 

 

 

 

2

2

 

2

р а

 

 

 

 

 

 

 

 

 

 

 

Рх

Ру

 

Рг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Рх,

Ру,

Рг — составляющие

импульса

молекулы;

р — модуль

импульса)

всегда

входит

в энергию

молекулы е как

независимое

слагаемое;

члены"

е в р

и

г к о л

в

выражении

(IV. 15)

от

перемен­

ных рх,

Ру, рг не зависят. Заметим, что даже в выражении для полной

энергии системы взаимодействующих

частиц энергия поступательного

движения молекул представляет сумму независимых слагаемых

 

 

 

 

 

 

 

 

РІІ

+ Р% +

РІІ

 

 

 

 

 

 

 

"(Р.Я)

 

=

=1

І

 

-+H'{p',q),

 

 

(IV.18)

 

где p

не зависит от составляющих

импульса

молекул

pxl,

р

ргі

(t =

l , ....

N).

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим далее

распределение

молекул

по скоростям.

Распре­

деление по скоростям

было впервые

выведено Максвеллом

(1860 г.)

на основании молекулярно-кинетического подхода. Здесь мы выведем распределение Максвелла из формул (IV. 10), (IV. 15), (IV.17). Энер­ гию молекулы идеального газа можем представить в виде суммы

+ PI +

- + •' ,

( I V . 19)

 

 

где е' = е' (р/, а) не зависит от рх, ру, pz (р' — набор обобщенных импульсов, исключая рх, pv, pz). Переменные рх, pv, pz, p', q полностью определяют положение и движение молекулы.

100

Из формул (IV. 10) и (IV. 19) вытекает следующее выражение для вероятности заданного состояния молекулы:

dw(px,

ру, рг, p',

q) =

 

2ткТ

Jtf—

 

= Ае

е

dpxdpydpzdp'dq»

(IV.20)

Плотность распределения вероятностей

р(рх, ру, рг,

р', а), как мы

видим, распадается на произведение двух независимых сомножите­ лей, один из которых определяется составляющими импульса моле­ кулы рх, ру, рг, другой — остальными обобщенными импульсами и координатами. Такой вид функции р является следствием разделения соответствующих переменных в выражении для энергии (IV. 19), которое представляет сумму двух независимых слагаемых. Приходим к выводу, что распределение по импульсам рх, ру, рг и распределение

по другим переменным

независимы,

причем

 

 

 

 

 

pj +

РІ+РІ

 

 

р(Рх.

Ру, Рг ) = Яе

 

2mkT

 

 

 

.

 

(IV.21)

где В не есть функция

импульсов.

 

 

 

 

Вероятность заданных значений рх,

ру,

рг не зависит

от того, ка­

ковы значения других

переменных, и определяется

выражением

 

 

„ 2

, 2 ,

2

 

 

dw(px,

ру,

pz)=Be

2mkT

dpxdpydpz.

 

(IV.22)

 

 

Значения импульсов

задаем в интервалах: от рх

до рх

+ dpx для

составляющей по оси х, от ру до ру + dpy для составляющей по оси у , от рг до рг + dpz для составляющей по оси г; значения других пе­ ременных при этом могут быть любые. Результат (IV.22) соответствует

интегрированию

выражения

(IV.20) по

всем возможным

значени-

• ям р' и q.

 

 

 

 

 

 

Определим постоянную В из условия нормировки и будем при этом

считать, что составляющие импульса могут меняться

в пределах от

— с о до оо, хотя, строго говоря, они не могут быть бесконечно

боль­

шими. Однако

расширение пределов интегрирования до ± о о не ска-

жется на

величине интеграла, поскольку величина

exp I

- — —

(а = X, у , г) быстро убывает с ростом абсолютного значения ра

и при

больших ра

практически равна нулю; большие значения ра не дают

ощутимого

вклада в величину

интеграла

j exp I —

• I

dpa.

 

101

Таким образом,

 

„2 , „2 , 2

оо со оо

рх + ри » + p z

1

г

г

Г

2 т А Г

 

— = J

j

j e

dpxdpydpz

=

 

—со —во •

 

 

 

 

 

 

-2

 

 

р?

 

 

 

 

 

оо

_^

 

 

оо

Л*,

 

СО

l'y

 

 

л

 

У

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р

2mAr

 

Л

2mkT

 

2 т £ Г

 

=

\

е

 

dp^

e

 

dp y

\ e

dp2 =

 

 

 

 

I

f

2mkT

dpx)

=

(2кткТ)

 

 

 

 

= I

\

e

 

[см. (П.1) J.

 

 

 

 

 

 

 

 

 

 

После подстановки значения 5 в выражения (IV.21) и (IV. 22)

найдем:

 

 

 

 

 

 

 

3

 

РІ+РІ+РІ

 

 

Р(Р*.

Ру,

Pz)

= (2*mkT)

2

 

2mkT

(IV.23)

 

е

 

;

 

 

 

 

 

 

2

 

2mkT

 

dw

х,

р у ,

pz)

= (2r.mkT)

 

е

 

dpxdpydpzl

(IV.24)

что представляет окончательные формулы распределения по состав­

ляющим

импульса

молекул.

 

 

 

 

 

 

 

Заметим, что выражение для плотности распределения вероятнос­

тей р х,

ру,

рг) является произведением трех сомножителей,

каждый

из которых зависит лишь от одной составляющей

импульса:

 

где

 

 

 

Р ІРх, Ру. Pz) =

Р (Рх) Р (Ру)

Р (Pz)>

 

 

(IV.25)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2mkT

 

 

 

 

 

 

(

р а

) = (2-кткТ)

е

(<*=*, у,

z).

 

(IV.26)

Это

говорит

о том, что распределения по составляющим

импульса

Рх,

Ру, Рг

независимы. Вероятность

dw(px)

того,

что составляющая

импульса

по оси х

имеет значение в интервале от р х до рх

+

dpx не

зависит от того, каковы значения других составляющих и определяет--

ся

выражением

 

 

 

 

 

 

 

 

 

 

 

Р2х

 

 

 

 

 

 

 

2

2mkT

 

 

 

 

 

dw (рх) = (2nmkT)

е

dpx.

(IV.27)

 

Из распределения (IV. 24) по импульсам можно получить распре­

деление

по

скоростям,

если

сделать

замену

переменных

ра = mva

= X,

у,

z):

 

 

 

 

 

 

 

dw (ѵх, ѵу, vg)

=\~çj-j-f)

e

 

dvxdvydvz.

(IV.28)

102

Функция распределения по составляющим скорости имеет вид

I

m

2kT

 

f(vx, Vy, Vz)

, -

(IV.29)

Формула (IV.29) представляет запись распределения

Максвелла.

 

Вероятность dw

х)

того,

 

что

составляющая

скорости

молекулы

вдоль оси X имеет значение в интервале от ѵх до

ѵх

+ dvx,

есть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

„ . . 2

 

 

 

 

 

 

 

 

 

 

 

dw (

 

/

от

 

2kT

 

 

(IV.30)

 

 

 

 

 

 

 

~

\

2ъкТ

du.

 

 

 

 

 

 

 

 

 

 

^

 

 

 

 

 

Аналогичные выражения

можно записать для составляющих ѵу и

ѵг.

 

Отметим

это

является

весьма

существенным

замечанием),

что

в

любой

реальной

систе­

 

 

 

 

 

 

 

 

ме

(газ, жидкость,

поверх­

 

 

 

 

 

 

 

 

ностный

слой

 

вещества)

 

 

 

 

 

 

 

 

распределение

по

скоро­

 

 

 

V'

 

 

стям центров инерции

мо­

 

 

 

 

 

лекул

представляет собой

 

 

 

 

 

распределение

Максвелла.

 

 

 

 

 

Такой результат

объясня­

 

 

 

 

 

 

 

 

ется тем,

что

в

функции

 

 

 

 

 

 

 

 

Гамильтона

системы всег­

 

 

 

 

 

 

 

 

да

выделяются

 

слагаемые

 

 

 

 

 

 

 

 

вида РІІ/2ШІ

индекс

 

 

 

0

 

 

Ѵх

частицы;

а

=

х,

 

у,

г) и

 

 

 

 

 

функция

Гамильтона

мо­

Рис.

17.

Распределение

по составляющей

жет

быть

представлена

в

 

 

 

скорости

ѵх

 

 

форме (IV. 18).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Распределение по составляющим скорости имеет вид нормального

распределения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2kT

 

 

(IV.31)

 

 

 

 

 

f K )

=

2nkT

 

 

 

(et = X,

y,

2 ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция f{va)

четная: f(va)

=

/(—va).

Функция

имеет максимум при

Va =

0,

причем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/(0)

=

2тЖ

 

 

(IV.32)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с повышением температуры высота максимума понижается. При і>а->- ->±оо кривая f(va) асимптотически приближается к нулю. Площадь, ограниченная кривой f (ѵа) и осью ѵа, равна единице, согласно условию

00

 

 

 

нормировки: jf{va)dva=

1. С

повышением температуры кривая ста-

00

 

 

 

новится более пологой

(рис.

17),

т. е. увеличивается вероятность

состояний с большими

значениями

| о а | .

103

Определим вероятность того, что модуль скорости молекулы

V =

~Ѵ ѵ2х + v2y + ѵ\

 

имеет

значение

в интервале от

ѵ до

ѵ +

dv.

При

этом

направление вектора

скорости может быть любым. Иначе

говоря,

 

мы

должны

определить

вероятность того,

что

конец

вектора ѵ в пространстве скоростей

ѵх,

ѵу,

ѵ2 лежит в шаровом

слое

радиуса

ѵ и толщины dv. Перейдем к сферическим координатам в рас­

сматриваемом

пространстве:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dvxdvydvz

 

д(ѵх,

 

Ѵу,

 

vz)

dvd%dy =

vi&\r\%dvdU<ç.

 

(IV.33)

 

 

 

=—\*[

 

;

 

*'

 

 

 

 

 

 

 

д(ѵ,

 

ѳ,

<р)

 

 

 

 

 

 

 

 

Согласно

(IV.28) и (ІѴ.ЗЗ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dw (v,

Ѳ,

<p) =

/I —m- —

\I

2

e

2 k T

i/2

sin ѲгіиШср,

 

(IV.34)

 

 

 

 

 

 

 

\

2шт

J

 

 

 

 

 

 

 

 

где

dw

(v,

Ѳ,ф) — вероятность

того,

что вектор скорости имеет мо­

дуль в

интервале от

t» до

v +

dv и заданное направление (углы в

интервале от Ѳ до Ѳ +

d9 и от ф до ф +

dtp)*. Для определения dw (v) —

вероятности заданного значения модуля скорости независимо от на­ правления вектора скорости следует выражение (IV.34) проинтегри­ ровать по всем возможным значениям Ѳ и ф:

3 mu'

dw (v) = f k

T

\ e

 

ѵЧѵ I dy

sin ѲгіѲ =

 

 

 

 

о

о

 

 

3

_

то'

 

= 4

л

Й г ) 2

е

2 k T ѵ Ч ѵ -

( І Ѵ - 3 5 )

Результат (ІѴ.35) мы могли бы получить просто умножением функции распределения (ІѴ.29) на объем сферического слоя 4nv2dv, поскольку все состояния с заданным модулем скорости имеют одну и ту же энер­ гию е п о с т и равновероятны (изображающие точки частицы в простран­ стве ѵх, ѵу, vz расположены на сфере). Функция распределения по модулю скорости дается выражением

 

 

 

 

3

_

тѵг

 

 

 

 

,(о)==4я(іаг) 2е

 

2kTѵК

(І Ѵ - 36)

 

Возможные

значения модуля

скорости

определены интервалом

(О,

оо), причем

/(0) = / ( о о ) = 0.

Кривая

f{v)

асимметричная, с мак-

 

* Выражение

(IV.34)

свидетельствует о

том,

что распределения

молекул

по

модулю скорости v и

направлению

вектора

скорости независимы.

Так как

энергия поступательного движения молекулы не зависит от направления век­ тора скорости, все направления вектора скорости являются равновероятными, распределение по составляющим Ѳ и <р вектора скорости беспорядочное.

104

симумом (рис. 18). Значение ѵ*, отвечающее

максимуму

функции

f(v)

 

(наиболее вероятное значение модуля скорости), находим из усло­

вия

djjv)

0. Дифференцирование выражения (IV.36) дает

dv

 

 

 

 

 

 

 

 

 

 

 

 

2кТ

 

mV

 

 

 

 

 

4тс 2nkT )

2ѵ*

о,

 

 

 

 

kT

 

откуда

 

 

 

 

 

 

 

 

 

 

 

2kT

 

 

( I V . 37)

 

 

 

 

 

 

 

 

 

Величина

ѵ* зависит

от массы

молекулы m

и температуры систе­

мы

 

Т. При

повышении

температуры

наиболее

вероятное

значение

Рис. 18. Распределение по модулю скорости

модуля скорости возрастает. Наиболее вероятные значения ѵ* двух разных молекул при заданной температуре относятся друг к другу как

 

 

(IV.38)

где mlv. пгг — массы молекул.

 

определяе­

Если в правую часть (IV.35) подставить значение ѵ*,

мое равенством (IV.37), то выражение для вероятности

принимает

более компактную форму:

 

 

dw (ѵ)

ѵЧѵ.

(ІѴ.39)

Из распределения (ІѴ.35) по модулю скорости легко получить рас­ пределение по энергии поступательного движения молекул. Произве-

105

дем в выражении (IV.35) замену переменных согласно соотношениям:

 

mi>2

/ 2е \ 2

1

/ 2 \ 2

s

2

(IV.40)

ê

= —

; о = —

; dy = —

de.

 

2

V m /

2

V m

 

 

 

Получим:

 

 

 

e_

_1_

 

 

 

2

W

2

 

 

dw (e)

=

e

s

de;

(IV.41)

 

 

1/ 71 (kT)3

E_

 

 

 

 

 

2

J

_

 

 

 

kT

2

 

 

/00

=

- =

e

e

.

(IV.42)

Кте (feT)3

§3. Средние значения некоторых функций скорости

поступательного движения частицы

Средние будем определять согласно общим формулам (1.32) и (1.33):

 

X =

J xf (х) dx;

А =

J Л (х) f (*) dx,

где

/ (х) — плотность

распределения

вероятностей для величины х;

А(х)

— некоторая однозначная

функция от х.

Найдем среднее значение компоненты скорости в положительном

направлении

оси

х.

Воспользовавшись

 

распределением

(IV.30) и

учтя только значения

ѵх > 0,

получим

 

 

 

 

 

 

 

 

 

-L

to

 

mvl

 

 

_L

ce

 

-L

 

 

 

 

2

 

х

 

 

2

 

2

-

(

m

 

\

С

vxe

2fcr

( m

\

kT'

Ç .

dt

kT

vx =

\

 

—-1

l

dvx

= I ——

e~l

2nm

x

2ккТ

/

J

 

.

V 2nkT

j

m

J

 

Величина vx характеризует среднее расстояние, которое проходит мо­ лекула за единицу времени в положительном направлении оси х. Такое же расстояние она проходит в отрицательном направлении оси. При усреднении по обоим направлениям получим ѵх = 0.

Среднее значение модуля скорости молекулы вычислим, учтя рас­ пределение (IV.39):

_ *

4

f

ѵ *2

V =

vf (v) dv — ——

\ v9e

dv,

J

Vu

o*3 J

 

2kT

где v* =УV" m . Интеграл в правой части выражения представ­ ляет интеграл Пуассона типа

со

 

/3 = \ х3е

~~ dx

о

2d*

 

106

(см. Приложение I), где а = 1 /ѵ*. Следовательно*,

4

1

2

. „ * 4 =

— о* 1,13Ü*.

(IV.43)

Среднее значение квадрата составляющей скорости vi находим, используя распределение (IV.30):

оо

2 оо

£ .

ОО —CO

Заметив, что интеграл в правой части есть

 

/ , =

J

x4~"'

dx =

~Л/

-^-

 

 

 

—ОО

 

 

 

'

 

 

при а =

, получим

=

 

. Очевидно,

 

 

 

 

 

 

 

feT

 

 

 

 

ѵ\

=

y 2

= [ J 2

=

^m_ .

 

(IV.44)

Среднее значение квадрата модуля скорости равно

 

 

~ф =

vi +

и?; + ѵ\ =

 

(IV.45)

 

 

 

 

"

m

 

 

 

Средняя квадратичная скорость имеет значение

 

 

У~&

=

у^Ш-

=

у^JL

и*.

(IV.46)

Относительное положение величин и*, о и \^ ѵ2, показано на рис. 18. Вычислим среднее значение кинетической энергии поступательно­

го движения частицы:

 

 

m 1 о

о

Л

тѵ2

3

 

 

 

«„ост =

— \fx

+ 0 ; + ѵ\) =

= т

kT.

(ІѴ.47)

 

 

(*

— од:2

 

 

 

 

 

* Интегралы типа

I

х"е

dx

удобно

вычислять

также с помощью

 

 

о

 

 

 

 

 

 

 

формул для Г-функции

(см. Приложение II) . Так,

 

 

оо

_ . И І _

 

 

ОО

 

 

 

 

 

\

ѵ*е

dv = v** \

e-'tdt=—

 

Г (2)

о

 

 

 

о

 

 

 

 

 

107

Величину е п о с т можно рассчитать и другими способами, например, следующим:

tn Г

епост=-Г- v*f(v)dv,

1 О

где функция f(v) дается выражением (IV.36), или же согласно формуле

œ

£пост = j E / ( e ) de,

О

где / ( е ) — функция, определяемая равенством (IV.42).

Значение энергии, приходящейся в среднем на одну степень сво­

боды поступательного

движения

частицы,

равно

_

кт

(а=х,у,г).

 

 

Епоста = —

(IV.48)

Таким образом, энергия в среднем равномерно распределяется по степеням свободы поступательного движения, и на каждую степень свободы приходится средняя энергия kT/2. Вклад поступательного движения в среднюю энергию моля газа составляет

£пост = < Г П 0 С Т = - | - Я 7 \

(IV.49)

где N0— число Авогадро, R — газовая постоянная. Вклад в молярную теплоемкость определится величиной

Споет = ^ f - T = - f R-

(ІѴ.БО)

§ 4. Число ударов молекул о единицу поверхности. Давление идеального газа

Найдем число молекул, ударяющихся за единицу времени о еди­ ницу поверхности. Обозначим эту величину (величину потока) через W.

Допустим, что площадка dS расположена перпендикулярно оси х, и молекулы налетают на эту плещадку слева. Выделим сначала мо­ лекулы, которые имеют составляющую скорости вдоль оси х в интер­ вале от ѵх до ѵх + dvx. Число таких молекул в единице объема равно

dn(vx)=—

dw(vx)=

\ - ^ r j е

dvx,

(IV.51)

где N— ебшее число

частиц в

сбъеме V; для вероятности

dw(vx)

использовано выражение (IV.30). За единицу времени до площадки dS долетят все молекулы с заданным значением ѵх, расположенные в

108

объеме vxdS, т. е. dn (vx)vxdS молекул*.Число ударов молекул с за­ данным значением составляющей ѵх о единицу поверхности в едини­ цу времени, следовательно, есть

Ni

m \

~1кт

(IV.52)

dW (ѵх) = vxdn (vx) = —

[

dvx.

Полное число ударов молекул о единицу поверхности за единицу времени получим, проинтегрировав выражение (IV.52) по всем поло­ жительным значениям ѵх (отрицательные значения составляющей отве­ чают движению молекул не к площадке, а от нее):

(IV.53)

о

С помощью выражения (IV.52) легко рассчитать давление идеаль­ ного газа. Сила, с которой идеальный газ действует на стенку сосуда, определяется исключительно упругими ударами молекул газа о стен­ ку, и эту силу можно приравнять изменению количества движений частиц в единицу времени (сила X время = изменение количества движения). Молекула, имеющая составляющую скорости ѵх (пола­ гаем, что стенка перпендикулярна оси х), при упругом ударе о стенку меняет знак этой составляющей на обратный и отдает стенке количест­ во движения 2тѵх. Сила, действующая на. единицу поверхности стенки (давление), равна изменению количества движения стенки в единицу времени, отнесенному к единице поверхности. Число ударов с задан­ ным значением составляющей ѵх за единицу времени о единицу по­ верхности определяется формулой (IV.52), так что

* То, что составляющие скорости ѵу

и ѵг

не равны нулю(

можно не

прини­

мать во

внимание: поскольку значения ѵу

и — ѵ у і ѵ г и — ѵ г

равновероятны,

столько

же молекул будет в среднем входить в выделенный объем v^dS

через

боковые

стенки, сколько и выходить из него. Число ударов о поверхность опре­

деляется

лишь составляющей скорости

ѵх,

перпендикулярной

к поверхности.

109

Соседние файлы в папке книги из ГПНТБ