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

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

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

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

Добавлен: 15.10.2024

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

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

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

После интегрирования в выражении (16.22) по заряду по­ лучим

Й = Й0 + АЙD+ — [

n%s — 2л —

п%3+

 

Р

L 2,/а

 

№ее

 

 

ч,

Z_

1 + In 2 \

 

06.26)

+ <X%KD

2

2*'.

)

'

4 ( 2 + 1 )

 

где a = e2j3xB. Этот результат

отличается

численными коэффи­

циентами от соответствующего выражения, полученного в ра­ боте [2], в которой допущены некоторые неточности. Эти не­ точности исправлены в приведенной формуле для й [9]. Отме­

тим, что член, пропорциональный а2, не

может быть

получен

с помощью введенных в этом параграфе

квантовых

корреля­

ционных функций, поскольку последние вычислены лишь с точ­ ностью до а. Выражение для этой поправки АЙ, однако, мож­

но получить из следующих

полуклассических соображений.

Запишем классическое выражение для бинарной

функции

корреляции в виде суммы двух членов

 

К к л ( r ) = 1 -

ехР( - И л г) + К л (г),

(16.27)

где второй член порядка а2. Тогда вклад в термодинамический потенциал, обусловленный К'кл (г), равен

е*

 

Ай = -у f d(? j -у- [ K i i + К е е - Wei].

(16.28)

о

 

Выражение K lin можно определить, воспользовавшись тем фак­

том, что на малых расстояниях (r<^lD) взаимодействие двух рассматриваемых частиц практически не искажается присутст­

вием остальных зарядов (т. е. роль экранирования мала)

и обя­

зано иметь больцмановский вид:

 

К кл (г) = exp (— etfSlr).

(16.29)

Если в экспоненте дописать дебаевский экранирующий множи­ тель, то выражения (16.29) и (16.27) будут совпадать с точно­ стью до членов первого порядка по а. Тогда с точностью до членов второго порядка по а можно записать

Ккл (г) = ехр ^---- [1 +

а2ф (г) +

О (а3)],

(16.30)

причем функция ф(г) при

может

иметь

особенности сла­

бее особенности типа г-2, ибо лишь в этом случае выражение

(16.30) стремится к распределению

Больцмана (16.29)

при

r-v0. Сравнивая выражения (16.30)

и (16.27), найдем

 

' ехр (—2хд г)

(1 6 .3 1 )

 

+ ф ( г )

2 {гу,^ 2

 

192


Подставляя это выражение соответственно для взаимодейст­ вий ее, ei, и в формулу (16.28), можно видеть, что первый член, пропорциональный г~2, приводит к логарифмически расходя­ щимся при г-*-0 интегралам, тогда как члены с функциями Ф(г), не содержащими особенностей типа г~2, конечны. Огра­ ничиваясь логарифмической точностью, последние можно не учитывать, и тогда получим

Нижний предел интегрирования следует выбирать по-разному в зависимости от типа взаимодействия. Так, ион-ионное взаи­

модействие (классическое) должно быть

ограничено на

рас­

стоянии

минимального

сближения г “ин ~ Z 2e2|3, в то

время

как

электрон-ионное [как

следует

из выражения (16.11)] — на

электронной длине волны ft, а электрон-электронное — на

при­

веденной длине волны %ее=%\'

2. Учитывая эти замечания,

вы­

ражение

(16.32) с логарифмической точностью можно записать

в виде

 

 

 

 

 

 

 

2л(е2Р)2 (z 4 l n - ^ - f

l n - ^ — 2Z2l n - ^ ) .

(16.33)

Интегрируя по заряду в соответствии с формулой

(16.28),

по­

лучаем

 

 

 

 

 

 

^

= 12(/а2;-1)Р ( Z2 1ПW

+ 1П ^

- 222 1П П Г ) ’

 

которое совпадает с результатом работы [2].

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

7

Зак. 635

193


такой плазмы с помощью теории возмущений представляется неоправданным. Необходимо преобразовать исходный гамиль­ тониан взаимодействующих электронов и ионов так, чтобы но­ вый гамильтониан в нулевом приближении соответствовал трехкомпопентной системе. Затем уже по отношению к преоб­ разованному гамильтониану можно применить теорию возму­ щений. В работе [3] выполнено каноническое преобразование гамильтониана, отвечающее этой цели. Метод позволяет после­ довательно учесть взаимодействие между всеми компонентами частично ионизованной плазмы и избежать расходимости ста­ тистической суммы атомов.

 

СПИСОК ЛИТЕРАТУРЫ

 

 

 

 

1.

Бонч-Бруевич М. А., Глас ко В.

Б.

«Оптика и спектроскопия». 1963,

т.

14,

2.

с. 495.

«Ж.

экспернм. и теор. физ.», 1959,

т.

36,

Веденов А. А., Ларкин А. И.

 

с. 1133.

 

 

 

 

3.Воробьев В. С., Хомкин А. Л. «Теор. и матем. физ.», 1971, т. VIII. с. 109.

4.Красников Ю. Г. «Ж. экспернм. и теор. физ.», 1967, т. 53, с. 2223.

5.Ландау Л. Д., Лифшиц Е. М. Статистическая физика. М., «Наука», 1964.

6.

Ландау

Л.

Д., Лифшиц Е. М. Квантовая механика. М., Физматгиз, 1963.

7.

Ларкин

А.

И. «Ж. экспернм. и теор. физ.», 1960, т. 38, с. 1896.

1957.

8.

Смирнов В.

И. Курс высшей математики. Т. IV. М.,

Гостехиздат,

9.

Трубников

Б. А., Елесин В. Ф. «Ж. экспернм. и теор.

физ.», 1964,

т. 47,

с. 1279.

10.Carley D. Phys. Rev., 1963, v. 131, p. 1406.

11.Kudrin L. P. Proc. of the VI Intern. Conf. on Ioniz. Phenomena in Gases,

SERMA, Paris, 1963, v. Ill, p. 431.

12. Montroll E., Ward J. Phys. Fluids, 1958, v. 1, p. 55.


Г л а в а в о с ь м а я

ВВЕДЕНИЕ В ТЕРМОДИНАМИКУ

ПЛОТНОЙ ПЛАЗМЫ

§ 17. УМЕРЕННО ПЛОТНАЯ КЛАССИЧЕСКАЯ ПЛАЗМА

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

■Пил = <те23/г0«

1

(17.1)

г] = е ^ у - о «

1 •

(17.2)

Последнее условие означает, что в дебаевской сфере находится мало частиц (в случае точного равенства единице одна части­ ца). Если плазма квантовая, то система является умеренно плотной, когда параметр rpm — 1.

Если г)<С1, термодинамические функции можно вычислить с заданной точностью. В другом предельном случае г|3>1 за­ дача тоже решается «точно», так как возможно разложение термодинамических величин по обратной величине этих пара­ метров. Это решение, которое будет рассмотрено в десятой гла­ ве, соответствует очень плотной системе частиц. Для системы фермионов, как было показано в первой главе, случай сверхвы­ соких давлений описывается выражениями, соответствующими идеальному газу, поскольку для rinn^l потенциальная энергия взаимодействия на частицу мала по сравнению с кинетической, т. е. взаимодействие является слабым. Особенно труден для теоретического рассмотрения случай, когда сравнимы потенци­ альная и кинетическая энергия, приходящиеся на одну частицу. При этом говорят, что задача не содержит малого пара­ метра. Выход из положения состоит в том, чтобы продвинуться

в область г|~ 1 как

со стороны малых плотностей, так и со

стороны больших, а

затем построить правдоподобную интер­

поляцию.

Рассмотрим подробно весьма интересный подход в термоди­ намике умеренно плотной плазмы, показывающей, как далеко можно продвинуться при анализе классической плазмы со сто­ роны малого параметра. Этот подход построен на диаграмм­ ном методе в классической статистике [7], который является обобщением изложенной в третьей главе графической техники Майера.

195

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

заряды одного

знака вносят определяющий вклад

в

к о р р е ­

л я ц и о н н у ю

э н е р г и ю п л а з м ы . Например,

в

астрофи­

зике существуют условия, когда в водородной плазме электро­ ны вырождены, причем энергия Ферми велика по сравнению с

потенциалом ионизации атома водорода.

Тогда

рассмотрение

однородного электронного фона

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

удовлетвори­

тельным приближением, так как

протоны

не

вырождены и

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

талле.

дебаевский радиус

1в= (4л^2пр) -1/2 едини­

Удобно считать

цей длины и ввести следующие безразмерные параметры:

R = г ф о , е = q ^ ! l D = (4л)4

( q b i u

f, )4 \

(17.3)

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

между

частицами

i

и /.

Таким

образом,

отношение кулоновской потенциальной

 

энергии

к

кинетиче­

ской порядка е для

зарядов

на расстоянии

и порядка е 2/3

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

ы(1, 2) = e/R.

Итак, если рассматривается система N одинаковых куло­ новских частиц в объеме V при температуре р-1, удобно ввести так называемый п о т е н ц и а л с р е д н и х с ил ш( 1,2) сле­ дующим образом:

ехр[— w (1,2)] = ZN'V2j . . . Jd3, d4, . . ., clNX

 

 

X

exp

S

/)

 

 

 

где

Z.v — конфигурационный

интеграл.

Это

выражение

совпа­

дает

с определением

бинарной

функции

корреляции,

если

u(i,

/)= и (г,;)Р , где

и(гц) — потенциал

парного

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

вия. Функцию до( 1,2)

можно представить графически (так на­

зываемые кустовые диаграммы) с помощью

введения вспомо-

1ательной функции

Г (1,2),

удовлетворяющей

интегральному

уравнению

 

 

 

 

 

 

 

 

 

Г(1,

2) =

/(1, 2) -}- /I ( е!3/ (1,

3) Г (2, 3),

(17.4)

где п — плотность

частиц;

/= ехр [—и (1,2)]—1 — введенная в

третьей главе функция Майера;

Г — скалярная

функция

един­

ственной векторной переменной ггд Тогда

 

 

 

 

 

в>(1, 2) =

м (1, 2) -1-/(1, 2) — Г (1,

2) — X (1, 2),

(17.5)

196


где Х(1,2) выражается через бесконечную сумму многократных интегралов, а каждый многократный интеграл соответствует 12 неприводимым диаграммам специального вида, подробно опи­

санным в работе [7]. Подынтегральное выражение каждого из

этих интегралов

содержит

некоторую

комбинацию

функций

Г и /.

 

 

функция

f

стремится

к —и (1,2)

=

Для фиксированного R

= e/R при в->0.

В этом

пределе функция

Г (1,2)

стремится

к

корреляционной

функции

в

приближении

Дебая — Хюккеля:

Г0 (1, 2) =

- (q2$/r12) exp ( - r M

=

-

(г/R) exp ( -

R). (17.6)

Если в достаточно мало, то кустовые диаграммы сводятся к известным уже нам кольцевым диаграммам в классической ста­ тистике и представляют последовательность членов по степе­ ням в. Порядок диаграммы, содержащей I линий и к точек, есть г’~к, так как каждая линия вносит множитель, пропорцио­ нальный в, а каждое интегрирование дает множитель порядка

l2Dn ~ £ ~ \ Не

будем обсуждать здесь вопрос быстроты

сходи­

мости членов

высших порядков в этом разложении.

Однако

ниже обсудим

вопрос о том, почему в области в<С^<С1 нет

серьезных трудностей и можно разложить Г в ряд по степеням е (первый член разложения Г0), но для /?^в этот подход не­ правомерен и необходимо модифицировать теорию. Ниже будет показано также, что Г сравнительно просто вычислить с точ­ ностью до второго порядка по в, включая область R ^ i в. Ока­ зывается возможным оценить также Г при в=1, привлекая численный счет.

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

 

 

Г (/с) = / (/с) [ 1 — л/

к)]-1.

 

(17.7)

Для

больших R сферически

симметричная

функция

f(R) =

= ехр[—«(Я)]—1 стремится к

—u(R).

 

Ее

фурье-компоненту

запишем в виде

 

 

 

 

 

 

 

 

/ (к) — — (4ле//с2) Ф (т]),

и =

(игк)'/’.

(17.8)

Подставляя Цк) в выражение

(17.7) и переходя к обратному

преобразованию Фурье, получаем

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

Г (/?) = -

2_

Ф(л)

sin (kR).

(17.9)

 

л

 

 

+

Ф (Г))

 

 

 

 

О

 

 

 

 

 

 

Эту

функцию можно протабулировать

для любых значений е

и R.

Такое вычисление было проделано в работе [6]. Результат

197