Учебники / Физика конденсированных сред
.pdf342 |
Глава 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) . |
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. Экстремальность квазиравновесного ансамбля |
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 dτ Sp{Pm |
|
(t)τ Pn |
|
(t)1−τ }. |
(6.44) |
||||
ρ |
ρ |
||||||||
0 |
|
|
|
|
|
|
|
|
Наконец, определяя скалярное двух операторов соотношением
(Pn, Pm)qt = |
1 dτ 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) можно получить