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

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

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

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

Добавлен: 15.10.2024

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

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

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

где пе, щ, па — плотности электронов, ионов и «атомов» соот­ ветственно. Для графиков типа 2 получаются аналогичные вы­ ражения, нужно лишь вместо множителя Z подставить фактор

( Z 2 + 1 ) .

Волновые функции фп(р) для водородоподобных систем мож­ но записать в виде

Фп, /, m (Р) = F *l (Р)

 

(0> ф);

 

Л ,.д р )=

2

(п— / — 1)!

■J'* п222(г+1) /! X

 

[-J

(П + /)!

(12.24)

 

 

 

 

X

П V

Г1+1 ( П2Р2 — 1

 

(nV +

1)'+2 п~г_1 \

п2р2 + 1

 

где CVN (,v) — функция Гегенбауэра (см., например, работу [2]).

Подставляя выражения (12.24) в формулу (12.23), получаем для s-состояний:

А Д п , О, О

5 (

Z — 1

у

a0x.D

при n =

1;

\

z

)

~ Г

 

(12.25)

3

Z — 1

\ 2

a0x D

 

при П

2, 3, 4

J

Z

)

~~J~П

 

 

где а0 — боровский п= 2 получим

 

> Гп to

о о

=

14

(

Z — 1 'V

3

\

 

Z

-/

радиус; п — главное квантовое число. Для

— 6

Z

1

^2

а0Хд

АДгд ,о ==

 

Z

)

 

 

 

¥ f l

;

А£г,1,1

= АДг, 1,-1 = АДг, 1,о- (12.26)

р

 

 

 

 

 

Для графиков типа 3 получается следующий вклад в энергети­ ческий сдвиг s-состояний:

о

9

(Z — I)2

----e2ai

_____ — [Z% +

2

0

z*

+ (Z— 1)3/га —

 

при п =

АДп,о,о — ■

 

 

(12.27)

е2а02 (2-^ -1— [23п/ +

+ (Z— 1)3п„ —

 

при п = 2 ,3 ,...

Расходящимися при малых q являются также графики типа 5 и 6, где заштрихованным квадратом обозначен «эффективный потенциал»

Г (р, р', g) — 2,

(2л)3 Р (Д' — еР') S

(Дк -

Sp) Фк (Р) Фк (Р') (12.28)

 

Д ' - Д к

140


Вычисление вклада от этих диаграмм (графиков «с лестницей») приводит к результату

А£„ =

4яег (Z — 1 )2 Хд

dq

 

X

(2л)3 р

IЯ1(Я2 + х%)

 

 

 

 

 

X

dpdp'$>„ (р) т|>* (р ') (Ек - вр) ^

(Р)

t k (Р ')

(Еп — е®) (ig4 —

 

 

(12.29)

 

— £ k +

Не +

Hi)

Нетрудно видеть, что при п=?^=к

АЕп = 0, а при п= к

АЕп = — (1/4) (Z — I)2 e2xD/e^.

(12.30)

Для выяснения физического смысла выражения (12.30) вычис­ лим собственную энергию свободного электрона

2 * 0 » ) = 7 Г ^ - У Л ^ С о(Р-9)[^*(Ч, <7*)-^(Ч)Ь 02.31)

(2л)3 р лаЛ J <?<

где

^*(q>

=

V3(q)ne(q. я*)

 

1+V(q)ne(q. Я*)

(12.32)

 

 

4 яе3 .

V(q)=-

Яг

Gq(P- q) = (ip* — ><7* + И- 8 p - q ) - 1 .

В сумме по <74 выделим член с <74 = 0 и рассмотрим подынтеграль­ ное выражение при малых q. Поскольку при <74 = 0 и q= 0, то

__

4 яе2

4 л е2 __

4 л е2х 2

 

q2 Х2

qi

qi (^2 |_ и2)

4ле

2х 2

dq

(12.33)

 

 

2 * (р) = (2л)3 Р J1 я2 (Я + * 2)

>Р« + — ер — ер_ „

Для свободного (невозмущенного) электрона полюса определя­

ются нулем выражения

\рц—ер+ ц = 0 . Поэтому

интеграл в вы­

ражении

(12.33)

следует

вычислять в с м ы с л е г л а в н о г о

з н а ч е н и я :

 

 

 

 

 

 

 

 

 

 

 

4 яе2х 2

 

dq

 

 

1

 

(12.34)

 

Уе(Р) =

 

J

q* (q* + х 3)

е р

~

е р—q

 

 

(2 я )3 р

 

Поскольку ер—e p_q = (l/2m) (2pq—q2), то

интеграл

 

 

dq

 

 

 

оо

 

 

я

sin Odd

 

 

 

 

 

4/тгя f

*2

f

 

- f

Я2 (<72 + *2)

 

ep-q

 

 

8 p

J

<?2 +

J

2pq cos 0 — я2

141


Далее

 

 

 

 

 

 

 

г

== —

f - dx

 

 

2p

 

sin 0rf0

 

 

Г

 

1 2pq cos 0

q*

2pq,JJ xx— q/2p

2pq

JJ xx—ql2pql:

+

0

 

 

—1

 

 

—1

 

r __L

Г

dx

_ _ L in 1 — qi2P

 

 

2pq

J

x — q/2p

 

2pq

1+ ql2p

 

■+ s

2P

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

00

In i -

g/2p

= 2я m p_1 j* •

dq

 

 

<7 (9 * + x!)

1 +

q/2Pj

oo

0

Подставляя этот результат в выражение (12.34), получаем

У4’ (р) = —е2к/4ерр.

Следовательно, А Е п в формуле (12.30) представляет собой поправку к собственной энергии «атома», т. е. иона с зарядом (Z—1) при его свободном движении. Действительно, усредняя выражение (12.30) по распределению Максвелла, получаем де­ баевскую энергию иона с зарядом (Z—1)

Д£ = — (1/2) (Z— 1)2е2к.

Этот сдвиг не зависит от состояния «атома». Отметим, что для обменных графиков, например диаграмм типа 4, можно полу­ чить выражения, сходящиеся при всех значениях q, и сдвиг имеет порядок (аох)2р-1, что согласуется с обычным вычисле­ нием обменного интеграла для гелия.

Из формул (12.25) и (12.26) следует, что сдвиг уровней про­ порционален корню из плотности и растет с температурой как Р-1/2. Сдвиг растет также с увеличением главного квантового числа. Эффект (в рассмотренном порядке по плотности) отсут­ ствует для водорода (Z = 1), что физически очевидно ввиду ко­ роткодействующего характера взаимодействия нейтрального ато­ ма водорода с плазмой. О величине эффекта можно судить из сравнения полученного сдвига с линейным штарковским рас­ щеплением

Дшт ~ пе2а0г0п2,

где п — главное квантовое число; г0— среднее расстояние меж­ ду частицами. Легко видеть, что

ДШт/Д£п ~ xr0 < 1.

142


Условие применимости теории возмущений имеет вид АЕп/Еп<С <С1. Следует отметить, что при выводе выражения для сдвига пренебрегалось eg~ P _1 по сравнению с Еп, поэтому получен­ ные формулы справедливы при достаточно низких температурах.

Напомним, что еще в § 5 делалась попытка интерпретиро­ вать результат (12.25) с помощью полуклассического рассмот­ рения эффекта динамического экранирования орбитального электрона. Конечно, столь грубое рассмотрение не позволило получить строгий количественный результат. Однако буквенное выражение (12.21) уже было получено на основе изучения флуктуаций плотности в приближении RPA [см. формулу (5.23)].

§ 13. ДВУХЧАСТИЧНАЯ ФУНКЦИЯ ГРИНА И УРАВНЕНИЕ СОСТОЯНИЯ ПЛАЗМЫ

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

Система частиц с короткодействующими силами

Рассмотрим систему одинаковых взаимодействующих частиц в объеме V в состоянии термодинамического равновесия при температуре р-1. Гамильтониан системы представим в виде

H = H 0 + H l t

(13.1)

где Н0— оператор системы невзаимодействующих частиц; Н {— гамильтониан взаимодействия:

Я! =

-у I г() (%х) Ф (х2) и (хх — х2)

ф (Хх) ф (х2) с1х^х2;

(13.2)

 

U (Хх

Х2) = U (Xj — х.2)

б (G /2)

 

— потенциал

парного

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

(мгновенного);

g' — кон-

 

 

 

 

л

станта взаимодействия. Отметим, что полевые операторы ф(х) и

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

ZN = Sp exp [—р (Я — рЯ)] = exp (—PQ),

то выражение £2 через двухчастичную функцию Грина имеет вид

Q = Q 0 „ .. ‘ -

f С2 (ft. ръ Ра, Pi) и (Pi Рг) dp,dp2dp3l (13.3)

2 (2я)9 р3 .)

g J

о

 

где Q0 — термодинамический потенциал идеальной системы.

143