Целью работы является изучение динамики развития рельефа, теплового свечения и отражательной способности поверхности кремния под действием лазерного излучения вблизи порога плавления; на этой основе сделан вывод о реальном механизме образования рельефа. Схема экспериментальной установки приведена на рис. 1. В опытах использовались механически полированные поверхности <111> кремния толщиной 4 мм. Разогрев образца осуществлялся p-поляризованным излучением лазера на неодимовом стекле (lambda=1.06 мкм) с длительностью импульса излучения на полувысоте 0.8 мс. Интенсивность теплового излучения нагретой поверхности регистрировалась в спектральном интервале 0.4-0.9 мкм. Механическая модуляция излучения зондирующего He-Ne лазера (f=~ 60 кГц) использовалась для разделения сигналов дифракции и теплового свечения. Динамика сигнала теплового свечения поверхности для различных плотностей мощности q воздействующего излучения приведена на рис. 2. При q=< 114 кВт/см² наблюдается приблизительно колоколообразная временная зависимость (рис. 2, а), значение максимума которой плавно растет с ростом плотности мощности излучения. При q<= 120 кВт/см² вблизи вершины кривой появляется точка перегиба (рис. 2, б), с ростом q переходящая в минимум (рис. 2, в), глубина которого, а также расстояние между максимумами растут с q. Динамика сигнала дифракции зондирующего излучения, типичная для структур с периодом d=2lambda (0=< theta=< 9^o), приведена на рис. 3, а. Рис. 3, б иллюстрирует динамику структур lambda, ранее не наблюдавшихся в остаточном рельефе [1] (остающемся после затвердевания поверхности). Обнаружено наличие остаточных структур lambda, занимающих лишь 5-10% зоны воздействия.Здесь и ниже для сокращения текста наряду с полным терминов "структуры с периодом d" используется сокращенный вариант "структуры d". Динамика структур lambda и 2lambda регистрировалась одновременно в одной фазе воздействия, так же как и структур d^- и 2d^-. На рис. 3, в, г приведена динамика структур d^- и 2d^- (9=< theta=< 30^o) [2] соответственно, полученная при q=123 кВт/см² и lambda=12^o. Как показали опыты, существование структур 2d^- связано с существованием структур d^-: при формировании последних формировались и структуры 2d^-. При низких q рельеф 2d^- занимал преимущественно центральную часть зоны облучения, а d^- - периферийную. В интервале 95-170 кВт/см² глубина рельефа 2d^- и занимаемая им площадь возрастали, причем структуры d^- и 2d^- возникали в перемежающихся пространственно ограниченных областях. Наблюдавшееся при q<= 170-190 кВт/см² уменьшение амплитуды сигнала дифракции, соответствующего остаточным структурам 2d^-, связано с уменьшением занимаемой ими площади - формированием в пределах кольца, внутри которого зона локально гладкого расплава. Динамика изменения отражательной способности поверхности на lambda=0.63 мкм изучалась путем регистрации интенсивности зеркально отраженного излучения, что позволило определить время начала плавления и оценить снизу долю образующегося на поверхности расплава (заметим, что часть отраженного излучения терялась вследствие дифракции). Изменение отражения обусловлено как температурной зависимостью оптических постоянных кремния [3], так и появлением расплава. Аппроксимируя данные [3], получим, что при T=T_melt рост коэффициента отражения при lambda=0.63 мкм составляет =~25%. Оценка показывает, что для рис. 4, б максимальная доля расплава поверхности не превышает 20%. Моменты начала процесса плавления, определенные по точке перегиба кривой теплового излучения (рис. 4, в) и по кривой отражения (рис. 4, б), совпадают при учете температурной зависимости оптических постоянных. [!tb] Схема экспериментальной установки. 1 - юстировочный He-Ne лазер; 2 - импульсный лазер на неодимовом стекле; 3 - поляризатор; 4 - зондирующий He-Ne лазер; S - образец кремния; D - диафрагма, проектируемая линзой L на образец S; C - измеритель энергии импульса ИМО-2Н; ФД₁- ФД₃ - кремниевые фотодиоды ФД-256; Ф₁- Ф₄ - фильтры; Ф₅ - фильтр ИКС-1. Глубина остаточного рельефа обычно составляла 0.1-0.2 мкм. Профиль рельефа в поперечном направлении имел относительно плоское "дно" и резко поднимающиеся края (рис. 5). Обсудим данные по динамике излучательной способности поверхности. Отметим, что излучательная способность твердофазного кремния выше, чем жидкофазного. Действительно, спектральный коэффициент излучения varepsilon(lambda) определяется уравнением баланса varepsilon(lambda)=1-rho(lambda)-tau(lambda), где rho(lambda) и tau(lambda) - спектральные коэффициенты отражения и пропускания соответственно. Для расплава кремния (металла) varepsilon₁(lambda)=~0.2, так как rho(lambda)=~0.8 и tau(lambda)=~0. В области фундаментального поглощения varepsilon₂(lambda)=~0.7, так как rho(lambda)=~0.3, tau(lambda)=0 (толстая пластина) и varepsilon₂/varepsilon₁=3.5. [!b] Форма импульса теплового свечения поверхности кремния, возникающего под действием импульса излучения неодимового лазера. q. кВт/см²: а - 113, б - 134, в - 142, г - 215 кВт/см². Появление излома и минимума на кривой теплового свечения есть следствие неоднородного плавления поверхности [4]. Действительно, с ростом q увеличивается доля расплава на поверхности, что и вызывает спад интенсивностив свечения. Увеличение интенсивности свечения при остывании (второй максимум на рис. 2, в) связано с кристаллизацией расплава, поскольку при этом растет varepsilon. Сравнение динамики сигналов дифракции, отражения и теплового свечения показывает, что поверхность, соответствующая первому максимуму, представляет собой сосуществование областей жидкой и твердой фаз. Значения интенсивности кривой теплового свечения содержат информацию о доле расплава на поверхности. В предположении, что температура поверхности не превышает T_melt, можно оценить долю площади расплава на поверхности (S_x/S) S_x/S=(1-varepsilon₂/varepsilon₁)^-1Delta /A, (1) где A - амплитуда второго максимума; Delta=A-B, где B - величина сигнала в момент существования частичного расплава, например минимума кривой свечения. Оценка по формуле (1) величины S_x/S при взятом из эксперимента значении величины Delta, соответствующей началу роста резонансных структур при низких q, дает значение S_x/S=~0.26-0.3. Обсудим механизм образования структур 2lambda и lambda/(1-sintheta). В динамике можно выделить две стадии развития рельефа. Первая характеризуется началом локального плавления, фиксировавшегося по трем признакам: моменту появления излома в кривой теплового свечения, увеличению отражательной способности поверхности и наличию остаточного микрорельефа. Образование рельефа на этой стадии обусловлено как увеличением плотности кремния при плавлении, так и действием капиллярных сил в зонах локального расплава. По мере роста q размеры областей локального расплава растут, рельеф сглаживается и появляется минимум в динамике сигнала дифракции [!tb] Изменение интенсивности дифракции излучения зондирующего He-Ne лазера на различных структурах. а - 2lambda, б - lambda (theta=~ 0^o), в - d^-, г - 2d^- (theta=12^o); излучение неодимового лазера с q=140 и 120 кВт/см² соответственно. Штриховая линия - сигнал теплового свечения поверхности. Вторая стадия характеризуется резким ростом высоты рельефа при кристаллизации расплава в условиях его продолжающегося облучения. На этой стадии (рис. 3, а) за времена =<50 мкм происходит значительное перераспределение материала при совместном действии капиллярных сил и фактора разницы в плотностях расплава и твердой фазы. Это и обусловливает наблюдавшееся образование резко поднимающихся гребней в профиле рельефа. Впервые влияние последнего фактора на рост рельефа из расплава было учтено в [5]. Судя по динамике интенсивности свечения, структуры 2d^- образуются лишь на стадии кристаллизации жидкой фазы, поэтому механизм их образования таков же, как и для структур 2lambda и d^-. Существенное различие сигналов дифракции на структурах lambda и 2lambda и преобладание в остаточном рельефе последних свидетельствуют о динамической конкуренции роста структур двух типов, пространственно сосуществующих в зоне воздействия, возможная причина которой обсуждалась в [1]. [!b] Изменение коэффициента отражения ( а) и теплового свечения поверхности ( б) под действием импульса лазерного излучения ( в) при q=130 кВт/см² во взаимной временной привязки. [!b] Профиль рельефа поверхности кремния, полученный при использовании интерференционного микроскопа (q=140 кВт/см², структура 2lambda). Из экспериментальных данных по динамике развития рельефа, теплового сигнала и изменения отражения поверхности непосредственно следует, что формирование упорядоченного микрорельефа (lambda, 2lambda, d^-, 2d^-) происходит на частично расплавленной поверхности. Поляритонная модель разрушения поверхности (интерференция поверхностных электромагнитных волн (ПЭВ) с падающей волной ) [6] предполагает наличие поверхностно-активной среды, в данном случае сплошного расплава полупроводника. В настоящее время теория существования ПЭВ на частично металлизированной поверхности отсутствует. Ориентация и период структур lambda/(1-sintheta) находятся в рамках поляритонной модели, а для структур lambda, 2lambda, 2lambda/(1-sintheta) выходят за них. Таким образом, в работе изучена динамика теплового свечения поверхности и рельефа, формирующегося под действием лазерного излучения в режиме неоднородного плавления поверхности. Интерпретирована динамика поведения сигнала теплового свечения поверхности. Показано, что образование структур 2lambda, lambda/(1-sintheta), 2lambda/(1-sintheta) происходит в основном на стадии кристаллизации расплава и связано с действием капиллярных сил и фактора изменения плотности материала при фазовом переходе "плавление". Обнаружено образование преимущественно динамических структур lambda. Сделан вывод о том, что образование периодического микрорельефа происходит на частично расплавленной поверхности. Детальное объяснение причин формирования рельефа 2lambda и 3lambda и особенностей динамики рельефа lambda требует более подробного изучения. Данная работа частично поддержана грантами N RNS300 и R3W300 Международного научного фонда и правительства России и грантом N 94-02-05802 Российского фонда фундаментальных исследований.

Подсчитывается количество просмотров абстрактов ("html" на диаграммах) и полных версий статей ("pdf"). Просмотры с одинаковых IP-адресов засчитываются, если происходят с интервалом не менее 2-х часов.

Дата начала обработки статистических данных - 27 января 2016 г.