При термическом окислении неоднородно легированного слоя константы окисления становятся функциями глубины прокисления кремния x = ay (a = 0.45), а через глубину и функцией времени окисления. Дифференциальное уравнение (5.10) становится нелинейным вида
dy = f ( y, t ), dt
его необходимо решать численными методами.
Расчет концентраций носителей заряда n и p в (5.12) проводится по концентрации примеси в легированных слоях с использованием выражений (5.5)–(5.9).
5.5. Перераспределение примесей при термическом окислении
При термическом окислении происходит перераспределение примесей между полупроводником и растущим оксидом. При невысоких скоростях окисления, меньших скорости продвижения диффузионного фронта, отношение концентраций примеси в оксиде C(SiO2) и кремнии C(Si) на границе раздела устанавливается постоянным, равным коэффициенту сегрегации:
m = C (SiO2 ).
C (Si)
Если энергия образования оксида примеси больше, чем оксида кремния, то m < 1, примесь не захватывается оксидом, а оттесняется в кремний, накапливаясь на границе раздела (примеси P, As, Sb в Si). Если энергия образования оксида примеси меньше, чем оксида кремния, то наоборот: m > 1, примесь захватывается оксидом, поверхность кремния обедняется примесью (примеси В, Al в Si). Значения коэффициентов сегрегации для основных легирующих примесей в кремнии приведены в прил. 5.
В системе с движущейся границей раздела Si–SiO 2 x = y – aw, где w –
толщина оксида, a = 0.45 – отношение толщины окисленного кремния к толщине образовавшегося оксида. При постоянном коэффициенте диффузии примеси в кремнии D = const, а в оксиде Dокc = 0 уравнение диффузии принима-
ет вид [14]
∂С = D |
∂2C + a |
dw |
∂C . |
(5.13) |
||
|
||||||
∂t |
∂x2 |
dt ∂x |
|
|||
Начальное и граничное условия имеют вид |
|
|
|
|
|
|
C ( x, y ) = C0 ( y ) при y ≥ 0, D |
∂C (0, t ) = ( m − a ) |
dw |
C (0, t ). |
(5.14) |
||
|
||||||
|
∂t |
|
|
dt |
|
|
16
При параболическом законе окисления ( w = 
Bt ) задача имеет точное решение в интегральной форме [15]:
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
( z − x − a Bt ) |
||||||||||||||
|
C ( x, t ) = |
|
|
|
|
|
|
∫ |
C |
( z ) exp − |
|
|
|
|
|
|
|
|
|
− |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
2 πDt 0 |
|
|
|
|
|
|
|
|
4Dt |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
m Bt − z |
|
|
|
za B D |
( z + x + a Bt ) |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||
− |
|
|
|
|
exp |
|
|
|
|
|
|
|
exp − |
|
|
|
|
|
|
|
|
dz. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
m Bt + z |
|
|
|
|
|
Dt |
|
|
4Dt |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
В случае однородного начального распределения легирующей примеси в
кремнии (C(x, 0) = C0 = const) и отсутствия начального оксида (w(0) = 0) рас-
пределение примеси в кремнии после термического окисления описывается аналитическим выражением [14]
|
|
|
C ( x, y ) = C + Aerfc |
|
|
x |
|
|
+ a |
|
|
|
, |
|
|
|||||
|
|
|
|
B D |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
||||||||||||||||
|
|
0 |
|
|
Dt |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
(a - m) |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
C |
|
|
B D |
|
|
|
|
||||||||||
A = |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
||
|
1 |
exp(-a2 B D) |
- (a - m) |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
× erfc |
a2 B D |
|||||||||||||
|
|
B D |
||||||||||||||||||
|
|
|
|
|
||||||||||||||||
|
|
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где x – координата, отсчитываемая от границы раздела Si–SiO 2. Распределе-
ние концентрации примеси в оксиде по глубине при этом постоянно:
Cокс = m C0 + Aerfc(a
B
D ) .
В случае неоднородного исходного распределения примеси по дополни-
тельной функции ошибок (после первой стадии диффузии при C(0, t) = Спов = = const и D0 = const)
C ( x, t ) = C |
erfc |
|
|
x |
|
|
|
2 |
|
|
|
||
D t |
|
|||||
пов |
|
|
||||
|
|
|
|
0 0 |
|
|
распределение примеси в кремнии после термического окисления описывается аналитическим выражением [14]
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
C ( x, t ) = |
|
|
Q |
|
|
|
|
x |
|
|
|
|
B |
|
|
|
|
||||||||||||||
|
|
|
|
|
exp − |
|
|
|
|
|
+ a |
|
|
|
|
|
|
− |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
πDt |
|
|
|
|
|
2 |
|
Dt |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
(m - a) |
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
+ a |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
B D Q erfc |
|
|
|
|
|
|
|
|
|
|
B D |
|
|
|
|||||||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
- |
|
|
|
|
|
|
|
|
|
|
|
|
|
Dt |
|
|
|
|
|
|
|
|
|
|
|
|
, |
||||||
1 - (m - a) |
|
|
|
|
|
|
exp(a |
|
|
|
|
|
|
)2 erfc(a |
|
|
|
|
) |
||||||||||||||
|
πB D |
|
B D |
|
B D |
|
|||||||||||||||||||||||||||
17
где Q – напомним, количество примеси в легированном слое, введенное на первой стадии диффузии:
Q = 2Cпов
D0t0 .
π
5.6.Решения уравнения диффузии при D = const
Вряде операций технологии ИМС используется диффузионное или ионное легирование с низкими концентрациями легирующих примесей (C < ni), напри-
мер при создании карманов и подлегировании каналов в технологии МОП ИМС, создании резисторов и базы в технологии биполярных ИМС. В этих случаях коэффициент диффузии легирующей примеси не зависит от концентрации (D = const) и уравнение диффузии имеет вид
∂C = D |
∂2C . |
(5.15) |
∂t |
∂x2 |
|
В общем случае при произвольном начальном распределении примеси C(x, 0) = C0(x) дифференциальное уравнение (5.15) имеет решение в инте-
гральной форме для бесконечного пространства (−∞ < x < + ∞) в виде
|
C ( x, t ) = |
|
|
1 |
|
|
+∞ |
0 |
( z )exp |
|
|
( z − x)2 |
|
|
||||||
|
|
|
|
|
|
|
∫ |
− |
|
|
dz, |
|
|
|||||||
|
|
|
|
|
|
|
|
C |
|
|
|
|
||||||||
|
|
|
|
|
2 |
|
πDt −∞ |
|
|
|
|
4Dt |
|
|
|
|||||
а для полубесконечного пространства ( 0 ≤ x < +∞ ) в виде |
|
|
||||||||||||||||||
|
1 |
|
∞ |
|
|
|
|
|
|
( z − x) |
2 |
|
|
( z + x) |
2 |
|
||||
C ( x, t ) = |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
∫ |
C |
( z ) exp − |
|
|
|
± exp − |
|
|
dz, |
||||||||
|
|
|
|
|
|
|
||||||||||||||
|
2 πDt |
|
0 |
|
|
|
|
4Dt |
|
4Dt |
|
|
||||||||
|
0 |
|
|
|
|
|
|
|
||||||||||||
где z – переменная интегрирования, знак плюс между экспонентами соответ-
ствует отражающей границе при x = 0 |
dC ( x, t ) |
= 0 |
, а знак минус – погло- |
|
|||
dx |
|
|
|
щающей границе при x = 0 (С( x, t ) = 0 ). |
|
|
|
В отдельных частных случаях уравнение диффузии (5.15) имеет точные решения в виде аналитических или известных табулированных функций: erf(z),
erfc(z), exp(–z 2), где z = |
|
x |
|
|
|
|
|
|||
|
|
|
|
. |
Функция ошибок определяется следующим |
|||||
|
|
|
|
|||||||
2 |
|
Dt |
|
|
|
|
|
|||
выражением: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
z |
|
|
erf ( z ) = |
|
∫ exp(− y2 )dy. |
|||||||
|
|
|
|
|||||||
|
π |
|||||||||
|
|
|
|
|
|
|
|
0 |
||
18
Дополнительная функция ошибок дополняет функцию ошибок до 1:
erfc( z ) = erf ( z ) – 1.
Дополнительная функция ошибок определяет распределение примеси по глубине:
1.При диффузии из бесконечного (постоянного) источника (C(0, t) =
=Cпов = const) на этапе загонки легирующей примеси:
С( x, t ) = C |
|
erfc |
|
|
x |
|
. |
||||||||
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|||||||||||
|
|
пов |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
2 |
|
Dt |
||||||
2. При диффузии из ограниченного источника с количеством примеси Q |
|||||||||||||||
dC ( x, 0) |
= 0 |
|
|
|
|
|
|||||||||
и отражающей границе при x = 0 |
|
|
|
|
|
|
на этапе разгонки: |
||||||||
|
|
dx |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
C ( x, t ) = |
|
Q |
|
|
|
|
|
x |
2 |
|
|
||||
|
|
|
exp |
|
|
|
. |
||||||||
|
|
|
|
|
|
|
|
|
|||||||
πDt |
|
|
|
|
|||||||||||
|
|
|
|
|
|
4Dt |
|||||||||
3. При диффузии в неограниченном теле из полубесконечного пространства с равномерным распределением примеси (C(x, 0) = C0 при x ≤ 0):
C ( x, t ) = |
С |
|
x |
||||
0 |
erfc |
|
|
|
|
. |
|
|
|
|
|
|
|||
|
2 |
2 |
|
Dt |
|||
Этот случай имеет место в планарно-эпитаксиальной технологии ИМС при выращивании слаболегированных эпитаксиальных пленок на сильнолегированных подложках или скрытых слоях.
4. При диффузии из полуограниченного равномерно легированного тела (C(x, 0) = C0 при x ≥ 0) наружу (испарении) со связывающей (поглощающей)
границей (C(0, t) = 0): |
|
|
|
|
|
|
|
C ( x, t ) = C |
erf |
|
|
x |
|
. |
|
|
|
|
|
||||
|
|||||||
0 |
|
|
|
|
|
||
|
|
2 |
|
Dt |
|||
5. При наличии активационного барьера для испарения примеси граничное условие на поверхности имеет вид
D ∂C (0, t ) = sC (0, t ),
∂x
где s – коэффициент испарения примеси. Решение уравнения диффузии для полуограниченного тела с равномерным исходным легированием имеет в этом случае вид [16]
19
|
|
|
x |
|
|
|
|
|
|
|
|
sx |
|
|
2t |
||||||
C ( x, t ) = C0 |
x |
|
|
|
t |
|
|
s |
|||||||||||||
erf |
|
|
|
|
+ erfc |
|
|
|
+ s |
|
|
exp |
|
+ |
|
|
. |
||||
|
|
|
|
|
|
|
D |
|
|
||||||||||||
|
|
|
|
|
|||||||||||||||||
|
|
|
2 Dt |
2 Dt |
|
|
|
D |
|
|
D |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Возможность потерь и перераспределения примесей вследствие испарения необходимо учитывать при отжигах полупроводниковых подложек и легированных слоев с незащищенной поверхностью в вакууме, нейтральных средах (Ar, N2), а также при предэпитаксиальной обработке подложек в водороде.
Параметры температурной зависимости коэффициентов испарения и определяемых совместно с ними коэффициентов диффузии легирующих примесей в кремнии приведены в прил. 6.
5.7.Диффузионная загонка
Втехнологии ИМС диффузия легирующих примесей осуществляется обычно в две стадии. Первая стадия – загонка – проводится из примесносиликатного стекла с постоянной поверхностной концентрацией, равной, как правило, предельной растворимости примеси в кремнии при температуре
диффузии: Cпов = Cпред (см. прил. 3). Решение уравнения диффузии из бесконечного источника при D = const и Спов = const имеет вид
С( x, t ) = C |
erfc |
|
|
x |
|
. |
(5.16) |
|
|
|
|
|
|||||
|
||||||||
пов |
|
|
|
|
|
|
||
|
|
2 |
|
Dt |
|
|||
Дополнительная функция ошибок erfc в (5.16) |
аппроксимируется следу- |
|||||||
ющими аналитическими выражениями: |
|
|
|
|
|
|
|
|
erfc( z ) ≈ exp(−( z + 0.3)2 )
с погрешностью не хуже 18 % при z ≤ 3 . Наиболее точная аппроксимация дополнительной функции ошибок с погрешностью не более 1.5 · 10–5 % при любых z дается в [17]:
|
|
|
erfc( z ) ≈ (at + bt2 + ct3dt4 + et5 )exp(−x2 ), |
где t = |
|
1 |
, p = 0.3275911, a = 0.2548296, b = –0.2844967, c = 1.4214137, |
|
|
||
|
|
||
1 |
+ pz |
||
d = –1.4531520, e = 1.0614054.
Во многих технологических операциях изготовления ИМС используется диффузионное или ионное легирование с высокой концентрацией легирующей примеси (C > ni), например при разделительной диффузии, создании скрытых слоев и областей эмиттера, подлегировании коллектора в техноло-
20