Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Князев Б.А. Низкотемпературная плазма и газовый разряд. 2000.pdf
Скачиваний:
123
Добавлен:
26.08.2013
Размер:
1.65 Mб
Скачать

Лекция 4

1.4.Излучательные процессы в низкотемпературной плазме

1.4.1.Роль излучения в низкотемпературной плазме, классификация переходов

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

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

En = E1/n2,

(1.4.1)

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

En(s) = I 1

1

.

(1.4.2)

n2

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

Рис. 21: Классификация атомных переходов.

56

Рассмотрим сначала чисто излучательные переходы (рис. 21). Их можно разделить на связанно-связанные, связанно-свободные, свободно-связанные и свободно-свободные. Примером последних являются тормозное излучение и тормозное поглощение при движении электрона вблизи атома. Связанно-свободные и свободносвязанные переходы происходят при захвате электрона и фотоионизации атома, соответственно. Связанносвязанные можно разделить на спонтанные и вынужденные. К последним относятся как поглощение фотона, так и его вынужденное излучение, индуцированное пролетающим рядом первичным фотоном. В последнем случае, как известно, излученный фотон имеет ту же фазу, что и первичный. В этом случае мы говорим о когерентном излучении.

1.4.2.Тормозное излучение и поглощение

Рассмотрим свободно-свободные переходы. Тормозное излучение "ответственно"за излучение плазмы в непрерывном спектре (континуум). Следуя [23] выполним качественный анализ спектра тормозного излучения. Последовательность действий будет следующей. (а) Сначала вычисляется высокочастотная асимптотика классического спектра рассеяния электронов, имеющих скорость v, на точечном кулоновском центре с зарядом Ze (формула Крамерса). (б) Вычисленный для всех ω спектр принимают за "эталонный".

(в) Затем вводят т.н. фактор Гаунта g(ω), который позволяет получить истинный спектр, учитывающий что электрон взаимодействует с "экранированным"ядром, и Zi << Zeff < Z.

Высоким частотам в спектре соответствуют, очевидно, близкие пролеты ρ << a, где a = Ze2/mv2

– характерная “кулоновская” длина. При этом излучаемая частота ωeff по порядку величины равна угловой скорости вращения электрона вблизи иона ωrot. Используя законы сохранения энегии и момента импульса

 

 

 

 

 

 

mv

2

 

 

 

2

 

 

 

 

 

 

 

Ze

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

mvmax

 

,

 

 

 

 

 

 

 

(1.4.3)

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

r0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mvρ = mvmaxr0 ,

 

 

 

 

 

 

 

 

 

 

(1.4.4)

и пренебрегая для близких пролетов левой частью в (1.4.3), получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vmax

=

ρv

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.4.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2v2

=

Ze2

 

,

 

 

 

 

 

 

 

 

 

 

 

 

(1.4.6)

 

 

 

 

 

 

 

 

 

 

 

2r02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r0 =

 

 

 

 

ρ2

 

 

 

 

 

ρ2

 

 

 

 

 

 

 

 

 

 

(1.4.7)

 

 

 

 

 

 

 

<< ρ .

 

 

 

 

 

 

 

 

2 · (Ze2/mv2)

2a

 

 

 

 

 

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

 

 

r0

 

r02

 

 

ρ4

 

 

 

 

ρ

3

· a

 

ρ

 

3

 

 

ωrot =

vmax

=

 

ρv

=

ρv

· 4a2

= 4

 

 

 

 

a

 

 

 

 

 

v

 

4

 

a

 

ω>>ω ,

(1.4.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

57

где

 

 

 

 

 

v

 

mv3

(1.4.9)

ω=

 

=

 

a

Ze2

представляет собой характерную частоту спектра излучения, соответствующую ρ = a.

Теперь нетрудно оценить энергию, излучаемую в дипольном приближении электроном, имеющим ско-

рость v и пролетающим мимо иона с прицельным расстоянием ρ

 

 

+2 e2

 

 

 

 

 

 

 

 

2e2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E(ρ) = −∞

 

 

 

 

w2(t)dt

 

 

 

 

(wmax)2 · (∆teff) ,

 

 

(1.4.10)

 

3

c3

3c3

 

 

где максимальное ускорение электрона равно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w2

 

 

=

vmax

=

 

 

Ze2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.4.11)

 

 

 

 

 

 

 

 

 

2mr02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max

 

 

 

r0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а характерное время пролета можно оценить как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(∆t)eff

r0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.4.12)

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя (1.4.5), (1.4.7), (1.4.11) и (1.4.12), преобразуем энергию к следующему виду

 

2

 

 

 

 

vmax

 

r0

 

 

vmax

 

 

 

 

 

 

ρ

3

v

3

 

 

 

 

16a v

3

 

 

(wmax) · (∆teff)

4

 

 

 

 

 

 

 

 

=

 

 

 

3

 

 

 

=

 

 

 

 

=

 

4

 

(1.4.13)

 

 

 

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

r02

 

vmax

 

 

 

r0

 

 

 

r04

 

 

 

ρ5

 

Домножив энергию на соответствующее дифференциальное сечение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ρ) =

16a4v3

 

 

2e2

=

 

2 · 16

 

e2a4v3

,

 

 

 

 

 

(1.4.14)

E

 

 

ρ5

 

· 3c3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

c3ρ5

 

 

 

 

 

 

 

 

получим величину

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ρ) = ∆

E

(ρ)

·

2πρdρ = 2π

·

 

16 · 2

 

e2a4v3

 

,

 

 

(1.4.15)

 

 

 

ρ4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

c3

 

 

 

 

 

 

 

 

имеющую размерность [эрг·см2]. Смысл ее ясен, если учесть, что умножив ее на поток электронов nev и

плотность инов ni, мы получим мощность тормозного излучения из единицы обьема dP (v, ρ) для электронов с прицельным расстоянием ρ и скоростью v.

Используя введенное выше приближение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

3

v

4a2v

 

 

ω(ρ) ωrot = 4

 

 

 

 

 

=

 

 

,

(1.4.16)

ρ

a

ρ3

нетрудно перейти от переменной ρ к ω

 

 

 

 

 

 

 

 

 

 

 

 

 

|dω| = 12

a2v

 

 

 

 

 

(1.4.17)

 

dρ .

 

 

 

ρ4

 

 

 

Тогда, приняв 12 4π, получим выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

(ω) =

16 e2a2v2

· dω = f(ω) ,

 

(1.4.18)

3

 

c3

 

 

58

не зависящее от частоты. Окончательно

 

 

 

 

 

 

 

 

 

(ω) =

16 e2(Ze2)2v2

=

16 Z2e6

dω .

(1.4.19)

 

 

 

 

 

 

3 c3(mv2)2

3 c3m2v2

 

 

 

 

Строгая классическая теория [25] дает для спектральной интенсивности ту же зависимость, но с допол-

 

 

 

 

 

 

 

 

 

 

нительным коэффициентом π/

3

:

 

 

 

 

 

 

 

 

 

 

 

(ω)

=

16π Z2e6

(1.4.20)

 

 

 

 

3

 

 

 

.

 

 

 

c3m2v2

 

 

 

3

Эта формула называется формулой Крамерса. Выражение (1.4.20) справедливо, как это следует из его вы-

вода, для v << vmax (см. (1.4.3)), то есть для частот

 

 

 

 

mv3

(1.4.21)

ω >>ω=

 

.

2

 

Ze

 

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

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

можно найти в [23]. В пределе малых частот (дальние, почти прямолинейные пролеты), когда ω <<ω= mv3/Ze2, гаунт-фактор, например, есть

g(ω) ln ωω . (1.4.22)

По порядку величины гаунт-фактор g(ω) для тормозного излучения лежит в интервале 1 ÷4 (см. рис. 22).

Рис. 22: Гаунт-факторы для свободно-свободных и свободно-связанных переходов [24].

Размерность (ω)/dω для дальнейшего удобно переписать в расчете на единичный спектральный интервал, выраженный в с1, а именно [(эрг·см2)/1)]. Проинтегрировав теперь это выражение по распределению электронов по скоростям, получим спектральную плотность мощности тормозного излучения

[эрг/с· с1·см3]

 

 

(ω)

 

 

 

dP (ω)

= 0

ni nevf(v) 4πv2 dv .

(1.4.23)

 

 

59

Akm1 .
|Dkm|2 gm ,

1.4.3.Линейчатое излучение. Вероятность перехода. Сила осциллятора

Вероятность перехода между дискретными уровнями k и m атомной системы (коэффициент Эйнштейна для спонтанного излучения) равна

4ω3

Akm 1] = 3 km (1.4.24)

c3

где gm – статвес конечного состояния, Dkm – матричный элемент дипольного момента атомной системы

N

D = ers ,

s=1

а rs - координата s-го электрона атома относительно ядра. Частота излучения, испускаемого или поглощаемого при переходе между этими уровнями, равна

ωkm =

Em − Ek

.

(1.4.25)

 

 

 

Время жизни уровня k по отношению к радиационному распаду находится суммированием обратных вероятностей излучения на все нижележащие уровни

τ =

m

Введем Fω - форму спектральной линии, пронормировав ее так, что ∫ Fω= 1. В отсутствие внешних возмущений Ее ширина ∆ω, в соответствии с принципом неопределенности, пропорциональна τ1 = Akm. Эта “естественная” ширина линии для разрешенных (дипольных-дипольных) переходов равна по порядку величины ∆ω 108 с1, что много меньше частоты перехода ω0 1015 с1. Линия имеет лоренцовский

контур

 

γ

1

 

 

 

Fωест = εω =

,

0

ε(ω) = 1 ,

(1.4.26)

 

 

 

2π

(ω − ω0)2 + (γ/2)2

где γ – постоянная затухания. Переходя к сечению поглощения фотона, для классического случая вынужденных колебаний упруго связанного электрона в поле электромагнитной волны E = E0 sin ωt найдем

σω =

πe2

 

γ

,

(1.4.27)

mc

 

(ω − ω0)2 + (γ/2)2

где γ = ∆ω – коэффициент затухания.

Для естественного уширения (1.4.27) γ = Akm, и сечение поглощения в центре линии равно

σmax 4πe2 109 см2,

ω0 mcγ

что много больше газокинетического сечения πa2. В классическом приближении площадь под спектральной линией осциллятора даже при наличии столкновений является константой, не зависящей от ширины линии

γ. В расчете на один осциллятор, в соответствии с общепринятой традиционной нотацией, она равна [2]

 

 

1

 

 

πe2

 

σω=

 

σω=

 

= 2.64 · 102 см2/с .

(1.4.28)

2π

mc

60

Действительно, при резонансном возбуждении осциллятора (ω = ω0) амплитуда колебаний растет пропорционально времени от последнего столкновения ra(t) t. Между столкновениями набирается энергия

W ra2 τc2, которая при следующем столкновении переходит в хаотическую. Тогда энергия отбираемая

от поля за секунду, которая пропорциональна величине σω, есть

 

 

W

 

ra2

 

τc2

= τc.

(1.4.29)

 

τc

τc

τc

Поскольку ∆ω τc1, то ∫ σωdν σωmax · ω (W/τc) · ω =const.

Если атом взаимодействует с окружающей средой, то это приводит к дополнительному уширению спектральной линии. Механизмы уширения мы рассмотрим в следующем разделе, а сейчас получим полезные соотношения, описывающие вероятность переходов для линий с некоторым реальным сечением σω. Энергия, поглощаемая на данном переходе в 1 см3 за 1 с, вне зависимости от характера уширения линии равна

4πnm Iωσωdω ,

где Iω [эрг/с см2 с1 ср] – “интенсивность” излучения, а S = Iω(Ω)dωdΩ – физически измеримая величина, называемая плотностью потока энергии. Если Iω const в пределах линии, то ее можно вынести из под интеграла и площадь линии служит мерой поглощательной способности вещества в линии.

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

ного излучения и распределение Больцмана для заселенностей уровней, получим

 

 

gk π2c2

 

 

gk πc2

(1.4.30)

σω =

 

 

 

Akm · Fω =

σω=

 

 

 

Akm .

gm

ω2

2gm

ω2

Пронормировав выражение (1.4.30) на площадь классического осциллятора (1.4.28) мы получим величину, называемую силой осциллятора, которая, по соглашению, для поглощения считается положительной

fmk =

(∫ σω)

=

gk mc3

Akm .

(1.4.31)

 

 

 

πe4/mc

2gm e2ωkm2

 

 

 

 

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

−gkfkm = gmfmk.

Величины fkm подчиняются правилу сумм

fkm = N ,

m

где N – полное число электронов в атомной системе, причем суммирование ведется по всем состояниям дискретного и непрерывного спектров. Силы осцилляторов практически всех линий протабулированы. Их можно найти, например, в [4, 26, 27, 28]. Для примера некоторые из них приведены в Таблице 6.

61

Элемент

H

H

N

He

Na

Hg

 

 

 

 

 

 

 

λ, A

6563(Hα)

 

 

5876

5890

5461

переход

2s − 3p

2p − 3s

2p − 3d

2p3p0 3d3D

3s − 3p

6s6p − 6s7s

f

0.416

0.014

0.695

0.62

0.98

0.08

 

 

 

 

 

 

 

Таблица 6: Силы осцилляторов некоторых переходов.

Рис. 23: Вероятность излучательных переходов с уровня k = 15 водородоподобного атома на нижележащие уровни m.

Вероятность перехода можно теперь записать в виде

 

 

 

 

 

 

 

A

=

2e2ωkm2 gm

f

 

= 3.79

·

109

(Ek − Em)2

f

km|

.

(1.4.32)

 

 

Ry2

km

 

mc3gk

mk

 

 

|

 

 

Для водородоподобных атомов fkm описывается квазиклассической формулой Крамерса, погрешность которой составляет 30% при всех k и m,

 

32

 

1

 

 

1

 

 

 

 

 

 

1

 

 

 

1

 

 

1

 

 

 

Ry

3

 

 

 

 

 

 

 

 

 

 

 

 

3.9

 

 

 

 

 

 

.

(1.4.33)

fkm

3π3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k5

m3

(1/k2

 

 

1/m2)3

k5

m3

 

k

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

− E

 

 

 

Подставив (1.4.33) в (1.4.32) и учитывая множитель Z4, получим для водородоподобных атомов

 

 

 

 

 

 

 

 

 

 

Akm =

 

1.6 · 1010Z4

.

 

 

 

 

 

 

 

 

 

(1.4.34)

Для k = 2, m = 1, Z = 1 получим

 

 

 

 

 

 

k3m(k2 − m2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

=

1.6 · 1010

 

 

=

160 · 108

= 5

·

108

c1 .

 

 

 

 

21

 

 

8

·

1(4

1)

 

32

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Согласно табличным данным [16], A21 = 6.26 · 108 с1 (g2/g1 = 6/2). Видно, что совпадение теории и эксперимента достаточно хорошее для оценок. Для переходов на нижние уровни (m << k)

Z4

Ak,m<<k 1.6 · 1010 k5m .

62

Вероятность максимальна при переходе в основное состояние (m = 1). Для перехода на соседнее нижнее состояние вероятность в два раза меньше

A

=

 

1.6 · 1010z4

 

0.8

 

1010

Z4

.

(1.4.35)

k3(k − 1)[k2 − k2 + 2k − 1]

·

 

k,k−1

 

 

 

k5

 

Экстремальное значение Ak, найдем, приняв нулю производную по m знаменателя формулы

 

[k3m(k2 − m2)] = k3[(k2 − m2) + m · (2m)] = k3(k2 3m2) .

(1.4.36)

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку вторая производная 6k

m < 0, при m = k/ 3 Akm принимает минимальное значение:

 

 

 

A

 

 

1.6 · 1010Z4

.

 

 

 

 

 

(1.4.37)

 

 

 

 

 

 

 

 

 

 

 

 

k,k/3

k5 · 0.38 k

 

 

 

 

 

 

Зависимость вероятности перехода от номера нижнего уровня, вычисленная из (1.4.34), приведена на рисунке (23). Если зафиксировать нижнее состояние, то вероятности переходов с высоких уровней значительно меньше, чем с ближайших, что демонстрируется рисунком 24). Суммарную вероятность радиацион-

Рис. 24: Вероятность излучательных переходов на уровень m = 1 водородоподобного атома с

верхних уровней k.

ного распада уровня k можно оценить (см. [3]) как

Ak =

Akm = 1.6 · 1010Z ln k3

− k .

(1.4.38)

 

 

 

 

 

 

 

 

 

 

k5

2

 

 

 

m<k

 

 

 

Приведенные выше выражения для Akm и fkm можно использовать и для неводородоподобных атомов, если использовать в них эффективные квантовые числа

k =

Z2

Ry

(1.4.39)

 

 

,

 

 

 

Ek

 

где Ek – реальная энергия связи данного уровня многоэлектронного атома.

63