Доклады БГУИР
2011 № 8 (62) УДК 621.315.593:539.216.2
УСЛОВИЯ ОБРАЗОВАНИЯ НАНОКЛАСТЕРОВ ПРИ ФОРМИРОВАНИИ НАНОСТРУКТУРИРОВАННЫХ ПЛЕНОК ПОЛИКРИСТАЛЛИЧЕСКОГО КРЕМНИЯ, ЛЕГИРОВАННОГО ГЕРМАНИЕМ
Н С. СТРОГОВА, А С. СТРОГОВА, А.А. КОВАЛЕВСКИЙ, В.М. БОРИСЕВИЧ, Д.В. ПЛЯКИН
Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 9 июня 2011
Проведены исследования топографии поверхности, кристаллической структуры пленок, поликристаллического кремния, легированного германием, полученных разложением смеси моносилан - моногерман - водород. Изменение структуры Ge и SiGe-кластеров анализировалось на основании данных спектроскопии комбинационного рассеяния света. Также проанализированы зависимости роста пленок от полного давления в реакторе и парциального давления моносилана.
Ключевые слова: нанокластеры, формирование наноструктурированных пленок кремния, самоорганизация нанокластеров.
Введение
Наноструктурированные пленки поликристаллического кремния, легированного германием, вызывают интерес с точки зрения возможности создания нанокластеров германия и твердого раствора кремний-германий, которые важны для приборных применений [1].
Методика эксперимента
Исследуемые образцы были получены методом пиролиза моногидридов Si и Ge в установке «Изотрон 4-150» при общем давлении газовой смеси 35-40 Па и объемном соотношении моногермана к моносилану 0,001-0,003. Температура осаждения при этом не превышала 680°С.
В качестве подложек использовались кремниевые пластины КЭФ 4,5, ориентированные в плоскости (100). На поверхность пластин наносился тонкий оксид кремния в процессе окисления в сухом кислороде, а оксиды диспрозия и иттрия формировались в процессе их нанесения и последующего окисления. Температура образцов в процессе осаждения контролировалась с помощью платино - платино - родиневой термопары.
Исследования топографии поверхности полученных структур проводились с помощью атомно-силового микроскопа Ntegra Prima. В качестве зондов применялись кремниевые канти-леверы с радиусом закругления вершин зондирующих пирамид 15-20 нм (SCNC12,NT-MDT), а также вискеры, отличительной особенностью которых является завершение зондирующей пирамиды вытянутым узким и более острым конусом из материала на основе углерода (NCC05, NT-MDT). Применение вискеров заметно улучшало контраст изображений островков Ge и SiGe. Кроме того, с помощью атомного силового микроскопа определялись плотность, размеры и форма кластеров.
Изменение структуры Ge и SiGe-кластеров анализировалось на основании данных спектроскопии комбинационного рассеяния света (КРС). Спектры КРС на оптических фононах
получали при облучении структур Ar лазером с длиной волны X = 514,5 нм и регистрировали с помощью спектрометра ДФС-52. Спектры регистрировались при комнатной температуре.
Результаты и их обсуждение
Проведены исследования топографии поверхности, кристаллической структуры пленок поликристаллического кремния, легированного германием (ППК [Ge]), полученных разложением смеси моносилан - моногерман - водород. Этот процесс в сочетании с пост окислением рассматривался в качестве альтернативной, более «чистой» технологии, позволяющей сформировать нанокластеры германия или сплава кремний-германий различной структуры и размеров. Исследовано влияние температуры на скорость осаждения пленки кремния, легированного германием и давления остаточных газов в ростовой камере установки «Изотрон 4-150» на морфологию поверхности тонких (40-400 нм) пленок кремния, осажденных при низких (540-670°С) температурах подложки с различной структурой. Установлено, что размер нано- и микрокристаллитов кремния, германия и твердого раствора кремний-германий существенно снижается с понижением температуры подложки в процессе осаждения и достигает ~ 10 нм при Тп = 540°С. При изменении общего (полного) давления в реакторе от 40 Па до 140 Па скорость осаждения пленок изменяется от 4,5 до 13,1 нм/мин, а размер нанокластеров уменьшается с 10 нм до 2,5 нм (рис. 1), причем по линейному закону. Разброс по размеру нанокластеров на пластинах по всей изотермической зоне реактора и по площади самих пластин, при зазоре между пластинами (И) 5 мм, не превышал ±15%.
V,, им/мин 13,0
12,0 11,0
10,0 9,0 8,0 7.0 6,0 5,0
30 40 50 60 70 80 90 100 110 120 130 140 150 160 Plot, Пэ
Рис. 1. Зависимость скорости G роста пленки поликристаллического кремния в LPCVD-процессе от полного давления в реакторе Plot
Исследование зависимости размера нанокластеров от величины общего давления в камере при скорости потока моносилана 7,0 + 0,5 дм3/час, при температуре 620 °С (893 К) показало, что размер нанокластеров германия и сплава кремний-германий при осаждении их в режиме выращивания ППК [Ge] сильно зависит от величины давления в реакторе.
Распределение нанокластеров по размеру сформированных в режиме осаждения ППК [Ge], при парциальном давлении SiH4 35 Па и 45 Па для скорости потока 7,0 дм3/ч и 10,5 м3/ч вдоль диаметра подложки (рис. 2) показывает, что радиальная неоднородность по размеру снижается при уменьшении давления. Если разброс по размеру нанокластеров, сформированных при давлении 45 Па, составляет + (12,5 - 13,0)%, то аналогичная величина для нанокластеров, сформированных при давлении 35 Па, снижается до + (7,5 - 10)% (см. рис. 2).
При осаждении пленок ПКК при общем давлении в реакторе 20 - 50 Па скорость осаждения (Vp) крайне мало зависит от давления в реакторе (120 - 100 А/мин). В это же время Vp сильно зависит от температуры осаждения (7П) и величины скорости потока и парциального давления SiH4. Это объясняется тем, что лимитирующей стадией процесса является скорость разложения моногидрида (МГ) на поверхности подложки. Уход в сторону более низких величин давлений в реакторе, чем 50 Па ощутимо влияет на структуру и морфологию пленок. С уменьшением общего давления при осаждении слоев интенсивность поверхностной диффузии адсорбированных атомов повышается.
0 2 4 6 8 10 1 2 1 4 1 6 18 20 22 24 26 28 30 32 34 36 38 40 42 Па
Рис. 2. Зависимость скорости Ур роста пленки поликристаллического кремния в LPCVD-процессе от парциального давления силана Р,ш
С уменьшением давления поверхность растущей пленки становится более «чистой» и более зеркальной, что обусловлено снижением количества адсорбированных на ней неконтролируемых примесей и побочных продуктов реакции, а основной процесс идет по гетерогенному механизму.
Перевод процесса осаждения слоев поликристаллического кремния в область, как можно более низких давлений МГ, представляет значительный практический интерес, поскольку способствует повышению стабильности процесса осаждения за счет перевода его в область, где Ур практически нечувствительна или малочувствительна к величине давления в объеме изотермической зоны. Более того, в таких условиях имеется возможность по управлению как скоростью процесса, так и структурно-морфологическими свойствами осажденных слоев.
Повышение однородности слоев по толщине, а, следовательно, должно быть и по электрофизическим параметрам, при уменьшении общего давления в изотермической зоне реактора при одной и той же величине скорости потока газов, обусловлено в первую очередь тем, что вероятность протекания реакций в газовой фазе, т.е. гомогенных, становится достаточно малой и, кроме того, практически исключаются полимеризация МГ и их производных с образованием кластеров кремния в газовой фазе, поскольку концентрация таких производных как 81Н и GeH4 пропорциональна квадрату парциального давления соответственно 81Н4 и GeH4.
Учитывая близость свойств 81Н с GeH4, интерес представляет именно комбинация этих компонентов. Причем, если GeH4 использовать как всего лишь легирующий компонент, его объемы по отношению к 81Н должны быть в пределах 0,001-0,003 об. частей.
Скорость роста пленок поликристаллического кремния, легированного германием, сильно зависит от величины парциального давления и скорости потока МГ в зоне реакции. И действительно выяснено, что изменение парциального давления 81Н в реакторе от 10 Па до 40 Па приводит к увеличению Ур от 50 до 110 А/мин при прочих равных условиях, а это, в свою очередь, определяет скорость формирования нанокластеров.
При скорости потока МГ и, в частности, 81Н^7 и 10 дм3/час адсорбция МГ диссоциативная, т.е. при адсорбции 81Н на поверхности подложки происходит диссоциация молекулы 81Н4 на 81Н2(ад) и Н2, GeH4 на GeH2 (ад) и Н2, каждая из которых занимает свой адсорбционный центр:
81Н (газ) ^ 81Н2(ад) + Н2(газ),
§Щ2(ад) ^ 81(Тв) + Н2(газ), GeH4(газ) >GeH2(ад) + Н2(газ), GeH2(ад) ^ Ge(тв) + Н2(газ).
(1) (2)
(3)
(4)
Если на растущей поверхности имеются радикалы 81Н2 и GeH2 или их аналоги типа 81Н3, GeH3, 81Н, GeH, они вступают во взаимодействие между собой, например:
81Н + ОеИ-
->81Н + ОеН + Н .
(5)
В итоге адсорбированными частицами является кремний и германий
Si, Ge Geад, Si, Geад^S1, Ge (тв), (6)
т.е. характер адсорбции не меняется при изменении скорости потока МГ, что подтверждает возможность сделанных предположений о возможном влиянии величины t на величину начальной вероятности «прилипания» молекул МГ, в частности БИ^ поскольку объем GeH4 ничтожно мал по отношению к объему БИ^ хотя в целом коэффициент адсорбции равен
9 = + ^е + 0Н + 98Ш2 + ^еН2 + 9пов.связей = 1.
В области низких давлений, когда Vp зависит от скорости потока МГ, модель предсказывает большую зависимость скорости осаждения V, = РРБ/Б0 от , что как показывают результаты работы [2], подтверждается на практике. Перевод процесса осаждения слоев ПКК^е) в область как можно более низких давлений МГ представляет значительный практический интерес. Это обусловлено тем, что низкие давления газовой смеси способствуют повышению стабильности процесса осаждения за счет перевода его в область, где V, практически нечувствительна или малочувствительна к величине давления в объеме изотермической зоны. Более того, в таких условиях имеется возможность по управлению как скоростью процесса, так и структурно-морфологическими свойствами осажденных слоев и, в частности, плотностью нанокла-стеров.
Повышение однородности слоев по толщине, и, следовательно, и по электрофизическим параметрам, при уменьшении общего давления в изотермической зоне реактора при одной и той же величине скорости потока газов обусловлено в первую очередь тем, что вероятность протекания реакций в газовой фазе, т.е. гомогенных, таких как (формулы 1-6) становится достаточно малой и, кроме того, практически исключаются полимеризация МГ и их производных с образованием кластеров кремния в газовой фазе, поскольку концентрация таких производных, как БЩ2 и GeH2, пропорциональна квадрату парциального давления соответственно и GeH4.
Как мы уже отмечали выше, добавки GeH4 приводят к довольно заметному повышению скорости осаждения. Добавление GeH4 к в количестве 0,002 об. частей по отношению к объему при прочих равных условиях приводит к увеличению скорости практически в два раза в центре изотермической зоны.
Эта зависимость позволила определить эффективный порядок реакции по GeH4 (п = AlgV,/AALgGeH ), величина которого п = 1,25±0,02 для температуры 893 К не изменяется по
всей зоне роста слоев.
Дробный порядок характерен для сложных химических реакций, а его положительное значение указывает на увеличение скорости реакции. В нашем случае это усиление гетерогенного маршрута процесса разложения БЩф Таким образом, при температуре 893К преобладает гетерогенный маршрут реакций. При этом гетерогенный маршрут является более эффективным поставщиком радикалов, БЩ2 для которых GeH4 является инициатором их образования. В итоге при постоянной температуре подложки, постоянных парциальных давлениях МГ и постоянных скоростях потока газов V, слоев можно рассчитать по формуле
V = N (031 +^е), (7)
"0
а содержание германия в пленке по формуле: 0
С = , (8)
" 0Ge +03, ,
где г - частота встраивания атомов в кристаллическую решетку, которая зависит от концентрации атомов кремния и германия, N - концентрация на поверхности подложки, Ы0 - концентрация в газовой фазе, - коэффициент адсорбции германия, 9Б1 - коэффициент адсорбции кремния.
Диапазон температур процесса осаждения слоев сверху ограничивается процессом преобразования структуры слоев в результате их рекристаллизации, а снизу очень малой скоростью осаждения.
Фундаментальные исследования процесса разложения SiH4, показали, что реакции его термического разложения являются многомаршрутными и многостадийными, зависящими от условий разложения. Как уже отмечалось выше, процесс разложения SiH4 и в том числе GeH4 как аналога SiH4, может протекать по двум маршрутам - гетерогенному и гомогенному, имеющим разные кинетические временные характеристики. Эти различия позволяют разделить составляющие процесса. Величины энергии активации гетерогенного (Е а(гет)) и гомогенного (Е а(гом)) процессов различны и зависят от условий разложения гидридов, и в первую очередь от температуры подложки. Зависимость скорости осаждения от температуры используется для определения энергии активации:
Е„ =-8,3 х 10"
V
2,3^ ^
' ВК
(
11
, кДж/моль
(9)
Т
К1 2
Т
1
где Е^ - энергия активации процесса, кДж/моль; V2 - скорость реакции осаждения при температуре подложки Т2; VI - скорость реакции осаждения при температуре подложки Т1.
По сути, энергия активации - это то избыточное количество энергии (по сравнению со средней величиной), которой должна обладать молекула в момент столкновения, чтобы быть способной к данному химическому взаимодействию. В большинстве своем увеличение ТП способствует усилению такого взаимодействия. При повышении ТП скорость химических реакций увеличивается весьма значительно.
Скорость реакции сильнее изменяется с температурой в тех реакциях (процессах), в которых энергия активации больше. И наоборот, если в данной реакции энергия активации незначительна, то скорость этой реакции (этого процесса) слабо изменяется с температурой. Многочисленные данные, посвященные процессу разложения SiH4 в реакторах атмосферного давления при температуре 550-900 °С свидетельствуют, что Е а(гет) ~ 200кДж/моль, а Е а(гом) ~ 230 235 кДж/моль [3-16]. При пониженных давлениях для Е а(гет) характерны величины 150170 кДж/моль, для Е а(гом)~200 кДж/моль [16-18]. По сути, результаты с использованием наших экспериментальных данных полностью согласуются с приведенными выше данными для процессов, происходящих при пониженных давлениях.
Эффективная энергия активации, рассчитанная нами по уравнению (9) с использованием экспериментальных данных (рис. 3) для случая, когда величина полного давления в реакторе 40 Па, парциального давления SiH4 35Па и скорость потока газов 7 дм3/ч, при отсутствии GeH4, составляет величину 175 кДж/моль, а при наличии GeH4 - 135кДж/моль в кинетической области и 33,5 и 28 кДж/моль, соответственно в диффузионной зоне.
Т , оС
103/Т , к-1
Рис. 3. Зависимость скорости осаждения пленок поликристаллического кремния, легированного германием, от температуры подложки: 1-2 - полное давление в реакторе 40 Па, парциальное давление моно-силана 35 Па, скорость потока газов 7,0 дм3 /ч; 3-4 - полное давление в реакторе 50 Па, парциальное давление моносилана 40 Па, 0еИ4^1И4=0,002 (1-4).
100
10
0,1
1,10
1,12
1,14
1,16
1,18
1,20
Снижение скорости потока газов в реакторе до 3,5 дм3/ч при прочих равных условиях приводит к увеличению энергии активации до 207 кДж/моль, а повышение ее до 10,5 дм3/час к снижению энергии активации до 188 кДж/моль. Таким образом, малые скорости потока газов способствуют усилению температурной зависимости V, слоев, и напротив повышение скорости потока газов снижает степень этой зависимости.
Увеличение парциального давления до 40 Па при общем давлении в реакторе 50 Па при прочих равных условиях при отсутствии GeH4 приводит к уменьшению энергии активации до 104 кДж/моль, а при наличии GeH4 до 133 кДж/моль в кинетической области и до 11,6 и 56 кДж/моль, соответственно, в диффузионной области. Анализируя эти результаты можно констатировать, что при больших, чем 35 Па, давлениях БЩ4 и общем давлении в реакторе больше 40 Па усиливается влияние наличия GeH4 на процесс осаждения с повышением температуры в изотермической зоне реактора.
Оценка эффективной энергии активации процесса разложения БЩ4 для условий наших экспериментов показала, что в отсутствие GeH4 в зоне реакции ее величина составляет 104 -207 кДж/моль. В этом случае сильный вклад на ее величину вносят парциальное давление и скорость потока газов в изотермической зоне реактора. Они являются важнейшими технологическими параметрами процесса осаждения ПКК, которые позволяют управлять как скоростью процесса осаждения, так и структурно-морфологическими свойствами осаждения пленок. При наличии GeH4 в зоне реакции в объеме 0,002 по отношению к объему при условиях: скорость газового потока 7,0 дм3/ч, парциальное давление БЩ4 35 Па и общее давление газовой смеси в реакторе 40 Па - энергия активации снижается с 175 до 137кДж/моль. А это говорит о том, что в этом случае V, слоев в меньшей степени зависит от температуры в изотермической зоне реактора. Однако увиличение парциального давления до 40 Па и выше, а общего давления до 50 Па и выше в реакторе, напротив, в присутствии GeH4 способствует увеличению эффективной энергии активации с 104 до 135 кДж/моль, т.е. приводит к большей степени зависимости процесса осаждения слоев от температуры.
К аналогичным результатам, но только в диффузионной области приводит увеличение скорости потока газов до 10,5 дм3 /ч. Повышение давления в реакторе выше 40 Па при температуре > 640 °С (913 К) приводит к тому, что появляется матовость поверхности слоев и его для этих условий следует считать критическим. Этот факт указывает на увеличение вклада гомогенного маршрута разложения гидридов. Во всех остальных случаях имел место гетерогенный маршрут разложения МГ. Это согласуется с выводами других исследователей о разложении БИ^ [17-21], что при сравнительно низких температурах разложения БИ^ преобладает гетерогенный маршрут реакций, а с повышением температуры скорости гомогенных реакций увеличиваются больше, чем скорости гетерогенных реакций. При этом гетерогенный маршрут является более эффективным поставщиком радикалов БЩ2 и GeH2; БЩ3 и GeH3; БЩ и GeH, для которых GeH4 является инициатором их образования.
Таким образом, величина скорости потока газов, величина общего давления в реакторе, величина парциального давления БЩ4 в реакторе и величина объемного соотношения GeH4/SiH4 как и ТП в изотермической зоне реактора (камеры), являются важнейшими технологическими параметрами процесса осаждения ультратонких пленок ПКК, которые позволяют управлять V
Заключение
Общее давление компонентов газовой смеси в камере, парциальное давление БЩ4, объемное соотношение GeH4/SiH4, скорость потока газов и ТП являются основополагающими факторами технологического процесса осаждения ультратонких наноструктурированных слоев ПКК и больше всего влияют на скорость осаждения. Дополнительное присутствие GeH4 в малых количествах (0,001-0,003 об. частей) в зоне разложения БЩ4 способствует увеличению V, в 1,5-2 раза.
Энергия активации процесса разложения БЩ4 без и в присутствии GeH4 в области температур 560-630°С, в зависимости от условий разложения МГ в кинетической области изменяется от 104 до 207, а в диффузионной области - от 11,6 до 72 кДж/моль.
Повышение общего давления в реакторе парциального давления SiH4 и ТП в процессе разложения SiH4 с поверхности подложки приводит к переходу в газовую фазу, т.е. гетероген-но-радикально-цепной маршрут, преобладающий при температуре 620±5°С, приобретает гете-рогенно-гомогенный цепной характер при температуре 630 °С.
CONDITIONS FOR THE FORMATION OF NANOCLUSTERS DURING THE PREPARATION OF NANOSTUCTURED FILMS OF POLYCRYSTALLINE SILICON DOPED WITH GERMANIUM
N.S. STROGOVA, AS. STROGOVA, A.A. KOVALEVSKII, V.M. BORISEVICH, D.V. PLYAKIN
Abstract
Investigations of surface topography, crystal structure of polycrystalline silicon films doped with germanium, obtained by decomposition of monosilane - monogermane - hydrogen mixture are presented. Changes in the structure of Ge and SiGe - clusters were analyzed by Raman spectroscopy. The dependencies of film growth on the total pressure in the reactor and the partial pressure of silane were also analyzed.
Литература
1. Ковалевский А.А., Борисевич В.М, Долбик А.В. // Материалы. Технологии. Инструменты. 2007. Т. 12, №4. С. 65-72.
2. Afanas'er V.V. // Microelectronic Engineering. 2001. Vol. 59. P. 335-339.
3. Ahmed W.D., Meakin B., Stoemenous J. // Proc. 10th Inc. Conf. Chem. Vapor Deposit. Pennington. 1987. P. 449-459.
4. Potapov A. V., Orlov L.K. // Proccedings of 6th Inc. Symp. Nanoctructures: Physics and Technology. St. Peterburg, Russia. 1998. P. 487-490.
5. Ковалевский А.А. // Материалы Междунар. науч.-техн. конф. Новые технологии изготовления многокристальных модулей. Минск. 2002. С. 110-114.
6. Ковалевский А.А., Крищенко В.П. // Материалы Междунар. науч.-техн. конф. Новые технологии изготовления многокристальных модулей. Минск. 2002. С. 44-48.
7. СухойМ.С. // Неорганические материалы. 1982. Т. 18, №7. С. 1077-1083.
8. Kovalevsky A.A., Borisenko V.E., Borisevich V.M. et al. // J. Inorganic Materials. 2005. Vol. 41, №12. P. 1260-1265.
9. Ковалевский А. А. // Материалы. Технологии. Инструменты. 2006. Т. 11, №1. С. 51-55.
10. Ковалевский А. А., Долбик А.В. // Материалы. Технологии. Инструменты. 2006. Т. 11, №2. С. 60-65.
11. Kovalevsky A A., Dolbik A V., Votech S.N. // Russian Microelectronics. 2007. Vol. 36, №3. P. 179-185.
12. Ковалевский А.А., Долбик А.В., Борисевич В.М. и др. // Материалы междунар. науч. конф. Нанострук-турные материалы. Минск. 2008. С. 138.
13. Ковалевский А.А., Долбик А.В., Борисевич В.М. и др. // Материалы 1-ой междунар. науч. конф Наност-руктурные материалы. Минск. 2008. С. 139.
14. HuangXinming et. al. // Jap. J. Appl. Phys. Pt.2. Letters. 2000. Vol. 39, №11B. P. L1115- L 1117.
15. Paul D.J. et al. // Advanced Materials. 1999. Vol. 11, №3. Р. 191-204.
16. Oliver Kermarrec et al. // Solid State Technology. 2002. Vol. 45, №3. P. 55.
17. Ковалевский А.А., Строгова А.С., Плякин Д.В. // Материалы. Технологии. Инструменты. 2009. Т. 14, №3. С. 90-95.
18. Eversteyn E.C. // Philips Res. Rep. 1974. Vol. 29. P. 45.
19. Liehr M. et al. // Appl. Phys. Lett. 1990. Vol. 56. P. 629.
20. Meakin D. // J. Appl. Physics. 1987. Vol. 67, №11. P. 5031.
21. Yu H.Y. et al. // Appl. Phys. Lett. 2002. Vol. 81. P. 376-378.