Parameters and Composition of Plasma in a Mixture of CF4 + H2 + Ar: Effect of the CF4/H2 Ratio
- Autores: Miakonkikh A.V.1, Kuzmenko V.O.1, Efremov A.M.2, Rudenko K.V.1
-
Afiliações:
- Valiev Institute of Physics and Technology of the Russian Academy of Sciences
- Ivanovo State University of Chemistry and Technology
- Edição: Volume 53, Nº 1 (2024)
- Páginas: 91-101
- Seção: ТЕХНОЛОГИИ
- URL: https://ogarev-online.ru/0544-1269/article/view/259608
- DOI: https://doi.org/10.31857/S0544126924010104
- ID: 259608
Citar
Texto integral
Resumo
A study of the electrophysical parameters of plasma and the kinetics of plasma-chemical processes in CF4 + H2 + Ar mixture by varying the CF4/H2 ratio were carried out. Using diagnostic methods and plasma modeling together, it was found that replacing tetrafluoromethane with hydrogen a) leads to a decrease in plasma density and an increase in electronegativity; and b) causes a disproportionately sharp drop in the concentration of fluorine atoms. The reason for the latter effect is an increase in the frequency of death of atoms in reactions of the form CHFx + F → CFx + HF, initiated by heterogeneous recombination according to the CFx + H → CHFx mechanism. A simultaneous increase in the concentration of polymer-forming radicals CHxFy (x + y < 3) indicates an increase in the polymerization load of the plasma on the surfaces in contact with it.
Palavras-chave
Texto integral
ВВЕДЕНИЕ
Фторуглеродные газы семейства CxHyFz нашли широкое применение в технологии интегральной микроэлектроники для структурирования поверхностей кремния и его соединений [1, 2]. Основным инструментом здесь являются процессы реактивно-ионного травления (РИТ), реализуемые в индукционных разрядах пониженного давления и сочетающие физические (ионное распыление) и химические (гетерогенные реакции нейтральных частиц) эффекты взаимодействия плазмы с обрабатываемой поверхностью [3]. Известно также, что выходные характеристики РИТ (скорость, анизотропия, селективность по отношению к материалу маски и нижележащего слоя) в значительной степени зависят от полимеризационной способности плазмообразующего газа, определяющей скорость осаждения и толщину фторуглеродной полимерной пленки на обрабатываемой поверхности [2, 3]. В частности, толстая сплошная полимерная пленка снижает скорость травления кремния (из-за уменьшения потока атомов фтора на границе полимер/поверхность), но увеличивает анизотропию процесса и селективность в системе SiO2/Si [3–6]. Последние два эффекта обусловлены ингибированием химических реакций на боковых стенках формируемого рельефа и меньшей толщиной пленки на кислородсодержащей поверхности [5, 6]. Таким образом, оптимизация выходных характеристик РИТ под задачи конкретного процесса требует целенаправленного регулирования концентраций различных типов активных частиц плазмы, определяющих кинетику травления и полимеризации.
Отличительной особенностью современной технологии РИТ является использование плазмообразующих смесей, комбинирующих фторуглеродный газ с одним или несколькими дополнительными компонентами. Так, по результатам предшествующих исследований [7–12] было установлено, что добавка O2 приводит к росту концентрации атомов фтора, снижает концентрацию полимеробразующих радикалов, а также вызывает окислительную деструкцию уже осажденной фторуглеродной полимерной пленки. Фактически, это обусловливает возможность применения газов с высокой полимеризационной способностью (например, C4F8 или CHF3) в процессах, требующих высоких скоростей травления и минимального загрязнения поверхности фторуглеродным полимером. Логично предположить, что расширение области применения возможно и для газов с низкой полимеризационной способностью (например, для CF4) за счет введения в плазмообразующую смесь компонента, стимулирующего полимеризацию. В частности, в работах по смеси CF4 + H2 [13–15] было показано, что увеличение доли водорода приводит к снижению концентрации атомов фтора, увеличивает концентрацию фторуглеродных радикалов и молекул HF, а также влияет на кинетику травления Si и SiO2 через изменение толщины и химического состава осаждаемой полимерной пленки. Хотя на качественном уровне многое объясняется протеканием процессов вида CFx + H → CHFx и CHFx + F → CFx + HF, полноценный кинетический анализ до сих пор отсутствует. Кроме того, явными недостатками модельных работ [14, 15] являются постулирование постоянной концентрации электронов (что в принципе не согласуется с существенными изменениями состава плазмы) и отсутствие сравнения расчетных концентраций частиц с экспериментальными данными. Такая ситуация не позволяет сделать однозначного вывода об адекватности кинетических схем, предлагаемых авторами для объяснения наблюдаемых эффектов.
В последнее время актуальность исследования полимеризации в плазме определяется интересом к процессу атомно-слоевого травления – перспективному методу прецизионного плазмохимического травления, обладающему достаточной точностью и селективностью для формирования перспективных приборов наноэлектроники [16]. Один из наиболее распространенных подходов к реализации процесса состоит в последовательном циклическом формировании тонкого модифицированного слоя на поверхности материала, например за счет осаждения тонких фторуглеродных пленок из плазмы, и последующей активации реакции между материалом, подвергаемом травлению, и фторуглеродной пленкой за счет бомбардировки ионами из плазмы химически нейтрального газа. Подходы к реализации атомно-слоевого травления с использованием осаждения фторуглеродных пленок были исследованы для полимерообразующих плазм C4F8, CHF3 [17], CF3I [18]. К процессу при этом предъявляются противоречивые требования: отсутствие спонтанного травления и при этом формирование достаточно тонкой пленки, реакция которой с поверхности может быть активирована ионной бомбардировкой без существенного распыления атомов поверхности. По этой причине исследование скорости формирования полимерной пленки и ее элементного состава становится важным условием разработки процессов атомно-слоевого травления. Идея прецизионного контроля пленкообразования за счет выбора соотношения потоков газов плазмообразующей смеси Ar/CF4/H2 для целей атомно-слоевого травления была предложена в работе [19]. Было установлено, что, изменяя состав смеси Ar/CF4/H2, можно регулировать баланс конкурирующих процессов травления и пленкообразования, а равновесие определяется составом поверхности, подвергающейся травлению. Предложенный подход позволил оптимизировать селективность травления Al2O3, HfO2, AlN по отношению к другим материалам.
В наших предшествующих работах было проведено детальное изучение физико-химических свойств плазмы в смесях CF4 + Ar [20, 21], CF4 + O2 [11, 12] и CF4 + Ar + O2 [11, 22] в условиях, типичных для РИТ кремния и его соединений. Совместное использование методов диагностики и моделирования плазмы позволило: а) установить взаимосвязи между начальным составом смесей и кинетикой плазмохимических процессов, определяющих стационарные концентрации активных частиц; б) провести верификацию кинетических схем, предпосылок и алгоритмов моделирования при сравнении расчетных и измеренных концентраций атомов фтора. Идея данной работы заключается в использовании аналогичного подхода для исследования смеси CF4 + H2 + Ar при варьировании соотношения CF4/H2. Соответственно, основное внимание было направлено: а) на формирование кинетической схемы, обеспечивающей адекватное описание кинетики плазмохимических процессов в присутствии водорода; б) выявление влияния добавок водорода на параметры электронного и ионного компонентов плазмы; в) идентификацию механизмов, определяющих влияние водорода на химическую активность и полимеризационную способность плазмы.
МЕТОДИЧЕСКАЯ ЧАСТЬ
2.1. Оборудование и условия эксперимента
Эксперименты проводились в условиях индукционного ВЧ (2.0 МГц) разряда, в реакторе PlasmaLab 100 (OIPT, UK). В качестве постоянных внешних параметров плазмы выступали скорость потока газа (q = 60 ст. см3/мин), давление газа (p = 10 мтор) и вкладываемая мощность (W = 1500 Вт), что соответствовало плотности мощности примерно 0.03 Вт/см3). Варьируемым параметром являлся начальный состав плазмообразующей смеси CF4 + H2 + Ar, задаваемый соотношением индивидуальных расходов ее компонентов. В частности, неизменное значение qAr = 50 ст. см3/мин обусловливало постоянную долю аргона yAr = qAr / q на уровне 83%, при этом варьирование qH2 в диапазоне 2–8 ст. см3/мин соответствовало yH2 = 3–13% за счет пропорционального снижения содержания тетрафторметана.
2.2. Диагностика плазмы
Данные по электрофизическим параметрам плазмы получали с помощью одиночного зонда Ленгмюра (Espion, Hiden Analytical, UK). Для минимизации искажений зондовых вольт-амперных характеристик (ВАХ) из-за пленкообразования на зондах проводилась их очистка в плазме 50% Ar + 50% O2 в течение примерно 2 мин перед и после каждого измерения. Эффективность данной процедуры подтверждается удовлетворительной воспроизводимостью результатов в серии последовательных измерений при неизменных условиях возбуждения плазмы. Неизменность размеров зонда до и после проведения диагностики свидетельствует об отсутствии травления его материала активными частицами плазмы. Обработка измеренных ВАХ базировалась на известных положениях зондовой теории для разрядов низкого давления [23]. Результатом обработки являлись данные по температуре электронов (Te), суммарной концентрации положительных ионов (n+) и концентрации электронов (ne).
Спектры излучения измеряли с помощью Ocean Insight HR4Pro в диапазоне длин волн 200–900 нм с оптическим разрешением 0.7 нм. При определении интенсивностей линий была учтена спектральная функция пропускания волокна. Стационарные концентрации атомов фтора и водорода определяли методом оптической актинометрии. В первом случае использовалась аналитическая пара F 703.8 нм/Ar 750.4 нм, в которой для обоих возбужденных состояний доминирует заселение прямым электронным ударом [24, 21]. Расчеты концентраций атомов фтора проводили на основе измеренных интенсивностей излучения IF и IAr по соотношению
, (1)
где nAr = yArN – концентрация аргона, N = p / kBTgas – общая концентрация частиц в реакторе при температуре газа Tgas; – актинометрический коэффициент, определяемый отношением соответствующих длин волн λ, констант возбуждения ke и вероятностей оптических переходов [25].
Сечения возбуждения и параметры оптических переходов для атомов фтора и аргона хорошо известны по литературным данным [25–27]. В работе [25] было показано также, что определенные таким образом концентрации атомов фтора удовлетворительно согласуются с результатами масс-спектральных измерений.
Что касается атомов водорода, использование аналогичного подхода может быть осложнено диссоциативным возбуждением молекул H2 + e → H + H* + e (kde, H2), которое составляет конкуренцию прямому процессу H + e → H* + e (ke, H) [28–30]. Отношение измеренных интенсивностей излучения атомов водорода и аргона в данном случае формально может быть записано как
, (2)
при этом переход к виду уравнения (1) требует . Расчеты констант скоростей c использованием известных сечений соответствующих процессов [31] показали, что kde, H2 << ke, HnH (например, kde, H2 ~ 5.3 × 10–13 см3/с и ke, H ~ 4.1 × 10–11 см3/с при Te = 3 эВ для линии H 656.4 нм) в силу высокой пороговой энергии диссоциативного возбуждения. Поэтому заметный вклад диссоциативного возбуждения в экспериментально определяемую величину IH можно ожидать лишь при малых (менее 5%) степенях диссоциации молекул H2, обеспечивающих nH / nH2 < 0.01. В то же время из имеющихся работ [29, 30] можно заключить, что характерные значения степени диссоциации H2 в индукционных разрядах низкого давления при сопоставимых с нашими значениями плотности вкладываемой мощности и концентрации электронов превышают 10%, что соответствует nH / nH2 > 0.2. Таким образом, определение концентрации атомов водорода соотношению вида уравнения (1) в исследуемом нами диапазоне условий является правомерным. Из данных работы [28] можно заключить также, что концентрации атомов водорода, определенные с помощью актинометрических пар Н 486.3 нм/Ar 750.4 нм и H 656.4 нм/Ar 750.4 нм, неодинаково согласуются с результатами независимых измерений. В частности, в первом случае лучшие результаты наблюдаются для области очень низких степеней диссоциации H2, а во втором – при степенях диссоциации более 5–8%. Поэтому для целей нашего исследования была выбрана пара H 656.4 нм/Ar 750.4 нм. Во всех расчетах использовалось значение Tgas = 600 К. Данная величина отражает результаты экспериментальных измерений данного параметра в плазме Ar [32], CF4 [33] и H2 [29], приведенные к актуальным значениям давления, плотности мощности и соотношения компонентов в плазмообразующей смеси.
2.3. Моделирование плазмы
Расчеты состава плазмы CF4 + H2 + Ar проводились на основе кинетической схемы, которая ранее неоднократно использовалась для моделирования плазмы CHF3 + Ar (табл. 1) [20, 34, 35]. Допустимость такого подхода обусловлена тем, что газовая фаза в обоих случаях представлена идентичным набором частиц. Отметим, что особенностями данной схемы являются: а) вынужденное игнорирование фторуглеродных компонентов вида СHxFy с x > 1 (из-за высокой неопределенности данных по механизмам и сечениям их диссоциации в плазме); б) связанное с этим отсутствие реакций частиц семейства метана, инициируемых образованием радикалов CH по R11. Тем не менее даже в рамках таких упрощений имеет место удовлетворительное согласие расчетных параметров и состава плазмы CHF3 с данными независимых экспериментов [35, 36]. Очевидным объяснением этого факта служат незначительные концентрации компонентов вида СHxFy с x > 1 и CHx, обуславливающие их пренебрежимо малое влияние на энергетическое распределение электронов и кинетику плазмохимических процессов. Что касается плазмы CF4 + H2 + Ar, то по первому пункту ситуация представляется еще более благоприятной. Причиной этого служит большее, по сравнению с плазмой CHF3, число промежуточных процессов, необходимых для образования “проблемных” частиц. Для оценки принципиальности второго пункта кинетическая схема была дополнена набором реакций для CHx (x = 1–4) из работ по моделированию плазмы метана [37–39]. Предварительные расчеты показали, что учет этих процессов не приводит к значимым (т. е. превышающим величины погрешностей входных параметров модели) изменениям концентраций основных типов активных частиц, определяющих физические и химические эффекты взаимодействия плазмы с поверхностью кремния и его соединений. Поэтому в табл. 1 соответствующие реакции не представлены.
Таблица 1. Основные процессы, определяющие состав плазмы в смеси CF4 + H2 + Ar
№ п/п | Реакция | k | № п/п | Реакция | k |
1 | CFx + e → CFx–1 + F + e | f (Te) | 21 | F2 + CFx → CFx+1 + F | f (Tgas) |
2 | CFx + e → CFx–1+ + F + 2e | f (Te) | 22 | CFx + F → CFx+1 | f (Tgas) |
3 | CFx + e → CFx–1 + F– | f (Te) | 23 | CHFx + F → CHFx+1 | f (Tgas) |
4 | F2 + e → 2F + e | f (Te) | 24 | F2 + H → HF + F | f (Tgas) |
5 | F2 + e → F2+ + 2e | f (Te) | 25 | H2 + F → HF + H | f (Tgas) |
6 | F2 + e → F + F– | f (Te) | 26 | CFx + H → CFx-1 + HF | f (Tgas) |
7 | F + e → F+ + 2e | f (Te) | 27 | CHFx + F → CFx + HF | f (Tgas) |
8 | H2 + e → 2H + e | f (Te) | 28 | CHFx + H → CHFx–1 + HF | f (Tgas) |
9 | H2 + e → H2+ + 2e | f (Te) | 29 | F– + X+ → F + X | — |
10 | H + e → H+ + 2e | f (Te) | 30 | CFx → CFx (s) | f (g) |
11 | CHFx + e → CHFx–1 + F + e | f (Te) | CFx(s) + F → CFx+1 | ||
12 | CHFx + e → CHFx–1+ + F + 2e | f (Te) | CFx (s) + H → CHFx | ||
13 | CHFx + e → CFx + H + e | f (Te) | 31 | F → F (s) | f (g) |
14 | CHFx + e → CFx+ + H + 2e | f (Te) | F (s) + CFx → CFx+1 | ||
15 | CHFx + e → CFx–1 + HF + e | f (Te) | F (s) + F → F2 | ||
16 | CHFx + e → CHFx–1 + F– | f (Te) | F (s) + H → HF | ||
17 | HF + e → H + F + e | f (Te) | 32 | H → H (s) | f (g) |
18 | HF + e → HF+ + 2e | f (Te) | H (s) + CFx → CHFx | ||
19 | HF + e → H + F– | f (Te) | H (s) + H → H2 | ||
20 | Ar + e → Ar+ + 2e | f (Te) | H (s) + F → HF |
Примечание: X+ в R29 — любой положительный ион.
Алгоритм моделирования плазмы базировался на системе типовых допущений, используемых при описании физики и химии индукционных ВЧ-разрядов низкого (p < 20 мтор) давления [3, 21, 22]. В качестве основных предпосылок модели с учетом специфики исследуемой системы выступали следующие.
Высокая степень ионизации газа (n+ / N > 10–4) определяет существенный вклад равновесных электрон-электронных соударений в формирование функции распределения электронов по энергиям (ФРЭЭ). Такая ситуация допускает использование максвелловской ФРЭЭ для описания кинетики процессов под действием электронного удара (R1–R20), при этом константы скоростей могут быть найдены по параметрическим соотношениям вида [30, 35].
- Уравнение квазинейтральности электроотрицательной плазмы имеет вид ne + n– = n+, где – концентрация отрицательных ионов, определяемая балансом скоростей диссоциативного прилипания и ион-ионной рекомбинации. Таким образом, концентрация электронов связана с концентраций положительных ионов соотношением
, (3)
где k29 ≈ 10–7 см3/с [8, 35] – эффективная (принимаемая условно равной для всех типов положительных ионов) константа скорости ион-ионной рекомбинации; – суммарная частота диссоциативного прилипания, формируемая соответствующими процессами для CFx (R2), F2 (R6), CHFx (R16) и HF (R19).
Константы скоростей для R3, R6 и R16 известны по литературным данным в виде параметрических функций температуры электронов [8, 30, 35]. Константа скорости R19 доступна только по результатам единичных измерений и/или оценок [26, 27, 40, 41], при этом единственным (т. е. не обеспечивающим возможности сравнения с аналогичными данными других авторов) источником сечения этого процесса является расчетная работа [42]. Тем не менее значение k19 ~ 5 × 10–11 см3/с при Te = 2–4 эВ, полученное при интегрировании данного сечения с максвелловской ФРЭЭ, адекватно отражает уменьшение эффективности прилипания электронов к молекулам галогенводородов в ряду HI–HBr–HCl–HF [40, 41].
- Гетерогенная рекомбинация атомов и радикалов происходит по первому кинетическому порядку относительно их концентрации в газовой фазе (механизм Или–Ридила). Справедливость данной предпосылки обусловлена тем, что стенки металлической разрядной камеры в условиях низких давлений газа (высоких длин свободного пробега частиц) характеризуются высокими степенями насыщения адсорбированными частицами. Соответственно, константа скорости рекомбинации для частицы молярной массы M и характеризующейся вероятностью рекомбинации [30, 34–38] может быть найдена как k ≈ (r + l) υTγ / 2rl [43], где – скорость хаотического теплового движения.
Выходными параметрами модели служили скорости процессов образования и гибели частиц, их концентрации и плотности потоков на обрабатываемую поверхность.
РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
На рис. 1 представлены экспериментальные и расчетные данные по влиянию содержания водорода на параметры электронного и ионного компонентов плазмы. Анализ этих результатов позволяет выделить ряд значимых свойств исследуемой системы, определяющих кинетику плазмохимических процессов и стационарный состав плазмы:
Рис. 1. Электрофизические параметры плазмы в смеси CF4 + H2 + Ar при p = 10 мтор и W = 1500 Вт: 1 — температура электронов; 2 — концентрация электронов; 3 — суммарная концентрация положительных ионов; 4 — относительная концентрация отрицательных ионов. Пунктирные линии представляют значения, полученные при моделировании плазмы.
- Температура электронов (рис. 1, а) практически не зависит от содержания H2 в плазмообразующем газе, демонстрируя слабую (в пределах погрешности эксперимента) тенденцию к росту в смесях, богатых водородом. Очевидно, что причиной этого является неизменность доминирующих каналов потерь энергии электронами и незначительная деформация ФРЭЭ в условиях, когда основными компонентами газовой фазы являются молекулярные частицы (рис. 2, а).
- Концентрация электронов (см. рис. 1, а) монотонно снижается с ростом доли водорода в смеси (в ~ 1.4 раза при 3–13% H2) на фоне практически постоянной суммарной частоты ионизации (1.2 × 104–1.3 × 104 с–1 при 3–13% H2). Последнее обусловлено тем, что основной вклад в νiz вносит реакция R20, протекающая при неизменной концентрации исходного компонента. Таким образом, определяющим фактором для зависимости является рост частоты гибели электронов в процессах диссоциативного прилипания νda (1.4 × 103–1.8 × 103 с–1 при 3–13% H2) по причине k19 > k3 и nHF > при > 5–6% (см. рис. 2, а). Хорошее согласие экспериментальных и расчетных значений свидетельствует о корректном описании кинетики этих частиц в рамках предложенной кинетической схемы.
- Концентрация положительных ионов (рис. 1, б) также снижается с ростом доли водорода в смеси, следуя изменению суммарной скорости ионизации νizne. Чуть более медленная тенденция по сравнению с аналогичной зависимостью концентрации электронов обусловлена увеличением абсолютной (2.4 × 1010–3.2 × 1010 см–3 при 3–13% H2) и относительной (0.14–0.25 при 3–13% H2) концентрации отрицательных ионов по причине роста νda.
Рис. 2. Концентрации нейтральных частиц (а) и скорости образования атомов фтора (б) в плазме смеси CF4 + H2 + Ar. Пунктирными линиями на рис. а показаны, которые выделяют водородсодержащие частицы. Условия возбуждения плазмы соответствуют рис. 1.
Из представленных выше данных можно заключить, что особенностью рассматриваемой системы является рост электроотрицательности и снижение плотности плазмы при увеличении содержания электроположительного компонента в плазмообразующей смеси. Это формальное противоречие обусловлено изменениями кинетики нейтральных частиц, приводящими к образованию и доминированию более электроотрицательных продуктов плазмохимических реакций. Рассмотрим соответствующие вопросы более детально.
Анализ кинетики нейтральных частиц показал, что основные свойства плазмы CF4 + H2 + Ar в условиях < 3% аналогичны свойствам бинарной смеси CF4 + Ar, известным по результатам предшествующих работ [20–22]. В частности, основными фторсодержащими компонентами газовой фазы являются исходные молекулы CF4, радикалы CF3 и атомы F, при этом выполняются условия nHF > и (см. рис. 2, а). Такая ситуация обеспечивается: а) ступенчатым характером образования фторуглеродных радикалов по R1, каждый элементарный акт которой приводит к высвобождению атома фтора; б) более высокими значениями вероятностей и констант скоростей гетерогенной рекомбинации для менее насыщенных радикалов (k30 ~ 125 с–1 для CF3 и ~ 315 с–1 для CF). Основной (около 80%) вклад в суммарную скорость образования атомов фтора вносит диссоциация CFx (x = 3, 4) по механизмам R1 и R2, при этом примерно 10% обеспечивается реакцией R4 (рис. 2, б). Заметная роль R4 обусловлена высокой константой скорости этого процесса (k4 ~ 1.4 × 10–9 см3/с при Te = 2.5 эВ, по сравнению с k1 ~ 2.3 × 10–11 см3/с и k2 ~ 3.0 × 10–11 см3/с для x = 4) и значительной концентрацией молекул F2 в плазме из-за их эффективного образования по R31. Соответственно, совокупность гетерогенных процессов F + CFx → CFx + 1 и F + F → F2 в составе R30 и R31 более чем на 90% формирует суммарную скорость гибели атомарного фтора.
Увеличение доли H2 в смеси обогащает газовую фазу атомами водорода (за счет роста скорости диссоциации молекул H2 по R8, примерно в 35 раз при 3–13% H2), что приводит к перераспределению параллельных каналов рекомбинации в составе R30 и R31 в пользу CFx + H → CHFx. Так как основным фторуглеродным радикалом при отсутствии водорода в смеси является CF3, в качестве доминирующего продукта R30 и R31 ожидаемо выступает CHF3 (см. рис. 2, а). Немонотонное (с максимумом при 6–8% H2) поведение концентрации этих частиц обусловлено аналогичным изменением скорости их образования по причине нарастающего дефицита CF3. Причиной этого является не только замедление генерации CF3 по R1 из-за снижения содержания тетрафторметана в плазмообразующем газе, но и рост частоты гибели этих частиц при взаимодействии с атомами водорода по R26. Кроме того, группа атомно-молекулярных процессов R24–R28 обеспечивает эффективную конверсию F2, CFx (x = 1–3) и CHFx (x = 1–3) в молекулы HF. Это приводит к выполнению условия nHF > и доминированию молекул HF уже при содержаниях водорода в смеси более 5–6%. Отметим также, что константы скоростей R27 и R28 увеличиваются примерно на порядок величины при переходе от более насыщенных к менее насыщенным исходным частицам (например, k27 ~ 3.0 × 10–12 см3/с для x = 3 vs. ~ 3.2 × 10–11 см3/с для x = 2 и 1). Таким образом, преимущественная генерация радикалов с низкими x существенно замедляет снижение концентрации CF2 (примерно в два раза при 3–13% H2) по сравнению CF3 (примерно в 20 раз при 3–13% H2) и приводит к росту концентрации CF. Еще одним важным моментом является тот факт, что частота гибели атомов фтора в паре R25 + R27 превышает совокупный эффект от группы гетерогенных процессов R30–R32. Поэтому добавка водорода приводит к резкому росту частоты гибели атомов F и снижению их концентрации (примерно в 8 раз при 3–13% H2, см. рис. 2, а) на фоне слабого изменения скорости образования этих частиц в процессах электронного удара (см. рис. 2, б). Последний эффект обусловлен тем, что снижение суммарной скорости R1 + R2 почти полностью компенсируется R17, в том числе за счет k17 > k1 + k2. Концентрация атомов водорода ожидаемо увеличивается с ростом доли H2 в смеси, однако характер этого изменения является более слабым по сравнению с (см. рис. 2, а). Причинами здесь являются снижение эффективности образования атомов по R8 из-за аналогичного измерения концентрации электронов, а также их расходование в атомно-молекулярных процессах R24, R26 и R28. Последние составляют реальную конкуренцию гетерогенной рекомбинации атомарного водорода по “классическому” механизму H + H → H2 в составе R32.
По результатам оптико-спектральных измерений было найдено, что замещение CF4 на H2 при yAr = const вызывает монотонное снижение интенсивности линии Ar 750.4 нм (примерно в 1.3 раза при 3–13% H2, рис. 3, а). Очевидной причиной этого является аналогичное изменение функции возбуждения ke, Arne (2.1–1.5 с–1 при 3–13% H2), обусловленное снижением концентрации электронов (см. рис. 1, а). Таким образом, результаты зондовой и спектральной диагностики плазмы находятся в согласии друг с другом. Из рис. 3, б можно видеть также, что имеет место хорошее согласие расчетных концентраций атомов фтора и водорода с результатами их актинометрического определения. Данный факт подтверждает адекватность кинетической схемы, использованной для моделирования плазмы, а также свидетельствует о корректности расчетных данных по составу газовой фазы.
Рис. 3. Интенсивности излучения (а) и концентрации атомов (б) в плазме смеси CF4 + H2 + Ar. На рис. а: 1 — Ar 750.4 нм; 2 — F 703.4 нм; 3 — H 656.4 нм. На рис. б: 2 — атомы фтора; 3 — атомы водорода. Пунктирные линии на рис. б представляют значения, полученные при моделировании плазмы.
Еще одним важным эффектом увеличения доли водорода в смеси является рост суммарной концентрации (примерно в 2.8 раз при 3–13% H2, см. рис. 2, а) и плотности потока (Γpol = 1.6 × 1016–5.7 × 1016 см–2с–1, или примерно в 3.6 раз при 3–13% H2) полимеробразующих радикалов СHxFy (x + y < 3). Противоположные тенденции изменения последнего параметра и плотности потока атомов фтора обеспечивают резкий рост отношения Γpol / ΓF (0.4–11.8, или примерно в 25 раз при 3–13% H2), отслеживающего изменение полимеризационной способности плазмы [20, 22, 34]. Ранее в работе [34] было показано, что аналогичные по абсолютной величине изменения отношения Γpol / ΓF достигаются в смесях CF4 + CHF3 + Ar и CF4 + C4F8 + Ar при более чем половинном замещении тетрафторметана на второй фторуглеродный компонент. Таким образом, комбинирование тетрафторметана с водородом обеспечивает более широкие диапазоны регулирования полимеризационной способности плазмы при меньшей концентрации газа-добавки. Данная особенность потенциально расширяет возможности плазмы CF4 в процессах РИТ, требующих высокой селективности травления в паре SiO2/Si.
ЗАКЛЮЧЕНИЕ
Проведено исследование электрофизических параметров и состава плазмы в смеси CF4 + H2 + Ar в условиях индукционного ВЧ (2 МГц) разряда. Установлено, что замещение тетрафторметана на водород при постоянном содержании аргона не приводит к существенным изменениям температуры электронов и суммарной частоты ионизации, но сопровождается снижением плотности и ростом электроотрицательности плазмы. Причиной последнего эффекта является рост суммарной частоты прилипания, обусловленный эффективным образованием молекул HF. При анализе кинетики нейтральных частиц найдено, что снижение концентрации атомов фтора с ростом доли водорода в смеси является значительно более резким по сравнению с изменением скорости их образования в плазме. Такая ситуация обусловлена увеличением частоты гибели атомов в реакциях вида CHFx + F → CFx + HF, протекание которых обеспечивается гетерогенной рекомбинацией по механизму CFx + H → CHFx. Удовлетворительное согласие расчетных и экспериментальных данных по концентрациям атомов фтора и водорода свидетельствует об адекватности кинетической схемы, использованной для моделирования плазмы. Еще одним важным эффектом увеличения содержания водорода в плазмообразующем газе является рост концентрации радикалов СHxFy (x + y < 3). Это свидетельствует о росте полимеризационной нагрузки плазмы на контактирующие с ней поверхности.
ФИНАНСИРОВАНИЕ
Исследование выполнено за счет гранта Российского научного фонда № 23-29-0077 1, https://rscf.ru/project/23-29-00771/.
КОНФЛИКТ ИНТЕРЕСОВ
Авторы заявляют, что у них нет конфликта интересов.
Sobre autores
A. Miakonkikh
Valiev Institute of Physics and Technology of the Russian Academy of Sciences
Autor responsável pela correspondência
Email: miakonkikh@ftian.ru
Rússia, Moscow
V. Kuzmenko
Valiev Institute of Physics and Technology of the Russian Academy of Sciences
Email: miakonkikh@ftian.ru
Rússia, Moscow
A. Efremov
Ivanovo State University of Chemistry and Technology
Email: miakonkikh@ftian.ru
Rússia, Ivanovo
K. Rudenko
Valiev Institute of Physics and Technology of the Russian Academy of Sciences
Email: miakonkikh@ftian.ru
Rússia, Moscow
Bibliografia
- Wolf S., Tauber R.N. Silicon Processing for the VLSI Era. Volume 1. Process Technology. New York: Lattice Press, 2000.
- Nojiri K. Dry etching technology for semiconductors. Tokyo: Springer International Publishing, 2015.
- Lieberman M.A., Lichtenberg A.J. Principles of plasma discharges and materials processing. New York: John Wiley & Sons Inc., 1994.
- Stoffels W.W., Stoffels E., Tachibana K. Polymerization of fluorocarbons in reactive ion etching plasmas, J. Vac. Sci. Tech. A. 1998. V. 16. P. 87–95.
- Standaert T.E.F.M., Hedlund C., Joseph E.A., Oehrlein G.S., Dalton T.J. Role of fluorocarbon film formation in the etching of silicon, silicon dioxide, silicon nitride, and amorphous hydrogenated silicon carbide, J. Vac. Sci. Technol. A. 2004. V. 22. P. 53–60.
- Schaepkens M., Standaert T.E.F.M., Rueger N.R., Sebel P.G.M., Oehrlein G.S., Cook J.M. Study of the SiO2-to-Si3N4 etch selectivity mechanism in inductively coupled fluorocarbon plasmas and a comparison with the SiO2-to-Si mechanism, J. Vac. Sci. Technol. A. 1999. V. 17. P. 26–37.
- Plumb I.C., Ryan K.R. A model of the chemical processes occurring in CF4/O2 discharges used in plasma etching. Plasma Chem. Plasma Process. 1986. V. 6. P. 205–230.
- Kimura T., Noto M. Experimental study and global model of inductively coupled CF4/O2 discharges, J. Appl. Phys. 2006. V. 100. P. 063303 (1–9).
- Efremov A.M., Bashmakova D.E., Kwon K.-H. Features of plasma composition and fluorine atom kinetics in CHF3 + O2 gas mixture, Chem. Chem. Tech. 2023. V. 66. No. 1. P. 48–55.
- Vasenkov A.V., Li X., Oehlein G.S., M.J. Kushner. Properties of c-C4F8 inductively coupled plasmas. II. Plasma chemistry and reaction mechanism for modeling of Ar/c-C4F8/O2 discharges, J. Vac. Sci. Technol A. 2004. V. 22. P. 511–530.
- Lim N., Efremov A., Kwon K.-H. A comparison of CF4, CHF3 and C4F8 + Ar/O2 Inductively Coupled Plasmas for Dry Etching Applications. Plasma Chem. Plasma Process. 2021. V. 41. P. 1671–1689.
- Baek S.Y., Efremov A.M., Bobylev A.V., Choi G., Kwon K.-H. On relationships between plasma chemistry and surface reaction kinetics providing the etching of silicon in CF4, CHF3, and C4F8 gases mixed with oxygen, Materials. 2023. V. 16. P. 5043 (1–18).
- Marra D.C., Aydil E.S. Effect of H2 addition on surface reactions during CF4/H2 plasma etching of silicon and silicon dioxide films, J. Vac. Sci. Technol. A. 1997. V. 15. P. 2508–2517.
- Knizikevicius R. Real dimensional simulation of SiO2 etching in CF4 + H2 plasma, Appl. Surface Sci. 2004. V. 222. P. 275–285.
- Gorobchuk A. Numerical modeling of silicon processing technology in CF4/H2 plasma, Proceedings of 2015 International Siberian Conference on Control and Communications (SIBCON). Omsk, Russia, 2015. P. 1–4.
- Kim D.S., Kim J.B., Ahn D.W. et al. Atomic Layer Etching Applications in Nano-Semiconductor Device Fabrication, Electron. Mater. Lett. 2023. V. 19. P. 424–441.
- Kim Y., Kang H., Ha H., Kim C., Cho S., Chae H. Plasma atomic layer etching of molybdenum with surface fluorination, Applied Surface Science. 2023. V. 627. P. 157309.
- Kim S.Y., Park I.-S., Ahn J. Atomic layer etching of SiO2 using trifluoroiodomethane, Appl. Surf. Sci. 2022. V. 589. P. 153045.
- Kuzmenko V., Lebedinskij Y., Miakonkikh A. et al. Selective atomic layer etching of Al2O3, AlNx and HfO2 in conventional ICP etching tool, Vacuum. 2023. V. 207. P. 111585.
- Efremov A.M., Murin D.B., Kwon K.-H. Concerning the Effect of Type of Fluorocarbon Gas on the Output Characteristics of the Reactive-Ion Etching Process, Russian Microelectronics. 2020. V. 49. No. 3. P. 157–165.
- Efremov A., Lee J., Kwon K.-H. A comparative study of CF4, Cl2 and HBr + Ar inductively coupled plasmas for dry etching applications. Thin Solid Films. 2017. V. 629. P. 39–48.
- Efremov A., Lee B.J., Kwon K.-H. On relationships between gas-phase chemistry and reactive-ion etching kinetics for silicon-based thin films (SiC, SiO2 and SixNy) in multi-component fluorocarbon gas mixtures. Materials, 2021. V. 14. P. 1432 (1–27).
- Shun’ko E.V. Langmuir Probe in Theory and Practice. Boca Raton: Universal Publishers, 2008.
- Herman I.P. Optical diagnostics for thin film processing. San Diego: Academic Press, 1996.
- Lopaev D.V., Volynets A.V., Zyryanov S.M., Zotovich A.I., Rakhimov A.T. Actinometry of O, N and F atoms, J. Phys. D: Appl. Phys. 2017. V. 50. P. 075202 (1–17).
- Christophorou L.G., Olthoff J. K. Fundamental electron interactions with plasma processing gases. New York: Springer Science & Business Media, 2012.
- Raju G.G. Gaseous electronics. Tables, atoms and molecules. Boca Raton: CRC Press, 2017.
- Skoro N., Puac N., Lazovi S., Cvelbar U., Kokkoris G., Gogolides E. Characterization and global modelling of low-pressure hydrogen-based RF plasmas suitable for surface cleaning processes, J. Phys. D: Appl. Phys. 2013. V. 46. P. 475206 (1–13).
- Iordanova S., Koleva I., Paunska T. Hydrogen degree of dissociation in a low pressure tandem plasma source, Spectroscopy Letters. 2011. V. 44. P. 8–16.
- Kimura T., Kasugai H. Properties of inductively coupled rf Ar/H2 plasmas: Experiment and global model, J. Appl. Phys. 2010. V. 107. P. 083308 (1–9).
- Lavrov B.P., Pipa A.V. Account of the fine structure of hydrogen atom levels in the effective emission cross sections of Balmer lines excited by electron impact in gases and plasma, Optics and Spectroscopy. 2002. V. 92. No. 5. P. 647–657.
- Celik Y., Aramaki M., Luggenholscher D., Czarnetzk U. Determination of electron densities by diode-laser absorption spectroscopy in a pulsed ICP, Plasma Sources Sci. Technol. 2011. V. 20. P. 015022 (1–12).
- Cunge G., Ramos R., Vempaire D., Touzeau M., Neijbauer M., Sadeghi N. Gas temperature measurement in CF4, SF6, O2, Cl2, and HBr inductively coupled plasmas, J. Vac. Sci. Technol. A. 2009. V. 27. No. 3. P. 471–478.
- Efremov A., Son H.J., Choi G., Kwon K.-H. On Mechanisms Influencing Etching/Polymerization Balance in Multi-Component Fluorocarbon Gas Mixtures, Vacuum. 2022. V. 206. P. 111518 (1–10).
- Ho P., Johannes J.E., Buss R.J. Modeling the plasma chemistry of C2F6 and CHF3 etching of silicon dioxide, with comparisons to etch rate and diagnostic data, J. Vac. Sci. Technol. B. 2001. V. 19. P. 2344–2367.
- Proshina O.V., Rakhimova T.V., Zotovich A.I., Lopaev D.V., Zyryanov S.M., Rakhimov A.T. Multifold study of volume plasma chemistry in Ar/CF4 and Ar/CHF3 CCP discharges. Plasma Sources Sci. Technol. 2017. V. 26. P. 075005 (1–26).
- Gogolides D., Mary D., Rhallabi A., Turban G. RF Plasmas in Methane: Prediction of Plasma Properties and Neutral Radical Densities with Combined Gas-Phase Physics and Chemistry Model, Jpn. J. Appl. Phys. 1995. V. 34. P. 261–270.
- Herrebout D., Bogaerts A., Yan M. et al. One-dimensional fluid model for an rf methane plasma of interest in deposition of diamond-like carbon layers, J. Appl. Phys. 2001. V. 90. P. 570–579.
- Semenova O.A., Efremov A.M., Barinov S.M., Svettsov V.I. Kinetics and concentration of active particles in nonequilibrium low temperature methane plasma, High Temperature. 2014. V. 52. No. 3. P. 348–355.
- Adams N.G., Smith D. Dissociative attachment reactions of electrons with strong acid molecules, J. Chem. Phys. 1986. V. 86. P. 6728–6731.
- Abouaf R., Teillet-Billy D. Fine structure in the dissociative attachment cross sections for HBr and HF, Chem. Phys. Letters. 1980. V. 73. No. 1. P. 106–109.
- Xu Y., Gallup G.A., Fabrikant I.I. Dissociative electron attachment to vibrationally and rotationally excited H2 and HF molecules, Phys. Rev. A. 2000. V. 61. P. 052705 (1–7).
- Chantry P.J. A simple formula for diffusion calculations involving wall reflection and low density, J. Appl. Phys. 1987. V. 62. P. 1141–1148.
Arquivos suplementares
