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

Учебники / Физика конденсированных сред

.pdf
Скачиваний:
50
Добавлен:
30.03.2022
Размер:
3.97 Mб
Скачать

§ 2. Экстремальность квазиравновесного ансамбля

341

(6.17), (6.18). Чтобы лучше понять смысл параметров, входящих в определение (6.22), сравним его с каноническим распределением Гиббса

 

1

 

ρ0 =

Z exp{−β(H − ζN )}.

(6.23)

В этом выражении Z−статистическая сумма, ζ− химический потенциал системы, H и N −операторы Гамильтона и числа частиц, β− обратная температура в энергетических единицах.

Из сравнений формул (6.22), (6.23) следует, что равновесное распределение – это распределение с заданным значением энергии и числа частиц, а квазиравновесное – это распределение с заданным значением средних Pn t . Величина φ(t) в выражении (6.22) носит название функционала Масье – Планка и, как и статистическая сумма Z , определяется условием нормировки

φ(t) = ln Sp{exp{− PnFn(t)}}.

(6.24)

n

 

Выбор параметров Pn и функций Fn(t) зависит от конкретной задачи. В частном случае гидродинамического режима, когда измеримыми величинами являются энергия системы, дрейфовый импульс и число частиц, набор операторов Pn и сопряженных им термодинамических функций Fn(t) может быть выбран следующим образом:

 

 

Операторы

H

 

N

 

 

 

P

 

 

 

Термодинамические функции

β(t)

 

β(t)ζ(t)

 

 

 

β(t)mV (t)

 

 

 

 

 

 

 

Здесь P −оператор суммарного импульса частиц системы,

V −их дрейфовая скорость, m−масса.

 

 

 

 

 

Перейдем теперь к построению термодинамики квазиравно-

весного распределения.

 

 

 

 

 

 

Используя определения (6.10) и (6.22), запишем выражение

для энтропии системы

 

 

 

 

 

 

S(t) = φ(t) + Pn tFn(t).

(6.25)

n

342

Глава 6. Метод НСО

Это уравнение можно рассматривать как преобразование Лежандра – переход от одного термодинамического потенциала к другому (от φ(t) к S(t) ) для неравновесной системы. Это становится совершенно очевидным, если произвести вариацию функционала Масье – Планка (6.24):

 

 

 

 

δφ(t) = δ ln Sp{exp{−

PnFn(t)}} =

 

n

 

 

 

Sp{PmδFm(t) ×

= [Sp{exp{−

PnFn(t)}}]1

 

n

 

m

 

× exp{−

 

 

 

 

PnFn(t)}} = − Pm tδFm(t). (6.26)

 

n

 

 

m

Последнее выражение в правой части формулы (6.26) записано

сучетом соотношений (6.6), (6.22), (6.24).

Сдругой стороны, используя определение энтропии (6.25) и явный вид квазиравновесного распределения (6.22), получаем

 

 

δS(t) = δφ(t) + (δ Pn tFn(t) + Pn tδFn(t)).

(6.27)

n

Подставляя в эту формулу значение δφ(t) , определяемое выражением (6.26), получаем

 

 

δS(t) = Fn(t)δ Pn t.

(6.28)

n

Соотношения (6.26), (6.28) можно интерпретировать следующим образом: при записи энтропии роль независимых переменных играют величины Pn t , а при записи функционала Масье Планка величины Fn(t) .

Полученные результаты позволяют обобщить соотношения Гиббса – Гельмгольца на случай неравновесной термодинамики. Вычисляя функциональную производную от функционала Масье – Планка и используя уравнение (6.26), имеем

Pm t =

δφ(t)

(6.29)

δFm(t) .

§ 2. Экстремальность квазиравновесного ансамбля

343

Подставляя этот результат в выражение для энтропии, получаем обобщение соотношений Гиббса – Гельмгольца на случай неравновесной термодинамики:

S(t) = φ(t)

m

δφ(t)

(6.30)

δFm(t) Fm(t).

 

 

 

 

Эта формула выражает энтропию системы через функционал Масье – Планка. Легко можно получить и обратное соотношение. Действительно, из выражения для вариации энтропии получаем

Fn(t) = δ Pn t .

 

 

(6.31)

 

 

δS(t)

 

 

 

 

Тогда формула для энтропии вновь дает

 

φ(t) = S(t)

 

 

δS(t)

Pn t.

(6.32)

 

 

 

 

 

 

 

δ Pn

 

t

 

m

 

 

 

 

 

 

 

 

Отличие этих соотношений от их равновесных аналогов сводится только к замене частных производных на функциональные.

Для понимания смысла квазиравновесного распределения ρ(t) очень важно выяснить, можно ли использовать это распределение для описания неравновесных процессов?

Вычислим производство энтропии в квазиравновесном состоянии. Усредняя оператор производства энтропии (6.13) по

квазиравновесному распределению, получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

˙

t

 

 

˙

˙

 

 

 

˙

(6.33)

S(t) q = Sp(t)[φ(t) +

PnFn(t) +

PnFn(t)]}.

 

 

 

 

 

 

 

n

 

 

 

n

 

 

Учитывая соотношение (6.26), получаем

 

 

 

 

 

 

 

 

 

 

t ˙

 

 

 

 

 

 

 

 

 

˙

 

 

 

 

 

 

 

 

 

φ(t) =

Pm Fm(t).

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

Подставляя этот результат в выражение (6.33), находим

 

˙

t

 

 

 

 

t

˙

˙

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

S(t) q = Sp(t)

[(Pn − Pn )Fn(t) + PnFn(t)]} =

 

= (Sp{

 

 

 

 

 

n

 

 

 

 

 

 

 

(t)Pn} − Pn t)F˙n(t) + Sp{

 

(t)iLSˆ(t)} = 0.

(6.34)

ρ

ρ

n

344 Глава 6. Метод НСО

При выводе последнего соотношения мы учли, что ρ(t) и опе-

ратор энтропии

ˆ

коммутируют между собой, и поэтому

S(t)

 

 

 

 

ˆ

 

 

 

 

 

 

Sp(t)iL S(t)} = 0.

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

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

В завершение параграфа найдем связь между вторыми функциональными производными от потенциалов S(t) и φ(t) и корреляционными функциями по квазиравновесному состоянию:

 

 

 

δ Pm t

=

 

 

δ2φ(t)

=

 

 

 

 

δFn(t)

δFn(t)δFm(t)

 

 

 

 

 

 

 

=

δ

 

 

 

 

 

 

 

Sp{Pm exp{−[φ(t) +

PkFk(t)]}}.

(6.35)

δFn(t)

 

 

 

 

 

 

 

k

 

 

Сделаем небольшое математическое отступление и вычислим производную по параметру от операторной экспоненты. Рассмотрим вначале более простой вопрос о разложении экспоненты

exp{(A + B)t}

в степенной ряд. Здесь A и B – не коммутирующие между собой операторы, а t−некоторый параметр. Введем обозначения:

§ 2. Экстремальность квазиравновесного ансамбля

345

exp{(A + B)t}, ≡ D(t) ≡ G(t) exp(At).

Составим уравнение движения для функции D(t) :

dD(t)

= (A + B)D(t) =

dG(t)

exp(At) + G(t)A exp(At). (6.36)

dt

 

 

dt

Оператор A коммутирует с операторной экспонентой exp(At) . Поэтому второе равенство в выражении (2.36) можно записать в виде

(A + B)D(t) = dG(t) exp(At) + D(t)A. dt

A + B тоже коммутирует с D(t) и поэтому

(A + B)D(t) = D(t)(A + B).

Подставляя этот результат в формулу (6.36) и сокращая одинаковые члены, получаем

D(t)B = dG(t) exp(At), dt

или

dG(t)

= exp{(A + B)t}B exp{−At}.

(6.37)

dt

Учитывая, что exp{(A + B)t} = G(t) exp(At) и используя последнее уравнение, получаем

d ln G(t1) = exp(At1)B exp(−At1) dt1.

Интегрируя это дифференциальное уравнение с учетом гранич-

ного условия G(0) = 1,

ln G(0) = 0, получаем

 

G(t) = exp

{

exp()B exp(

}

(6.38)

t

 

) dλ ,

0

t

exp{(A + B)t} = exp{ exp()B exp(−Aλ)dλ} exp(At).

0

346

Глава 6. Метод НСО

Если оператор B мал (малость оператора понимается как малость соответствующих матричных элементов) и можно обойтись первыми членами разложения, то вместо (6.38) получаем

G(t) = 1 + t

exp()B exp(

) ;

exp (A + B)t

}

=

0

 

 

 

{

 

 

 

 

 

 

 

 

= exp(At) + t

exp()B exp(−Aλ)exp(At).

 

(6.39)

 

0

 

 

 

 

 

 

На основании этой формулы выведем правило дифференцирования операторной экспоненты по параметру. Используя определение производной, имеем

 

 

 

 

d

exp(P1λ1 + P2λ2) =

 

 

 

 

 

 

 

1

 

2

=

lim

[exp(P1λ1 +P2λ2 +P2 λ2) exp(P1λ1 +P2λ2)].

 

λ2

λ20

(6.40)

Считая, что P2 λ2 есть малый оператор, и полагая t = 1 , на основании формулы (6.39) получаем

exp(P1λ1 + P2λ2 + P2

λ2) = exp(P1λ1 + P2λ2) +

+ 1 exp[(P1λ1 + P2λ2)λ]P2

λ2 exp[(P1λ1 + P2λ2)λ] dλ.

0

 

(6.41)

С учетом этого результата в итоге имеем

d

2

exp(P1λ1 + P2λ2) =

1

=exp[(P1λ1 + P2λ2)λ]P2 exp[(P1λ1 + P2λ2)(λ − 1)] dλ.

0

(6.42)

§ 2. Экстремальность квазиравновесного ансамбля

347

Возвратимся снова к формуле (6.35). С помощью выражения (6.42), найдем функциональную производную:

 

 

 

 

 

δ

 

 

 

 

 

 

 

 

 

 

exp{−

PkFk(t)} =

 

 

 

δFn(t)

 

 

 

 

 

 

 

k

 

 

 

 

 

= 1 exp[

 

 

 

 

 

 

 

 

 

 

PkFk(t)τ ]Pn exp[

k

 

PkFk(t)(τ − 1)] dτ. (6.43)

0

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее, действуя аналогично с учетом того, что

 

exp(−φ(t)) = [Sp{exp(

k

PkFk(t))}]1,

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ

exp(

φ(t)) =

Sp{Pn exp(

k PkFk(t))}

.

 

δFn(t)

 

 

 

 

 

 

 

[Sp exp(

 

k

PkFk(t)) ]2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{

 

}

 

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

 

δ

 

 

 

Sp{Pm

 

(t)} =

 

 

 

ρ

 

 

δFn(t)

 

= Pn t Pm t 1 Sp{Pm

 

(t)τ Pn

 

(t)1−τ }.

(6.44)

ρ

ρ

0

 

 

 

 

 

 

 

 

Наконец, определяя скалярное двух операторов соотношением

(Pn, Pm)qt =

1 Sp (Pm

 

 

Pm

 

t)ρ(t)τ (Pn

 

Pn

 

t)ρ(t)1−τ

,

 

{

 

 

 

 

 

 

 

 

 

 

}

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ Pm qt

 

 

 

 

 

t

 

 

 

 

 

 

 

получаем

 

 

= (Pm, Pn)q .

 

 

 

 

 

 

(6.45)

 

δFn(t)

 

 

 

 

 

 

Подведем некоторые итоги. Исходя из принципа экстремальности информационной энтропии построено выражение

348

Глава 6. Метод НСО

для квазиравновесного статистического оператора (6.22). Смысл этого распределения состоит в том, что оно описывает только что приготовленный ансамбль неравновесных систем, в котором еще не началась эволюция и не развились потоки.

Ключевым для понимания метода НСО является соотношение (6.6), устанавливающее равенство средних значений базисных операторов Pn , вычисленных с использованием неравновесного и квазиравновесного распределений. Истолковать это соотношение можно следующим образом. К моменту времени, когда сформировался квазиравновесный ансамбль, единственным набором величин, измеримых в неравновесной системе, был набор переменных Pn . В дальнейшем эволюция системы происходит так, что новых медленно меняющихся динамических переменных не появляется, и средние значения Pn t операторов Pn медленно эволюционируют благодаря зависимости от времени сопряженных термодинамических сил Fn(t) .

Что касается термодинамических сил Fn(t) , то они формируются в ходе реальной эволюции системы и будут зависеть от неравновесных процессов, протекающих в системе. Нахождение термодинамических сил Fn(t) будет темой подробного обсуждения в параграфе, посвященном линейным релаксационным уравнениям в методе НСО.

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

§ 3. Граничные условия и уравнение Лиувилля для НСО

Рассмотрим неравновесную систему, состояние которой на достаточно больших временах описывается набором макроскопических переменных Pn t . Как уже неоднократно отмечалось, это означает, что только эти величины являются измеримыми в данной системе и что сделанное предположение не нарушает общности рассмотрения. Чаще всего набор величин

§ 3. Граничные условия и уравнение Лиувилля для НСО 349

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

Будем предполагать, что в момент времени t0 , который для удобства будет отнесен на отрицательную бесконечность (конечно, имеется в виду «физическая бесконечность», т. е. времена, значительно большие, нежели некоторое характерное для данной системы время размешивания, за которое «вымирают» несущественные для дальнейшей эволюции корреляции), приготовлен квазиравновесный ансамбль систем, описываемый квазиравновесным распределением ρ(t).

Сформулируем начальное условие для неравновесного статистического оператора ρ(t) . Будем полагать, что в момент времени t0 неравновесный и квазиравновесный статистические операторы совпадают.

Сформулируем теперь условие, позволяющее записать неравновесный статистический оператор в виде некоторого функционала от квазиравновесного распределения. Мы уже отмечали, что квазиравновесный статистический оператор ρ(t) не удовлетворяет уравнению Лиувилля и под действием оператора эволюции будет трансформироваться в отличие от неравновесного распределения ρ(t) , которое является интегралом движения.

Будем считать, что если приготовить квазиравновесное распределение, а затем предоставить системе возможность эволюционировать, то квазиравновесное распределение ρ(t) через некоторое время порядка времени размешивания трансформируется в неравновесное распределение ρ(t) .

На языке математики это последнее условие и сформулированное выше граничное условие для НСО с учетом введенных ранее определений (6.2) – (6.4) можно записать в виде

lim exp(it1L)

 

(t + t1, 0) =

lim exp(it1L)ρ(t + t1, 0). (6.46)

ρ

t1→−∞

t1→−∞

Уравнение (6.46) не только позволяет выразить НСО ρ(t) через квазиравновесное распределение ρ(t) , но и вносит необратимость в поведение величины ρ(t) . Действительно, достаточно в этом уравнении устремить t1 +, чтобы теория описывала

350

Глава 6. Метод НСО

не возрастание, а убывание энтропии в системе. Причина этого понятна. В уравнении (6.46) квазиравновесное распределение, сформированное в момент времени t0 = −∞, в ходе эволюции трансформируется в неравновесное распределение при t > t0 . Иначе говоря, направление спонтанно текущего процесса задано и меньшему значению времени соответствует более упорядоченное состояние. Если положить t0 = +, то система с течением времени будет переходить из менее упорядоченного в более упорядоченное состояние, что и соответствует уменьшению энтропии с течением времени. Применяя теорему Абеля, согласно которой

 

0

 

lim f (t) = lim

exp( t)f (t)dt,

(6.47)

t→−∞ →0 −∞

если этот предел существует, перепишем уравнение(6.46) в следующем виде:

lim 0

exp( t1)

 

(t + t1, t1)dt1

= lim 0

exp( t1)ρ(t + t1, t1)dt1.

ρ

0 −∞

 

 

 

0 −∞

 

(6.48) Уравнение (6.48) допускает интересную интерпретацию. По существу, формула (6.48) утверждает, что сглаженные (усредненные) по достаточно большому промежутку времени статистические операторы ρ(t + t1, t1) и ρ(t + t1, t1) равны между собой. Часто сглаживание, определяемое формулой (6.48), называют взятием инвариантной части. Очевидно, что ρ(t + t1, t1) =

= ρ(t) , поэтому

lim

0

exp( t1)

 

(t + t1, t1)dt1 = ρ(t).

(6.49)

ρ

0

−∞

 

 

 

 

Из уравнений (6.48), (6.49) следует, что в ходе эволюции квазиравновесное распределение трансформируется в неравновесное распределение. В этом, собственно, и состоит физический смысл уравнения (6.48). Результат (6.49) можно получить