Физика плазмы, 2020, T. 46, № 8, стр. 730-745

Корреляция между зондовыми и спектроскопическими измерениями температуры в индуктивно связанной плазме смеси неона с аргоном

Z. I. Khattak a*, A. W. Khan b, A. U. Rahman c, M. Shafiq a**

a Department of Physics, Quaid-i-Azam University
45320 Islamabad, Pakistan

b Department of Physics, Gomal University
29050 D. I. Khan, Pakistan

c Department of Physics, International Islamic University
44000 Islamabad, Pakistan

* E-mail: zikhattak@phys.qau.edu.pk
** E-mail: ii.mshafiq@qau.edu.pk

Поступила в редакцию 23.09.2019
После доработки 21.01.2020
Принята к публикации 26.03.2020

Полный текст (PDF)

Аннотация

Исследуются параметры плазмы, включая температуру возбуждения $({{T}_{{exc}}})$, электронную температуру $({{T}_{e}})$, толщину скин-слоя $({{\delta }_{p}})$, потенциал плазмы $({{V}_{p}})$ и функцию распределения электронов по энергиям (ФРЭЭ) в двух модах поддержания стимулированного магнитным полюсом индуктивного ВЧ-разряда (MaPE-ICP) в смеси неона с аргоном. Исследования проводятся с помощью зонда Ленгмюра (LP) и оптической эмиссионной спектроскопии (OES). Основная цель работы – изучить возможность замены зондового измерения электронной температуры $(T_{e}^{{LP}})$ на спектроскопические измерения температуры возбуждения $(kT_{{exc}}^{{OES}})$ и электронной температуры ($kT_{e}^{{OES}}$). Измерения выполняются при различных значениях вкладываемой в разряд энергии и разном содержании аргона в смесях. Показано, что температура, определенная разными методами, уменьшается с ростом энерговклада и содержания аргона в смесях. При этом значения, полученные спектроскопически, меньше температуры, измеренной зондовым методом. Проведенное исследование может быть полезно в технологиях с химически активной плазмой, где использование зондовых методов имеет многочисленные ограничения.

Ключевые слова: плазма в смесях Ne–Ar, зонд Ленгмюра, оптическая эмиссионная спектроскопия, электронная температура, температура возбуждения, толщина скин-слоя, ФРЭЭ

1. ВВЕДЕНИЕ

В физике и химии плазмы наиболее важным, информативным и необходимым параметром является температура электронов, поскольку они не только участвуют в процессах возбуждения и ионизации нейтральных частиц, но и влияют на химические реакции в плазме [1]. Условия локального термодинамического равновесия (ЛТР) обычно не достигаются в низкотемпературной плазме, которая обозначается как неравновесная плазма (плазма без ЛТР). В такой плазме температура тяжелых частиц, как правило, слишком мала, чтобы влиять на химические процессы. Поступательную температуру свободных электронов часто связывают с температурой возбуждения, характеризующей связанные электроны в атомах и молекулах, поскольку свободные электроны контролируют процессы возбуждения, приводящие к распределению по возбужденным состояниям [2]. Температура возбуждения часто измеряется в плазме при высоких давлениях, и электронная температура постоянно оценивается в плазме атмосферного давления на основе измерения температуры возбуждения, что обусловлено их близкой связью [3, 4].

Для анализа параметров плазмы, включая ее потенциал $({{V}_{p}})$, плотность плазмы $({{n}_{e}})$, электронную температуру $\left( {{{T}_{e}}} \right)$ и функцию распределения электронов по энергиям (ФРЭЭ), используются различные диагностические методы. Определение параметров плазмы с целью изучения и оптимизации процессов в низкотемпературной плазме играет важную роль в плазменной обработке материалов [5]. Для этих целей часто используются электростатические зонды Ленгмюра, которые привлекают своей простотой, несмотря на сложную интерпретацию результатов и воздействие на саму плазму [6, 7]. В этом методе в плазму вводится тонкий проволочный электрод, к которому прикладывается дополнительное напряжение относительного референсного электрода, и измеряется ток через этот дополнительный электрод. Из измеренных характеристик ленгмюровского зонда можно определить различные параметры плазмы. Зондовый метод позволяет получать более надежные данные по электронной температуре в плазме без химических реакций (инертные газы) благодаря хорошо развитой теории зондов по сравнению с оптической эмиссионной спектроскопией, основанной на многих предположениях [2]. Зонды широко используются для диагностики плазмы, хотя при этом могут возникать определенные трудности в эксперименте, например, вызванные загрязнением зондового электрода или связанные с вызванными зондом возмущениями плазмы [6, 8].

Метод оптической эмиссионной спектроскопии, в отличие от зондового метода, позволяет исследовать плазму без ее возмущения. С его помощью можно определять несколько параметров плазмы [9, 10]. В этом методе измеряется интенсивность излучения на различных длинах волн при снятии возбуждения частиц в плазме. Часть информации извлекается с помощью численных кодов спектрального анализа, а часть – с помощью базы данных NIST [11]. С помощью метода больцмановских построений можно определить температуру возбуждения, а при помощи модифицированной схемы построения, изложенной в [12, 13], можно вычислить и электронную температуру.

Индуктивно связанная плазма в стимулированном магнитным полюсом ВЧ-разряде (MaPE-ICP) может существовать в двух модах: E-моде (известной также как электростатическая или емкостная мода) и H-моде (также известной как электромагнитная или индуктивная мода). Эти моды отличаются разными разрядными характеристиками [14, 15]. Разряд переходит в индуктивную или H-моду при увеличении подводимой ВЧ-мощности. Для емкостной или Е-моды характерны слабая интенсивность излучения света, высокая температура электронов, низкая плотность плазмы и ее высокий потенциал. В случае Н-моды плазма характеризуется более высокой интенсивностью излучения, высокой плотностью плазмы, низким потенциалом и небольшой температурой электронов [15, 16].

Во многих исследованиях оценивались электронная температура и температура возбуждения и изучалась их корреляция для разрядов в отсутствие локального термодинамического равновесия. Значения этих температур измерялись в плазме емкостного ВЧ-разряда, индуктивного ВЧ-разряда и плазме магнетронного разряда постоянного тока в широких диапазонах давления газа и подводимой мощности [2]. Обе температуры уменьшались с ростом давления и не сильно зависели от увеличения мощности. При этом температура $T_{{exc}}^{{OES}}$ стремилась к $T_{e}^{{LP}}$ с увеличением мощности и давления. Связь между $T_{{exc}}^{{OES}}$, определенной с помощью оптической эмиссионной спектроскопии, и $T_{e}^{{LP}}$, измеряемой методом зондов с ВЧ-компенсацией в емкостном ВЧ-разряде в аргоне, изучалась в [17]. Изменения в $T_{{exc}}^{{OES}}$, как и в $T_{e}^{{LP}}$, показали тенденцию увеличения доли высокоэнергичных электронов при изменении мощности и давления газа. Температура электронов также определялась в индуктивном и вспомогательном разрядах в неоне с помощью метода (прямого и модифицированного) больцмановских построений [18]. В [19] выполнен сравнительный анализ температур $T_{e}^{{LP}}$ и $T_{{exc}}^{{OES}}$ в двух разных модах (E- и H-модах) для MaPE-ICP-разряда в аргоне. Исследовались зависимости $T_{e}^{{LP}}$, $T_{{exc}}^{{OES}}$ и их отношения $T_{e}^{{LP}}{\text{/}}T_{{exc}}^{{OES}}~$ от приложенной мощности (5–50 Вт) и давления газа (15–60 мТорр). На основе полученных данных проверялась оправданность замены $T_{e}^{{LP}}$ на $T_{{exc}}^{{OES}}$ в обеих модах MaPE-ICP-разряда. В настоящее время различные параметры плазмы исследовались при переходе между модами MaPE-ICP-разряда в неоне [20]. Было показано, что в обеих модах разряда плотность электронов растет, а электронная температура (определенная и зондовым, и спектроскопическим методами) падает с увеличением подводимой мощности и давления газа. Однако мало что известно о связи внешних условий как с электронной температурой, определенной методом зондов и модифицированным методом больцмановских построений, так и с температурой возбуждения, найденной из спектра излучения, для плазмы низкого давления в смесях неона с аргоном в случае обеих мод MaPE-ICP-разряда.

Целью данной работы является сравнительное изучение электронной температуры, определенной методом зондов и модифицированным методом больцмановских кривых, и температуры возбуждения, измеренной методом оптической эмиссионной спектроскопии в двух модах (E- и H-моды) для разрядов низкого давления в смесях инертных газов (Ne–Ar). Поведение этих параметров исследуется при изменении внешних условий, к которым относятся подводимая ВЧ‑мощность и концентрация аргона в его смеси с неоном для MaPE-ICP-разряда. Также изучается переход между E- и H-модами в связи с параметрами плазмы и внешними условиями. В отличие от обычных исследований таких разрядов, проводится анализ зависимостей потенциала плазмы, толщины скин-слоя и ФРЭЭ от ВЧ-мощности и содержания аргона в смеси для обеих мод разряда.

2. ЭКСПЕРИМЕНТАЛЬНАЯ УСТАНОВКА

На рис. 1 представлена система (с запатентованным ферритовым сердечником), используемая для инициирования MaPE-ICP-разряда и описанная в [20]. Высота сделанной из нержавеющей стали цилиндрической камеры составляла 250 мм при внутреннем радиусе 170 мм. Разряд инициировался с помощью подключенного через согласующую цепь ВЧ-генератора (частота 13.56 МГц) в конфигурации с ферритовым сердечником. Исследования проводились в автоматическом режиме, позволяющем минимизировать отражение подаваемого импульса напряжения от нагрузки. Вкладываемая в разряд мощность находилась по разнице между величинами подводимой мощности и отраженной. С помощью турбомолекулярного и роторного насосов реактор откачивался до давления 10–3 Па. Рабочее давление в реакторе контролировалось с помощью Pirani- и Penning-датчиков. Скорость напуска неона и аргона в камеру контролировалась и измерялась посредством Teledyne Hastings измерителя массового потока.

Рис. 1.

Схема источника стимулированной магнитным полюсом индуктивно связанной плазмы (MaPE-ICP) и системы зондовой диагностики.

3. ДИАГНОСТИКА ПЛАЗМЫ

3.1. Зонд Ленгмюра

Как описано в [20], в работе для определения параметров плазмы использовался ВЧ компенсированный одиночный цилиндрический зонд Ленгмюра. Зонд, сделанный из вольфрамовой проволоки диаметром Dp = 0.2 мм с выступающей частью длиной Lp = 8.5 мм, помещался почти в центр камеры. Зонд был соединен с блоком управления и компьютером по аналогии с [21], что позволяло осуществлять сбор данных. В этой системе присутствовал надежный механизм чистки зонда до блеска посредством его бомбардировки электронами. Типичная IV-характеристика одиночного зонда Ленгмюра и ее вторая производная приведены на рис. 2. Потенциал плазмы определяется по положению точки, где вторая производная IV-характеристики обращается в ноль.

Рис. 2.

IV кривая, измеренная одиночным цилиндрическим зондом, и ее вторая производная.

Температура электронов. Вольтамперная характеристика (I–V) и ФРЭЭ определялись зондовым методом. Характеристика I–V интерпретировалась с помощью численного кода, поставляемого вместе с оборудованием. Этот код обрабатывает зондовую I–V-характеристику на основе теории Laframboise [22] для бесстолкновительного слоя. В рамках данной теории, с использованием ряда предположений был предложен метод для определения параметров плазмы по току, собираемому электрически проводящим сферическим или цилиндрическим зондом из бесстолкновительной максвелловской покоящейся плазмы. Зондовая характеристика, связывающая ток на зонд и его потенциал, определялась для обеих геометрий в случае, когда радиус зонда мог быть больше радиуса дебаевского экранирования горячими заряженными частицами вплоть до 100 раз, для всего диапазона отношений ионной температуры к электронной и для отношения потенциала зонда к тепловой энергии горячих частиц в диапазоне от –25 до +25 [22].

ФРЭЭ находилась методом Дравестейна из второй производной кривой I–V по формуле [8]

(1)
${{F}_{p}}(\varepsilon ) = \frac{1}{{\sqrt \varepsilon }}{{f}_{e}}(\varepsilon ) = \frac{{2\sqrt {2{{m}_{e}}} }}{{{{e}^{3}}{{S}_{p}}}}\frac{{{{d}^{2}}I}}{{d{{U}^{2}}}}.$

Здесь $~{{f}_{e}}\left( \varepsilon \right)$ – ФРЭЭ, ${{S}_{p}}$ – площадь поверхности зонда и I – собираемый зондом ток в зависимости от потенциала $~U = \left( {V - {{V}_{p}}} \right)$, где ${{V}_{p}}$ – потенциал плазмы.

Концентрация (${{n}_{e}}$) и температура (${{T}_{e}}$) электронов определялись при интегрировании ФРЭЭ

(2)
${{n}_{e}} = \int\limits_0^\infty {\sqrt \varepsilon ~{{F}_{p}}\left( \varepsilon \right)d\varepsilon } $
и

(3)
${{T}_{e}} = \frac{2}{3}\left\langle \varepsilon \right\rangle ~ = \frac{2}{{3{{n}_{e}}}}\int\limits_0^\infty {{{\varepsilon }^{{3/2}}}{{F}_{p}}\left( \varepsilon \right)d\varepsilon } {\kern 1pt} {\kern 1pt} .$

Здесь $\left\langle \varepsilon \right\rangle $ – средняя кинетическая энергия электронов.

3.2. Оптическая эмиссионная спектроскопия

Спектр излучения разряда в смеси неона с аргоном восстанавливался с помощью программного обеспечения спектрометра “Ocean Spectrometer (HR4000)”, соединенного с компьютером через USB оптоволоконный кабель, как описано в [20]. Показанный на рис. 3 спектр излучения снимался через окно, расположенное вблизи отверстия для зонда, чтобы получать параметры плазмы при изменении внешних условий: ВЧ мощности и содержания аргона в смеси.

Рис. 3.

Спектр излучения разрядной плазмы в смеси Ne–Ar при фиксированной ВЧ-мощности, разных концентрациях аргона и давлениях 1 Па (a), 5 Па (б) и 10 Па (в).

Для оценки температуры возбуждения и электронной температуры с помощью зондового метода и модифицированного метода больцмановских построений были выбраны приведенные в таблице 1 линии Ne–I в диапазоне длин волн 300–800 нм. Этим линиям соответствуют большие разности между энергиями верхних и нижних уровней, что дает наилучшую аппроксимацию прямыми линиями для экспериментальных точек. Необходимая спектроскопическая информация для неона, в частности, – величины ${{A}_{{ij}}}$, ${{g}_{i}}$ и ${{E}_{i}}$ – брались из базы данных [11].

Таблица 1.

Выбранные линии Ne–I для определения $T_{{exc}}^{{OES}}$ и $T_{e}^{{OES}}$ с помощью метода больцмановских построений и его модификации, соответственно

${{\lambda }_{{mn}}}$, нм ${{E}_{m}}$, эВ ${{g}_{m}}$ ${{A}_{{mn}}}$, 107 с–1 $\mathop \sum \nolimits_{m > n} {{A}_{{mn}}}$, 107 с–1
352.05 20.36886 1 0.93 1.86
359.35 20.29728 5 0.099 0.099
363.37 20.25918 1 0.11 0.48
585.25 18.96596 1 6.82 6.91
607.43 18.71138 1 6.03 6.06
609.62 18.70407 5 1.81 5.27
614.31 18.63679 5 2.82 4.98
626.64 18.69336 3 2.49 5.08
640.22 18.55511 7 5.14 5.14
650.65 18.57584 5 3.00 4.90
667.83 18.70407 5 2.33 5.27
703.24 18.38162 3 2.53 3.71

3.2.1. Температура возбуждения. При получении значений температуры возбуждения предполагалось, что верхние энергетические уровни выбранных атомарных переходов находятся в состоянии локального термодинамического равновесия; то есть, плотность заселения этих уровней подчиняется больцмановскому распределению. Температура возбуждения ($kT_{{exc}}^{{OES}}$) измерялась методом больцмановских построений [13] с помощью спектроскопических данных для выбранных излучающих линий Ne–I и соотношения

(4)
$\ln \left( {\frac{{{{I}_{{ij}}}\lambda }}{{{{A}_{{ij}}}{{g}_{i}}}}} \right) = - \frac{{{{E}_{i}}}}{{{{k}_{B}}T_{{exc}}^{{OES}}}} + C.$

Здесь ${{I}_{{ij}}}$ – интенсивность атомарной линии излучения при переходе с i-го уровня на j-й уровень, λ   – длина волны соответствующего перехода, ${{A}_{{ij}}}$ – коэффициент Эйнштейна для спонтанного излучения, ${{g}_{i}}$ и ${{E}_{i}}$ = статистический вес и энергия i-го уровня. Значения $\ln \left( {{{I}_{{ij}}}\lambda {\text{/}}{{A}_{{ij}}}{{g}_{i}}} \right)$ определялись при измерении интенсивности излучения и длины волны. Далее, они строились в зависимости от ${{E}_{i}}$, и значения $kT_{{exc}}^{{OES}}$ вычислялись из наклона прямых, аппроксимирующих экспериментальные точки, как показано на рис. 4a.

Рис. 4.

Больцмановское построение для расчета температуры возбуждения (а) и модифицированное больцмановское построение для расчета температуры электронов (б).

Главным источником ошибок при получении значений $k{{T}_{{exc}}}$ из выражения (4) является неопределенность в знании коэффициентов ${{A}_{{ij}}}$, неточность измерения интенсивности излучения и выбор переходов с недостаточно большой разницей между энергиями верхнего и нижнего уровней. Однако использование логарифмических зависимостей позволяет снизить эти ошибки. Например, ошибка на уровне 20% для отношения ${{I}_{{ij}}}{{\lambda }_{{ij}}}~{\text{/}}{{g}_{i}}{{A}_{{ij}}}$ приводит лишь к ошибке в 5% для логарифма от этой величины [13].

3.2.2. Электронная температура. В большинстве случаев плазма не находится в локальном термодинамическом равновесии, особенно это относится к лабораторной плазме низкого давления. Отклонения от локального термодинамического равновесия связаны с большой разницей в массах электронов и положительных ионов, с сильными градиентами в плазме и с влиянием диффузии. Поскольку локальное термодинамическое равновесие может быть подтверждено только для высоких энергетических уровней, то в этом случае можно говорить о частичном равновесии. Даже в случае этого равновесия значения $k{{T}_{e}}$ могут сильно отличаться от оцененных значений $k{{T}_{{exc}}}$, давая информацию только о низкоэнергичной части ФРЭЭ.

В неравновесной плазме в рамках коронального приближения плотности заселения возбужденных состояний получаются из баланса между их возбуждением электронным ударом из основного состояния атомов (процесс возбуждения) и спонтанным излучением (процесс снятия возбуждения). В результате температура $k{{T}_{e}}$ будет отличаться от температуры возбуждения $k{{T}_{{exc}}}$ и не может быть измерена в рамках простого метода больцмановских построений. Изучаемая плазма не находится в состоянии локального термодинамического равновесия, поскольку она низкотемпературная и находится при низком давлении, когда баланс возбужденных состояний соответствует корональной модели, что проверялось для выделенных линий неона в [23]. Способ получения электронной температуры для неравновесной плазмы в случае коронального баланса основан на модифицированном методе больцмановских построений [13]. Этот метод использовался для оценки электронной температуры ($kT_{e}^{{OES}}$) по снятым спектрам с помощью спектроскопических данных выбранных линий для Ne-I на основе соотношения

(5)
$\ln \left( {\frac{{{{I}_{{ij}}}\mathop \sum \nolimits_{i > j} {{A}_{{ij}}}}}{{h{{f}_{{ij}}}{{A}_{{ij}}}{{b}_{{1i}}}}}} \right) = - \frac{{{{E}_{i}}}}{{kT_{e}^{{OES}}}} + C.$

Здесь ${{I}_{{ij}}}$ – интенсивность линий, ${{A}_{{ij}}}$ – коэффициент Эйнштейна, $\mathop \sum \nolimits_{i > j} {{A}_{{ij}}}$ – сумма коэффициентов Эйнштейна по всем переходам из верхнего состояния, $h{{f}_{{ij}}}$ – разность энергий для i-го и j-го состояний, ${{b}_{{1i}}}$ – коэффициент в экспоненциальной аппроксимации, который можно определить при обработке экспериментальных точек, ${{E}_{i}}$ – энергия i-го уровня и $kT_{e}^{{OES}}$ – температура электронов в энергетических единицах.

Величины $\ln \left( {{{I}_{{ij}}}\mathop \sum \nolimits_{i > j} {{A}_{{ij}}}{\text{/}}h{{f}_{{ij}}}{{A}_{{ij}}}{{b}_{{1i}}}} \right)$ вычислялись по интенсивностям выделенных спектральных линий и спектроскопическим константам и далее представлялись на графиках в виде зависимости от ${{E}_{i}}$. Значения $kT_{e}^{{OES}}$ можно оценить из наклона прямых линий, аппроксимирующих экспериментальные точки, как показано на рис. 4б.

Основные ошибки при получении экспериментальных данных, необходимых в методе больцмановских построений, возникают из-за неопределенности в величинах ${{A}_{{ij}}}$ и интенсивностях излучения, а также из-за неопределенностей при численной подгонке коэффициентов возбуждения ${{k}_{{1i}}}$ с помощью выражения (7) из [13]; то есть, из-за ошибок численного определения ${{b}_{{1i}}}$. Мы определили пределы ошибок в виде ±5% для метода больцмановских построений с учетом неточности при определении выражения под логарифмом в формуле (5). Кроме того, интенсивные линии излучения обычно наблюдаются для переходов с малой (1 эВ и меньше) разницей между энергиями верхнего и нижнего уровней перехода, что приводит к увеличению погрешности при определении температур возбуждения и электронов. Поэтому в данной работе выбирались линии с наибольшей разностью энергий, хотя они были и слабыми. Из-за разброса в экспериментальных точках, обрабатываемых модифицированным методом больцмановских построений, полученным данным по электронной температуре была приписана погрешность ±20% [13].

4. РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

На рис. 5–7 приведены измеренная зондовым методом электронная температура и температуры возбуждения и электронов, измеренные методом оптической эмиссионной спектроскопии, в зависимости от ВЧ-мощности. Кривые приведены для различных давлений и содержания аргона в смесях. При увеличении ВЧ-мощности наблюдается уменьшение всех температур. Снижение электронной температуры можно объяснить тем, что высокоэнергичные электроны теряют свою энергию в неупругих столкновениях с атомами, в результате чего увеличивается скорость ионизации [24]. Также увеличение вкладываемой в разряд ВЧ-мощности приводит к интенсификации ионизационных процессов и, как следствие, к росту концентрации электронов в разряде. При этом увеличивается частота столкновений электронов, что приводит к уменьшению температуры электронов и температуры возбуждения [25]. Из данных на рис. 5–7 следует, что при 1 Па и 65 Вт наблюдается резкий спад температуры (переход между E- и H-модами). Мощность, при которой наблюдается этот переход, снижается с ростом содержания аргона в смеси. Резкий спад температуры при переходе между E- и H-модами можно объяснить резким увеличением концентрации электронов из-за более эффективного поглощения ВЧ-энергии при наличии ферритового сердечника в MaPE-ICP-разряде. Согласно рис. 5–7 температуры, определенные зондовым методом и спектроскопически, имеют одинаковые спадающие тренды с ростом мощности при постоянном давлении. Однако температура, измеренная зондами, оказывается выше температуры, полученной из спектроскопических измерений. Это различие можно приписать тому, что возбуждение атомов осуществляется высокоэнергичными электронами, которые не могут эффективно передавать энергию атомам в неупругих столкновениях, поскольку значительная доля этих электронов быстро диффундирует на стенки, где они рекомбинируют. Кроме того, как следует из критерия Грима, в данном эксперименте концентрация ${{n}_{e}}$$(\sim {\kern 1pt} {{10}^{8}}{\kern 1pt} - {\kern 1pt} {{10}^{{11}}}\;{\text{с}}{{{\text{м}}}^{{ - 3}}})$ меньше критической концентрации, и отклонение от локального термодинамического равновесия также может привести к занижению спектроскопической температуры [26]. Наблюдаемое различие в температурах может также быть связано с тем, что спектроскопическая температура соответствует усредненным по объему значениям $kT_{e}^{{OES}}$, а “зондовая” температура $T_{e}^{{LP}}$ относится к локальным значениям на оси разрядной трубки [27]. Помимо этого, дополнительная неопределенность в определении $kT_{e}^{{OES}}$возникает из-за использования коэффициентов возбуждения для линий излучения. Поэтому оптическая эмиссионная спектроскопия не в состоянии дать электронные температуры, совпадающие с зондовыми измерениями [28].

Рис. 5.

Температура электронов, определенная методом зондов, в зависимости от ВЧ-мощности при разных концентрациях аргона и давлении 1 Па (a), 5 Па (б) и 10 Па (в).

Рис. 6.

Температура электронов, определенная методом оптической эмиссионной спектроскопии, в зависимости от ВЧ-мощности при разных концентрациях аргона и давлении 1 Па (a), 5 Па (б) и 10 Па (в).

Рис. 7.

Температура возбуждения, определенная методом оптической эмиссионной спектроскопии, в зависимости от ВЧ-мощности при разных концентрациях аргона и давлении 1 Па (a), 5 Па (б) и 10 Па (в).

На рис. 8–10 показано изменение $T_{e}^{{LP}}$, $~kT_{e}^{{OES}}$ и $kT_{{exc}}^{{OES}}$ с ростом содержания аргона в смеси при различных значениях ВЧ-мощности и давления газа. Температуры, определенные зондовым методом и спектроскопически, уменьшаются при увеличении содержания аргона. Такое поведение можно объяснить на основе сечений и потенциалов ионизации. Атомы неона имеют больший потенциал ионизации и меньшее сечение ионизации по сравнению с атомами аргона. Поэтому атомы аргона более эффективно ионизуются. В этих условиях электронная температура снижается из-за неупругих столкновений электронов с частицами в плазме [25].

Рис. 8.

Температура электронов, определенная методом зондов, в зависимости от концентрации аргона при разных значениях ВЧ-мощности и давлении 1 Па (a), 5 Па (б) и 10 Па (в).

Рис. 9.

Температура электронов, определенная методом оптической эмиссионной спектроскопии, в зависимости от концентрации аргона при разных значениях ВЧ-мощности и давлении 1 Па (a), 5 Па (б) и 10 Па (в).

Рис. 10.

Температура возбуждения, определенная методом оптической эмиссионной спектроскопии, в зависимости от концентрации аргона при разных значениях ВЧ-мощности и давлении 1 Па (a), 5 Па (б) и 10 Па (в).

Потенциал плазмы ${{V}_{p}}$ – ее важный параметр, который управляет энергией ионов в пристеночном слое. Любые изменения ${{V}_{p}}$ влияют на электрическое поле Е, в котором свободные электроны набирают энергию, необходимую для продолжения ионизационных процессов в плазме. На рис. 11 показан потенциал плазмы в зависимости от подводимой ВЧ-мощности при различных значениях давления и содержания аргона в разряде. Сначала, когда разряд существует в Е-моде, потенциал плазмы растет с увеличением мощности. Но потом, после перехода в Н-моду, потенциал снижается. Такое поведение потенциала можно объяснить следующим образом. В Е-моде ${{V}_{p}}$ увеличивается, что свидетельствует о малом количестве электронов в объеме плазмы из-за низкой скорости ионизации при малой мощности и эффективной диффузии высокоэнергичных электронов на стенки реактора. Однако с ростом ВЧ-мощности скорость ионизации и концентрация электронов увеличиваются в объеме плазмы, что приводит к падению потенциала. Другой причиной уменьшения потенциала плазмы может быть эффективное удержание плазмы в объеме при растущей напряженности создаваемого магнитным полюсом магнитного поля [29]. Из рис. 11 также следует, что потенциал плазмы уменьшается с увеличением содержания аргона в плазме. Это может быть объяснено ростом концентрации электронов при увеличении доли аргона в смеси.

Рис. 11.

Потенциал плазмы в зависимости от ВЧ-мощности при разных концентрациях аргона и давлении 1 Па (a), 5 Па (б) и 10 Па (в).

На рис. 12 приведена толщина скин-слоя в зависимости от ВЧ-мощности при различных давлениях и концентрациях аргона. Концентрация плазмы и интенсивность испускаемого излучения существенно меняются при переходе между E- и H-модами. В E-моде концентрация очень низкая (ne ~ 109 см–3), в то время как в H-моде она значительно выше (ne ~ 1010–1011 см–3). В [30] также был приведен аналитический критерий устойчивости Н-моды. Согласно этому критерию, чтобы данная мода была устойчива в разряде низкого давления, толщина скин-слоя $({{\delta }_{p}})$ должна быть меньше критического размера $\sqrt {2{\text{/}}3} R$, где R – радиус реактора. Толщина скин-слоя в бесстолкновительной плазме вычисляется согласно

(6)
${{\delta }_{p}} = {{\left( {\frac{{{{m}_{e}}}}{{{{n}_{e}}{{e}^{2}}{{\mu }_{o}}}}} \right)}^{{1/2}}}.$
Рис. 12.

Толщина скин-слоя в зависимости от ВЧ-мощности при разных концентрациях аргона и давлении 1 Па (a), 5 Па (б) и 10 Па (в).

Здесь $~{{m}_{e}} = 9.11 \times {{10}^{{ - 31}}}$ кг – масса электрона, e = $ = 1.19 \times {{10}^{{ - 19}}}$ К – абсолютная величина его заряда, и µ0 = 1.26 × 10–6 м кг с–2 А–2 – магнитная проницаемость вакуума. При подстановке этих величин и измеренной концентрации электронов $\left( {{{n}_{e}}} \right)$ в (6) можно вычислить толщину скин-слоя. В рассматриваемом случае критический размер равен 13.8 см. Как видно на рис. 12, получаемая толщина скин-слоя в H-моде меньше 13.8 см при внешней ВЧ-мощности 65–80 Вт и 0% аргона, а также при 20–60 Вт и 20–100% аргона. Это является следствием того, что толщина скин-слоя резко падает при переходе разряда в Н-моду, поскольку ${{n}_{e}}$ увеличивается на один-два порядка величины по сравнению с концентрацией электронов в Е‑моде.

На рис. 13–15 представлены измеренные ФРЭЭ в области перехода между модами для разных концентраций аргона в смеси и давлений газа (1, 5 и 10 Па). В E-моде ФРЭЭ является немаксвелловской, а в Н-моде эти функции близки к максвелловской. На рис. 13а виден бимаксвелловский характер ФРЭЭ в Е-моде, что является отличительной чертой нелокальной кинетики электронов в плазме при низких давлениях [31]. В разрядах низкого давления, где частота столкновений электронов с нейтральными частицами $({{\nu }_{{e - n}}})$ значительно меньше частоты ВЧ-источника $({{\omega }_{{RF}}})$, обычно преобладает бесстолкновительный нагрев электронов при их взаимодействии со скин-слоем вблизи рамки антенны. Энергичные электроны могут преодолевать барьер амбиполярного потенциала и достигать скин-слоя. В результате энергичные электроны приобретают дополнительную энергию, что приводит к ускорению ионизации и возбуждения нейтральных частиц при неупругих столкновениях с такими электронами. В то же время электроны с низкой энергией не способны преодолеть амбиполярный потенциальный барьер и достичь скин-слоя, где происходит бесстолкновительный нагрев [32]. На рис. 14a и 15a показано, что ФРЭЭ превращается из бимаксвелловской в дравестейновскую в режиме Е-моды разряда. Переход к ФРЭЭ типа дравестейновских при высоких давлениях (5 и 10 Па) связан с тем, что механизм нагрева электронов вместо бесстолкновительного (стохастического) становится столкновительным (омическим). Дравестейновский характер ФРЭЭ формируется прежде всего благодаря эффективному нагреву низкоэнергичных электронов в столкновительном режиме [33].

Рис. 13.

ФРЭЭ в E- и H-модах при разных концентрациях аргона, фиксированном давлении (1 Па) и ВЧ-мощности 10 Вт (a) и 65 Вт (б).

Рис. 14.

ФРЭЭ в E- и H-модах при разных концентрациях аргона, фиксированном давлении (5 Па) и ВЧ -ощности 10 Вт (a) и 65 Вт (б).

Рис. 15.

ФРЭЭ в E- и H-модах при разных концентрациях аргона, фиксированном давлении (10 Па) и ВЧ-мощности 10 Вт (a) и 65 Вт (б).

Как показано на рис. 13б–15б, в Н-моде ФРЭЭ являются максвелловскими для всех исследованных давлений. Это можно приписать возросшей концентрации электронов в разряде при увеличении ВЧ-мощности и концентрации аргона, в результате чего выросла частота электрон-электронных (e–e) столкновений, что привело к максвеллизации ФРЭЭ. Однако с ростом содержания аргона в смеси происходит обеднение высокоэнергичного хвоста ФРЭЭ. Этот эффект можно объяснить следующим образом. Увеличение доли аргона приводит к большей длине свободного пробега для низкоэнергичных электронов из-за эффекта Рамзауэра и к росту концентрации высокоэнергичных электронов. Потенциал ионизации аргона (15.7 эВ) больше аналогичной величины для неона (21.5 эВ). Поэтому присутствующие в плазме высокоэнергичные электроны участвуют в неупругих столкновениях с атомами аргона, ионизуя их. В результате истощается высокоэнергичный хвост ФРЭЭ, и увеличивается концентрация низкоэнергичных электронов.

5. ВЫВОДЫ

В данной работе исследуются параметры плазмы в смесях неона с аргоном для двух различных мод разряда. В частности, изучается корреляция между электронной температурой, полученной из зондовых измерений, и температурой электронов, и температурой возбуждения, определенных с помощью оптической эмиссионной спектроскопии. Все измеренные температуры демонстрируют аналогичные спадающие тренды при увеличении ВЧ-мощности и содержания аргона в смесях. Это свидетельствует о том, что спектроскопическое определение температуры может быть хорошей альтернативной измерению температуры зондовым методом. В то же время из измерений следует, что спектроскопическая температура меньше температуры, определенной с помощью зондов, ($kT_{e}^{{OES}}$, $kT_{{exc}}^{{OES}}$< $T_{e}^{{LP}}$); это различие приписывается отклонению состояния плазмы от локального термодинамически равновесного. По аналогии с интенсивностью излучения, концентрация электронов, толщина скин-слоя и температуры, определенные спектроскопически и методом зондов, испытывают резкий скачок при переходе между Е- и Н-модами. Это указывает на то, что температуры $T_{e}^{{LP}}$, $kT_{e}^{{OES}}$ и $kT_{{exc}}^{{OES}}$ могут быть хорошими индикаторами перехода между Е‑ и Н-модами в плазме в смеси неона с аргоном. Также анализируется поведение потенциала плазмы при переходе между модами. Из измерений следует, что потенциал растет в Е-моде, но показывает тенденцию к спаду при переходе к Н‑моде. При этом выполняется полученный в [21] аналитический критерий стабильности Н‑моды, поскольку вычисленная толщина скин-слоя меньше критического размера разрядного реактора (13.8 см) для ВЧ-мощности в диапазоне 20–60 Вт (H-мода). При низких ВЧ-мощностях, когда разряд существует в Е-моде, измеренные Ф-РЭЭ оказываются немаксвелловскими. С ростом ВЧ-мощности наблюдается переход между Е- и Н-модами, при котором происходит скачкообразное снижение $T_{e}^{{LP}}$, $kT_{e}^{{OES}}$, $kT_{{exc}}^{{OES}}$ и толщины скин-слоя, а ФРЭЭ становится максвелловской.

Список литературы

  1. Grill A. Cold plasma in materials fabrication. IEEE Press, New York, 1994.

  2. Park H., Choe W. Current Applied Physics. 2010. V. 10. P. 1456.

  3. Moon S.Y., Choe W., Uhm H.S., Hwang Y., Choi J. Phys. Plasmas. 2002. V. 9. P. 4045.

  4. Jovićević S., Ivković M., Konjević N. Spectrochimica Acta Part B: Atomic Spectroscopy. 2001. V. 56. P. 2419.

  5. Onga C., Chin O., Nisoab M., Paosawatyanyangc B. Jurnal Fizik Malaysia. 2004. V. 25.

  6. Godyak V., Piejak R., Alexandrovich B. Plasma sources sci. technology. 1992. V. 1. P. 36.

  7. Kang N., Oh S.-G., Ricard A. J. Phys. D: Applied Phys. 2008. V. 41. P. 155203.

  8. Godyak V., Demidov V. J. Phys. D: Applied Phys. 2011. V. 44. P. 233001.

  9. Fuller N., Donnelly V.M., Herman I.P. J. Vacuum Science & Technology A: Vacuum, Surfaces, and Films. 2002. V. 20. P. 170.

  10. Zhu X.-M., Chen W.-C., Li J., Pu Y.-K. J. Phys. D: A-pplied Phys. 2008. V. 42. P. 025203.

  11. https://www.nist.gov/.

  12. Griem H.R. Principles of plasma spectroscopy. Cambridge University Press, 2005.

  13. Gordillo-Vázquez F., Camero M., Gomez-Aleixandre C. Plasma Sources Sci. Technology. 2005. V. 15. P. 42.

  14. Hopwood J. Plasma Sources Sci. Technology. 1992. V. 1. P. 109.

  15. Singh S. J. Applied Phys. 2008. V. 103. P. 083303.

  16. Park I.-S., Kang H.-J., Kim K.-H., Chung C.-W. Phys. Plasmas. 2017. V. 24. P. 123506.

  17. Park H., You S., Choe W. Phys. Plasmas. 2010. V. 17. P. 103501.

  18. Tanisli M., Sahin N., Demir S. arXiv preprint arXiv:1706.05015. 2017.

  19. Khan A., Jan F., Saeed A., Zaka-ul-Islam M., Abrar M., Khattak N., Zakaullah M. Current Appl. Phys. 2013. V. 13. P. 1241.

  20. Khattak Z.I., Shafiq M., Khan A.W. IEEE Transactions Plasma Sci. 2019. V. 47. P. 2665.

  21. Lee M.-H., Jang S.-H., Chung C.-W. J. Appl. Phys. 2007. V. 101. P. 033305.

  22. Laframboise J.G. Theory of spherical and cylindrical Langmuir probes in a collisionless, Maxwellian plasma at rest. Toronto Univ Downsview (Ontario) Inst. for Aerospace Studies, 1966.

  23. Rehman N., Zakaullah M., Khan F., Naseer S. J. Appl. Phys. 2008. V. 104. P. 123304.

  24. Jan F., Khan A., Saeed A., Zakaullah M. Contributions Plasma Phys. 2013. V. 53. P. 492.

  25. Orazbayev S., Ussenov Y.A., Ramazanov T., Dosbola-yev M., Utegenov A. Contributions Plasma Phys. 2015. V. 55. P. 428.

  26. Yanguas-Gil A., Cotrino J., González-Elipe A.R. J. Appl. Phys. 2006. V. 99. P. 033104.

  27. Kano K., Suzuki M., Akatsuka H. Plasma Sources Sci. Technology. 2000. V. 9. P. 314.

  28. Tanışlı M., Şahin N., Younus M., Rehman N.U., Demir S. Phys. Plasmas. 2017. V. 24. P. 102123.

  29. Younus M., Rehman N., Shafiq M., Naeem M., Zaka-ul-Islam M., Zakaullah M. Plasma Phys. Rep. 2017. V. 43. P. 588.

  30. Lee M.-H., Chung C.-W. Phys. Plasmas. 2006. V. 13. P. 063510.

  31. Lee H.-C., Lee M.-H., Chung C.-W. Appl. Phys. Lett. 2010. V. 96. P. 041503.

  32. Lee H.-C., Chung C.-W. Phys. Plasmas. 2012. V. 19. P. 043505.

  33. Chung C., Chang H.-Y. Appl. Phys. Lett. 2002. V. 80. P. 1725.

Дополнительные материалы отсутствуют.