CC BY
11УДК 621.315.592:546.286
Т.Н. Смыслова, С.В. Ивин, Л.К. Орлов
АДСОРБЦИОННЫЕ И КИНЕТИЧЕСКИЕ ХАРАКТЕРИСТИКИ МОЛЕКУЛ
МОНОСИЛАНА НА РОСТОВОЙ ПОВЕРХНОСТИ SI(100) В УСЛОВИЯХ ФОРМИРОВАНИЯ ДВОЙНЫХ СВЯЗЕЙ РАДИКАЛОВ С ПОВЕРХНОСТЬЮ
Институт физики микроструктур РАН
Получены аналитические выражения, связывающие скорость встраивания атомов кремния в растущий кристалл с характерной частотой пиролиза молекул моносилана на поверхности кремния в диапазоне температур, соответствующих эпитаксиальному росту пленок Si(100).
Ключевые слова: моносилан, поверхность, скорость роста, пиролиз, частота диссоциации.
Решение проблемы получения высококачественных приборных кремниевых структур при гомо- и гетероэпитаксии в методе вакуумной химической эпитаксии требует осуществления в процессе роста прецезионного контроля процессов, протекающих на ростовой поверхности. Разработка соответствующих методов контроля, однако, невозможна без детального понимания физико-химических явлений, протекающих на эпитаксиальной поверхности роста. Необходимый минимум знаний предполагает наличие знания о конкретных путях распада молекул, т.е. типе, скоростях и энергиях активации реакций, протекающих на поверхности, значениях концентраций и временах жизни продуктов пиролиза на горячей поверхности.
Ясное понимание кинетики поверхностных реакций и знание соответствующих кинетических коэффициентов позволяет проводить расчет скоростей протекания физико-химических процессов и анализ концентраций продуктов пиролиза на ростовой поверхности. Данные, полученные в результате проведения кинетического анализа и расчета энергетического спектра радикалов молекул, адсорбируемых поверхностью, можно использовать как для прогнозирования толщины и состава получаемых в процессе эпитаксиального роста слоев, так и с целью повышения скорости протекания поверхностных физико-химических реакций, а соответственно и скорости роста пленок. Последнее является весьма актуальным при проведении низкотемпературного ростового процесса в методе вакуумной эпитаксии из газовой фазы.
Многочисленные исследования особенностей ростового гидридного процесса, проведенные разными методами, показали, что значительное уменьшение скорости роста пленок кремния при пониженной температуре эпитаксии связано исключительно с эффектом заполнения поверхностных связей продуктами распада гидридов, в частности, водородом. В конечном счете, падение скорости роста пленок связано с особенностями пиролиза молекул гидридов на ростовой поверхности. Разрабатываемые кинетические модели позволяют описывать физико-химические процессы на ростовой поверхности и определять области допустимых значений параметров, отвечающих за рост пленок. Развиваемые модели, основываясь на использовании конкретных данные технологических экспериментов, позволяют уточнять, а в ряде случаев, и определить значения кинетических коэффициентов, например, таких как скорость распада молекул, префакторы и энергии активации процессов десорбции и адсорбции, коэффициенты кристаллизации, имеющих фундаментальное значение.
Далее мы обсудим наиболее часто рассматриваемые схемы распада моносилана, базирующиеся, в общем случае, на решении системы уравнений баланса [1, 2], но, в отличие от ранее выполненных работ [2-6], рассмотрим ситуацию, когда адсорбируемые атомы и молекулы образуют не одиночные, как предполагалось ранее, а двойные связи с поверхностью.
© Смыслова Т.Н., Ивин С.В., Орлов Л.К., 2010.
Основываясь на структуре поверхностных связей атомов в кристалле, можно предположить, что рассматриваемая в работе ситуация вполне реализуема на поверхности кремния, вырезанном в направлении (100), где каждый атом кремния на поверхности имеет две свободные болтающиеся связи [7].
Развиваемая кинетическая модель процессов роста базируется на общепринятых представлениях и включает в себя различные стадии ростового процесса в условиях молекулярного течения газов: адсорбцию молекул ростовой поверхностью, их пиролиз и процессы десорбции, процессы кристаллизации и сегрегации. Введение временных производных в исходную систему уравнений позволяет рассматривать нестационарную кинетику роста и исследовать динамику протекания физико-химических процессов на ростовой поверхности. Развиваемые модели позволяют активно использовать данные технологических экспериментов и реальные технологические зависимости для вычисления кинетических коэффициентов, таких как частоты распада радикалов молекул, адсорбируемых поверхностью, энергии активации процессов пиролиза и десорбции, коэффициенты кристаллизации.
Решение поставленной ранее задачи проводилось путем решения системы кинетических уравнений, связывающих безразмерные концентрации 9; атомов и молекул, адсорбируемых ростовой поверхностью, с атомарными и молекулярными потоками веществ, поступающих к подложке, и с температурой последней. Для достижения поставленных целей наибольший интерес представляет определение поверхностных концентраций продуктов распада гидридов, а также определение либо уточнение констант, входящих в кинетические уравнения и ответственных за скорость протекания того или иного процесса на поверхности роста. Окончательный выбор рабочей системы уравнений, используемой для анализа процесса пиролиза и кристаллизации, в общем случае определяется составом либо видом используемых газов, а также выбором конкретной схемы их распада и способа десорбции водорода с поверхности. Общая система уравнений, описывающая кинетику распада гидридов на поверхности эпитаксиальной пленки, подчиняется общей схеме, по которой записываются кинетические уравнения в физической химии [8]. В самом общем случае система содержит большое число уравнений, задаваемых типом поверхностных реакций с чаще всего неизвестными кинетическими коэффициентами, в свою очередь, требующих доопределения методами квантовой химии либо непосредственно из эксперимента.
Возможные обсуждаемые в литературе схемы пиролиза молекул силана на ростовой поверхности рассматривались в работах [1, 3, 5]. Далее мы рассмотрим только простейшие модели распада моносилана. Наиболее часто при анализе ростовых процессов предполагается, что лимитирующим время распада молекулы процессом является пиролиз молекулы силена БШз [2, 4, 6], вследствие чего всеми остальными процессами и взаимодействиями на поверхности роста пренебрегается.
Если основным радикалом силана, адсорбируемым поверхностью, выступает молекула Б1Н2 с двумя свободными ковалентными связями, то она вполне может находиться на поверхности кристалла, занимая не одну, как рассматривалось ранее [1-6], а две позиции [7]. В настоящем работе это предположение выбрано в качестве рабочей гипотизы, на основе которой записана система кинетических уравнений. Последняя строится на базе следующих физико-химических реакций для ] = 2:
^Н^) + (к+2) ^ §1Н2 + кН + ((2-к)/2)И2Т, ЖН + п ^ & + пН + ((2-п)/2)Н2Т, (1)
где символ Б! (БМ?) обозначает адатом кремния (радикал молекулы силана), связываемый с поверхностью двумя ковалентными связями, индексы к, п = 0 ^ 2 - показывают, на какой ступени процесса пиролиза и с какой эффективностью происходит захват водорода поверхностью. Соответствующая система кинетических уравнений для безразмерных поверхностных концентраций продуктов в стационарном пределе в этом случае принимает вид
0 = (к+2)!(8^зщ4/п8)(0пы)к+2 - (п+2)!у§ш2(е§щ2)2(е&)п, 0пы = 6& - 2е§1и2, 0 = k•(k+2)!(Ss1н4Fs1н4/ns)(еnbl)k+2 + п^п +2)^Ш2^щ2)2(е&)п - х^еи , (2)
о = (п+2)^ш2^ш2)2(е&)п - 2^(е^)2, 2е^ш+2еs1+ен + е& = 1,
где еS1H2, еS1, ен - безразмерные (отнесенные к концентрации поверхностных связей пз = 6.78х 1014 см-2) концентрации радикалов SiH2, адатомов кремния и водорода на поверхности роста; у^н2 - скорости распада радикалов силана на поверхности, у - коэффициент десорбции водорода; г® - коэффициент встраивания атомов кремния в решетку кристалла;, еы , епЬ1 - концентрация блокированных и неблокированных свободных е& связей. Фактори-альные множители соответствуют числу перестановок, то есть числу возможных конфигураций мест, на которые можно посадить адсорбируемые молекулы (атомы).
Исключая из системы уравнений (2) поверхностные концентрации, нетрудно найти соотношение, связывающее частоту распада радикала моносилана со скоростью кристаллизации г&:
VS1H2 = 2п+2 а / |(п+2)![у - (а/2rsl) 1/2]2 [у + 2р1/(к+2) - 3(a/2rsl)1/2]n}. (3)
Коффициенты в (3) равны:
а = /1о; Р = Впы(к+2) = а/((к+2)! Бвш4 Р^щ); у = (1 - р1/(к+2) - 0Я), (4)
где УёГ (А/с), 10 = 1.23 А. Безразмерные значения концентрации продуктов распада и плотность свободных поверхностных связей соответственно равны
= (а/2^)1/2; еsш2 =(у - (а/2^)05)/4; 0& = 1 - 2esш2 - 29^ - 0н. (5)
Из уравнения (3) нетрудно установить асимптотики поведения кривой у^н2(г&). В реальной ситуации коэффициент кристаллизации велик. Если стремится к ю, то нетрудно найти минимальное значение частоты распада молекул:
VslH2 = 2п+2а/{(п+2)! у2 [у + 2р1/(к+2)]п}. (6)
Величина VS1H2 является функцией всех параметров задачи, включая выбор модели распада и характер заполнения поверхности водородом. Формула (6) показывает зависимость скорости роста пленки от величины VS1H2 и позволяет оценить значение скорости распада моносилана на поверхности роста в условиях высокой скорости встраивания в решетку кристалла дифффундирующих по поверхности адатомов кремния.
Чтобы понять, какие решения реально допустимы в эксперименте, рассмотрим характерный вид зависимостей коэффициента встраивания и безразмерных поверхностных концентраций атомов кремния и радикалов молекулы силана от скорости ее распада для ряда наиболее интересных моделей системы. Расчет проведем, используя формулы (3,5) при температуре роста пленки 450°С для различных значений параметров (/=2,к,п), соответствующих разным способам захвата поверхностью высвобождающегося из молекулы водорода. Ситуация с (к,п) = (2,0) описывает адсорбцию водорода на первой стадии процесса пиролиза, ситуация с (к,п) = (0,2) - адсорбцию водорода на второй стадии процесса. Ситуация с (к,п) = (1,1) соответствует тому, что атомы водорода в результате разрыва связей в молекуле не сразу переходят в свободное состояние, а сначала переходят в связанное состояние на поверхности, а затем в процессе диффузии они объединяются в молекулу водорода, десорби-руя с поверхности.
Такая ситуация может оказаться более предпочтительной, так как энергия указанного перехода атома водорода существенно меньше энергии связи его в молекуле. Решение, соот-
ветствующее захвату поверхности всех атомов водорода (к,п) = (2,2), отсутствует. Характерный вид зависимостей г&(у08Щ2) представлен на рис. 1, а. На рис. 1, б, в приведены соответствующие значения поверхностных концентраций 98; и 98щ2. При высоких значениях коэффициента встраивания, соответствующего низкому уровню заполнения поверхностных связей адатомами кремния, а именно такая ситуация наиболее близко соответствует реальному эксперименту.
в)
Рис. 1. Зависимость коэффициента встраивания ^ (я), поверхностных концентраций в81 (б) и в8щ2 (в) от префактора частоты распада молекул моносилана для моделей 211, 220 и 202
при температуре роста Т^ = 45°°С, = 0.007 и давлении газа в реакторе Р8щ4 = 0.3 мТор
Нетрудно оценить значение скорости распада молекул силана на поверхности роста слоя кремния у8щ2 = у°8щ2 ехр(-Еа/кТёГ) для Еа = 0.08 эВ [10].
Выражение (3) при заданной скорости роста пленки и степени ее покрытия водородом устанавливает однозначным образом взаимосвязь между коэффициентом кристаллизации и скоростью распада моногидрида на поверхности слоя. Знание значений параметров и легко позволяет найти поверхностные концентрации атомов кремния бзь радикалов молекулы силана 08щ2 и концентрацию свободных связей на поверхности 0^.
На рис. 2 на примере одной модели (2,0,2) показаны характерные зависимости скорости встраивания и поверхностных концентраций 98ь 98щ2 от скорости пиролиза молекул для более широкого диапазона температур. Области допустимых значений концентраций 0 < вэь 98щ2, 9& < 1 определяют нижний предел значения частоты распада молекул силана. Это значение, как видно из рис. 1 и рис. 2, зависит как от температуры роста, так и от характера
протекания реакций на поверхности, описываемого той или иной кинетической моделью. Из приведенных рисунков можно оценить минимальное значение, которое принимает префак-тор скорости распада молекул на поверхности роста пленки при заданных значениях параметров задачи. Из рис. 1 следует, что при низкой (450°С) температуре пиролиза префактор скорости распада молекул лежит в диапазоне 40-60 с-1. Скорость распада достаточно велика
Л! 1
и по порядку величины близка к значению Узщ2 ~ 10с- . С повышением температуры роста до 600-700°С значение префактора уменьшается до значений порядка 1, скорость распада падает до нескольких долей секунды.
0,15-
I
с/э
Ф
0,10-
0,05-
С
3
10
100 1000
0 „-1 V , С
в)
Рис. 2. Зависимость коэффициента встраивания г81 (я), поверхностных концентраций 081 (б) и 68щ2 (в) от префактора частоты распада молекул моносилана у°81ш для модели 202 и значений температур роста Т^(°С) = 600 (1), 550 (2), 500 (3) при £81н4 = 0.007 и давлении газа
в реакторе ^¡га = 0.3 мТор
Таблица 1
Значения скорости роста пленок кремния из моносилана от температуры роста при давлении в реакторе ^¡га = 3 мТор и соответствующая данным температурам степень покрытия поверхности 81 адсорбируемым водородом
т °с 600 550 500 450
А/с 0.237 0.136 0.045 0.009
бы = Пы П 0.091 0.253 0.544 0.839
Таким образом, расчеты показывают, что температурная зависимость скорости распада молекул моносилана определяется не активационным множителем, а температурной зависимостью префактора v°síh2. Необходимо также отметить (рис. 1) сильную зависимость частоты распада молекул от того, на какой ступени процесса распада происходит захват атомов водорода с молекулы поверхностью. Распад происходит с наибольшей скоростью, если молекулы водорода адсорбируются на первой ступени пиролиза (кривая 1 на рис. 1). Скорость пиролиза при данной температуре роста минимальна, если молекулы водорода адсорбируются только на второй ступени пиролиза (кривая 3 на рис. 1).
Из формулы (3) можно также оценить минимальное значение, которое может принимать коэффициент кристаллизации rSi для подстилающей поверхности, характеризуемой низким коэффициентом растворения атомов кремния. Предполагая скорость распада силана vSiH2 большой величиной, для коэффициента кристаллизации rSi получаем значение rSi = а/2у2. Второе решение всегда является лишним, не соответствующим реальным значениям поверхностных концентраций. Минимально допустимые значения коэффициента встраивания атомов кремния, соответствующие предельно высокой скорости распада молекул силана на поверхности роста, зависят только от параметров (k, n), то есть от характера десорбции водорода с поверхности в диапазоне низких температур роста.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований по проектам № 08-02-00065 и № 08-02-97017p - "Поволжье ".
Библиографический список
1. Greve, D.W. Materials Science & Engineering. B18, 22 (1993).
2. Potapov, A.V. Thin Solid Films / A.V. Potapov, L.K. Orlov, S.V. Ivin. - 1999. V. 336. - P. 191-196 (1999).
3. Орлов, Л.К. Поверхность, Рентгеновские, синхротронные и нейтронные исследования / Л.К. Орлов, Т.Н. Смыслова. - 2005. №11. С. 99-105.
4. Орлов, Л.К. Физика и техника полупроводников / Л.К. Орлов, Т.Н. Смыслова. 2005. № 39. С. 1320-1324.
5. Орлов, Л.К. Физика и техника полупроводников / Л.К. Орлов, Т.Н Смыслова. 2006. Т. 40. № 45. С.1427-1432.
6. Орлов, Л.К. Физика и техника полупроводников / Л.К. Орлов, Т.Н. Смыслова. 2007. Т. 41. № 56. С.1517-1522.
7. Зенгуил, Э. Физика поверхности / Э. Зенгуил. - М.: Мир. 1990. - 536 с.
8. Товбин, Ю.К. Теория физико-химических процессов на границе газ-твердое тело / Ю.К. Тов-бин. - М. Наука. 1990. - 287с.
9. Appl. Phys.Lett / K.J. Kim [at all.] 1993. V. 62. P. 3461. - 3465.
1°. Vittadini, A. Phys. Rev. Lett. Ser. A / A. Vittadini // Selloni.1995. V. 75. P. 4756-4759.
Дата поступления в редакцию 07.07.2010
T.N. Smislova, S.V. Ivin, L.K. Orlov
ADSORPTION AND KINETIC CHARACTERISTICS OF MONOSILAME MOLECULES ON GROWING SURFACE SI (100) AT DOUBLE BOUNDING FORMATION BETWEEN
RADICALS AND SERFACE
The analytic forms linking the silicon atom embedding rate with the typical frequency of silane molecular py-rolysis on the silicon surface were obtained. The temperature range corresponds to the Si(100) films epitaxial growth.
Keywords: monosilane, surface, growth rate, pyrolysis, dissociation frequency.
