- •1. Расчет коэффициентов ускорения для основных механизмов отказов.
- •1.1 Электромиграция.
- •1.2 Коррозия.
- •1.3 Зависящий от времени пробой диэлектрика (tddb)
- •2. Моделирование отказов в мэмс переключателях, обусловленных накоплением заряда
- •2.1. Моделирование отказов, обусловленных инжекций и захватом носителей заряда в ловушечный центры в диэлектрике.
- •Туннельная инжекция.
- •Моделирование деградации порогового напряжения мдп–транзистора вследствие захвата горячих носителей.
- •Порядок работы с программой
Туннельная инжекция.
При приложении к затвору МДП-структуры больших величин напряжений в подзатворном диэлектрике возникает высокое значение напряжённости электрического поля. В этом случае величина потенциального барьера оказывается сравнима с дебройлеровской длиной волны электрона и возможна туннельная инжекция cвободных носителей из полупроводника в зону проводимости диэлектрика. Как и ранее выделим в фазовом пространстве {k,t} элемент объёма dK=dKx, dKy, dKz вблизи значений волнового числа k, тогда число состояний dz в этом объёме будет:
(9)
число электронов dn в этом объёме
(10)
Будем считать, что туннельный ток обусловлен электронами, имеющими компоненту скорости Vx обладающими значениями энергии E и имеющими отличную от нуля при этом значении энергии величину коэффициента прозрачности D(Е)
Тогда
(11)
Подставляя (10) в (11) и пользуясь тем, что D(Е) зависит только от компоненты Vx получим:
(12)
Учитывая:
Получаем:
(13)
К
30
.
(14)
Подставляя это в выражение для тока и интегрируя получаем:
(15)
Здесь у выражает понижение потенциального барьера за счет эффекта Шоттки:
t(y) , (y) – функции параметра у слабо от него зависящие и табулированные в [z].
Выражение (15) часто записывают в виде:
(16)
г
29
Для контакта Si–SiO2 высота туннельного барьера Ф0=3,2 эВ, ,m*=0.42me , F0=2,3108 В/см , А=10–7 А/В2.
31
Моделирование деградации порогового напряжения мдп–транзистора вследствие захвата горячих носителей.
Физическая и математическая модель деградации порогового напряжения.
При протекании постоянного тока через МДП структуру происходит накопление заряда на ловушках в диэлектрике, а также генерация новых ловушек. В связи с этим возникает необходимость определения пределов изменения порогового напряжения при известных параметрах центров захвата: количества типов ловушек, эффективной плотности, сечения захвата.
Для случая постоянного тока через структуру действительна кинетическая модель, в соответствии с которой захваченный в диэлектрике заряд равен:
Q =
q
0i{1–exp[–i
(jt/q)]}.
(1)
Захват заряда приводит к сдвигу напряжения плоских зон в соответствии с формулой:
VFB= –(qtox/0) 0i{1–exp[–i (jt/q)]}, (2)
где N–количество типов ловушек; n0i – эффективная плотность данного типа ловушек (1/см2); i – сечение захвата данного типа ловушек; (см2); j – плотность тока (А); t – время инжекции (сек); dox – толщина диэлектрика (см).
32
Программа моделирования позволяет рассчитать изменения напряжения плоских зон в соответствии с формулой (2) и представляет результаты в виде графика на экране монитора.
Порядок работы с программой
Включить компьютер и загрузить программу «Проводник». Запустить программу расчета МОП.exe двойным нажатием левой кнопки мыши.
Ввести данные и нажать «Enter». При этом появится значение напряжения плоских зон VFB, а в правой части экрана будет представлен график VFB(t).
Изменяя входные данные, исследовать зависимость VFB(t) для нескольких типов ловушек, при различных параметрах ловушек и уровнях тока инжекции (см. таблицу вариантов).
Результаты
моделирования представить в виде графика
VFB(t)=f(t)
при j=const.
По кинетическим кривым для случая
ловушек одного типа можнорассчитать
сечение захвата, предварительно
представив результат в координатах
от t.
N 1 2 3 4 5 6 7 8 9 10 11 12
J
t
N0
I
ni0
tox, 10 20 30 40 50 60 70 80 85 90 95 100
ox=3,8, 0=8,8610–14 Ф/см, q=1,610–19 Кл.
Выводы
Обобщая результаты ускоренных испытаний ИС на долговечность, можно сделать следующие выводы:
1. Уравнение Аррениуса является базовым для описания зависимости долговечности ИС от температуры. Отказы, возникающие в условиях повышенной нагрузки, обусловлены механизмами отказов, которые, как правило, не подчиняются зависимости Аррениуса и не проявляются в ИС, спроектированных с учетом надежности.
2 Ускоренные испытания показывают, что долговечность ИС носит бимодальный характер. Зависимость долговечности ИС от температуры хорошо описывается логарифмически нормальным распределением. Температура в сочетании с повышенным напряжением является универсальным воздействующим фактором, позволяющим выявлять механизмы отказов элементов конструкции БИС (кристалл - оксид - металл).
3. В производственных условиях ведущие фирмы-изготовители интенсивность отказов для периода нормальной работы ИС на «ваннообразной» кривой интенсивности отказов оценивают по результатам квалификационных испытаний с использованием распределения хи-квадрат.
4. Ускоренные испытания позволяют с определенной достоверностью за относительно короткий срок испытаний оценить долговечность ИС, выражающуюся сотнями тысяч лет, т. е. позволяют изучить период «старения» ИС через взаимосвязь механизмов отказов с периодом времени их проявления.
