Вычислительные технологии
Том 12, № 5, 2007
ВЛИЯНИЕ ВЧ-РАЗРЯДА НА ПРОЦЕСС ПЛАЗМОХИМИЧЕСКОГО ТРАВЛЕНИЯ КРЕМНИЯ В CF4/O2*
А. Г. Горовчук, Ю. Н. Григорьев Институт вычислительных технологий СО РАН, Новосибирск, Россия
e-mail: al@ict.nsc.ru
The effect of RF-discharge on silicon etching process in tetrafluoromethane — oxygen mixture in a plasma-chemical reactor is investigated numerically. It was shown that the decrease of average density of energy electrons in the reactor due to the oxygen inflow can decrease the etching rate by up to 30%.
Введение
Одним из важнейших технологических процессов, применяемых в производстве изделий микроэлектроники, является плазмохимическое травление полупроводниковых материалов. Травление поверхности обрабатываемого материала производится химически активными частицами, генерируемыми в низкотемпературной плазме высокочастотного (ВЧ) разряда. Скорости объемных процессов ионизации, генерации и рекомбинации частиц, участвующих в обработке полупроводниковых материалов, определяются внутренними характеристиками ВЧ-разряда — концентрациями и функциями распределения электронов и ионов по энергиям, а также сечениями соответствующих реакций. Изменение концентраций электронов и ионов в разряде способно существенно повлиять на основные характеристики процесса травления.
Математическое моделирование плазмохимического травления предоставляет широкие возможности для поиска взаимосвязей факторов, определяющих качество и скорость обработки образцов. В результате исследований [1-3] разработана адекватная численная модель реакторного процесса на основе уравнений многокомпонентной гидродинамики с учетом конвективно-диффузионного переноса отдельных компонент смеси, позволяющая исследовать тонкие физические эффекты процесса плазмохимическо-го травления. На основе данной модели исследовались особенности интенсификации травления кремния в плазме CF4/O2. В [1-3] показано, что добавленный к CF4 кислород, образуя углеродсодержащие соединения, способствует возрастанию концентрации активных атомов фтора, повышая в несколько раз скорость травления образца. Было
* Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 05-01-00359, № 07-01-00315), президентской программы поддержки ведущих научных школ РФ (грант № НШ-9886.2006.9).
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2007.
найдено, что максимум средней скорости травления достигается в области с относительно низким содержанием кислорода. При добавках 02, больших 25%, наблюдается образование хемосорбционного слоя кислорода на кремнии, приводящие к эффекту гистерезиса на графике зависимости средней скорости травления от концентрации фтора при различном содержании 02 [4]. Исследование данного эффекта в зависимости от интенсивности хемосорбции кислорода на кремнии выполнено в [3]. В этой работе показано, что абсолютное значение и положение максимума средней скорости травления в зависимости от процентного содержания 02 в смеси полностью определяется процессом хемосорбции кислорода на кремнии.
Вместе с тем влияние добавки кислорода неоднозначно и может оказать также и отрицательное воздействие. Согласно данным спектральных измерений [4-7], процесс травления кремния в смеси СГ4/02 по мере возрастания содержания кислорода сопровождается уменьшением концентрации электронов при постоянной мощности, вкладываемой в разряд, и постоянном давлении. Это приводит к существенному снижению скорости генерации активных частиц и, как следствие, к уменьшению их концентрации. Дальнейшие исследования [5] показали, что снижение электронной плотности не зависит от конструкционных размеров реакционной камеры и параметров рабочего режима. Кроме того, в [6] было показано, что падение концентрации электронов не зависит от их энергии, по крайней мере в диапазоне электронных энергий 8 — 20 эВ. Следовательно, эксперименты [4-6] показывают, что влияние добавки кислорода в СГ4 на плотность электронов в ВЧ-разряде носит достаточно универсальный характер.
Таким образом, присадка кислорода в СГ4 вызывает прямо противоположные процессы: с одной стороны, увеличивается производство активных частиц Г одновременно с обеднением смеси радикалами СГХ; с другой — идет образование хемосорбционного слоя на кремнии, затрудняющего доступ активных частиц к образцу, снижение плотности электронов и, как следствие, эффективности диссоциации молекул исходной смеси. До настоящего времени на основе математического моделирования эффект снижения электронной плотности в СГ4/02 практически не рассматривался. Известны лишь отдельные работы, в которых количественные оценки его влияния на характеристики процесса травления отсутствуют. В этой связи представляет интерес определить степень снижения производства активных частиц в плазме ВЧ-разряда при учете данного эффекта.
В работе исследуется влияние снижения электронной плотности в ВЧ-разряде при добавке 02 в СГ4 на скорость и однородность травления образцов. Рассматривается зависимость формы кривой гистерезиса от производства активных частиц при разной интенсивности хемосорбции атомарного кислорода на кремнии. Исследуется эффективность диссоциации исходной газовой смеси СГ4/02 в условиях пониженной электронной плотности в ВЧ-разряде.
Теоретический анализ
В работе моделируется полностью загруженный плазмохимический реактор радиальной схемы (рис. 1). На рисунке введены обозначения: г — радиус реакционной камеры реактора, Ь — межэлектродное расстояние, г0 — радиус нижнего электрода, г — радиус выходного отверстия, Н — длина выходной трубки. Течение газовой смеси обозначено стрелкой и направлено к центру реактора. Торцевые поверхности реакционной камеры
реактора служат электродами, между которыми инициируется высокочастотный разряд. На нижнем электроде равномерно размещены полупроводниковые пластины кремния (г г < г < г0, г = Н). Газовая смесь поступает в реакционную камеру через кольцевую поверхность с распределенными отверстиями (г0 < г < г^, г = 0), расположенную на периферии нижнего электрода. В разрядной зоне реактора в результате реакций диссоциации исходной газовой смеси электронным ударом образуются химически активные частицы, участвующие в травлении полупроводниковых пластин. Продукты реакций и непрореагировавший газ удаляются в аксиальном направлении через выпускное отверстие в центре реактора (0 < г < гг, г = 0). Конструкционные размеры и режимы работы реактора брались из работы [8].
Как было показано в [4 - 7], влияние добавки 02 в СР4 на электронную плотность не зависит от выбора конструктивной схемы и режимов работы реактора. Это позволило
Рис. 1. Схема радиального плазмохимического реактора: 1,2 — ВЧ-электроды; 3 — обрабатываемые пластины кремния; 4 — зона плазменного разряда; 5 — подача исходного газа; 6 — откачка прореагировавшей смеси
Рис. 2. Зависимость нормированной электронной плотности от процентного содержания О2 в
СР4/О2 [4]
использовать в данной работе экспериментальные данные по снижению электронной плотности [4]. Зависимость электронной плотности от содержания 02 в СГ4 описывалась соотношением пе = пе0 • (хо2). Здесь пе0 — однородная по объему реактора плотность электронов в отсутствие кислорода, хо2 — молярная доля кислорода, $е(жо2) — графическая зависимость из работы [4] (рис. 2).
Численное моделирование проводилось на основе двумерной математической модели плазмохимического реактора травления [1-3], в которой уделялось особое внимание многокомпонентной химической кинетике газофазных и гетерогенных реакций. В расчетах использовались 14 газофазных химических реакций из работы [9], которые при моделировании процесса травления кремния в смеси СГ4/02 удовлетворительно предсказывали концентрации реагентов Г, С0, С02, С0Г2. Для описания процессов адсорбции радикалов СГ2, СГ3 на кремнии брались гетерогенные реакции с участием Г, СГ2, СГ3, СГ4, 0 из работ [8, 10]. В результате кинетика химических реакций включала следующие процессы:
СГ4 + е- —^ке1 СГз + Г + е-, (1)
СГ4 + е- —>ке2 СГ2 + 2Г + е-, (2)
02 + е- — >ке3 0 + 0 + е-, (3)
С0Г2 + е- —>ке4 С0Г + Г + е-, (4)
С02 + е- —>кеБ С0 + 0 + е-, (5)
СГз + СГз + М - С2Гб + М, (6)
Г + СГз + М - СГ4 + М, (7)
Г + СГ2 + М - СГ3 + М, (8)
0 + СГз - С0Г2 + Г, (9)
0 + СГ2 - С0Г + Г, (10)
0 + СГ2 - С0 + 2Г, (11)
0 + С0Г —С02 + Г, (12)
Г + С0Г + М — ^ С0Г2 + М, (13)
Г + С0 + М - С0Г + М, (14)
Г + Г + М —Л10 Г2 + М, (15)
Г 2 + М —^ к.и р + Г + М, (16)
СГз —^кв1 СГз (в), (17)
СГ2 —СГ2 (в), (18)
Г + СГ2 (в) —СГз, (19)
Г + СГз (в) —>к*4 СГ4, (20)
СГз + СГз (в) —^6, (21)
СГ2 (в) + 0 - —С0 + 2Г, (22)
СГз (в) + 0 - —С0 + 3Г, (23)
0 —Л8 0 (в), (24)
0 (в) + Г —>км 0 + Г, (25)
4Г + —Й1Г4 т . (26)
Набор реакций содержит 16 газофазных и 10 гетерогенных реакций на образце, в которых принимают участие 12 реагентов — Г, Г2, СГ2, СГ3, СГ4, С2Г6, О, 02, СО, С02, СОГ, С0Г2. Реакции (1)-(5) описывают диссоциацию молекул тетрафторметана и кислорода электронным ударом с образованием химически активных атомов фтора и кислорода, ке1 — ке5 — соответствующие константы скоростей диссоциации; реакции (6)-(16) представляют объемную рекомбинацию атомов и радикалов с участием третьего тела М, — к^ц — соответствующие константы скоростей рекомбинации; реакции (17)-(23) — гетерогенные реакции адсорбции-десорбции радикалов СГ2 и СГ3 на поверхности кремния, кз1 — кз7 — константы скоростей гетерогенных реакций; реакции (24) и (25) — реакции хемосорбции фтора и кислорода на поверхности кремния, кз8, кз9 — константы реакций хемосорбции; (26) — реакция спонтанного травления образца, к3 — соответствующая константа скорости реакции. Обозначение (в) использовано для радикалов и частиц, адсорбированных на поверхности образца. Значения констант реакций взяты из работ [8-12].
Распределение концентраций компонент находилось из системы уравнений конвективно-диффузионного переноса:
V -УСг = V- (аА-т(Ужг)) + Сг{Сг,С3), г,з = 1,..., 12, (27)
где Сг, хг — молярная концентрация и молярная доля частиц сорта г; С — молярная концентрация газовой смеси; Пг-т — коэффициент диффузии г-й компоненты в многокомпонентной смеси; Сг — источник частиц г-го сорта, г = Г, Г2, СГ2, СГ3, СГ4, С2Г6,
О, 02, СО, С02, СОГ, С0Г2.
Правая часть системы содержит набор определяющих газофазных реакций (1) —(16), которые устанавливают сложные взаимозависимости скоростей генерации частиц, записываемые в виде закона действующих масс:
i
аг = Х^Т — ¿гт )[*> П х7г — кьг\\ Х/г Ь
т=1 j=1 j=1
где игт, игт — стехиометрические коэффициенты прямых и обратных реакций для вещества г в реакции г = 1, 2,...,1; kfт, кьт — соответственно константы прямой и обратной реакции г. Структура источникового слагаемого такова, что вносит в уравнения степенную нелинейность относительно концентраций частиц. Гетерогенные реакции (химического травления, рекомбинации, адсорбции и др.) учитывались в краевых условиях на образце. Краевые условия на границах расчетной области записывались в виде локального баланса потоков соответствующей компоненты.
Физические характеристики среды — коэффициенты вязкости и теплопроводности смеси вычислялись по аппроксимационным формулам Вильке и Мэзона — Саксены, известным из кинетической теории газов [13]. Коэффициент диффузии многокомпонентной смеси вычислялся по закону Бланка [14]:
х х -1 Г г = Г, СГ2, СГз, Г2, О, СО, СО2, СОГ, СОГ2;
Пг-т =(1 — хг) £ П- ,
^ [ 3 = СГ4, О2,
где в качестве основных компонент смеси брались СГ4 и О2. Здесь П^ — коэффициенты бинарной диффузии. Формулы учитывают зависимость диссипативных коэффициентов
от температуры и характеристик чистых газов. Межмолекулярное взаимодействие описывалось потенциалом Леннарда — Джонса.
Важное значение в процессе травления имеет пассивация поверхности кремния адсорбентами. Конкурирующие процессы адсорбции радикалов СР2, СЕ3 и хемосорбции атомарного кислорода на кремнии закрывают доступ атомам фтора к образцу и препятствуют спонтанному травлению кремния. Обозначим доли поверхности кремниевой пластины, покрытой СР2, СЕ3 и хемосорбированным кислородом, как "ср2 , "ср3 , "о соответственно. Равновесные балансы массовых потоков на образце для компонент СР2, СЕ3, О дают следующие соотношения на неизвестные "ср2, "ср3, "о [8, 10]:
к«2ЖСР2/(к*Жр + кзбЖо) = "сР2/(1 - "ср2 - "срз - "о),
ЖсрЗ/(А^Жр + А^ЖсРз + к^жо) = "срЗ/(1 - "ср2 - "срЗ - "о), а^жо/жр = "о/(1 - "ср2 - "срз)• Решением системы линейных уравнений являются:
к«2ЖсР2
"СР2 АСР2
к«2 ЖСР2 + к«3 ЖР + к«6Жо + АСР2 А1,СР2 + А2,СР2 >
а1,СР2 а2,СР2
"СРз АСРз а1,СРз А2,СРз
"о
А
А
= к51ЖсРз
жо
2,0
к«3 ЖР + к«6Жо
к^жр + к5бЖср3 + к«7Жо
а— (к^зжр + к^бЖо);
жр
к«1ЖсРз
к«1жсРз + к^жр + к^бжсРз + Асрз
+ А
1,СРз + А2,СРз > к«4ЖР + к«5ЖСР3 + к«7Жо
к^з жр + к^бЖо
к«2ЖСР;
а — (к^4Жр + к^5Жср3 + к«7Жо);
жр
а жо
1,0 = ЖР
а«жо + жр + А1,о + А2,о; к«2ЖсР2
жр
к«зЖр + к^бЖо'
_к^1ЖсРз_
к^4Жр + к^5Жср3 + к«7Жо
(28)
(29)
Полученные формулы имеют простой вид, в них неизвестные характеристики адсорбционных слоев выражены через константы реакций. Здесь параметр а8 = к^/к^д пропорционален отношению двух кинетических процессов — хемосорбции кислорода на кремнии (реакция 24), при котором образуется пленка О, предотвращающая процесс травления кремния, и его десорбции (реакция 25) [11, 12]. Значение параметра а характеризует интенсивность хемосорбции кислорода на кремнии.
Частными случаями приведенных формул являются следующие. Если в (28) положить хср2 = хо = 0, то получим выражение для $ср3 из модели трехкомпонентной химической кинетики [10]:
к^ХсЕз
^СЕз "
(к«1 + к5б)хсрз + к«4хр
Если в (29) пренебречь компонентами хср2,хсрз С хр,хо, то получим $о из работы [12]:
_ а«хо
"о ~-:-•
а«хо + хр
Первое выражение отражает только адсорбцию радикалов СГ3, второе — чистую "конкуренцию" процессов хемосорбции атомарных фтора и кислорода на кремнии.
Численные расчеты показывают, что $ср2, "$се3 С 1 [1, 2], и неизвестный параметр а« можно приближенно оценить по соотношению адсорбированных на образец частиц атомов фтора и кислорода:
к«8 ^о 5о /ШР кТ .
а« = тг = -, ^ = \ п-, З = F, 0, (30)
к«9 ¿р у шо у 2пШj
где Шj, Fj — молекулярная масса и интенсивности молекулярных потоков частиц З-го сорта, £р,$о — коэффициенты прилипания фтора и кислорода на кремнии [12].
Доли поверхности кремниевого образца, покрытые адсорбированными СГ2, СГ3 и О, являются важными параметрами задачи, поскольку входят в краевые условия на образце и существенно влияют на скорость травления. Скорость спонтанного травления, ее среднее значение, А/мин, и индекс неоднородности травления вычислялись по следующим формулам [8, 10]:
V« = 1.81 ■ 1010 (1 — ^о — ^СР2 — ^СЕз) к«Ср, (31)
Ко
— Щ!
К
= Vs2 — Vsl
1гП
2г« '
где к« — константа скорости реакции спонтанного травления, см/с; Ср — молярная концентрация фтора, моль/см3; г«2, г51 — максимальная и минимальная скорости травления. Из-за адсорбции СГ2, СГ3, О на образце скорость травления является сложной нелинейной зависимостью от компонент смеси — Г, СГ2, СГ3, О. Линейная зависимость скорости травления от концентрации фтора может быть использована в диапазоне добавок О2, где $се2 ,^се3 , С 1.
Решение сформулированной задачи осуществлялось численным методом, кратко изложенным в [15]. Для решения системы уравнений (27) использовался метод конечных разностей. Расчетная область по пространственным переменным покрывалась равномерной сеткой П^. Число узлов в ее верхней части (0 < г < г^, Н < г < Н + Ь) составляло 70 х 45. Искомое решение находили методом стационарных итераций с постоянным шагом т. В сеточной области П каждое из уравнений для относительных молярных
концентрации хг аппроксимировалось неявной итерационном конечно-разностной схемой расщепления со стабилизирующей поправкой:
к+1/2 ь
Х__ ГХхк+1/2 + щх+ ^
-— Щ х, + щ х, + ^
к
хк+1 - хк+1/2
^г_ - _ ГХ (хк+1 — )
(хг хг ) ,
Т
где символы хк обозначают сеточные функции решения хк п т на к-временном слое. Нижние индексы г,п,т в хк п т относятся к частицам г-го сорта в (п,т)-м узле сетки. Аппроксимируемые операторы , Щ1 включали в себя первые и вторые производные от искомых функций по направлениям г и г соответственно. Для построения разностных операторов Щх, ЬХ во внутренних узлах сетки П использовались центральные разности второго порядка аппроксимации производных по пространственным переменным. Полученные конечно-разностные уравнения аппроксимировали исходные дифференциальные уравнения (27) с точностью 0(т + Л-1 + Л|), где Л,2 — шаги в направлениях г и г. Каждое из конечно-разностных уравнений, записанное в дробных шагах на трехточечном шаблоне, представлялось на сетке П в виде системы линейных алгебраических уравнений с трехдиагональной матрицей. Такая система линейных уравнений эффективно решалась методом прогонки, являющейся пространственно устойчивой.
Отметим особенности организации вычислительного алгоритма. Решение уравнений движения вязкого теплопроводного газа находилось независимо от решения системы уравнений переноса (27). Уравнения Навье—Стокса в приближении Буссинеска рассчитывались в переменных "функция тока—вихрь" совместно с уравнением баланса энергии. Краевые условия для вектора скорости записывались в этих же переменных. При решении уравнения переноса вихря на каждом новом временном шаге Т решались уравнения для функции тока и температуры. По найденным распределениям скоростей потока газа и температуры из системы уравнений конвективно-диффузионного переноса (27) находились концентрации компонент смеси. Далее определялись характеристики спонтанного травления — скорость и индекс неоднородности травления.
Рассмотрим более подробно схему построения итерационного процесса вычисления относительных концентраций компонент хг. Установившиеся значения компонент вектора скорости и температуры подставляются в систему уравнений (27) и строится итерационный процесс вычисления относительных концентраций компонент хг. Пусть в момент времени в узлах сетки известны значения сеточных функций хк. Для перехода от момента к моменту времени ¿к+1 — + т последовательно решаем конечно-разностный аналог уравнения (27) для каждой г-й компоненты. Все значения хк+1 на новом временном слое вычисляются одновременно. Заметим, что источниковые члены уравнений (27) и краевые условия на границах расчетной области зависят от искомых величин хк нелинейным образом. Для разрешения этой нелинейности значения хк в ис-точниковых членах и в краевых условиях на границах брались с предыдущей итерации к. Итерационный процесс продолжается до установления приближенного стационарного решения. Стационарное решение считается найденным, если выполняется следующее условие:
тах
г ,
хк+1 _ хк гг
Тхк+1
< 10-10.
Для оценки точности численного алгоритма были выполнены тестовые расчеты радиального изотермического реактора для расчетных данных, совпадающих с данными работы [16], и сопоставление результатов. Расчеты проводились в рамках трехком-понентной химической кинетики на регулярной сетке размерами 200 х 150. Отмечено качественное совпадение сопоставляемых результатов. Для количественной оценки различия скоростей травления рассчитывались максимальные локальное ||б||с и среднеквадратичное \\^\\ь2 отклонения. Найдено, что их значения не превосходят величин 0.03873 и 0.03131 соответственно. Моделирование неизотермического плазмохимическо-го реактора радиальной схемы непрерывного действия с трехкомпонентной химической кинетикой [10] на регулярной сетке размерами 200 х 150 показало аналогичные результаты. Верификация модели плазмохимического реактора с 12-компонентной химической кинетикой [8] выполнялась на качественном уровне на основе сравнения с экперимен-тальным фактом существования гистерезиса на графике зависимости скорости плазмо-химического травления от концентрации фтора. В целом тестовые расчеты показали вполне удовлетворительное совпадение результатов, полученных для радиального плаз-мохимического реатора в рамках различных независимых численных моделей.
Результаты и обсуждение
С целью сокращения изложения полученных результатов рассмотрим влияние понижения электронной плотности в реакторе на примере рабочего режима р = 0.5 Торр, Q = 200 см3/мин, Т* = 300 X. На приведенных ниже рисунках данные соответствуют выбранному режиму. При других рассчитанных в работе режимах качественные показатели практически не изменяются.
Рис. 3. Зависимость средней скорости травления от процентного содержания О2 в СР4/02: 1 - = 1, а* = 0.3; 2 - §е = 1, а* = 1; 3 - = ), а* = 0.3; 4 - §е = §е(хо2), а* = 1
Добавка кислорода позволяет существенно повысить концентрацию активных частиц и скорость травления кремниевых подложек. Это свойство иллюстрируется зависимостью средней скорости спонтанного травления от процентного содержания 02 в СГ4/02 (рис. 3, кривые 1, 2). Вместе с тем по мере увеличения содержания кислорода в тетрафторметане при постоянных прочих условиях концентрация электронов монотонно уменьшается. Это снижает скорости диссоциации и ионизации электронным ударом, что приводит к уменьшению производства активных атомов фтора. Одновременно снижается средняя скорость травления (рис. 3, кривые 3, 4). Как видно из графиков, заметное расслоение кривых наступает в области максимальных скоростей травления при концентрациях кислорода, больших 30%.
Для характеристики снижения средней скорости травления и эффективности диссоциации молекул исходного газа вычислялись относительные отклонения:
=
г'8 - г''8
Н> Н"
С г — С г
С 'г
г = Г, СР4, 02,
где г'8, С 'г — средние скорость травления и концентрации фтора, тетрафторметана, кислорода соответственно при = 1; г''8, С''г — средние значения с учетом зависимости "$е(хо2). В таблице приведены ^, , %СР , ^с0 в зависимости от добавки 02. С увеличением процентного содержания 02 в СГ4 в пределах до 60 % падение средней скорости травления может составлять до 30 %. Однако качественных изменений в распределении средней скорости травления от процентного содержания 02 не наблюдается. Положение максимума средней скорости травления сохраняется и зависит только от интенсивности хемосорбции атомарного кислорода на кремнии — параметра (рис. 3, кривые 1, 2 и 3, 4). В расчетах рассматривались два значение этого параметра — 0.3 и 1. Первое (базовое значение) соответствует выбору коэффициентов прилипания для атомарных фтора и кислорода на основе анализа одной из схем химических реакций в сравнении с экспериментальными данными [12], второе — более реалистичному сценарию процесса травления [3]. С увеличением параметра относительные отклонения средней скорости травления, связанные с уменьшением электронной плотности, сохраняются.
Сравнение кривых 1,2 с кривыми 3, 4 на рис. 4, показывает, что с уменьшением концентрации электронов в ВЧ-разряде производство активных частиц падает. Количественные изменения средней концентрации фтора на образце в диапазоне добавок кислорода 0-40 % совпадают с изменениями средней скорости травления. Последнее свидетельствует о практически линейной зависимости скорости травления от концентрации фтора и слабом эффекте адсорбции частиц на образце в указанном диапазоне. Положение максимального значения концентрации фтора в зависимости от процентно-
Относительные отклонения средних скорости травления, концентраций фтора, тетрафторметана и кислорода в зависимости от содержания О2 в СР4 при а8 = 0.3
г
8
%, О2 ^, % % р, % %СР4 , % %02 ,%
10 8.95 8.95 2.05 25.98
20 14.79 14.74 8.41 42.07
30 21.22 21.07 22.25 58.75
40 24.88 24.44 44.15 66.49
50 26.66 26.86 68.82 74.99
60 28.90 29.66 81.61 84.58
Рис. 4. Зависимость средней концентрации фтора на образце от процентного содержания О2 в СР4/О2: 1 - ^е = 1, а* = 0.3; 2 - §е = 1, а* = 1; 3 - ^ = ^е(жо2), а* = 0.3; 4 - §е = ^е(жо2), а» = 1
го содержания 02 сохраняется. С ростом параметра а8 относительные изменения средней концентрации фтора на образце, связанные с уменьшением электронной плотности, также сохраняются.
Изменение скоростей объемных процессов ионизации, диссоциации и рекомбинации частиц отражается на расчетных концентрациях компонент смеси (рис. 5). В целом у поверхности образца отмечается снижение средних концентраций Г, Г2, СО, С02 с сохранением положения максимумов компонент в зависимости от добавки 02. Значения концентраций СГ2, СГ3 на 1-3 порядка ниже, чем концентрации других компонент смеси, и монотонно спадают по мере увеличения содержания 02. Для атомарного кислорода, дающего существенный вклад в конкурирующие процессы на образце, наблюдается заметное снижение концентрации, но в диапазоне добавок 02, больших 40%. Таким образом, баланс конкурирующих процессов на образце не изменяется, а следовательно, не изменяется и сценарий процесса травления. Единственное качественное, но несущественное отличие наблюдается для компоненты С0Г2, для которой в диапазоне добавок до 40 % 02 концентрация падает со смещением ее максимума с 25 % содержания 02 к 20 %. Изменение концентраций компонент смеси затрагивает и компоненты исходной газовой смеси — концентрации СГ4 и 02 увеличились, что свидетельствует об их менее интенсивном расходе в реакциях диссоциации электронным ударом (1)-(3). Относительные изменения концентраций СГ4 и 02, представленные в таблице, показывают, что понижение электронной плотности в ВЧ-разряде, связанное с добавкой 02, ведет к увеличению доли компонент исходной газовой смеси. При больших добавках кислорода эффективность диссоциации СГ4/02 падает практически вдвое.
Расчет концентраций реагентов показывает (рис. 5), что концентрация Г у образца превосходит концентрации радикалов СГ2 и СГ3 не менее, чем на 2 и 3 порядка соответственно. Концентрация атомарного кислорода по мере увеличения содержания 02 в
Рис. 5. Зависимость средних концентраций компонент смеси на образце от процентного содержания 02 в СР4/О2 при §е = $е(жо2) и а* = 0.3: 1 - Р; 2 - СР4; 3 - Р2; 4 - О; 5 - О2; 6 -СО; 7 - СО2; 8 - СОР2
СГ4/02 растет и становится того же порядка, что и концентрация Г. В связи с этим можно ожидать, что среди конкурирующих процессов адсорбции радикалов СГ2, СГ3 и хемосорбции 0 на кремнии будет преобладать последний. Расчеты показывают, что доли поверхности кремния, покрытые адсорбентами СГ2, СГ3, не превосходят 0.2 и 0.01 % соответственно. С увеличением добавки кислорода 02 значения $ср2 , "$ср3 монотонно падают, что объясняется расходованием СГ2, СГ3 в реакция рекомбинации (6)-(8) и реакциях взаимодействия кислорода с углеродом (9)-(12), в результате которых образуются С0, С02, С0Г, С0Г2. При добавках 02, больших 50%, $ср2, "$ср3 практически равны нулю. Последнее связано с общим обеднением смеси радикалами типа СГХ.
Доля поверхности кремния, занятая хемосорбированным кислородом, растет пропорционально росту концентрации атомарного кислорода (рис. 6). При небольшом содержании 02 сравнительно небольшая часть поверхности кремния связана кислородом (в пределах 1-2%). Значение $о становится заметным при добавках 02, больших 35%. Начиная с 45%-й добавки 02, параметр $о быстро возрастает с ростом доли 02, устремляясь к 1 при 100%-м содержании 02. Последнее связано с балансом двух процессов — увеличением доли атомарного кислорода и обеднением смеси атомами фтора, осуществляющих реакцию травления, в результате чего все большая часть поверхности кремния становится пассивной из-за хемосорбированного кислорода. С ростом параметра а3 увеличивается интенсивность хемосорбции кислорода на кремнии и ее действие становится заметным уже при малых добавках 02. При экстремальном значении а3 = 100 параметр $о растет линейно с увеличением доли 02 в смеси и при 50 %-м содержании кислорода практически вся поверхность образца насыщена хемосорбированным кислородом.
Конкуренция процессов хемосорбции атомарных фтора и кислорода на кремнии приводит к гистерезисной зависимости скорости спонтанного травления от концентрации
Рис. 6. Зависимость доли поверхности кремния, покрытой хемосорбированным кислородом, от процентного содержания О2 в СР4/О2 при де = 1: 1 — а8 = 0.3; 2 — а8 = 1; 3 — а8 = 5; 4 — а, = 100
Рис. 7. Зависимость средней скорости травления от средней концентрации фтора на кремнии при различном содержании О2 в СР4/О2. Маркерами на кривых обозначена доля О2 в СР4/О2, растущая в направлении стрелки от 30 до 60% с шагом в 5%: 1 — де = 1, а8 = 0.3; 2 —
де = 1, а = 1; 3 — де = $е(жо2), «« = 0.3; 4 — де = де(хо2), = 1
фтора (рис. 7) [4]. В предыдущих расчетах [1, 3] были получены зависимости формы кривой гистерезиса и положения начала закручивания кривой от параметра Учет зависимости (хо2) не приводит к качественным изменениям полученных ранее
■■и 0.12
0.11
0.1
0.09
0.08
0.07
0.06
и 05 >
- с : / 5 ©
. <2> - / / с )
15 [ Г в \
Я [Г г ) \ \ —□— 1 - -О- - 2
1 о \ \
: ш \ ч Р
-1 1 1 II III с к , Ип ^ и ч! г е г* Г* J1 1 1 1
10 20 30 40 50 60 70 80 90
%,о2
Рис. 8. Зависимость индекса неоднородности травления от процентного содержания О2 в ОР4/С2 при ав = 0.3: 1 - 0е = 1; 2 - 0е = $е(жо2)
результатов. Понижение электронной плотности в реакторе вызывает одновременное пропорциональное уменьшение амплитуд средних скорости травления и концентрации фтора на образце, что ведет к стягиванию кривой гистерезиса к началу координат без заметного изменения ее формы. Количественно это подтверждает таблица относительных отклонений средних скорости травления и концентрации фтора. Положение начала закручивания кривой относительно процентного содержания кислорода в смеси сохраняется.
На рис. 8 представлена зависимость индекса неоднородности травления от процентного содержания 02 в СГ4/02. Индекс неоднородности травления, рассчитанный при постоянной средней электронной плотности в реакторе = 1, находится в пределах 5.5-10 %. Это связано с тем, что в области г = 23 — 27.5 см с ростом добавки кислорода в пределах 0-25 % наблюдается наибольшие изменения скорости травления. На практике образцы обрабатываются при небольших добавках 02 < 30 % и располагаются на некотором расстоянии от внешнего и внутреннего краев нижнего электрода, поэтому резким падением скорости травления в этой зоне можно пренебречь. В оставшейся области обеспечивается травление, близкое к однородному. Индекс неоднородности травления, рассчитанный с учетом зависимости $е(жо2), оказывается выше на 1-5%, чем в реакторе без учета данного эффекта. Неоднородность травления составляет более 10 % при любом содержании кислорода. Наибольшие изменения скорости травления при добавках кислорода до 25 % локализуются на внешнем краю образца. При добавках кислорода, больших 25 %, неравномерное травление имеет место в достаточно широкой области пластины, что не соответствует требованиям к обработке образцов.
Заключение
Выполнены расчеты процесса травления кремния в плазме CF4/O2 с учетом экспериментальных данных по снижению электронной плотности в зависимости от концентрации 02. Сравнение полученных результатов с результатами предыдущих расчетов [1, 3], где этот эффект не учитывался, позволяет сделать следующие выводы.
1. Снижение электронной плотности в ВЧ-разряде, вызванное добавкой кислорода, приводит к падению скорости спонтанного травления в пределах до 30 %. При этом снижение в диапазоне оптимальных скоростей травления, достигаемых при концентрациях кислорода 25-40%, превышает 20%. Тем не менее основной эффект повышения скорости травления в газовой смеси CF4/O2 по сравнению с чистым тетрафторметаном по порядку величины сохраняет свое значение.
2. Понижение электронной плотности в реакторе приводит к одновременному пропорциональному уменьшению амплитуд средних скорости травления и концентрации фтора на образце с сохранением формы кривой гистерезиса. Качественные изменения в сценарии процесса травления отсутствуют.
3. Для исследования физических механизмов влияния электронной компоненты на процессы в плазмохимическом реакторе необходимо математическое моделирование ВЧ-разряда.
Список литературы
[1] Grigoryey Yu.N., Gorobchuk A.G. Numerical Simulation of plasma-chemical reactors // Comput. Technol. 2003. Vol. 8. Special Issue. Pt. 2. P. 53-73.
[2] Григорьев Ю.Н., Горовчук А.Г. Оптимизация состава смеси для травления Si в CF4/O2 // Тр. Междунар. конф. RDAMM-2001. Новосибирск, 2001. Ч. 2. C. 217-224. (Электронная публикация).
[3] Григорьев Ю.Н., Горовчук А.Г. Особенности интенсификации травления кремния в плазме CF4/O2 // Микроэлектроника. 2007. Т. 36, № 5. (в печати).
[4] Mogab C.J., Adams A.C., Flamm D.L. Plasma etching of Si and SiO2 — The effect of oxygen additions to CF4 plasmas // J. of Applied Physics. 1978. Vol. 49, N 7. P. 3796-3803.
[5] Coburn J.W., Chen M. Optical emission spectroscopy of reactive plasmas: A method for correlating emission intensities to reactive particle density // J. of Appl. Phys. 1980. Vol. 51, N 6. P. 3134-3136.
[6] d'Agostino R., Cramarossa F., De Benedicts S., Ferraro G. Spectroscopic diagnostics of CF4 —O2 plasmas during Si and SiO2 etching processes // J. of Appl. Phys. 1980. Vol. 52, N 3. P. 1259-1265.
[7] Дерюгин А.А., Словецкий Д.И. Моделирование механизма химических реакций в тлеющем разряде в смесях тетрафторметана с кислородом // Химия высоких энергий. 1983. Т. 17, № 4. С. 358-367.
[8] Venkatesan S.P., Trachtenberg I., Edgar T.F. Modeling of silicon etching in CF4/O2 and CF4/H2 plasmas // J. of the Electrochemical Society. 1990. Vol. 137, N 7. P. 2280-2290.
[9] Plumb I.C., Ryan K.R. A model of the chemical processes occurring in CF4/O2 discharges used in plasma etching // Plasma Chemistry and Plasma Processing. 1986. Vol. 6, N 3. P. 205-230.
[10] Park Syng-Kyu, Economou D.J. A mathematical model for etching of silicon using CF4 in radial flow plasma reactor // J. of the Electrochemical society. 1991. Vol. 138, N 5. P. 1499-1508.
[11] Schoenborn Ph., Patrick R., Baltes H.P. Numerical simulation of a CF4/O2 plasma and correlation with spectroscopic and etch rate data // J. of the Electrochemical Society. 1989. Vol. 136, N 1. P. 199-205.
[12] Kopalidis P.M., Jorine J. Modeling and experimental studies of a reactive ion etcher using SF6/O2 chemistry // J. of the Electrochemical society. 1993. Vol. 140, N 10. P. 3037-3045.
[13] Ферцигер Дж., Капер К. Математическая теория процессов переноса в газах. М.: Мир, 1976.
[14] Reid R.C., Prausnitz J.M., Poling B.E. The Properties of Gases and Liquids. N.Y.: McGraw-Hill, Inc. 1986.
[15] Григорьев Ю.Н., Горовчук А.Г. Численная оптимизация планарных реакторов индивидуального плазмохимического травления // Поверхность. 1996. № 2. С. 47-63.
[16] Dalvie M., Jensen K.F., Graves D.B. Modelling of reactors for plasma processing. I. Silicon etching by CF4 in radial flow reactor // Chemical Engineering Sci. 1986. Vol. 41. P. 653-660.
Поступила в редакцию 14 июня 2007 г., в переработанном виде — 30 июля 2007 г.