Файл: Кудрин, Л. П. Статистическая физика плазмы.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 156

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

где статистическая

сумма большого канонического

ансамбля

Zn =

2

(ехР(РН'ЛО 2 ехР (— £*л/Р)]-

(7.8)

 

N

п

 

Здесь данное состояние системы определяется не только набо­ ром квантовых чисел к, но и числом частиц, которое является переменной величиной.

Более компактно статистические суммы ZG и ZN можно записать через след соответствующих операторов *:

Zq =

Sp [exp (—рЯ)]; ZN =

Sp [exp )P (рЯ — Я))].

(7.9)

Такая форма

записи отражает чрезвычайно

 

важное

обстоя­

тельство: д л я

в ы ч и с л е н и я

с л е д а

м о ж н о

п о л ь ­

з о в а т ь с я л ю б о й п о л н о й

с и с т е м о й

в о л н о в ы х

ф у н к ц и й .

 

 

 

 

 

 

Введем матрицу плотности

 

 

 

 

р4= ехр[Р(рЯ — Я)] = exp [Р (рЯ — Я,,)] 5(Р).

(7.10)

Разобъем гамильтониан системы Я на две части:

 

 

 

Я = Я0-|-Я;,

 

 

(7.11)

где Яв — оператор, описывающий

движение свободных

частиц,

У»Ч

 

взаимодействия. Тогда, как

в этом легко мо­

а Н \ — оператор

жет убедиться читатель с помощью непосредственного вычис­

ления, мацубаровская матрица плотности 5(Р) удовлетворяет уравнению

(Р)/дР = — Ях (p/S (Р).

(7.12)

Граничное условие очевидно:

 

5(0 )= 1,

(7.13)

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

Нг (т) = ехр [ - (рЯ - Я0) г] ЯI ехр [(рЯ - Я0) т] (7.14)

/* ч

соответствует оператору Н{ в представлении взаимодействия. Отметим, что каноническое представление (7.14) отлично от

* В соответствии с общими правилами под ехр (—(ЗЯ) понимают опера­

тор, собственные функции которого совпадают с собственными функциями

гамильтониана Я, а собственные значения есть ехр (—Р£„).

3 Зак. 635

65


обычного гейзенберговского представления операторов через

шредингеровский оператор (в данном случае H i ) тем, что вместо

полного гамильтониана Н в «обкладки» входит лишь оператор

системы свободных частиц.

(электронный газ)

Рассмотрим однокомпонентную плазму

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

«простран-

ственно-временная» точка х = (г, т), тогда

оператор

H i в пред­

ставлении взаимодействия можно записать в виде

 

#! (т) = е2j tjj+ (х) и (х) i| (x) dx,

(7-15)

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

ф+ (х) = exp [— (\iN Нп)т] ф* (г) exp [(p/V — Н0) т];

ф (х) = ехр [— (pJV — Н0) т] ф (г) ехр [(рЛ/ — Н0) т].

Поскольку среднее от матрицы S(p) по ансамблю невзаимодей­ ствующих частиц

<ГS ^ — SPS (Р) ехР [НрЛГ — #о)1

SPP _

ZN

Sp exp [p (рЛ' — H0)}

Spf)„

l n

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

pQ = pQ0ln < S ( p ) > 0,

(7.16)

где

 

 

Й0 -

In Zjv0)

(7.17)

естьтермодинамический потенциал системысвободных

ча­

стиц.

Очень существенно, что введенное представление взаимодей­ ствия позволяет вычислять средние величины по ансамблю не­ взаимодействующих частиц. Таким образом, если бы удалось

выразить величину < 5 (р )> о в виде ряда по параметру взаи­ модействия е2(точнее, по отношению амплитуды рассеяния к среднему расстоянию между частицами), то, во всяком слу­ чае, в рамках теории возмущений можно было бы решить за­ дачу о вычислении термодинамического потенциала, если, ко­ нечно, при этом не появится расходимость в отдельных членах ряда теории возмущений. Как будет показано ниже, расходи­ мость действительно имеет место в случае системы кулоновских частиц.

66


Возникает вопрос, зачем нужно строить способ описания системы по теории возмущений, если уже давно существует хорошо известная термодинамическая теория возмущений (см., например, работу [2]). Дело в том, что вычисление каждой следующей поправки к термодинамическим функциям по стан­ дартной теории возмущений связано с громадными математиче­ скими трудностями. Построение же теории возмущений мацубаровским методом чрезвычайно просто. Введение диа­ граммного метода, аналогичного методу графиков Фейнмана в квантовой электродинамике, позволяет выделить наиболее су­ щественные члены ряда теории возмущений и провести выбо­ рочное суммирование ряда диаграмм, отвечающих разным по­ рядкам по параметру взаимодействия.

Как известно из теории

поля,

решение

уравнения

(7.12) с

граничным условием (7.13)

можно представить в виде

 

S(P) = Гехр

 

 

 

 

 

 

 

• f dxn X

 

 

 

 

 

 

 

 

 

О

 

X Т [Я4 (Tj),

ЯДт2), .

.

., ЯД т„)],

 

(7.18)

где

Т — оператор

упорядочения

по

 

«мнимому

времени»; т==

— it.

Оператор Вика Т действует на

произведение операторов

так,

что операторы,

относящиеся

к «более

раннему

моменту

времени» т, располагаются справа,

а

операторы,

относящиеся

к «большему» т, располагаются слева.

Так,

упорядоченное по т

«*■4 /N

произведение операторов ф+ и ф есть

ф + (та) ф (Т2)

Т {ф+ (т^фДг)) =

ф (т2) ф+ (tj)

при T j > T 2

^ 19^

при тх < т2.

Задача состоит в том, чтобы построить ряд теории возмущений

по е2 для величины < S(f})> 0. Для этого удобно ввести одноча­ стичную и двухчастичную функции Грина в представлении взаи­ модействия:

Gi(х4, л-2) = < Т ф+ (хх)ф(х2) S (Р)> 0 = < Гф+ (хх)ф(х2)> ;

| ^ ^

С2 (хь х2; Хд, х4) = <Тф + Ы Ф+ (х2) ф (х3) ф (х4) S (3) > 0,

]

где индекс 0 означает по-прежнему усреднение по ансамблю не­ взаимодействующих частиц.

з* 67


6

Представим | H\(x)dx в виде

о

 

j Нг (* ) dx = ~ J ф+ (a'j) \ jx (х2)и (хг — х2) ф (xs) ф (х4) X

 

X б (ati х3) б (х2 — х4) dx1dx2dx3dxi,

(7.21)

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

U(Л4 — Л'.2) -> hi (Xj — А'о).

(7.22)

Продифференцируем отношение Zx/ Z ^ = <5(|3) > 0 по парамет­

ру Я,:

#)=<:?РХе

 

 

 

 

d

^ H 1(x)dxyQ. - L .

(7.23)

 

 

 

dX

 

 

 

 

 

 

 

 

 

 

Разделим

обе части равенства (7.23) на < S (P )> 0.

При

этом

получается,

что

 

 

 

 

 

d in (ZN

 

= —

f dx1dx2dx3dx4u (x4 — x2) б (x4— x3) X

 

dk

 

2k

J

 

 

 

X 6 (x2 — x4) < Г ф+ (a'j)ф+ (x2) ф (x3) ф (a-4) S (P)> 0 =

 

1

c

2 (Xu

Аг’ *3’

X*>u ^ ~~ x^ 6 (Xl ~

6 ^

~ Xi) X

= ~2X

)

 

 

 

X dx4 dx« dx3 dx4.

 

 

(7-24)

Следовательно, можно установить простую связь между тер­ модинамическим потенциалом и двухчастичной функцией Грина G2:

ег

Q — Q0 = — С — С Go (*х, х2; х3, х4) и (х1— х>) б (х4 — х3) X

2 J е* J

 

о

(7.25)

X б (х2 —■х4) dx4 dx2 dx3 dx4.

В импульсном представлении, в котором обычно удобнее рабо­ тать,

 

ег

 

 

Аа “ 1г “ Q" =

f0f

- j’ (p-

"*> x

X U(Pi — p2)dpI dp2 dpa;

p4 = Pi + p2 — p3.

(7.26)

Отметим, что интегрирование по «координатам» и «импуль­ сам» записано условно, так как в технике температурных функ-

68


ций Грина координата

х=(х, х4) =

(х,

т)- Аналогично

импульс

р также четырехмерен,

т. е. р = (р ,

р4).

Поэтому, если

в даль­

нейшем будет записано интегрирование по х и по р, то это озна­ чает, что в обычном смысле следует интегрировать по вектор­ ным составляющим этих величин. Что же касается четвертой компоненты импульса, то ввиду конечного интервала изменения т(0, р) по р4 проводится не интегрирование, а суммирование по

дискретным значениям

 

 

 

 

р4 = (2п + 1)

л/Р,

п = 0,1,2, .

. .

(7.27)

Технику суммирования

пока

не будем

рассматривать здесь,

а остановимся на физическом смысле введенных формально од­ ночастичной и двухчастичной функций Грина. Одночастичная функция Грина G (1, Г) описывает распространение возмуще­ ния, при котором одна частица добавляется или удаляется из многочастичной равновесной системы. Например, когда оператор рождения, действуя первым, вызывает возмущение, до­ бавляя частицу в пространственно-временной точке (iv, т4-). Затем происходит распространение возмущения, которое про­ должается до «момента времени» п , когда частица удаляется в точке (гь t\), в результате чего возмущение снимается и систе­ ма возвращается в равновесное состояние. При Ti<Ti' первым

действует оператор ф. Распространение возмущения, которое

возникает теперь

из-за удаления частицы в точке (гь /i), про­

должается до «момента времени» т г ,

когда оно снимается в ре­

зультате добавления частицы в точке

iv . Аналогично, двухча­

стичная функция

Грина G2 (1, 2; 1',

2') описывает различные

возмущения, вызываемые удалением или добавлением двух ча­ стиц, производимыми в различной последовательности. Напри­

мер,

если Ti и т2 наступают «позже», чем тг, т2',то

G2 описы­

вает

добавление двух частиц,

за которым следует

удаление

двух частиц. Если же ti и т г

наступают «позже», чем т2 и т2->

то двухчастичная функция Грина описывает возмущение и воз­ вращение к равновесию в результате добавления и удаления одной частицы. Разумеется, столь непосредственный четкий фи­ зический смысл имеют так называемые временные функции Грина. Поэтому не случайно ставим слова «позже» и «рань­ ше» в кавычки, когда говорим о температурных функциях Грина, для которых мерилом «времени» является обратная температура |3 (точнее, ip).

Перейдем к графическому изображению отдельных членов ряда теории возмущений по константе взаимодействия е2. Вообще метод графиков можно представить себе просто: в виде рисунков изображают различные физические процессы, которые происходят с частицами. Распространение, например, кванта света или плазмона можно представить пунктирной линией, а распространение свободной частицы можно изо­

69