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

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

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

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

Добавлен: 15.10.2024

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

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

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

где а = ре2игЛ3/2<С1.

В разложении

(13.23)

по а сохраним чле»

ны с отрицательными степенями а

и

члены, пропорциональные

In а. В результате интегрирования по К получим

 

- PAQ = (2/3) (яр3)‘/! | 3/а -

(1/3) (е2Р)312 In §e*xD). (13.24)

Рассмотрим теперь многокомпонентную систему из атомов,

электронов и ионов. При этом

 

 

 

 

AQ = — j ~ ~ j

dptdp2dp3 и (р3 — рх) ^

g ip 2 (Pi, р2; Рз, р^,

о

 

 

г/

 

(13.25)

 

 

 

 

 

где P4= P i+ P 2—Рз;

gij — константа

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

i-й и /'-й

компонент, а функции Грина

определяются из

уравнений

типа (12.1). Формулу (13.25) дополняют условия равновесия и квазинейтральности плазмы:

+ ^ < 7 ^ = 0. <13-26>

где pi — заряд частиц i-й компоненты.

Рассмотрим случай, когда в плазме присутствуют электроны е, ионы i с зарядом Z и (Z—1) раз ионизованные атомы. По­ следние представляют собой одноэлектронную связанную си­ стему, и в дальнейшем будем условно называть их просто ато­ мами а. При этом в плазме имеются следующие виды взаимо­ действий: ее, ii, ie, аа, еа, ia. Взаимодействие с потенциалом отталкивания учитывается квазиклассически совершенно так же, как это делалось при выводе формулы (13.25). В результа­ те получим, например, для взаимодействия И:

— рдп,

т 2

+

л{33 I_________ ^_________

 

Ke+&ti + ( z - 1)2- Ы

 

+

(№Z*f g?In фе°2°хГ\

(13.27)

*a = W l h + Z % + ( Z - l Ш 'Ч

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

= - w Zexp 10 ^ +1*'» (-wP х

x .ff j* dpdp'u (р — р') exp (— 0££е) (р') ф£е (р) (13.28)

о

149



к квазиклассическому пределу, начиная с уровня £n„>i легко получить сходящееся выражение

 

_

Р (AQJKB=

 

оf ^

оIuexp (е~а71) X

 

X /

x ( l - e r f y ^ i a l ^ - J - j + ^ e x p H a O X

 

exp (— at)

 

оо

 

 

Л

 

 

 

 

 

 

 

it

+ j

[exp (е~а70 -

1 ] tdt\,

(13.29)

 

 

 

где

 

 

 

 

 

 

 

 

 

*». =

26nl/V; б =

 

«

1;

а =

«

1.

Разложение по а с указанной выше точностью приводит к сле­ дующему выражению, не зависящему от По:

- р (АЙ = - y

VWe%teZ* Це + Z%i + ( Z - l f

1аГ Ч‘~

-

^ ( p e 2Z)3 U eIn (PеЧх)~\

(13.30)

Аналогично вычисляется вклад от квазиклассической части взаимодействия еа. Вклад от основного состояния атома Е0, ко­ торое в первом приближении можно считать кулоновским, уже учтен в й 0- Вклад в уравнение состояния от сдвига основного уровня можно получить, пользуясь построенной нами четырех­ мерной теорией возмущений [3]. При этом рассмотренное выше адиабатическое приближение уже не является правильным, а взаимодействие не является потенциальным, поскольку и зави­ сит ОТ <74-

Построение термодинамической теории возмущений

Нетрудно видеть, что Ай можно выразить через функцию Г, определенную формулой (13.14). Эту функцию иногда называют

о б о б щ е н н о й в е р ш и н о й , или

ч е т ы р е х п о л ю с н и к о м ,

или э ф ф е к т и в н ы м в з а и м о д е й с т в и е м :

dpldg’ )13'1^

AQ= (2л)®к оf^’Р^l(PГIФ

Р’

ег

 

 

р

 

 

причем Г удовлетворяет уравнению

 

 

Г (р, p',g) = (2я)3 РU (р - р’) + ^

J U (р -

Р') Ф (Р, g) X

Х Г (p,pf,g)dp'.

(13.32)

150


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

лоновского

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

м(р) и Ди(р, р4) — малого возму-

о

а- о 0 -

+ 0

6

 

 

 

+ *»*+ © +

+

яс

+ - - -

Рис. 13. Диаграммы лестничного типа.

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

U = и + Ли.

(13.33)

Эффективное взаимодействие Г также представим в виде

Г(Р, р', g) = Г0(р, р', g).+ Г,(р, р \ g),

(13.34)

причем Г0 изображается «лестницей» с кулоновским потенциа­ лом и (рис. 13). Подставляя функцию (13.34) в уравнение (13.32) и пренебрегая членами второго порядка малости, полу­ чаем уравнение для Ti

Г1 (Р. Р', g) = (- ^ - р J и (Р— Р>) Ф (Ри g) Гх (ръ р \ g) dPl +

+

j Ли (p —pi) <p (pi, g) Г0 (plt p', g) dp,.

(13.35)

В первом приближении учтем второй член в правой части этого уравнения, т. е.

Г*1' =

2 j d4Au (Ч. Яд Ф — <7. g) Г0 Я, р \ g)- (13.36)

'р я*

Легко видеть, что Г{0) не зависит от р\. Решая уравнение

(13.25) методом итераций, убеждаемся, что Г, также не зави­ сит от р'. Выражение (13.36) для Г,(1) с Ли в виде петли можно

151


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

Рис. 14. К построению четырехмерной термодинамической тео­ рии возмущений.

Аналогичный график ж получается для собственно энергетиче­ ской части.

Пусть

X(Р. Р', g) = Гх (р, р', g) Ф (р, g).

(13.37)

Тогда уравнение для %имеет вид уравнения Шредингера с пра­ вой частью:

© (Р, Р. g ) = 72Э Т "

2

Ч>( Р ’ Я ) j d4Au (4.fc) Ф(Р — Я ’

Я ) X

v Р

ЯАуРа

 

 

X

Г0(р — q, р', g),

(13.38)

так что решение для %выражается через ф-функции относитель­ ного движения двух частиц. Записывая в качестве Аи цепочку вместо петли, выполняя суммирование по четвертым компонен­

152

там импульсов,

получаем с помощью выражений (13.31),

(13.37)

и (13.38) при q—>-0

 

 

 

 

 

*

 

о" f

f dpdpxdgexp (р (ц, +

ц, — в*)) X

 

 

(2я)«р

J е*

J

 

 

 

 

 

 

 

О

 

 

 

 

X 5

]

 

(Р)

(Р) (Е'к -

ер') < (Pi) %

(Pi) X

 

х

п , к

 

 

 

(£„ — eP) exP (— P£n)

(13.39)

( 4

- s p) e x p ( — P £ k)

 

 

 

 

 

 

 

 

 

(£k - £p,)2

(£; - v ) s

 

С точностью до

 

высших

степеней

отношения

£ n/p-Cl наиболь­

ший вклад в AQ дает член с n = k. Выделяя вклад основного со­

стояния атома,

имеем

 

 

 

z j

 

| 1>„<Р) I * -

 

- е х р ( - Р Е „ ) ф f

 

 

J

£0-- £Р1

В результате простого вычисления получим

 

AQ^ 5а0>ша/22Р,

(13.40)

где па— плотность «атомов».

 

к Ай от собственно

Аналогично

вычисляется поправка

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

го вычисления получается выражение типа

(13.40), где

вместо

коэффициента 5/2Z стоит —5(Z2+1)/4Z2.

Вычисление

графи­

ков типа е й ж (см.

рис. 14)

приводит

к

дебаевскому

члену,

обусловленному изменением энергии

иона

с зарядом

(Z—1)

при его движении

как

целого.

Этот

член

уже учтен

ранее.

С учетом поправки

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

для

термодинамическо­

го потенциала рассматриваемой системы с точностью до чле­

нов порядка | 2

имеет вид

 

- S

5'— S T

■+ f

1 “

 

~

т ( Ч ± У а°у11а-

(13-41)

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

153