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

Классическая статистика ч2 - Сотский А.Б

..pdf
Скачиваний:
25
Добавлен:
18.04.2015
Размер:
955.53 Кб
Скачать

E E0 (a) = δE → 0 .

(37.17)

Тогда

 

Γ(E, a) = Ω(E , aE + OE2 ) = Ω(E, aE + OE2 ) .

(37.18)

0

 

Подставив (37.18) в (37.16) и опустив величины второго порядка малости,

находим: δE = kT .

Для макроскопических

систем, энергия которых

E ~ NkT отношение

δE ~

 

 

 

 

1

→ 0 .

(37.19)

 

 

 

E

N

энергии dE физически

Выражение (37.19)

означает,

что интервал

бесконечно мал, что оправдывает представление (37.18).

Итак, мы нашли, что энтропия изолированной равновесной системы может быть взята в виде (37.15). Выражение (37.15) называют формулой Больцмана. Согласно этой формуле, энтропия названной системы пропорциональна логарифму фазового объема в Γ − пространстве, заключенного между двумя физически бесконечно близкими гиперповерхностями постоянной энергии (37.9).

По поводу постоянной d в (37.15) можно сделать следующее замечание. Из квантовой механики известно, что классические понятия координаты и импульса можно ввести только в квазиклассическом приближении. При этом, в соответствии с соотношением неопределенностей Гейзенберга DpDq ³ h , где h - постоянная Планка,

минимальный объем фазовой ячейки в Γ − пространстве одной частицы,

совершающей трехмерное движение, равен h3 ,

а в Γ − пространстве N

таких частиц он равен h3 N . Таким образом,

величина h3 N

является

естественной единицей измерения фазового

объема. Выбрав

в (37.15)

d = h3 , мы можем отождествить отношение

Γ / h3 N с числом различных

состояний макроскопической системы в фазовом объеме Γ . Тогда энтропию изолированной системы

 

Γ

 

 

S = k ln

 

 

 

(37.20)

 

3 N

h

 

 

 

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

H ( X , a) = E .

Следует отметить, что статистические определения энтропии (28.10) и (37.20) не совпадают. Дело в том, что эти выражения относятся к физически различным системам: к системе в термостате, энергия которой может флуктуировать и к замкнутой системе со строго фиксированной энергией. Соответствие формул (28.10) и (37.19) мы рассмотрим в квантовой статистике.

51

§ 38. Формула Эйнштейна для вероятности макроскопических параметров

Выделим некоторую малую подсистему в большой изолированной равновесной системе с фиксированной энергией E . Найдем вероятность того, что некоторый макроскопический параметр подсистемы F ( X , a)

находится в диапазоне

F < F ( X , a) < F + dF .

(38.1)

Неравенству (38.1) соответствуют системы фазового ансамбля, изображающие точки которых в Γ − пространстве заключены между гиперповерхностями

F ( X , a) = F ,

F ( X , a) = F + dF .

(38.2)

С другой стороны, эти точки должны принадлежать гиперповерхности постоянной энергии

H ( X , a) = E .

(38.3)

Таким образом, изображающие точки рассматриваемых систем располагаются в 6N −1 − мерной полосе σ , которая находится на гиперповерхности (38.3) и заключена между линиями, образованными пересечением гиперповерхностей (38.2) и (38.3). Качественно эта ситуация

σ

H(x,a)=E

F(x,a)=F F(x,a)=F+dF

Рис. 38.1

изображена на рис. 38.1, где полоса σ заштрихована. Очевидно, что площадь dσ полосы σ будет зависеть от значения F , которое входит в уравнения (38.2), то есть dσ = dσ (F ) .

Поскольку большая система описывается микроканоническим распределением (37.1), фазовая плотность вероятности w( X , E, a)

одинакова для всех точек гиперповерхности (38.3). Тогда из

52

геометрических соображений очевидно, что вероятность dυ события, состоящего в попадании изображающей точки ансамбля в указанную выше полосу, будет пропорциональна ее площади:

dυ ~ dσ (F ) .

(38.4)

Ввиду того, что гиперповерхности (37.9) при условии (37.17) физически бесконечно близки одна к другой, (38.4) можно записать также в виде

dυ ~ ΔΓ(F ) , (38.5) h3N

где ΔΓ(F ) − объем фазового пространства, примыкающий к площадке dσ

и заключенный между гиперповерхностями (37.9).

Сопоставим состояниям большой системы, которым соответствуют изображающие точки, расположенные в пределах полосы σ , энтропию

 

ΔΓ(F )

 

S (F ) = k ln

h

3 N

.

(38.6)

 

 

 

 

Выражение (38.6) является естественным обобщением определения энтропии равновесного ансамбля (37.15) на случай неравновесного ансамбля, изображающим точкам которого доступна не полная гиперповерхность H ( X , a) = E , а лишь расположенная на ней полоса σ .

Определение энтропии (38.6) называют постулатом Больцмана. Из (38.5) и (38.6) заключаем, что

 

S (F )

dυ ~ exp

 

.

 

 

k

Вероятность dυ , очевидно, будет пропорциональна и дифференциалу dF . Таким образом,

 

S (F )

 

dυ = w(F )dF = C exp

 

dF ,

(38.7)

 

 

k

 

где w(F ) − плотность вероятности для параметра F ,

 

 

S (F )

−1

C =

exp

 

dF

 

 

 

k

 

− нормировочная постоянная.

Если рассмотреть не один, а несколько макроскопических параметров подсистемы F1, F2 ,..., то рассуждения, подобные

представленным выше, позволяют записать вероятность обнаружения данных параметров в диапазонах (F1, F1 + dF1 ),(F2 , F2 + dF2 ),... в виде

 

 

S (F1, F2

,...)

 

 

w(F1, F2

,...)dF1dF2

... = C exp

 

 

 

 

dF1dF2

...,

(38.8)

 

 

 

 

 

 

 

k

 

 

 

 

 

где w(F1, F2 ,...) − совместная плотность

 

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

случайных

величин F1, F2 ,...,

 

 

 

 

 

 

 

 

 

53

 

dF1

dF2

 

S (F , F ,...) −1

 

C =

...exp

1

2

 

(38.9)

 

 

 

 

 

 

k

 

 

 

− нормировочная постоянная.

Выражение (38.8) называется формулой Эйнштейна для вероятности макроскопических параметров.

§ 39. Термодинамическая теория флуктуаций

Формулы (38.7) и (38.8) получены исходя из микроканонического распределения плотности вероятности для фазового ансамбля равновесных систем. Поскольку ансамбль рассматривается в фиксированный момент времени (см. § 23), эти формулы дают вероятность мгновенных значений макроскопических параметров малой подсистемы, расположенной внутри большой изолированной равновесной системы. Но сейчас нас интересуют не мгновенные, а термодинамические, то есть средние по времени значения этих параметров, которые и наблюдаются экспериментально.

Рассмотрим зависимость некоторого макроскопического параметра подсистемы F от времени. График этой зависимости схематически

представлен на рис. 39.1. Измеряемые значения

F (t) даются интегралом

(23.3). Но в зависимости от протяженности интервала усреднения τ

этот

интеграл допускает различную интерпретацию.

 

 

Введем в рассмотрение времена релаксации для малой подсистемы и

для большой системы, равные τ r

и Τr , соответственно. Из опыта известно,

F

 

 

 

_

 

 

 

_ F

 

 

 

F(t)

 

 

 

T

τ

 

 

 

 

t

 

 

Рис. 39.1

 

 

 

54

 

 

что время релаксации системы возрастает пропорционально ее размерам. Поэтому выполняется неравенство Τr >> τ r .

Если τ >> Τr , то усреднение (23.3) даст константу F , не зависящую

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

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

τ r << τ ,

(39.1)

τ << T << Τr .

(39.2)

Здесь T − характерный временной масштаб ощутимого изменения параметра F (рис. 39.1). При условии (39.1) подсистема в течение времени регистрации τ может рассматриваться как равновесная, а в силу

неравенств (39.2) среднюю величину F (t) вида (23.3) можно отождествить

с мгновенным значением F и трактовать это значение как зависящий от времени термодинамический параметр подсистемы. Эволюцию этого параметра можно рассматривать как равновесный случайный процесс с характерным временем T , вызванный воздействием на подсистему со стороны остальной части большой системы.

Таким образом, при условиях (39.1), (39.2) мы можем использовать для вероятности термодинамических параметров подсистемы формулу Эйнштейна (38.8) и считать, что приращения этих параметров подчиняются дифференциальным соотношениям равновесной термодинамики. Это открывает возможность построения термодинамической теории флуктуаций.

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

Обозначим через V1, E1, S1 (E1,V1 ) объем, энергию и энтропию

подсистемы, а через V2 , E2 , S2 (E2 ,V2 )

аналогичные величины для

термостата. Будем считать выполненными равенства

E = E1 + E2 ,

V = V1 + V2 ,

где E и V − энергия и объем полной системы, значения которых фиксированы. Заметим, что в силу аддитивности энтропии

S = S1 (E1,V1 ) + S2 (E2 ,V2 ) = S1 (E1,V1 ) + S2 (E E1,V V1 ) = S (E1,V1 ) . (39.3)

Пусть энергия и объем подсистемы испытывают флуктуации, то есть отклоняются от своих средних, или равновесных значений E1 и V1 ,

достигаемых в пределе при τ → ∞ в (23.3). Обозначим приращения этих параметров как

E1 = E1

 

 

 

 

(39.4)

E1, V1 = V1 V1 .

55

Разложим энтропию (39.3) в ряд Тейлора по степеням приращений (39.4). Считая эти приращения малыми, ограничимся в этом разложении величинами нулевого, первого и второго порядков:

S = S (

 

 

 

) + L (

E ,

V ) + L (

E2

,

V 2

,

E

V ) ,

(39.5)

E

,V

1

1

1

1

1

2

1

 

1

 

1

1

 

где L1 и L2 − некоторые линейные функции указанных аргументов.

Поскольку в состоянии равновесия энтропия большой системы достигает максимума (см. § 16),

 

 

 

 

 

 

 

 

 

L 1=

 

 

 

S

 

 

E1

+

 

S

 

 

V1 = 0 ,

 

 

 

 

 

 

 

 

(39.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, E1 =

 

 

 

 

 

 

 

 

 

 

 

 

 

где производные вычисляются при V1 = V1

E1 .

 

 

 

 

 

 

 

 

 

 

Для слагаемого L 2 в (39.5) имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2 =

1

 

2 S

 

 

 

 

 

2

 

+

1 2 S

 

 

 

 

 

2

+

 

2S

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

E2

 

 

 

 

E1

 

 

V 2

 

 

 

V1

 

 

 

 

 

 

E1 V1

 

 

 

 

 

2

 

 

 

 

2

 

 

E V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

V

 

 

 

 

 

 

 

 

 

 

 

 

1

 

E

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

(39.7)

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2 S2

 

 

 

 

 

 

1

 

2 S2

 

 

 

 

 

 

 

 

 

 

2 S2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

E2

 

 

E1

+

 

 

 

V 2

 

 

 

V1

 

+

 

 

 

 

 

E1 V1,

 

 

 

 

 

 

 

2

 

2

 

 

 

E V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

V

 

 

 

 

 

 

 

 

 

 

 

2

E

2

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где производные

 

вычисляются

 

 

при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V1 = V1

,

 

 

E1 = E1

 

и

учтено,

что

V2 = −

 

V1 ,

 

E2 = − E1 . Сумму трех первых слагаемых в (39.7) запишем в

виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψ =

1

 

 

 

 

2 S

 

E +

 

 

2 S

 

 

 

 

+

 

1

 

 

 

 

2 S

 

 

 

 

+

2 S

 

 

 

 

 

E

 

 

1

 

 

 

 

 

 

 

1

 

V

 

 

 

 

V

 

 

 

1

 

 

V

1

E

.

 

 

 

E2

 

E V

 

 

 

 

 

V

 

 

 

E V

 

2

 

 

1

 

 

 

 

1

 

 

 

 

 

1

 

 

 

2

 

1

2

 

 

1

 

1

 

 

 

 

 

 

 

 

 

1

V1

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

E1

 

 

 

 

1 1

 

 

Поскольку малая подсистема может рассматриваться как равновесная, в соответствии с основным равенством термодинамики имеем

S

1

 

S

 

 

P

 

 

1

 

 

 

 

 

1

 

1

 

 

 

 

=

 

,

 

V

 

=

 

,

E

T

T

 

1

V

1

 

1

E

1

 

 

 

1

 

 

 

 

 

1

 

 

 

где P

и T − давление и температура в подсистеме. Тогда приращение

1

1

 

производной (S1

E1 )V , вызванное приращениями параметров E1 и V1 ,

 

 

1

составит

где

где

 

S1

 

 

2 S1

 

 

2S1

 

 

1

 

 

T1

 

 

 

=

E2

 

E1 +

 

V1 =

 

 

 

= −

T 2 .

E

E V

T

 

1

V

 

1

V

 

1 1

 

 

1

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

T1 = T1 T , T − равновесное значение температуры. Аналогично

S

2 S

 

 

2 S

 

P

T P P T

 

 

1

 

=

1

 

 

V1 +

1

E1 =

1

=

1

1

,

 

 

 

V 2

 

 

 

 

 

 

 

 

V

E V

T

T 2

 

 

1

E

 

1

E

 

1 1

 

1

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

P = P P , P − равновесное значение давления. Следовательно,

1

1

56

 

E T

V T P P T

 

 

T ( E + P V )

 

P V

Ψ = −

1

1

+

1

 

 

 

1

1

= −

1

1

1

+

1 1 .

2T 2

 

2

 

 

 

T 2

 

 

2T 2

 

 

 

 

 

 

 

 

 

 

 

 

2T

Но в малой подсистеме

S1 = (

 

E1 + P V1 ) /T . Поэтому

 

 

 

 

 

 

 

Y =

 

P

V

T S

 

 

 

 

 

 

 

 

 

 

1

1

1

1

.

 

 

 

(39.8)

 

 

 

 

 

 

 

2T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оставшиеся три слагаемых в (39.7) имеют порядок O ( N1 / N2 ) , где

N1 и N2 − числа

частиц в

 

подсистеме

и

термостате,

соответственно.

Действительно, рассмотрим, например, первое из указанных слагаемых. Поскольку приращения (39.4) малы, термостат близок к состоянию равновесия. Поэтому приближенно верна формула равновесной

термодинамики

(S2

E2 )V

= 1 T2 ,

 

 

где

 

 

T2 − температура

термостата.

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2 S

 

 

 

E12 =

1

 

 

 

1

 

 

 

 

 

 

 

 

 

E22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E12 .

 

 

 

 

2

 

2

E

T

 

 

 

 

 

 

 

2

 

V

 

 

 

 

 

 

 

 

2

 

 

2

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

=

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E2 T2

 

N2 Eɶ2 T2

 

 

 

где Eɶ

-удельная энергия, производная по которой от 1/T

конечна. За

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

оценку

квадрата

приращения

энергии DE2

 

 

 

примем

дисперсию этой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

величины D( E1 ) . Поскольку величина E1

 

аддитивна относительно числа

частиц в подсистеме, в соответствии с §

24

 

имеем D(

E1 ) ~ N1 . Таким

образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

S22

 

 

 

 

 

 

 

 

 

 

N1

 

 

 

 

 

 

 

 

 

 

 

DE12 = O

 

.

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E2 V

 

 

 

 

 

 

 

 

 

N2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогичные оценки могут быть получены и для двух оставшихся из трех рассматриваемых слагаемых. Таким образом, при N1 / N2 → 0 данные

слагаемые пренебрежимо малы по сравнению с Ψ .

Итак, сохранив в (39.7) только сумму трех первых слагаемых вида (39.8), для экспоненциального множителя в выражении (38.8) имеем

 

S

 

DPDV - DT DS

 

S (

E

,

V

)

exp

 

 

= exp

 

 

× exp

 

 

 

 

 

.

 

2kT

 

k

 

 

k

 

 

 

 

 

 

 

Здесь для краткости индексы «1» опущены и подразумевается, что все приращения величин относятся к малой подсистеме, которую в дальнейшем будем именовать просто системой. Тогда для рассматриваемой моновариантной системы с двумя степенями свободы выражение (38.8) приобретает вид

57

w(F , F )dF dF = C exp

 

P V

T S

dF dF ,

(39.9)

 

 

 

1 2 1 2

2kT

 

 

1 2

 

 

 

 

 

 

 

где множитель exp S (

 

 

 

) / k

 

включен в нормировочную постоянную C ,

E

,V

 

 

 

 

 

 

 

 

 

 

которая в соответствии с (38.9)

равна

P V T S −1

 

 

 

 

 

C =

d (

 

F1 ) d ( F2 )exp

 

.

 

2kT

−∞

 

 

−∞

 

 

 

Здесь F1 = F1 F1 , F2 = F2 F2 , величины с чертой обозначают средние

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

Выражение (39.9) составляет основу для вычисления флуктуаций основных термодинамических величин. Проиллюстрируем это на примерах.

Выберем в качестве независимых термодинамических параметров

системы объем и температуру ( F1 = V ,

 

F2 = T ). Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P =

 

P

 

 

V

+

P

 

T ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V T

 

 

T V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

S

 

 

 

 

 

 

 

 

P

 

 

 

 

CV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S =

 

 

 

 

 

V +

 

 

 

 

 

T =

 

 

 

 

 

 

 

 

V +

 

T .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V T

 

 

 

 

 

T

 

V

 

 

 

 

 

 

T V

 

 

 

 

 

 

Согласно (39.9),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V 2 P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(V ,T ) = C exp

 

 

 

 

 

2

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

V

,

 

 

 

 

(39.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2kT V

T

 

 

 

 

 

 

 

 

 

2kT

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CV

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C = exp −

 

 

 

 

 

 

 

T

 

d (

 

T )

×

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2kT

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(39.11)

 

 

1

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

CV (P V )

 

 

 

×

exp

 

 

 

 

 

 

 

 

 

 

 

V

 

 

d (

V )

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

.

 

 

 

 

 

 

 

 

 

 

 

 

 

kT

 

 

T

 

 

 

−∞

2kT V T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пользуясь

 

 

выражениями

 

(39.10),

 

 

 

(39.11),

найдем моменты

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V 2 ,

T 2 . Выполнив

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V T ,

 

 

 

 

 

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

 

интегрирование в

определении

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

,

 

= d ( V ) d ( T ) w(V ,T )(

V T , V 2 , T 2 ),

 

 

 

 

 

 

 

V 2

T 2

 

 

V T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(39.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

58

 

 

 

V

 

V

2

,

 

= −kT

 

 

 

 

P T

 

2

T 2 = kT . CV

(39.13)

(39.14)

Исходя из выражений (39.12) – (39.14), легко рассчитать дисперсию энергии системы. Для этого следует воспользоваться установленным в курсе термодинамики соотношением

E =

 

E

V +

 

E

T =

T

 

P

P

 

V + C T .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

V

 

 

T

 

 

 

T

 

 

 

 

 

 

T

 

 

V

 

 

 

V

 

 

Усредняя квадрат правой части в (39.15) и учитывая (39.12) находим

 

= −kT

 

V

T

 

P

P

2

+ C kT 2 .

E2

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

P

T

 

T

 

 

 

 

 

 

 

 

 

V

 

 

(39.15)

– (39.14),

(39.16)

 

 

 

 

 

 

 

 

Для определения момента

V P запишем выражение

 

 

P

 

P

 

 

P =

 

 

V +

 

 

T .

(39.17)

 

 

 

V T

 

T V

 

 

Умножая обе части (39.17) на V , усредняя результат и пользуясь выражениями (39.12), (39.13), получаем

 

 

 

V

P

 

= −kT .

 

 

(39.18)

Аналогично, умножая обе части равенства (39.17) на

T и пользуясь

(39.12), (39.14), находим

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

kT

2

 

 

 

 

 

 

 

 

 

 

 

P T =

 

 

 

 

 

.

(39.19)

 

 

 

 

 

 

 

 

 

 

 

T

C

 

 

 

 

 

 

 

 

 

 

 

V

V

 

 

 

Из (39.12) – (39.15)

нетрудно также заключить, что

 

 

 

 

 

 

= kT 2 .

 

 

(39.20)

 

 

 

E

T

 

 

Корреляционные моменты, содержащие энтропию, проще всего рассчитать, если в качестве независимых термодинамических параметров выбрать давление и энтропию ( F1 = P , F2 = S ). Тогда

 

V

 

P +

V

 

 

V

 

 

 

T

S ,

V =

 

 

 

 

 

S =

 

 

 

P +

 

 

 

P S

 

 

 

S P

 

 

P S

 

 

P S

 

T

 

 

 

 

T

 

 

 

T

 

 

T

 

T =

 

P

+

 

 

S =

 

 

P +

 

 

S .

 

CP

 

P

S

 

 

 

S P

 

 

P

S

 

 

В результате плотность вероятности в формуле (39.9) приобретает вид

 

 

 

 

1

 

V

 

 

2

 

 

2

1

 

 

w(P, S ) = C exp

 

 

 

 

 

P

 

S

 

 

,

(39.21)

 

 

 

 

 

 

 

 

 

2kT P

S

 

 

 

 

2k CP

 

где

59

 

 

1

 

V

 

 

 

 

2

 

 

 

−1

 

 

 

2

1

 

−1

 

C =

exp

 

 

 

 

 

 

P

 

d ( P)

 

 

exp −

S

 

 

d (

S )

=

 

 

 

 

 

 

 

−∞

 

2kT

 

P S

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

2k CP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2π k

 

 

CPT

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(V P)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисляя интегралы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

,

P2

,

S 2

) = d ( P) d ( S ) w(P, S )(

P S, P2 , S 2 ),

 

 

P S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(39.22)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P2 = −

 

 

 

 

 

,

 

 

 

 

 

(39.23)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(V P)

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= k C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S 2

P

.

 

 

 

 

 

 

 

 

 

(39.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

отметить, что согласно (39.12) и (39.22) случайные величины

T ,

V и

P,

S статистически независимы. Этого нельзя сказать о парах величин

V ,

P ,

P,

T и

 

E, T (см. выражения (39.18) – (39.20)). Существенно

также, что поскольку параметры V ,CV ,CP , S, E являются экстенсивными (они пропорциональны числу частиц в системе N1 ), а T , P − интенсивными

(они не зависят от N1 ), из выражений (39.13), (39.14), (39.16), (39.23), (39.24) можно заключить, что флуктуации (корни квадратные из

дисперсий)

экстенсивных

параметров пропорциональны

N1 ,

 

а

 

 

 

 

 

флуктуации

интенсивных

параметров обратно пропорциональны

N1 .

Еще одним следствием рассмотренной теории флуктуаций являются так называемые условия устойчивости термодинамических систем. Для их вывода заметим, что при флуктуациях параметров малой подсистемы

энтропия большой замкнутой системы

испытывает отклонение S от

равновесного значения, которое согласно (39.5) и (39.8) равно

S =

P

V

T

S

1

1

1

1

.

 

2T

 

 

 

 

 

 

 

Поскольку равновесное значение энтропии большой замкнутой системы

максимально (см. §16), выполняется неравенство

S ≤ 0 . Следовательно,

P V T S ≤ 0 ,

(39.25)

где все величины относятся к малой подсистеме, а индекс «1» у них для краткости записи опущен.

Выберем в качестве независимых термодинамических параметров малой подсистемы энтропию S и объем V . Тогда

60