ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА
2024 • ФИЗИКА • Вып. 2
УДК 538.971, 533.924 PACS 52.40.Hf
Исследование процесса вторичной конденсации LiF в ходе реактивного ионного травления тонкопленочного ниобата лития во фторсодержащей плазме
А. А. Козлов ^, У. О. Салгаева, И.В. Петухов, В. А. Журавлев, А. Б. Волынцев
Пермский государственный национальный исследовательский университет, Пермь, Россия t andreykozlov10031994@gmail.com
Работа посвящена исследованию процесса вторичной конденсации LiF на поверхности тонкопленочного ниобата лития (НЛ) X-среза в ходе реактивного ионного травления (РИТ) в комбинированной плазме ICP/CCP разрядов в газовой смеси Ar + SFe. Процесс РИТ НЛ может быть представлен как топохимическая реакция (ТХР), в ходе которой происходит процесс вторичной конденсации LiF из пересыщенной паровой фазы на поверхность НЛ, описываемый уравнением Колмогорова-Аврами-Ерофеева. Показано, что процесс конденсации характеризуется индукционным периодом, который зависит от скорости РИТ и дефектности кристаллической структуры тонкопленочного НЛ. Для исследования процесса конденсации LiF использовались методы оптической спектральной рефлектометрии, рентгеновской энергодисперсионной спектроскопии (EDS), рентгеновской фотоэлектронной спектроскопии (XPS), сканирующей электронной микроскопии (СЭМ). Для изучения плотности дислокаций в тонкопленочном НЛ использовались методы ямок травления и рентгенографический метод. Впервые показано, что длительность индукционного периода ТХР на поверхности НЛ обратно пропорциональна скорости РИТ. Также процесс РИТ характеризуется критическими скоростями травления, выше которых наблюдается рост тонкой пленки LiF по механизму двумерного зародышеобразования Фольмера-Вебера. Экспериментально продемонстрировано, что критические скорости уменьшаются с ростом давления процесса. Впервые была экспериментально оценена плотность дислокаций в тонкопленочном НЛ методом ямок травления. Проведено экспериментальное сравнение индукционных периодов ТХР на поверхности объемного и тонкопленочного НЛ. Показано, что повышенная плотность дислокаций в тонкопленочном НЛ значительно сокращает длительность индукционного периода. Также в работе определено влияние модификации структуры НЛ при помощи протонного обмена на кинетику РИТ. Продемонстрировано, что замещение ионов Li+ на H+ приводит к уменьшению степени пересыщения паровой фазы LiF над поверхностью НЛ и, как следствие, к значительному увеличению индукционного периода ТХР и увеличению скоростей РИТ в 5-9 раз по сравнению с не модифицированным НЛ.
Ключевые слова: тонкопленочный ниобат лития; реактивное ионное травление; индукционный период; плотность дислокаций; протонный обмен
Поступила в редакцию 23.04.2024; после рецензии 03.05.2024; принята к опубликованию 03.05.2024
The investigation of the LiF redeposition process during reactive ion etching of thin-film lithium niobate with fluorine plasma
© Иванов И. И., Петров П. П., Сидоров С. С. и др., 2021
распространяется на условиях лицензии Creative Commons Attribution 4.0 International (CC BY 4.0).
A. A. Kozlov ^, U. O. Salgaeva, I. V. Petukhov, V. A. Zhuravlev,
A. B. Volyntsev
Perm State University, Perm, Russia t andreykozlov10031994@gmail.com
This paper investigates LiF redeposition during reactive ion etching (RIE) of X-cut thin-film lithium niobate (TFLN) with Ar/SF6 ICP/CCP plasma. RIE of TFLN is a topochemical reaction (TCR), described by the Kolmogorov-Avrami-Erofeev equation, during which there occurs LiF redeposition from the supersaturated vapor phase onto the sample surface. It is shown that the redeposition process is characterized by an induction period, which depends on the etch rate and the TFLN crystal structure imperfections. To study the LiF redeposition process, we used spectroscopic reflectometry, energy-dispersive X-ray spectroscopy (EDS), X-ray photoelectron spectroscopy (XPS), and scanning electron microscopy (SEM). The dislocation density of TFLN was investigated with the help of etch pit method and X-ray diffraction. This study is the first to show that the duration of the induction period of TCR on the TFLN surface is inversely related to the etch rate. The RIE process is characterized by critical etch rates, above which a thin LiF film grows according to the two-dimensional Volmer-Weber mode. It was experimentally shown that the critical etch rates decrease with increasing the process pressure. For the first time, the dislocation density of TFLN was experimentally evaluated by means of the etch pit method. An experimental comparison of the induction periods of TCR on the surface of bulk and TFLN was carried out. It has been revealed that the increased dislocation density of TFLN significantly reduces the duration of the induction period. The impact of the modification of the TFLN structure with proton exchange on the kinetics of RIE was studied. It was demonstrated that the replacement of Li+ with H+ leads to decrease in the degree of LiF vapor saturation above the TFLN surface and, as a consequence, to a significant increase in the induction period of TCR and an increase in the etch rates by 5-9 times compared to unmodified TFLN.
Keywords: thin-film lithium niobate; reactive ion etching; induction period of a topochemical reaction; dislocation density; proton exchange
Received 23 April 2024; revised 03 May 2024; accepted 03 May 2024 doi: 10.17072/1994-3598-2024-2-33-55
1. Введение
Ниобат лития НЛ) уже более 50 лет ис-
пользуется как материал-платформа для изготовления интегрально-оптических (ИО) устройств ра-диофотоники [1]. Использование диффузионных технологий для формирования ИО волноводов на НЛ накладывает ряд ограничений на топологию и параметры создаваемых элементов [2], например:
1. Из-за низкого контраста показателя преломления (ПП) (типичная разность между ПП подложки и волновода не превышает величину Ап — 0.02), поле моды в волноводе имеет размеры 8-9 мкм. В связи с этим при создании электрооптических модуляторов зазор между электродами (в случае использования Х-среза НЛ) составляет 1820 мкм, что, как следствие, приводит к необходимости применения высоких (7-9 В) напряжений для управления характеристиками оптического излучения.
2. Типичные критические радиусы кривизны волноводов, изготовленных по диффузионным технологиям, составляют 1-3 см, что не позволяет изготавливать миниатюрные ИО волноводные элементы, например, компактные кольцевые резонаторы.
С помощью технологии управляемого скола ^таг^и) изготавливаются пластины тонкопленочного НЛ. На основе тонкопленочного НЛ могут быть изготовлены волноводы с контрастом ПП Ап — 0.7 [2-4], что позволяет преодолеть ограничения, описанные выше. В данной работе рассматриваются пластины с конфигурацией слоев, представленной на рис 1.
Рис. 1. Схематическое изображение поперечного сечение пластин тонкопленочного НЛ
Ключевым технологическим процессом формирования профиля ИО волновода (рис. 2) на пластинах тонкопленочного НЛ является процесс «сухого» (плазменного) травления по фоторезистивной или металлической маске.
Для проведения данного процесса наиболее часто используются методы: ионного травления в плазме инертных газов (Аг, Хе) [5-7] и реактивное
ионное травление (РИТ) в плазме газовой смеси реакционного (SF6, CF4, ВС13) и инертного газов [810].
В процессе ионного травления удаление материала, не закрытого маской, происходит под воздействием потока бомбардирующих поверхность вы-сокоэнергетичных ионов инертных газов. В случае РИТ материал подвергается одновременному воздействию как химически активных частиц (ХАЧ) реакционного газа, генерирующихся в плазменном разряде, так и потоку бомбардирующих поверхность высокоэнергетичных ионов.
Рис. 2. Поперечное сечение волновода на
тонкопленочном НЛ
Процесс РИТ позволяет формировать волноводы с углом наклона боковой стенки 65-80 , низкой шероховатостью поверхностей (единицы нанометров) и характеризуется относительно высокой селективностью процесса (скорость травления удаляемого материала значительно выше скорости травления маски). Далее в данной работе будет рассматриваться только процесс РИТ тонкопленочного НЛ. Аналогично публикациям [8, 11] акцент будет сделан на процесс РИТ в смеси газов Аг + SF6 как одной из наиболее часто используемых для данного процесса.
Как было показано в работах В. В. Гуляева [ 12, 13], процесс травления НЛ во фторсодержащей плазме характеризуется образованием на поверхности пластины нелетучего продукта реакции LiF. Гуляев В. В. справедливо классифицирует такой тип реакций как топохимическую (ТХР). ТХР - это тип твердофазных реакций, который протекает на границе раздела «твердый продукт реакции - твердый реагент». Однако, в своих работах В. В. Гуляев не учитывает влияние реальной структуры НЛ (в частности, неравновесных дефектов, например, дислокаций) на протекание такого рода реакций.
В ряде зарубежных работ также отмечается, что в ходе РИТ и плазмохимического травления (ПХТ) НЛ во фторсодержащей плазме происходит процесс вторичной конденсации (переосаждения, с англ. redeposition) LiF на поверхности, который приводит к повышению шероховатости поверхности НЛ и, как следствие, к росту оптических потерь в волноводах, изготовленных в ходе такого процесса РИТ, до значений вплоть до 60 дБ/см [2,6,8,9]. Однако на момент написания данной работы не было найдено исследований, нацеленных на изучение процесса вторичной конденсации LiF.
Целью данной работы является исследование процесса вторичной конденсации LiF в ходе реактивного ионного травления тонкопленочного НЛ в плазме газовой смеси Аг + SF6.
2. Влияние реальной структуры кристалла на особенности кинетики реактивного ионного травления тонкопленочного LiNbOз
2.1. Особенности ТХР, протекающей в ходе РИТ тонкопленочного НЛ во фторсодержащей плазме
В. В. Гуляевым было показано, что плазмохими-ческое травление (ПХТ) НЛ в плазме фторсодержа-щих газов является ТХР, которая протекает на границе раздела твердый продукт реакции (LiF) -твердый реагент (УКЬОз) [12, 13].
Классический механизм ТХР включает в себя несколько стадий протекания реакции, которые показаны на интегральной кривой преобразования п (соответствует зависимости отношения суммарной площади областей, на которых протекла реакция, к площади непрореагировавшей фазы во времени; здесь и далее в работе рассматривается только 2D случай), рис. 3. Скорость травления - это первая производная dп/dt.
Рис. 3. Интегральная кривая преобразования [14]
Участок кривой АВ - это индукционный период или период скрытого протекания реакции, который предшествует стадии роста твердого продукта реакции на поверхности. Инструментально преобразование поверхности реагента в продукт реакции еще не фиксируется. На данной стадии происходит образование ядер-зародышей на поверхности НЛ.
Участки кривой ВС, CD, DE отвечают стадии роста скорости ТХР и связаны с ростом уже образовавшихся ядер продукта и их последующим перекрытием. Считается, что ускорение реакции связано с тем, что происходит когерентный рост продукта реакции с периодом решетки отличным от периода решетки реагента. Это приводит к самоускорению реакции за счет возникновения упругих напряжений в месте сопряжения кристаллических
решеток (автокаталитический процесс, параметры кристаллических решеток: аьш = 0.403 нм, аш = 0.515 нм). Начиная с участка EF происходит торможение реакции с полным перекрытием поверхности реагента продуктом реакции, дальнейшая скорость реакции лимитируется диффузионными процессами через слой продукта реакции.
Необходимо заметить, что процесс ПХТ НЛ, показанный В. В. Гуляевым в работах [13, 14], характеризуется высокими давлениями - 0.1-1 мбар, образец находится под плавающим потенциалом плазмы [15]. Как следствие, для частиц плазмы характерны низкие энергии (не превышают 5-10 эВ). Такие энергии недостаточны для эффективного распыления нелетучего LiF с поверхности (пороговые энергии ионного распыления LiF составляют величину 20 эВ для Аг+, типовые значения коэффициента распыления Y - 1 молекула LiF на 1 ион Аг+ с энергией 100 эВ [17]). Также процесс травления характеризуется высокими энергиями активации и с формальной точки зрения является классической ТХР.
В процессе РИТ в плазме газовой смеси Аг + SF6, как было показано в работе [18], происходит стимулированное ионной бомбардировкой протекание химической реакции (энергиями Аг+ нельзя пренебречь, ЕАг+ в зависимости от режима РИТ могут достигать значений 250 эВ). В ходе реакции образуются следующие продукты реакции: LiF, NbF5, OF2, О2. Продукты реакции десорбируются с поверхности как самопроизвольно ввиду низкой температуры испарения при давлении процесса Р = 0.005 мбар (характерно для NbF5, OF2, О2), так и вследствие ионного распыления (LiF). Схема процесса РИТ НЛ показана на рис. 4, а). Продукты реакции откачиваются вакуумной системой, и при увеличении расстояния г от пластины концентрация продуктов реакции падает (рис. 4, б), в процессе РИТ устанавливается некоторый динамический профиль концентрации с максимумом у поверхности пластины.
OF2 имеют температуру кипения ниже рабочей температуры процесса (Трав. - 200 °С, Тшп^5 = -21 °С, Ткип^2 = - 215 °С [18], пары NbF5, OF2 ненасыщенные и испаряются с поверхности самопроизвольно), LiF при обозначенных параметрах процесса находится над поверхностью LiNbOз в пересыщенном состоянии (Тшпьш = 558 °С> Трав) [18].
Степень пересыщения LiF может быть выражена формулой:
О = ^парц/^ (1)
где Рпарц - парциальное давление LiF, Рб - давление насыщенного пара при данной температуре.
Здесь и далее по тексту под Трав подразумевается средняя температура, которую достигает образец в ходе нагрева плазмой при выходе на стационарную температуру по истечении 3 мин; в момент инициации плазмы образец имеет температуру Тнач = 20 °С, равную температуре хладагента, охлаждающего столик-держатель образца, и выходит на Тстац спустя 3 мин [18].
Для грубой оценки парциального давления LiF воспользуемся уравнением Менделеева-Клапейрона (более строгое рассмотрение требует использования уравнения Ван-дер-Ваальса). Тогда можно записать выражение для парциального давления LiF (без учета работы системы вакуумной откачки):
(2)
Рпарц =
тшКГ
Рис. 4. Схема процесса РИТ НЛ (а) и концентрационный профиль паров LiF вблизи поверхности (б)
Вблизи поверхности каждый продукт реакции имеет свое парциальное давление. При рабочем давлении Р = 0.005 мбар продукты реакции NbF5,
парц у •
г-ЬьР ^кам
где - это масса LiF, которая образуется в ходе РИТ, ^ыр - молярная масса LiF, Т - температура процесса, Укам - объем вакуумной камеры. При этом масса LiF может быть выражена как:
тш = тшьоз% = PыNboз^ЫNboз^, (3)
с учетом площади образца и толщины пленки НЛ, удаленной в процессе травления, запишем выражение для шыр :
тш = (4)
где рыв- - плотность LiF, УРИТ - скорость РИТ, £ -площадь образца, который подвергается травлению; УРИТ t - толщина пленки НЛ, которая удалена в процессе травления за время травления V, £ - коэффициент, определяющий массу продукта реакции LiF, образующегося по схеме LiNbOз+8F = LiF +NbF5+OF2+O2, при реакции НЛ объемом Уыыьоз = БУрит t.
В итоге, наблюдаем линейную зависимость степени пересыщения от УРИТ:
= Й^Рьшьоз^ит^ (5)
Данное выражение дает грубую оценку и не учитывает работу системы вакуумной откачки.
Для грубой численной оценки Рпарц возьмем следующие значения констант: = 25.9 • 10-3 кг/моль, Т = 473 К, Укам = 0.064 м3, римо = 4700 кг/м3; Б = 2.25 • 10-4 м-2, Урит = 5 • 10-10 м/с, t = 120 с, £ = 0.087.
Тогда Рпарцх^ = 1.3 • 10-2 Па или 1.3 • 10-4 мбар,
для сравнения рабочее давление Р = 5 • 10-3 мбар. Полученное значение по порядку величины согласуется с типовыми значения парциальных давлений, наблюдающихся при вакуумном напылении материалов [19, 20].
Поскольку рабочее давление процесса РИТ Р = 0.005 мбар много ниже давления тройной точки для ЫБ (~103 Па), то в ходе РИТ происходит процесс конденсации ЫБ на поверхности НЛ. Впервые экспериментально это продемонстрировал Н. Nagata в работе [21] при травлении НЛ в плазме
Таким образом, в процессе РИТ протекает не классическая ТХР с формальной точки зрения, а процесс вторичной гетерогенной кристаллизации ЫБ, который переосаждается из пересыщенной газовой фазы на поверхность НЛ. Пересыщение газовой фазы ЫБ является движущей силой процесса конденсации.
Как показано в работе [18], процесс РИТ НЛ характеризуется схожими кинетическими закономерностями, которые присущи классическим ТХР (например, ускорение скорости травления при росте ЫБ на поверхности - автокаталилический процесс).
Рост тонкой пленки ЫБ по поверхности НЛ может быть описан классической для ТХР моделью -
геометрико-вероятностной моделью изотермической кристаллизации Колмогорова-Аврами-Ерофе-ева (КАЕ). Степень превращения во времени описывается следующим уравнением (6) [16]:
^ = 1 — е
-лпу2ь2
(6)
где V - это скорость роста новой фазы LiF на поверхности, X - среднее число зародышей на единицу площади.
При выводе уравнения Колмогоровым сделано несколько предположений, важнейшие из которых: 1) рост новой фазы происходит на фиксированных центрах зародышеобразования (рассматривается для простоты двумерный случай); 2) рост новой фазы происходит равномерно во все стороны с постоянной скоростью (характерно для изотермического протекания процесса); 3) процесс образования зародыша на потенциальном активном центре описывается распределением Пуассона, возникновение зародышей на поверхности случайно.
Используемая в данной работе экспериментальная установка Этна-100 имеет водоохлаждаемый столик без гелиевого поддува, на котором располагается образец в ходе травления. В процессе РИТ условия не могут считаться изотермическими ввиду нагрева образца плазмой и низкой эффективности теплоотвода. Более строго процесс неизотермической кристаллизации описывается моделью Одзавы и моделью Мо, однако данные модели требуют знания завимости V = Т) [22]. В то же время ТРИТ является функцией мощностей, давления и потока газов. Переход к моделям неизотермической кристаллизации с практической точки является избыточным усложнением задачи.
В наблюдаемых условиях эксперимента рост
температуры со временем в процессе РИТ (типовой температурный график показан в работе [18]) приводит к изменению пересыщения. При этом в зависимости а = Pпарц/Ps знаменатель растет намного быстрее, чем числитель. Так для Т = 20 °С ~ 10-62 Па, а при Т = 300 °С Р^ы* ~ 10-33 Па [23]. Степень пересыщения айв- уменьшается с момента инициации плазмы в ходе роста температуры.
Таким образом, уравнение КАЕ является оценкой «сверху». Это означает, что реально наблюдаемые в эксперименте степени преобразования поверхности Пэксп меньше, чем птеор. С практической точки зрения такая оценка является достаточной.
При использовании процесса РИТ для формирования ИО волноводов ФИС рост LiF на поверхности является нежелательным из-за повышения степени шероховатости поверхностей. Поэтому необходимо понимание механизмов роста LiF на поверхности НЛ. В зависимости от степени пересыщения паровой фазы рост кристаллов происходит по трем различным механизмам, приведенным ниже (ранжирование по степени понижения пересыщения [24]):
1. путем образования и роста двумерных зародышей (механизм Фольмера-Вебера);
2. дислокационный рост по механизму Бартона-Кабреры-Франка (БКФ) [25, 26];
3. путем нормального послойного роста (характерно для крайне низких пересыщений).
Степень пересыщения паровой фазы ЫБ в процессе РИТ достаточно велика (с не ниже 1030) для осуществления роста ЫБ на поверхности по механизму двумерного зародышеобразования.
При росте пленки из паровой фазы продолжительность индукционного периода зачастую является наиболее длительной стадией [27]. Поэтому для обеспечения контролируемости процесса РИТ требуется понимание механизмов, которые влияют на длительность индукционного периода роста ОБ на поверхности.
Если полагаем, что Пинд - это степень преобразования поверхности НЛ, которая соответствует окончанию индукционного периода, по завершению которой ОБ на поверхности детектируется инструментально, то соответствующая ей длительность индукционного периода из уравнения (6) может быть выражена следующим образом:
т =
'-инд
N
1п(1 — Линд 1
\
(7)
Фактически, под понятием индукционный период в рамках данной работы подразумевается время задержки нуклеации 1нукл (время, которое требуется для роста дозародышевого образования ада-томов до размера термодинамически устойчивого, сверхкритического) плюс время, в ходе которого зародыш достигает размеров, детектируемых инструментально 4нст. По уравнению (7) видно, что длительность индукционного периода уменьшается с увеличением скорости роста ОБ на поверхности.
Поскольку мы полагали, что Тинд = tн:
то
справедливо рассматривать скорость роста LiF на поверхности как зависящую от частоты зародыше-образования и скорости роста зародышей-островков до размеров, детектируемых инструментально.
Частота зародышеобразования растет с ростом
■ф2
пересыщения как I = "ж, где коэффициент у1 зависит от свойств среды, скорости доставки частиц к поверхности зародыша, а у2 - от температуры и поверхностного натяжения.
Скорость роста кристалла есть линейная функция от пересыщения при двумерном механизме зародышеобразования: VLiF = уз&ив, при этом константа у3 также зависит от материалов и условий роста [28].
В процессе РИТ, как было показано выше, степень пересыщения газовой фазы оив- зависит от скорости РИТ Урит (скорости выхода продуктов реакции) линейно.
Тогда в соответствии с показанным выше из уравнения (7) получим выражение, демонстрирующее, что продолжительность индукционного периода образования LiF на поверхности обратно пропорциональна скорости РИТ:
т ~!/и (8)
нд Чит'
Таким образом, в процессе РИТ НЛ должно наблюдаться сокращение индукционного периода образования LiF на поверхности с увеличением скорости РИТ.
Увеличение скорости РИТ приводит к повышению концентрации продукта реакции LiF над пластиной НЛ и, как следствие, к повышению степени пересыщения. Сокращение индукционного периода с ростом пересыщения является общей закономерностью зародышевого механизма роста кристаллов.
В ходе работы данная зависимость также была проверена экспериментально.
2.2. Влияние плотности неравновесных дефектов на кинетику ТХР
При рассмотрении влияния реальной структуры кристалла на кинетику ТХР ограничимся описанием только неравновесных дефектов (собственные точечные дефекты наравне с температурой и давлением могут считаться параметрами состояния кристалла [29]).
Под кинетикой в данной работе будем понимать не только механизм изменения скорости протекания ТХР, но и особенности перехода «индукционный период - период роста скорости ТХР» (в частности, длительность индукционного периода, по истечении которого на поверхности детектируется рост продукта реакции LiF).
Как было показано выше (уравнение (7)), продолжительность индукционного периода уменьшается при росте числа зародышей на единицу площади по закону 1/Х. При этом необходимо заметить,
что для высоких степеней пересыщения поверхностная плотность центров зародышеобразования Я Ф Nduсл (Nдисл плотность дислокаций), поскольку возможно появление зародышей надкритического размера флуктуационным путем. Так, согласно работе [27], типовая X ~ 109 см-2 (что значительно превышает типовые значения Nduai, характерные для кристаллов полупроводниковой промышленности).
Дислокации вносят вклад в рост кристалла как центры образования зародышей (экспериментально показано в [30]). В зародышевом механизме роста кристаллов считается, что области с выходом дислокаций снижают энергетический порог образования зародышей [31]).
Академик В. В. Болдырев в своих трудах развил идею, связанную с тем, что такие неравновесные дефекты, как дислокации в кристаллах, играют ключевую роль при зародышеобразовании в ходе ТХР [32-34]. Такое явление В. В. Болдырев назвал автолокализацией - рост зародыша на потенциально активном центре порождает дислокации несоответствия, что приводит к самоускорению реакции (автокатализ, химическая реакция протекает с положительной обратной связью).
Ранее в работах [34, 35] было показано, что радиационное воздействие на кристалл (поток электронов, гамма-квантов, ионов) - это один из способов управляемого изменения кинетики ТХР; P. van der Put в монографии [36] продемонстрировал, что увеличение плотности дислокаций в результате радиационного воздействия приводит к уменьшению длительности индукционного периода ТХР.
Показанные выше закономерности подчеркивают важность изучения реальной структуры тонкопленочного НЛ (в частности, плотности дислокаций) для контролируемого протекания процесса РИТ. В данной работе авторы намерены провести экспериментальное сравнение длительности индукционного периода ТХР объемного и тонкопленочного НЛ.
2.3. Изменение реальной структуры тонкопленочного НЛ методом протонного обмена как способ контроля кинетики ТХР
В работах [9, 37-39] показана возможность уменьшения степени переосаждения LiF и повышения скорости плазмохимического травления НЛ во фторсодержащей плазме. Это достигается за счет проведения предварительного протонного обмена (ПО) и снижения концентрации Li в кристаллической решетке НЛ в ходе ПО.
Увеличение скорости травления НЛ обеспечивает повышение селективности процесса РИТ (увеличение отношения скорости травления НЛ к скорости травления маски). Такое повышение селективности позволяет снизить требования к
маске для травления (выполнить переход от металлических масок к полимерным, уменьшить толщину маски) и тем самым улучшить качество формируемых волноводных структур.
В процессе ПО происходит замена ионов Ы+ в кристаллической решетке НЛ на ионы Н+ в соответствии со следующей схемой:
LiNbOo +хН+ = + xLi+
(9)
Процесс ПО является не чистой заменой ионов
на Н+, а сопровождается внедрением Н+ в меж-доузельные позиции. Так, после ПО в адипиновой кислоте (Тпо = 213 С, t = 60 мин) отношение концентраций Ы+:Н+ = 1:1.16 [40].
Формирование дефектов по типу междоузель-ных ионов Н+ с концентрацией пм.у. и дефектов замещения с концентрацией пз. ассоциируется с увеличением свободной энергии кристалла на величину АВ = АВм.у + АВ3.
Известно, что мерой химической активности твердых тел является избыток их свободной энергии по отношению к нормальному состоянию (характеризующемуся равновесной концентрацией дефектов) за счет накопления неравновесных дефектов в кристаллической решетке [41] Таким образом, создание дефектов в кристаллической решетке НЛ должно сказываться на скоростях РИТ.
Однако остаются малоизученными вопросы: 1) как зависит скорость РИТ от стехиометрии кристалла Ы№03 (степени замещения), в том числе для тонкопленочного НЛ; 2) изменится ли скорость РИТ с расстоянием от исходной поверхности вглубь пленки? Ввиду диффузионной природы процесса ПО встает вопрос о возможности травления модифицированного протонным обменом слоя НЛ с одинаковой скоростью по глубине. Этот вопрос крайне важен с практической точки зрения при разработке технологического процесса травления тонкопленочного НЛ для создания ФИС. Авторы статьи планируют провести ряд экспериментов, позволяющих исследовать описанные выше особенности.
Также необходимо отметить, что так как процесс ПО приводит к снижению концентрации О в кристаллической решетке НЛ, это, в свою очередь, может сказаться на степени пересыщения паров ЫБ над поверхностью пластины в процессе РИТ. Можно предположить, что ПО может приводить к увеличению длительности индукционного периода тинд. В ходе работы будут выполнены эксперименты, направленные на сравнение Тинд исходных образцов тонкопленочного НЛ и модифицированных ПО.
3. Дефекты структуры тонкопленочного НЛ
3.1.Описание технологии управляемого скола (БтаНСШ)
Как было показано в разделе 2.2, реальная структура кристалла НЛ оказывает существенное влияние на протекание ТХР. Поэтому важно проанализировать технологию изготовления исследуемых подложек тонкопленочного НЛ.
На данный момент пластины тонкопленочного НЛ изготавливаются по двум технологиям: технологии управляемого скола (БтаНСШ) [2-4, 42, 43] и технологии бондинга пластины НЛ на подложку из другого материала с последующим хемомеханиче-ским утонением НЛ до требуемых толщин [43, 44]. Объектом исследования в данной работе является тонкопленочный НЛ Х-среза, изготовленный по технологии управляемого скола (NanoLN, Китай).
На рис. 5 приведена схема ключевых этапов изготовления пластин тонкопленочного НЛ с помощью технологии управляемого скола.
На первом этапе (рис. 5, а) осуществляется процесс ионной имплантации ионов Не+ в пластину-донор объемного НЛ. При ионной имплантации остановка Не+ происходит на некоторой проективной глубине, определяемой начальной энергией ионов. Вторым важным параметром является флюенс ионов - в рассматриваемой технологии Ф — 1016 -1017 см-2.
Рис. 5. Схема этапов технологии изготовления пластин тонкопленочного НЛ методом управляемого скола
На рис. 6 приведен профиль концентрации дефектов в НЛ после процесса ионной имплантации при типовых режимах, используемых при изготовлении пластин тонкопленочного НЛ.
В слое НЛ, в котором происходило торможение ускоренных ионов при ионной имплантации, можно условно выделить три области, каждая из которых характеризуется различной степенью повреждения кристаллической структуры за счет упругих соударений [45]:
1. Приповерхностный слой с относительной концентрацией дефектов до —0.7 на Х-срезе НЛ (значительно выше, чем для 2-среза). В данном слое преимущественным типом дефектов являются вакансии типа №□ (перестановка № в вакантные октаэдрические пустоты).
2. Переходная область, в которой формируются 2D и 3D кластеры из точечных дефектов.
3. Аморфная область с полным нарушением периодичности расположения атомов в кристаллической решетке НЛ.
Рис. 6. Концентрационные профили дефектов в НЛ для ЕНе = 325 кэВ и флюенсе Ф = 4 х 1016 см-2 [4]
На втором этапе (рис. 5, б) изготовления тонкопленочного НЛ выполняется подготовка пластины-акцептора - кремния (реже НЛ) с изолирующим слоем SiO2 на поверхности с толщиной от 2 мкм.
Далее (рис. 5, в) производится бондинг пластины-донора к пластине-акцептору при комнатной температуре. Для обеспечения адгезии может проводиться предварительная функционализация поверхностей ОН-группами [46].
На следующем этапе (рис. 5, г) осуществляется отжиг пластин (Т = 150-250 °С), в ходе которого внутри аморфной области образуются пузыри-скопления гелия. В результате чего происходит откалывание остаточной части пластины-донора, а тонкая пленка НЛ остается на поверхности пластины-акцептора.
На финальном этапе (рис. 5, д) выполняется отжиг при Т = 400-500 °С в течение нескольких часов для аннигиляции точечных дефектов и восстановления кристаллической структуры в пленке НЛ. После отжига производится хемомеханическая полировка поверхности, в ходе которой полностью удаляется аморфный слой с поверхности, а также большая часть переходного слоя [4].
3.2. Механизм увеличения плотности структурных дефектов в тонкопленочном НЛ
Для изготовления пластин тонкопленочного НЛ в качестве пластин-доноров используются пластины объемного НЛ, которые характеризуются собственной плотностью дефектов кристаллической структуры (так, например, плотность дислокаций в объемном НЛ — 104 см-2 [47]).
При изготовлении тонкопленочного НЛ повышается плотность структурных дефектов в тонких пленках НЛ в ходе ионной имплантации и финального отжига.
Так, в работе [48] впервые показано, что после ионной имплантации ускоренных ионов Не+ и последующего отжига НЛ происходит увеличение плотности дислокаций, которое влечет за собой экстремальное изменение ПП материала.
Gischkat Т. с соавторами в работе [49] показали, что после ионной имплантации Не+ в НЛ эффективная аннигиляция точечных дефектов при отжиге наблюдается только при низких относительных концентрациях дефектов (до Потнос —0.4). При высоких концентрациях дефектов (потнос> 0.4) в процессе отжига происходит агломерация вакансий в вакансионные диски с дальнейшим их схлопыва-нием в дислокационные петли. Наиболее отчетливо этот эффект проявляется на Х-срезе НЛ. Более подробное описание механизма образования дислокационных петель после ионной имплантации дано в источниках [50, 51].
При изготовлении тонкопленочного НЛ после отжига методом полировки удаляются аморфная и переходная области, однако, если проанализировать рис. 6, можно отметить, что на поверхности тонкой пленки НЛ после хемомеханической полировки остается, область с потнос> 0.4, в которой потенциально могут образовываться дислокационные петли. Таким образом, можно ожидать, что поверхность тонкопленочного НЛ будет характеризоваться плотностью дислокаций, повышенной относительно пластин объемного НЛ.
В 2023 г. в работах [52, 53] впервые проведено целенаправленное исследование изменения структуры НЛ после ионной имплантации и термической обработки при режимах, используемых в коммерческих технологиях изготовления пластин тонкопленочного НЛ. С помощью методов диффузного рассеяния рентгеновского излучения и просвечивающей электронной микроскопии было показано, что отжиг НЛ после ионной имплантации приводит к восстановлению кристаллической структуры, однако в процессе происходит перестройка кластеров точечных дефектов в дислокационные петли.
Несмотря на наличие экспериментальных свидетельств образования дислокационных петель в ходе отжига НЛ, подвергшегося имплантации ионов Не+, в литературе не приведены количественные оценки или измеренные значения плотности дислокаций в тонкопленочном НЛ для сравнения с данными для объемного НЛ.
На данный момент существует ряд работ, в которых методами дифракционного структурного анализа показано относительное уширение пика кривой качания тонкопленочного НЛ относительно объемного НЛ [43, 54, 55], что может быть косвенным признаком повышенной плотности дислокаций.
Авторы данной работы намерены провести экспериментальное измерение плотности дислокаций на поверхности пластин тонкопленочного НЛ прямым методом - методом ямок травления.
4. Методика эксперимента, описание образцов и оборудования
4.1. Описание образцов
Исходные образцы: 1) тонкопленочный НЛ (X-срез, структура пластин приведена на рис.
1; 2) объемный ниобат лития (Х-срез, CQT, толщина пластин - 1 мм). Пластины для проведения экспериментов были разделены на кусочки площадью 15 х 15 мм2 методом дисковой резки.
4.2. Методика эксперимента по определению зависимости длительности индукционного периода ТХР от скорости травления и плотности дислокаций
Процесс РИТ производился на установке плаз-мохимического осаждения и травления МБТ Этна-100, принципиальная схема которой приведена на рис. 7.
Рис. 7. Схема установки для проведения процесса РИТ
Подробное описание принципа работы установки приведено в работе [18].
Для определения зависимости длительности индукционного периода от скорости травления Тинд = /(Урщ) следующие параметры процесса РИТ поддерживались неизменными: Р = 0.005 мбар, №сер = 80 Вт, <Аг: QsF6 = 38: 2 см3/мин. Для изменения скорости травления производилась вариация WICP в диапазоне от 100 до 600 Вт. Как было показано ранее, зависимость Утрав от WICP линейная [18].
Скорость РИТ вычислялась при помощи выражения (10):
^травл
^нач ^Ki
(10)
где dm и dtroti — это начальная и конечная толщины пленки НЛ, t - время травления.
Измерение толщин пленок НЛ производилось при помощи спектрального рефлектометра Wuhan EOptics SR-C. Толщины пленок определялись с погрешностью 0.2 %, но не менее 2 нм.
При каждом значении мощности WiCp изменялось время травления с шагом 1 мин (на низких мощностях - 100 и 200 Вт - с шагом 4 мин). Время увеличивалось до тех пор, пока не начинали появляться следующие признаки окончания индукционного периода (образцы исследовались по окончании времени травления): 1) на поверхности НЛ методом энергодисперсионной спектроскопии (EDS) детектировался F, что косвенно свидетельствовало о наличии соединения LiF на поверхности; 2) коэффициент достоверности аппроксимации R2 в модели расчета спектрального рефлектометра становился ниже 0.95, что интерпретировалось как необходимость вводить в модель помимо слоев SiO2 и LiNbO3 слой LiF на поверхности; 3) при помощи сканирующего электронного микроскопа (СЭМ) на поверхности определялись области с пористой морфологией с типичными размерами пор от единиц до десятков нанометров в диаметре. После окончания индукционного периода все указанные признаки наблюдались одновременно. Более подробно методика определения перехода «индукционный период - период роста скорости ТХР» описана в работе [18]. Каждому исследуемому режиму соответствовал один образец.
Для сравнения длительностей индукционного периода ТХР объемного и тонкопленочного НЛ производилось травление при следующем режиме РИТ: P = 0.005 мбар, Wicp = 500 Вт, WCCP = 80 Вт, QAr: QSF6 = 38: 2 см3/мин. Время изменялось также с шагом 1 мин.
4.3. Методика определения плотности дислокаций в НЛ
Определение плотности дислокаций в объемном и тонкопленочном НЛ было проведено при помощи двух альтернативных методов: 1) метод ямок травления, 2) рентгенографический метод.
4.3.1. Рентгенографический метод
Рентгенографический метод применялся для качественного сравнения совершенства структуры объемного и тонкопленочного НЛ. Сравнение производилось на основании сопоставления полуширин дифракционных максимумов. Съемка дифракционных кривых производилась на рентгеновском двухкристальном дифрактометре ДРОН-УМ1 (рис. 8).
Ускоряющее напряжение рентгеновской трубки 1 U = 30 кВ, ток через рентгеновскую трубку I = 10 мА. Материал анода Co, Kß- линия, ÄKß = 1.62073 Ä.
травление НЛ. Поэтому был собран экспериментальных стенд для выявления ямок травления в парах (Рис. 9). За идею были взяты работы И. И. Сорокиной и Б. Б. Педько. [57, 58].
Рис. 8. Принципиальная схема дифракто-метра ДРОН-УМ1 (обозначения в тексте)
В качестве кристалла-монохроматора 2 использовался бездислокационный Si, установленный в положение брэгговского максимума, соответствующего Kp-линии от плоскости (111) в первом порядке отражения. Использование такой конфигурации позволяло считать угловую расходимость первичного пучка малой и полагать, что ширина спектральной линии определяется только микродеформациями образца.
После кристалла-монохроматора устанавливались щели 3: 1.0 мм - ограничение вертикального размера пучка, 0.2 мм - ограничение ширины пучка. Исследуемый образец 5 и счетчик 7 располагались на платформе гониометра 4. Перед счетчиком устанавливалась вертикальная щель 6 шириной 0.05 мм. В установке использовалась фокусирующая схема по Брэггу-Брентано.
Интенсивность определялась числом импульсов на счетчике, регистрируемых за 10 с. Для всех образцов снимались кривые первого порядка от плоскости (110). Дифрактограммы снимались при последовательном вращении образца на угол 0 при неподвижном счетчике (кривые качания, ю-кри-вые). Шаг изменения угла - 0.0025
Для сравнения дефектности кривые качания для тонкопленочного и объемного НЛ были построены на одном графике с нормированием на максимум интенсивности. Количественная оценка не производилась. Выделение вкладов физического ушире-ния кривых отражения для субмикронных пленок является не тривиальной задачей.
4.3.2.Метод ямок травления
Метод ямок травления, напротив, применялся для количественной оценки плотности дислокаций на поверхности изучаемого НЛ. Методика определения плотности дислокаций методом ямок травления в смеси KOH +NaOH в объемном НЛ приведена в статье [56].
Данная методика не применима для изучения плотности дислокаций в тонкопленочном НЛ, поскольку пластиной-носителем является Si, который растворяется с высокой скоростью в большинстве неорганических соединений, в которых проводится
Рис. 9. Схема экспериментального стенда для выявления ямок травления в парах
Выявление ямок травления производилось в парах смеси кислот HF: HNO3 = 2:1. Рабочая температура смеси +70 °C. Образец был обращен пленкой НЛ к раствору. Температура раствора контролировалась при помощи бесконтактного пирометра Fluke. Оптимальное расстояние от образца до поверхности раствора — 1 см. Время обработки в парах - 1 ч 30 мин. После обработки в парах ямки травления исследовались при помощи СЭМ Tescan Mira3.
Расчет средней плотности дислокаций по ямкам травления производился по следующей формуле:
N =
^дисл
N,
ср 5 '
(11)
где Ыср - среднее количество ямок в поле СЭМ, Б -площадь поля зрения в см2.
4.3.3. Методика изучения профиля плотности дислокаций в тонкопленочном НЛ по глубине
Как было показано выше, процесс ионной имплантации оказывает существенное влияние на формирование дислокационных петель в ходе изготовления тонкопленочного НЛ [53]. Определение профиля распределения плотности дислокаций по глубине пленки НЛ является самостоятельной актуальной задачей, поскольку плотность дефектов оказывает прямое влияние на контролируемость процесса травления.
Изначально планировалось проводить контроль плотности дислокаций по глубине методом ямок травления. Однако, как показали эксперименты, применение метода ямок травления в парах некорректно использовать при плотности дислокаций ниже ~ 106 см-2. При снижении плотности дислокаций наблюдалось не только увеличение времени их выявления, но и, как следствие, появление большого числа поверхностных дефектов и отслаиваний в пленке НЛ.
В качестве альтернативы был использован косвенный метод, основанный на идее, предложенной в статье [59]. Г. В. Гордеевой и коллегами на примере ионного распыления стали 0Х15Н16М3Б ионами Аг+ с энергией 10 кэВ было экспериментально продемонстрировано, что образующиеся от дислокаций поля упругих напряжений влияют не только на развитие каскада упругих столкновений атомов в тонком поверхностном слое, но и на скорость распыления. По мере увеличения плотности дислокаций происходило уменьшение коэффициента распыления стали.
Так, для изучения изменения плотности дислокаций по глубине была отработана методика, в ходе которой производилось пошаговое удаление слоя за слоем НЛ при помощи метода ионного травления на установке Этна-100 при следующем режиме: Р = 0.005 мбар, WIcp = 500 Вт, Wccp = 100 Вт, <2лг = 20 см3/мин, t = 5 мин.
В экспериментальном режиме ионного травления, описанном выше, ионы Аг+ обладают энергией, не превышающей 0.2 кэВ. Согласно теории распыления твердых тел П. Зигмунда [60], при энергиях падающих ионов менее 1 кэВ нельзя считать полноценным развитие линейного каскада упругих соударений. При энергии ионов 0.2 кэВ можно говорить о так называемом режиме единично выбитых с поверхности атомов (е.в.а), схема которого приведена на рис. 10.
\ Е > Ц ^У
о о \ о о о
' ' \ о о '.о / о* о
» / * /
о о о о о о о о\ о о
Рис. 10: Схема движения иона Аг+ в приповерхностном слое в режиме единично выбитых атомов
В таком режиме ион Аг+ испытывает ряд упругих соударений с атомами, входящими в кристаллическую решетку, и возможно возникновение ситуации, когда обратно отраженный ион будет иметь кинетическую энергию большую, чем энергия связи поверхностного атома (ЕШн> и). Это приводит к распылению приповерхностных атомов.
В режиме е.в.а. Аг+ при перемещении в приповерхностных слоях кристаллической решетки может быть представлен как несобственный междо-узельный атом (н.м.а). Ионный радиус Аг+ равен 157 пм [23]. Н.м.а. создает вокруг себя области гидростатического сжатия [61]. В то же время краевые дислокации (отдельные участки дислокационной
петли имеют чисто краевую компоненту) имеют области гидростатического сжатия и растяжения по разные стороны от экстраплоскости. Таким образом, упругие поля, создаваемые дислокацией, будут оказывать влияние на динамику движения ионов Аг+. Авторы полагают, что закономерность будет аналогичной приведенной в работе [59]: чем выше плотность дислокаций, тем сильнее влияние упругих полей на движение ионов и ниже коэффициент распыления.
Исходя из приведенных выше рассуждений, методика эксперимента была модифицирована до нижеследующей: производилось пошаговое ионное травление по режиму, указанному выше; на каждом шаге производилось измерение скорости ионного травления; строилась зависимость Vион.трав = /От), где Т - это глубина, на которой измерялась скорость (если считать за Т = 0 исходную поверхность пленки НЛ). Эксперимент был проведен для тонкопленочного НЛ и объемного НЛ при одних условиях ионного травления. На объемном НЛ предварительно создавалась Сг маска для травления с целью определения скорости травления при помощи СЭМ по геометрии профиля, сформированного в ходе травления.
Метод является косвенным, исходя из чего определить функцию распределения плотности дислокаций Ыдисл = /Т) по полученным экспериментальным данным не удается, однако позволяет делать выводы о пространственном распределении области с повышенной концентрацией дислокационных петель, а также провести сравнительный анализ тонкопленочного и объемного НЛ.
4.4. Протонный обмен в тонкопленочном НЛ
Процесс ПО производился в закрытом циркониевом реакторе в расплаве бензойной кислоты (СбИ5СООН) квалификации «ЧДА». Формирование протонобменного слоя на тонкопленочном НЛ производилось при Т = 175 °С. Были получены образцы с несколькими временами ПО: 30 мин, 2 ч и 4 ч.
Для того чтобы оценить влияние ПО на процесс РИТ тонкопленочного НЛ, необходимо перейти от технологического параметра - времени ПО к результирующему отношению Li/Nb в кристаллической решетке либо молярному процентному содержанию Li2O (сы2в) аналогично [62], которое устанавливается в результате замены ионов Li+ на Н+. Одним из возможных методов оценки сшо является косвенный метод, основанный на экспериментальном измерении степени анизотропии Ап = пе - п0 [63].
Функция, описывающая зависимость си2о от Ап для длины волны X = 1550 нм и комнатной температуры, имеет вид:
сЫ2о = 37.2024 - 147.9245ДП. (12)
Данная функция получена W. Schlarb и К. Betzler в работе [62] на основе обобщенного уравнения Зель-мейера, описывающего показатель преломления НЛ как функцию длины волны и стехиометрии. Так, для конгруэнтного НЛ Ап на длине волны X = 1550 нм - - 0.075 и сШо — 48.3 %.
Для экспериментального измерения величины Ап для всех образцов тонкопленочного НЛ после ПО применялся метод призменного ввода на установке МеМсоп 2010/М. Данный метод основан на возбуждении волноводных мод в тонкой пленке. Подробное описание метода приведено в статьях [64, 65]. Схема измерительного блока приведена на рис. 11.
Рис. 11. Принципиальная схема измерительного блока установки ЫеМсоп
Точность измерения ПП составляет ± 0.0001. Для определения анизотропии Ап производилось измерение ПП пе и п0. 1111 пе измерялся при распространении света с ТЕ поляризацией вдоль кристаллографической оси Y, ПП по - при распространении света с ТЕ поляризацией вдоль оси кристаллографической Ъ (для Х-среза тонкопленочного НЛ).
Поскольку толщина пленки НЛ (й = 600 нм) удовлетворяет условию существования только фундаментальной ТЕ0 моды, то полноценно восстановить профиль ПП по глубине по общепринятому методу Вентцеля-Крамерса-Бриллюэна не представляется возможным. Будем считать, что результатом измерения являются некоторые усредненные по толщине пленки ПП.
4.5. Методика эксперимента по изучению влияния ПО на кинетику РИТ тонкопленочного НЛ
На первой стадии в ходе эксперимента определялось влияние степени замещения ионов Li+ на Н+ при ПО на скорость РИТ. Травление производилось по следующему режиму: Р = 0.005 мбар, WIcp = 200 Вт, Wccp = 80 Вт, д^: QsF6 = 38: 2 см3/мин. Для всех образцов строились профили скорости РИТ по глубине (Урит = ДЬ•)) по методике, аналогичной приведенной в п. 4.3.3.
Влияние степени замещения на длительность индукционного периода определялось согласно методике и режиму травления, приведенному в разделе 4.2.
5. Результаты экспериментов и их обсуждение
5.1. Зависимость длительности
индукционного периода от скорости РИТ
На рис. 12 представлена экспериментальная зависимость длительности индукционного периода ТХР тинд от мощности на 1СР источнике плазмы для тонкопленочного НЛ.
с > о Переход не наблюдался • Окончание инд. периоде
I
< >
< >
0
100 200 300 400 500 600
¡А/к?, Вт
Рис. 12. Зависимость длительности индукционного периода ТХР от мощности на 1СР источнике
При увеличении мощности на 1СР источнике плазмы наблюдалось сокращение длительности индукционного периода.
На рис. 13 показана временная динамика степени преобразования тонкопленочного НЛ после окончания индукционного периода (режим травления: Р = 0.005 мбар, WIcp = 500 Вт, Wccp = 80 Вт, длг'- QsF6 = 38:2 см3/мин).
По окончании индукционного периода на поверхности НЛ детектируются области с пористой морфологией, соответствующие областям, где наблюдался рост LiF на поверхности (на рис. 13 выделены контурами).
Степень превращения поверхности считалась как отношение суммы площадей областей, на которых детектируется LiF, к общей площади снимка на СЭМ): а) при г = 2 мин п - 0.02, б) при г = 3 мин П — 0.20, в) при г = 4 мин п = 1.
Переход от индукционного периода к периоду роста скорости ТХР наблюдался в ходе эксперимента только для мощностей от 300 Вт включительно и выше. При более низких мощностях на 1СР источнике плазм продолжительность индукционного периода велика, поэтому тонкая пленка НЛ стравливается быстрее, чем заканчивается индукционный период и начинается детектируемый инструментально рост LiF на поверхности НЛ.
Рис. 13. Временная динамика преобразования поверхности тонкопленочного НЛв ходе ТХР: а) t = 2 мин (окончание индукционного периода), б) t = 3 мин, в) t = 4 мин
30 25
I 20 * 15 10
с о Переход не наблюдался • Окончание инд. периода
о
>о
• •
• • •
20 30 40 50 60 70 80 90 100 Угрза, нм/мин
Рис. 14. Зависимость индукционного периода ТХР от скорости травления
На рис. 14 показана зависимость длительности индукционного периода ТХР на поверхности от скоростей РИТ, которые соответствовали режимам при мощностях WICP на рис. 12. Дополнительно была произведена вариация Wicp вблизи WICP = 300 Вт с шагом в 30 Вт для уточнения зависимости.
Экспериментальная зависимость Тинд = /(Урит) качественно согласуется с выражением (8): при росте скорости РИТ (росте пересыщения паровой фазы LiF) наблюдается сокращение продолжительности индукционного периода.
При этом стоит заметить, что на кривой выколотыми точками отмечены те скорости, при которых произошло полное удаление пленки НЛ с поверхности без образования LiF. Таким образом, только закрашенные точки соответствуют истинному индукционному периоду.
Экспериментально наблюдается, что в диапазоне скоростей от 54 нм/мин до 59 нм/мин находится некоторое критическое значение скорости РИТ Укрит. При этом, если Урит < Укрит, то процесс РИТ протекает на стадии индукционного периода. LiF не детектируется на поверхности, возможно травление на полную глубину в непрерывном режиме травления.
Трактовка наличия таких скоростей может быть следующая:
1) это есть реальные критические скорости, которые соответствуют критическим пересыщениям паровой фазы;
2) при снижении WICP происходит понижение средней температуры процесса и, как следствие, увеличение Тинд до таких значений, что Тинд > tполн, где 4олн - время для полного стравливания пленки НЛ.
Подробнее обсудим первый вариант. Как было показано ранее, пересыщение газовой фазы LiF над поверхностью пластины НЛ зависит от скорости выхода продуктов реакции в ходе РИТ (выражение 5). Отсюда можно сделать вывод, что критическим скоростям Укрит соответствуют акрит, ниже которых рост тонкой пленки LiF не наблюдается. Наличие критических пересыщений является характерной особенностью зародышевого механизма роста пленки LiF (механизм Фольмера-Вебера) [19, 28].
При пересыщениях ниже критических вероятность образования зародышей с надкритическими размерами мала и процесс зародышеобразования не детектируется инструментально. Как было показано в работе [18], на стадии индукционного периода LiF детектируется только методом XPS, что говорит о крайне низких аналитических толщинах. Вероятно, что при а < акрит рост LiF реализуется по механизму дислокационного роста и не может быть обнаружен используемым инструментарием (кроме XPS).
Гипотезу о реализации механизма двумерного зародышеобразования можно косвенно подтвердить при помощи дополнительного эксперимента по измерению зависимости критических скоростей от давления Укрит = /P). Как было сказано ранее, при Укрит достигаются критические пересыщения
акрит.
Согласно уравнению Фольмера-Флада для двумерного зародышеобразования критическое пересыщение может быть выражено следующим образом [64]:
^(^крит)
N
1дА
"4Ш
(13)
где М - молярная масса вещества, осаждающегося из пара; р - плотность, е - поверхностное натяжение (аналогично плотности берется значение, равное объемным величинам), Л - предэкспоненциаль-ный множитель в выражении, описывающем число зародышей, достигающих надкритического раз-
мера в единицу времени: J = Ае кг . Коэффициент А~Р/\l2nmkT - число столкновений молекул пара с единицей площади поверхности зародыша (P - давление, m - масса молекулы пара) [66].
В итоге, из уравнения (13) (для простоты полагаем T = const; в реальности процесс РИТ можно считать изотермическим только после выхода на Тстац по истечении 3 мин с момента инициации плазмы) может быть выражена теоретическая зависимость критического пересыщения от давления:
Ig(^Kpnx) = const
N
1
IgP'
(14)
Предполагается, что зависимость Укрыт от давления будет иметь схожую форму.
Для построения экспериментальной зависимости производилась вариация давления в диапазоне от 0.005 до 0.16 мбар, фиксированными оставались следующие параметры: Wicp = 500 Вт, Wccp = 80 Вт, QAr: QSF6 = 38:2 см3/мин. Для измерения критических скоростей на каждом давлении сначала определялся режим РИТ, при котором детектировался LiF на поверхности. Для этого производилось увеличение времени травления от tm4 = 2 мин (соответствует длительности индукционного периода на данном режиме РИТ для P = 0.005 мбар) с шагом 1 мин. Затем пошагово (каждый шаг на новом образце) снижалась скорость травления (путем уменьшения времени травления с шагом 10 с) до тех пор, пока не фиксировалась Укрыт, при которой LiF на поверхности инструментально не детектировался (достижение Укрыт «сверху»).
На рис. 15 представлена экспериментальная зависимость Укрыт от давления P.
5.2. Результаты исследования
дислокационной структуры в тонкопленочном НЛ
5.2.1. Рентгенографический метод
На рис. 16 изображены нормированные ю-кри-вые тонкопленочного и объемного НЛ. Кривые нормированы на максимум интенсивности 1тах, за начало координат принято положение дифракционного пика.
1.0
<0.6
0.4
0.2
0.0
/ / 11 • L —- U 1 \ \\ Тонк. НЛ ■ Объем. НЛ
1 f *i 1 1 1
1 / i / / / Т \ I 1 и \ \
/ 1 • I 1 / 1 • 1 > * \ \ ^
\Л »
-0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 Да, °
Рис. 16: ю-кривые исходного тонкопленочного и объемного образцов LiNbOз, снятые с порядка отражения п = 1 (индекс отражения (110))
Полуширина дифракционного пика втонк составила величину 0.023°, воб = 0.012° соответственно. Сравнивая полуширины ю-кривых можно сделать качественный вывод о том, что тонкопленочный НЛ характеризуется повышенной степенью несовершенства кристаллической структуры, которая, вероятно, может быть связана с образованием дислокационных петель в объеме пленки НЛ в ходе ионной имплантации. Количественные оценки плотности дислокаций давать затруднительно ввиду сложности учета всех факторов, которые влияют на уширение пика для субмикронной пленки НЛ.
5.2.2. Метод ямок травления
На рис. 17 показано СЭМ изображения поверхности тонкопленочного НЛ после выявления ямок травления в парах ОТ + HNOз.
Рис. 15. Зависимость Укрит от давления процесса РИТ
Форма кривой на рис. 15 по форме схожа с кривой, описываемой выражением (14), что является еще одним косвенным свидетельством того, что LiF на поверхности НЛ растет по механизму двумерного зародышеобразования.
Рис. 17. Поверхность тонкопленочного НЛ после выявления ямок травления в парах
Ямки травления имеют форму искаженных ромбоэдров, что, согласно [67], характерно для X-среза. Оси ромбов параллельны друг другу, отношение длинной оси к короткой c/a — 1.45-1.55. Если обратить внимание на рис. 13, а, то видно, что для овальных областей, на которых фиксируется начало образования новой фазы LiF характерно аналогичное отношение c/a. Это еще раз подтверждает факт того, что процесс зародышеобразования новой фазы LiF на поверхности НЛ начинается преимущественно с термодинамически благоприятных мест - мест выхода дислокаций на поверхность.
Усредненная по пяти снимкам плотность дислокаций на исходной поверхности тонкопленочного НЛ составила величину Ытонк = 5.7-107 см-2.
Плотность дислокаций для объемного НЛ, определяемая методом ямок травления, составляет величину Ыоб = 5.0 104 см-2 [56]. Таким образом, метод ямок травления демонстрирует увеличение плотности дислокаций в тонкопленочном НЛ по сравнению с объемным на 3 порядка.
5.2.3. Профиль ионного травления тонкопленочного НЛ по глубине
Как было показано выше, плотность дислокаций на исходной поверхности тонкопленочного НЛ Nrnom = 5.7 • 107 см-2. После удаления 100 нм НЛ методом ионного травления по режиму P = 0.005 мбар, Wicp = 500 Вт, Wccp = 80 Вт, Qa- = 20 см3/мин плотность дислокаций, выявленная методом ямок травления в парах, составила величину Ытонк = 9.7 • 106 см-2. Однако, при дальнейшем последовательном удалении 100 нм НЛ выявить ямки травления уже не удавалось. При утонении пленки, по-видимому, снижалась плотность дислокаций и росло время выявления ямок в парах. При длительном выдерживании в парах (более 3 ч) поверхность значительно деградировала из-за эффекта похожего на блистеринг (рис. 18), который мог, вероятно, возникать из-за накопления водорода в приповерхностном слое НЛ.
где Т - глубина от исходной поверхности тонкопленочного НЛ (рис. 19). Аналогичный профиль был построен для объемного НЛ.
Рис. 18. Поверхностные дефекты, образующиеся в ходе выявления ямок травления в парах, на поверхности тонкопленочного НЛ
Для изучения профиля дефектов были построены зависимости скорости ионного травления по режиму, описанному в разделе 4.3.3, Уион.трав = /(Т),
Рис. 19. Профиль скорости ионного травления НЛ по глубине
По графику видно, что профиль скорости ионного травления тонкопленочного НЛ имеет сигмо-идальную форму. Скорость ионного травления растет при удалении приповерхностных слоев и на глубине Т — 300 нм начинает приближаться к скорости ионного травления объемного НЛ Х-срез (CQT). Ситуация наблюдалась при двух различных мощностях WCCP (энергии ионов, соответствующие Wccp = 100 Вт и Wccp = 120 Вт, равны Еаг+ — 200 эВ и Еаг+ — 250 эВ соответственно).
Профиль скорости травления может быть аппроксимирован с некоторой погрешностью следующей функцией:
пю =
b + с exp(-kh2)'
(15)
Такой вид функции был выбран из следующих соображений: согласно данным в работе [69], описывающим распределения различных ионов с энергиями 200-400 кэВ после ионной имплантации в LiNbOз, распределение ионов Не+ после ионной имплантации единственное в ряде ионов Н, Не, Ве, В, F, Р и Ga является гауссовым, остальные элементы имеют распределение Пирсона IV типа. Распределение ионов имплантируемого материала коррелирует с распределением плотности индуцированных в кристаллической решетке точечных дефектов, и, как следствие, неравновесных дефектов, возникающих в ходе отжига, следующего за ионной имплантацией).
Коэффициенты в распределении, предположительно, имеют следующий физический смысл: А -коэффициент пропорциональности между скоростью ионного травления и плотностью дислокаций; Ь - это плотность дислокаций объёмного НЛ, использованного для изготовления пластин тонкопленочного НЛ; сехр(— кк2) - функция распределения дислокаций, дающих добавку к плотности дислокаций объемного НЛ после процесса ионной имплантации, характерную для тонкопленочного
НЛ (см. рис. 5, п. 3.1); коэффициент k определяется пространственным распределением структурных дефектов, которые возникают в ходе ионной имплантации.
Выше приведенные результаты являются косвенным свидетельством того, что плотность дислокаций максимальна на исходной поверхности тонкопленочного НЛ (по методу ямок травления Ытонк = 5.7Т07 см-2) и спадает до величин, близких к плотности дислокаций объемного НЛ (по методу ямок травления Ыоб = 5.0 104 см-2), использованного для изготовления пластин тонкопленочного НЛ, на глубине — 300 нм и далее остается постоянной.
5.3. Влияние плотности дислокаций на длительность индукционного периода ТХР
Количественная оценка плотности дислокаций на исходной поверхности объемного и тонкопленочного НЛ, полученная методом ямок травления, показывает, что плотности дислокаций различаются почти на 3 порядка.
Как было показано в разделе 2.2, увеличение плотности дислокаций должно приводить к сокращению индукционного периода ТХР.
Эксперимент, выполненный по методике, описанной в п. 4.2, показал, что длительность индукционного периода на режиме P = 0.005 мбар, Wicp = 500 Вт, Wccp = 80 Вт, Qai-: Qsf6 = 38: 2 см3/мин для тонкопленочного НЛ составила величину ттонк = 2 мин, а для объемного НЛ Тоб = 13 мин, соответственно.
Если бы зародыши образовывались только на дислокациях (в приближении X = Ыдисл в выражении 7), то при Ытонк = 5.7 • 107 см-2 и Кб = 5.0 • 104 см-2 индукционный период Ттонк должен быть в ~ 33 раза короче, чем Тоб. В реальности Тоб и Ттонк различаются только в 6.5 раз на заданном режиме. Вероятной причиной расхождения является тот факт, что Ыдисл вносят вклад в общую поверхностную плотность центров зародышеобразования X (который не превышает 1% для тонкопленочного НЛ, если принять за X ~ 109 см-2 [27], и пренебрежимо мал для объемного НЛ), но подавляющая часть зародышей образуется флуктуационным путем.
Еще одной вероятной причиной расхождения может быть точность оценки плотности дислокаций методом ямок травления, поскольку по теории Кабреры эффективно выявляются только винтовые дислокации [69]).
5.4. Влияние протонного обмена на кинетику ТХР в ходе РИТ тонкопленочного НЛ
В ходе работы был проведен ПО в расплаве бензойной кислоты на трех образцах тонкопленочного НЛ по режимам: Тпо = 175 °С, ti = 0.5 ч, t2 = 2 ч, t3 = 4 ч.
На каждом из образцов после ПО были измерены усредненные по толщине пленки НЛ коэффициенты анизотропии Ап = пе -п0 по методике, приведенной в п. 4.4, по которым вычислены молярные процентные содержания Li2O ^20 по выражению (12).
На рис. 20 показана теоретическая зависимость
Ап =
Рис. 20. Зависимость коэффициента анизотропии от cu2o на длине волны X = 1550 нм
После процесса ПО получены значения An для образцов, согласно которым по графику на рис. 20 определены молярные процентные содержания
Li2O.
Необходимо иметь в виду, что на графике отмечены некоторые усредненные значения молярных процентных содержаний Li2O для разных образцов, так как процесс ПО проходил на всю глубину пленки НЛ.
Для всех образцов были получены профили скорости РИТ по режиму: P = 0.005 мбар, Wicp = 200 Вт, Wccp = 80 Вт, Qat : Qsf6 = 38:2 см3/мин (рис. 21).
• cLho = 36.2 % ■ Сц2о = 28.3 % т си,о = 22.3 %
X ИСХОДНЫЙ
120
100
600
Рис. 21. Профиль скорости РИТ по глубине для разных степеней замещения
Из рис. 21 видно, что скорость травления имеет некоторый профиль по глубине, который определяется вероятно диффузионным профилем концентрации ионов Н+, внедренных в ходе ПО. Скорость
РИТ образца после ПО й = 0.5 ч у поверхности существенно превышает скорость травления исходного тонкопленочного НЛ и составляет —90 нм/мин. На глубине около —280 нм скорость РИТ снижается до величины исходного тонкопленочного НЛ. При этом, согласно оценкам авторов, глубина 280-300 нм соответствует глубине диффузии ионов Н+ в ходе процесса ПО в течение 0.5 ч при Т=175 °С. Схожий характер изменения скорости травления по глубине наблюдается и для образца после ПО в течение 2 ч (с учетом предполагаемой авторами глубина диффузии ионов Н+ - около 600 нм).
Стоит обратить внимание на экспериментальные точки для образца с временем ПО Ь = 4 ч. Модификация пленки, как полагают авторы, произошла на глубину, значительно большую, чем физическая толщина пленки НЛ. Это привело к тому, что скорость РИТ практически не изменяется по всей глубине пленки НЛ. Такой профиль позволяет подбирать время для обеспечения любой глубины травления (поскольку скорость равномерна по всей глубине), что крайне важно для формирования заданного профиля канального волновода ФИС.
Однако подчеркнем одну характерную особенность длительного процесса ПО в тонкопленочном НЛ. Она связана с растрескиванием пленки, типичный вид которой представлен на рис. 22.
Рис. 22. СЭМ изображение поверхности тонкопленочного НЛ после ПО в течение 4 ч
Такое растрескивание характерно для механизма релаксации растягивающих напряжений, наличие которых в пленках НЛ на пластинах-носителях из Si подтверждается экспериментами, приведенными в ряде исследований, например, [70].
Для дальнейших исследований влияния модификации тонкопленочного НЛ с помощью ПО на длительность индукционного периода было выбрано время ПО t2 = 2 ч. В качестве экспериментального режима был выбран режим травления P = 0.005 мбар, Wicp = 500 Вт, Wccp = 80 Вт, Qa-: Qsf6 = 38:2 см3/мин. Напомним, что для исходного тонкопленочного НЛ длительность индукционного периода на данном режиме составляет величину т = 2 мин.
Для исходного тонкопленочного НЛ скорость РИТ на данном режиме составляла 86 нм/мин, на образце после ПО t2 = 2 ч она достигла значений 250 нм/мин.
После процесса травления до остаточных толщин 100 нм LiF на поверхности НЛ не детектировался методом EDS. Дополнительно было проведено исследование химического состава поверхности при помощи метода рентгеновской фотоэлектронной спектроскопии (XPS) по методике, описанной в предыдущей работе [18] (рис. 23). Данная методика позволяет провести «поверхностный» анализ химического состава и определять наличие на поверхности даже ультратонких пленок.
0 200 400 600 800 1000
Есвязи' эВ
Рис. 23. XPS спектр тонкопленочного НЛ после ПО и РИТ по режиму, описанному выше
На поверхности образца методом XPS детектировались следующие элементы и соответствующие им энергии связи: Elîiba = 54.47 эВ, Elîibb = 56.52 эВ, Efis = 685.71 эВ, ENb3d = 208.17 эВ, EoisA = 531.08 эВ, Eoisb = 532.62 эВ. Пики Fe, Ni, C принадлежат загрязнениям на поверхности образца. На поверхности процентное соотношение элементов следующее: Liis: Nb3d: Ois: Fis: Cis = 30.8: 20.9: 35.2: 5.3: 7.7 ат.%.
По экспериментальному спектру необходимо отметить, что детектируются пики Li с двумя энергиями, среди которых ELi1sA = 54.47 эВ соответствует LiNbO3, а пик ELi1sB = 56.52 эВ - LiF. По зависимости средней глубины выхода фотоэлектронов от энергии электронов (универсальная кривая, рис. 25.27 в [23]) можно сделать вывод, что для энергий — 50 эВ глубина выхода не превышает 5 Â. Это означает, что возможный верхний предел толщины LiF на поверхности, при котором еще возможно детектирование пика, Elusa от НЛ не превышает 5 Â. При этом процентное содержание F равно 5.3 ат. % при практически трехкратном росте скорости РИТ относительно исходного тонкопленочного НЛ. Это может быть свидетельством того, что после процесса ПО пар LiF демонстрирует снижение пересыщения.
Снижение пересыщения паровой фазы LiF после процесса ПО приводит к значительному росту продолжительности индукционного периода, по окончании которого может наблюдаться LiF на поверхности, так что Тинд » W, где t„олн - это время полного травления пленки НЛ при наблюдаемых
скоростях. Этот вывод был проверен экспериментально для всего диапазона значений мощностей WICP до 600 Вт, Wccp до 300 Вт, возможных для использования на установке Этна-100. По крайней мере, для выделенного диапазона мощностей данный вывод является справедливым.
Таким образом, протонный обмен является эффективным способом увеличения скорости РИТ тонкопленочного НЛ во фторсодержащей плазме. Вместе с тем предварительный ПО позволяет значительно увеличивать продолжительность индукционного периода образования LiF, что дает возможность травить НЛ на полную глубину в непрерывном режиме.
Травление исходного тонкопленочного НЛ, напротив, требует соблюдения условия tpит <Тинд, для того, чтобы избежать повышенной шероховатости за счет образования LiF на поверхности.
6. Заключение
В ходе работы было представлено исследование процесса вторичной конденсации LiF в ходе реактивного ионного травления тонкопленочного нио-бата лития в плазме газовой смеси Аг + SF6.
Было показано, что длительность индукционного периода образования LiF на поверхности пластины НЛ сокращается с ростом скорости РИТ. Впервые определено, что процесс РИТ НЛ во фтор-содержащей плазме характеризуется критическими скоростями РИТ, выше которых наблюдается рост LiF на поверхности по механизму двумерного заро-дышеобразования Фольмера-Вебера.
Также было показано, что критические скорости снижаются с ростом давления процесса.
Впервые было экспериментально определена плотность дислокаций в тонкопленочном НЛ, изготовленном по методу управляемого скола. Также было показано экспериментально влияние повышенной плотности дислокаций на длительность индукционного периода образования LiF по сравнению с объемным НЛ.
Дополнительно было продемонстрировано влияние предварительного ПО на кинетику РИТ.
7. Благодарности
Авторы выражают благодарность ПАО «Пермская научно-производственная приборостроительная компания» (ПАО «ПНППК) и лично А.Г. Андрееву и В.К. Струку за оказанную помощь при проведении исследований.
Авторы благодарны УЦКП «Современные нанотехнологии» и лично В.Я. Шуру и Е. В. Шишкиной за помощь в проведении исследований методом XPS.
Также авторы выражают признательность Т.Д. Ратманову, А.Ю Клочкову за ценные советы по проведению экспериментов и интерпретации полученных данных.
Список литературы
1. Агрузов П. и др. Ниобат лития как платформа для создания интегрально-оптических устройств радиофотоники // Фотон-Экспресс. 2023. № 6. C. 124-125.
2. Zhu D. et al. Integrated photonics on thin-film lithium niobate // Advances in Optics and Photonics. 2021. Vol. 13. N. 2. P. 242-352.
3. Poberaj G. et al. Lithium niobate on insulator (LNOI) for micro-photonic devices // Laser Photonics Rev. 2012. Vol. 6. N. 4. P. 488-503.
4. Jia Y., Wang L., Chen F. Ion-cut lithium niobate on insulator technology: Recent advances and perspectives // Appl. Phys. Rev. 2021. Vol. 8, 011307.
5. Kaufmann F. et al. Redeposition-free inductively coupled plasma etching of lithium niobate for integrated photonics // Nanophotonics. 2023. Vol. 12. N. 8. P. 1601-1611.
6. Ulliac G. et al. Argon plasma inductively coupled plasma reactive ion etching study for smooth side-wall thin film lithium niobate waveguide application // Optical Materials. 2016. Vol. 53. P. 1-5
7. Wang C. et al. Second harmonic generation in nano-structured thin-film lithium niobate waveguides // Optics Express. 2017. Vol. 25. N. 6. P. 6963-6973.
8. Kozlov A. et al. Reactive ion etching of X-cut LiNbO3 in an ICP/TCP system for the fabrication of an optical ridge waveguide // Appl. Sci. 2023. Vol. 13. N. 4, 2097.
9. Hu H. et al. Plasma etching of proton-exchanged lithium niobate // Journal of Vacuum Science and Technology A. 2006. V. 24. N. 4. P. 1012-1015.
10. Baghban M. A. et al. Low-loss ridge waveguides in thin film lithium niobate on-insulator (LNOI) fabricated by reactive ion etching // OSA Technical Digest. 2016, ITu1A.2.
11. Bahadori M. et al. High performance fully etched isotropic microring resonators in thin-film lithium niobate on insulator platform // Optics Express. 2019. Vol. 27. N. 15. P. 22025-22039.
12. Козлов А. и др. Физико-химические особенности реактивного ионного травления тонкопленочного ниобата лития для формирования канальных волноводов фотонных интегральных схем // Мат. 26-й Всерос. конф. «Актуальные проблемы физической и функциональной электроники». Ульяновск: изд-во Ульяновского филиала Института радиотехники и электроники им. В. А. Котельникова РАН, 2023. С. 222-223
13. Гуляев В. и др. Высокоскоростное плазмохими-ческое травление ниобата лития // Конденсированные среды и межфазные границы. 2010. Т. 12, № 4. C. 360-368.
14. Гуляев В. В. Исследование кинетики и механизмов взаимодействия газоразрядной фторсодер-жащей плазмы с поверхностью LiNbO3: дис. канд. техн. наук. Воронеж: изд-во Воронежского гос. ун-та, 2011. 136 с.
15. Волков А. В. Разновидности систем плазменного травления и их особенности: учеб.-метод. пособие. Самара: изд-во Самарского гос. аэрокосм. ун-та, 2013. 17 с.
16. Воробьев А. Х., Богданов А. В. Кинетика нелинейных химических систем. М.: изд-во Московского ун-та, 2019. 158 с.
17. Neidhart T. et. al. Potential sputtering of lithium fluoride by slow multicharged ions // Phys. Rev. Lett. 1995. Vol. 74. N. 26. P. 5280-5283.
18. Козлов А.А. и др. Исследование кинетики реактивного ионного травления тонкопленочного ниобата лития во фторсодержащей плазме. // Вестник Пермского университета. Физика. 2024. № 1. С. 56-71.
19.Хасс Г., Тунн Р., Франкомб Г. Физика тонких плёнок. Т. 1. М.: Мир, 1968. 343 с.
20. Данилин Б. С. Вакуумное нанесение тонких пленок. М.: Энергия. 1967. 304 с.
21. Nagata H. et al. Growth of crystalline LiF on CF4 plasma etched LiNbO3 substrates // Journal of Crystal Growth. 1998. Vol. 187. N. 3-4. P. 573-576.
22. Ситникова В. Е., Пономарева А. А., Успенская М. В. Методы термического анализа: практикум. СПб: изд-во ун-та ИТМО, 2021. 152 с.
23. Григорьев И. С., Мейлихов Е. З. (ред.) Физические величины: справочник. М.: Энергоатомиз-дат, 1991. 1232 с.
24. Рябцев Н. Г. Материалы квантовой электроники. М.: Советское радио, 1972. 384 с.
25. Burton W. K., Cabrera N., Frank F. C. The growth of crystals and the equilibrium structure of their surfaces // Philosophical Transactions of the Royal Society of London A. 1951. Vol. 243. N. 866. P. 299358.
26. Чернов А. А. и др. Современная кристаллография. Т. 3. Образование кристаллов. М.: Наука, 1980. 401 с.
27. Дубровский В. Г., Цырлин Г. Э. Кинетика роста тонких пленок при зародышевом механизме формирования слоев // Физика и техника полупроводников. 2005. Т. 39. №. 11. С. 1312-1319.
28. Мочалов И.В. Выращивание оптических кристаллов. СПб: изд-во ун-та ИТМО, 2012. 75 с.
29. Широков Ю.Г. Механохимия. Теоретические основы: учебное пособие. Иваново: изд-во Ивановского ун-та, 2015. 214 с.
30. Krasulin Y. L. Dislocations as active centers in top-ochemical reactions // Theoretical and Experimental Chemistry. 1969. Vol. 3. N. 1. P. 31-35.
31. Леонюк Н. И. и др. Кристаллография: зарождение, рост и морфология кристаллов. М: изд-во «Юрайт», 2017. 152 с.
32. Болдырев В. В. Топохимия термического разложения твердых веществ // Успехи химии. 1973. V. 42. N. 7. P. 1161-1183.
33. Boldyrev V. V. The effects of various kinds of preliminary treatments on topochemical reactions //
Proc. of Int. Conf. "Reaction Kinetics in Heterogeneous Chemical Systems". Amsterdam: Elsevier, 1975. P. 159.
34. Болдырев В. В. Управление химическими реакциями в твердой фазе // Соросовский образовательный журнал. 1996. №. 5. С. 49-55.
35. Brown M. E. The Prout-Tompkins rate equation in solid-state kinetics // Thermochimica Acta. 1997. Vol. 300. N. 1-2. P. 93-106.
36. Van der Put P. J. The Inorganic Chemistry of Materials: How to Make Things out of Elements. New York: Springer, 2013. 406 p.
37. Ulliac G. et al. Ultra-smooth LiNbO3 micro and nanostructures for photonic applications // Microelectronic Engineering. 2011. Vol. 88. N. 8. P. 24172419.
38. LiX. P., Chen K. X., Hu Z.F. Low-loss bent channel waveguides in lithium niobate thin film by proton exchange and dry etching // Optical Materials Express. 2018. Vol. 8. N. 5. P. 1322-1327.
39. Aryal A. et al. High-quality dry etching of LiNbO3 assisted by proton substitution through ^-plasma-surface treatment // Nanomaterials. 2022. Vol. 12. N. 16, 2836.
40. Nekvindovâ P. et al. Annealed proton exchanged optical waveguides in lithium niobate: differences between the X-and Z-cuts // Optical Materials. 2002. Vol. 19. N. 2. P. 245-253.
41. Гусева А. Ф., Нейман А. Я., Анимица И. Е. Твердофазные реакции при получении и эксплуатации неорганических материалов. Екатеринбург: изд-во УрФУ, 2005. 76 с.
42. Boes A. et al. Status and potential of lithium niobate on insulator (LNOI) for photonic integrated circuits // Laser Photonics Rev. 2018. Vol. 12. N. 4, 1700256
43. Li Q. et al. Characterizations of single-crystal lithium niobate thin films // Crystals. 2022. Vol. 12. N. 5, 667.
44.Xue G. et al. Integrated fabrication and ferroelectric domain adjustment of lithium niobate single crystal films based on silicon substrate // Materials and Design. 2022. Vol. 215, 110447.
45. Wesch W., Wendler E. Ion beam modification of solids. Cham: Springer, 2016. 534 p.
46. Suni T. Direct wafer bonding for MEMS and microelectronics: Dr. Sci. dissertation. Tietoe: VTT Research Center, 2006. 140 p.
47. Ponomarev R. et al. Investigation of near-surface layer dislocation density of lithium niobate single crystal wafers using chemical etching // Journal of Physics: Conference Series. 2021. Vol. 2086. N. 1, 012031.
48. Brown S. D. et al. Dislocation distributions in ion implanted optical waveguides // Nuclear Instruments and Methods in Physics Research Section B. 1994. Vol. 94. N. 4. P. 417-423.
49. Gischkat T. et al. Annealing behavior of lithium niobate irradiated with He-ions at 100 K // Nuclear
Instruments and Methods in Physics Research Section B. 2009. V. 267. N. 8-9. P. 1492-1495.
50. Martynenko Y. V., Moscovkin P.G. Dislocation structure formation due to ion implantation // Radiation Effects and Defects in Solids. 1994. Vol. 129. N. 3-4. P. 193-198.
51. Butcher P. N., March N. H., Tosi M. P. (eds.). Crystalline Semiconducting Materials and Devices. New York: Springer, 2013. 664 p.
52. Lang E. et al. Defect structures as a function of ion irradiation and annealing in LiNbO3 // Thin Solid Films. 2023. Vol. 768, 139719.
53. Dolabella S. et al. Correlated disorder by defects clusters in LiNbO3 single crystals after crystal ion-slicing // Materials and Design. 2023. Vol. 231, 112001.
54. Han H., Cai L., Hu H. Optical and structural properties of single-crystal lithium niobate thin film // Optical Materials. 2015. Vol. 42. P. 47-51.
55. Hu X. et al. Size-controlled polarization retention and wall current in lithium niobate single-crystal memories // ACS Applied Materials and Interfaces. 2021. Vol. 13. N. 14. P. 16641-16649.
56. Сосунов А. В., Петухов И. В. и др. Влияние предварительного отжига пластин ниобата лития на характеристики протонообменных волноводов // Кристаллография. 2022. Т. 67. № 6. С. 982989.
57. Сорокина И. И., Педько Б. Б. Способ травления монокристаллов ниобата лития. Патент СССР: №1 583 479 A1, 1990.
58. Педько Б. Б., Сорокина И. И., Бойцова К. Н. Способ травления оксидных нелинейных диэлектрических и оптических монокристаллов. Патент РФ: №2 558 898 C1, 2014.
59. Гордеева Г. В. и др. Распыление и радиационный блистеринг стали 0Х15Н16М3Б после холодной деформации и рекристаллизации // Атомная энергия. 1990. Т. 68. №4. С. 244-247.
60. Sigmund P. Theory of sputtering. I. Sputtering yield of amorphous and polycrystalline targets // Physical Review. 1969. Vol. 184. N. 2, 383.
61. Новиков И. И. Дефекты кристаллического строения металлов. М.: Металлургия, 1975. 208 с.
62. Schlarb U., Betzler K. Refractive indices of lithium niobate as a function of wavelength and composition // Journal of Applied Physics. 1993. Vol. 73. N. 7. P. 3472-3476
63. Атучин В. В. Зависимость показателей преломления ниобата лития от химического состава кристалла // Вестник СибГУТИ. 2009. №. 3. С. 39-45.
64. Tien P. K. Light waves in thin films and integrated optics // Applied Optics. 1971. Vol. 10. N 11. P. 2395-2413.
65. Ulrich R., Torge R. Measurement of thin film parameters with a prism coupler // Applied Optics. 1973. V. 12. №. 12. P. 2901-2908.
66. Сутугин А. Г. Спонтанная конденсация пара и образование конденсационных аэрозолей // Успехи химии. 1969. Т. 38. №. 1. С. 166-191.
67. Nassau K., Levinstein H. J., Loiacono G. M. Ferroelectric lithium niobate. 1. Growth, domain structure, dislocations and etching // Journal of Physics and Chemistry of Solids. 1966. Vol. 27. N. 6-7. P. 983-988.
68. Wilson R. G. et al. Ion implantation in LiNbO3 // Proc. Symp. "Defect Properties and Processing of High-technology Nonmetallic Materials". Boston, Massachusetts, USA. 1984. Vol. 24. P. 181.
69. Ван Бюрен X Г. Дефекты в кристаллах. М.: изд-во иностранной литературы, 1962. 610 с.
70. Takigawa R. et al. Residual stress in lithium niobate film layer of LNOI/Si hybrid wafer fabricated using low-temperature bonding method // Micromachines. 2019. Vol. 10. N. 2. P. 136.
References
1. Agruzov P. et al. Niobat litiia kak platforma dlia sozdaniia integral'no-opticheskikh ustroistv radiofotoniki [Lithium niobate as a platform for the creation of integrated optical devices of radiophoton-ics]. Foton-Ekspress [Photon-Express], 2023, no. 6, pp. 124-125. (In Russian).
2. Zhu D. et al. Integrated photonics on thin-film lithium niobate. Advances in Optics and Photonics, 2021, vol. 13, no. 2, pp. 242-352.
3. Poberaj G. et al. Lithium niobate on insulator (LNOI) for micro-photonic devices. Laser Photonics Rev, 2012, vol. 6, no. 4, pp. 488-503.
4. Jia Y., Wang L., Chen F. Ion-cut lithium niobate on insulator technology: Recent advances and perspectives. Appl. Phys. Rev., 2021, vol. 8, 011307.
5. Kaufmann F. et al. Redeposition-free inductively coupled plasma etching of lithium niobate for integrated photonics. Nanophotonics, 2023, vol. 12, no. 8, pp. 1601-1611.
6. Ulliac G. et al. Argon plasma inductively coupled plasma reactive ion etching study for smooth side-wall thin film lithium niobate waveguide application. Optical Materials, 2016, vol. 53, pp. 1-5.
7. Wang C. et al. Second harmonic generation in nano-structured thin-film lithium niobate waveguides. Optics Express, 2017, vol. 25, no. 6, pp. 69636973.
8. Kozlov A. et al. Reactive ion etching of X-cut LiNbO3 in an ICP/TCP system for the fabrication of an optical ridge waveguide. Appl. Sci, 2023, vol. 13, no. 4, 2097.
9. Hu H. et al. Plasma etching of proton-exchanged lithium niobate. Journal of Vacuum Science and Technology A, 2006, vol. 24, no. 4, pp. 1012-1015.
10. Baghban M. A. et al. Low-loss ridge waveguides in thin film lithium niobate on-insulator (LNOI) fabricated by reactive ion etching. OSA Technical Digest, 2016, ITu1A.2.
11. Bahadori M. et al. High performance fully etched isotropic microring resonators in thin-film lithium niobate on insulator platform. Optics Express, 2019, vol. 27, no. 15, pp. 22025-22039.
12. Kozlov A. et al. Fiziko-khimicheskie osobennosti reaktivnogo ionnogo travleniia tonkoplenochnogo niobata litiia dlia formirovaniia kanal'nykh volnovodov fotonnykh integral'nykh skhem [Physico-chemical features of reactive ion etching of thin-film lithium niobate for the fabrication of channel waveguides of photonic integrated circuits]. Ak-tual'nye problemy fizicheskoi i funktsional'noi el-ektroniki [Current Problems of Physical and Functional Electronics]: Conference proceedings. Ulyanovsk, 2023, pp. 222-223. (In Russian).
13. Gulyaev V. V. et al. High-rate plasma-chemical etching of lithium niobate. Condensed Matter and Interphases, vol. 12, no. 4, pp. 360-368 (In Russian).
14. Gulyaev V. V. Issledovanie kinetiki i mekhanizmov vzaimodeistviia gazorazriadnoi ftorsoderzhashchei plazmy s poverkhnost'iu LiNbO3 [Study of the kinetics and mechanisms of interaction of gas-discharge fluorine-containing plasma with the LiNbO3 surface] Cand. tech. sci. thesis. Voronezh, 2011. 136 p. (In Russian).
15. Volkov. A Raznovidnosti sistem plazmennogo trav-leniya i ix osobennosti [Types of plasma etching systems and their features]. Samara: Samara State University, 2013. 17 p. (In Russian)
16. Vorobyev A. Kh., Bogdanov A. V. Kinetika nelineinykh khimicheskikh sistem [Kinetics of Nonlinear Chemical Systems]. Moscow: Moscow State University, 2019. 158 p. (In Russian).
17. Neidhart T. et. al. Potential sputtering of lithium fluoride by slow multicharged ions. Phys. Rev. Lett., 1995, vol. 74, no. 26, pp. 5280-5283.
18. Kozlov A. A., Salgaeva U. O., Zhuravlev V. A., Volyntsev A. B. The study of the kinetics of thin-film lithium niobate reactive ion etching in a fluorine-containing plasma. Bulletin of Perm University. Physics, 2024, no. 1, pp. 56-71.
19. Hass G., Tunn R., Frankomb G. Physics of Thin Films, vol. 1. Oxford: Academic Press, 1963. 350 p.
20. Danilin B. Vakuumnoe nanesenie tonkikh plenok [Vacuum Deposition of Thin Films]. Moscow, 1967. 304 p. (In Russian).
21. Nagata H. et al. Growth of crystalline LiF on CF4 plasma etched LiNbO3 substrates. Journal of Crystal Growth, 1998, vol. 187, no. 3-4, pp. 573-576.
22. Sitnikova V. E., Ponomareva A. A., Uspenskaya M.V. Metody termicheskogo analiza [Thermal Analysis Methods]. St. Petersburg: ITMO University, 2021. 152 p. (In Russian).
23. Grigoriev I. S., Meilikhov E. Z., Radzig A. A. Handbook of Physical Quantities. Boca-Raton: CRC Press, 1997. 1548 p.
24. Ryabcev N. G. Materialy kvantovoi elektroniki [Quantum Electronic Materials]. Moscow: Soviet Radio, 1972. 384 p. (In Russian).
25. Burton W. K., Cabrera N., Frank F. C. The growth of crystals and the equilibrium structure of their surfaces. Philosophical Transactions of the Royal Society of London A, 1951, vol. 243, no. 866, pp. 299358.
26. Chernov A. A. et al. Sovremennaia kristallografiia [Modern Crystallography] : in 4 vols. Vol. 3. Crystal formation. Moscow: Nauka, 1980. 401 p. (In Russian)
27. Dubrovsky V. G., Cirlin G. E. Growth kinetics of thin films formed by nucleation during layer formation. Semiconductors, 2005, vol. 39, no. 11, pp. 1267-1274.
28. Mochalov I.V. Vyrashchivanie opticheskikh kristal-lov [Growing of Optical Crystals], St. Petersburg: ITMO University, 2012. 75 p. (In Russian).
29. Shirokov Yu.G. Mekhanokhimiia. Teoreticheskie osnovy [Mechanochemistry. Theoretical Foundations]. Ivanovo: Ivanovo University, 2015. 214 p. (In Russian).
30. Krasulin Y. L. Dislocations as active centers in top-ochemical reactions. Theoretical and Experimental Chemistry, 1969, vol. 3, no. 1, pp. 31-35.
31. Leonyuk N. I. et al. Kristallografiia: zarozhdenie, rost i morfologiia kristallov [Crystallography: Nu-cleation, Growth and Morphology of Crystals] Moscow: "Urait" Publishing House, 2017. 152 p. (In Russian).
32. Boldyrev V. V. Topochemistry of the thermal decomposition of solid substances. Russian Chemical Reviews, 1973, vol. 42, no. 7, pp. 515-528.
33. Boldyrev V. V. The effects of various kinds of preliminary treatments on topochemical reactions. Proc. of Int. Conf "Reaction Kinetics in Heterogeneous Chemical Systems". Amsterdam: Elsevier, 1975. P. 159.
34. Boldyrev V. V. Upravlenie khimicheskimi reaktsi-iami v tverdoi faze [Control of chemical reactions in the solid phase]. Sorosovskii obrazovatel'nyi zhurnal [Soros Educational Journal], 1996, no. 5, pp. 49-55 (In Russian).
35. Brown M. E. The Prout-Tompkins rate equation in solid-state kinetics. Thermochimica Acta, 1997, vol. 300, no. 1-2, pp. 93-106.
36. Van der Put P. J. The Inorganic Chemistry of Materials: How to Make Things Out of Elements. New York: Springer, 2013. 406 p.
37. Ulliac G. et al. Ultra-smooth LiNbO3 micro and nanostructures for photonic applications. Microelectronic Engineering, 2011, vol. 88, no. 8, pp. 2417-2419.
38. Li X. P., Chen K. X., Hu Z.F. Low-loss bent channel waveguides in lithium niobate thin film by proton exchange and dry etching. Optical Materials Express, 2018, vol. 8, no. 5, pp. 1322-1327.
39. Aryal A. et al. High-quality dry etching of LiNbO3 assisted by proton substitution through ^-plasma-surface treatment. Nanomaterials, 2022, vol. 12, no. 16, 2836.
40. Nekvindova P. et al. Annealed proton exchanged optical waveguides in lithium niobate: differences between the X- and Z-cuts. Optical Materials, 2002, vol. 19, no. 2, pp. 245-253.
41. Guseva A. F., Neiman A. Ya., Animitsa I. E. Tver-dofaznye reaktsii pri poluchenii i ekspluatatsii ne-organicheskikh materialov [Solid-phase reactions in the production and use of inorganic materials] : a study guide. Yekaterinburg: Ural Federal University, 2005. 76 p. (In Russian).
42. Boes A. et al. Status and potential of lithium niobate on insulator (LNOI) for photonic integrated circuits. Laser Photonics Rev., 2018, vol. 12, no. 4, 1700256.
43. Li Q. et al. Characterizations of single-crystal lithium niobate thin films. Crystals, 2022, vol. 12, no. 5, 667.
44. Xue G. et al. Integrated fabrication and ferroelectric domain adjustment of lithium niobate single crystal films based on silicon substrate. Materials and Design, 2022, vol. 215, 110447.
45. Wesch W., Wendler E. Ion Beam Modification of Solids. Cham: Springer, 2016. 534 p.
46. Suni T. Direct wafer bonding for MEMS and microelectronics: Dr. Sci. dissertation. Tietoe: VTT Research Center, 2006. 140 p.
47. Ponomarev R. et al. Investigation of near-surface layer dislocation density of lithium niobate single crystal wafers using chemical etching. Journal of Physics: Conference Series, 2021, vol. 2086, no. 1, 012031.
48. Brown S. D. et al. Dislocation distributions in ion implanted optical waveguides. Nuclear Instruments and Methods in Physics Research Section B, 1994, vol. 94, no. 4, pp. 417-423.
49. Gischkat T. et al. Annealing behavior of lithium ni-obate irradiated with He-ions at 100 K. Nuclear Instruments and Methods in Physics Research Section B, 2009, vol. 267, no. 8-9, pp. 1492-1495.
50. Martynenko Y. V., Moscovkin P.G. Dislocation structure formation due to ion implantation. Radiation Effects and Defects in Solids, 1994, vol. 129, no. 3-4, pp. 193-198.
51. Butcher P. N., March N. H., Tosi M. P. (eds.). Crystalline Semiconducting Materials and Devices. New York: Springer, 2013. 664 p.
52. Lang E. et al. Defect structures as a function of ion irradiation and annealing in LiNbO3. Thin Solid Films, 2023, vol. 768, 139719
53. Dolabella S. et al. Correlated disorder by defects clusters in LiNbO3 single crystals after crystal ion-slicing. Materials and Design, 2023, vol. 231, 112001.
54. Han H., Cai L., Hu H. Optical and structural properties of single-crystal lithium niobate thin film. Optical Materials, 2015, vol. 42, pp. 47-51.
55. Hu X. et al. Size-controlled polarization retention and wall current in lithium niobate single-crystal memories. ACS Applied Materials and Interfaces, 2021, vol. 13, no. 14, pp. 16641-16649.
56. Sosunov A. V., Petukhov I. V. et al. Influence of thermal pretreatment of lithium niobate plates on the characteristics of proton-exchange waveguides. Crystallography Reports, 2022, vol. 67, no. 6, pp. 980-986.
57. Sorokina I. I., Ped'ko B. B. Sposob travleniia monokristallov niobata litiya [Method for etching lithium niobate single crystals]. USSR Patent no. 1 583 479 A1, 1990. (In Russian).
58. Ped'ko B. B., Sorokina I. I., Boytsova K. N. Sposob travleniya oksidnykh nelinejnykh dielektricheskikh i opticheskikh monokristallov [Method for etching oxide nonlinear dielectric and optical single crystals]. RU Patent no. 2 558 898 C1, 2014. (In Russian).
59. Gordeeva G. V et al. Raspylenie i radiatsionnyi blistering stali 0X15N16M3B posle kholodnoi de-formatsii i rekristallizatsii [Sputtering and radiation blistering of steel 0X15H16M3E after cold deformation and recrystallization]. Atomic Energy, 1990, vol. 68, no. 4, pp. 244-247. (In Russian).
60. Sigmund P. Theory of sputtering. I. Sputtering yield of amorphous and polycrystalline target. Physical Review, 1969, vol. 184, no. 2, pp. 383.
61. Novikov I.I. Defekty kristallicheskogo stroeniia metallov [Defects in the crystal structure of metals]. Moscow: Metallurgiia Publ., 1975. 208 p. (In Russian).
62. Schlarb U., Betzler K. Refractive indices of lithium niobate as a function of wavelength and composition. Journal of Applied Physics, 1993, vol. 73, no. 7, pp. 3472-3476.
63. Atuchin V.V. Zavisimost' pokazatelei prelomleniia niobata litiia ot khimicheskogo sostava kristalla [Dependence of the refractive indices of lithium ni-obate on the chemical composition of the crystal]. Vestnik SibGUTI [Bulletin of Siberian State University of Telecommunications and Information Science], 2009, no. 3, pp. 39-45. (In Russian).
64. Tien P. K. Light waves in thin films and integrated optics. Applied Optics, 1971, vol. 10, no. 11, pp. 2395-2413.
65. Ulrich R., Torge R. Measurement of thin film parameters with a prism coupler. Applied Optics, 1973, vol. 12, no. 12, pp. 2901-2908.
66. Sutugin A. G. The Spontaneous Condensation of Vapour and the Formation of Condensation Aerosol. Russian Chemical Reviews, 1969, vol. 38, no. 1, pp. 79-93.
67. Nassau K., Levinstein H. J., Loiacono G. M. Ferroelectric lithium niobate. 1. Growth, domain structure, dislocations and etching. Journal of Physics and Chemistry of Solids, 1966, vol. 27, no. 6-7, pp. 983-988.
68. Wilson R. G. et al. Ion implantation in LiNbO3. Defect Properties and Processing of High-technology Nonmetallic Materials: Symposium held November
1983 in Boston, Massachusetts, USA. 1984, vol. 24, p. 181.
69. Van Buren H. G. Imperfections in Crystals. Amsterdam: Horth-Holland, 1961. 676 p.
70. Takigawa R. et al. Residual stress in lithium niobate film layer of LNOI/Si hybrid wafer fabricated using low-temperature bonding method. Micromachines, 2019, vol. 10, no. 2, p. 136.
Просьба ссылаться на эту статью в русскоязычных источниках следующим образом: Козлов А. А., Салгаева У. О., Петухов И. В., Журавлев В. А., Волынцев А. Б. Исследование процесса вторичной конденсации LiF в ходе реактивного ионного травления тонкопленочного ниобата лития во фтор-содержащей плазме // Вестник Пермского университета. Физика. 2024. № 2. С. 33-55. doi: 10.17072/19943598-2024-2-33-55
Please cite this article in English as:
Kozlov A. A., Salgaeva U. O., Petukhov I. V., Zhuravlev V. A., Volyntsev A. B. The investigation of the LiF redeposition process during reactive ion etching of thin-film lithium niobate with fluorine plasma. Bulletin of Perm University. Physics, 2024, no. 2, pp. 33-55. doi: 10.17072/1994-3598-2024-1-33-55
Сведения об авторах
1. Козлов Андрей Андреевич, аспирант кафедры нанотехнологий и микросистемной техники, Пермский государственный национальный исследовательский университет, ул. Букирева, 15, Пермь, 614068
2. Салгаева Ульяна Олеговна, кандидат физ.-мат. наук, научный сотрудник кафедры нанотехнологий и микросистемной техники, Пермский государственный национальный исследовательский университет, ул. Букирева, 15, Пермь, 614068
3. Петухов Игорь Валентинович, кандидат хим. наук, доцент кафедры физической химии, Пермский государственный национальный исследовательский университет, ул. Букирева, 15, Пермь, 614990
4. Журавлев Владимир Андреевич, магистр кафедры нанотехнологий и микросистемной техники, Пермский государственный национальный исследовательский университет, ул. Букирева, 15, Пермь, 614068
5. Волынцев Анатолий Борисовчи, доктор физ.-мат. наук, профессор, заведующий кафедрой нанотехно-логий и микросистемной техники, Пермский государственный национальный исследовательский университет, ул. Букирева, 15, Пермь, 614068
Author information
1. Andrei A. Kozlov, Postgraduate Student, Department of Nanotechnologies and Microsystems Engineering, Perm State University; 15, Bukireva st., Perm, 614068, Russia.
2. Uliana O. Salgaeva, Candidate of Physical and Mathematical Sciences, Research Assistant, Department of Nanotechnologies and Microsystems Engineering, Perm State University; 15, Bukireva st., Perm, 614068, Russia.
3. Igor V. Petukhov, Candidate of Chemical Sciences, Associate Professor, Department of Physical Chemistry, Perm State University; 15, Bukireva st., Perm, 614068, Russia.
4. Vladimir A. Zhuravlev, Master's Student, Department of Nanotechnologies and Microsystems Engineering, Perm State University; 15, Bukireva st., Perm, 614068, Russia.
5. Anatoliy B. Volyntsev, Doctor of Physical and Mathematical Sciences, Professor, Head of the Department of Nanotechnologies and Microsystems Engineering, Perm State University; 15, Bukireva st., Perm, 614068, Russia.