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

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

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

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

Добавлен: 15.10.2024

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

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

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

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

о

-

- о

и

о с г ю

 

 

о

-

- о

Рис. 5. Кольцевые диаграммы.

 

 

рис. 5. Нетрудно видеть, что члены

наинизшего порядка по е2,

т. е. члены наинизшего порядка по

амплитуде рассеяния / =

= е2р, суть кольцевые диаграммы (интегралы)

(3*0.

С точностью

до соответствующего комбинаторного

множителя,

определяе­

мого числом способов построения кольца, неприводимые коль­ цевые интегралы имеют вид:

Р*, о =

(1/2V) J" /12/23 • •

 

, , dr,dr, . . .dr

“+•’

)

(6.37)

Pi,о =

№ ( г),

 

/к+1,1 1 а

 

 

 

 

 

 

 

где f(r) задано выражением

(6.1)

с кулоновским

потенциалом

взаимодействия

 

 

 

 

 

 

 

и (г) — е2/У.

 

 

(6.38)

Как было показано во второй главе, f(r) представляет со­ бой не что иное, как бинарную функцию корреляции. Тогда, если плазма слабо неидеальна, т. е. амплитуда рассеяния ре2 мала по сравнению со средним расстоянием между частицами, бинарную функцию корреляции можно приближенно записать

в виде

/ (г) «

р« (г) + ± ри2 (г), M

l /кТ.

(6.39)

Тогда

 

 

 

р, 0 = у р2 j

и2 (г) dr + о (евР3) *

j UuWncMr,.

(6.40)

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

вие требования квазинейтральности плазмы.

В первом приближении можно учесть поправку к давлению идеального газа, считая амплитуду рассеяния очень малой, т. е. формально устремляя / = ре2 к нулюТогда для кольцевого груп­

60


пового интеграла рк>0, соответствующего диаграмме с (7с+1) вершинами, имеем выражение

Рк о = (—

^и ( ^

и (г2з)

•и (Гк+ i,i) с1ггdr2 .

. . drK+{.

 

 

 

 

 

 

 

 

(6.41)

Поскольку фурье-компонента

кулоновского потенциала взаи­

модействия имеет вид

 

и (q) =

4ле2/<72.

 

(6.42)

то можно написать

 

 

 

 

 

 

 

 

 

 

(~ Р)*+1

 

 

 

 

(6.43)

 

Рк ,0 ~~

(2я)з

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

5 (р) =

1-----рЗ» Г -Л 'Ш * .-

(6.44)

 

 

 

2 (2я)3

Kl

J 1+рр«(?)

 

Подставляя формулу

(6.42) в выражение

(6.44), получаем

 

 

5 (р) (5-1 = е3 (лрр)‘/г.

 

(6.45)

Из выражений (6.34) и (6.45)

 

получим

уравнение

состояния

вида

 

 

 

 

 

 

 

 

 

Р =

np~l — -i- (л$)'ие3п ‘и.

(6.46)

Второй член в правой части представляет собой дебаевскую по­ правку к давлению идеального газа.

Учет следующих членов разложения по малому параметру приводит к дальнейшим поправкам по плотности п к термоди­ намическим функциям системы.

§ 7. О ДИАГРАММНОМ МЕТОДЕ В ТЕРМОДИНАМИКЕ ПЛАЗМЫ. АНАЛОГИЯ С ДИАГРАММАМИ ФЕЙНМАНА

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

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

методов

кванто­

вой статистики многих частиц

(в частности,

и для

системы

многих кулоновских частиц—

плазмы). Этот

метод в

различ­

ных вариантах интенсивно разрабатывался как советскими, так и зарубежными теоретиками. Подмеченная в 1955 г. японским физиком Мацубара [5] аналогия между матрицей рассеяния S в квантовой электродинамике и матрицей плотности р в кван­ товой статистике сыграла большую роль в бурном развитии статистической физики бозе- и ферми-систем. Указанная анало­ гия позволила перенести основные элементы математического аппарата, хорошо разработанного в теории поля, для описания статистического ансамбля многих частиц. В частности, в кван­

61


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

диаграммную технику, являющуюся «двойником»

известного

в квантовой электродинамике графического метода

Фейнмана.

При квантовомеханическом исследовании систем, состоящих из очень большого числа одинаковых, произвольно взаимодейст­ вующих частиц, удобным способом описания системы является м е т о д в т о р и ч н о г о к в а н т о в а н 'и я. Смысл его заклю­ чается в том, что для описания свойств системы используется волновая функция, зависящая не от координат каждой частицы системы, а характеризуемая числами заполнения, т. е- числом частиц, находящихся в данном состоянии (N\, N2, ...), а именно волновой функцией фдг,, лг2••• (основные формулы см. в при­ ложении) .

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

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

информацию.

 

 

системы,

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

состоящей из большого

числа

одинаковых частиц *,

удобно

описывать

с помощью

полевых

г е й з е н б е р г о в с к и х

опе ­

р а т о р о в

р о ж д е н и я и у н и ч т о ж е н и я частиц-

Эти

операторы в отличие от соответствующих операторов в представ­

лении Шредингера зависят от времени. Оператор ф+(г, t), дей­ ствуя слева направо на какое-либо состояние системы, добав­ ляет одну частицу к этому состоянию в пространственно-вре-

/ Ч

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

через произведение операторов ф+ и ф **. Например, плотность частиц в точке (г, t)

« (г, t) = ф+ (г, ОФ(г, 0-

(7-1)

* Одинаковых опять-таки из соображений простоты.

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

62


Поскольку акт удаления и последующего немедленного восста­ новления частицы в точке (г, t) служит мерой плотности частиц в этой точке, оператор полного числа частиц можно предста­ вить в виде

N (t) ■-=f dnj>+ (г, /) -ф (г, t).

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

 

Л2 /•

л

(г>

+

Н (t) —

2^Г J dr^ + (г’ *)

+ -J-1 drdr'y+

(г, t) зр+ (г',

t)v{ \ г — г '|) i|) (r'f

t) ф (г, t). (7.2)

При этом первый член в правой части описывает одночастич­ ный оператор кинетической энергии системы, а второй — опера­ тор двухчастичного взаимодействия — представляет потенци­ альную энергию системы. Если система частиц находится во

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

/ч /ч

ниже, приобретает другой вид (равно как и операторы -ф и ф+).

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

онных соотношений следует, что ф2(г, t)=Q. Это и есть выра­ жение принципа Паули, запрещающего нахождение в одной и той же пространственно-временной точке двух идентичных фермионов.

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

нее значение любого оператора А, характеризующего систему, может быть получено усреднением по большому каноническому ансамблю:

 

X « ' I А | i> exp [— р (Ei—iiNi)]

 

< А >

i

 

(7.3)

 

1i. exp [— P (Ei — pNt)]

 

где | t > — нормированная на

единицу волновая

функция со­

стояния системы с энергией

и числом частиц А,-; р — хими­

ческий потенциал.

Суммирование проводится по

всем состоя-

63


ниям системы с любым возможным числом частиц. Более ком­ пактно среднее значение (7.3) может быть представлено в виде

S p [e x p { -p (tf-^ A /)} A ]

Splexp { - Р ( « - р Л 0 } ]

где символ Sp характеризует след оператора. Термодинамическое состояние системы определяется пара­

метрами р и р. Основному состоянию системы соответствует нулевая температура, т. е. р_1 = 0.

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

гих

кулоновских частиц. Если мы хотим

изучать

свободную

энергию системы F(р, V), то она выражается через

реальную

плотность частиц

в системе n = N/V. Тогда

уравнение

состоя­

ния

получается

в

явном виде

дифференцированием

свободной

энергии F(р, V)

по объему при постоянном числе частиц N. Если

же

нужно

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

потенциал

 

системы

Q(p,

V, р),

т. е.

рассмотреть

систему с

переменным

числом

частиц, то Q выражается не через реальную плотность п, а че­

рез

величину

 

 

 

 

 

 

Тогда плотность п необходимо также выразить Е(р), и для по­

лучения уравнения состояния в явном виде

нужно исключить

из полученного уравнения

и из уравнения

(7.4) параметр р.

Отметим, что параметр £,

который в случае

вириального раз­

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

Описанные выше подходы к рассмотрению термодинамиче­ ских функций соответствуют двум альтернативным описаниям системы: каноническому ансамблю Гиббса и большому кано­ ническому ансамблю. Поэтому и статистическая сумма по со­ стояниям системы записывается по-разному в этих двух слу­ чаях. А именно:

f

=

_ p - 4 n Z G,

(7.5)

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

 

 

 

ZG=

2 e x p ( - p £ J

(7.6)

 

 

К

 

выражает суммирование

по состояниям системы к

с фиксиро­

ванным числом частиц N.

 

 

 

Во втором случае

 

 

 

П =

- р - Ч п ZN,

(7.7)

64