Материал: TRfDAZqnfW

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

При термическом окислении неоднородно легированного слоя константы окисления становятся функциями глубины прокисления кремния 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(aB 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

Источник: https://studfile.net/preview/16438825/