Расчет профиля диффузии сурьмы в кремнии

АННОТАЦИЯ



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

Полученные результаты хорошо согласуются с известными данными из литературы.

Полученные таким методом данные могут широко использоваться при расчете необходимой глубины залегания примеси.

Объем расчетно-пояснительной записки составляет ___ листа, она содержит 19 рисунков, 2 таблицы, приложение. В работе использовано 7 литературных источников.




СОДЕРЖАНИЕ



	ВВЕДЕНИЕ

		ТЕОРЕТИЧЕСКАЯ ЧАСТЬ

		Феноменологические модели

		Механизмы диффузии атомов в полупроводниках

		Коэффициент диффузии примесного атома

		Зарядовые состояния точеных дефектов

		ПРАКТИЧЕСКАЯ ЧАСТЬ

	ВЫВОДЫ

	СПИСОК ЛИТЕРАТУРЫ

	ПРИЛОЖЕНИЕ




ВВЕДЕНИЕ



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

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

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




1 ТЕОРЕТИЧЕСКАЯ ЧАСТЬ



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



1.1 Феноменологические модели



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



 (1)



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



 (2)




Закон Фика можно обосновать теоретически методами термодинамики. Именно, плотность потока  примеси А в общем случае даётся произведением её концентрации  на обобщённую скорость , которая, в свою очередь, есть произведение подвижности  на обобщённую термодинамическую силу ,



 (3)



Подвижность связана с коэффициентом диффузии соотношением Эйнштейна



 (4)



Термодинамическая сила равна градиенту химического потенциала примеси  с противоположным знаком,



 (5)



Сам же химический потенциал может быть записан в виде



 (6)



где  - равновесная растворимость примеси в матрице. В однородном полупроводнике , и из (4), (5) и (6) получаем



 (7)

Уравнение непрерывности следует из интегрального баланса диффундирующих частиц в замкнутом объёме , ограниченном поверхностью  (рис. 1).





Рис. 1. Баланс диффундирующих частиц



Полное число частиц в объёме  может измениться вследствие 4 причин:

- вытекания их из объёма  через поверхность . Скорость этого процесса равна 

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

- генерации частиц внутренними источниками с локальной скоростью , тогда во всём объёме за 1с возникнет  частиц.

- генерации частиц внешним источником с локальной скоростью , тогда в объём вводится за 1с  частиц.

Уравнение баланса имеет вид





Оно справедливо для любого объёма V. Преобразовав поверхностный интеграл в объёмный с помощью теоремы Остроградского-Гаусса, получим уравнение непрерывности в виде



  (8)



Или, после подстановки (1) в (8) ,



   (9)



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



 (10)



 (11)



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



 (12)




 (13)



с начальными условиями 



 (14)



 (15)



Таким образом, в этом случае следует решать систему уравнений (8), (12), (13). Конкретный вид слагаемых, описывающих изменения концентраций внутренних источников и стоков, определяется физическими механизмами взаимодействия их с диффундирующими частицами.

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





Уравнения (8) и (9) принимают вид



(16)



(17)



Если внутренние источники и стоки подвижны, то в уравнения (16) и (17) следует добавить слагаемые, описывающие их диффузию и дрейф,



 (18)

 (19)



Уравнения (9), (18) и (19) необходимо дополнить граничными условиями на свободной границе с вакуумом и границах раздела контактирующих сред различной физической природы.

Если через свободную поверхность  производится диффузия примеси в объём («загонка примеси»), то на ней задаётся постоянная концентрация,



 (20)



Если производится диффузия примеси из заданного начального распределения  («разгонка примеси»), то свободная поверхность считается непроницаемой для примеси, что эквивалентно равенству нулю производной от  по нормали к ,



 (21)



На границе раздела  сред 1 и 2 задаются условие непрерывности потока дифузанта



 (22)



и условие соответствия концентраций



, (23)


где  - коэффициент распределения, равный отношению предельных растворимостей диффузанта в средах 1 и 2.[1]



1.2 Механизмы диффузии атомов в полупроводниках



В кристаллических полупроводниках диффундируют собственные и примесные атомы. Диффузия собственных атомов называется самодиффузией, диффузия примесных атомов - примесной диффузией.

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

Примеси внедрения малого радиуса диффундируют по прямому междуузельному механизму. Это значит, что они совершают переход или акт миграции, или скачок, непосредственно из одного междуузлия в другое (рис. 2).





Рис. 2. Прямой междуузельный механизм диффузии



Коэффициент диффузии  по прямому междуузельному механизму,

 (24)



где , - энергия активации диффузии примеси внедрения.

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





Рис. 3. Вакансия  и собственный междуузельный атом 



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


(25)



При  , .

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

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

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



 (26)



Исчезновение комплекса происходит по реакции



 (27)



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





Рис. 4. Вакансионный механизм диффузии



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

Выделяют три разновидности междуузельного механизма диффузии:

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



 (28)



- механизм вытеснения описывается квазихимической реакцией



 (29)



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

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



 (30)



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



Рис. 5. Междуузельный механизм диффузии



Исчезновение примесно-междуузельного комплекса происходит по реакции



(31) [1]





1.3 Коэффициент диффузии примесного атома



Если бы все примесные атомы  входили в состав одного из комплексов, то коэффициент диффузии примеси  был бы равен коэффициенту диффузии соответствующего комплекса,

либо 

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



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





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



Тогда (32)



В состоянии термодинамического равновесия



(33)



Относительный вклад каждого механизма характеризуют параметрами



(34)



(35)



Они связаны естественным соотношением


 (36)



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



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







 (37)



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

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

2. В IV группе наблюдается переход от междуузельного механизма к вакансионному с возрастанием атомного радиуса:



3. Та же тенденция наблюдается в V группе:



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



 (38)



где  - обобщённый индекс дефекта.



1.4 Зарядовые состояния точеных дефектов



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

Рассмотрим сначала вакансии. Установлено, что вклад в диффузию примесей вносят в основном пять зарядовых состояний – нейтральноу , однократно и двукратно положительно заряженные  и , однократно и двукратно отрицательно заряженные  и . Состояние  реализуется при захвате нейтральной вакансией одного электрона, состояние  - при захвате ею двух электронов. Такие состояния являются акцепторными. Аналогично состояния  и  реализуются, если вакансия отдаёт один или два электрона соответственно. Такие состояния являются донорными. Нейтральная вакансия, следовательно, должна обладать двумя электронами, которые могут быть отданы. Таким образом, можно считать, что зарядовые состояния вакансии создаются четырьмя электронами.

Равновесные концентрации заряженных вакансий определяются положениями их энергетических уровней  и уровня Ферми .[1]

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

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

Можно получить



 (39)



 (40)



 (41)



 (42)



 (43)




 есть уровень Ферми в собственном полупроводнике. Все каналы дают аддитивный вклад в результирующий коэффициент диффузии ,









 (44)[1]




2 ПРАКТИЧЕСКАЯ ЧАСТЬ



Задание: Методом численного решения одномерного уравнения диффузии необходимо исследовать зависимость профиля концентрации Sb в кремний при диффузии из поверхностного источника постоянной концентрации N0 от парциальных вкладов различных зарядовых состояний точечных дефектов в коэффициент диффузии в интервалах температур и концентраций:

850 ˚С ≤ Т ≤ 1200 ˚С

1015 см-3 ≤ N0 ≤ 1021 см-3.

Для расчетов будем использовать программу Mathcad 11.

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

Диффузия сурьмы будет идти практически стопроцентно по вакансионному механизму [1], так как сурьма – элемент пятой группы и относительный вклад междоузельного механизма составляет порядка 1%. Коэффициент диффузии будем рассчитывать по уравнению (44).

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

Таким образом, выражение для коэффициента диффузии примет вид:



(45)



Согласно уравнению Аррениуса:



 (46)



, (47)



Т.о. (48)



Аррениусовские параметры приведены в табл. 1 [1]:



Таблица 1. Аррениусовские параметры Sb в Si[3]



Примесь





Sb

,

0,214



3,65

,

15,0



4,08



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



(49)



В данном уравнении вкладом слагаемого  можно пренебречь.

Введем обозначение .

Концентрации электронов и дырок можно выразить через концентрацию собственных носителей.



 (50)



 (51)



Концентрации заряженных дефектов также являются функциями температуры:



 (52)



 (53)



 (54)



 (55)



Аналогично находим:



(56)




(57)



.(58)



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



(59)







Совершенно аналогично для междоузельных атомов:





 (60)





Также находим :



 (61)



(62)




 (63)



 (64)



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



(65)



Из него мы определяем х.

Коэффициенты А, В, С и P имеют вид:



(66)



(67)



 (68)



, где(69)



, (70)




, (71)



, (72)



, (73)



, (74)



, (75)



, (76)



, (77)



, (78)

. (79)



В эти выражения входят равновесные концентрации вакансий и междоузлий, которые можно определить из (25).

Температурную зависимость для ширины запрещенной зоны определяет соотношение Варшни:



.(80)



Для кремния ,  эВ/К,  К [4]





см2/с

К-1

Рис. 6. Зависимость коэффициента диффузии от обратной температуры



Решая это уравнение и записывая коэффициент диффузии как

(81)



где EF – решение уравнения, получаем зависимости рис. 6.

Увеличивая масштаб, получим:





см2/с

К-1

Рис. 7. Зависимость коэффициента диффузии от обратной температуры в увеличенном масштабе





К-1

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

Зависимость в полулогарифмическом масштабе от обратной температуры приведена на рис. 8.

Увеличивая масштаб:





К-1

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





Рис. 10. Зависимость коэффициентов диффузии различных примесей в Si от обратной температуры[5]

Полученные зависимости подтверждают предполагаемую многими авторами зависимость D(T) = D0exp(-ΔE/kT).

Сравнивая данную зависимость с практическими результатами (рис. 10), например из [5], видим, что сходство теоретических и практических данных очень хорошее. Данные для сурьмы здесь приведены для максимальной поверхностной концентрации 1019 см-3 по работе [6].

Зависимость от концентрации в логарифмическом масштабе (рис. 11):





см2/с

см-3

Рис. 11. Зависимость коэффициента диффузии от концентрации доноров в логарифмическом масштабе



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

Для определения начала нелинейного участка, найдем производные коэффициента диффузии по концентрации (рис. 12).

Как видно из полученных зависимостей коэффициент диффузии не зависит от концентрации доноров до следующих значений концентрации (табл. 2).





см-3

см-3аб



см-3

см-3вг



см-3д

Рис. 12. Зависимость производной коэффициента диффузии от концентрации доноров при: а – 1123 К, б – 1223 К, в – 1323 К, г – 1423 К, д – 1473 К



Таблица 2. Пределы начала нелинейности в зависимости коэффициента диффузии

Т, К

1123

1223

1323

1423

1473

1573

Nd, см-3













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

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

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

Предельная растворимость сурьмы в кремнии составляет  см-3 при температуре 1573 К [7] и снижается до 2*1019 при 1123 К, таким образом можно утверждать, что коэффициент диффузии будет постоянен при любой концентрации доноров, то есть можно использовать решение уравнения Фика только для случая постоянного коэффициента диффузии.

Максимальная погрешность такого решения по коэффициенту диффузии составит (для 1473 К):



,(82)




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

Отметим сразу, что собственная концентрация при выбранном диапазоне температур превышает , в то время как максимальная предельная растворимость <  см-3, то есть концентрация сурьмы в кремнии никогда не превысит собственную концентрацию при равновесном процессе.

Пусть диффузия ведется на глубину 10 мкм. Зададим четыре температуры процесса (1173 К, 1273 К, 1373 К и 1473 К) и сравним полученные профили концентрации в этом случае на пятом временном слое. Выбранное t составляет полчаса.

Для этого решим уравнение диффузии:



(83)



где - безразмерные переменные.

Полный текст программы приведен в приложении 1.

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



см-3

Рис. 13. Профиль распределения примеси при разных температурах диффузии (время процесса – 2,5 часа, ND = 1019 см-3, температуры - 1173 К,  1273 К,  1373 К,  1473 К)



Как видно из полученных распределений (рис. 13) увеличение температуры значительно ускоряет диффузию, что проявляется в увеличении глубины залегания примеси. Из графика видно, что при температуре 1173 К концентрация примеси на глубине 0,1 мкм достигает значения 1010 см-3 за 2,5 часа, в то время как при 1473 К та же концентрация за то же время наблюдается на глубине 6 мкм. То есть целесообразно повышать температуру диффузии, это дает существенную экономию времени.

Посмотрим, как влияет на диффузию изменение поверхностной концентрации (рис.14). Пусть температура процесса Т=1473 К. Тогда возможные распределения концентрации (поверхностные концентрации 1018 см-3, 1019 см-3, 2*1019 см-3):



см-3

Рис. 14. Профиль распределения примеси при разных поверхностных концентрациях (время процесса – 0,5; 2,5; 5 часов, ND = 1018 см-3, ND = 1019 см-3, ND = 2*1019 см-3, температура 1473 К)



Обозначение на графике:

- - 0,5 часа, 1473 К, ND = 1018 см-3;

-- 2,5 часа, 1473 К, ND = 1018 см-3;

-- 5 часов, 1473 К, ND = 1018 см-3;

-- 0,5 часа, 1473 К, ND = 1019 см-3;

-- 2,5 часа, 1473 К, ND = 1019 см-3;

-- 0,5 часа, 1473 К, ND = 2*1019 см-3;

-- 2,5 часа, 1473 К, ND = 2*1019 см-3;

-- 5 часов, 1473 К, ND = 2*1019 см-3.

Очевидно, что чем больше поверхностная концентрация в источнике, тем большая концентрация доноров будет в распределении, и тем глубже будет лежать p-n-переход. Исходя из зависимостей, при поверхностной концентрации 2*1019 см-3 за 5 часов концентрация 1010 см-3 достигается на глубине 7,7 мкм, а за 2,5 часа на глубине 6,2 мкм.

Концентрация же в 1 см-3 при 5-часовом процессе и поверхностной концентрации 2*1019 см-3 достигается на глубине 13,3 мкм (рис. 15).





см-3

Рис. 15 Профиль распределения примеси при разных поверхностных концентрациях (время процесса – 0,5; 2,5; 5 часов, ND = 1018 см-3, ND = 1019 см-3, ND = 2*1019 см-3, температура 1473 К) в уменьшенном масштабе



Построим зависимость для температуры 1473 К и поверхностной концентрации 1019 см-3 и посмотрим, как примесь диффундирует во времени (рис. 16).

Из зависимости видно, что при данных условиях, для достижения1010 см-3 на требуемой глубине в 10 мкм требуется проводить процесс более 10 часов. За 10 часов на глубине в 10 мкм будем иметь концентрацию доноров в 2*106 см-3, а за полчаса – 10-3 см-3.



см-3

Рис. 16. Эволюция профиля распределения примеси (время процесса – 0,5;  1,5;  2,5;  7,5;  10 часов, ND = 1019 см-3, температура 1473 К)



Покажем максимальную погрешность при использовании постоянного коэффициента диффузии и учета его непостоянства (рис. 17):





см-3

Рис. 17. Погрешность определения профиля диффузии

Здесь  - концентрация 2*1019 см-3, 1473 К, время – 5 часов, постоянный коэффициент диффузии.  - то же, непостоянный коэффициент диффузии.

Видно, что при проведении процесса в течении 5 часов погрешность на почти 10 мкм при концентрации в 10-6 см-3 составит не более 60 Å, что допустимо даже для самого современного уровня производства. Это говорит о том, что для данной примеси можно вести расчет в предположении постоянства коэффициента диффузии.

Также интересно посмотреть на влияние на профиль вклада в коэффициент диффузии отрицательно заряженных дефектов. Для этого исключим из уравнения (45) второе слагаемое. Тогда получим (рис. 18):





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




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





Рис. 19. Сравнение профилей диффузии с учетом и без учета отрицательно заряженных дефектов (время 2,5 часа, температура 1473 К, ND = 2*1019 см-3)



Как видно, ошибка при таком вычислении при 1473 К, 2,5-часовом процессе и ND = 2*1019см-3 достигает почти 50%, то есть влиянием этого вклада пренебрегать никак нельзя.




ВЫВОДЫ



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

Проведено сравнение полученных температурных зависимостей коэффициента диффузии с практическими. Сходимость результатов найдена очень хорошей. Это говорит о довольно точном значении аррениусовских коэффициентов.

Определены пределы начала нелинейности коэффициента диффузии в зависимости от поверхностной концентрации. Было выяснено, что максимальная предельная растворимость сурьмы в кремнии ( см-3 при температуре 1573 К) не превышает этих значений. Определенная максимальная погрешность как в профиле (60 Å), так и в самом коэффициенте диффузии () при пренебрежении его непостоянством подтвердила эти выводы.

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

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

Определены точные значения времени процесса для проникновения примеси на необходимую толщину. Показано, что для ускорения процесса необходимо выбирать как можно большие значения поверхностной концентрации и температуры (2*1019 см-3 при 1473 К).

Также показано, что для сурьмы концентрация примеси ни при каких температурах не может превысит собственную концентрацию (например, , а максимальная предельная растворимость ).

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




СПИСОК ЛИТЕРАТУРЫ



1. Головатый Ю.П., Косушкин В.Г. Ионная имплантация и диффузия. Методические указания к выполнению курсовой работы по дисциплине «Математическое моделирование технологических процессов». – Калуга, 2007

2. Fahey P. M., Griffin P. B. and Plummer J. D. Rev. Mod. Phys.. - v. 61, 1989. - p.289

3. Van Vechten J.A. Handbook on Semiconductors, edited by T. S. Moss, vol. 3: Materials, Properties and Preparation, ed. S. P. Keller. - North-Holland, New York, 1980. - p. 1

4. http://ioffe.ru

5. Мазель Е.З., Пресс Ф.П. Планарная технология кремниевых приборов. – М.: Энергия, 1974. – 384 с.

6. Fuller, Ditzenberger // J. Appl. Phys, 1954, v. 25, p. 1439-1440; 1956, v. 27, p. 544-553

7. Бургер Р, Донован Р. Окисление, диффузия, эпитаксия. – М.: Мир, 1969