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

У першому доданку в правій частині цього виразу усереднюється адитивна величина, а в другому – бінарна. Такі середні розглядалися в підрозділі 6.2. Використовуючи формули (6.9) і (6.12), придатні для будь-яких адитивних і бінарних величин, знаходимо

 

n r1 n r2

 

n d

3

 

 

r2

 

 

 

 

 

 

 

 

 

r r1

r

 

r

 

 

 

n

2

d

3

 

d

3

 

 

 

 

r2 r

 

 

 

r

 

r

 

 

 

 

 

 

 

 

 

 

 

r

 

r r1 r

 

 

 

 

 

(7.50)

 

n r1 r2 n 2

 

r1 r2

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

де – введена в підрозділі 6.1 радіальна функція розподілу. Підставляючи (7.50) у (7.49), одержуємо

C

 

r1 r2

 

n r1 r2 n v

 

r1 r2

 

,

(7.51)

 

 

 

 

де

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.52)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

r

 

n

 

r

 

1 .

 

 

 

 

 

 

 

 

 

 

 

При r

r1 r2

флуктуації

 

густини в точках

r1 і r2

стають статистично незалежними, функція v прямує до нуля, а – до одиниці.

Дисперсія числа частинок N0 2 у виділеному в си-

стемі об’ємі V0 дорівнює

N0 2 d 3r1 d 3r2 n r1 n r2 ,

де інтегрування виконується за об’ємом V0 . Підставляючи сюди (7.51), знаходимо

N0

2

 

 

d

3

 

(7.53)

 

n V0 1

 

rv r .

260

Вхідний у цю формулу інтеграл виразимо через ізотермічну стисливість. З цією метою використовуємо отриманий у підрозділі 7.4 вираз для дисперсії числа частинок у виділеному об’ємі V0 . Тоді (7.53) можна записати у вигляді

d 3rv r 1 kT

N

 

V

 

 

 

 

0

 

0

.

(7.54)

V

2

 

 

 

 

P

 

 

0

 

 

 

T

 

Якщо кореляція між положеннями

частинок

відсутня

v 0 , з формули (7.53) випливає отриманий у підрозділі 7.1 вираз для дисперсії числа частинок в ідеальному газі:

N0 2 nV0.

Розкладемо флуктуючу густину (7.48) у ряд Фур'є:

 

 

 

 

 

 

 

 

 

n r

1

nqeiqr ,

 

(7.55)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

q

 

 

 

 

 

 

 

 

 

 

де V – обмежений об’єм,

nq

компонента Фур'є. Вона

дорівнює

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nq d 3r n r e iqr .

 

(7.56)

Оскільки функція n r дійсна, з (7.56) одержуємо

 

 

 

 

 

 

 

 

 

 

n

n

 

 

.

 

 

 

 

 

 

(7.57)

 

 

 

 

 

 

 

 

 

 

q

 

 

q

 

 

 

 

 

 

 

 

Виразимо

 

 

n

 

 

2

через функцію v :

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

2

 

 

n n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

q

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

e iq r1 r2

 

 

1

 

 

 

 

 

2

 

 

 

(7.58)

 

 

 

 

 

 

 

 

 

 

 

d3r

 

d3r

n

r

 

n

 

r

 

 

 

N 1

v q ,

261

де N – число частинок в об’ємі V ,

v q d 3rv r e iqr

– компонента Фур'є функції v r . Тут ми врахували (7.51) і перейшли від інтегрування за r1 і r2 до інтегрування за

r1 r2 і r2 .

Нехай повний об’єм системи V і її температура T фіксовані. Тоді мінімальна робота, необхідна для відхилення системи від стану рівноваги, дорівнює

Rmin FП ,

де FП – зміна повної вільної енергії системи. Згідно

з (7.13) імовірність такого відхилення в результаті флуктуацій дорівнює

 

 

F

 

(7.59)

w exp

П .

 

 

kT

 

 

Зміну вільної енергії тіла, обумовлену флуктуаціями густини, представимо у вигляді

2

 

 

 

d 3r2

 

 

n r1 n r2 ,

 

FП

1

 

d 3r1

 

r1

r2

(7.60)

 

 

 

де – друга варіаційна похідна від FП за n. Переходя-

чи в (7.60) до компонентів Фур'є і враховуючи (7.57), знаходимо

FП

1

q

 

nq

 

2 ,

(7.61)

 

 

2V

 

q

 

 

 

де q – фур'є-компонента функції r .

Підставляючи (7.61) у (7.59), одержуємо імовірність флуктуації у вигляді

262

 

1

q

 

2

 

w exp

nq

.

2kTV

 

 

q

 

 

 

 

 

 

Звідси видно, що флуктуації різних компонентів Фур'є nq статистично незалежні і є гауссівськими. Згідно з (7.57),

 

2 входить у суму (7.61) двічі (з q

кожен квадрат

n

 

q

 

і q ), тому розподіл імовірностей флуктуацій nq дається виразом

exp

 

1

q

 

n

 

2 .

(7.62)

 

 

 

 

 

 

 

 

q

 

 

 

 

 

kTV

 

 

 

 

 

 

Якщо врахувати, що

nq 2 Re nq 2 Im nq 2

– сума квадратів двох незалежних величин, і порівняти формулу (7.62) з (7.5), одержимо дисперсію флуктуючої величини nq :

n

 

2

 

kTV

.

(7.63)

 

 

 

 

q

 

 

q

 

 

 

 

 

 

 

Порівнюючи цей вираз з (7.58), знаходимо

v q

kT

1.

(7.64)

n q

Ця формула дозволяє одержати радіальну функцію розподілу (7.52), якщо відома функція q .

Зміну вільної енергії тіла представимо у вигляді інтеграла за об’ємом:

F

 

 

d 3r f r

f

,

(7.65)

П

 

 

 

 

 

263

де f r – флуктуюча густина вільної енергії, f – її середнє значення. Воно постійне вздовж тіла. При малих флуктуаціях розкладемо f r f в ряд за ступенями відхилення (7.48) густини від середнього значення. Лінійний член розкладу пропорційний n r . При інтегруванні за об’ємом він дає нуль, тому що

d 3rn r n V .

Квадратичний член дорівнює

1

 

2 f

 

n r

2

 

 

 

 

 

2

 

 

,

(7.66)

 

 

 

2

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

де друга похідна від густини вільної енергії обчислюється у стані рівноваги. Оскільки f стосується одиниці об’єму

речовини, то

fn

– хімічний потенціал. Отже,

 

2 f

 

 

 

 

n2

 

 

.

 

 

 

n

 

 

T

 

T

Але при T const маємо

d 1n dP.

Тому

 

2 f

 

 

 

 

1

P

 

 

n2

 

 

 

 

 

 

.

(7.67)

 

 

 

 

n

 

n

n

 

 

 

T

 

T

 

 

 

T

 

264

Підставляючи (7.66) з урахуванням (7.67) у (7.65) і переходячи до компонентів Фур'є (7.56), одержуємо

FП

1

 

P

1

 

 

nq

 

2

.

(7.68)

 

 

 

 

 

 

 

 

 

 

 

 

2n

n T V

q

 

 

Порівняння цієї формули з (7.61) дає

q

1

 

P

3

 

P

,

(7.69)

 

 

 

 

 

 

 

n

 

n T

 

 

T

 

 

де 1n – питомий об’єм. Компонента Фур'є q у розглянутому наближенні виявилася незалежною від q. Це означає, що у формулі (7.61) виконаний розклад q в ряд за ступенями q і в цьому розкладі залишений головний член 0 (7.69). Іншими словами, вирази (7.65) і (7.66) дозволяють розглянути лише довгохвильові флуктуації густини. Члени порядку q2 при розкладі q бу-

дуть ураховані нижче. Отримані в цьому підрозділі співвідношення будуть використані в теорії фазових переходів.

7.9. Узагальнена сприйнятливість

Щоб пов'язати флуктуації з дисипацією енергії в системі, розглянемо її реакцію на змінне зовнішнє поле. Гамільтонову функцію системи представимо у вигляді

H H0 V ,

(7.70)

де H0 – гамільтонова функція за

відсутності поля, а V

враховує взаємодію системи із заданим зовнішнім полем. Нехай V залежить від часу за законом

V t xf t ,

(7.71)

265

де x – деяка величина, яка характеризує систему, f t – задана функція часу. У квантовій теорії H (7.70) необхідно

замінити оператором Гамільтона ˆ системи, а – ерміто-

H x

вим оператором xˆ. Наприклад, розглянемо взаємодію системи заряджених частинок із змінним електричним полемt . Гамільтоніан взаємодії дорівнює

ˆ

P t ,

де ˆ – оператор електричного дипольного моменту систе-

P

ми. Гамільтоніан взаємодії магнетика із змінним магнітним полем H t дорівнює

ˆ

MH t ,

де ˆ – оператор магнітного моменту.

M

Будемо вважати, що в момент t система знаходилася у стані рівноваги, а середнє значення x дорівнювало нулю. У цей момент включається зовнішнє поле. Воно виводить систему зі стану рівноваги. З'являється середнє зна-

чення x t , індуковане зовнішнім полем. Воно є відгуком

системи на зовнішнє поле. Якщо поле слабке, зв'язок між відгуком x і причиною f , яка викликає цей відгук, буде

лінійним. Ми розглянемо випадок, коли цей зв'язок однорідний. Якщо f 0, відгук x відсутній. Крім того, цей

зв'язок повинен задовольняти принцип причинності. Причина повинна передувати наслідку. Це означає, що зв'язок між x і f повинен бути таким:

 

t

 

 

.

 

 

 

 

 

(7.72)

x t dt t t

f t

266

У цьому співвідношенні вважається, що зв'язок x з f нелокальний. Значення x у момент t визначається значеннями f в усі попередні моменти часу. У наступному підрозділі буде показано, що функція , яка входить у (7.72),

зображується рівноважним середнім типу (7.27). Вона залежить від різниці часів t t . Принцип причинності вимагає,

щоб t t 0. Переходячи в (7.72) до нової змінної інте-

 

 

грування t t , можемо відгук (7.72) представити так:

 

 

xt d f t .

(7.73)

0

 

Величини xt , t , f t , які входять у (7.72), предста-

вимо у вигляді інтегралів Фур'є. Тоді з (7.72) випливає, що фур'є-компоненти цих функцій пов'язані співвідношенням

x f ,

(7.74)

де

 

 

 

dt t ei t .

(7.75)

0

 

Ця функція називається узагальненою сприйнятливістю.

Оскільки функція t

дійсна, з (7.75) одержуємо

 

 

 

.

(7.76)

 

Величина комплексна:

 

i ,

(7.77)

де і – дійсна і уявна частини . Підставляючи

(7.77) у (7.76), одержуємо

267

, . (7.78)

Таким чином, – парна функція частоти, а –

непарна.

У співвідношенні (7.75) змінна дійсна. Однак корисно вийти за межі дійсної осі і розглядати аналітичні властивості узагальненої сприйнятливості в комплексній площині частоти. Оскільки i , одержуємо

ei t ei t e t .

Тоді з (7.75) і обмеженості t при t 0 випливає, що функція регулярна у верхній напівплощині комплек-

сної частоти. Цей результат є наслідком інтегрування в (7.75) за напіввіссю t 0 , тобто принципу причинності. У нижній -напівплощині вираз (7.75) непридатний, тому що має там особливості. Щоб одержати у ниж-

ній напівплощині, необхідно здійснити аналітичне продовження інтеграла (7.75) з області Im 0 в область

Im 0.

З регулярності узагальненої сприйнятливості в області Im 0 (або принципу причинності) випливають важливі

співвідношення між функціями і . Щоб одержати їх, розглянемо функцію

,

де – фіксована дійсна частота. Ця функція регулярна усередині контуру C у площині , показаного на рис. 7.3.

268

Рис. 7.3. Контур інтегрування в інтегралі (7.79)

У силу теореми Коші інтеграл від цієї функції дорівнює нулю:

 

 

 

 

0.

(7.79)

 

d

 

С

Він залишається рівним нулю й у тому випадку, коли радіус малого півкола на рис. 7.3 прямує до нуля, а радіус r великого півкола прямує до нескінченності. Коли r ,

інтеграл за великим півколом онулюється, оскільки 0,

тому прямує до нуля швидше, ніж 1 . Отже, при 0, r з (7.79) одержуємо

 

 

 

 

 

 

 

 

 

 

 

lim

 

d

 

 

 

d

 

d

 

 

0. (7.80)

 

 

 

0

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

Тут d – інтеграл за малим півколом. У цьому інтегралі

С

покладемо

ei

іперейдемо до інтегрування за кутом . Тоді з (7.80)

випливає

 

1

 

 

 

 

 

P

d

,

(7.81)

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

269

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