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

РОЗДІЛ 3. РОЗПОДІЛ ГІББСА

3.1.Канонічний розподіл

Уцьому підрозділі ми одержимо функцію розподілу

підсистеми в термостаті. Припустимо, що підсистема і термостат утворюють замкнуту систему з енергією E0 .

Між підсистемою і термостатом відбувається обмін енергією, а обмін частинками заборонений. Уся система перебуває в стані рівноваги. Макроскопічний стан підсистеми характеризується її об’ємом V і температурою T. Сукупність таких підсистем з однаковими V і T утворює канонічний ансамбль Гіббса.

Будемо вважати, що частинки системи підкоряються законам класичної механіки. Нехай q і p – набір узагаль-

нених координат і імпульсів підсистеми, а q і p – термо-

стата. Тоді

H q, p H q , p E0 ,

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

Функцію

розподілу замкнутої

 

системи позначимо

 

 

. Вона нормована умовою (1.22):

q, p, q , p

 

 

 

 

 

1,

 

 

d d q, p, q , p

де

 

 

 

 

 

 

 

 

d

dqdp

 

 

 

 

 

N !h3N

 

 

 

 

 

 

 

– число станів підсистеми в елементі об’єму її фазового простору, а

111

d dq dp N !h3N

– число станів термостата. Тут N і N – числа частинок підсистеми і термостата відповідно. У цих виразах враховується тотожність частинок підсистеми і термостата окремо, однак частинки підсистеми можуть відрізнятися від частинок термостата.

Як і в підрозділі 1.8, зручно вважати, що енергія всієї

системи

не фіксована точно, а перебуває в проміжку

E0 , E0 E0 . Тоді функція розподілу всієї системи, згідно

з (1.54), дорівнює:

 

 

 

 

1

, E0 H H E0 E0 ,

 

 

 

 

 

 

E0

q, p, q , p W0

 

 

 

 

 

 

0, H H E0 , H H E0 E0 ,

(3.1)

 

 

 

 

де W0 – статистична вага всієї системи. Вона залежить від

E0 , E0

і об’єму системи. Функція розподілу підсистеми

q, p

може бути отримана з (3.1) шляхом інтегрування

за фазовими координатами термостата:

q, p d q, p, q , p .

Підставляючи сюди (3.1), знаходимо:

 

 

1

 

 

 

 

 

 

 

q, p

 

 

 

 

 

dq dp

 

.

 

E0

 

 

 

N !h

3N

W0

 

 

 

 

 

 

 

 

 

 

E0 H q, p H q , p

 

 

 

 

 

 

 

 

E0 H q, p E0

 

 

 

 

 

 

Інтеграл, який сюди входить, дорівнює статистичній вазі термостата W при енергії E0 H q, p . Отже,

112

 

 

W E0 H q, p

 

q, p

 

 

 

 

 

 

 

 

 

 

 

.

(3.2)

 

 

W0

E0

 

 

 

 

 

 

 

 

 

 

Згідно з (1.62), функції W0 і W можуть бути виражені че-

рез ентропію всієї системи S0

та ентропію термостата S :

W

E

exp

 

1

S

 

E

 

,

 

 

 

 

0

 

 

0

0

 

 

 

 

0

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

W H

exp

 

S H

 

,

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

де k – постійна Больцмана. У результаті функція розподілу (3.2) дорівнює:

q, p exp

1

S E

H q, p

1

S

 

E

 

. (3.3)

 

 

 

0

 

0

 

k

 

0

 

 

 

k

 

 

 

 

 

 

 

Будемо вважати, що енергія підсистеми мала в порівнянні з енергією всієї системи. Тоді ентропія S може бути розкладена в ряд за ступенями H . Обмежуючись лінійним членом розкладу, одержуємо:

S E H q, p

S E

H q, p

S E0

 

.

 

 

0

 

0

 

E0

 

 

 

 

 

 

 

 

(3.4)

Похідна обчислюється при постійному об’ємі і дорівнює:

 

E0

 

 

1

 

 

S

 

,

(3.5)

E

 

T

 

 

 

 

 

0

 

 

 

 

 

де T – температура термостата, яка збігається з температурою підсистеми. Підставляючи (3.4) і (3.5) у (3.3), одержуємо:

113

1

 

 

 

S E

 

 

H q, p

 

 

 

 

 

 

 

 

S E

 

 

q, p ek

 

0

0 0

e

 

kT .

Цей розподіл називається канонічним розподілом Гіббса. Запишемо його у вигляді:

q, p

1

 

 

H

 

q, p

 

 

 

 

exp

 

 

.

(3.6)

Z

 

kT

 

 

 

 

 

 

 

 

 

Постійна Z може бути знайдена з умови нормування

(1.22):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dqdp

 

H q, p

 

 

 

Z V ,T

 

 

 

 

 

 

 

 

 

e

 

kT

 

 

.

(3.7)

h3N N !

 

 

 

 

 

 

 

 

 

 

 

 

Вираз (3.7) називається статистичним інтегралом.

 

Ясно, що у випадку квантової підсистеми в термостаті функція розподілу n , введена у підрозділі 1.3, має вигляд:

 

 

 

1

 

 

En

 

 

 

 

 

 

 

 

 

n

 

 

e kT ,

(3.8)

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

де En – рівні енергії підсистеми, а

 

Z e

 

En

 

 

kT

(3.9)

 

n

– статистична сума (статсума). Вона отримана з умови нормування (1.23). Вираз (3.8) являє собою власне значення статистичного оператора

 

 

 

ˆ

 

 

1

 

H

 

ˆ

kT

 

 

e

,

Z

 

 

 

 

(3.10)

де ˆ – оператор Гамільтона підсистеми.

H

Із формул (3.6) і (3.8) видно, що температура рівноважної системи повинна бути позитивною. У протилежному випадку імовірність стану зростала б із збільшенням енер-

114

гії, що безглуздо. Видно також, що розподіл (3.6) збігається з формулою (1.52), отриманою з рівняння Ліувілля. Функція розподілу за енергією (1.55) відрізняється від функції розподілу (3.6) у фазовому просторі додатковим множником, рівним густині станів підсистеми.

3.2. Канонічний розподіл Гіббса і термодинаміка

Маючи функції розподілу (3.6) і (3.8) для системи в термостаті, ми можемо термодинамічні величини такої системи обчислювати за формулами (1.27) і (1.28). Однак кращим виявляється інший метод.

Замість нормувальної константи Z , яка входить у формули (3.6) і (3.8), введемо іншу константу F , зміст якої необхідно з'ясувати:

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

1

e kT .

 

 

 

 

(3.11)

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тоді розподіли (3.6) і (3.8) набувають вигляду:

 

 

 

 

 

 

F H q, p

 

 

 

q, p

e

 

 

 

 

 

 

 

 

 

kT

,

(3.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F En

 

 

 

 

 

 

n

e kT

.

 

 

(3.13)

 

 

 

 

 

 

 

 

 

 

 

 

Ентропія системи (1.57) дорівнює:

 

 

 

 

 

 

 

 

 

 

 

 

 

S k

 

k

F En

 

 

F E

 

ln

,

 

 

 

 

 

 

 

kT

 

 

 

T

 

де E En – внутрішня енергія системи. Ми одержали

термодинамічне співвідношення (2.36), яке пов'язує вільну енергію з внутрішньою. Отже, величина F у (3.11) являє собою вільну енергію системи. Вона дорівнює:

F V ,T kT ln Z V ,T .

(3.14)

115

Це співвідношення справедливе як у квантовій статистиці, так і в класичній. Воно дозволяє одержати вільну енергію, якщо відомі статистичний інтеграл (3.7) або статистична сума (3.9). Інші величини обчислюються так, як описано в розділі 2. Зокрема, внутрішня енергія системи пов'язана із статсумою співвідношенням:

E kT 2

 

ln Z ,

T

 

 

де похідна береться при постійному об’ємі і постійному числі частинок.

3.3. Великий канонічний розподіл

Розглянемо рівноважну відкриту підсистему в термостаті. Вона обмінюється з термостатом не тільки енергією, але і частинками. Її макростан характеризується об’ємом V , температурою T і хімічним потенціалом . Сукупність

таких підсистем утворює великий канонічний ансамбль Гіббса. Оскільки число частинок у підсистемі фігурує як додаткова змінна, визначення функції розподілу, введеної в підрозділі 1.3, потребує уточнення.

Нехай Nn – імовірність того, що квантова система

містить N частинок і перебуває у квантовому стані n. Умова нормування для цієї функції розподілу має вигляд:

 

 

Nn 1.

(3.15)

N 0 n

Внутрішня енергія системи і середнє число частинок в ній рівні:

 

 

Nn ENn ,

(3.16)

E

 

 

Nn

 

116

 

 

Nn N ,

(3.17)

N

 

 

Nn

 

де ENn – рівень енергії системи N частинок, які

перебувають у стані n.

У класичній статистиці

N q, p dqdp h3N N !

являє собою імовірність того, що система містить N частинок, а її фазова точка міститься в елементі об’єму dqdp її фазового простору. Ця функція нормована умовою:

 

 

dqdp

 

 

 

 

N q, p 1.

(3.18)

 

 

N 0

h3N N !

 

Середні в класичній статистиці обчислюються за допомогою формули

 

 

 

 

dqdp

 

 

 

 

 

 

N q, p GN q, p ,

(3.19)

G

 

 

 

 

N 0

h3N N !

 

де GN q, p – динамічна змінна системи N частинок. Функція розподілу Nn рівноважної відкритої підсисте-

ми в термостаті може бути отримана методом, використаним у підрозділі 3.1. Ми використаємо інший метод. Він придатний і у випадку канонічного ансамблю Гіббса.

Ентропія підсистеми дорівнює:

S k Nn ln Nn .

(3.20)

Nn

 

З'ясуємо, яка функція розподілу приводить до максимуму ентропії при додаткових умовах (3.15)–(3.17). Умови (3.16) і (3.17) фіксують середню енергію системи і середнє число частинок в ній. Ця задача розв’язується методом Лагранжа.

117

Варіюючи ентропію (3.20) і умови (3.15)–(3.17) за функцією розподілу Nn , прирівнюючи ці варіації до нуля,

одержуємо:

 

ln Nn 1 Nn 0 ,

(3.21)

Nn

 

Nn 0,

(3.22)

Nn

 

ENn Nn 0,

(3.23)

Nn

 

N Nn 0.

(3.24)

Nn

 

Домножимо (3.22)–(3.24) на множники Лагранжа , ,і складемо з (3.21):

ln Nn 1 ENn N Nn 0.

Nn

 

Звідси випливає

 

ln Nn 1 ENn N 0.

(3.25)

Система рівнянь (3.25) разом із (3.15)–(3.17) дозволяє визначити Nn і множники Лагранжа , , . Із (3.25)

знаходимо:

Nn e 1 e ENn N .

Такий розподіл називається великим канонічним розподілом Гіббса. Запишемо його у вигляді:

Nn

1

e

ENn N

.

(3.26)

 

 

 

 

 

 

 

Цій функції відповідає статистичний оператор

118

 

1

ˆ ˆ

 

 

ˆ

e H N

,

(3.27)

 

 

 

 

 

ˆ

де N – оператор числа частинок у підсистемі. Константаможе бути знайдена з умови нормування (3.15):

e

ENn N

.

(3.28)

 

Nn

Ця величина називається великою статистичною сумою. Якщо обмін частинками між підсистемою і термостатом

відсутній, розподіл (3.26) переходить у канонічний розподіл (3.8). Це означає, що множник Лагранжа пов'язаний

зтемпературою співвідношенням 1kT .

Утому випадку, коли підсистема підкоряється законам класичної механіки, розподіл (3.26) переходить у

N q, p

1

 

 

HN q, p N

 

 

 

e

 

,

(3.29)

 

 

 

 

 

 

 

 

 

де HN q, p – гамільтонова функція, а

 

 

 

 

dqdp

 

 

H N q, p N

 

 

 

 

 

 

e

 

 

(3.30)

 

 

 

 

 

h3N N !

 

 

 

 

N

 

 

 

 

 

 

– великий статистичний інтеграл. Вирази (3.28) і (3.30) можуть бути записані у вигляді:

e N ZN ,

(3.31)

N

 

де ZN – статистичний інтеграл (3.7) або статистична сума (3.9) системи N частинок.

3.4. Великий канонічний розподіл і термодинаміка

Щоб з'ясувати зміст константи у (3.26), покладемо

119

1

e ,

(3.32)

 

 

 

де – нова константа. Тоді функції (3.26) і (3.29) набувають вигляду:

 

Nn e

ENn N

,

 

 

(3.33)

 

 

 

 

 

 

 

N q, p e

HN q, p

N

(3.34)

 

 

 

.

Ентропія відкритої системи дорівнює:

 

 

 

 

S k

 

k

 

 

 

.

 

ln

 

E

N

 

Ця рівність збігається з (2.91), якщо – великий потенціал, а – хімічний потенціал підсистеми.

Із формули (3.32) випливає:

 

V ,T , kT ln V ,T , .

(3.35)

Ця формула дозволяє обчислити великий потенціал відкритої системи, якщо відомі велика статистична сума (3.28) або великий статистичний інтеграл (3.30).

3.5. Ізобарично-ізотермічний ансамбль

Наприкінці підрозділу 1.1 був згаданий ще один тип контакту системи з термостатом – система в циліндрі під поршнем. Її температура T , тиск P і число частинок фіксовані, а об’єм флуктує. Відповідний ансамбль Гіббса називається ізобарично-ізотермічним.

Функція розподілу nV квантової системи в циліндрі визначається так, що nV dV являє собою імовірність того, що об’єм системи міститься в інтервалі V ,V dV і вона перебуває в квантовому стані n . Ця функція нормована умовою:

120

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