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

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

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

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

Добавлен: 10.04.2024

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

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

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

Я-теорема Больцмана подверглась ожесточенной критике со сторо­

ны многих

ученых, в частности, за ее кажущуюся несовместимость

с законами

механики.

Одно из возражений, высказанное Лошмидтом, состояло в следую­

щем.

Допустим, газ,

находившийся

в начальном состоянии

Л0 , через

время

t

оказался

в

состоянии

At.

В

соответствии

с Я-теоремой

Н( - < Я 0

(S,

S0 ).

Если обратить направления скоростей, то, согласно

законам механики,

за время t

система

должна

вернуться в

исходное

состояние А0,

которому

отвечает

значение Я

= Я 0 .

Но

Я-теорема

требует, чтобы Я 0 <^ Ht\

Можно не ставить вопроса о том, как обратить

скорости молекул; важно лишь, что с механической точки зрения подобное состояние системы не исключено.

Другое возражение (парадокс Цермело) основывалось на воз­ вратной теореме Пуанкаре. Если система в момент tx находится в неравновесном состоянии Ах, то согласно Я-теореме в системе проис­ ходят процессы, приближающие ее к состоянию равновесия и сопро­ вождающиеся ростом энтропии. В момент времени t2 > tt система будет находиться в некотором состоянии А2 , которому отвечает энтро­ пия S2 > Sx (Sx и S2 — значения энтропии в моменты времени tx и t2 соответственно). Однако, как показывает возвратная теорема Пуан­

каре,

в

некоторый

момент времени t3> t2

система должна

оказаться

в

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

с состоянием Аъ

так что

Ss

~

Sx.

Очевидно,

переход от состояния

А2

к состоянию

As

должен

сопровождаться уменьшением энтропии, что противоречит Я-теореме.

Как согласовать обратимость механических процессов на молеку­ лярном уровне с наблюдаемой на опыте необратимостью макроскопи­ ческих процессов, было указано самим Больцманом в-его более позд­

них работах: разрешение кажущихся противоречий

состоит в том,

что Я-теорему следует понимать как статистическую

закономерность.

Многое для уточнения смысла Я-теоремы дали также работы П. Эренфеста и Т. Эренфест.

Статистический характер закона возрастания энтропии

вытекает

из самого определения энтропии (III.70), связывающего эту

функцию

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

макроскопического параметра X*

точностью до интервала

АХ),

она почти наверняка придет к этому

состоянию; если же система

уже

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

Статистический подход позволяет объяснить приведенные выше парадоксы. Действительно, система может вернуться в исходное

81


неравновесное состояние, как того требует возвратная теорема Пуан­ каре, и такой процесс будет сопровождаться уменьшением энтропии. Однако вопрос в том, насколько вероятен подобный процесс. Неболь­ шие флуктуации будут происходить довольно часто. Значительным же отклонениям от равновесия отвечает фазовый объем, составляю­ щий ничтожно малую долю всего объема энергетического слоя; изобра­ жающая точка может попасть в данную область лишь по прошествии огромного промежутка времени и будет находиться в этом объеме ничтожно малое время. Ниже в качестве примера приводятся данные, показывающие, как быстро возрастает среднее время возвращения к исходному неравновесному состоянию с ростом относительных

отклонений от

равновесия. Рассматриваются

флуктуации плотности

в объеме 1 см3,

составляющем часть большого объема газа при нормаль­

ных условиях.

 

 

 

 

 

 

Относительное

отклонение

Время

возвращения

 

плотности

от

равновесной

 

 

 

 

2

10-ю

 

4- Ю - 3 сек

 

3

10-ю

 

1

»

 

4

10-ю

 

21

мин

 

5

10-ю

 

5

мес.

 

6

10-ю

 

3-Ю4 лет

 

7

Ю-ю

 

2 - Ю 1 0 »

Действительно, прав был Больцман, ответив на возражение Цермело о необходимости возвращения систем в исходное, далекое от равновесного, состояние: «Долго же Вам придется ждать!».

Легко понять, что вероятность зафиксировать в макроскопическом опыте самопроизвольное сжатие газа практически равна нулю. Ве­ роятность того, что молекулы идеального газа соберутся все в одной половине сосуда при N = 102 3 есть 2~і023 ~ Ю - 3 ' 1 0 " . Эта вероятность не нулевая, как того и требует механика, но столь мала, что событие следует практически отнести к невозможным.

Подведем итог сказанному. Итак, переход системы из равновесного в неравновесное состояние допустим, но вероятность значительных отклонений от равновесия, связанных с заметным уменьшением энтро­ пии изолированной системы, практически нулевая. В то же время небольшие отклонения от равновесия происходят очень часто; в ка­ кие-то моменты времени энтропия системы уменьшается. Статистичес­ кая интерпретация энтропии, следовательно, раскрывает смысл вто­ рого начала термодинамики и указывает границы его применимости: закон возрастания энтропии в изолированной системе (и постоянства энтропии при равновесии) справедлив лишь, если пренебречь флуктуационными процессами.

§ 9. Каноническое распределение

Каноническое распределение Гиббса — статистическое распределе­ ние для систем, имеющих заданное число частиц N, заданный объем V и способных обмениваться энергией с окружением. В общем случае, если помимо потенциала, создаваемого стенками сосуда, имеются

82


другие

внешние силовые поля, задается набор

внешних координат

аъ..., as,

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

энер­

гии системы не наложено никаких ограничений, и в этом отличие сис­

темы канонического ансамбля от системы микроканонического

ансам­

бля. Система канонического ансамбля находится

в жесткой,

непро­

ницаемой для частиц оболочке, но эта оболочка

неадиабатическая,

теплопроводящая, что и делает возможным обмен энергией

между

системой и окружением. Энергетическое взаимодействие системы с окружением предполагается, однако, достаточно слабым, чтобы сис­ тему можно было считать статистически независимой (энергия взаи­ модействия пренебрежимо мала по сравнению с полной энергией сис­ темы). Предполагается, что макроскопическое состояние окружения практически неизменно, и ставится задача найти плотность распреде­ ления вероятностей в фазовом пространстве для различных состояний системы. Иначе говоря, требуется решить вопрос о виде функциональ­ ной зависимости р(Я). Рассмотрим наиболее простой вариант вывода.

Допустим, что окружение интересующей нас системы 1 составляет

очень

большая система 2. Системы 1 и 2 взаимодействуют слабо, так

что в выражении

для полной энергии совокупности «система 1 + сис­

тема

2» можем

пренебречь

членом,

связанным с

взаимодействием:

 

 

Я ( р , q)=Hl(p{l),

<7( 1 ) ) +

Н2{2),

<7 ( 2 ) ),

(III.101)

где величины без индекса относятся к совокупности; функция Гамиль­ тона Я х для системы 1 зависит только от микропараметров рх и qx , относящихся к данной системе; аналогичное справедливо для системы 2. Будем считать, что совокупность систем 1 и 2 в целом изолирована:

Я ( р , с ) = £ = const.

Если системы 1 и 2 находятся в равновесии внутри себя и между собой, то совокупность «1 + 2» является в целом равновесной и к ней можно применить выражение ( I I 1.39):

 

Pi = pi(tf0;

Р = Р

);

Р

= Р ( Я ) =

Р

1

+

Н

) .

(III.102)

 

 

2

2

2

 

 

 

 

 

2

 

Поскольку системы

1 и 2, по условию,

квазинезависимы, норми­

рованные

плотности вероятности

р = N!hfNp

 

 

должны

удовлетворять

равенству

( I I 1.67):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

!Рі № ) н № ) = Р ( Я )

или

1пр(Я) = 1пР і (ЯО + ln Р 2 2 ) .

( I I I . 103)

Дифференцирование обеих частей уравнения ( I I I . 103) дает

d In р dH = d In pi dH,

+ d ln p.'2 dH.

dH

~dîû

83


Учитывая, что dH = dHx + dH2 — 0 (по условию изоляции со­

вокупности систем 1 + 2) и dHx

— —dH2,

получаем

fL^ÎL^Éhh^^

(III.104)

dH\

dH%

 

где a — некоторый параметр, общий для систем 1 и 2 и не зависящий

от микроскопических

переменных. Равенство ( I I I . 104) должно вы­

полняться при любых

значениях Нх. Используем

теперь

предполо­

жение о том, что система 2 является очень большой.

Обмен

энергией

с системой 1 изменит состояние системы 2, только очень мало (при бес­ конечно большом размере системы 2 изменения ее состояния за счет обмена энергией с системой 1 конечного размера будут бесконечно ма­ лыми). Систему 2 можно назвать термостатом. Так как состояние систе­ мы 2 практически неизменно, параметр а, относящийся к ней, следует считать постоянным. Получается так, что система 2 (окружение) как бы задает значение параметра a для интересующей нас системы 1, и мы рас­ сматриваем вероятности всевозможных состояний систем 1 при опре­

деленном, заданном значении a

(кроме того, очевидно, при заданных

V и N). Интегрируя выражение

( I I I .

104) при a = const, находим вид

функциональной зависимости Рі(#і):

 

 

 

 

7і = Л е о Я і

 

( I I I .

где А — постоянная для заданных

условий.

 

Поскольку в дальнейшем нас будут интересовать только состоя­

ния

системы 1, а наличие термостата

будет учитываться

через па­

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

 

 

q) =Ae"iF*

ч ) ,

( I I I . 10

где все переменные относятся к системе; параметр a задается

окружени­

ем

(термостатом).

 

 

 

 

Вероятность того, что микропараметры системы имеют

значения

в интервале от р до р + dp для импульсов и от q до q + dq для коор­ динат, определим согласно соотношению ( I I 1.68):

dw (р, q) = — , Ае"и ( р ' "> dp dq.

Nlh

По условию нормировки

JJе«я (Р. ч) dpdq=l,

 

1

( I I I . 1

где интегрирование проводится повеем состояниям, совместимым с за­ данными значениями V, N (и при любых значениях энергии). Если возможные значения энергии сверху не ограничены, что имеет место

для всех

систем, то условием сходимости интеграла в левой части

( I I I . 107)

является

84


Положим

 

a < 0 * .

 

( I I I . 10

 

 

 

 

 

 

 

а = - 4 " .

где Ѳ > 0,

( I I I . 10

 

 

о

 

 

 

и вместо

( I I 1.105) запишем:

 

 

 

 

 

 

Я (р, <7)

 

 

 

Р ( р ) 9 ) = Л е

 

;

( Ш Л Ю )

 

 

 

 

И (р, д)

 

 

 

РІР. Я)=—*~т

Ае

9 .

(III.111)

 

 

NW

 

 

 

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

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

ностей в фазовом пространстве есть

экспоненциально убывающая

функция

Я . Формула ( I I I . I l l ) представляет запись

канонического

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

Гиббса.

 

 

 

Величину Ѳ Гиббс назвал модулем

канонического

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

Условием равновесия между системами, находящимися в энергети­

ческом контакте, является

равенство величин Ѳ для этих систем: если

Ü 2 E £ =

_ _ L и

i ! £ £ _ _ _ _ L f

( і 1 і . 1 1 2 )

то при равновесии, согласно ( I I I . 104)

 

 

Ѳ 1 = Ѳ 2

=

Ѳ.

(III.113)

Основываясь на равенстве

( I I I . 113),

можем утверждать, что параметр

Ѳ имеет смысл температуры

(вспомним, что в основе определения тем­

пературы лежит принцип равенства при равновесии температуры всех тел, находящихся в энергетическом контакте). Величину Ѳ называют статистической температурой. Эта величина положительна и согласно ( I I I . ПО) должна иметь размерность энергии, чтобы показатель степени Я/Ѳ был безразмерным. Позднее мы покажем, что

 

Ѳ = М \

 

(III.114)

где Т — абсолютная термодинамическая

температура,

k — постоян­

ная Больцмана. Таким образом, задание

параметра Ѳ равносильно за­

данию температуры системы.

 

 

Теперь мы можем уточнить задание

параметров,

определяющих

макроскопическое

состояние системы канонического

ансамбля: это

параметры T, N,

V. Рассматривается система в жесткой, непроницае­

мой для частиц теплопроводящей оболочке, помещенная в термостат.

В общем случае, когда,

помимо давления, на систему действуют

другие внешние силы

и система многокомпонентная, задается набор

параметров

 

 

T,

Ni

NK, at, ... , as.

* При a>0 функция eaH является непрерывно возрастающей функцией энергии и, если значение энергии не ограничено сверху, интеграл $$ еdpdq

расходится, условия нормировки выполнены быть не могут. Положительному значению а отвечал бы случай, когда состояния с большей энергией имеют боль­ шую вероятность (см. гл. V I I , § 7).

85