- •1. Моделирование процессов диффузии
- •1.1. Законы процесса диффузии
- •1.2.1. Модель диффузии из неограниченного источника
- •1.2.2. Модель диффузии из ограниченного источника
- •1.2.3. Примеры расчета этапов диффузии
- •1.3. Модель диффузии примесей из легированных оксидов
- •1.4. Модель диффузии примесей в кремнии с одновременным термическим окислением
- •1.5. Диффузия при высокой концентрации примесей
- •2. Моделирование процессов ионной имплантации
- •2.1. Модель Линдхарда, Шиотта и Шарфа
- •2.2. Модель Пирсона
- •3. Расчет сопротивления слоя
- •4. Моделирование структуры биполярного транзистора
- •4.1. Конструктивно-технологические особенности биполярных транзисторов интегральных микросхем
- •4.2. Расчет параметров транзисторной структуры при заданных режимах технологического процесса
- •4.3. Выбор режима технологического процесса для заданных параметров слоев структуры
- •Список РЕКОМЕНДУЕМОЙ литературы
- •ОГЛАВЛЕние
МИНОБРНАУКИ РОССИИ
––––––––––––––––––––––––––––––––––––––––––––––––––––
Санкт-Петербургский государственный электротехнический университет «ЛЭТИ» им. В. И. Ульянова (Ленина)
–––––––––––––––––––––––––––––––––––––––––––
И. И. ЗЯТЬКОВ А. Н. КРИВОШЕЕВА
БАЗОВЫЕ ПРОЦЕССЫ ПЛАНАРНОЙ ТЕХНОЛОГИИ
Учебное пособие
Санкт-Петербург Издательство СПбГЭТУ «ЛЭТИ»
2018
1
УДК 621.382.075(07)
ББК З85-04-02я7
З-99
Зятьков И. И., Кривошеева А. Н.
З-99 Базовые процессы планарной технологии: учеб. пособие. СПб.: Изд-во СПБГЭТУ «ЛЭТИ», 2018. 63 с.
ISBN 978-5-7629-2380-4
Рассмотрены одномерные аналитические модели диффузии и ионной имплантации, позволяющие рассчитать профили распределения концентрации примесей в полупроводнике при различных граничных условиях на поверхности полупроводниковых пластин, соответствующих различным технологическим процессам. Приведены примеры расчетов параметров технологических процессов формирования полупроводниковых структур, а также расчет элементов интегральных схем, получаемых методами диффузии и ионной имплантации.
Предназначено для подготовки бакалавров, обучающихся по направлению 200100.62 «Электроника и наноэлектроника».
УДК 621.382.075(07)
ББК З85-04-02я7
Рецензенты: кафедра нано- и микроэлектроники ФГБОУ ВО «Пензенский государственный университет»; д-р физ.-мат. наук, вед. науч. сотр. ФТИ им. А. Ф. Иоффе РАН Е. В. Калинина.
Утверждено редакционно-издательским советом университета
в качестве учебного пособия
ISBN 978-5-7629-2380-4 |
© СПбГЭТУ «ЛЭТИ», 2018 |
2
1. МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ДИФФУЗИИ
1.1. Законы процесса диффузии
Процесс диффузии в полупроводниковой технологии представляет собой направленный перенос атомов примесей, обусловленный тепловым хаотическим движением при наличии градиента концентрации примесей.
Математическое описание диффузионных процессов в виде двух законов, основанных на уравнениях теплопроводности, впервые было предложено в 1855 г. немецким ученым А. Фиком применительно к идеальным газам и растворам. Эти законы и лежат в основе макроскопической теории диффузии
примесей в полупроводниках. |
|
В большинстве случаев диффузия в полу- |
|
проводниках может быть представлена как од- |
|
номерный процесс. На практике это означает, |
|
что глубина диффузии (например, вдоль оси x) |
|
значительно меньше размеров площади, на ко- |
Рис. 1.1. Полупроводниковая |
торой она происходит (в направлении осей |
структура с диффузионным |
слоем: w – толщина |
|
y и z) (рис. 1.1). |
кристалла; xj – глубина |
Рассмотрим уравнение Фика в одномер- |
залегания p-n-перехода |
|
|
ном приближении. |
|
Первый закон Фика характеризует скорость проникновения атомов од-
ного вещества в другое: |
|
|
Jx = −D dN |
, |
(1.1) |
dx |
|
|
где Jx – плотность потока или число атомов вещества, переносимых в единицу времени через единичную площадь, (с·см2)–1; D – коэффициент диффузии,
см2/с; dNdx – градиент концентрации диффундирующих атомов, см–4. Знак
«минус» в правой части (1.1) показывает, что диффузия происходит в направлении убывания концентрации N (совпадает с направлением оси x), тогда как градиент концентрации направлен в сторону увеличения концентрации (противоположен направлению оси x) (рис. 1.2). Необходимо обратить внимание на то, что плотность потока определяется как разностью концентраций dN в точках с координатами x1 и 0, так и расстоянием dx между ними.
3
Второй закон Фика определяет скорость накопления диффундирующей примеси в любой плоскости, перпендикулярной направлению диффузии, как функцию времени. Для одномерного случая закон имеет вид:
dN |
|
d |
dN |
|
|
|
dt |
= |
|
D |
|
, |
(1.2) |
|
||||||
|
dx |
dx |
|
|
||
где |
dN – изменение концентрации диффундирующей примеси во времени. |
|
|
dt |
|
|
|
Коэффициент диффузии численно равен плотности |
|
|
потока атомов вещества при единичном градиенте кон- |
|
|
центрации. Так как диффузионный поток атомов веще- |
|
|
ства идет в направлении выравнивания перепада концен- |
|
|
трации, то коэффициент диффузии D является мерой |
|
|
скорости, с которой система способна при заданных |
Рис. 1.2. Профиль |
условиях выровнять разность концентраций. Коэффици- |
|
|
распределения |
ент диффузии в общем случае может зависеть от кри- |
|
концентрации |
сталлографического направления. Однако свойства таких |
|
примесей |
|
материалов, как кремний, германий, арсенид галлия, имеющих кубическую кристаллическую решетку, можно считать изотропными.
Уравнения (1.1) и (1.2) справедливы в макроскопическом приближении без рассмотрения детальных взаимодействий между атомами примесей и полупроводника. Однако количественные параметры этого взаимодействия учитываются коэффициентом диффузии, представляемым обычно соотношением
D = D0 exp(−∆E / kT ), |
(1.3) |
где ∆E – энергия активации процесса диффузии, эВ; |
k = 8,625 10–5 эВ/К – |
постоянная Больцмана; T – температура, К.
Как следует из (1.3), коэффициент диффузии D сильно зависит от температуры процесса диффузии T. Удобно изображать эту зависимость в полулогарифмическом масштабе, т. е. ln D = f(1/T) (рис. 1.3):
ln D = ln D0 − ∆kTE .
Предэкспоненциальный множитель D0 – значение D при 1/T = 0, что
соответствует интерполяции D в область бесконечно больших температур. Энергия активации процесса диффузии пропорциональна тангенсу угла наклона α прямой (рис. 1.3). Если α1 > α2, следовательно, ∆E1 / kT > ∆E2 / kT.
4
Значение энергии активации для наиболее часто используемых примесей составляет единицы электронвольт. Как коэффициент D0, так и энергия актива-
ции процесса диффузии |
∆E определяются физико- |
|
химическими свойствами вещества, в котором проис- |
|
|
ходит диффузия, а также физико-химическими свой- |
|
|
ствами диффундирующих примесей. В табл. 1.1 при- |
|
|
ведены аппроксимации |
коэффициентов диффузии |
Рис. 1.3. |
примесей замещения в кремнии. |
Температурные |
|
Уравнение (1.2) предполагает зависимость коэф- |
зависимости |
|
коэффициентов |
||
фициента диффузии D от концентрации примесей. По- |
диффузии примесей |
|
лучить решение этого уравнения, т. е. найти зависимость распределения концентрации примесей N (x, t) от координаты за время t, в аналитическом виде невозможно. Для расчета профиля распределения примесей необходимо использовать численные методы с применением ЭВМ.
|
Таблица 1.1 |
Коэффициенты диффузии примесей в кремнии |
|
|
|
Вид примеси |
Коэффициент диффузии |
Алюминий Al |
4,8·exp (–3,44/kT) |
Мышьяк As |
68,6·exp (–4,25/kT) |
Бор B |
0,544·exp (–3,425/kT) |
Галлий Ga |
3,6·exp (–3,49/kT) |
Сурьма Sb |
5,6·exp (–3,92/kT) |
Фосфор P |
3,85·exp (–3,66/kT) |
Если коэффициент диффузии D считать постоянным (такое допущение справедливо в большинстве случаев диффузии в полупроводниках, а для примесей в кремнии (табл. 1.1) это соответствует N ≤ 1019 см–3, то (1.2) примет вид:
dN |
= D |
d 2N |
. |
(1.4) |
|
dt |
dx2 |
||||
|
|
|
Однако и для этого уравнения получить аналитическое решение удается только в некоторых, хотя и очень важных с практической точки зрения, частных случаях.
1.2.Простейшие модели распределения концентрации примесей
вполупроводниковых структурах, изготовленных методом диффузии
Решение уравнения (1.4) в зависимости от начальных и граничных условий дает разные виды распределения концентрации примесей в полупровод-
5
нике. Наибольший практический интерес представляют два частных случая этого решения, которые характеризуют соответственно две стадии диффузии в технологии изготовления полупроводниковых структур. Первый случай представляет собой диффузию из неограниченного источника (загонка примесей), а второй носит название диффузии из ограниченного источника (разгонка примесей).
1.2.1.Модель диффузии из неограниченного источника
Вслучае проведения диффузии из неограниченного источника примесей (загонка примесей) на поверхности полупроводниковых пластин (x = 0) поддерживается постоянная концентрация примесей, т. е. N (0, t) = Ns = const при
любом времени проведения диффузии. Отметим, что в литературе этот способ проведения диффузии имеет различные названия: диффузия при постоянной концентрации примесей на поверхности, диффузия из источника неограниченной мощности и др.
Профиль распределения концентрации примесей по глубине кристалла после проведения процесса диффузии в течение времени t1 при постоянном коэффициенте диффузии D1 (постоянной температуре процесса) описывается с помощью распределения по дополнительной функции ошибок (erfc-распре- деления):
N (x,t) = Ns erfc |
|
x |
|
. |
(1.5) |
|
|
|
|||
|
|
||||
1 |
2 |
D1t1 |
|
||
|
|
||||
Символом erfc обозначена функция дополнения интеграла ошибок до единицы (сокращение от error function complementary). Между нею и функцией ошибок erf существует зависимость:
erfc y = 1 – erf y.
В табл. 1.2 даны значения erfc y, а на рис. 1.4 она представлена в полулогарифмическом и линейном масштабах.
Глубина залегания p-n-перехода xj, иначе говоря, металлургическая граница p-n-перехода, определяется как расстояние от поверхности в глубь полупроводника, на котором концентрация введенных примесей N(x; t) равна концентрации примесей в исходном полупроводнике Nисх, т. е. при x = xj
N(xj; t) = Nисх.
6
|
|
|
|
|
|
|
|
|
|
|
Таблица 1.2 |
|
|
|
|
Значения функции erfc y |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
y |
erfc y |
|
y |
erfc y |
|
|
y |
erfc y |
|
y |
erfc y |
0 |
1,00000 |
|
1,00 |
0,15730 |
|
2,00 |
|
0,00468 |
|
3,00 |
0,00002209 |
0,10 |
0,88754 |
|
1,10 |
0,11980 |
|
2,10 |
|
0,00298 |
|
3,10 |
0,00001165 |
0,20 |
0,77730 |
|
1,20 |
0,08969 |
|
2,20 |
|
0,00186 |
|
3,20 |
0,00000603 |
0,30 |
0,67137 |
|
1,30 |
0,06599 |
|
2,30 |
|
0,00114 |
|
3,30 |
0,00000306 |
0,40 |
0,57161 |
|
1,40 |
0,04772 |
|
2,40 |
|
0,000689 |
|
3,40 |
0,00000152 |
0,50 |
0,47950 |
|
1,50 |
0,03390 |
|
2,50 |
|
0,000407 |
|
3,50 |
0,000000743 |
0,60 |
0,39614 |
|
1,60 |
0,02365 |
|
2,60 |
|
0,000236 |
|
3,60 |
0,000000356 |
0,70 |
0,32220 |
|
1,70 |
0,01621 |
|
2,70 |
|
0,000134 |
|
3,70 |
0,000000167 |
0,80 |
0,25790 |
|
1,80 |
0,01091 |
|
2,80 |
|
0,000075 |
|
3,80 |
0,000000077 |
0,90 |
0,20309 |
|
1,90 |
0,00721 |
|
2,90 |
|
0,000041 |
|
3,90 |
0,000000035 |
Рис. 1.4. Нормализованные распределения концентрации для функции Гаусса
идополнительной функции ошибок (erfc) в масштабах: а – полулогарифмическом;
б– линейном
Количество диффузанта Q, см–2, введенное в полупроводниковую пластину за время проведения диффузии t1 при температуре диффузии Т1 и отнесенное к единице площади диффузионного слоя, определяется выражением:
Q = 2N |
|
|
D1t1 |
|
=1,13N |
|
|
N |
|
x |
. |
(1.6) |
|
s |
|
s |
D t |
s |
|||||||||
|
|||||||||||||
|
π |
1 1 |
|
j |
|
|
|||||||
|
1 |
|
1 |
|
|
1 |
|
|
|
||||
Рассматриваемый случай соответствует таким технологическим процессам, когда примесь на поверхности полупроводника постоянно восполняется взамен продиффундировавшей в глубь полупроводника. Такая ситуация реализуется, например, при проведении диффузии методом открытой трубы, когда диффузант вводят в газ-носитель в виде пара или газа, либо источник примеси предварительно наносят на поверхность полупроводника, как это происходит при диффузии из неорганических пленок, резистов-диффузантов. При этом источник не успевает истощиться за время проведения диффузии. В реальных структурах Ns1 может достигать при данной температуре предельной
7
растворимости, равной 1019…1021 см–3, в то время как исходная концентрация примесей в полупроводнике Nисх = 1014…1016 см–3. Поэтому распределение примесей N(x, t) удобнее изображать не в линейном (рис. 1.4, б), а в полулогарифмическом (рис. 1.4, а) масштабе.
Из уравнения (1.5) видно, что распределение примесей N(x, t) зависит как от температуры T1 (от T1 зависит коэффициент диффузии D1), так и от времени диффузии t1. О скорости диффузии примесей удобно судить, например, по глубине залегания p-n-перехода xj (при этом Na = Nd). Если положить постоянной температуру диффузионного процесса T1 (D1 = const), то по мере увеличения времени диффузии t1 примеси будут все глубже проникать в глубь полупроводника (рис. 1.5, а). Аналогично будет выглядеть распределение примесей, если изменять T1 при постоянном времени проведения диффузии t1 (рис. 1.5, б). Иначе говоря, при изменении любого параметра T1 или t1 получаем качественно одно и то же распределение. Это позволяет сделать следующие выводы.
Во-первых, совершенно очевидно, что будет справедливо и распределение, представленное на рис. 1.5, в. Таким образом, получен важный параметр, характеризующий диффузионный процесс, – произведение D1t1, или пара-
метр 
D1t1 , имеющий размерность длины.
Рис. 1.5. Профили распределения концентрации примесей при загонке: а – постоянная температура; б – постоянное время диффузии; в – постоянное произведение 
D1t1
Данный параметр называют «длиной диффузии», подразумевая, что на этой длине концентрация примесей уменьшится по сравнению с поверхностной Ns1 в е раз (е = 2,72) за время диффузии t1 (по аналогии с диффузионной
длиной неосновных носителей заряда в полупроводниках).
Во-вторых, распределение концентрации примесей не изменится, если, например, D1 увеличится в 2 раза, а t1 уменьшится во столько же раз, или
8
