Научная статья на тему 'Методы моделирования тепловых повреждений полупроводниковых приборов'

Методы моделирования тепловых повреждений полупроводниковых приборов Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
336
85
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Чумаков Владимир Иванович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Чумаков Владимир Иванович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Simulation of radioelectronic devices thermal failures

The reviewe of thermal failure models of radioelectronic devices is given. Linear heat model of pulse failures is developed.

Текст научной работы на тему «Методы моделирования тепловых повреждений полупроводниковых приборов»

УДК 621.382

МЕТОДЫ МОДЕЛИРОВАНИЯ ТЕПЛОВЫХ ПОВРЕЖДЕНИЙ ПОЛУПРОВОДНИКОВЫХ ПРИБОРОВ

ЧУМАКОВ В.И.

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

1. Введение

Исследования теплового механизма деградаций радиоэлементов и компонентов являются одной из актуальных задач при проектировании РЭА. Это объясняется тем, что тепловые перегрузки являются наиболее тяжелым режимом работы элементов РЭА. Кроме того, деградационные эффекты, связанные с превышением допусков теплового режима, приводят к необратимым повреждениям элементов и являются наиболее частой причиной катастрофических отказов РЭА. Возникновение тепловых перегрузок связано обычно с нарушением электрических режимов функционирования полупроводниковых элементов и компонентов РЭА, наиболее восприимчивых к воздействию внешних факторов. При этом деградационные эффекты могут проявляться в результате весьма кратковременных импульсных электрических перегрузок. Ужесточение требований к обеспечению теплового режима полупроводниковых приборов вызвано миниатюризацией и высокой степенью интеграции элементов (особенно цифровых ИС), переходом к субмикронным технологиям и повышением широкополосности и быстродействия РЭА и средств вычислительной техники (ВТ) [1-5].

Один из аспектов проблемы анализа деградационных эффектов в элементах РЭА при импульсных энергетических нагружениях заключается в ответе на вопрос: какова будет реакция объекта на воздействие токов и напряжений данной интенсивности. Ответ на этот вопрос зависит в первую очередь от теплофизических параметров самого объекта, т.е. от того, какое количество энергии будет рассеяно в структуре объекта в течение длительности воздействия и достаточно ли этой энергии для реализации деградационных эффектов.

В работе рассмотрены различные модели анализа теплового механизма деградации полупроводниковых приборов при импульсных электрических кон-дуктивных нагружениях либо перегрузках, возникающих в результате действия электромагнитного излучения. Приведенные результаты могут быть использованы для исследования деградационных эффектов в полупроводниковых элементах РЭА, возникающих в результате электрических перегрузок, которые вызываются нештатными аварийными режимами РЭА, а также для анализа повреждений, возникающих при облучении РЭА и средств ВТ электромагнитными полями, что характерно для проблем ЭМС и особенно для обострившейся в

настоящее время проблемы электромагнитной стойкости РЭА.

2. Линейные тепловые модели деградационных эффектов в РЭА

При тепловом механизме деградации анализ процесса диссипации энергии в объекте проводится на основе уравнения теплопроводности

£L = v<mvt» + 0, (і)

где T — температура в точке с координатой гв момент времени t; С — линейный дифференциальный оператор; kT, Cp, р — удельная теплопроводность, теплоемкость и плотность полупроводника; qv(t) — объемная плотность мощности источников (стоков) тепла в структуре полупроводника.

Граничные и начальные условия задаются структурой объекта и условиями его эксплуатации.

Моделирование процессов деградации в элементах и компонентах РЭА проводится традиционно в рамках линейной тепловой модели, в которой величины kT, Cp, считаются не зависящими от температуры, что приводит к упрощению анализа.

2.1. Модель W-B

Классическая тепловая модель (модель W-B) деградации дискретных полупроводниковых приборов (диодов и транзисторов) в результате электрических перегрузок приведена в работе [1]. Рассматривается модель плоского p-n перехода, в котором мгновенный плоский источник тепла мощностью Q создает тепловой поток в направлении нормали к плоскости перехода. На основании анализа однородного уравнения теплопроводности для плоского одномерного p-n перехода показано, что зависимость пороговой мощности Рп импульсной электрической перегрузки, при которой возникают деградации в p-n переходе, от длительности импульса электрической перегрузки описывается выражением

Рї/S = B] (Tn-TH) Г0’5 , (2)

где S — площадь перехода; Bj = ^лкт pCp — постоянная величина, определяемая теплофизическими константами материала; TH, Тп — начальная и конечная температура материала соответственно.

Под деградацией можно понимать любое явление, приводящее к изменению характеристик полупроводникового прибора при достижении температуры перехода Тп. Так, в качестве конечной температуры может рассматриваться температура плавления полупроводника, тогда деградационный эффект проявляется в виде проплавления перехода и необратимого отказа полупроводникового прибора [1, 5].

В качестве характеристики восприимчивости полупроводниковых приборов различных типов к тепловым перегрузкам вводится величина

wB =B3(Tn -Th)i (3)

называемая постоянной повреждений (постоянная W-B).

Зависимость (2) выполняется в диапазоне длительностей электрической перегрузки 10-7У2*10-5 с, что подтверждается результатами приведенных в [1] экспериментов.

Физическая интерпретация зависимости (2) основана на предположении о неравномерном распреде-

РИ, 1999, № 2

31

лении плотности тока, протекающего через p-n переход в диапазоне исследуемых длительностей импульса, локализации его в отдельных каналах, что эквивалентно уменьшению эффективной площади перехода, т.е. уменьшению зоны энерговыделения. При этом температура, локальных областей возрастает быстрее, чем происходит теплообмен с прилегающими областями полупроводника, и по достижении температуры плавления Тп возникает тепловое повреждение.

При укорочении длительности и увеличении мощности электрической перегрузки область энерговыделения в объеме полупроводника практически ограничивается объемом токовых каналов, и передачи энергии в прилегающие области не происходит. В этом случае первым слагаемым в правой части (1), которое описывает рассеяние мощности в структуре полупроводникового прибора в результате теплопроводности, можно пренебречь. Тогда зависимость мощности Рп от длительности импульса электрической перегрузки имеет вид:

Рї/S = B2(Tn-TH) t-1 , (4)

где В2 также определяется характеристиками полупроводника.

В режиме длинноимпульсных электрических перегрузок происходит относительно медленный разогрев всего объема полупроводника в результате теплообмена между его отдельными областями, т.е. наблюдается режим, близкий к стационарному. При этом в уравнении (1) можно пренебречь выражением в левой части и получить зависимость для Рп в виде Рї/S = В3(Тп-Тн) t0 , (5)

т.е. мощность повреждений не зависит от длительности импульса. Константа В3 в выражении (5) представляет собой плотность потока мощности, необходимой для разогрева структуры на 10 С. На рис. 1 приведена качественная зависимость пороговой мощности электрической перегрузки, при которой возникают деградации, от длительности импульса, построенная на основании выражений (2), (4) и (5) [2].

Построив качественную зависимость пороговой энергии импульсных электрических перегрузок, приводящей к деградации полупроводниковых приборов, Wn = Pnt, можно видеть, что минимальные значения пороговой энергии располагаются в короткоимпульсной области t< 10-7 c (рис .2). В классической тепловой модели зависимость (2) наблюдается при двух противоположных значениях полярности напряжения импульсной электрической перегрузки. При этом причины возникновения нештатного режима функционирования прибора и электрические

процессы установления тока перегрузки не рассматриваются.

На рис. 3, а приведена типовая структура диффузионного p-n перехода. Области [X1,xp] и [xn,X2] соответствуют квазинейтральным областям полупроводникового кристалла, [xp,xn] — обедненная область полупроводника. Эквивалентную электрическую схему такой структуры можно представить в виде рис.3, б, где Rb — сопротивление квазинейтральной области, Rj — сопротивление области обеднения.

При приложении импульса электрической перегрузки положительной полярности относительно области p переход смещается в прямом направлении и соотношение между величинами сопротивления имеет вид Rb >> Rj, поэтому диссипация энергии происходит в квазинейтральной области.

В случае отрицательной полярности электрической перегрузки относительно области p переход смещается в обратном направлении (запирается) и выполняется неравенство Rb << Rj. Через переход

протекает малый обратный ток, величина которого не зависит от приложенного напряжения и существенно зависит от температуры полупроводника. При этом в отсутствие теплоотвода и теплообмена с окружающей средой происходит относительно медленный разогрев полупроводника до температуры Тп. Максимальная температура достигается в области x=0. Данный режим характерен для диапазона длинноимпульсных перегрузок (рис. 1, 2), в котором пороговая мощность, приводящая к деградациям, не зависит от времени и описывается соотношением (5).

При значительных напряжениях обратной полярности на переходе происходит лавинный пробой вследствие ударной ионизации в обедненной области, сопротивление которой резко падает, что приводит к увеличению тока через переход, и процесс энерговыделения происходит аналогично прямосмещенному p-n переходу.

32

РИ, 1999, № 2

В режиме двуполярной электрической перегрузки происходит детектирование входного тока [6]. Результирующий ток представляет собой последовательность униполярных импульсов, соответствующих временным интервалам, в течение которых переход смещен в прямом направлении.

В общем случае выражение для мгновенной мощности, поглощенной полупроводниковым кристаллом, можно представить как

P(t) = i(t)Uj + i2(t)Rb,, (6)

где Uj — напряжение на обедненной области.

В модели W-B учет температурной зависимости теплопроводности полупроводника заключается во введении в полученные выражения усредненной величины теплопроводности. При этом постоянная повреждений оказывается зависящей только от разности температур (Tn-TH).

2.2. Объемная тепловая модель

Дальнейшее развитие модели W- B для трехмерного случая осуществлено в работах [7, 8]. Прирост температуры в полупроводниковом кристалле с характерными размерами a, b, c описывается с применением функции Грина G0, которая для трехмерного уравнения теплопроводности (1) имеет следующий вид:

Go(r, r', t, т)

1

|4- г)]3 2

exp

(r - г)2 4c(t - г)

. (7)

Здесь r’ — координаты мгновенного точечного источника тепла вида b(t,r), действующего в момент времени t; r, t — соответственно координаты точки и время определения температуры; a=kT/Cpp — коэффициент температуропроводности полупроводника (коэффициент тепловой диффузии).

Расчет критической температуры, при которой возникают деградации, в случае произвольной временной зависимости входной мощности P(t) осуществляется на основании интеграла Дюамеля:

1 t d

Т<0 = т"+c^lP{T)d(r^H('-r)ldr- (8>

Функция H(x) имеет смысл переходной характеристики и представляет собой температурную реакцию полупроводника на действие источника электрической перегрузки, мощность которого описывается единичной функцией Хевисайда 1(x). Таким образом, функции G0 и H связаны соотношением

H(') = JG0(r,v)dv . (9)

0

Для перегрузки в виде прямоугольного импульса мощности длительностью т получаем временную зависимость температуры:

Т(t) = Тн + Po[H(r,t) -H(r,t -r)]. (10)

Длительность электрической перегрузки tf, при которой возникают деградации (например, расплавление полупроводника при достижении температуры Tn), можно определить следующими критериями. В первом случае из (10) имеем

т (tf) = Тп = Тн + PoH (r, tf). (11)

Подставляя (11) в (8), приходим к выражению для нормализованной температуры в произвольный момент времени t:

f (t)=ТТ^-ТТн = J P (г)-

0

d

1

Po(t -Д

dr

(12)

Тп - Тн 0 d(t-г)

Здесь P0(x) представляет временную зависимость пороговой мощности деградаций при действии электрической перегрузки в виде прямоугольного импульса. В частности, используя выражения (2) и (3) модели W-B, получаем

tf d f о =f p (r) iTt—V

0 d(tf - t)

1

P0(tf -0

dr =

2Sw

tf

J P(0

B 0

V

-It

tf ~T

(13)

Подставляя в (13) критическую температуру, определяем критерий W-B оценки времени tf, необходимого для реализации тепловой перегрузки:

tf 1

- J P(г) , = dx = 1.

1

2Sw

B0

Ф

tf ~т

(14)

Другой критерий определения величины tfПOЛу-чается из уравнения

1

1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

df (tf) dt

(15)

При этом необходимо выполнение условия maxf(t)}>1. Подставляя значение tf в (15), получаем выражение для максимального тока электрической перегрузки.

Типичные зависимости нормализованной температуры f(t) для импульсов тока электрической пере-

грузки равной энергии Eи = J i2(t)dt прямоугольной — 1, экспоненциальной с постоянной времени т/ 5 — 2, синусоидальной — 3 и треугольной — 4 форм приведены на рис. 4, 5.

Такие формы импульсов рассматриваются при анализе процессов импульсной теплопроводности [9]. Мгновенная мощность рассчитывалась согласно (6) при Uj = 100IR (рис. 4), что соответствует режиму обратно смещенного р-n перехода и при Uj=I0)Rb (рис. 5), что характеризует переход, смещенный в прямом направлении. Для прямоугольного импульса функция f(t) является монотонной и максимальная температура достигается в конце импульса, в то время как при неоднородном разогреве импульсами сложной формы максимум температуры может достигаться в начале либо в середине импуль-

РИ, 1999, № 2

33

(*> д

Ь: С '

0 L^_l__________________________

О 0,2 0,4 0,6 0,8 t/r

Рис. 5

са. В случае экспоненциального импульса, основная энергия которого сосредоточена вблизи переднего фронта, темп разогрева в начале значительно выше, чем при других формах импульса, однако достижение высоких температур ограничено быстрым спадом импульсной мощности с течением времени. В режиме прямого смещения достигаются более высокие температуры, максимум температуры достигается раньше и, кроме того, темп разогрева и охлаждения (кривые 2, 3, 4) также более высок.

Для области короткоимпульсных электрических перегрузок t < 10-8 с ядро в выражении для нормализованной температуры (12) рассчитывается на основании (4). При этом для обобщенной температуры получаем

ДО = A J р{т)йт , (16)

0

первую очередь потому, что слоистые структуры типа p -i-n диодов являются элементами входных трактов СВЧ аппаратуры и наиболее часто подвергаются воздействию электромагнитных перегрузок. Кроме того, слоистые МДП и МОП структуры характеризуются достаточно низкими критериальными уровнями повреждений. В [10, 11] приведен анализ процесса разогрева многослойных полупроводнико -вых структур при протекании непрерывных и импульсных СВЧ токов. Использована теория развития

Л

х

яуу/у і

Ш:-::

ш

Рис. 7

теплового процесса в слоистой структуре (рис. 7), в которой распространяется поток тепла P(t) в положительном направлении оси x [ 11 ]. Предполагается, что торцевые поверхности образца металлизированы, поэтому температура этих поверхностей постоянна и равна Тн, теплообмен через боковую поверхность образца отсутствует, распространение тепла происходит в направлении оси у. В зависимости от смещения на p-i-n диоде функция источников тепловой мощности предполагается либо равномерной, либо синусоидальной вдоль оси х.

Получены выражения для определения максимального перегрева полупроводниковой структуры и предельной СВЧ мощности, рассеиваемой в приборе, в непрерывном

где A = -

1

T =

ми

4d2

- — постоянная величина. Таким

SB2(Tn -TH)

образом, температура эквивалентна энергии импульса электрической перегрузки. Отметим, что при короткоимпульсных перегрузках из (16) следует, что максимальная температура всегда достигается на заднем фронте импульса, что является результатом адиабатическо -го процесса разогрева (рис. 6 построен для режима обратного смещения p-n перехода).

2.3. Модель повреждений слоистых структур

Моделирование повреждений многослойных полупроводниковых структур имеет важное значение в

T =-------Р

мн „ , с1 мн

1 1 d

8 kp S

и импульсно-периодическом режиме:

(17)

Р1

J ми А

=1 жъ(2ш -1)3 Sd CpkT

m

1,5

1,0

0,5

aP'A'j

0,2 0,4 0,6 0,8

Рис. 6

t/r

1 - exp тг2(1ш -1)2аг

_ d2 _

1 - exp тт2(2ш -1)2ав

d2 _

sin(2m-1)-.(18)

P = P

1 MU 1 MH

Здесь обозначено: d, S — толщина и площадь структуры соответственно; 9 — период повторения импульсов; m — целое число.

Выражение для мощности при этом имеет вид 1 - ехр(-д£)

1 - ехр(-У) , ( )

где §=т/тг — относительная длительность импульса;

Я" kT S T м ак с тт

q — скважность импульсов; Рмн =-------iiatLs. для

8 d

кремния Рмн = 2,67Tмак S / d [Вт].

Величина ттназывается постоянной времени тепловой релаксации и зависит от теплофизических констант материала и характерных размеров d области энерговыделения:

тT = (d / Я-)2Cppl kT . (20)

Таким образом, постоянная времени %т связана с постоянной W-B соотношением

34

РИ, 1999, № 2

R

R

wB = kTJ ятт (TK - TH)/ d . Характерный размер области энерговыделения d при %т = 10-8 c для кремния составляет ~ 3,4 мкм. В цилиндрической системе, когда рассматривается перенос тепла в радиальном направлении в результате действия источника, расположенного параллельно оси цилиндра, величина d соответствует радиусу области энерговыделения.

Переходя к пределу при q ^ ю в выражении для импульсной мощности (19), получаем, что в случае одиночного импульса пороговая мощность дефектообразования описывается выражением

р , = р ----1----

1 и1 1 ми л ,

1 - exp(-Е) ’

график которого приведен на рис. 8.

Раскладывая функцию Ри1 в ряд при §< 1, приходим к выражению Ри1 = Рми/%~1, аналогичному (4) модели W-B, т.е. постоянная времени t тдает оценку верхней границы короткоимпульсной области, в которой справедлива аппроксимация (4). Кроме того, из рис. 8 видно, что при |>1 пороговая мощность практически не зависит от х, т.е. имеет место режим стационарного разогрева, соответствующий зависимости (5) модели W-B (ср. рис. 1).

Подставим tT в выражение для пороговой мощности повреждений (2); получим

( \_0’5

рп = кт№(Tk - Тн )^-J ,

т.е. мощность, необходимая для возникновения повреждения в области W-B, равна pWb = i3ooVS [мкВт], где площадь S выражена в мм2.

2.4. Метод эквивалентной линии

Метод эквивалентной линии основан на сходстве процессов распространения тепла в изотропной среде и распространения волны напряжения в RC длинной линии (рис. 9) [12, 13]. В идеальном случае предполагается, что скорость распространения волны в линии и скорость распространения тепла в среде равны да. Уравнение для напряжения в линии имеет вид

cUJ (x, t) 1 d2U (x, t)

---------_------------“----H Fo

З RC fix2 0

(21)

аналогичный одномерному уравнению теплопроводности, которое можно получить из (1).

Между величинами, входящими в (21) и (1), имеет место следующее соответствие: напряжение U ^ температура Т; погонное сопротивление R ^ тепло-

F0

Ф

C

т—Т

C

ис. 9

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

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

RT =——,

kT S ’

равное величине,

обратной плотности теплового потока, который создается в среде перепадом температуры в 1 К; погон-

ная емкость С ^ тепловая емкость

г - C m С - Crj,

равная теплоемкости образца массой m, имеющего характерный линейный размер l; постоянная времени RC температуропроводность среды а.

2.5. Моделирование повреждений пленочных элементов

Пленочные технологии широко используются при создании полупроводниковых и гибридных ИС, а также элементов и функциональных узлов малогабаритной, в частности, бортовой РЭА. Так, в полупроводниковых ИС средней степени интеграции только общая площадь металлизации (т.е. проводящей металлической пленки) составляет до 30% общей площади кристалла, поэтому анализ процессов в пленочных элементах микроэлектроники широко освещен в литературе (см., например, [12, 14, 15]).

Пробой диэлектрических пленок обусловлен рядом элементарных электрических процессов, которые происходят в пленках в сильных электрических полях.

1. Эмиссия электронов из катода в диэлектрик с вершин микроострий при напряженности электрического поля на микроостриях ~ 107 - 108 МВ/см. Типичная вольт-амперная характеристика (ВАХ) тонкопленочной структуры Al-SiO-Al имеет три участка. На первом (в области сравнительно слабых полей) ВАХ совпадает с законом, рассчитанным по формуле Фаулера- Нордгейма для автоэлектронной эмиссии в результате туннелирования электронов в диэлектрик под действием сильного поля вблизи острия.

2. Размножение электронов в диэлектрике в результате ударной ионизации. При этом рост результирующего тока несколько снижается за счет снижения коэффициента ударной ионизации диэлектрика.

3. Образование и разрушение области отрицательного объемного заряда (ООЗ). Образование ООЗ обусловлено наличием многочисленных энергетических уровней в запрещенной зоне диэлектрика в результате разупорядочения кристаллической структуры, являющихся ловушками для свободных носителей заряда. Такие уровни способны образовывать целые зоны, в которых концентрация ловушек изменяется в пределах 1014-1020 см-3. Электроны, захваченные ловушками, образуют в прикатодном слое диэлектрика толщиной не более 10-5 см область ООЗ, которая оказывает экранирующее действие на микроострия, уменьшая напряженность поля вблизи них. Тем самым ограничивается рост тока на ВАХ тонкопленочной структуры (участок 2).

РИ, 1999, № 2

35

Дальнейшее увеличение напряженности поля приводит к разрушению ООЗ вблизи острий и образованию ООЗ в более глубоких слоях диэлектрика. При этом экранирующее действие ООЗ снижается, что приводит к более сильному росту тока на ВАХ структуры (участок 3) [13].

Приведенные рассуждения отражают сценарий образования источника начального тока, протекание которого через диэлектрик приводит к интенсивно -му энерговыделению и возникновению тепловой деградации.

По электрической проводимости пленки можно разделить на 3 группы:

— проводящие (все элементы металлизации, межсоединения, контактные площадки и др.);

— резистивные (интегральные резисторы);

— изолирующие (защитные покрытия, изолирующие слои, подложки и др.).

Деградационные эффекты, вызванные тепловыми процессами в первых двух группах, обусловлены превышением предельных допустимых величин плот-[ і noe a ea Jnp, либо мощности рассеяния Рдоп, которые связаны соотношением

Рдоп~ Jпр°2r - P0SR ,

где a — сечение пленки; R — сопротивление пленочного элемента; Р0 — удельная мощность рассеяния материала пленки; SR — площадь поверхности пленочного элемента.

Если повреждение пленочного элемента вызвано электрической перегрузкой в результате электромагнитного воздействия [5], то критическая напряженность электрического поля, приводящего к возникновению деградационного эффекта, связана с величинами Рдоп и Jnp соотношением

Eкр _ JпрaRKпр ,

где Кпр—коэффициент преобразования электрического поля, который зависит от геометрии элемента и определяется, исходя из его антенной модели.

В [12] на основании анализа уравнения теплопроводности построена модель процесса неоднородного разогрева пленок SiO2. Показано, что локализация тепла в узком канале протекающего тока приводит к возникновению тепловой неустойчивости, аналогичной режиму с обострением [16, 17]. При этом вольт-амперные характеристики канала пробоя имеют S-образную форму с участком отрицательной проводимости. Получены выражения и оценки величин времени задержки пробоя:

f—a exp(a/T) / o'

tD ~----~ -----:----Ei\

°oEl L a/T

T

где a = (є-pEb)/2k ; e, p — параметры материала пленки; Eb—напряженность приложенного электрического поля; s о — удельная электрическая проводимость материала; Ei(x)—интегральная показательная функция; к — постоянная Больцмана.

Оценка времени переключения, т.е. времени перехода из непроводящего состояния в состояние с высокой проводимостью составляет t^=10-9 c.

2.6. Моделирование процессов деградации сложных структур

Применение тепловой модели позволяет анализировать деградационные процессы в сложных структурах [18]. Рассмотрим структуру металлический про-

водник сечением S — контактная площадка металлизации толщиной к] на поверхности диэлектрика (SiO2), имеющего толщину к2. Через структуру протекает импульс тока i(t) длительностью т (рис.10). Отметим, что такие структуры в ИС малой и средней степени интеграции располагаются по периферии кристаллов и занимают значительную часть их повер -хности, т.е. располагаются именно в тех областях, в которых деградационные эффекты проявляются наиболее часто [5].

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

Т

w п =J i 2(t)R п dt ,

0

1 l

где R п =---— сопротивление проводника длиной

U п S

l, удельная проводимость которого равна s„.

Считая, что эта энергия затрачивается на разогрев проводника до температуры плавления Тгш от начальной температуры То, получим выражение для квадрата действующего значения тока:

12 =

C пРп AT nS 2^

(22)

Здесь АТп = Тпл — То- — перегрев проводника; Сп и р„—теплоемкость и плотность материала проводника соответственно.

Рассуждая аналогично относительно энергии, отводимой в диэлектрик толщиной к2 из контактной площадки размерами ab за счет теплопроводности, получим следующее выражение:

12 . kTd ^ь2„ , (23)

где кТд — теплопроводность диэлектрика; ДТд — разность температур на границах слоя диэлектрика.

Из сравнения (14) и (15) находим критическую длительность импульса, при которой происходит деградационный эффект:

Здесь Рп

ткр ~ тТп

kTп ЕТп °п Sh2 kTд АТд &к b\ '

(24)

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

tt п = —п^п — постоянная времени тепловой релак-

kT п

сации проводника. Если длительность импульса т>ткр, то происходит расплавление проводника, если г<хкр, то вначале происходит расплавление площад-

36

РИ, 1999, № 2

ки. При одинаковых материалах проводника и кон-татктной площадки получаем из (16)

_ _ _ kTп Sh2

тк р тТ п , 2

kTd Ь2hl '

3. Заключение

Проведенный анализ теплового механизма повреждений полупроводниковых приборов показывает, что линейная тепловая модель достаточно точно описывает временные зависимости температуры и мощности повреждений. Критериальными параметрами являются постоянная W-B и время тепловой релаксации полупроводникового кристалла хт,, которые зависят от макроскопических характеристик полупроводника (теплоемкость, теплопроводность, плотность). При укорочении длительности воздействия характер временной зависимости мощности Pn(t) определяется уменьшением эффективных размеров области энерговыделения, что приводит к уменьшению пороговой энергии повреждений. Показано, что при одинаковых уровнях энергии импульсные перегрузки являются более опасными для функционирования РЭА, чем квазистационарные.

Исходными причинами, которые приводят к возникновению тепловых перегрузок, как правило, являются быстропротекающие процессы электрического пробоя, которые не рассматриваются в классических моделях, что ограничивает использование тепловых моделей в области малых времен. В ряде случаев деградационные эффекты в полупроводниковых элементах могут возникать на стадии развития электрического процесса, еще до наступления тепловых перегрузок. Подобные эффекты рассмотрены в [5] при анализе воздействий короткоимпульсных электромагнитных излучений на изделия электронной техники. Показано, что подпороговые воздействия приводят к отклонению параметров полупроводниковых приборов от номинальнах значений и интерпретируются как степень деградации II.

Физическая картина процессов дефектообразования в полупроводниковых приборах при воздействии импульсов перегрузки ультракороткой длительности может быть выяснена в рамках теории электрических процессов в полупроводниках, характерные времена которых значительно меньше времени тепловой релаксации

Дальнейшее уточнение сценария теплового процесса дефектообразования заключается в анализе нелинейного уравнения теплопроводности, аналогичного (1), с учетом температурных зависимостей тепло- и электрофизических параметров полупроводника [16, 19].

Литература. 1. Wunsch D.S.,, Bell R.R. Determination of threshold failure level of semiconductors diodes and transistors due to pulse voltage // IEEE Trans. on Nuclear

Sci., 1968. Vol.15, N6. Р.244-259. 2. R.J. Antinone. How to prevent circuit zapping // IEEE Spectrum. 1987, Vol.4, N24. Р.34-38. 3. A comparison of DC and RF pulse susceptibilities of UHF transistors // IEEE Trans. on Electromagnetic Compatibility, 1977. Vol.19, N9. Р.49-56. 4. Мырова Л.О., Чепиженко А. З. Обеспечение стойкости аппаратуры связи к ионизирующим и электромагнитным излучениям. М.: Радио и связь, 1988. 296с. 5. Блудов С.Б., Гадецкий Н.П., Магда И.И. и др. Еенерирование мощных СВЧ-импульсов ультракороткой длительности и их воздействие на изделия электронной техники // Физика плазмы, 1994. Т.20. Вып.7,8. С.712-717. 6 Бригидин А.М., Листопад Н.И., Титович Н.А., Ясюля Г.И. Исследование восприимчивости полупроводниковых приборов и интегральных схем к воздействию ВЧ и СВЧ помех // Радиотехника и электроника. Минск.: Вышэйшая школа. 1990. Вып.19. С. 115-119. 7. Dwyer V.M., Franklin A.J., Campbell D.S. Thermal failure in semiconductor devices // Solid State Electronics. 1990. Vol.33. Р.553-560. 8. Dwyer V.M., Franklin A.J., Campbell D.S. Electromagnetic discharge thermal failure in semiconductor devices // IEEE Trans. on Electr. Dev. 1990. Vol.37, N11. Р.2381-2387. 9. Камья Ф.М. Импульсная теория теплопроводности: Пер. с франц./ Под ред. А.В. Лыкова. М.: Энергия. 1972. 272с. 10. СВЧ устройства на полупроводниковых диодах. Проектирование и расчет / Под ред. И.В.Мальского и Б.В. Сестрорец-кого. М.: Сов. радио, 1969. 580с. 11. Mortenson K.K. Analysis of the temperature rise in PIN diodes caused by microwave pulse dissipation // IEEE Trans. on Electr. Dev., 1966. Vol.13, N3. Р.305-313. 12. Shousha A.H.M. Negative differential conductivity due to electrothermal instabilities in thin amorphous films // J.Appl.Phys. 1971. Vol.42, N12. Р.5131-5236. 13. Lange T.J., Hjellen G.A. A comparison oftest and model-predicted RF-pulse susceptibilities of UHF transistors // IEEE Trans. on Electromagnеtiс Compatibility. 1978. Vol.20, N.4. Р.513-514. 14. ВоробьевГ.А., МухачевВ.А. Пробой тонких диэлектрических пленек. М.:Сов. радио. 1977. 72с. 15. WongK.L. Effects of electromagnetic interference for electromagnetic pulses incident on microstrip circuits / IEE Proceedings. 1990. Vol.137, N1. Р.75-77. 16. Вирченко Ю.П., Водяницкий А.А., Ковтун Г.П. Локализация тепла и становление структуры теплового пробоя. Обзор: Харьков, ХФТИ, 1992. 32с. (препринт). 17. Галактионов В.А., Курдюмов С.П., Посашков С.А., Самарский А.А. Квазилинейное параболическое уравнение со сложным спектром неограниченных автомодельных решений. В кн.: Математическое моделирование. Процессы в нелинейных средах. М.: Наука, 1986. С.142-182. 18. Smith J.S. Electrical overstress failure analysis in microcircuits // 16 Annual Proceedings Reliability Physics. IEEE Proceedings Reliability Physics. 1978. Р.41-46. 19. МаллерР., Кеймингс Т. Элементы интегральных схем: Пер. с англ. М.: Мир, 1989. 630 с.

Поступила в редколлегию 05.05.99 Рецензент: д-р техн. наук Шокало В.М.

Чумаков Владимир Иванович, канд. физ.-мат. наук, доцент кафедры ОРТ ХТУРЭ. Научные интересы: генерации мощных электромагнитных излучений и их применение. Адрес: Украина, 61726, Харьков, пр. Ленина, 14, тел.: (0572) 47-29-55.

РИ, 1999, № 2

37

i Надоели баннеры? Вы всегда можете отключить рекламу.