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

Язык труда и переводы:
УДК:
538.971
Дата публикации:
23 ноября 2022, 13:28
Категория:
Математическое моделирование физических процессов и технических систем
Авторы
Аннотация:
Программный код открытого доступа NASCAM, реализующий метод статистического моделирования роста покрытий, задействован для прогнозирования структурных характеристик металлических покрытий (на примере алюминия), осаждаемых при работе магнетронной распылительной системы с жидкофазной мишенью. Роль факторов высокой плотности потока осаждаемых частиц и большого количества испаренных атомов в осаждаемом потоке по сравнению с традиционным случаем магнетронного распыления изучена. Получены особенности формирования шероховатости, пористости и плотности моделируемых пленок.
Ключевые слова:
статистическое моделирование, магнетронные распылительные системы, тонкие пленки, структура, морфология
Основной текст труда

Введение

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

v_{dep}=v_{ev}+v_{sput},

где v_{ev} и v_{sput} – составляющие скорости осаждения за счет испарения и распыления соответственно.

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

Методы

В качестве инструмента моделирования использовался пакет NASCAM [3]. В программном коде NASCAM реализован метод статистического моделирования, что позволяет рассчитывать системы, численность которых достигает сотен тысяч атомов [4].

Входные параметры программы подбирали на основе эксперимента [2] с жидкофазной магнетронной распылительной системой. Данные эксперимента представлены ниже:

режим работы… МРС со среднечастотным источником питания

осаждаемое вещество… Al

мощность МРС, Вт… 3500

рабочее напряжение U, В… 558

температура мишени T, K… 1493

рабочий газ… самораспыление

рабочее давление, Па… 0,18–0,2

скорость осаждения пленки v_{dep} , нм/c… 50,6 

скорость испарения v_{ev} , нм/c… 42,8 

скорость  распыления v_{sput} , нм/c… 7,8 

Расчет скорости осаждения, а также величин потоков распыляемых и испаряемых частиц, происходил на основе моделей, представленных в [5].

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

Результат расчета энергетического спектра представлен на рис. 1.

Рис. 1. Энергетический спектр осаждаемых атомов меди при распылении испаряющейся мишени МРС

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

При магнетронном осаждении для уплотнения структуры пленок для улучшения их функциональных характеристик используют ионную бомбардировку растущей пленки. Для реализации этого приема на подложку подается отрицательный потенциал U_{bias} , называемый потенциалом смещения, в результате чего часть ионов из плазмы магнетронного разряда вытягивается из плазмы и ускоряется в сторону подложки. Доля таких ионов в общем потоке, согласно экспериментальным наблюдениям, варьируется в пределах:

{\frac {F_{ion}}{F_{dep}}}=0,1\ldots 1.

Этот эффект также был учтен при моделировании. Доля ионов в общем потоке была задана, равная  0,1F_{dep} , где F_{dep}  — суммарная плотность потока осаждаемых частиц. Энергия ионов по абсолютному значению в эВ примерно равна напряжению смещения U_{bias} , выраженному в В.

Для оценки влияния скорости осаждения и энергетических спектров частиц на свойства получаемых покрытий были получены результаты для четырех режимов работы МРС:

  • А — испарение и распыление; высокая скорость осаждения;
  • Б — распыление; низкая скорость осаждения;
  • В — испарение; высокая скорость осаждения;
  • Г — испарение; низкая скорость осаждения.

Характеристики потоков частиц, осаждаемых в этих режимах, задавались исходя из экспериментальных данных, представленных выше. Данные для моделирования режимов А–Г содержатся в таблице.

 Входные данные для моделирования режимов А–Г

Режим

vdep, нм/с

vdep, МС/c

Доля испаряемой компоненты в общем потоке wev

А

50,6

224,6

0,846

Б

7,8

34,7

0

В

42,8

189,9

1

Г

7,8

34,7

1

Монослоем [МС] называется один плотноупакованный слой атомов. Зная атомную плотность вещества n [шт./м3] и толщину покрытия [м], можно рассчитать количество плотноупакованных слоев в этом покрытии:

{\text{[МC]}}={\text{[м]}}\cdot n^{1/3}.

Поскольку vdep состоит из двух компонент, эти компоненты можно представить только значением v_{dep} и долей одной из компонент в общем потоке w_{ev} . Таким образом

v_{ev}=v_{dep}\cdot w_{ev};

v_{sput}=v_{dep}\cdot \left(1-w_{ev}\right).

При исследовании влияния напряжения смещения на свойства тонких пленок сравнивались показатели шероховатости и пористости покрытий в двух режимах осаждения:

  • 2.А) U_{bias} = 0;
  • 2.Б) U_{bias}  = –100 В.

Осаждаемый поток для режимов 2.А и 2.Б был выбран на основе режима А.

Другие входные данные для моделирования в NASCAM были заданы следующими:

  • размер подложки 25×100 атомов;
  • количество моделируемых осаждаемых частиц 200 000 шт.;
  • температура подложки 25 °C;
  • материал подложки Cu;
  • кубический тип решетки;
  • ориентация подложки к потоку осаждаемых частиц – (100).

Результаты и обсуждение

На рис. 2 представлена зависимость шероховатости покрытия от количества осажденных атомов для режимов осаждения А–Г.

Рис. 2. Зависимость шероховатости пленки от количества осаждаемых атомов

Шероховатость выше при более низкой скорости осаждения. Режимы, при которых средняя энергия осаждаемых частиц выше, т. е. режимы А и Б, характеризуются более высокими значениями шероховатости. Рассчитанная шероховатость коррелирует с экспериментальными значениями [7].

Рис. 3. Изображение поверхностей осажденных пленок (на шкале указана высота покрытия от уровня подложки в единицах постоянных решетки): l.u. — постоянная решетки

Самый маленький характерный размер микроструктур на поверхности соответствует низкоскоростному термическому осаждению. Наибольший размер микроструктур получен для МРС с холодной мишенью.

Изменение шероховатости образцов в течение процесса осаждения покрытия при наличии потенциала смещения (режимы 2.А–2.Б) представлено на рис. 4.

Рис. 4. Зависимость шероховатости пленки от времени осаждения

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

Рис. 5. Изображение поверхностей осажденных пленок: l.u. — постоянная решетки

Количество микроструктур отрицательно коррелирует с шероховатостью поверхности. В режиме (2.Б) наблюдается их наибольшее количество. Распределение пор по размерам представлено на рис. 6.

Рис. 6. Распределение пор в покрытии согласно их размерам

Режим (2.Б) отличается плотноупакованной структурой и малым количеством пор больших размеров. С увеличением доли высокоэнергетических ионов, попадающих на подложку, граница образования открытых пор смещается в сторону поверхности покрытия. Численные значения плотности покрытия в относительных единицах представлены на рис. 7.

Рис. 7. Эволюция плотности покрытия

Заключение

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

Для систем осаждения с высокой скоростью осаждения и широким спектром энергетического распределения частиц характерны значения шероховатости ниже по сравнению со стандартными МРС с холодными мишенями.

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

Литература
  1. Sidelev D.V., et al. High-rate magnetron sputtering with hot target. Surface and Coatings Technology, 2016, vol. 308, pp. 168–173.
  2. Блейхер Г.А., Юрьева А.В., Сиделёв Д.В., Юрьев Ю.Н., Грудинин В.А. Исследование фундаментальных аспектов режима самораспыления при работе магнетронных распылительных систем с испаряющимися мишенями и его влияние на свойства осаждаемых покрытий. Итоговый отчет по НИР. Рег. № НИОКР AAAA-A18-118011090005-2. Томск, ТПУ, 2019, 39 с.
  3. NASCAM (NAnoSCAle Modeling). Université de Namur. URL: https://www.unamur.be/sciences/physique/ur/larn/logiciels/nascam (дата обращения 24.10.2022).
  4. Lucas S., Nita F., Moskovkin P. NASCAM user’s manual for version 4.6.2. University of Namur, 2018, 62 p.
  5. Блейхер Г.А., Кривобоков В.П. Эрозия поверхности твердого тела под действием мощных пучков заряженных частиц. Новосибирск, Наука, 2014, 248 с.
  6. Ekpe S.D., Dew S.K. Theoretical and experimental determination of the energy flux during magnetron sputter deposition onto an unbiased substrate. Journal of Vacuum Science & Technology A: Vacuum, Surfaces, and Films, 2003, vol. 21, no. 2, pp. 476–483.
  7. Bleykher G.A., et al. The properties of Cu films deposited by high rate magnetron sputtering from a liquid target. Vacuum, 2019, vol. 169, art. no. 108914.
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.