Файл: Смирнова Н.А. Методы статистической термодинамики в физической химии учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 188

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

верхности dS (вектор dS по величине равен площади dS элементарной поверхности; направление его совпадает с направлением положитель­ ной нормали к площадке). Так как число фазовых точек ансамбля постоянно (системы ансамбля не возникают и не исчезают), то убыль фазовых точек в элементе объема Д Г равна потоку фазовых точек через замкнутую поверхность AS, ограничивающую данный элемент объема:

д С РаТ = (В РѴ dS.

(III.16>

dt

(ДГ) (AS)

Согласно теореме Гаусса — Остроградского поток вектора через зам­ кнутую поверхность равен объемному интегралу от дивергенции (расхождения) вектора*, так что

ф Р VdS=

J div(P V) dT.

( I I I . 19

(AS) (ДГ)

Поскольку операции дифференцирования по времени и интегрирова­ ния по объему в левой части уравнения (III.16) независимы, порядок их может быть изменен. Учитывая это и равенство ( I I I . 19), можем и» уравнения ( I I I . 16) получить следующее:

— ^ -^- dT = J" div (Р V) dT.

(III.20)

(ДГ) (ДГ)

Уравнение (III.20) выполняется для произвольного объема Д Г , по которому проводится интегрирование. Следовательно, условием спра­ ведливости уравнения является равенство подынтегральных выраже­ ний в обеих частях:

( ^ А

= - d i v ( P K ) .

(111.21)

\ Ot

Jp, q

 

Производная ( —

)

характеризует изменение плотности фазовых

\ àt

Jp,

q

точек в единицу времени в окрестности некоторой фиксированной точки фазового пространства. Уравнение (III.21) есть известное в

гидродинамике уравнение неразрывности, записанное в

применении

к движению фазовой жидкости. Это уравнение является

следствием

непрерывности движения и постоянства числа фазовых точек ансамбля.

Используя формулы векторного

анализа, находим:

 

 

 

 

div(P V) = P d i v V+

KgradP,

 

(111.22)

*

Дивергенция трехмерного вектора а =

х, ау,

 

az) в системе

координат

X, у,

г есть, по

определению,

 

 

 

 

 

 

 

 

да г

даѵ

да,

.

 

,ТІІ

,„ч

 

 

d i v c = - _ L

ду

dz

 

(III.17)

 

 

дх

 

 

 

 

В общем случае r-мерного вектора а (а\

аг),

определенного

в системе

координат хі,

хг.

 

 

 

 

 

 

 

 

S ьт^да;aas

 

 

 

( І І І Л 8 )

50


где

 

 

 

(=i

V

d<7i

dpi

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

дР

дР

'

\

 

 

 

 

 

 

I

 

 

 

 

 

S b a - i q i + ^ : p i

) -

( І І І - 2 4 )

 

 

 

 

1=1

 

 

 

 

 

 

Согласно уравнениям (11.28)

при движении по фазовым

траекториям

 

 

dg,

_

_

др;

_ д*Н

 

 

(111.25)

 

 

dqt

 

 

dpi

dpidqi

 

 

 

 

 

 

 

 

 

так что

 

 

 

 

 

 

 

 

 

 

 

 

dqt +

dpi

_ Q .

d i v K

=

o.

 

(111.26)

Из

уравнений

(III.21) — (III.26)

следует:

 

 

 

 

 

 

дР \

F

I

дР •

дР

\

 

 

 

 

V 1

 

 

где

изменения

переменных

pt

и

<ft отвечают

движению

по фазовой

траектории (по любой из фазовых траекторий систем ансамбля). Уравнение (III.27)—одна из форм аналитической записи теоремы

Лиувилля. Следует подчеркнуть, что при выводе

его были

учтены

уравнения движения. Если принять во внимание

(11.38) и (11.28),

уравнение ( I I 1.27) можно

 

записать,

используя скобки Пуассона:

 

 

 

 

дР

 

 

 

 

 

 

 

 

 

= -{Р,Н].

 

(111.28)

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

( I I 1.27).

Полный

дифференциал функции

Р

(р, q, t)

 

 

 

 

,

àP

 

 

 

I

dP

dP

\

 

dP = (

 

 

 

.S.(*r*, + sr*')î

 

 

Ot

/р.

q

 

полная производная

по

времени

 

 

 

 

dP

I

дР\

 

F

дР '

дР .

 

 

 

\ i I

 

(III.29)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(=1

 

 

 

 

dP

Производная — . характеризует скорость изменения плотности

dt

фазовых точек при движении их по фазовым траекториям, т. е. изме­ нение плотности в непосредственной окрестности произвольно выбран­ ной движущейся фазовой точки. Из уравнений (III.27) и (III.29)

51


вытекает следующее:

~ - = о ,

(111.30)

at

утверждающее, что плотность фазовых точек при движении их по фазовым траекториям остается постоянной. Движение фазовых точек аналогично движению несжимаемой жидкости. Уравнение (III.30) является записью принципа сохранения плотности «фазовой жидкости» и наряду с уравнением ( I I 1.27) выражает сущность теоремы Лиувилля .

Найдем изменение во времени объема Д Г ,

занимаемого AL

точка­

ми, при движении

точек

по фазовым траекториям.

Полагаем, что

Д Г — очень малый

объем,

так что

внутри

него

плотность фазовых

точек приближенно можем считать постоянной: Д / ,

=

РД Г. Мы фикси­

руем границы фазового

объема Д Г таким образом, что в этом

объеме

находится заданное число точек Д L, и следим за движением этого объе­

ма. Очевидно,

d\L

 

Р d Д Г +

 

dP

 

 

 

 

 

=

Д Г

= 0

 

 

 

 

dt

 

dt

 

dt

 

 

 

 

и, в силу условия (III.30),

 

 

 

 

 

 

 

 

 

 

- * - А Л = 0 ,

 

 

 

 

( Ш . 3 1 >

 

 

 

dt

 

 

 

 

 

 

Таким образом, если Д L точек при своем движении

по фазовым

траек­

ториям перешли из

объема

Д Г в объем Д Г '

(рис. 8), то

 

дг = д г ,

хотя по форме эти элементы объема могут отличаться. Иначе говоря, объем, занимаемый фазовыми точками совокупности изолированных систем, при изменении состояния систем остается постоянным и может изменяться только по форме. Этот вывод можно распространить на объем любого размера, движущийся в энергетическом слое, так как всякий объем можно представить как сумму малых объемов. Итак, всякий фазовый объем, занятый заданным числом фазовых точек, при своем движении в энергетическом слое соответственно изменению состояния систем ансамбля остается неизменным по величине. Данная формулировка теоремы Лиувилля может быть названа принципом сохранения фазового объема.

Вывод уравнений (III.27), ( I I I . 30), (III.31), выражающих сущность теоремы Лиувилля, был основан на использовании канонических уравнений движения при описании поведения ансамбля изолирован­ ных систем. Выведенные соотношения справедливы только при ис­ пользовании канонических переменных, т. е. при описании движения изображающих точек ансамбля в пространстве обобщенных координат

и импульсов. Для пространства переменных qt и qt аналогичные об­ щие соотношения, в частности принцип сохранения фазового объема, выведены быть не могут. Этим объясняется то предпочтение, которое в статистической физике оказывают каноническим переменным.

52


Поскольку плотность фазовых точек связана с плотностью распре­ деления вероятностей соотношением ( I I I . 8 ) (P = pL, где L = const),

то теорема Лиувилля определяет

изменение р для произвольно выбран­

ной системы ансамбля. Вместо

уравнений ( I I I . 2 7 ) , ( I I I . 2 8 ) и ( I I I . 3 0 )

можем записать:

 

(•fL+s(*r;'+t''H

і=і

 

d t ] р Г - [ ? , Щ ;

(ІІІ.ЗЗ)

- ^ - = 0 .

(111.34)

Эти уравнения описывают изменение р при движении вдоль фазовой траектории данной системы. Дальнейшие утверждения о виде функцио­ нальной зависимости р(р, q) в большой степени основываются на уравнении ( I I I . 3 2 ) и вытекающих из него следствиях. Именно таким образом в статистической физике учитывают механические уравнения движения.

Теорема Лиувилля справедлива как для равновесных, так и для неравновесных ансамблей. Если ансамбль находится в статистическом

равновесии, то { ^ \

= 0, и уравнения

( I I I . 3 2 )

и ( I I I . 3 3 ) принимают

вид:

F

 

 

 

 

dp

• \

 

 

dp •

 

 

- ^ - Ч І + - Г - Р І ) = 0;

(III.35)

 

. dqt

apt

J

 

 

{ p , # } = 0 .

 

(III.36)

Равенства должны

выполняться

для

любого

энергетического слоя

при различных значениях р и Н. Отсюда следует, что в случае равно­ весного ансамбля плотность распределения вероятностей должна зависеть от р и q только через интегралы движения. Действительно, без ограничения общности можем предположить, что для равновес­

ного

ансамбля функция

р может

быть представлена

в форме

=

 

 

т

гл

Е

Ф І

=

Ф І ( Р >

а

) некоторая функция обобщенных

P

р(фі. •••> ф )>

-

 

 

 

координат

и

импульсов

 

(/ =

1,

т).

Полагаем,

что

величины

Ф І , ...,фо т

независимы. После подстановки в уравнение

( I I 1.35) находим:

 

 

F

p /

TJX

 

 

 

 

\

 

/

Т71

\

 

 

 

 

m

 

F

 

 

 

 

 

 

 

 

m

 

 

 

 

/=і

 

і=і

 

 

 

 

 

 

 

і=і

 

 

Так

как ~ -

Ф 0, то равенство при независимости величин ф г , ...,фт

53