Статы Экзамен / Статы / statukr
.pdfрозподілу макросистеми мультиплікативна, тому ентропія
(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