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

nonequilibrium_statistical_operator

.pdf
Скачиваний:
13
Добавлен:
10.02.2015
Размер:
697.54 Кб
Скачать

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

t + Pn

t δFn (t))+ δφ(t) .

(1.50)

n

 

 

 

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

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

t .

(1.51)

n

 

 

Соотношения (1.49), (1.51) можно интерпретировать следующим образом: при записи энтропии роль независимых переменных играют величины Pn t , а

при записи функционала МасьеПланка величины Fn(t).

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

< P >t = −

δφ(t)

.

(1.52)

 

m

δFm (t)

 

 

 

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

S(t) = φ(t)

δφ(t)

Fm (t) .

(1.53)

δF (t)

m

 

 

 

m

 

 

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

F (t) =

δS(t)

.

(1.54)

δ< P

>t

n

 

 

 

n

 

 

 

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

21

φ(t) = S(t)

δ

δS(t)

t Pn (t) t .

(1.55)

n

P (t)

 

 

 

 

n

 

 

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

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

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

S(t)

t

 

ρq (t)

(Pn Fn (t) + Pn Fn (t))+ φ(t)

 

 

q

= Sp

.

(1.56)

 

 

 

 

n

 

 

 

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

φ(t) = −Fm (t) Pm

t .

(1.57)

m

 

 

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

S(t)

t

 

ρq (t)

(

Pn Fn (t) +(Pn

Pn

t )Fn (t)

 

=

 

= Sp

 

 

q

 

 

 

n

 

 

)

 

 

= (Sp(ρq (t)Pn )

Pn

t

 

 

ˆ

.

(1.58)

 

)Fn (t) +Sp(ρq (t)iLS(t))= 0

n

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

ˆ

 

оператор энтропии S(t) коммутируют между собой, и поэтому

 

ˆ

(1.59)

Sp(ρq (t)iLS(t))= 0 .

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

22

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

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

переменных не появляется, и средние значения Pn t операторов Pn медленно

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

1.3.Граничные условия и уравнение Лиувилля для неравновесного статистического оператора

23

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

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

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

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

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

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

lim exp(it1L)ρq

(t +t1

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

(1.60)

t1→−∞

 

t1→−∞

 

Уравнение (1.60) не только позволяет выразить неравновесный статистический оператор ρ(t) через квазиравновесное распределение ρq(t), но и вносит необратимость в поведение величины ρ(t). Действительно,

24

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

 

 

0

 

 

lim f (t) = limε

exp(εt) f (t)dt,

(1.61)

t→−∞

ε→0

 

 

 

 

−∞

 

 

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

limε→0

ε 0

exp(εt1 )ρq (t +t1,t1 )dt1 =limε→0

ε 0

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

(1.62)

 

−∞

 

−∞

 

 

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

 

0

 

 

lim ε

−∞

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

(1.63)

ε→0

 

 

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

 

0

 

 

 

 

 

 

 

 

 

lim ε

exp(εt )exp(iLt )ρ

q

(t +t ,0)dt

(t,0) lim exp(εt )ρ(t +t ,t )

ε→0

1

1

1

1

t

→−∞

1

1 1

 

−∞

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25

 

 

 

 

0

 

 

 

lim exp(εt1 )exp(iLt1 )

 

+iL ρ(t +t1,0)dt1 .

(1.64)

t1

ε→0 −∞

 

 

 

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

0

 

 

lim exp(εt1 )exp(iLt1 )

 

+iL

t1

ε→0 −∞

 

 

Далее,

lim exp(εt1 )ρ(t +t1,t1 )

t1→−∞

ρ(t +t1,0)dt1 . (1.65)

= 0 ,

поскольку величина в этой формуле является конечной и должна стремиться к нулю после выполнения термодинамического предела и вычисления средних. Поэтому выражение (1.64) по существу является определением неравновесного статистического оператора:

 

ρ(t,0) = ε 0

exp(εt1 )exp(iLt1 )ρq (t +t1,0)dt1 .

 

(1.66)

 

−∞

 

 

 

 

 

Найдем теперь уравнение движения, которому удовлетворяет

неравновесный

статистический

оператор

(1.66).

Для

этого

продифференцируем уравнение (1.66) по времени t:

∂ρ(t,0)

= ε 0

exp(εt1 )exp(iLt1 )

d

ρq (t +t1,0)dt1 =(интегрируем по частям)=

t

dt

−∞

 

 

0

= εexp(εt1 )exp(iLt1 )ρq (t +t1,0) 0−∞ −ε (ε+iL)exp(εt1 )exp(iLt1 )ρq (t +t1,0)dt1 =

−∞

26

=(учитываем (1.66)) = εexp(εt )exp(iLt )ρ

q

(t +t ,0)

 

0

−ερ(t,0) iLρ(t,0). (1.67)

 

1

1

1

 

−∞

 

 

 

 

 

 

 

Учитывая, что при t1 → −∞, exp(εt1) 0, получаем уравнение Лиувилля, содержащее бесконечно малый источник в правой части:

∂ρ(t,0) +iLρ(t,0)

t

= −ε(ρ(t,0) −ρq (t,0)) .

(1.68)

Необходимо отметить, что равенство нулю выражения (1.65) выполняется, в чем легко убедиться (см. задача 1.2), если вспомнить формулу (1.62).

Следует сказать несколько слов о смысле бесконечно малых источников в правой части уравнения движения для неравновесного статистического оператора (1.68). Как известно, уравнение Лиувилля (1.20) является обратимым во времени. Вместе с тем мы хорошо знаем, что в реальных системах имеется спонтанное нарушение симметрии динамических уравнений относительно операции обращения времени. Таким образом, в исправленных с учетом второго закона термодинамики в динамических уравнениях должно быть снято вырождение состояний, связанное с симметрией относительно операции обращения времени.

Более последовательно интерпретировать возникновение источников в правой части уравнения (1.68) в духе идеологии квазисредних Н. Н. Боголюбова [21]. Очевидно, что с этих позиций все средние, которые вычисляются при использовании метода неравновесного статистического оператора, являются квазисредними, а член ε(ρ(t,0) − ρq(t,0)), снимающий вырождение уравнения Лиувилля относительно операции обращения времени, в некотором идеализированном виде учитывает контакт системы с термостатом, приводящий к релаксации неравновесного распределения, если систему предоставить эволюционировать самой себе. Тогда величина ε может быть интерпретирована как обратное время релаксации неравновесного распределения к квазиравновесному.

1.4.Интегральные уравнения и теория возмущений для неравновесного статистического оператора

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

27

H = H0 +V , L = L0 + LV .

(1.69)

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

Рассмотрим уравнение Лиувилля с источниками (1.68) для неравновесного статистического оператора

t +iL0 +iLV ρ(t,0)

= −ε(ρ(t,0) −ρq (t,0)) .

(1.70)

Преобразуем уравнение (1.70) в эквивалентное ему интегральное уравнение. Вычитая из правой и левой частей уравнения (1.70) выражение

t +iL0 ρq (t,0) ,

приведем его к следующему виду:

 

+iL

+ε δρ(t,0)

= −

+iL

ρ (t,0) +iL ρ(t,0))

 

,

(1.71)

 

 

 

 

 

 

0

 

 

0

 

q

V

 

 

 

t

 

 

 

t

 

 

 

 

 

 

 

где δρ(t,0) = ρ(t,0) −ρq (t,0).

Вводя

оператор

эволюции

exp(itL0 ) со

свободным гамильтонианом H0 и умножая (1.71) на интегрирующий множитель exp(εt +itL0 ) , представим левую часть (1.71) в виде полной производной по времени

d

e

εt

e

itL

δρ(t,0)

= −e

εt

e

itL

+iL

 

 

 

(t,0)

+iL ρ(t,0))

 

 

 

 

0

 

0

 

 

 

ρ

q

.

(1.72)

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

t

0

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полагая, что

28

lim eεt eitL0 δρ(t,0) = 0,

t→−∞

проинтегрируем уравнение (1.72) по времени от −∞ до t. Имеем

t

 

 

 

 

 

eitL0 δρ(t,0) = − dt1eε(t1 t )eit1L0

 

 

+iL0

ρq (t1

,0) +iLV ρ(t1

,0))

(1.73)

t1

−∞

 

 

 

 

 

 

или, положив t1 t t1 , окончательно

0

 

 

 

 

 

 

 

 

ρ(t,0) = ρq (t,0) dt1eεt1 eit1L0

 

+iL0

ρq (t +t1,0)

+iLV ρ(t +t1

,0))

 

. (1.74)

 

 

 

 

t1

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Если взаимодействие V не входит явно в выражение для операторов Pn, что мы ниже и будем предполагать, то в уравнении (1.74) первые два члена под знаком интеграла зависят от V лишь неявно, через параметры Fn(t). Поэтому можно считать, что эти члены описывают термические возмущения. Третий член под знаком интеграла в (1.74) явно зависит от взаимодействия V, и можно считать, что он описывает механические возмущения.

Уравнение (1.74) можно записать и в другой форме:

ρ(t,0) = ρ0 (t,0) i 0

dt1eεt1 eit1L0 LV ρ(t +t1,0) ,

(1.75)

−∞

 

 

где

0

 

 

 

 

ρ0 (t,0) q (t,0) dt1eεt1 eit1L0

 

 

+iL0

ρq (t +t1,0) =

 

t1

 

−∞

 

 

 

 

=(интегрируем по частям)=ε 0

dt1eεt1 eit1L0 ρq (t +t1,0)

(1.76)

−∞

29

есть статистический оператор, не содержащий явной зависимости от взаимодействия V. Согласно (1.76) выражение ρ0 (t,0) представляет собой инвариантную часть квазиравновесного статистического оператора по отношению к эволюции со свободным гамильтонианом H0. Следует отметить, однако, что ρ0 (t,0) зависит от точных значений функций Fn(t) или < Pn >t , определяемых через обобщенные кинетические уравнения, включающие взаимодействия.

Преобразуем выражение (1.76) для ρ0 (t,0) к другому виду. Используя тождество Кубо и правило дифференцирования операторной экспоненты получим

iL0ρq

(t,0) =

1

[eS (t,0) , H0 ]

= −1

dτe−τS (t,0)

1

[S(t,0), H0 ]e(τ−1)S (t,0)

,

i

i

 

 

 

 

 

 

 

0

 

 

 

 

 

 

ρq (t,0) =

eS (t,0)

= −1

dτe−τS (t,0)

S(t,0) e(τ−1)S (t,0) .

 

 

 

t

 

 

 

 

 

 

 

t

0

 

 

 

t

 

Принимая во внимания эти соотношения и действие оператора exp(itL0 )

exp(itL0 )A = exp(itH0 / )Aexp(itH0 / ) ,

перепишем выражение (1.76) к виду

ρ0 (t,0) = ρq (t,0) + 0

dt1eεt1 1

dτeiH0t1 / e−τS (t+t1 ,0) S(t +t1,0)e(τ−1)S (t+t1,0)eiH0t1 / ,

−∞

0

 

 

 

(1.77)

где

S(t,0) =

S(t,0)

+

1

[S(t,0), H0 ].

t

i

 

 

 

Делая замену

ρ(t,t) = eiH0t / ρ(t,0)eiH0t / , ρ0 (t,t) = eiH0t / ρ0 (t,0)eiH0t / , V (t) = eiH0t / VeiH0t / ,

30

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]