Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
80
Добавлен:
20.05.2015
Размер:
8.64 Mб
Скачать

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

(1.57) адитивна. У випадку

чистого стану n0 ентропія

системи дорівнює нулю. Дійсно, у цьому стані

 

1, n n0 ,

 

n nn

 

 

0

0, n n0.

 

 

 

 

Тоді

 

 

 

nn ln nn

ln n n 0.

0

0

0

0

n

Переконаємося в тому, що коли ентропія системи максимально можлива при заданих зовнішніх умовах, то функція розподілу від індексу стану n не залежить. Це відповідає повній невизначеності інформації про стан системи. Відповідний ансамбль є найбільш хаотичним.

З'ясуємо, яка функція розподілу приводить до максимуму ентропії при одній додатковій умові:

n 1.

n

Це – відома задача на умовний екстремум. Вона розв’язується методом множників Лагранжа. Необхідною умовою екстремуму ентропії є рівність нулю її варіації за функцією розподілу:

ln n 1 n 0.

(1.59)

n

 

Звідси не можна зробити висновок, що вираз в дужках дорівнює нулю, тому що величини n не незалежні. Вони

пов'язані однією умовою:

 

n 0.

(1.60)

n

 

За цією умовою одна з величин, наприклад,

n , може

 

1

51

бути виражена через інші незалежні величини:

n1 n2 n3 ...

Домножимо рівність (1.60) на множник Лагранжа , що не залежить від n, і складемо його з (1.59):

ln n 1 n 0.

n

Виділимо в цій сумі доданок з n n1 :

ln n1 1 n1 ln n 1 n 0.

n n1

Підберемо так, щоб перший доданок став рівним нулю:

ln n1 1 0.

Оскільки у

величини n незалежні, маємо:

n n1

 

 

ln n 1 0, n n2, n3,...

У результаті ми одержали замкнуту систему рівнянь для

n і :

ln n 1 0, n n1, n2,...,

n 1.

n

З цієї системи випливає

n e 1 ,

тобто функція розподілу дійсно не залежить від індексу стану n.

З формули (1.59) випливає, що варіаційна похідна ентропії за функцією розподілу дорівнює:

52

1 ln n.n

Отже,

2 1 0.n2 n

Ентропія стану, у якому n не залежить від n, макси-

мальна.

Формула (1.57) справедлива і для класичної системи. Тільки середнє визначається тепер формулою (1.27):

 

 

d

q, p ln q, p .

(1.61)

 

 

 

 

 

h3N N !

 

Видно, що ентропія (1.58) і (1.61) безрозмірна.

Застосуємо формулу (1.58) до мікроканонічного ансамблю. Підставляючи функцію розподілу (1.56) у вираз (1.58), одержуємо

lnW ,

(1.62)

де W – статистична вага макростану системи. Формула (1.62) має ім'я Больцмана.

Як приклад використання формули (1.62) обчислимо ентропію дворівневої системи. Розглянемо N незалежних частинок, кожна з яких має внутрішню ступінь свободи. Останній відповідають два рівні енергії: 0 , 0 . Одержи-

мо внесок внутрішніх ступенів свободи частинок в ентропію. Для цього необхідно знайти число мікростанів W , що реалізують макростан системи. Макростан характеризується числом частинок N N N і енергією системи

E 0 N N 0M ,

де N – числа частинок на верхньому і нижньому рівнях.

53

Мікростан системи – певний розподіл частинок за рівнями. Оскільки переставляння частинок на одному рівні не приводять до нових мікростанів, знаходимо:

 

 

W C N

 

N !

.

 

 

 

 

 

 

 

N

 

N !N !

 

 

 

 

З огляду на

 

 

 

 

 

 

 

N

1

N M ,

 

N

1

N M ,

2

 

 

 

 

 

2

 

 

одержуємо для ентропії дворівневої системи вираз:

ln

 

 

N !

 

 

.

 

 

 

1

1

 

 

 

 

N M !

 

N M

!

 

 

 

2

2

 

 

 

Ентропія залежить від числа частинок у системі і її енергії.

1.10. Ентропія ідеального газу

Визначимо внесок поступальних ступенів свободи частинок в ентропію класичної системи. Частинки системи вважаємо незалежними. Тоді гамільтонова функція системи залежить тільки від імпульсів частинок:

3N

2

 

 

H p

pi

.

(1.63)

 

i 1

2m

 

Тут m – маса частинки, індекс і нумерує компоненти імпульсів частинок. Поверхнею постійної енергії Е в імпульсному просторі є сфера

pi2 2mE

i

радіусом 2mE .

Статвага макростану системи (1.14) пропорційна об’єму кульового шару у фазовому просторі між близькими

54

гіперповерхнями постійної енергії. Оскільки підінтегральна функція в (1.14) не залежить від координат частинок, інтегрування за координатами кожної частинки дає об’єм системи V , а інтеграл за координатами усіх частинок

дорівнює V N . Залишається

обчислити об’єм кульового

шару між гіперсферами

H p E і

H p E E

в імпульсному просторі.

Переконаємося в тому, що при великому числі ступенів свободи системи майже весь об’єм багатомірної кулі зосереджений у тонкому приповерхневому шарі.

Відомо, що об’єм k -мірної кулі дорівнює

V C Rk ,

(1.64)

k k

 

 

де R – її радіус, а

 

 

 

 

 

Ck

k

2

.

 

 

 

 

k

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

Тут – гамма-функція Ейлера. Ми сподіваємося, що позначення гамма-функції і числа станів однією буквою не призведе до непорозумінь. Об’єм кульового шару між сферами з радіусами R і R R дорівнює:

C Rk C

 

R R k V

 

 

1

R

k

k

1

.

k

k

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

Цей вираз можна записати так:

 

 

 

 

 

 

 

 

R

k

 

 

 

k

 

 

 

 

 

 

R

 

 

Vk

1 1

 

 

 

.

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

55

При k 1 він приблизно дорівнює

 

 

R

 

V

1 e k

R

V .

k

 

k

 

 

 

 

Це означає, що статвагу у формулі Больцмана (1.62) для ентропії можна замінити числом станів E під поверхнею постійної енергії:

ln .

(1.65)

Використовуючи формулу (1.64) для об’єму кулі радіусом

R 2mE 12 у 3N -мірному імпульсному просторі,

одер-

жуємо для числа станів вираз

 

 

 

E,V

3N 2V N 2m 3N

2 E3N 2

.

(1.66)

 

3N

3N

 

 

 

 

 

 

h

 

N !

 

1

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

1

 

 

 

 

 

 

 

При

N 1 з урахуванням

z 1 z z ,

 

 

 

 

 

 

 

 

2

 

 

 

звідси випливає отриманий раніше вираз для одночастинкової густини станів (1.20).

Оскільки N 1, можна у (1.66) використати формулу

Стірлінга:

 

 

 

 

N N

 

N !

 

.

(1.67)

 

 

e

 

Тут e 2,7 – основа натуральних логарифмів. Тоді ентропія ідеального газу дорівнює:

 

 

4 em

3

V

E

3

 

 

 

2

2

 

 

E,V N ln e

 

 

 

 

 

 

 

 

.

(1.68)

3h2

 

 

 

 

 

 

 

 

N

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

56

Цей вираз називається формулою Сакура – Тетроде.

З формули (1.68) видно, що ентропія адитивна. З ростом

числа

частинок при незмінних питомих величинах

V

 

 

N

і E

вона зростає пропорційно N . Відзначимо, що ади-

N

 

 

тивність ентропії досягнута завдяки врахуванню N !

у зна-

меннику (1.10). Якби тотожність однакових частинок не була врахована, адитивність ентропії була б загублена. Ця ситуація відома під назвою парадокса Гіббса.

Розв’язуючи рівняння (1.68) відносно енергії, одер-

жуємо:

 

 

 

 

 

 

 

 

 

 

E ,V N

3h2

 

 

N 23

2

 

 

 

 

 

 

 

exp

 

 

1 .

(1.69)

 

23

 

 

 

4 me

V

3 N

 

 

Енергія ідеального газу виражена через параметри V і , що характеризують його макростан. Із (1.69) легко знаходимо частинні похідні

E

 

2 E

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

3 V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.70)

 

E

 

2 E

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

3 N

 

 

 

які будуть використані надалі.

1.11. Закон зростання ентропії

Ентропія системи, що перебуває в стані статистичної рівноваги, визначається співвідношеннями (1.58) і (1.61). Функція розподілу, що входить у ці формули, явно не залежить від часу. Природно припустити, що ентропія нерівноважної системи дається тими ж формулами (1.58) і (1.61), у які необхідно підставити нерівноважну функцію

57

розподілу. Тоді ентропія нерівноважної системи буде залежати від часу. У класичній статистиці вона дорівнює:

 

t d q, p,t ln q, p,t ,

(1.71)

де d

dqdp

 

– число станів в об’ємі dqdp .

 

h3N N !

 

 

 

 

Відповідно до нульового принципу термодинаміки нерівноважна система, полишена на саму себе, прагне перейти в стан рівноваги. У цьому стані функція розподілу і макроскопічні величини від часу не залежать. Перехід до рівноваги супроводжується зростанням ентропії системи. Закон зростання ентропії замкнутої системи

d t

0

(1.72)

dt

 

 

є одним із формулювань другого принципу термодинаміки. Знак рівності в (1.72) стосується рівноважної системи, ентропія якої не змінюється, а знак – нерівноважної системи. Процеси, при яких ентропія замкнутої системи не змінюється, називаються оборотними, а процеси, що супроводжуються зростанням ентропії, називаються необоротними.

Спробуємо одержати нерівність (1.72), використовуючи формулу (1.71). Для цього порівняємо ентропію нерівноважної системи в момент t 0

0 d 0 q0 , p0 ,0 ln q0 , p0 ,0

з ентропією (1.71) у пізніший момент t . Нехай точкиq0 , p0 і q, p лежать на фазовій траєкторії системи.

Згідно з теоремою Ліувілля,

q0 , p0 , 0 q, p,t , d 0 d .

58

Отже, ентропія нерівноважної системи, визначена співвідношенням (1.71), не залежить від часу: t 0 . Цей

результат можна одержати шляхом диференціювання ентропії (1.71) за часом з використанням рівняння Ліувілля (1.46). Таким чином, рівняння Ліувілля не суперечить детермінізму класичної механіки.

Щоб одержати закон зростання ентропії, Пауль і Тетяна Еренфести запропонували замінити точну функцію розподілу у (1.71) «огрубленою». Вони усереднили функцію

розподілу

за малими комірками об’єму

i

у фазовому

просторі системи:

 

 

 

 

 

i t

1

d q, p,t .

 

(1.73)

 

 

 

 

 

i

 

 

 

 

 

i

 

 

Об’єм i

вибирається не дуже малим, щоб він мав багато

точок ансамблю. З іншого боку, об’єм i не повинен бути

занадто великим, щоб цілком не згладити функцію розподілу. Функція (1.73) стала в межах комірки i . Вона нормо-

вана так само, як і функція . Дійсно,

d i i i i

1

d d d . (1.74)

 

i

 

 

 

i

i

i

 

 

 

 

 

 

 

 

i

i

Підставимо (1.73) у формулу Гіббса (1.61):

t

 

1

 

 

d q, p,t ln q, p, t

 

 

 

 

h

3N

N !

 

 

 

 

 

 

(1.75)

 

1

 

i i ln i ,

 

 

 

 

 

 

 

 

 

 

h3N N ! i

 

 

 

 

 

 

де q, p – координати елемента об’єму i . Переконаємося в тому, що ентропія замкнутої системи, визначена співвід-

59

ношенням (1.75), зростає з часом. Припустимо, що в початковий момент часу огрублена функція розподілу збігається з точною:

 

q0 , p0 ,0 q0 , p0 ,0 .

 

(1.76)

Порівняємо

ентропію

(1.75)

у момент

t з ентропією

в початковий момент часу:

 

 

 

0

1

 

d

0

q0 , p0 ,0 ln q0 , p0 ,0 ,

(1.77)

 

 

h3N N !

 

 

 

 

 

 

 

де точки

q0 , p0

 

і

q, p

лежать на

одній

фазовій

траєкторії. З огляду на d 0 d і початкову умову (1.76), одержуємо

0

 

1

 

 

d

q0 , p0 ,0 ln q0 , p0 ,0

.

 

 

 

 

 

h3N N !

 

 

 

 

 

Згідно з теоремою Ліувілля,

 

 

 

 

q0 , p0 , 0 q, p,t .

 

Тоді

 

 

 

 

 

 

 

 

0

 

1

 

 

 

d q, p,t ln q, p, t .

(1.78)

 

 

 

 

 

h

3N

N !

 

 

 

 

 

 

У співвідношенні (1.75) можна зняти знак ~ над функцією розподілу, яка стоїть перед логарифмом. Дійсно,

t

 

 

 

1

 

 

i i ln i

 

1

 

 

 

ln i d

 

 

 

 

 

 

 

 

 

 

 

h

3N

N !

h

3N

N !

 

 

 

 

 

 

 

i

 

 

 

i

i

(1.79)

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

d ln i

 

 

 

d ln .

 

 

 

 

 

 

 

 

 

 

 

h

3N

N !

h

3N

N !

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

Вираз ln i

сталий усередині комірки i , і тому він вне-

60

Соседние файлы в папке Статы