2007
Доклады БГУИР
октябрь-декабрь
№ 4 (20)
УДК 532.23: 533.9
ПЕРИОДИЧЕСКИЕ СТРУКТУРЫ, СФОРМИРОВАННЫЕ НА ПОВЕРХНОСТИ МОНОКРИСТАЛЛИЧЕСКОГО КРЕМНИЯ ПРИ ВОЗДЕЙСТВИИ КОМПРЕССИОННОЙ ПЛАЗМЫ
Н.Т. КВАСОВ1, Ю Г. ШЕДКО1, В В. УГЛОВ2, В.М. АСТАШИНСКИЙ3
'Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
2Белорусский государственный университет пр. Независимости 4, 220080, Минск, Беларусь
3Институт молекулярной и атомной физики НАН Беларуси пр. Независимости, 70, 220072, Минск, Беларусь
Поступила в редакцию 28 марта 2007
Рассмотрено образование периодических структур под воздействием компрессионного плазменного потока на поверхности монокристаллической кремниевой пластины как результат зарождения и развития гидродинамической неустойчивости типа неустойчивости Кельвина-Гельмгольца, вызванной течением плазмы вдоль поверхности пластины. Предложено описание кремниевого расплава с помощью уравнений магнитной гидродинамики для идеальной несжимаемой жидкости. Для анализа стабильности системы плазма-расплав использована линейная теория устойчивости, адаптированная для случая изменяющейся со временем глубины расплавленного слоя. Получено аналитическое выражение, которое связывает период модулированной поверхности пластины с величинами, характеризующими течение плазмы.
Ключевые слова: периодические структуры, кремний, плазма, компрессионный поток, неустойчивость Кельвина-Гельмгольца.
Введение
В настоящее время важнейшую роль в технике играют материалы с целенаправленно модифицированными поверхностными слоями, поскольку именно свойствами поверхностных слоев определяются эффективность работы и долговечность компонентов сложных компактных систем. Поэтому поиск новых материалов и разработка новых методов модификации свойств уже существующих материалов остаются неизменно актуальными. С данным направлением современного материаловедения тесно связана и одна из основных задач микро- и нанотехнологии — проблема формирования низкоразмерных структур на поверхности или в приповерхностных слоях твердых тел. Несмотря на внедрение новых (часто экзотических и дорогостоящих) материалов, чрезвычайно важно, что еще не исчерпаны возможности материалов, традиционно применяемых в микроэлектронике. Подтверждением этому служат эксперименты по воздействию компрессионных плазменных потоков (КПП) на монокристаллические кремниевые пластины (МКП).
Экспериментально показано [1-3], что при воздействии КПП на поверхности МКП образуются субмикронные структуры, состоящие из параллельных друг другу цилиндров
длиной 50-100 мкм, диаметром 0,1-0,8 мкм. Цилиндры располагаются с периодом 1-2 мкм на поверхности МКП. В процессе обработки МКП поглощает энергию 5-25 Дж/см2, проплавляется на глубину 6-10 мкм и затем вновь кристаллизуется. Поэтому толщина расплава возрастает от нуля (в момент начала плавления) до максимального значения, потом убывает обратно до нуля к моменту завершения кристаллизации. При проведении исследований КПП получали в газоразрядном плазменном ускорителе типа магнитоплазменный компрессор [4]. Исследования велись в режиме остаточного газа (откачанная вакуумная камера компрессора заполнялась азотом до давления 400 Па). Концентрация заряженных частиц плазмы в области максимального сжатия достигает (5-10)-1017 см-3, температура — 1-3 эВ. Скорость плазмы составляет (4-7)-106 см/с.
В [5] предпринята попытка объяснить формирование субмикронных структур на МКП как реализацию минимума поверхностной энергии в рамках предположения о ее сильной анизотропии. В [6] проведено численное моделирование процесса формирования цилиндрических структур на основе уравнений Эйлера. В [5] также рассмотрены термоупругие напряжения, возникающие в твердой части МКП. Однако эти попытки не смогли объяснить наблюдаемую экспериментально морфологию поверхности. Не была установлена связь периода модулированной (возмущенной) поверхности с параметрами плазмы, не учитывались зависимость толщины расплава от времени и течение плазмы вдоль МКП и расплава.
Течение плазмы вдоль кремниевой пластины и неустойчивость Кельвина-Гельмгольца
Неустойчивость Кельвина-Гельмгольца (НКГ) возникает в случае плоскопараллельного течения сред, отделенных плоской границей раздела [7, 8]. При этом граница претерпевает модуляцию вследствие резкого роста первоначально малых возмущений. Возмущение, имеющее максимальную скорость роста, называют модой максимальной нестабильности (ММН). Длина волны ММН имеет порядок а/К0, где а«1 Дж/м2 — поверхностное натяжение жидкого кремния [9], К0~106-107 Дж/м3 — плотность кинетической энергии плазмы [1-3]. Легко видеть, что длина волны ММН в нашем случае составляет 10-10 мкм, поэтому необходимо учесть возможность возникновения НКГ и исследовать стабильность системы расплав-плазма по отношению к малым возмущениям, длины волн которых лежат в данном диапазоне.
Кроме течения плазмы вдоль МКП необходимо учесть и другие факторы. Известно [10], что электропроводность жидкого кремния имеет порядок 106 Ом-1-м-1, что является характерным для металлов (например, ртути). Согласно [10], другие свойства жидкого кремния также отвечают свойствам типичных металлов. Это касается, в частности, кинетических коэффициентов. Описание течения кремниевого расплава будем вести, как это часто делают при описании течения жидких металлов, с помощью уравнений магнитной гидродинамики (МГД). С целью упрощения используем уравнения МГД для случая идеальной несжимаемой среды [11]:
р| ~ + (-V) | = -У(р + В ■ Н/2) + (Н -У)В,
д?
(1)
&у V = 0,
(2)
дВ
= 10
(г?х В),
(3)
В = 0 .
(4)
Здесь V — скорость течения среды; р — давление; р — плотность; В = ц0Н — индукция магнитного поля; Н — напряженность магнитного поля; ц0 — магнитная проницаемость вакуума. Силы тяжести не учитываются, так как на микронных масштабах их
влияние незначительно на фоне других факторов. К уравнениям (1)-(4) нужно еще добавить уравнение состояния и уравнение сохранения энтропии [7, 11].
Известно, что одной из отличительных особенностей КПП является мощное вмороженное магнитное поле (см. [1-3] и ссылки там). Чтобы учесть взаимодействие расплава с полем плазмы и вместе с тем не усложнять чрезмерно задачу, мы используем (1)-(4) совместно с уравнением сохранения энтропии и соответствующим уравнением состояния для описания динамики плазмы. В дальнейшем под плазмой мы будем понимать не КПП в целом, а плазму, находящуюся у поверхности расплава.
К сожалению, найти аналитическое решение (1)-(4) в общем виде не удается, однако в этом нет необходимости при исследовании на устойчивость. Мы имеем систему из двух текучих сред, граничащих друг с другом и разделенных первоначально плоской границей раздела. Введем декартову систему координат так, чтобы плоская граница раздела совпадала с плоскостью ХОУ. Рассматривая НКГ, положим, что среды могут иметь отличную от нуля скорость течения V(г), параллельную плоскости ХОУ. Плотность р(г), давление р0(г),
скорость V(г), магнитное поле В(г) и удельную энтропию также будем считать зависящими только от координаты г. Магнитное поле будем считать параллельным плоскости ХОУ. Плотность, давление и удельная энтропия связаны между собой в каждой точке уравнением состояния. То, что все указанные величины считаются функциями только одной координаты, не будет являться существенным ограничением общности, если они мало меняются на расстояниях порядка длины волны ММН (т.е. порядка 1 мкм) вдоль плоскости ХОУ .
Для выяснения характера устойчивости начального стационарного состояния рассмотрим малое отклонение от него, которое численно выразим в виде малых возмущений в распределениях всех величин, характеризующих систему. Так, давление теперь будет р0(г) + р , скорость V0(г) + и , магнитное поле В (г) + Ь , где р , и, Ь — малые по величине произвольные возмущения. С учетом движения жидкостей новое положение займет и граница раздела. Будем считать, что форма поверхности теперь задается функцией г = х, у, ,). Так как начальное положение границы раздела задается уравнением г = 0, то х, у, ,) представляет собой возмущение и считается малым. Уравнениям МГД теперь будут удовлетворять уже возмущенные значения физических величин. Сохранив в (1)-(4) только члены не выше первого порядка по возмущениям, можем записать
р(г) (д, + V(г) • V)и = -V (р + Н(г) • Ь ) + (Н(г) • V) Ь + Ь2Й'(г) - п2р(г)V'(г) + , (5)
и = 0, (6)
(д, + V(г) • V)Ь = (((г) • V) + Ь Г(г) -и2 В'(г), (7)
^у Ь = 0, (8)
д
где д = —; е — орт вдоль оси 02, а штрих над векторными величинами означает
' д,
дифференцирование всех их компонент по координате г .
В (5) добавочный член 5 = а5(гучитывает поверхностное натяжение на границе раздела сред. Ось 02 направлена из расплава в плазму. Плазме отвечает область г > 0, расплаву -к(,) < г < 0 . к(,) > 0 — глубина расплава.
Обычно при исследовании гидродинамической стабильности различных физических систем применяют следующий прием [7, 8]. Проводят фурье-преобразования по х, у, , всех возмущений в уравнениях (5)-(8) и граничных условиях, отражающих в нашем случае невозможность смешения плазмы с расплавом и отсутствие протекания расплава в твердую часть пластины [7]. Однако из-за того, что к зависит от ,, нам следует проводить преобразования только по х и у. Действительно, проведение преобразования Фурье
по времени эквивалентно представлению решения в виде суперпозиции членов, содержащих временную зависимость в виде exp(/p t). Подобную временную зависимость можно ожидать, когда толщина расплава неизменна. В случае изменяющейся со временем толщины расплава преобразование Фурье по времени не эффективно, т. к. не приводит к упрощению задачи. В дальнейшем будем обозначать преобразованные величины теми же самыми символами, что не вызовет недоразумений. Отметим также, что в результате проведения фурье-преобразований по двум координатам появится новый параметр к — волновой вектор, лежащий в плоскости XOY.
Далее необходимо различать два случая. В первом случае к • B(z) = 0, и путем несложных, но несколько громоздких преобразований удается исключить из (5) все неизвестные кроме uz, bz, £(t). Интегрируя по z, чтобы избавиться от дельта-функции, получаем еще одно граничное условие [8]. Добавляя к ним соотношения, необходимые для определения bz, £(t), получаем
к2р(z)Dzuz = 5z (ik • H'(z) bz - ik • V\z)p(z) uz + p(z)Dz3zuz ), (9)
Dzbz = 0, (10)
s(ik • H'(z) bz - ik • V'(z)p(z) uz +p(z)Dzd uz ) = a k4^(t), (11)
±0 Ç(t ) = u z |z=o+0, (12)
- - - д где введены обозначения Dz =dt +ik •V(z), д z =— и 8(y) = y(z + 0)-y(z - 0).
z dz
Во втором случае к • B(z) ф 0 . Аналогичным образом приходим к уравнениям к2L zy = d L zdz y, (13)
s(( z d z |) = ik 4a^(t), (14)
D0±0 () + />|z=0±0)= 0. (15)
где L z=p(z)Dz2 + (( •H(z))(( • B(z)) и посредством соотношения Dz| = iuz введена
вспомогательная величина |.
Соотношения (9)-(12) или (13)-(15) весьма громоздки и содержат частные производные по времени и координате z. Нужно помнить, что решаются эти уравнения для двух областей пространства с изменяющейся геометрией, и различать случаи к • B(z) = 0 и к • B(z) ф 0, что неудобно для расчетов. Рассмотрим частный случай, когда задача заметно упрощается. Предположим, что H, B, p, V постоянны в пределах каждой из областей, где z > 0 или -h(t) < z < 0 . Тогда можно показать, что и при к • B(z) = 0 и при к • B(z) ф 0 выполняется
(a к3 + p+ (+ + ^2 ) + p_ ( Д2_ + та - )(( coth (h )Dz- ))(t ) = J (t), (16)
где та± =^(к • H±)(к • B±))p± , Dz± = d + ih • V± , а D-1 означает оператор, обратный оператору
D . Знаками "+" и "-" отмечены величины, относящиеся к плазме и расплаву соответственно. Неоднородная часть J (t) зависит сложным образом от времени, от вида функции h(t) и функций, определяющих начальные условия для | или uz. Кроме того, Dz-Ç = 0 в те моменты
времени, когда к = 0 (т. е. в момент начала плавления и в момент завершения кристаллизации). Функция 3 (,) обращается тождественно в нуль, например, в тех случаях, когда течение среды потенциально. Ограничимся для упрощения случаем потенциального течения, 3(,) = 0 .
Несмотря на это, (16) в общем виде решить не удается, однако можно легко получить информацию о его асимптотических свойствах. При 3(,) = 0 в (16) явная зависимость от времени будет лишь через функцию к(,) , стоящую в аргументе гиперболического котангенса. Учитывая свойства котангенса, можно заключить, что когда к(,) близка к нулю, соЛ (кк) « кпри этом для решения (16) будет важен конкретный вид зависимости к от ,. Наоборот, когда кк >> 1, будет соШ (кк)«1. При этом детали поведения функции к(,) уже не сыграют заметной роли. На практике можно с хорошей точностью взять соШ(кк) «1 уже при кк > 1. Для таких волн выполняется уравнение, которое было бы справедливо для расплава неограниченной толщины:
(ак3 + р+ ((+ +ш+) + Р- ((- + ®-))(0 = 0 . (17)
Решением этого ОДУ является, как и должно быть для неограниченно толстого слоя расплава, ) ~ ехр(/С). Тогда из (17) следует
ак3 +((-Н+)((-В1+)-р+(со + к?-V;)2 + (к•#)) •В-)-р-(о + к •К)2 = 0. (18)
Из (18) можно выразить две ветви частот о (к), которые будут либо вещественными одновременно, либо комплексно-сопряженными с ненулевой мнимой частью. В последнем случае для одной из ветвей Яе(/га) будет положительной, и £(,) ~ ехр (/о,) будет
экспоненциально расти со временем. Такое поведение означает нестабильность системы плазма-расплав.
Чтобы проследить особенности расчетов, связанные с данными дисперсионными соотношениями о (к), рассмотрим частный случай. Пусть расплав и плазма текут в одном направлении (можно считать, что расплав увлекается потоком плазмы), выберем это направление в качестве оси ОХ. Пусть также В- = В+ = В (т. е. нет токов, текущих по границе плазма-расплав). Тогда
о = -ки0 (19)
Здесь введены следующие обозначения: и, =(рЛ+р_К)/(р,+р_), (20)
К = р,р. (П-^ , (2,)
р+ + р- 2
Ж(кх,ку) = ак3 + 2|-1 (к • В)2 - 2к_2К0. (22)
Величины К0 и и0 характеризуют, очевидно, движение сред относительно друг друга. Обозначение Ж вводится для упрощения дальнейшего анализа. Учитывая вид (19), легко заметить, что мнимая часть частоты появится тогда и только, когда будет отрицательной величина Ж, определяемая соотношением (22). Анализируя свойства функции Ж(кх, ку),
приходим к следующим выводам.
1. Стабилизировать систему, обеспечив неотрицательность Ж для всех к , может только магнитное поле величины В > К0 = В0, направленное вдоль направления течения
сред. Легко оценить, что B0 « 3 Тл. Таким образом, для стабилизации требуются очень большие значения поля.
2. Существует ММН. Длина волны моды максимальной нестабильности (а значит, и период модулированной поверхности) определяется как
= 3™- 1 , (23)
K0 1 -Г + ^/(Г-1)2 + 4Гsin2 ф
где Г = Б2/ц0K0 — безразмерный параметр; ф — угол между векторами B и V±. В тех
случаях, когда магнитное поле слишком слабое, чтобы стабилизировать систему, т.е. при Б < Б0, будет выполнено неравенство Г < 1 и Л0 будет вещественной неотрицательной величиной.
3. Важной характеристикой является также угол в между скоростью течения сред и
волновым вектором ММН, т. е. угол между OX и направлением, вдоль которого возникнет периодическая модуляция поверхности. Он определяется соотношением
18(29.) = (24)
cos 2ф-Г
Учитывая, что sin2 ф изменяется от 0 до 1, легко показать, что 1 -Г+^(Г-1)2 +4rsin2 ф
пробегает значения от 0 до 2. Максимальное значение этой функции дает, очевидно, наименьшую длину волны возмущения, которое может приводить к неустойчивости системы плазма-расплав при заданной величине магнитного поля и скорости течения. Максимальное значение достигается при ф = + п/2, т.е. когда поле ортогонально скорости течения сред. При этом зависимость от величины поля исчезает, а 9c = 0, так что периодичность проявится вдоль
направления течения сред. Тогда минимально возможная длина волны определяется соотношением
Л = — . (25)
2K0
Используя известные значения величин, которые входят в (25), и пренебрегая течением расплава, получаем, что Л лежит в интервале 0,2-2,5 мкм. Можно сразу отметить, что данный диапазон охватывает все наблюдающиеся экспериментально периоды рассматриваемых субмикронных структур [1-3]. Важно подчеркнуть, что диапазон значений Л, а не строго одно значение, как это следует из формулы (25), возникает из-за того, что входящие в (25) параметры известны с некоторой погрешностью.
Заключение
Исследована проблема образования периодических структур на поверхности пластины монокристаллического кремния при воздействии компрессионной плазмы как результат зарождения и развития неустойчивости. Предложена модель образования периодических структур, учитывающая взаимодействие расплава кремния с магнитным полем компрессионной плазмы, растекание плазмы вдоль поверхности кремниевой пластины и высокую электропроводность жидкого кремния. В рамках предложенной модели показано, что расплав кремния с плоской свободной границей неустойчив по отношению к сколь угодно малым возмущениям. Расчет показал, что модам максимальной нестабильности отвечает модуляция свободной поверхности с периодом, лежащим в диапазоне 0,2-2,5 мкм. Данный диапазон охватывает все наблюдающиеся экспериментально периоды субмикронных структур [1-3].
PERIODIC STRUCTURES FORMED ON THE SURFACE OF SILICON MONOCRYSTAL TREATED WITH COMPRESSION PLASMA FLOW
N.T. KVASOV, Y.G. SHEDKO, V.V. UGLOV, V.M. ASTASHYNSKI
Abstract
Formation of periodic submicron structures on the surface of silicon wafer under the action of compression plasma flow is considered to be a manifestation of Kelvin-Helmholtz hydrodynamic instability caused by the plasma flows along the melt free surface. The melt flow is described with the help of motion equations for an ideal non-compressible conducting medium. Linear stability theory modified to include the time-dependent melt depth is employed to analyze stability of plasma-melt system. The expression, which relates the modulated surface period and plasma parameters, is obtained.
Литература
1. Углов В.В., Анищик В.М., Асташинский В.В. и др. // Письма в ЖЭТФ. 2001. Т. 74, № 4. C. 234-236.
2. Uglov V. V., Anishchik V.M., Astashynski V. V., et al. // Surface and Coatings Technology. 2002. Vol. 158-159. P. 273-276.
3. Astashynski V.M., Ananin S.I., Askerko V.V. et al. // Surface and Coatings Technology. 2004. Vol. 180-181. P. 392-395.
4. Асташинский В.М., Ефремов В.В., Костюкевич ЕА. и др. // Физика плазмы. 1991. Т. 17, № 9. С. 1111-1115.
5. АнищикВ.М., УгловВ.В., Пунько А.В., и др. // Перспективные материалы. 2003. № 5. С. 5-11.
6. Anishchik V.M., Uglov V.V., Punko A.V. et al. // XXXIII Междунар. конф. по физике взаимодействия заряженных частиц c кристаллами: Тез. докл. Moscow. 2003. P. 116.
7. Ландау Л.Д., Лифшиц ЕМ. Теоретическая физика. Т. 6. Гидродинамика. М., 1986.
8. Chandrasekhar S. Hydrodynamic and Hydromagnetic Stability. Oxford: University Press, 1961.
9. Ниженко В.И., Флока Л.И. Поверхностное натяжение жидких металлов и сплавов (одно- и дЕухкомпонентные системы). Справочник. М., 1981.
10. Глазов В.М., Чижевская С.Н., Глаголева Н.Н. Жидкие полупроводники. М., 1967.
11. ЛандауЛ.Д., Лифшиц Е.М. Теоретическая физика. Т. 8. Электродинамика сплошных сред. М., 2003.