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

Следовательно,

Если считать, что одна из частиц, например а, находится в на­ чале координат, то выражение (4.30) совпадает с первыми дву­ мя членами разложения (4.20).

Связь между функциями корреляции различной кратности можно описать системой уравнений, отличной от системы (4.23). Простой вывод этой «зацепляющейся» системы уравнений мож­ но найти, например, в работе [7]. Именно в таком виде эта система уравнений была впервые получена Н. Н. Боголюбовым [1]. Не будем выписывать последовательную систему уравнений Боголюбова («боголюбовские цепочки»), а приведем лишь урав­

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

В том же приближении можно считать, что даже пары частиц не находятся настолько близко друг к другу, чтобы Каъ было существенно отлично от единицы. Тогда, вводя малые величины

в правой части член, содержащий лишь f ac. Остальные же чле­ ны обращаются в нуль вследствие изотропии газа. В первом

31

Если взять дивергенцию от обеих частей этого уравнения, помня, что

uab = ZaZ„e2/r, г = | г„ — га |,

которое является уравнением типа Дебая — Хюккеля для само­ согласованного поля, уже рассмотренным для однокомпонент­ ной плазмы в § 1. Подставив решение для f(г) в формулы (4.37) и (4.33), опять получим выражение (4.30) для бинарной корреляционной функции в дебаевском приближении.

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

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

§ 5. МЕТОД КОЛЛЕКТИВНЫХ КООРДИНАТ

Одним из интересных подходов к описанию поведения систе­ мы кулоновских частиц явился развитый Пайнсом и Бомом [9] в 50-е годы метод коллективных координат, который сыграл большую роль в понимании многих явлений, имеющих место в низкотемпературной плазме, а также в квантовом электронном

32

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

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

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

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

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

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

Прежде всего рассмотрим вопрос о движении совокупности частиц в методе коллективных координат и выясним область применимости теории при таком подходе. Для простоты, как и раньше, объем системы положим равным единице. Тогда в слу­ чае точечных частиц их плотность

где гг— координата i-го электрона.

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

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

При этом символ S' означает, что член с к = 0 следует опустить.

Выразив v ,= п через уравнение движения для t-й частицы и подставив получившееся выражение

к. /

34

в формулы (5.4), получим уравнение, описывающее изменение «к во времени, т. е.

Разобьем вторую сумму на две части, соответствующие об­ ластям суммирования с к = к' и к^=к'. При этом сумма с кт^к' содержит фазовые множители exp[i(k'— к)г*], которые зависят от положения данной частицы. При усреднении эти фа­ зовые множители должны быть малы, так как в отсутствие сильного взаимодействия имеется большое число частиц, распре­ деление которых в пространстве близко к хаотическому. Если считать, что результат такого усреднения равен нулю, то члены с кф к ' следует опустить, что соответствует приближению хао­ тических фаз, или RPA (Rand Phase Approximation). Тогда вместо уравнения (5.6) имеем

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

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

Из неравенства (5.10) следует более наглядное выражение

Лк —1~Ыр Лк —0

2 35