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

книги / Математические модели элементов интегральной электроники

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

тывающцх изменение геометрических размеров на технологических операциях [15]:

1 = —2*ф— 1.6JC/, г = г т+ 2хф,

(4.5 5

где LT, ZT—топологические (закладываемые в фотошаб­ лон) размеры канала; x j—глубина р—л-перехода диф­ фузионных областей; ОД*^— глубина «боковой» диффу­ зии; Хф^\—1,3 мкм — изменение размера на операциях фотолитографии и тоавления.

Рис, 4.9. Топология интегрального МДП-транзнстора:

-----изменение размеров после проведения определяющих технологических

операций.

На рис. 4,10 показаны экспериментальные и смоде­ лированные характеристики МДП-транзисторов с дли­ ной канала 6 и 17.мкм при разных температурах. Точ­ ность моделирования статических характеристик при ис­ пользовании рассмотренной модели приближается к 10%.

В [16] при построении физико-топологической модели МДП-транзистора используется другой вид аппроксима­ ции зависимости подвижности от поля

Р*рэ = щ/(1 -\~Ех/ EKXI Ey/EKffi),

(4.56)

где Ex=[Ex(fi) +Ех(хк)]12 —усредненное по сечению ка­ нала поперечное поле; E1{xi и Е1{уi — постоянные аппрок­ симации. Нетрудно показать, что усредненное поле £*= =[2£*(0) —Ср/8пео)]/2. Используя (4.25) и (4.20), получаем

Е - — и'а — У— ВкУ Т .

(4.57)

БП

2*д

'

'

Используя остальные Допущения предыдущего разде­ ла, подставляя аппроксимацию (4.57) в (4.56) и затем в выражение для тока (4.30) и интегрируя его по длине канала с учетом граничных условий (4.29), получим сле­ дующее выражение, неявно определяющее /с:

т___________ \^CAZUKy_________ w

с " “ L[UKy + U c - U a - I c(Rc + RB)] А

ис-*с*с

 

I

U , - U t + UB - U - B o t V i , + U

 

Unci + U t - U t + U B - и - В о е V b + u

X

 

с^н

 

и с 'с*с

 

 

dU

-d.u

, (4.58)

J

Uo-\-Ug— U — Вое V^o^U

*М+'Л

 

где £/Ka:i—2£к*1еп/ед; Uliv=LEvyi. Для численного реше­ ния этого уравнения может быть использован итерацион­

ные. 4.10. Экспериментальные

(О О О )

и смоделированные (-----------)

по уравнениям (4.31), (4.40),

(4.45) выходные характеристики М Д П -

транзисторов с параметрами лгд=1100

А, Л ^ = 5 *1 0 15 см-3 и с разной

длиной канала:

 

 

а)

L -I7

мхм, Z=13,5 мкм, ЛГП0В- 811 см -а, ------ 25*С,---------100°С;

б)

L-6

мкм, Z=50 мкм, ЛГП0В= 6- 10" см -2.

 

ный метод Ньютона [см.

(4.33) — (4.35) и (4.37)]. Пола­

гая /?с=/?и=0 и интегрируя (4.58), получаем следующее явное выражение для тока стока в крутой области ха­ рактеристик:

т

P OC Z Z U KX I U KU чу

/ с _ _

£------------ А

X —

*UcУд Ч~4^ос(га —Zi) 4-

 

2С/к*» Г

Zz —Сг

— Сг In

Z2--- C \ ]

Уку +

C t In

Zi —Сг

Z i-C T j

 

Т ^ с Г ['

 

+ с . - с г

[ (6* +

biC,)ln Z - c \ — <Ь г + biC ^ ,п ? , - е г ]

 

 

 

 

(4.59)

где

Zi =

у<р9-\-Цн'

гГ ~ У ъ + У с\

 

 

U0

UB-}-Уо -f- UKXU

 

C I — ■ 1 /2 ( Boc- j -

у В2ос - | -

4 a 2) ;

 

C2 =

1/2 (Boc -

]/j32oc +

4a2);

 

b\ =

UKX\ -|- 2B2oc; bz =

2B0Ca2.

В пологой области характеристик для описания за­ висимости тока стока от управляющих напряжений ис­ пользуется выражение (4.40) /с (^с> t/отс) =/*с/ (1— —/отс/В), учитывающее уменьшение эффективной длины

канала за^ счет расширения

Истк

3т6ор

Спшн

обедненной

области отсечки

 

 

 

/0тспри

увеличении напря­

 

 

 

жения на стоке. Для расчета

 

 

 

/отсв [16] получены форму­

 

 

 

лы, учитывающие влияние по­

 

 

 

движных

носителей тока

в

 

 

 

области

отсечки на

ширину'

 

 

 

этой области (рис. 4.11) при

 

 

 

следующих

допущениях:

в

Рис. 4.11. Форма области отсеч­

области

отсечки подвижные

ки, принятая в [16] для вычис­

носители

 

имеют

постоян­

ления /отс.

состав­

ную.скорость Унас;

под действием

нормальной

ляющей поля носители выталкиваются из поверхностной области и поглощаются стоковым р—л-переходом глуби­ ной Xj; во всей обедненной области отсечки носители имеют пост/отоянную плотность; распределение потенциала в области может быть описано одномерным уравне­ нием Пуассона. С учетом этих допущений можно записать

<4®»

где /Wc/(£*oc); *ос=(Xj—Xi(/ioTc )у — глубина обед­ ненной области отсечки, линейно увеличивающаяся с ро­ стом у (рис. 4.11). Граничные условия для этого урав­ нения следующие: у левой границы области отсечки‘на­ пряжение равно £/ото; с увеличением у напряжение воз-

203

растает и у границы со стоком становится равным Uc

U(y = 0 ) = U o r c ,

dU

Vi\ас

U(y = L) = Uz.

dy у

- И

 

(4.61)

 

 

 

При этих условиях решение уравнения Пуассона (4.60) записывается в виде

_ // _ ^Д. /2

1 J _________ f i £ ±

L _ V

 

U' • — 1/отс — _ _ «* о*

[ ‘ yNfl.Zt)Hac(x>j — Лк)2*

 

епе0

 

х ( ш ^ -

+ V )] + /0TCJ&C’

 

(4.62)

 

 

где Ec=vuacl\io —электрическое поле в точке

отсечки

(#=0). Выражение (4.62)

по сравнению с

(4.46)

учиты­

вает влияние глубины стокового рп-перехода на 10Тс, а следовательно, и на величину выгодного сопротивления МДП-транзистора в пологой области. Видно, что при увеличении Xj для заданного напряжения на стоке ши­ рина / 0тс также увеличивается и выходное сопротивление уменьшается. Выходное сопротивление максимально при Xj=x1{, поэтому для его увеличения желательно умень­ шить глубину стокового р—/i-перехода, используя для этого, например, ионное легирование.

Обычно х ^ х к

 

о

и выраже­

4 мкм, x„«100 А)

ние (4.62) можно упростить

 

 

и 'с - и °*—

/г°« [ 1

2/с

l nX

qNuZvmcXj

 

 

 

(4.63)

Выражая 10тс из (4.40) и подставляя его в (4.63), полу­ чаем следующую неявную формулу для тока стока в по­ логой области характеристики:

U'с и ол

. qNjJJ

2ens0

 

Х[‘+ « я н г г К * - ')]+ “ •(' - тг>

(4.64)

где 1*с — значение тока стока на границе крутой и по­ логой областей. Экспериментальные и смоделированные характеристики МДП-транзисторов с разными длинами

204

каналов приведены на рис. 4.12, 4..13. Как видно, точ­ ность моделирования составляет 5—10%. Модель прове­

рялась

на приборах,

имеющих L= 3—25 мкм;

дгд=

= 900—1500 А;

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

под­

ложки

р=2—5

Ом-см;

jCj= 3,5 мкм; 1>иас = 6-10е

см/с;

£ с = 3-10*В/см; Л'к= 100 А; £«*»= Ю6 В/см; р0=200см2/Вс.

Рис. 4.12. Экспериментальные (ООО) и смоделированные (---------

)

по уравнениям (4.58), (4.64) выходные характеристики

МДП-транзи-

сторов с разной длиной канала при Z= 100 мкм, р=5

Ом-см:

 

a)3,7 мкм; 6) 5,8 мкм; в) L - Ю,8мкм; г) 1=20,9 мкм.

Система определения параметров этой модели такая же, как для модели, определяемой соотношениями (4.47) —(4.52). Отличие состоит только в формуле для вычисления параметра Euyi, которая получается следую­ щим образом. Предварительно по выражению (4.59) вы­ числяются значения 1**с=1с(Ецу1-+00), т. е. без учета зависимости 1лрэ(Еу). Определяя разность между значе­ ниями 1**о и экспериментальными значениями тока /с можно вычислить величину EKyi:

2С/кхт

 

ZzCz

Cz)

Zi — Cz ) -

- C , In

f Zz — C i

(4.65)

 

Zi —Cl

■)}

Две рассмотренные модели МДП-транзистора учи­ тывают., практически все основные эффекты, свойствен­ нее'работе прибора. Они .применимы в широком интер­ вале изменения структурных параметров, управляющих напряжений и температуры. В крутой области характе­ ристик обе модели имеют примерно одинаковую точ­ ность; в пологой области точность второй модели выше

г, кОм

г,кОм

Рис. 4.13. Экспериментальные (ООО) и смоделированные (---------) зависимости дифференциального выходного сопротивления от режи­

ма

при Z [00 мкм:

 

Т-20,9 мкм.

a)

L=3.7 мкм; б) Т-»5,8. мкм; в) £-10,8 мкм; г)

для МДП-транзисторов с

глубокими р—«-переходами

( X j > 2 —3 мкм), а точность

первой

выше при большие

напряжениях на затворе, так как в ней учтено влияние поля затвор—сток на ширину обедненной области от­ сечки.

Перейдем к рассмотрению модели Фромана — Бентчковского и Вадаша. Эта модель также построена отно­ сительно физико-топологических параметров и предна­ значена для машинного проектирования МДП ИС [17]. При описании характеристик прибора в ней учитываются влияние заряда обедненного слоя, зависимость подвиж­ ности от поперечного поля, влияние паразитных сопро­ тивлений истока и стока, конечная проводимость стока в пологой области,-зависимость характеристик от темпе-

206

ратурЫ; т. е. в ней учтены те же физические эффекты, кроме эффекта |Лр?(£у)> что и в первых двух моделях. Учет модуляции подвижности поперечным полем осу­ ществляется в усредненной форме. Поверхностная эф­ фективная подвижность рРо усредняется по длине канала и считается независимой от координаты у:

Р*рэ = Рч» (£кха/<Екср)

(4.66)

где ЕхСр—усредненное некоторым образом поле Ех(у)\ ЕКх2 и CI —коэффициенты аппроксимации; ро— подвижность при слабых лоля,х, Усредненное поле определяется следующим образом:

Я*р= (У. - */. ■+ Uв- 0 Ж ) вдЦ'пхЯ.

(4.67)

Значения коэффициентов Екх2 и Cit определенные экспериментально, равны: Екх2—6 • 104 В/см; С\ — посто­ янная, не зависящая от толщины диэлектрика и концен­ трации примеси в подложке (Ci=0,15 для р-канальных и 0,36 для n-канальных приборов). Подставляя аппрок­ симацию (4.66) в общее выражение для тока (4.30) и интегрируя его, получаем

г _ UpoCnZ

f

^КХ2епХд

с ~

L

[ е д (С /з -^ о + У а - 0 ,5 У с ) J

X VC- JI CRC (Us-Uo + UB- U

- B ocV9o-\-U)dU. (4.68)

ип+‘А

 

 

 

Ток /с из этого неявного выражения вычисляется итера­ ционным методом. Если паразитные сопротивления исто­ ка и стока равны нулю, то (4.68) приводит к явному выражению для тока стока в крутой области характе­ ристик:

t

P'0CnZ

|"

£|С*2еП*Д

1^4/

с

L

[*д(Уз — Уо + У0 +

0,5УС) J

Л

X {(С/, - У 0+ У в -

и„) (Ус -У н ) —

J- (Ос -

Ун)* -

 

з" Вое [(9о +

Ус)3'2 — (fo — Уи)3/2]

(4.69)

В пологой области уменьшение эффективной длины канала учитывается по формулам (4.40), (4.45). Напря­ жение отсечки по стоку, определяющее границу между

крутой и пологой областями, выбирается из условия не­ прерывности производной dIJdU0 при переходе через эту границу. Проведенные экспериментальные исследова­ ния показали хорошую точность модели и однозначную связь ее параметров со структурными параметрами МДП-транзистора. В табл. 4.2 приведены две группы

 

 

 

^ Т а б л и ц а 4.2

 

Параметры прибора

Эксперимен­

Теоретические

 

тальные

ЛГД,

см-1

8-1014

м о и

хД| мкм

0,105

0,110

iVnoB. см"2

5-1011

4,7.10п

Ро. СМ2/В-С

175

170

и

В

3,6

3,65

Zb

мкм

28

26

Lu

мкм

3,8

4

Z2,

мкм

10

12,5

1г.

мкм

46,5

50

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

Рис. 4Л4. Ключевая схема (а), экспериментальные (ООО) и рассчи­ танные (--------- ) с помощью модели Фромана—Бентчковского и Вадаша выходные характеристики активного 77 и нагрузочного Т2 транзисторов (б).

ментально из непосредственны^ электрофизических из­ мерений. Хорошее согласие между теоретическими и экспериментальными значениями параметров позволяет использовать эту модель для проектирования МДП ИС.

208

На рис. 4.14 приведены экспериментальные и смодели­ рованные вольт-амперные характеристики активного 77 и нагрузочного Т2 транзисторов ключевой схемы.

Модель Фромана — Бентчковского и Вадаша имеет ряд недостатков. В уравнениях не учтена зависимость поверхностной эффективной подвижности от продольно­ го поля, что может внести заметную ошибку при моде­ лировании МДП-транзисторов с короткими каналами (L<5 мкм). Зависимость рьрэ(Вх) учтена в усредненной форме (не учтено изменение по длине канала), что снижает точность расчетов при сильной деформации ка-

а

В

Рис. 4.15. Структура

МДП-транзистора (а) и его активная

область {б).

 

нала (при большом изменении потенциала канала по его длине), которая имеет место, если напряжения затвор— исток и затвор—сток значительно отличаются друг от друга. Процедура вычисления тока по уравнениям этой модели содержит два внутренних итерационных цикла: один —для расчета границы между крутой и пологой областями, второй —для вычисления тока с учетом па­ разитных сопротивлений. Это требует больших затрат машинного времени. Аппроксимация зависимости ЦрЭ(£*ср) степенной функцией также требует большого объема машинных вычислений. Перечисленные недостат­ ки снижают эффективность модели и ограничивают ее применение, особенно для МДП-траизисторов с коротки­ ми каналами.

Динамический режим. Активная область МДП-тран- зистора (рис. 4.15) наиболее полно может быть пред­ ставлена эквивалентной схемой с распределенными па­ раметрами (рис. 4.16), соответствующей основным урав­ нениям (4.19), (4.20). Эквивалентная схема содержит

ёмкостные элементы диэлектрика CAdy, обеДненйого слой подложки C00dy и элемент проводимости канала gn/dy, которые определяются следующими выражениями:

/*»

еД £ °

Г»

dQoc

1 /Л ГВос

/ А 7 г\\

 

 

Loc—

 

— 2"Сдт т

4,/и

 

£к =

Сд ||1рэ [(С/'з *Р---босV"?) ~\~

 

 

+ b - g f (и '* - f -

j.

 

Рис. 4.16. Представление активной области МДП-транзистора экви­ валентной схемой с распределенными параметрами.

Распределенная эквивалентная срсема учитывает собст­ венную инерционность прибора и может быть использо­ вана для анализа его высокочастотных и низкочастотных характеристик. Исследование высокочастотных характе­ ристик £18—20] показывает, что инерционность прибора характеризуется временем 1—10 нс; длительность же переходных процессов даже в самых быстродействующих цифровых МДП-транзисторнырс схемах определяется пе­ резарядом паразитных и нагрузочных емкостей схемы и составляет 50—100 нс. Поэтому при расчете подав­ ляющего большинства схем собственной инерционностью МДП-транзистора можно пренебречь и для моделирова­ ния использовать низкочастотные схемы.

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

210

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