книги / Математические модели элементов интегральной электроники
..pdfтывающцх изменение геометрических размеров на технологических операциях [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В- 8*Ю11 см -а, ------ 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С/кхт |
|
Zz— Cz |
— 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