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

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

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

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

Добавлен: 15.10.2024

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

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

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

Разложим 6Ni в ряд Фурье

с периодическими

граничными

условиями в кубе единичного объема:

 

6.V,- = V

лкеехр (ikr,-),

(5.33)

k+ 0

 

 

причем «1 =п-к. Поскольку величина лк комплексна, ее удобно

записать в виде

 

лк = rkexp(icpk),

(5.34)

где гк = г_к', фк = —ф—к* Отметим, что независимой здесь являет* ся только половина величины лкНетрудно видеть, что

bNfiNj

 

 

лк |2 =

2яе2^ - ^ - Г к .

(5.35)

2 Z j I Ч — ri

I

 

 

 

 

it

 

t+i

 

 

 

 

 

Тогда

 

 

 

1

 

 

W ~

 

exp

2ne2Pn

2

(5.36)

 

 

Гк

где л — средняя плотность электронов.

Найдем теперь элемент объема в пространстве Гк и фк, рас­ сматривая N{ как непрерывную переменную и предполагая, что в выражении (5.32) 8Ni = l. Тогда это выражение можно умно­ жить на элемент объема в ^-пространстве (dNь dN2, ..., dNit...), который также равен единице. Элемент объема в пространстве Гк >фк получим с помощью якобиана J при преобразовании от iVj-пространства:

X [^Фк_с/фк2 . . .^Фк,1-

(5.37)

Вероятность WJ есть произведение функций, каждая из кото­ рых зависит от единственного значения к. Физически это озна­ чает, что в рассмотренном приближении статистические флук­

туации,

связанные с различными значениями

волнового

числа

к,

независимы.

учитывая,

что

г-к

Вычислим теперь среднее значение <Гк2> ,

~ г к

и величину г 2 в выражении (5.36)

можно

заменить

на

2 г 2.

Тогда суммирование необходимо проводить

по

поло­

вине возможных значений к. В результате получим

 

< r 2 ^ $Wrk [ •

d (rk) • • • ]( • •

-^Фк • • •) =

___ П»___

] > [ . .

. d ( r l ) . . . ] ( . .

. d < t k . . . )

М +

40


или

<Гк> =

п

(5.38)

lD2

 

Если /г2» х 2 = /д2, то <г£ >-*-п

и соответствует хаотическому

распределению частиц в системе. В обратном предельном слу­ чае (62<Cx2), получим

<г1> = пкЧ1.

(5.39)

Забавно, что этот результат можно себе представить с по­ мощью следующей аналогии. Пусть по достаточно большой пло­ щади в разных направлениях движутся автомобили, причем их движение довольно хаотично. Конечно, водители стараются из­ бежать столкновений. Милиционеры также стараются коррели­ ровать движение. Но увы! Аварии случаются! Пусть число этих аварий и есть флуктуация плотности автомобилей на площади. Рассмотрим два случая. Пусть погода настолько плоха (туман), что водитель может видеть соседнюю машину лишь на малом расстоянии %~а, где а — размер автомобиля. Тогда, поскольку

ft-СId, движение

автомобилей

хаотично

(или

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

независимо) и число аварий пропорционально

их

плотности п.

В ясную погоду (ft» /D) водители хорошо видят друг друга,

видят они и милиционеров. Число аварий

падает

в результате

«осмысленного»

коррелированного

движения

в

отношении

(1о/к)2, где / д — характеризует

некоторую

«площадь взаимо­

действия». Этот результат и описывается

формулой (5.39), из

которой следует

важный вывод, что

на больших

расстояниях

к у л о н о в с к и е с и л ы п р и в о д я т к с и л ь н о м у с н и ­

ж е н и ю ф л у к т у а ц и й .

Поскольку вклады отдельных компонент Фурье, как это было отмечено выше, независимы, можно написать среднюю кулонов­

скую энергию для

данного k:

 

< ик> =

< 2netnl/k2> = (1/2) [P(l + k 4 l ) ] - \

(5.40)

При &<Сх <ык> = 1/(2р). Если вспомнить, что в случае длинных волн компоненты «к ведут себя почти коллективно и колеблются гармонически, то предыдущий результат легко понять физически на основе закона о равнораспределении энер­ гии по степеням свободы, согласно которому средняя потен­ циальная энергия гармонического осциллятора равна 1/(2р).

При /г» х

< « к > = 2гоге2/&2.

(5.41)

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

41


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

~ р - 1= и \

Подчеркнем в заключение, что приближение, в котором Ni считается непрерывной величиной, справедливо, строго говоря, лишь для малых значений волнового числа k:

где г0— среднее расстояние

между частицами,

определяемое

выражением (1.2). Тем не

менее для больших

k

это прибли­

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

значение < г *>,

а именно п.

Это означает, что выражение (5.38) практически

 

годится

для

любых значений k.

 

физически к методу Бома и

Несколько иной, но близкий

Пайнса подход к описанию роли флуктуаций плотности в тер­

модинамике был предложен Н. Н. Боголюбовым

[1] в

50-е

годы (а также В. М. Галицким и А. Б. Мигдалом). Речь идет

о так называемом методе «лишних переменных»

(коллективных

переменных).

Рассмотрим однокомпонентную систему кулоновских частиц,

находящуюся в термодинамическом равновесии.

Запишем кон­

фигурационный интеграл для системы N частиц:

 

 

ы(| г* — Г] |)1сМг2 . . .drN,

(5.42)

где ы|г*—iy |) — потенциальная энергия взаимодействия

куло­

новских частиц.

частиц гу

как

Введем наряду с координатами отдельных

бы лишние переменные п к [см. выражение (5.1)], в качестве

которых выступают компоненты Фурье от плотности частиц п(г). Разложим и(т) в ряд Фурье*:

и (г) = 2>(k) exp(ikry).

к

 

Энергию взаимодействия выразим через

переменные п^\

/1*7 а

 

 

(5.43)

Из выражения (5.42) следует, что

(5.44)

v (к) = 4ле2//г2.

Для простоты, как и раньше, положим объем системы равным единице.


Если это так, то второй член в правой части уравнения (5.43). бесконечен. Однако не нужно забывать о наличии компенсирую­ щего равномерно распределенного заряда, который для k = 0 приводит к тому, что этот член следует опустить. Физически это

означает выполнение требования

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

Выражение для энергии

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

(5.43)

допускает

факторизацию конфигурационного

интеграла

Zn.

Перейдем

в этом

интеграле от переменных

гь г2, ..., rN к

переменным

гь г2, ...,

rN, nk], «k2. ..., п kp

т. е. будем как бы приписывать си­

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

ZN = exp

n2ftv (0)

. Р

 

X

 

 

 

 

 

 

X jexpj—| ^ J n v ( k ) n kn_k| D ( . . . nk . . .) . . . (dnkdnk)

. . .,

где

 

 

 

 

 

(5.45)

 

 

 

 

 

 

D (. . . nk . . . ) = j П

6 [”k — y f

cos (kr j) ] X

 

X 6 flk -y=r ^

Sin (kr/)

drldrZ

• • drJV>

(5.46)

nk = n? + inf;

 

n£ =

ney\

n* =

— щ к.

(5.47)

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

ными. Функция D(... nk...) играет роль якобиана для перехода

от переменных rj

к n k и может быть

представлена

в виде

D ( . . . nk . . .) =

) exp

in 2

«к

X

 

 

 

 

к

 

 

X jexp

— ^

®кехР( 'кг) ^rJ

(d(Sl)>

(5.48)

где

 

 

 

 

 

 

(do) = П (dcokdcokj;

о)к =

сок — цок :

 

 

сок = со_к;

со_к = — 0

 

 

43


При этом использовано известное интегральное представление для 6-функции:

оо

6 (х) = j exp (2лт.ш) day.

Можно явно проинтегрировать выражение (5.48) по dr, если

разложить экспоненту, стоящую под интегралом в ряд по Ограничившись тремя членами разложения

ехр

 

 

 

 

 

 

 

 

(5.49)

получим

 

 

 

 

 

D (. . .пк . . .) =

3! Уп X

X

д3

 

А>(. . . nk . . .),

dnv дпк дпк

где

ki-bk2-j-k3—О *'1 2

 

 

 

D0( . . . пк . . .) — J' ехр

2

сок«к — я2 2

©к(о_к\ (dto). (5.50)

 

I

k

к

J

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

D0( . . , nk . . .) = ехр

Y ^ (ЯкЯ-к + In л)

(5.51)

Тогда

ZN = ехр j— -^ - v (0 ) + -|- ^jTjnv

/X j* ехр | ---- —l ^

[nPv (k) + 1] rtk«-:< + In я |

. . . (dnk)'

(5.52)

или

 

 

 

 

 

ZN = exp

v (0) 2Г nv (k) pj IT J exp 1— [(nv (к) P +

 

+ l) ((«к)2 +

(«к)2) + In я]) . . . dnidnl.

(5.53)

Вычислив интеграл no

d

n получим

 

 

Z„ = exp { -

(0) +

j ' «* « (■} IT

(5'54)

 

 

 

к

 

 

44