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

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

.pdf
Скачиваний:
124
Добавлен:
22.10.2023
Размер:
15.13 Mб
Скачать

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

Для изучения статистического поведения термодинамических систем • Гиббс предложил метод ансамблей. По Гиббсу ансамбль — это доста­ точно представительный ряд микросостояний, соответствующих, т. е. совместимых, с данным макросостоянием. Иначе, ансамбль можно представить как очень большое число (в пределе стремящееся к бес­ конечности) аналогичных термодинамических систем, находящихся, однако, в различных стадиях развития, т. е. различных микросостоя­ ниях.

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

Микроканонический ансамбль соответствует изолированной сис­ теме, которая не может обмениваться с окружением ни веществом, ни

энергией. Она характеризуется постоянством v, U, N

(или рядом Nv

N2,

если система многокомпонентна).

 

 

Канонический ансамбль соответствует закрытой

изотермической

системе, которая характеризуется значениями у, Л; и Т. Такая система находится в термическом равновесии с окружающей средой и может обмениваться с ним энергией. Канонический ансамбль оказался наи­ более удобным для целей статистической термодинамики.

Оба перечисленных ансамбля вместе называют малыми каноничес­ кими ансамблями.

Большой канонический ансамбль соответствует открытой изотерми­ ческой системе, характеризуемой объемом v, температурой Т и хими­ ческим потенциалом \i. Он находится как в тепловом, так и материаль­ ном равновесии с окружающей средой и может обмениваться с ним . и энергией, и веществом.

Развитие во времени микроскопического динамического состояния каждого члена малого канонического ансамбля будет описываться траекторией его изображающей точки в 2УУ/-мерном фазовом простран­ стве. Если в том же самом пространстве помещены изображающие точки всех других членов ансамбля, то получится как бы облако движущихся

изображающих точек. Мгновенная

плотность

облака в данной точке

может быть охарактеризована функцией распределения

 

p(qi...qf

P l . . . P f t ) ,

'

(VI.13)

или в сокращенной записи р (q, р, t). Функция распределения норми­ руется следующим образом:

\...\dqi...dptp{q,

 

р,

0 =

1,

 

(VI.14)

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

долю

от

общего

числа си- >

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

на

единицу гиперобъема

фазового

пространства в области точки ql

...

pf в

момент t.

Подынтегральное

— 180 —

выражение в (VI. 14) — также доля систем, но в элементе гиперобъема dqr ... dpf. По-видимому, р выражает также вероятность, что изобра­ жающая точка наудачу взятой системы из ансамбля окажется в той или иной единице объема фазового гиперпространства.

§ 4. Теорема Лиувилля

Движение «облака» точек в фазовом пространстве (говорят также о жидкости вероятности) определяется законами механики, т. е. уравнениями Гамильтона (VI. 11), имеющими непрерывные решения. Отсюда следует, что соседние изображающие точки принадлежат сосед­ ним траекториям. Элементы гиперобъема, определенные содержанием

того или

иного числа

точек, деформируются

н е п р е р ы в н о

(в смысле

непрерывности

функций)

в течение

развития во

времени.

В этом параграфе будет показано, что величина

гиперобъема,

занима­

емого этими точками, тем не менее

о с т а е т с

я

п о с т о я н н о й .

А это эквивалентно утверждению, что распределение вероятности мож­ но изобразить в виде несжимаемой жидкости.

Как уже говорилось, фазовые траектории суть непрерывные кри­

вые — они

не могут

начинаться или обрываться. Другими словами,

изображающие точки

не

могут зарождаться или

уничтожаться, т. е.

«жидкость

вероятности»

неуничтожима. Поэтому внутри

каждого

с т а ц и о н а р н о г о

 

элемента гиперобъема Aw скорость увели­

чения числа изображающих точек должна в

точности

равняться

общему числу точек,

покидающих Асо. Если

беспредельно умень­

шать Асо, получим аналог уравнения течения неуничтожимой жид­ кости, т. е.

 

 

-g- + div(p^) = 0,

*

(VI. 15)

где v — вектор

скорости

течения.

 

 

В трехмерном пространстве член с дивергенцией

(расхождением

вектора) имеет

вид

 

 

 

 

div (pw) =

~ (pvx) + -~ (pvy) + ~ (pvz).

 

(VI. 16)

Основываясь на тождестве (VI. 16), запишем в развернутом виде обоб­ щение уравнения (VI. 15) для 2Nf измерений

 

Nf

 

+aiM= 0 -

 

 

 

(VLI7)

Дифференцируя

произведения в (VI. 17), получим

 

-2-+2

w

 

 

Применяя уравнения Гамильтона (VI. 11) и замечая, что для непре­ рывных функций порядок дифференцирования безразличен, убежда-

181 —

емся в равенстве нулю второй суммы в (VI. 18). В самом деле каждый член этой суммы принимает согласно (VI.5) следующий вид:

dpidqi ~ dqtdpi '

однако с разными знаками.

Далее, так как р рассматривается в качестве функции qit pt и t, остающиеся члены (VI. 18) представляют собой полную производную

по

времени,

т. е.

 

 

 

 

т

 

 

что

можно

также получить, деля

полный

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

р (q, р, t) на dt.

 

 

 

Следует

обратить внимание на

различие

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

dp /dt и dp /dt. Первая представляет собой изменение плотности вероят­

ности во времени по соседству со

с т а ц и о н а р н о й

точкой

qtpi фазового пространства, тогда

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

dp /dt

является соответствующим изменением по соседству с д в и ж у щ е й-

ся изображающей точкой, следующей своей траектории в соответствии

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

Учитывая уравнения (VI. 18) и (VI. 19), приходим к теореме Лиувилля (1838)

 

 

 

^ - = 0,

 

 

 

(VI.20)

утверждающей, что по соседству с движущейся в фазовом

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

изображающей точкой

плотность

изображающих

точек

остается

постоянной.

 

 

 

 

 

 

 

Поэтому, как утверждалось в начале этого параграфа,

элемент

гиперобъема Асо может с течением

времени менять

свои очертания,

но не расширяться и не сжиматься.

 

 

 

 

Рассмотрим простой пример, связанный с падением материальной

точки массы т в поле

земного

тяготения. Импульс

(количество дви­

жения)

и координата

точки

выразятся

следующими

функциями:

 

 

p = mgt,

q = ^gt*,

 

 

(VI.21)

где g — ускорение силы тяжести.

 

 

 

 

 

Для

определения фазовой траектории

исключаем из (VI.21) время

и получаем уравнение

энергии

 

 

 

 

 

 

 

Я =

 

 

 

 

(V I -22)

Следовательно, в этом случае фазовой траекторией является парабола (рис. VI.2). Пусть одновременно с точой / нас интересует движение трех

182 —

других точек, движущихся в том же поле тяготения, но отличающихся от точки / другими начальными значениями р0 и q0, как это записано в табл. 26.

 

 

 

 

Т а б л и ц а

2 6 .

 

<?0

р

я

 

точки

Ро

 

1

0

0

mgt

W

 

2

dpa

0

mgi+dpu

 

3

0

dq„

mgt

^g*2 + dq0

4

dpi,

dqu

mgt + dpQ

 

Из этих данных таблицы и рис. VI.2 видно, что четыре рассматрива­ емые точки, которые в' начале движения образовали на фазовой плос­

кости р — q прямоугльник

dpudq0,

 

с течением времени изменяют свое

 

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

расположение. При

 

этом,

однако, точки /

и 3,

а

также

 

2 и

4

будут

оставаться

все время

 

на одинаковой высоте и расстоя­

 

ние между ними dqu остается неиз­

 

менным. В то же время верхняя

 

сторона

четырехугольника

1342

 

движется

параллельно

 

оси q

 

быстрее,

чем

нижняя,

вследствие

 

наличия

у точек 2 и 4

начального

 

импульса dp0.

Поэтому прямоуголь­

Рис. VI.2. К теореме Лиувилля (Фа­

ный

четырехугольник

превращает­

зовая площадь для падающих мате­

ся с

течением

времени

в

паралле­

риальных точек остается постоянной).

лограмм. И тем не менее

площадь

 

фигуры, образуемой четырьмя точками на фазовой плоскости, и во время движения остается неизменной, как того и требует теорема Лиувилля.

§ 5. Эргоидная гипотеза

иосновные постулаты классической статистической механики

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

183 —

эргоидной гипотезой, впервые предложенной Больцманом и Максвел­ лом: изображающая точка изолированной системы (v, N, U постоянны) перед возвращением в исходное положение проходит все достижимые, т. е. совместимые с заданными условиями точки фазового простран­ ства.

Иными словами, через достаточно длительное время механическая система вернется в исходное состояние, пройдя все другие достигаемые состояния, причем достижимость в данном случае ограничена соблю­ дением закона сохранения энергии *:

U (или E) = H(q, р) = const.

(VI.23)

Эргоидная гипотеза, которую Максвелл называл «принципом непре­

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

гипотезой.

Более того, она нуждается, как показали П. Эренфест и Т.

Афанась-

ева-Эренфест, в уточнении, заменяющем «возвращение в исходное

положение»

«приближением в

конечном

счете

сколь

угодно близко

к каждой достижимой

точке фазового пространства».

Необходимость

уточнения

эргоидной

гипотезы

связана

также

и с

существованием

динамических систем, изображающие точки которых никогда не прони­ кают в достижимые области их фазового пространства. Такие системы можно назвать неэргоидными и привести в качестве примера планетную систему, в которой планеты будут, вероятно, всегда оставаться в плос­ кости эклиптики, хотя орбиты, расположенные, скажем, в перпенди­ кулярной плоскости, энергетически вполне «достижимы». Таким образом, эргоидная гипотеза или квази-эргоидная гипотеза (по Эрен­ фест ам) налагает известные ограничения на системы, подвергающиеся

изучению с точки зрения статистической механики, а именно

т а к и е

с и с т е м ы

д о л ж н ы

б ы т ь

э р г о и д н ы м и . Как

выясня­

ется,

эргоидные системы

состоят из достаточно большого числа доста­

точно

взаимодействующих частиц.

 

 

Эргоидная

гипотеза

совместно

с теоремой Лиувилля

приводит

к основным положениям статистической механики, которые иногда при­ нимают постулативно. Во-первых, это — постулат равной вероятности:

для

изолированной

системы все достижимые области

фазового

прост­

ранства имеют равные априорные

вероятности.

 

 

 

Во-вторых, теорема о средних значениях: среднее по времени

(доста­

точно длительному)

значение физически наблюдаемой величины F (q,p)

для

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

среднему-значению этой величины

по ансамблям.

 

Первое из этих

утверждений

вытекает из того, что

изображающая

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

* Д л я рассматриваемых консервативных систем функция Гамильтона выражает полную энергию системы.

184 —

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

некоторой динамической переменной

6 (р, а). Таким образом, исполь­

зуя верхнюю черту для обозначения

среднего по времени,

имеем

6 = 6 [q(t), p{t)}.

(VI.24)

Например, давление газа соответствует средней скорости переноса ко­

личества движения

на единицу

поверхности сосуда.

 

С р е д н е е п о

в р е м е н и

динамической переменной

задается

выражением

 

 

 

 

 

X

 

 

fl~M = 7

\ 4q(t), P(t)]dt,

(VI.25)

 

б

 

где т — некоторое время, «достаточно долгое», чтобы сделать возмож­ ным имеющим физическое значение измерение рассматриваемого термодинамического свойства. Среднее по ансамблю определяется дру­ гим соотношением:

| 6 (q,

р) | = \... \ 6 (q, p)p{q, р,

t) dqx...

dqfdpx

...dpf

и представляет

собой среднее значение

8 среди

всех систем ансамбля

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

Ш 7 ) = \В(Я, р)\.

(VI.26)

При этом следует понять, что равенство (VI.26) является основным для связи механического описания системы с термодинамическим и имеет смысл только для эргоидных систем.

§ 6. Подсчет микросостояний по Больцману

В системах, состоящих из большого числа одинаковых молекул, например одного моля какого-либо химически чистого газа, для опи­ сания механического состояния проще применять фазовое |л-про- странство 2/ измерений, если / — число степеней свободы молекулы (см. § 2 гл. VI). Как уже говорилось, точка в таком пространстве будет точно определять координаты (qx ... qf) и импульсы х ... р}) данной

185 —

единственной молекулы. Число измерений в у-пространстве, применя­ емом для описания состояния всей системы в целом, будет в N раз боль­ шим, т. е. равно 27V/, если N —число молекул в системе. Если между молекулами системы отсутствуют (значительные) силы взаимодей­ ствия *, то у-пространство системы можно представить совокупностью индивидуальных ^.-пространств. Местонахождение точки в у-прост- ранстве системы описывает положение изображающей точки каждой

молекулы

в ее собственном (^-пространстве

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

(м и к р о с о с т о я н и е)

всей системы

в

целом. Очевидно также,

что мгновенное микросостояние системы N молекул будет характери­

зоваться

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

изображающих

 

точек в ^-пространствах.

Как уже говорилось ранее, с классических

позиций данному равновес­

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

Существенно ограничить число микросостояний, чтобы оно не было бесконечно большим. Для этого достаточно условиться определять координаты и импульсы частиц не с абсолютной точностью, а считать,

что эти характеристики

лежат в пределах

от qt до qt + dqi и от p t

до Pi + dpi. Иными словами, чтобы сделать

множество

макросостоя­

ний счетным, фазовое

^-пространство разбивают на

элементарные

ячейки объемом **

 

 

 

Дсо = dqi... dqf dpi...

dpf

 

и считают состояние молекулы достаточно определенным, «ели соответ­ ствующая ей фигуративная точка находится в данной элементарной ячейке.

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

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

Пусть система состоит из ./V молекул, а ее макросостоянию соответ­ ствует Nx фигуративных точек в первой ячейке, N2 точек во второй и т. д. и в общем случае Ni точек в t'-той ячейке. Так как перестановки

* Но все-таки некоторое взаимодействие между молекулами должно иметь место, так как иначе будет исключен обмен энергией между молекулами и окажется невозможным установление равновесного распределения — система, лишенная любого взаимодействия между молекулами, не будет эргоидна.

** В квантовой статистике, используя принцип неопределенности Гейзенберга dqdp 5= h, где h — постоянная Планка (точнее — деленная на 2л), считают элемен­ тарную ячейку Дсо = h. В классической статистике объем ячейки остается нео­ пределенным.

-186 —

( V L 2 7 )

внутри ячейки не считаются — они не дают новых микросостояний, — число последних находится как число перестановок с повторениями:

G " W ' . W -

Проиллюстрируем метод следующим примером. Условимся рас­ сматривать распределение молекул только в обычном трехмерном

пространстве. Ограничимся

шестью

молекулами

и тремя ячейками:

1,

2,

3

 

 

4,

5,

6

 

 

I I

1, з

4,

6 5,2

 

На приведенной выше схеме показаны два

«макросостояния».

I все молекулы находятся в одной ячейке и I I молекулы распре­

делены равномерно по всем ячейкам. Число «микросостояний», соответ­

ствующих

первому

«макросостоянию»,

равно единице: все молекулы

в первой

ячейке,

а перестановки внутри ячейки не считаются. Для

второго «макросостояния» подсчитаем

G по формуле (VI.27):

 

 

G = 2I 21 2! =

= 9 0 ,

т. е. второе «макросостояние» согласно классической статистике может быть получено 90 способами.

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

§7. Термодинамическая вероятность

иэнтропия

Из постулата равной вероятности (см. § 5 гл. VI) следует, что все микросостояния системы, совместимые с заданными условиями (на­ пример, с условием постоянства энергии), математически равно­ вероятны. На первый взгляд представляется, что это не может отвечасть реальной действительности. В самом деле, сравним два микро­ состояния одного моля газа. Пусть в одном он занимает весь объем сосуда, скажем 10 л, и молекулы его движутся хаотически. Будем считать, что такое микросостояние соответствует равновесному мак­ росостоянию. В другом же все молекулы собрались в объеме 1 см3 и движутся параллельно. Представить себе самопроизвольный пере­

ход первого

равновесного

состояния во второе неравновесное дейст­

вительно трудно.

Однако

гипотеза

равных

вероятностей

приводит

к правильным следствиям

и,

по-видимому,

справедлива.

Все дело

в том, как

часто

встречаются

те или

иные

микросостояния.

Для уяснения сути дела вообразим себе колоду из 52 карт, а также машину, которая один раз в секунду всю колоду перетасовывает. В ходе этой процедуры может получиться любое расположение карт, в том числе и упорядоченное расположение по старшинству и мас­ тям. Вероятности любого расположения («микросостояния») одина-

— 187 —

ковы. Но упорядоченное расположение единственное, а неупорядо­ ченных, т. е. прочих,

52! «^ 103 7 .

При частоте смены комбинаций 1 сек'1 каждая из комбинаций, в том числе и упорядоченная, будет встречаться 1 раз в 1030 лет. Следова­ тельно, хотя все'расположения и равновероятны, появление упорядо­ ченного расположения чрезвычайно мало вероятно, так как оно един­ ственное из 1037 всевозможных расположений.

Равным образом можно было бы сравнить и повторяемость или время ожидания упоминавшихся выше микросостояний газа с тем лишь отличием, что для неравновесного состояния оно невообразимо велико, что его трудно выразить с помощью обычно принятых спосо­ бов написания чисел, а написать, к примеру, 10з в , возведенное в сте­ пень 1030, можно, но это мало что даст.

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

состояния (макросостояния).

 

W = const G,

(VI.28)

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

пропорциональна

числу микросостояний, соответствующих данному макросостоянию. Более того, можно положить коэффициент пропорциональности рав­ ным единице и рассматривать W как относительную вероятность, т. е. вероятность данного состояния по отношению к вероятности ка­ кого-то стандартного состояния.

Можно также рассматривать W как числитель математической вероятности, в знаменателе которой стояло бы число всех мыслимых для данной системы микросостояний, совместимых с законами со­

хранения. Итак,

 

W = G^l,

(VI.29)

причем для обычных молекулярных систем в равновесном состоянии W очень велико.

Статистический смысл второго закона термодинамики состоит согласно Л. Больцману в том, что изолированная система эволюцио­ нирует преимущественно в направлении большей термодинамической вероятности. Иными словами, второй закон термодинамики не имеет такого абсолютного значения, как первый закон, ибо самопроизволь­ ное образование неравновесных состояний, например самоуплотне­ ние газа, не является абсолютно невозможным. Просто термодинами­ ческая вероятность таких состояний очень мала: им соответствует малое число микросостояний, и поэтому они должны реализовываться очень редко.

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

183 — ,

сумма термодинамических вероятностей всех возможных неравновес­ ных состояний:

^равнов. сост ^

^неравное, сост.

всех

Итак, система эволюционирует в сторону увеличения термодинами­ ческой вероятности, а равновесное состояние соответствует наиболь­ шему, максимальному значению W. Но именно так изменяется и энт­ ропия изолированной системы, поэтому эти величины должны быть взаимно связаны, т. е.

S = f(W).

(VI.30)

Для установления вида этой функции рассмотрим две независи­ мые системы с вероятностями состояний Wx и W2. Если составить из этих систем новую сложную систему, то ее термодинамическая ве­ роятность, согласно теореме умножения вероятностей, выразится произведением

 

W = WjW,.

 

Энтропии отдельно взятых систем согласно

(VI.30)

S ^ / 0 ^ ) и St

= f(Wt),

а энтропия суммарной системы равна

сумме этих энтропии, так как

S представляет собой экстенсивное,

аддитивное свойство системы

или

S = S1

+ S2

 

 

 

 

 

f(W)

= f(Wl)

+ f(Wt).

(VI.31)

Равным образом согласно уравнению (VI.30)

f(W1Wa)

= f(Wd

+ f(Wt).

Это выражение дифференцируем сначала по W{.

W2f

(W1Wt)

= f (Wi)

и вторично дифференцируем по W2:

 

 

 

/' (W^t)

+ WjWJ"

(WXW2)

= 0

или

 

 

 

 

/' (W) + Wf (W) = 0.

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

i + w у щ - - и или f l { W ) aw ~ — - w - .

Интегрируя, находим

\nf'(W) = — lnW + lnk = ln ~ ,

причем In к — константа интегрирования.

169 —

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