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

Но бесконечное сгущение уровней имеет место лишь в модели уединенной частицы. В действительности электрические поля, создаваемые окружающими частицами, не позволяют реализоваться квантовым состояниям, начиная с некоторой степени возбуждения. Кроме того, как отмечено в разделе 14.1, конечная ширина спектральных линий приводит к их слиянию вблизи границы серии. На рис.14.1.4 видно, как линии, постепенно сближаясь, в конце концов, сливаются, и дискретный спектр переходит в континуум еще до границы ионизации.

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

17.2 Критерий Инглиса-Теллера

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

 

 

где <r> в данном случае - среднее расстояние между ионами. Оно определяется концентрацией частиц N:

 

 

Напряженность поля можно выразить как функцию N :

 

 

Число однозарядных ионов вследствие электронейтральности плазмы равно числу электронов. Поэтому в данном случае справедливо

 

^.

 

Атом водорода на энергетическом уровне с главным квантовым числом n будем рассматривать как электрический диполь, момент d которого равен

 

^.

 

В поле с напряженностью E атом приобретает потенциальную энергию

 

 

Индекс 'St' напоминает об эффекте Штарка. Спектральные линии перекрываются, если потенциальная энергия U достигает размеров энергетической щели ΔE_n,n+1 между двумя соседними уровнями. В случае больших значений n асимптотическая формула для ΔE_n,n+1 получается аналогично формуле (14.1.7) для Δλ:

 

Из условия

 

 

вытекает окончательная оценка для главного n_St, начиная с которого уровни энергии перекрываются:

 

^.

Также перекрываются и линии с верхним уровнем n>n_S. Формулу (2.2) называют критерием Инглиса-Теллера. Она определяет сдвиг границ спектральных серий.

17.3 Плазменная частота

Выше, оценивая влияние на атом окружающих частиц, мы неявно предполагали, что их можно описывать в рамках модели идеального газа. Хотя это предположение в большинстве случаев оказывается верным, но заранее оно совсем не очевидно. Действительно, кулоновский потенциал очень медленно уменьшается с расстоянием, и естественно ожидать, что много заряженных частиц непрерывно влияют друг на друга. Но в модели идеального газа частицы должны взаимодействовать только попарно, причем в течение короткого промежутка времени, пока происходит столкновение. Это, на первый взгляд, серьезное противоречие решается при более внимательном анализе поведения плазмы. Оказывается, что газ, состоящий из заряженных частиц, в целом электронейтрален. Кроме того, в нем возникают флуктуации заряда, в результате которых частицы противоположных зарядов экранируют друг друга, поэтому силы электростатического взаимодействия быстро падают с расстоянием. Время установления электронейтральности определяется ленгмюровской частотой ω₀, а пространственный масштаб - дебаевским радиусом l_D.

В первой главе мы написали формулу (1.1.2) для плазменной, или ленгмюровской частоты из соображений размерности. Теперь выполним расчеты в рамках модели плоского конденсатора. Хотя эта модель далека от реальных условий, она позволяет наиболее простым способом вычислить величину ω₀.

Пусть в пространстве между двумя плоскостями находится электронейтральная плазма. Как и выше, предположим, что плазма состоит из электронов и однозарядных ионов, поэтому справедлива формула (2.1). Хотя электрическое поле с одинаковой силой действует на ионы и электроны, но легкие электроны значительно более подвижны. Поэтому в первом приближении мы будем принимать во внимание только движение электронов. Предположим, что в некоторый момент все электроны, находящиеся в узком слое толщиной x, сместились относительно положительных зарядов и расположились на границе слоя, как это показано на рис.17.4.1.

 

 

Между плоскостями образуется плоский конденсатор с поверхностной плотностью заряда

 

^,

 

внутри которого возбуждено электрическое поле

 

^.

 

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

 

^.

 

По второму закону Ньютона она вызывает ускорение

 

^.

 

Мы получили типичное уравнение колебаний, решением которого является гармоническая функция. Ее частоту мы уже написали в формуле (1.1.2) первой главы при анализе размерностей. Приведем ее еще раз здесь:

 

 

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

17.4. Дебаевский радиус

Плазменная частота ω₀ является динамической характеристикой плазмы. Существует также статическая характеристика l_D, имеющая размерность длины. Она называется дебаевским радиусом, или дебаевской длиной экранирования и определяет масштаб, начиная с которого плазма становится электрически нейтральной. Электронейтральность плазмы как системы зарядов имеет следующий смысл. Внутри объема, занимаемого плазмой, каждый ион или электрон окружен 'шубой' из противоположно заряженных частиц. Суммарный заряд шубы в точности равен заряду частицы (с противоположным знаком), а ее размер (по порядку величины) - дебаевскому радиусу. Таким образом, система 'заряд-шуба' электронейтральна на расстояниях, превышающих длину экранирования.

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

 

^.

 

Частица может его преодолеть, если ее тепловая энергия, пропорциональная температуре T, сравнима с величиной барьера:

 

.

 

Отсюда следует оценка расстояния между обкладками:

 

^.

 

Точные расчеты подтверждают эту оценку, а именно, дебаевский радиус оказывается равным

 

 

Мы случайно получили точную формулу, хотя все рассуждения имеют силу только как оценки по порядку величины. Смысл параметра l_D заключается в следующем. При усреднении по достаточно большим расстояниям плотность числа носителей заряда N можно считать постоянной величиной, но на малых масштабах возможны флуктуации заряда. Сказанное иллюстрируется рис.17.4.1:

 

 

Одинаковое число знаков '+'и '-' во всех строках и столбцах матрицы имитирует квазинейтральность плазмы на масштабах порядка стороны квадрата. Но на небольших расстояниях знаки распределены неравномерно, что указывает на флуктуации заряда. Каждый заряд внутри квадрата окружен 'шубой', в которой число противоположных знаков на единицу больше.

Дебаевский радиус называют также 'длиной экранирования'. Смысл этого термина поможет уяснить решение следующей задачи. На левой границе x=0 полупространства, занятого однородной плазмой, находится электрод, создающий в этой плоскости положительный потенциал φ₀.

 

(4.2) φ(0) = φ ₀.

 

Под влиянием его электрического поля в области x>0 происходит поляризация плазмы, схематически изображенная на рис.17.4.2.

 

 

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

Определим масштаб длины экранирования, предполагая плоское одномерное распределение зарядов. Для этого запишем уравнение Пуассона для потенциала φ:

 

^,

 

создаваемого электронами и однозарядными ионами

 

^.

 

По предположению, они находятся в тепловом равновесии. Следовательно, их распределение по координате описывается формулой Больцмана:

 

 

Посредством N_i,e(x) мы обозначили плотность числа ионов и электронов как функцию координаты, а N, как и выше - некоторое 'фоновое', усредненное значение. Подставляя выражения для N_i(x) и N_e(x) в уравнение Пуассона, получим

 

^.

 

Однако в таком нелинейном виде пользоваться уравнением нельзя. Вспомним, что распределение Больцмана дает усредненное по времени значение концентрации. Если φ - соответствующее среднее значение потенциала, то в правой части (4.3) средние значения функций N_i(x) и N_e(x) оказывается замененным функцией среднего значения. Такая операция допустима исключительно для линейных функций, следовательно, уравнение (4.3) приобретает смысл только после линеаризации. Другими словами, оно справедливо только в случае малого по модулю отношения eφ/T. Раскладывая экспоненты в ряд с точностью до линейных слагаемых, приходим к уравнению

 

^,

 

в котором введено обозначение

^.

 

В случае плоской симметрии уравнение (4.4) имеет вид

 

^.

 

Его решением является линейная комбинация

 

 

Из требования конечности потенциала определяем константу С₂:

 

^,

 

а из граничного условия (4.2) - константу С₁:

 

^.

 

Окончательно получаем:

 

^.

 

В рассматриваемой модели длина экранирования равна

 

^.

 

Она в раз меньше, чем (4.1). Различие обусловлено тем, что в упрощенной модели принималось во внимание экранирование только электронами, а в (4.4) учтены также ионы.

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

 

^,

 

где α нумерует компоненты плазмы. Параметр Z_α равен заряду частиц сорта α в единицах элементарного заряда e. Например, для электрона . В новых обозначениях уравнение Пуассона (4.3) принимает вид

 

^.

После линеаризации с учетом (4.6) снова получаем уравнение (4.4) для самосогласованного поля, но с другим выражением для κ:

 

^.

 

Решение уравнения (4.4) в случае сферической симметрии выполняется методом, аналогичным тому, который был применен в разделе (16.5). Для функции

 

 

уравнение самосогласованного поля приводится к виду:

 

^.

 

Его убывающее на бесконечности решение равно

^.

 

Константа C теперь определяется из условия, что потенциал вблизи частицы, имеющей заряд q, должен стремиться к кулоновскому:

 

^.

 

Для электрона

 

^,

 

а для иона

 

^.

 

Итак, приходим к окончательному выражению

 

 

Последняя формула дает математическое описание 'шубы' размера l_D, состоящей из частиц с противоположным зарядом и схематически изображенной на рис.17.4.1. На расстояниях меньше дебаевского радиуса поле каждой частицы определяется кулоновским потенциалом, а при r > l_D оно практически полностью экранируется 'шубой'.

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

 

17.5 Снижение потенциала ионизации

Электрон, находящийся в основном состоянии, несомненно, принадлежит своему атому. Но по мере перехода на все более возбужденные состояния он удаляется от ядра, и на него начинают влиять соседние заряженные частицы. Электроны, окружающие атом, все сильнее экранируют поле ядра, и когда радиус электронной орбиты превышает масштаб экранирования, сила притяжения атомного электрона к ядру начинает падать с расстоянием по экспоненциальному закону (4.9).

Согласно разделу (16.4), бесконечное число уровней вблизи границы ионизации существует только в том случае, если потенциал падает с расстоянием не быстрее, чем 1/r². В поле (4.9) при r > l_D связанные состояния не реализуются, и электрон становится свободным при энергии возбуждения, меньшей потенциала ионизации изолированного атома. Происходит как бы снижение потенциала ионизации. Главное квантовое число последнего возбужденного состояния можно оценить из условия

 

^,

 

или

 

^.

 

В объектах таблицы 17.5.1 штарковское уширение уровней приводит к большему сдвигу границы серии, чем понижение потенциала ионизации: n_St<n_D. Оценим плотность газа, при которой оба эффекта дают один и тот же результат. Из условия

 

 

следует уравнение

 

^.

 

Его решением является плотность числа частиц

 

^.

 

Столь большая величина N_equal реализуется только для весьма ограниченного круга объектов.

Численные оценки

Выполним численные оценки описанных выше параметров плазмы для нескольких конкретных объектов: установки лазерного термоядерного синтеза (ЛТС), управляемого термоядерного синтеза (УТС) с магнитным удержанием плазмы, солнечной короны (СК) и солнечного ветра (СВ).

 

Таблица 17.5.1.

	ЛТС	УТС	СК	СВ
TэВ	104	104	100	10
Ne(см-3)	1022	1014	108	10
nSt		30	180	1500
nD	10	103	104	3ћ105
ω0, 1/с	6ћ1015	6ћ1011	6ћ108	2ћ105
2π/ω0, с	10-15	10-11	10-8	3ћ10-5
τ, с	10-9	1	103	8ћ104÷3ћ105
lD, см	10-6	10-2	1	103
L, см	10-2	102	1010 1011	1091013

 

В последней строке таблицы приведены пространственные масштабы систем, а τ - характерное время длительности процессов. Хорошо видно, что во всех случаях плазменная частота настолько велика, что флуктуации заряда за промежуток времени τ многократно сменяют знак и, тем самым, приводят к 'фоновому' значению плотности заряда. Сравнение двух последних строк показывает, что в перечисленных объектах представление об электронейтральной плазме соответствует действительности.

 

17.6 Ионизационное равновесие

В условиях термодинамического равновесия степень ионизации газа однозначно определяется его температурой и плотностью. Построение теории, необходимой для вывода уравнений ионизации и диссоциации, было проведено химиком ван'т Хоффом. В 1919 году Эггертом эти идеи были применены к веществу звездных недр, а в 1920 году Сахá установил, что формула для равновесной ионизации дает ключ к пониманию спектральной последовательности звезд. Исключительная важность этой работы Сахá для астрономии привела к тому, что уравнение носит его имя. Последовательный вывод формулы Сахá проводится методами химической термодинамики и выходит за рамки настоящего курса. Здесь мы покажем, как ее можно получить путем не вполне строгих, но физически обоснованных рассуждений.

Равновесные населенности невырожденных уровней N_n и N₀ описываются формулой Больцмана:

 

^.

 

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

 

 

В левую часть входит элемент объема в фазовом пространстве, введенный в формуле (1.2.3) первой главы. Выпишем его еще раз:

 

^.

 

Подставив (1.2.3) в (6.1), получим распределение свободных электронов по импульсам:

 

^.

 

Проинтегрируем это выражение по импульсам и координатам:

 

^.

 

В пространстве импульсов удобно перейти к сферической системе координат. С учетом изотропии интегрирование по угловым переменным сводится к умножению на 4π:

 

^.

 

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

 

 

Подставляя (6.3) в (6.2), получим результат интегрирования по пространственным переменным:

 

^.

 

Интеграл в правой части выражается через гамма-функцию. Напомним известную в математике формулу

 

 

В данном случае ν = 2 и α = (2m_eT)^-1. Следовательно,

 

^.

 

Учитывая, что

 

^,

 

окончательно получим следующую формулу:

 

^.

 

Изначально мы опирались на формулу Больцмана для невырожденных уровней. Для учета вырождения необходимо каждую из трех величин: N_e, N_i и N₀ - разделить на ее статистический вес. Это выполняется с помощью следующих замен:

 

^.

 

Здесь g_i и g₀, - статистические веса иона и атома в основном состоянии, а статистический вес электрона полагаем равным двум - числу ориентаций его спина. Произведя замены, получим формулу Сахá:

 

^.

 

Перепишем ее, выразив множитель в скобках через длину волны де Бройля (1.2.8):

 

 

Здесь D_D=λ_D/(π^1/2). Произведение определяет число квантовых состояний, связанных с пространственным перемещением свободных электронов. Возможны два случая:

1)      Это означает, что длина волнового пакета, значительно меньше среднего расстояния между электронами. В такой среде новые свободные электроны легко находят свое место, и степень ионизации газа может оказаться значительно больше экспоненты в правой части (6.2).

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

Присутствие параметра отличает формулу Сахá от формулы Больцмана. Возможна заметная степень ионизации химического элемента при сравнительно низких температурах. Например, водород заметно ионизован в атмосферах звезд класса А0, хотя температура там почти в пятнадцать раз меньше его потенциала ионизации.

17.7 Сумма по состояниям

 

Приведем без вывода точное выражение для формулы Сахá, полученное методами химической термодинамики:

 

^.

 

В знаменателе левой части стоит N_a - полная плотность числа атомов, просуммированная по всем его возможным состояниям. Величина N_i теперь равна сумме числа ионов во всех состояниях. В правой части статистические веса иона и атома в определенном (основном) состоянии заменены статистическими суммами:

 

^.

 

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

Расходимость суммы по состояниям

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

 

^.

 

Если температура низка в сравнении с потенциалом ионизации, то слагаемые суммы (7.2) быстро уменьшаются по величине и после немногих первых членов становятся пренебрежимо малыми. Может возникнуть впечатление, что сумма по состояниям сходится. Но это не так, в чем легко убедиться, заменив в (7.2) все энергии уровней на потенциал ионизации. Сумма ряда, если она существует, от такой операции может только уменьшиться. Следовательно, справедливо неравенство

 

^,

 

из которого вытекает расходимость суммы по состояниям.

Этот результат справедлив для всех ионов и атомов, так как их возбужденные состояния вблизи границы ионизации хорошо описываются в модели атома водорода. Причина расходимости ряда (7.2) ясна из предыдущего материала. Бесконечное множество уровней существует только в уединенном атоме, следовательно, сумма (7.2) содержит только конечное число слагаемых.

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

В задаче об излучении плазмы связанными следует считать электроны, испускающие дискретные линии, а свободными - электроны, дающие непрерывный спектр. Если электрон находится на одном из верхних возбужденных уровней, то у него есть определенная вероятность либо 'свалиться' на какой-нибудь из нижних уровней с испусканием кванта дискретного спектра, либо подвергнуться ионизации, в результате которой он окажется свободным и сможет испустить квант непрерывного спектра. Наиболее последовательным представляется решение вопроса методами физической кинетики. Находятся вероятности переходов между дискретными уровнями, вероятность излучения с переходом на все нижележащие уровни и вероятность ионизации с учетом всех возможных процессов. К ним относятся удары электронов и быстрых ионов, а также ионизация тепловым и резонансным излучением, воздействие флуктуирующих микрополей. С учетом всех перечисленных процессов вычисляются заселенности высших возбужденных уровней. В расчет входят величины, характеризующие уже не равновесие, а вероятности различных процессов. В последующих главах мы рассмотрим эти процессы, а также покажем сравнительно простой способ вычисления всех необходимых вероятностей.

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

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