Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Микроструктуры интегральной электроники

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
20.63 Mб
Скачать

Му рассчитаем зависимость а (а*) от 77, воспользовавшись форму­

лами для надбарьерных токов для структур

без и с р+-слоем и

»ащав функцию распределения

N(1) в виде

N(l) = D ~ l (2 я)"1/2Х

Х ехр[— (7—J)2/2D2g\. В этом

выражении 7

— среднее

значение

Толщины; Dg — дисперсия, характеризующая

величину и

полуши­

рину максимума зависимости N(1). При Z)g<C 1 структура одно­

родна с

 

7, а при Dg^>\ вклад

элементов

площади с

различны­

ми 7 в

полный ток

равновероятен, т. е. структура

максимально

• неоднородна.

В рамках данной

модели

учет

последовательного

Сопротивления

Гб элемента

площади

структуры

сложен, поэтому

: расчет а (77) проводили при

7/<770, когда для каждого

элемента

площади справедлива

экспоненциальная

зависимость 7(77).

Со­

гласно

оценкам

770=ф *о— (Т/Т0) In (ТоГ^Т)

(7’0= 3 0 0 К) •

Выраже­

ния для

а (а*)

от 77 получили из

(2.77)

для

токов

над барьерами

Фо и ф°тах

(структуры

без и с р"1-слоем

соответственно)

 

 

с(-Ч„+(/)Г0/Г

 

' f n N {l)dH _

J ax е- Фщах Т«!Т £ (

3

N (I) di

 

а =

 

 

 

 

 

 

 

 

min

 

 

 

du

 

 

 

 

 

 

 

 

‘min

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е—(Фо+Н)То/г

 

 

 

j"

е

4>ma-xT°/T N (l)dl

 

 

 

J " N (i)d(l)+

 

при 7Х= 0,

 

 

 

 

 

 

 

 

min

 

 

 

 

(2.78)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

"фтах r °/r

d ( — cpjnax)

»,

 

 

 

 

 

 

 

 

 

 

 

 

 

--------------- N (/) d/

 

 

 

 

 

 

а =

J _________________àU

 

 

 

при

 

 

(2.79)

 

 

*max

 

m

T ,T

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

j

 

«'max Г»/Г дг (/) dl

 

 

 

 

 

 

 

a * (-7 7 ) , | -

a

(

-

77),

 

 

 

 

 

 

 

 

(2.80)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^_Й?тах)

_

7C___

7 j _ ________ ?________ \

 

 

 

 

dU ~

 

K - - ï

 

V

V2 (ф,

-

l/) +

K /r

j '

 

 

 

 

Эти зависимости, полученные численным интегрированием

.(2.78),

(2.79)

для

структур

при

различных

Dg, К, 7, фо, к, 7тах и

N d = 2-1016

см -3, изображены па

рис. 2.10,

2.11. На рис. 2.10

за­

висимости

а (77)

рассчитаны

при Dg<C 1

(р+-слой, близкий к одно­

родному), Dg^> 1

(р+'-слой максимально

неоднородный)

и для

од­

нородного р+-слоя при 7 = 7min, 7=7 .и 7=7тахВидно, что при Dgc 1 и для однородного р^-слоя с 7=7 зависимости а от 77 ие отли­

чаются.

Это объясняется тем, что при

1 число элементов пло­

щади с

7=7 велико и определяет вид зависимости а (77). При 7)?> 1

и для однородного р+-слоя с 7=7min зависимости а от 77 отлича­

ются незначительно. Из

зависимостей

а (77) для

промежуточного

случая

неоднородности

р+нслоя,

когда

Dg= 0 ,5

(рис. 2.10), видно,

что изменение параметров к, 7т ах, 7

в пределах

7= 0,5 —2

(рис.

2.10),

7тах= 2 —2,5

(рис. 2.11,в)

и

7i =

0,5— 1,5

(рис. 2.11,6)

вли­

яет в

основном на

величину a,

a

параметров К

и фо в пределах

К = 1 0 —50 (рис. 2.11,0) и фо— 0,3—0,6 (см. рис. 2.11,с) — на за­ висимость а от U. Так, .изменяя К, фо, получаем а — 38— 39 и 35— 39,7 В-1 соответственно. Такие же закономерности характерны и для однородной структуры. Учет неоднородностей уменьшает влияние К, 1 на а (17).

Таким образом, воспользовавшись определенным набором па­ раметров Dg, К, Т, Фо, h, /max, можно смоделировать эксперимен­ тальную зависимость а от U. Очевидно, что в случае незначи­ тельной неоднородности по толщине р+-слоя эту зависимость мож­ но смоделировать подбором К, I и фо для однородной модели. Поэтому для подтверждения (или опровержения) существования неоднородности р+-елоя необходимо исследовать, кроме зависимо­

сти

а (17),

другие

характеристики

структур. Этими

характеристи­

ками могут быть, например, температурные зависимости

а (17,

Т)

и C~2—f\U ).

Д ля рассмотрения

температурных

зависимостей

a (U, Т) провели

численные

расчеты d a /d (llf)

для

однородной и

неоднородной структур и показали, что они

совпадают

с

точ­

ностью до

1%, т. е. температурная зависимость

а (17, Т)

не

может

служить критерием неоднородности структуры.

 

 

 

 

 

 

 

 

Емкость структур с неоднородным по толщине р+-слоем пред­

ставим как сумму емкостей элементов

площади

S C — eS/4nD.

Поэтому полная емкость структур

 

 

 

 

V—1

 

 

 

 

 

 

eS

 

 

N (/) dt

 

max

 

 

 

 

 

 

 

 

 

J,

 

 

.......N (1)-dl ,

-

 

 

 

(2.81)

 

V

.................................. ......J

 

 

 

 

 

h

2 (ф0 — £/) + * /»

h

 

 

 

 

 

 

 

 

 

 

На

рис. 2.12 построены зависимости C~2— f(U )

с параметрами

К,

I и

фо, яелученными из а (17) для

однородной

структуры,

и с

па­

раметрами Dg, К, 1, фо, /| и /m-х, полученными

из а (17)

для

неод­

нородной

структуры. Видно,

что хотя а (U)

примерно совпадают,

зависимости C~2= f ( U ) для

однородной

и

неоднородной

структур

значительно отличаются.

 

 

 

 

 

 

 

 

 

 

 

 

 

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

полупроводником. Моделирование

зависимостей

/ (17) и а (17)

при

 

C~z-t0~f пФ

 

 

 

U<.Ur> для

структур

с неодно­

 

 

 

 

 

 

родным

распределением

Nd в

 

 

 

 

 

 

полупроводнике

вдоль

 

грани­

 

 

 

 

 

 

цы с металлом, определяю­

 

 

 

 

 

 

щим,

главным

образом,

неод­

 

 

 

 

 

 

нородность

 

ф°тах,

проведено

 

 

 

 

 

 

также

для

 

нормального рас­

 

 

 

 

 

 

пределения

7V (ф°тах)

[38].

ИЗ

 

 

 

 

 

 

(2.77)

следует,

что

в

случае

 

 

 

 

 

 

Рас. 2.12. Зависимости С -2

от

U,

 

 

 

 

 

 

рассчитанные для структур с одно­

 

 

 

 

 

 

родным (а) и неоднородным (б) по

 

 

 

 

 

 

толщине р+-слоем:

 

 

6—К-10;

 

 

 

 

 

 

а —К=10;

 

г,-0,7;

фо-0,35 эВ.

 

 

 

 

 

 

/1=0,5;

г» 0,6;

/гаах=2; фо=0,35

эВ

при

Г>=Ш0

яадбарьерного тока через различные участки структуры с неодно­ родным по концентрации легирующей примеси полупроводником влияние участков структур с NdB>Nd «a токи через структуру больше, чем участков с Ndn< .N d. Если характерные размеры не­ однородности полупроводника по Nd d меньше размеров структу­ ры, но больше толщины ОПЗ (£><Cd<S1/2), то, пользуясь общим

правилом

нахождения среднего, для средней величины

прямого

тока ТВр можно

записать [4]: 7Пр = J* Inp(Nd)P (Nd)dN d[ J*P (Nd) X

X dN d]-1,

где P (N d) — функция распределения

плотности вероят­

ности по

Nd. В первом приближении, считая,

что

n ^ f ( U ) ,

т. е.

n = c o n st,

ВАХ

для надбарьерного

тока можно

представить в

виде

 

 

 

 

 

 

 

 

7„p =

^ S «

p ( - | = r ) e Xp ( | § ) .

 

 

 

 

(2.82)

где ф'8 — средняя высота барьера; п — средний

(усредненный по

площади)

коэффициент идеальности

ВАХ. При

U —0

из

(2.82)

находим

 

 

 

 

 

 

 

 

 

t /тч

J* р Мг) ехр [ — (<р'в (Nd)/kT)\dNd

 

 

 

 

(2.83)

фа — —R i Ш------------------------------------------ .

 

 

 

 

 

 

J P (Nd)dNd

 

 

 

 

 

 

Усредненный дифференциальный коэффициент

идеальности

 

 

$P(N d)l(N d)dN d

 

 

 

 

 

(2.84)

J

P (Nd)I(Nd)ln(Nd)\-id N d

 

 

 

 

 

 

 

 

 

 

 

Анализ влияния неоднородности полупроводника по Nd на параметры струк­ туры М-р4 -п-И проведем для трех случаев распределения плотности вероятности

Р(Кл)

в интервале от Nd min до Ndm,\z•

 

 

 

 

а)

нормальное распределение от Namm до Ndmax

 

 

 

f ____ i _

- exp

 

(Nd -

Nd)

при Ndmia < N d< N dmax’

 

 

2a(Nd)*

P (N<0 = { о (Nd) У 2л

 

 

при Nd< A'dmin. Nd> N dmax!

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.85)

б) равномерное распределение в интервале от Nd min до Nd mux

 

 

 

г шах

Ndnita)

*

при

Nd min <C Nj <C Ndшах»

 

(2.86)

 

fM u

 

 

 

 

 

 

 

P (AW

 

 

 

 

при Nd< Nd min> Nd> N dmax!

 

 

lo

 

 

 

 

в) дискретное распределение

 

 

 

 

 

 

P (Nd) = ba (Nd — Ndl) +

(1 -

i) 0 (Nd - Wdi),

 

^ (2-87)

где a(Nd) — дисперсия;

b —отношение участка площади структуры о

А'<г=Л<и

Я общей площади структуры.

 

 

 

 

 

 

 

Очевидно,

что распределение

(2.85)

при cr(iVd)<#<* перехолнт в PlN d)^

«=s0(tfd), что

соответствует

однородной

по

Nd структуре с

Nd==Nd, а ПРИ

tf(Nd)<Nd в (2.86). На рис. 2.13,а показана зависимость I/I(Nd)

от напряжения

для распределения плотности

вероятностей

(2.85) для различных

значений

Рис. 2.13. Рассчитанная зависимость отношения тока 7, протекающего че­ рез неоднородную струкгуру М-р4- п-П, к току 7 через эту однородную структуру от U для различных

o(Nd)

Рис 214 Рассчитанная по (2105) зависимость отношения площади участка структуры М-рл-п—П с кон­ центрацией примеси в полупроводни­ ке Ndu при которой через него про­ текает половина суммарного тока, к общей площади структуры от кон­ центрации Ndi

o(Nd) и 3-1016 см“3. Как видно из рисунка, при o(Nd) >2- 101в см-3 от­ ношение J// (fid) с увеличением напряжения уменьшается, т. е. неоднородность по Na приводит к ослаблению зависимости тока от напряжения. Кривая, соот­

ветствующая

o(Nd) = Ю18,

совпадает с кривой, рассчитанной для распределе­

ния

(2.87).При <y(Afd) =

l*1016 отношение 7/7(#<0»7. Рассчитанные зависимости

от

a(Nd) отклонения параметров <р'ь, fid

для распределения (2 85) при

JV<j =

= 3

- 1016 от соответствующих параметров однородных структур

с Nd~fid

сви­

детельствуют

о том, что

эти отклонения

возрастают с ростом

o(Nd)

и

при

o(iVd) >

Ю18

стремятся

к

насыщению, так как при этих a(fid)

распределение

P (Nd)

(2.85)

переходит

в

распределение

(2.86), не зависящее от (У(Na).

 

 

Дискретное распределение плотности вероятности по Nd (2.87) соответствует наличию двух однородных участков: bS с концент­ рацией Ndi и (1—b)S с концентрацией примеси Nd2- В этом слу­ чае

I = b I (Ndl) +

(1 — b) I (Nd2),

(288)

фв = фв (Nd) +

Афв (Ndl, Nd2),

(2.89)

где Дфв = k T In {b + (1 — Ь) exp [ — Е (Дфё (Ndl, Nd2)/kT]}.

n — b l (Ndl) + (1 - Ь ) I (Nd2) {(b/n,) / (Ndl) H-

+ 1(1 - b ) / n 2] / ( t f d2)} -‘,

(2.90)

где tiu n2 — дифференциальные коэффициенты идеальности ВАХ однородных участков структуры с концентрациями примеси в по­ лупроводнике Ndi и Nd2 соответственно. Как видно из (2.88) —

(2.90), при Ь—>-1

ф'ь-ир'ь (N<u) > п-+п\ и при b-vO, J-*

84

-+ I(N d2), ф'ь-мр'ь (-Мгг), n-*-n2. Если параметры однородных уча етков с Ndi и N,12 таковы, что их вклад в / одинаков, т. е.

b /

(АТ*) = (1 - b) I

(Nd2) =

1/27,

 

(2.91)

то

 

 

 

 

 

 

 

n =

2 («J n2)/(nx + na).

 

 

(2.92)

Значение b = bi/2,

при

котором выполняется

условие (2.91),

зависит от

U, Ndi, Nd2, и при ехреЩ п\Ш « e x p eUln2kT

 

*»/2 = {1

h exp [ - (Дфз (Ndl, Nd2)/k Г)]}-1.

 

(2.93)

Зависимость bi/2= f(N di), рассчитанная по (2.93)

при

/= 1 0 -6

см,

р+ = 2-1018

см-3, JVd2 = 3-1016

см-3, приведена на рис.

2.14. Из

ри­

сунка видно, что для того чтобы ток, протекающий через участок структуры с концентрацией доноров Ndi в полупроводнике, давал вклад в общий ток, протекающий через структуру М-р+-л-П, больший, чем ток, протекающий через участок структуры с кон­ центрацией Nd2, его площадь должна быть больше половины пло­ щади структуры, если Ndi< iN d2 (при Nd2= 3-101®см-3 и JV<ii = 1015 см~3 необходимо условие (?>0,84), и меньше половины площади структуры, если Ndi> N d2 (при Nd2 = 3-1016 см~3 и Ndl = W 7 см~3 достаточно & >0,15).

Следовательно, ток через участки с Ndi > N d2 небольшой пло­ щади по сравнению с площадью структуры может существенно влиять на вид ВАХ. Влияние участков с повышенной концентра­

цией легирующей

примеси Nd2 и пониженной

высотой барьера

<р*о согласно (2.6)

на емкость структуры C0(Nd)

при U = 0, кото­

рую можно представить как сумму параллельных емкостей уча­ стка структуры и структуры, соответствующей остальной площа­

ди, можно учесть, если считать, что Со ~ (ф*о—р )1/2■ Из этой формулы следует, что наличие участков с пониженной высотой барьера должно приводить к увеличению емкости неоднородной

структуры

Со (Nd) по сравнению

с емкостью однородной структу­

ры С0 (N,i).

Однако измеряемая

мостовым методом емкость Сизм

структуры с участками повышенной проводимости будет умень­ шаться, так как

СИЗм = С„ ((1 + r 6/#o)2 + <»2q r | ] - ' ,

(2.94)

где i?o — дифференциальное сопротивление структуры при U = 0;

и — частота переменного сигнала.

Из

(2.94)

видно, что

при

уменьшении R0 величина CH3Mтакже уменьшается.

 

Моделирование зависимостей I(U )

и

a ((/)

при U < .U 0

для

структур с неоднородным распределением глубоких уровней в по­ лупроводнике вдоль границы с металлом позволило оценить ми­ нимальные значения N t.H ^N d для элементов площади структуры, на которых генерационно-рекомбинационный ток является доми­ нирующим и определяет ход a(U ) на рис. 2.7 [35].

Результаты эксперимента. Д ля определения механизма пере­ носа носителей заряда в структурах Al-p+-n-Si, в которых р+-слой неоднороден по толщине, сравнивали экспериментальные зависи­

мости a(U ) и С"2(U) с

моделированными

зависимостями

a(U )

для надбарьерного тока

(рис. 2.10, 2.11)

и C~2(U) (рис.

2.12)

для различных случаев неоднородности р+-слоя. Структуры А1-р+- n-Si площадью 5 —2 2 -10- 5 мкм2 изготовлены на n-Si (111) с р = 1 0 ‘1-см. Из рис. 2.11 видно, что зависимости а (£/) характеризуют­ ся двумя участками. В интервале напряжений от 0,1 до 0,2 В на­ блюдается незначительное уменьшение a, a при напряжениях £ /> > 0 ,2 В — более резкая зависимость а от U. Значение высоты барьера q°raax при U = 0 определяли по току ia : ф°тах=0,53—0,54 эВ, значение высоты барьера <р0 — по току ts для структур Al-n-Si, полагая, чго после образования р+-слоя значение ф0, определяе­ мое Г,С кремния, практически не изменится. Рассчитанное значе­

ние <| о -

0,42 эВ.

Сравнение экспериментальной зависимости <x(U)

с оаесчитанными

для однородного и

неоднородного по толщине

р-' -слоя

(рис. 2.10, 2.11) показывает,

что экспериментальная за ­

висимость а (U)

может быть описана моделированной для одно­

родного

р+-слоя

при 1 - 9 - 10-7 см, /(= 1 7 , фо=0,35 эВ; а для не­

однородного р+-слоя при Dg—0,5, /= Ы 0 ~ 6 см, /(= 10, ф0 = 0,35 эВ, lmax = 4-10“6 см, /1=6-10-7 см. Чтобы определить, какой из этих случаев имеет место для надбарьерного тока через структуру, проверили, удовлетворяют ли значения I и К для однородного р+- слоя экспериментальной ВФХ (рис. 2.12). Д ля этого сравнили ф*0,=0,82—0,83 эВ, определенное по отсечке прямой C~2—f(U ) на оси напряжений, с рассчитанным значением ф*оР для однород­ ного р+-слоя: ф*оР= ф 0 + /(/2/2. Расчет дает ф*ор» 0 ,5 эВ, что не совпадает со значением ф*оэ. Согласно теоретическим оценкам ф*о для структур с моделированным неоднородным по толщине р+- слосм ф*0р » 0 ,7 эВ, т. е. этот случай достаточно хорошо описыва­ ет результаты эксперимента: a —f(U ) и C~2—f(U ). Следователь­ но. для структур Al-p+-n-Si с неоднородным по толщине I р+-сло­ ем, когда тип неоднородности I учтен введением функции распре­ деления величины / в виде нормального распределения при запи­ си ВАХ (2.77), основным является надбарьерный ток, как и для структур с однородным р+-слоем.

Экспериментально влияние тока через участки структуры с повышенной концентрацией примеси в полупроводнике на вид ВАХ структур Al-p+-n-Si изучали, моделируя неоднородное распределение Sb в Si (Ад=3• 10'6 см-3) ме­ тодом электрохимической имплантации Sb в Si [40].

Наблюдали выход оже-электронов Sb на участках поверхности Si пло­

щадью —0,3% от общей площади структуры. Создавали структуры Al-p^-n-Si.

Емкость однородной структуры при (/=0

составляла C0(Nd) — 11 пФ, а неод­

нородной Co(Pd) —7 пФ. В структурах с

моделированной неоднородностью по

Ns наблюдается увеличение токов при небольших прямых напряжениях и при

обратных

U^ 8 В. Анализ

ВАХ проведем, используя (2.83) и (2.84). Как сле­

дует из

этих формул, для

структуры с моделированной неоднородностью по

Nd в Si при наличии двух однородных участков площадью bS с Ndi н (1—b)S с Nd2 надбарьерные токи через эти участки при прямых напряжениях £/wp<

Сср'п отличаются из-за разницы площадей и высот барьера

ср'вОУ^) и

(A'd2):

:

< ф /в(Л^(и). Поэтому при небольших прямых

напряжениях

(U ^

^

0,35 В) доминирующим является ток через участок с повышенной концентра­

цией примесей Nd2. При Unv^q>'B(Nd2) рост тока через участок структуры Л^а замедляется, так как U(Md2) = (Ыаг)Гб (гб — последовательное сопротивление объема полупроводника) и при r6>>R(Mdi) доминирующим становится ток, про­ текающий по остальной площади структуры с Ndi в Si. На этом участке ВАХ структуры с моделированной неоднородностью по Nd незначительно отличается от ВАХ структуры с однородной концентрацией Nd. По ВАХ для однородных структур определили <p'b=<0,79 эВ. Для участка структуры с Nd%высота барь­ ера (р'в—ОДТ эВ. Увеличение обратных токов для структуры обусловлено по­ явлением участка с повышенной проводимостью.

Экспериментально наблюдаемое уменьшение Co(Nd) для структур с моде­ лированной неоднородностью, как следует из (2.94), обусловлено уменьшением дифференциального сопротивления структур iRo(Nd) из-за наличия участка структуры с повышенной проводимостью.

Таким образом, экспериментально подтверждены результаты математического расчета параметров ВАХ при моделировании не­ однородности по концентрации примесей в я-полупроводнике структуры металл-р+-/г-полупроводник, имеющий участок с по­ вышенной концентрацией Nd2, подключенный параллельно струк­

туре, соответствующей

остальной площади.

Д ля найденных

<р*о(Мп) и

Ф *о(^о) из

рис. 2.14 находим, что площадь участка

с Nd2>N d\

составляет

не менее 0,25% от общей

площади струк­

туры. Результаты расчетов параметров ВАХ при других законах распределения концентрации примесей в полупроводнике — нор­ мальном и равномерном от Nmin до А/щах (рис. 2,14) — могут быть использованы при изучении причин отклонения ВАХ от первона­ чальных при внешних воздействиях на структуру (радиация и др.).

2.7. ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ОПЗ В М-п(п+)-п+(п)-П СТРУКТУРАХ

Параметры ОПЗ полупроводника зависят от закона распреде­ ления легирующей примеси в ОПЗ (§ 1.1). Поэтому прежде всего рассмотрим законы распределения примесей в п (пf) -слоях на по­ верхности п+(п) -полупроводника.

Законы распределения примесей в полупроводнике. Различные распределения концентрации легирующих примесей в полупроводнике получают при примене­ нии следующих методов: эпитаксиального выращивания; диффузии; переразделения примесей на границе полупроводник—внешняя фаза за счет гетерирования внешней фазой (испарения, термического окисления); ионного внедрения [1]. Например, при эпитаксиальном выращивании «-слоя на /^--полупроводни­ ке (подложке) распределение примесей определяется их диффузией из сильно­ легированной подложки, испарением примесей и некоторыми другими фактора­ ми. Если скорость роста слоя постоянна и коэффициент диффузии примеси ио-

стоянен и одинаков для

подложки

и и-слоя, то примесь в структуре

п-слой-

пл - полупроводник распределяется по закону

 

п (х) =s

ехр а2 crfc а

Ь + с-ехр (b— a)2erfc (b— а)-|-

 

Ь +

exp (b— a +

c) erfc (b -a-\-c)

(2.95)

Здесь no — концентрация примесей

в подложке; a~xj-2~\/ Qi\ b-

 

— W\/Qt\

c~ht~\/ Qt (G —

скорость роста эпитаксиальной пленки, h

- ско­

рость испарения примесей; Q

-- коэффициент диффузии примеси в подложке

и пленке). Концентрация примесей на границе эпитаксиального слоя полупро­

водника ns может

быть вычислена из выражения

(2 95)

при х = W.

 

Распределение

примеси в и-н+-структуре при

h— 0,

построенное

с исполь­

зованием (2 95) для разных отношений Wj'ŸQt,

показано на рис.

2.15 [1].

На этом рисунке показано также распределение примеси в и-и1-структуре при разных значениях концентрации примеси в эпитаксиальном стое nt, определя­ емых легированием растущего эпитаксиального слоя только извне, т е. без уче­

та диффузии примеси

из

подложки в

эпитаксиальный

слой. Видно, что

при

и (W ïY ql)^

1 распределение

примеси

с достаточной

точностью

ап­

проксимируется линейным

законом. При

 

и

Q /)» l

распределение

примеси с хорошей точностью аппроксимируется ступенчатым законом

 

п(х) = п1 при 0

 

< W; п(х) = п0 при

x>W .

 

(2.96)

Для промежуточного слоя распределение примеси можно аппроксимировать не­ сколькими прямыми.

Таким образом, если ОПЗ в п-п^ -структуре расположена в эпитаксиальном

я-слое, то для расчета

ее ширины при выполнении соотношений ns>ni,

(WVVQO ^

1 можно воспользоваться линейным, а при выполнении соотноше­

ний n ^ tiu

[W rY Q t)^\

— равномерным распределением примеси. Если обед­

нение примесями приповерхностной области полупроводника достигается в ре­ зультате высокотемпературного окисления, при котором примесь из полупро­ водника переходит в окисел, то закон распределения примеси в этой области

может быть представлен как

 

 

 

 

 

 

п (*) =

n0 -f (ns— n0) erfc

~УК21 + х]

( Л ^ \ х/2 I - 1

х

xjerfe

Л (Л * V '2

 

» [ т

\

Q !

J

(2.97)

 

 

 

 

 

 

Г

г

 

окисления

 

(dд

- толщина пленки

Здесь Kz — постоянная

в законе

 

окисла), b

постоянная,

которая

при окислении

кремния равна 0,44 [1]. За­

висимость

отношения п(х)1п0

от значения xl~]/nQt

при и8< п 0

при небольших

значениях п(х)/щ может быть аппроксимирована прямой. При линейном при­ ближении ПС х представим как

n(x) = ns -\-x(n0 — ns)iA при х ^ А ; п(х)^=п0 при х> А .

(2.98)

При /?<?<<по

(обедненный слой) справедливы следующие выражения для толщи­

ны обедненного слоя А и пв:

 

А =

l /л

Q t ехр (62 Кг/4Q) erfc [(6/2) (АГг/(?)1/2 ];

(2.99)

ns =

n0{l -h ( jt!/г/АС)ехр(&2/C2/4i?) erfc[(6/4) (AT3/Q>1/21/2 (/C.2/C?)l/2]}—1 •

(2 100)

lit-

 

 

 

Рис. 2 15 Распределение примесей в структу­ ре эпитаксиальный высокоомный слой толщи­ ной W на низкоомном полупроводнике

Из (2.100) видно, что при обеднении пу­ тем термического окисления толщина слоя А растет с увеличением времени окисления t и отношения /G/Q. Концентрация примеси в при­ поверхностном слое полупроводника ns не за­ висит от времени окисления и уменьшается с уменьшением коэффициента К

Таким образом, для рассмотренных мето­ дов изменения концентрации примеси в ns-слое

«в-л+-лолупрово1Дпика распределение примеси может быть аппроксимировано ли­ нейным (2.98) или ступенчатым (2 96) законам.

Параметры ОПЗ в я(п +)-п+(м) полупроводнике. Если концен­ трация примеси в полупроводнике задается законом (2.98), то в результате интегрирования уравнения Пуассона получаем

F (х) = (е/е0 г) {(х2/2 А) (п0 — na) + nsx — [(D2/2 А) (п0 — па) +

+ ns D]},

 

 

 

 

(2.101)

<р (х) = (еа/е0 е) [(п0 — ns)/Ds + *3)/6 A + ns (ЕР + х2)/2 —

 

 

- D 2 (п0 — па) х/2А — па D J,

 

 

(2.102)

D = (2/3) (D/Л)3 (1 —па/п0) + D *na/n0.

 

(2.103)

Если слой Шотки лежит в области х < А , то при D/A^>nsln0 и,

следовательно, ns<C«o

 

 

 

 

 

D-= (3 ЛМ/2)1/3.

 

 

 

(2.104)

При выполнении того же условия

и при D/A<^.(n^lno) (1—ns/no)~l

справедливо соотношение

 

 

 

 

D = (п0/па)112 L.

 

 

 

 

(2.105)

Из (2.104) и (2.105)

видно, что

D > L

при ns< n 0 и

D < L

при

«8> rt0. Величины

Е°2

и ф0—eU

можно

получить из

выражений

(2.101) и (2.102)

при * = 0 .

 

 

 

 

Зависимости напряженности поля Е (х ) и потенциала ф(х)

для

рассматриваемого случая показаны штриховыми линиями на рис. 2.16.

Если

концентрация примесей

в

области

полупроводника

0 <

< л г < /

распределена

по

произвольному

закону,

а

в области

К

< .x<zD

является постоянной

п(х)=П о,

то

Е°

и

фо—eU можно

представить как

 

 

 

 

 

 

 

 

 

 

 

£0 __ _ еп0

D - 1

+ J

 

dx 1

,

 

 

 

(2.106)

 

808

п0

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

Фо— eU --- е2 п0

D2 — /2

Г

^

Xdx

 

 

 

(2.107)

 

 

+

 

 

 

 

 

808

 

 

 

о

"0

 

 

 

 

 

 

Рис. 2.16. Напряженность поля и ход потен­ циала в слое Шотки:

------

равномерное распределение примеси;

линейное распределение примеси

При

выполнении

неравенств

I

'

 

I

[п (х)/n0]x d x > 1212

J \n (x )/n 0]d x '> l и |

о

 

 

о

 

 

значения Е°2 и фо—eU при заданных

величинах

D

будут

превышать соот-

в х ветстъующие значения для равномер­

ного распределения примесей в облас­

ти 0 < x < D . Если выполняются

обратные

неравенства,

то соотно­

шение между вышеприведенными величинами будет противопо­

ложным.

Это означает,

что при средних значениях

концентрации

п (х )< п п

наблюдается

увеличение толщины слоя

Шотки

(обед­

ненный примесями слой у поверхности), а при ti(x )> n 0 ею

умень­

шение (обобщенный примесями слой). Кроме того, из (2.107) видно, что r рассматриваемом случае изменяется зависимость ширины слоя Шотки от фо—eU.

2.8. ВОЛьГ-ФАРАДЯЫЕ ХАРАКТЕРИСТИКИ М-п^-п-П СТРУКТУР ПРИ Ш'ХКНГШШ РАСПРЕДЕЛЕНИЯ ПРИМЕСИ В ПОЛУПРОВОДНИКЕ

PaccMoipiiM ВФХ структур МП для случая усиления зависи­ мости С от U в заданном интервале напряжений, который реали­ зуется в М-Я‘-я-Н структурах [1]. Усиление зависимости С от U необходимо для оптимизации параметров СВЧ-диодов (§ 2.10)

Если распределение примеси в я'-слое линейно, то

С яе0 е а2 {I 2 +

(А213) 11 — (я8/я0) ] 1/2

(2.1081

при О й Л и

 

 

С — яе0 s a2JL (я)

 

(2.109)

при D > ;4 , л3»Л о

(а — радиус структуры).

 

Из приведенных выражений видно, что для обеспечения мак­ симального изменения емкости с напряжением необходимо вы пол • нение условия D й А. Например, при Л2Яз./Зяо->фо/2яе2яо получаем более сильную зависимость емкости от напряжения, чем при меньших прямых напряжениях.

При ступенчатом распределении примеси в полупроводнике усиление зависимости С от U достигается при W ^ L . Зависимость емкости от напряжения для такого закона распределения имеет вид:

С -п& два2 {L2 — W2 [(я8/я0) - 1]}-» /2.

(2.110)

Если Лз/лоЭ-1, то W2nsJn0 приближается

к значению фо/2ле2яо

и емкость МП структуры С сильнее зависит от напряжения, чем для равномерного распределения примеси.

Соседние файлы в папке книги