Неравновесные фазовые переходы индуцированные ионной бомбардировкой, на примере силикатных систем

Личный вклад автора:

- разработка, совместно с научным руководителем, теоретической модели, позволяющей описывать процессы, протекающие при термообработке систем, полученных ионной имплантацией;

- разработка, совместно с научным руководителем, алгоритма определения функции плотности свободной энергии Гиббса как функции температуры и объёмной концентрации тройных систем , ;

- разработка алгоритма численного решения квазилинейных уравнений параболического типа, возникающих в рассматриваемой задаче моделирования, обладающего безусловной устойчивостью по начальным данным и саморегуляризацией, связанной с переменным шагом по времени;

- разработка программного обеспечения для реализации вычислительных алгоритмов, связанных с применением разработанной теоретической модели;

- моделирование, с использованием разработанного подхода, процессов перераспределения компонентов, протекающих при термообработке систем , , , полученных ионной имплантацией;

- анализ результатов, полученных в ходе моделирования.

Публикации. По теме диссертации опубликовано 11 печатных работ, список которых приведён в конце автореферата.

Объем и структура диссертации. Диссертация состоит из введения, пяти глав, заключения и списка цитируемой литературы. Общий объем работы составляет 133 страницы, в том числе 118 страниц основного текста, 54 рисунка и список литературы из 114 наименований.

Краткое содержание работы

Во введении описывается тема работы, актуальность решаемой задачи, формулируются цели работы, даётся её структура и содержание.

В первой главе приводится аналитический обзор работ, посвящённых ионному синтезу скрытых диэлектрических слоёв. В данный момент хорошо известны технологии синтеза скрытых слоёв диоксида (SIMOX), нитрида (SIMNI) и оксинитрида кремния. Наибольшее практическое применение из них имеет SIMOX - процесс, который, как указывалось выше, не свободен от недостатков. Это обусловлено рядом его характерных особенностей.

Применяемые при ионном синтезе дозы составляют порядка (SIMOX, SIMNI) и несколько меньше при формировании оксинитрида, при энергии имплантации . С увеличением используемой в процессе энергии возрастает и величина дозы примеси. Снижение количества имплантируемой примеси - одна из задач, стоящих перед ионным синтезом. На данный момент разработаны варианты упомянутых технологий, позволяющие это делать. В частности, это имплантация в горячую подложку, имплантация с варьированием энергии ионов, синтез с помощью повторяющихся циклов имплантация - отжиг, стимулирование фазообразования предварительным созданием центров преципитации или дальнодействующим ионным облучением. Но окончательно данная задача не решена, поскольку дозы имплантации остаются достаточно высокими и, вместе с этим, качество формируемых скрытого и рабочего слоёв при уменьшении доз снижается.

Вторым вопросом, стоящим перед ионным синтезом, является снижение температуры стадии отжига. Для SIMOX - процесса она составляет . Термический отжиг при данной температуре проводится в течение нескольких часов и ведёт к значительным финансовым затратам и расходу материала рабочего слоя на термическое окисление. Снижение используемых доз реактивной примеси влечёт, как правило, увеличение указанной температуры стадии отжига. Таким образом, оптимизация параметров процесса ионного синтеза является важной задачей. При этом ясно, что наиболее эффективным способом оптимизации является теоретическое моделирование.

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

Альтернативой описанным методам может быть модель, в основе которой лежит классическая теория фазовых переходов Гинзбурга-Ландау [1]. Модель предполагает, что для любой системы можно найти “параметр порядка”, изменение которого отражает фазовый переход, претерпеваемый системой, и для которого справедливо релаксационное уравнение Халатникова - Ландау.

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

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

(1)

где - концентрация атомов -го компонента, находящихся в каком - либо из возможных для данного компонента состоянии, - время, - кинетическая постоянная, - функционал (2), аналогичный функционалу свободной энергии Гиббса.

(2)

где - коэффициенты диффузии, - объём, а функция на стадии постимплантационного отжига представляет собой функцию плотности свободной энергии.

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

- предполагается, что для каждого химического элемента существует полный дискретный набор состояний в данном материале;

- предполагается, что все атомы данного элемента, находящиеся в одном состоянии, считаются неразличимыми и характеризуются определенным значением концентрации в любой точке материала;

- предполагается, что если некоторый набор атомов входит в какое-то объединение, например, в состав какого-либо комплекса, то для этого комплекса существует ряд уравнений связи, определяемый составом комплекса и позволяющий уменьшить число переменных;

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

Суммируя уравнения системы (1) по всем состояниям данного компонента и вычисляя вариационные производные с учетом (2), систему уравнений (1) можно свести к системе, число уравнений в которой совпадает с количеством компонентов:

(3)

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

Основной задачей, которую необходимо решить при использовании уравнений вида (3), является задача определения вида функции плотности свободной энергии. В данной работе при решении этой задачи использовался феноменологический подход, основанный на анализе фазовых диаграмм. Преимущество данного подхода состоит в том, что он позволяет получать искомые выражения при минимальном количестве сведений о физической природе и составе моделируемой системы. При анализе твёрдых растворов можно ограничиться учётом лишь плотности свободной энергии смешения и рассматривать только конфигурационную часть энтропии [2]. Это приближение также использовалось в данной работе.

Решение системы уравнений в частных производных вида (3) с помощью аналитических и приближённо-аналитических методов представляет собой в общем случае сложную задачу. При произвольном выборе вида функции G более эффективным является применение численных методов. Исследование известных численных схем привело к необходимости разработки метода, обладающего лучшими характеристиками применительно к решаемой задаче. Разностная схема, с такими показателями, строилась для уравнения вида (3), с возмущающей функцией в виде степенного полинома. Используя разностную аппроксимацию для членов и заменяя частную производную на полную, а нелинейный член его приближением в виде разложения в ряд Тейлора , получили выражение для разностной схемы:

(4)

откуда в первом приближении:

(5)

В (4), (5) индекс обозначает сорт компонента; индексы и обозначают номер узла аппроксимирующей сетки, по координатным осям и соответственно; - параметр схемы, представляющий собой сумму множителей при ;- шаг по переменной ; - сумма слагаемых из (4), не содержащих в качестве множителя .

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

В третьей главе предлагаемый подход применяется к расчёту динамики системы .

Энтальпии смешения для бинарных оксидных систем были представлены в виде:

(6)

где x_i - мольная доля - го компонента, - зависящие от температуры коэффициенты, подлежащие определению.

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

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

(7)

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

С использованием полученных результатов, было выполнено моделирование динамики системы при различных исходных концентрациях кислорода и при различных температурах. Рассматривались концентрационные профили со стехиометрической и тремя значениями достехиометрической концентрации кислорода в максимуме. На рис. 2 показан результат моделирования с использованием достехиометрической концентрации кислорода при температуре .

Рис. 1. Функция свободной энергии смешения для системы Si-O

Рис. 2. Изменение концентрационных профилей кислорода, имплантированного в кремний, под влиянием термообработки. Температура . 1 - исходное распределение (концентрация в максимуме ), 2 - длительность термообработки , 3 - длительность термообработки

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

Четвёртая глава посвящена моделированию динамики системы в течение отжига.

В работе построение термодинамических функций тройных систем проводилось по следующему алгоритму. На первом этапе для получения функции свободной энергии смешения для тройной системы последовательно вычисляли все функции избыточной энтальпии вида (6), для всех бинарных систем, которые возможно выделить из состава рассматриваемой тройной системы. При определении функций для системы учитывалось, что в ней, при процессе ионного синтеза реализуемы два равновесных состояния. Первое из них соответствует кремнию, содержащему сравнительно низкую концентрацию бора (см^-3). Второе располагается в области, где концентрации бора и кремния одного порядка и соответствует соединению . Полученная функция была нормирована к стандартной энергии образования соединения , найденной по известным из литературы зависимостям [4],[5]. График полученной функции показан на рисунке 3. В качестве единственного равновесного состояния системы рассматривалось соединение . График полученной функции показан на рис.4.

Далее в работе был найден вид полинома, представляющего собой модельную функцию для избыточной энтальпии смешения для трёхкомпонентной системы. Коэффициенты полинома определялись по известным данным о величине энтальпии смешения для псевдобинарной системы [6] и значениям определённых функций бинарных систем. Общий вид избыточной энтальпии смешения:

 (8)

(9)

где - коэффициенты, функции температуры; - мольные доли.

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

Затем, исходя из соотношений линейной неравновесной термодинамики, были определены недостающие коэффициенты пропорциональности для системы уравнений вида (3).

Полученные результаты были применены при расчёте динамики системы . Во всех расчётах концентрация кислорода в максимуме имплантационного профиля была достехиометрическая. Исходные профили концентрации кислорода и бора моделировались с помощью программы TRIM. Энергия имплантации во всех случаях была для ионов кислорода и для ионов бора. Концентрационная зависимость коэффициента диффузии кислорода в системе аппроксимировалась линейной функцией по его значениям в и . Зависимостью от концентрации бора

Рис. 3. Функция свободной энергии Гиббса для системы при температуре

Рис. 4. Функция свободной энергии Гиббса для системы при температуре

пренебрегали, поскольку его массовое содержание в рассматриваемой системе недостаточно велико для оказания существенного влияния. При расчёте также учитывалась концентрационная зависимость коэффициента диффузии бора в системе от концентрации кислорода и бора . Его зависимость от концентрации кислорода представляет собой линейную аппроксимацию:

(10)

где коэффициенты и определялись по известным значениям коэффициентов диффузии бора в и .

Один из результатов расчётов можно видеть на рис. 5. На рисунке заметно стягивание профиля кислорода, что можно расценивать как формирование скрытого слоя новой фазы. Также видно формирование сложного профиля бора, с чётко выраженными экстремумами. При этом основной максимум локализуется в области высокого содержания кислорода.

Рис. 5. Изменение концентрационных профилей кислорода и бора, имплантированных в кремний, в результате термообработки. Температура отжига . Длительность отжига . 1 - исходное распределение кислорода (концентрация в максимуме ), 2 - исходное распределение бора (концентрация в максимуме ), 3 - профиль кислорода после термообработки, 4 - профиль бора после термообработки, 5 - область локального минимума на профиле бора

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

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

Результаты моделирования подтверждаются экспериментальными данными.

В пятой главе проведено моделирование динамики системы . Построение функции свободной энергии Гиббса для тройной системы осуществлялось по тому же алгоритму, что и для системы . Для построения искомого выражения были получены функции свободной энергии Гиббса бинарных систем и в дополнение к уже найденной для системы . При ионном синтезе для реализуемы два равновесных состояния. Первое из них соответствует кремнию, содержащему концентрацию фосфора на уровне _. Второе соответствует соединению . Стандартная энергия образования этого соединения, использованная для нормировки функции свободной энергии Гиббса, известна из литературы [7].

Функция свободной энергии Гиббса для системы рассчитывалась по аналогии с системой . Это является следствием того, что в обеих системах, при рассматриваемых температурах, не существуют области с низким содержанием кислорода. В качестве единственного равновесного состояния рассматривалось состояние, соответствующее соединению . Это сделано на основании того, что оксиды другого состава имеют значительно более низкую температуру плавления. Графики полученных функций , показаны на рисунках 6 и 7.

Коэффициенты полинома (9) были определены, опираясь на известные данные о величине энтальпии для псевдобинарной системы и

Рис. 6. Функция свободной энергии Гиббса для системы при температуре

Рис. 7. Функция свободной энергии Гиббса для системы . Температура

информацию об определённых ранее функциях для систем и . Коэффициенты пропорциональности определялись так же, как и для системы .

Полученные результаты были использованы при расчёте динамики системы . Профили концентрации кислорода и фосфора моделировались с помощью программы TRIM. Энергия имплантации во всех случаях была для ионов кислорода и для ионов фосфора. При расчётах для коэффициента диффузии кислорода учитывалась использованная ранее зависимость от концентрации кислорода в системе ; его зависимостью от концентрации фосфора пренебрегали. Коэффициент диффузии фосфора в системе рассматривался как функция концентрации кислорода и фосфора . Его зависимость от концентрации кислорода предполагалась линейной и описывалась выражением, аналогичным (9). Значения коэффициентов диффузии определялись по известным значениям коэффициентов диффузии фосфора в и .

Один из результатов проведённых расчётов можно видеть на рис. 8. Как можно заметить, изменения концентрационного профиля кислорода в результате термообработки имеют тот же характер, что и при формировании скрытого слоя в присутствии бора. На участке между скрытым слоем и поверхностным окислом концентрация кислорода заметно снижается по сравнению с исходным профилем. Профиль фосфора в результате непродолжительной термообработки приобретает сложную форму. Область наибольшей концентрации локализуется в формируемом скрытом слое. При этом форма распределения в этой области становится достаточно симметричной и подобной кислородному профилю за счёт смещения максимума распределения фосфора в направлении середины скрытого слоя. Кроме основного экстремума, на сформированном распределении фосфора наблюдается также два дополнительных максимума. Один из них, менее выраженный, расположен со стороны объёма подложки и существует лишь на начальной стадии отжига. В дальнейшем, в течение термообработки, наблюдается его деградация под действием диффузионных процессов. Второй локальный максимум локализован в растущем слое термического окисла.

Рис. 8. Изменение концентрационных профилей кислорода и фосфора, имплантированных в кремний, в результате термообработки. Температура отжига . Длительность отжига .1 - исходное распределение кислорода (концентрация в максимуме ), 2 - исходное распределение фосфора (концентрация в максимуме ), 3 - профиль кислорода после термообработки, 4 - профиль фосфора после термообработки

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

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

Результаты моделирования подтверждаются экспериментальными данными.

Основные результаты и выводы

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

2. Для системы установлено, что критическое значение объёмной концентрации кислорода, ниже которого, в данной области, формирование скрытого слоя не происходит, определяется максимумом функции плотности свободной энергии.

3. Установлено для систем и , что:

- скорость формирования скрытого слоя лимитируется скоростью диффузии кислорода в системе, полученной в результате имплантации;

- в результате термообработки происходит перераспределение кислорода и бора (кислорода и фосфора), в результате которого максимумы концентрационных профилей локализуются в одной области, а в области границ синтезируемого слоя наблюдается снижение концентрации бора (фосфора);

- при использовании в качестве дополнительно имплантируемого к кислороду компонента, как бора, так и фосфора, процесс ионного синтеза скрытого слоя происходит при более низких температурах по сравнению с классическим SIMOX - процессом.

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

5. Взаимодействие между атомами кислорода и бора (кислорода и фосфора) обусловливает локализацию максимумов концентрационных профилей в одной области.

Список цитируемой литературы

1. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Статистическая физика. - М.: Наука, 1976. Т. V., ч.1. - 584 с.

2. Суворов С.А., Сёмин Е.Г., Гусаров В.В. Фазовые диаграммы и термодинамика оксидных твёрдых растворов. -Л.:Из-во ЛГУ, 1986. - 140 с.

3. Ishimaru M., Tsunemori T., Harada S., Arita M., Motooka T. Microstructural evolution of oxygen implanted silicon during annealing processes // Nuclear Instruments and Methods in Physics Research B. - 1999. - №148. - P. 311-316.

4. Стромберг А.Г., Семченко Д.П. Физическая химия. - М.: Высшая школа, 1988. - 141 с.

5. Пригожин И., Кондепуди Д. Современная термодинамика. От тепловых двигателей до диссипативных структур - М.: Мир, 2002.- 461 с.

6. Baret G., Madar R., Bernard C. Silica-based Oxide Systems. Experimental and Calculated Phase Equilibria in Silicon, Boron, Phosphorus, Germanium, and Arsenic Oxide Mixtures // J. Electrochem. Soc. - 1991. - V.138. №9. - P. 2830-2835.

7. Равдель А.А., Пономаревская А.М. Краткий справочник физико-химических величин. -Л.: Химия, 1983. - 232 с.

Перечень работ, опубликованных по теме диссертации:

1. Кривелевич С.А., Денисенко Ю.И., Маковийчук М.И., Паршин Е.О., Цырулев А.А. Ионно - синтезированные силикатные слои - новое направление КНИ - технологии // Радиационная физика твёрдого тела: Тр. XIII международного совещания, Севастополь, 1 - 6 июля, 2002. - С. 315 - 319.

2. Кривелевич С.А., Цырулев А.А. Расчет функции свободной энергии и моделирование эволюции системы кремний - имплантированный кислород // Кремний - 2002: Тез. док. совещания по росту кристаллов, плёнок и дефектам структуры кремния, Новосибирск, 9 - 12 июля 2002 г. - Новосибирск, 2002.- С. 86.

3. Кривелевич С.А., Цырулев А.А. Моделирование ионного синтеза скрытых диэлектрических слоёв // Взаимодействие ионов с поверхностью: Cб. материалов XVI международной конференции, Звенигород, 25 - 29 августа 2003 г. - 2003. Т.2. - С. 94-97.

4. Krivelevich S., Buchin E, Denisenko Yu., Tsyrulev A. Ion synthesys of silicate glasses - simulation, process engineering and applied aspects // Micro - and nanoelectronics - 2003 (ICMNE-2003): Abs. Int. conf., Russia, Moscow - Zvenigorod, 6 - 10 October 2003.- P1 - 16.

5. Кривелевич С.А., Цырулев А.А. Моделирование ионного синтеза слоёв диоксида кремния в рамках теории Гинзбурга - Ландау // Радиационная физика твёрдого тела: Тр. XIII международного совещания, Севастополь, 30 июня - 5 июля 2003 г. - С. 393 - 397.

6. Бучин Э.Ю., Денисенко Ю.И., Кривелевич С.А., Паршин Е.О., Цырулев А.А. Диффузия бора и фосфора в слоях кремния, обогащённых имплантированных кислородом // “Кремний - 2004”: Тез. док. Иркутск, Листвянка, 5 - 9 июля 2004 г. - Иркутск, 2004. - С. 93.

7. Кривелевич С.А., Бачурин В.И., Цырулев А.А. Моделирование ионного синтеза скрытых диэлектрических слоёв // Поверхность. Рентгеновские, синхротронные и нейтронные исследования. - 2004.-№ 5.- С. 40-43.

8. Krivelevich S.A., Denisenko Yu. I., Tsyrulev A.A. Ionic synthesis of silica - based glasses: prospect, simulation and applied aspects // Proceedings of SPIE - May 2004. - V. 5401 - P. 119 - 128.

9. Кривелевич С.А., Цырулев А.А. Моделирование синтеза скрытых диэлектрических слоёв сложного состава, получаемых с помощью ионной имплантации // Вестник Нижегородского университета. Серия “Физика твёрдого тела”. - Вып. 1 (8). - 2005. - C. 91- 96.

10. Кривелевич С.А., Цырулев А.А. Моделирование формирования скрытых слоёв комбинированного состава // Тонкие плёнки в электронике: материалы XVII международного симпозиума, Москва, МГТУ им. Н.Э. Баумана, 2005, 8-10 сентября, C.393-397.

11. Krivelevich S.A., Tsyrulev A.A. Simulation of ion - beam synthesis of silica - based compounds// Micro - and nanoelectronics - 2005 (ICMNE-2005): Abs. Int. conf., Russia, Moscow - Zvenigorod, 3 - 7 October 2005.- P1 - 49.

Размещено на Allbest.ru