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

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

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

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

Добавлен: 15.10.2024

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

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

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

изменения интенсивности в максимумах Д /= (/Син—/кр)//син и соответствующие экспериментальные значения, а также сдвиг максимумов относительно центра линии (АХ) [7].

Зависимость А/ от плотности я* в интервале 1016—1017 см~3 оказывается удобной для практического использования моно­ тонной функцией. Следует отметить, что относительное измене­ ние интенсивности в максимумах крайне слабо зависит от спо­ соба усреднения. Поэтому метод измерения плотности по асим­ метрии линии Нр может оказаться более перспективным, чем метод диагностики плазмы по ширине линий, учитывая неточ­ ность измерения последней и несовершенство существующих теорий для ширины линий при умеренных плотностях плазмы. В работе Грима [15] делается попытка объяснить асимметрию водородных линий с помощью квадратичного штарк-эффекта. Очевидно, что квадратичный штарк-эффект дает вклад следую­ щего порядка по отношению a0IR0 по сравнению с квадрупольным членом и не является основным эффектом.

Изложенная теория не позволяет судить о центре линии, но это нас не интересует, поскольку максимумы интенсивности рас­ положены достаточно далеко от ее центра, так что отношение

сдвига максимума к

ударной

ширине у

при Т ~ 1 эв

и яг~

— 1017 смг3

 

 

 

 

 

 

 

 

 

 

А т

___ _ 6

R 0

^ g o vv

|

 

 

 

 

Яу

п2

 

 

 

 

Дальнейшее

улучшение теории

должно

касаться

функции

распределения

И7 (со)

и учета

электронного

уширения

линии.

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

Двухчастичная функция Грина и ширина спектральных линий атомов и ионов в плазме

В настоящем разделе распространим метод функций Грина на случай вычисления ширины спектральных линий в плазме с учетом корреляции заряженных частиц [8]. Такой подход, повидимому, наиболее корректен. Удобство метода состоит в воз­ можности суммирования бесконечного числа существенных чле­ нов, необходимость учета которых обусловлена наличием дальнодействующих кулоновских сил.

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

ниан плазмы Нпл можно записать в представлении вторичного квантования в виде

121


Яп

/ V

 

/ Ч

/ ХЧ

V 4 .

VN

*4 . -'Ч .

»"Ч \

 

Я 0

 

Н q =

1 8р

flp

Лр ~ h бр ^ р

^ р ' >

 

 

 

 

р

 

 

 

 

 

 

 

Hi =

4 ' S 1/qI ] ^ ^

 

a„-_q ap+q +

 

} 0 1.13)

 

 

 

PP'

 

 

 

 

 

 

 

 

*ч /s

%Zctp OCp'^Xp'— q £Zp_|_q ) ,

 

~\~ 2^ 0tp

CCp'OCp'—q ^ p + q

 

где ap— оператор

поглощения

электрона;

ap~— оператор его

рождения;

Л -I-

Л

 

 

 

операторы

для

ионов с

а р , а р— соответствующие

импульсом

р; Z — заряд ионов, вер

и г р — операторы

кинетиче­

ской энергии электронов и ионов соответственно.

Энергия взаимодействия атома с заряженной частицей +2с, находящейся на расстоянии R от атома, в дипольном прибли­

жении имеет вид

—- Ze2r cos Q/R2,

 

 

и (г, R) =

 

 

где г — координата атомного

электрона.

Плотность

энергии

взаимодействия атома с плазмой записывается аналогично:

и = er cos 0 [Ztii (R, t) ne(R, t)]

U (R),

(11.14)

 

dR

 

где Пи ne— плотность ионов

и электронов

в плазме;

U(R) =

= e/R — кулоновский потенциал. Для дебаевской плазмы U мож­ но считать возмущением. Вероятность безрадиационных перехо­ дов атома А при его взаимодействии с плазмой вычислим в борновском приближении, критерий применимости которого имеет вид [9]

| и (R) | < hv/R.

В случае дипольного взаимодействия атома с плазмой это нера­ венство можно записать следующим образом:

e2/hv < R/ant

где R — характерное расстояние между частицами; ап — размер

атома.

заряженных частиц в плазме я-~1017 смтг,

При плотности

(3_1~ 5 эв имеем

для электронов: е2///щ—4; Rja0 — 250. Следо­

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

Вероятность того, что под действием возмущения плазма пе­ рейдет из состояния п в т , а атом из состояния к в / , есть

w = — |< т , / | и | п, к> |2 б (ЕтЕп -| - ек — + e6g'), (11.15)

П

122


где б-функция выражает закон сохранения энергии; egg' — изме­ нение движения центра масс.

Матричный элемент перехода

Г . dU

М = <ш , / I и I п, k > = J фг -Q^V/-4i — ne]mindr,

4>к

где ф* и фк — волновые функции атома в состояниях / и k; d —

дипольный момент атома. Записывая плотность п во вторичном квантовании, получим в ^-представлении

 

wa =

 

Z 2 3

ata.

 

■d r t‘ (It

JП

X

 

 

 

 

р'

 

 

 

 

 

 

 

 

X 6 ( £ m Ea - f

ek -|

e, +

Egg'),

 

(11.16)

где

~

1

f

dx$\dv\ dx — проекции

дипольного

момента d

M;,k=—

l

на

направление

q; Дч = 4 ne2/q2. Для

получения

полной вероят­

ности перехода ek->-e; нужно просуммировать

выражение

(11.16) по

всем

конечным состояниям плазмы

и усреднить по

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

 

 

 

 

 

п = exp [Р (й + |л(- V; + |xeNe— Н)],

 

(11.17)

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

где

^ = ^ п

2q№ V < M ® )>

(11.18)

 

ttp'OCp'—q Cl]+ aB'а I

 

Фч(со)=

Z,

У\

 

 

Р'

 

 

х exp [р (Й +

ntNi +

neNe— £„)] б (Е„ Еп

со + egg'),

 

 

СО=

8к — 8;.

(11.19)

При интегрировании W q по q получаем полную вероятность пе­ рехода к-»-/ или обратное время жизни в состоянии к по отно­ шению к такому переходу, что и определяет ширину линии

(к—>- /):

Yw = fUVq.

(11.20)

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

123


члены, учитывающие изменение состояния ионов, в выражениях (11.16) — (11.20) следует опустить и можно пренебречь изме­ нением энергии поступательного движения атома egg' вследст­ вие электронного удара. Для нахождения функции Фч(со) можно выразить ее через двухчастичную временную функцию Грина

/С(гх, t\\

r2, t2). Поскольку внешнее поле отсутствует, К зави­

сит лишь

от разностей г= Г]—r2, t = t\_i2. Поэтому Фд(со) вы­

ражается через корреляционную двухчастичную функцию, ко­

торая представляет собой частный вид гриновской

функции, у

/S

которой координаты и времена операторов ф и ф+ попарносовпадают. Вместо К удобно рассматривать запаздывающую функцию ^(q, о)), аналитическую в верхней полуплоскости пе­ ременной (о. Тогда [10]

К (q, со)

( 11-21)

O q N = П [1 — ехр (—р<о)]

Функция К является аналитическим продолжением темпе­ ратурной корреляционной функции /С(q, т) на верхнюю полу­ плоскость переменной и. Зная ^(q, со) на действительной оси, можно по формуле (11.21) найти Фд(со). Отметим, что способ вычисления вероятности перехода, при котором взаимодействиеатома с плазмой рассматривается как парные столкновения электронов с атомом, соответствует тому, что при вычислении функции К учитывается лишь простейший график типа петли

n fa. (°) = -

n p + q / 2 — ” p - q / 2

(11.22)

(pq) — co — i8

 

где Hp=exp[p(p—ep)] в случае статистики Больцмана.

Для правильного учета «парных» взаимодействий в случае кулоновских сил следует суммировать бесконечную цепочку из

электронных петель П (q, со). При этом

 

 

 

АГЧ= Пч/(1 — б'дПд).

 

(11.23)

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

 

 

 

 

 

262q (Mbu)^* т_

П (q,

со)

(11.24)

U 4 —

1 — ехр (— рсо) Im

1 - 6 qn

(q,

со)

Ширина ук, /, соответствующая

переходу к о п р е д е л я е т с я

формулами (11.20)

и (11.24). Следовательно,

 

 

8е4

Vk, i 1 — ехр (— рсо) (Afk. ifJ,

где

П (q , со)

= f q2dq Im

q"- — 4ле2П (q, co)

(11.25)

(11.26)

124


Разобьем область интегрирования на две части:

(0, q\)

и

(qi, оо),

причем

Тогда для второй области

можно пре­

небречь

членом

4яе2П в знаменателе выражения

(11.26)

по

сравнению с q2. Это соответствует приближению чисто парных

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

выражением

для

мнимой части П (Ч,со) [10]

 

 

 

 

 

 

 

Im П (q, со) = п

(2^

 

Р) - [1— ехр (—0со)] X

 

 

(

Р

(

а>т'!г

?

М

,

отч

X ехр |

2

\

q--------2 т */.;Ь

<П -27>

ТО

 

 

 

 

 

 

 

I = (2ятр)'/2л (1 — е_рм) epto/2K0

+

 

 

+ 4я Г

 

П 2(Ч 'Ю>

^

,

(11.28)

.)

 

 

дЧпе2П (q,

со)

 

 

 

где Ко — функция Бесселя;

л плотность

электронов. Вторым

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

yki, ж 8яе4 (2т$)и (л//г2) exp (ftcok> $ /2) X

 

X К0(Й<ок., р/2) (Мк,,)2.

(11.29)

Для получения полной ширины уровня ук следует просумми­ ровать выражение (11.29) по всем конечным состояниям I, для

которых матричные элементы Мг_к отличны от нуля:

ук = 8яе4

(2лгР)'/г

(MkJ)2ехр (7гсйк,,;Р/2) X

 

 

I

 

 

 

 

 

 

ХК0(Йсок, ,р/2) gt ехр (—РЕ,),

 

(11.30)

где g(exp(—рEi) — статистический вес

 

состояния

/. В

случае

большой энергии

возбуждения

атома

(Йсок,г> р - ')

получим из

выражения (11.29)

 

 

(

Я )2 1 .

к

 

ук , = 8яе4(2лг)‘/2— •

(11.31)

 

 

 

(Afflk. l)'1’

 

 

что согласуется с известной формулой Эльверта [14].

В обратном предельном случае, когда Йсок, г'Ср-1, выраже­ ние (11.29) логарифмически расходится и следует учитывать экранирование заряженных частиц в плазме, т. е. корреляцию

125