МОДЕЛИРОВАНИЕ РАДИАЦИОННЫХ ВОЗДЕЙСТВИЙ НА МИКРОЭЛЕКТРОННЫЕ ИЗДЕЛИЯ
УДК 621.382:621.396.6
МОДЕЛЬ ФОРМИРОВАНИЯ ТОКОВ УТЕЧКИ ДИЭЛЕКТРИКОВ МОП-СТРУКТУР ПРИ ВОЗДЕЙСТВИИ ТЗЧ © 2015 г. А. В. Согоян, В. А. Полунин
Национальный исследовательский ядерный университет МИФИ Институт экстремальной прикладной электроники E-mail: vapln@spels.ru Поступила в редакцию 11.08.2014 г.
Деградация подзатворного диэлектрика МОП-структур является одним из основных факторов, определяющих долговечность работы современных КМОП СБИС. Экспериментальные данные свидетельствуют о том, что воздействие тяжелых заряженных частиц приводит к снижению характеристик надежности тонких (менее 10 нм) диэлектриков, возникновению токов утечки. При этом эффекты, обусловленные воздействием тяжелых заряженных частиц, оказываются "латентными". В работе предложена модель формирования токов утечки диэлектрика, обусловленных воздействием тяжелых заряженных частиц. Модель основана на решении системы кинетических уравнений и позволяет преодолеть ряд ограничений известного перколяционного подхода.
Б01: 10.7868/8054412691501010Х
Деградация подзатворного диэлектрика МОП-структур является одним из основных факторов, определяющих долговечность работы современных КМОП СБИС. Она может вызываться электрическими нагрузками (инжекцией горячих носителей из подложки и высокополевой инжек-цией) и действием ионизирующего излучения (ИИ), в частности — тяжелых заряженных частиц. Влияние воздействия тяжелых заряженных частиц на характеристики тонкого окисла неоднократно наблюдались экспериментально [1—4]. Так, в работе [2] установлено существенное влияние предварительного воздействия тяжелых заряженных частиц на электрическую прочность окисла. Пленки окисла толщиной 2.2—3.3 нм, облучавшиеся тяжелыми заряженными частицами, демонстрировали рост токов утечки и снижение пробивных напряжений. Эффект усиливался с ростом ЛПЭ тяжелых заряженных частиц и уменьшением толщины диэлектрика. Облучение диэлектрика толщиной 3.3 нм ионами Аи приводило к снижению критического заряда на два порядка. В работе [3] наблюдалось снижение характеристик надежности окисла (уменьшение времени до начала пробоя при электрическом стрессе 1вв) после воздействия ионов 81, Вг, I и Аи. При этом снижение 1во отмечалось также в приборах, в которых не наблюдалось существенной утечки сразу после облучения. Данный результат соответствует очевидному предположению о том, что причиной деградации окислов являются струк-
турные дефекты, генерируемые при прохождении тяжелых заряженных частиц.
В работе [4] экспериментально исследовалось возникновение токов утечки в окисле МОП-структур с плавающим затвором. Результаты эксперимента свидетельствовали о возникновении перманентных каналов утечки подзатворного окисла в результате воздействия ионов. Величина утечки коррелировала с ЛПЭ ионов.
Таким образом, имеющиеся экспериментальные данные свидетельствуют о том, что воздействие тяжелых заряженных частиц приводит к снижению характеристик надежности тонких (менее 10 нм) диэлектриков, возникновению токов утечки. При этом эффекты, обусловленные воздействием тяжелых заряженных частиц, оказываются "латентными", проявляясь в виде отказов приборов спустя некоторое время после воздействия.
В целом, дефектообразование в ПП и диэлектриках при прохождении тяжелых заряженных частиц обусловлено действием ПВА и термической модификацией структуры вблизи оси трека [5, 6].
Образование отдельных дефектов структуры первичными выбитыми атомами (ПВА) можно проанализировать с помощью численной модели [7]. На рис. 1 показаны расчетные распределения числа дефектов в пленке окисла толщиной 10 нм (показана вероятность генерации числа дефектов, превосходящего заданную величину в расчете на 1
5
65
Число дефектов
Рис. 1. Распределения числа вакансионных дефектов в слое SiO2 толщиной 10 нм (в расчете на 1 первичный ион).
0 20 40 60 80 100
г, нм
Рис. 2. Радиальные распределения плотности вакансионных дефектов в слое SiO2 (в расчете на 1 первичный ион).
первичный ион). Распределения нормированы с учетом событий, соответствующих отсутствию генерации дефектов. На рис. 2 показаны радиальные распределения плотности дефектов, генерированных в слое 8Ю2 при прохождении ионов. Представленные результаты расчета показывают, что средняя плотность структурных повреждений, обусловленных ПВА, быстро спадает при удалении от оси трека. Вместе с тем, как показывает анализ, при реализации событий с большим количеством повреждений имеет место образование кластеров дефектов в соответствии с традиционными представлениями.
Вторым механизмом образования повреждений при воздействии тяжелых заряженных частиц является термическая модификация структуры материала [5, 6]. Ионизационные потери энергии ОЯЧ передаются решетке посредством электрон-фононного взаимодействия, что может приводить к повышению температуры вещества выше точки испарения вблизи оси трека. Для анализа этого процесса принципиальное значение имеет пространственно-временной характер передачи энергии от излучения решетке.
На рис. 3 показана временная зависимость темпа передачи энергии от электронов решетке. Расчет проводился для случая прохождения иона
МОДЕЛЬ ФОРМИРОВАНИЯ ТОКОВ УТЕЧКИ ДИЭЛЕКТРИКОВ
67
о
В
т
300 X 100 250х 100 200 х 10° 150х 100 100х 100 50 х 100 00
10
10
100 101 фс
102
103
Рис. 3. Кинетика передачи энергии от электронов решетке. Расчет проводился для случая прохождения иона Ag с энергией 275 МэВ (57.3 МэВ см2/мг).
Ag с энергией 275 МэВ (57.3 МэВ см2/мг) через слой 8Ю2 толщиной 20 нм в рамках системы моделей, описанной в [8]. Представленные зависимости соответствуют пространственному усреднению по области трека. Результаты свидетельствуют, что основной вклад в передачу энергии вносят процессы рекомбинации электронов и дырок за времена менее 100 фс.
Анализ процессов термической модификации проводился средствами молекулярной динамики. Межатомное взаимодействие описывалось модифицированной моделью Тегюй" [9]. Моделировалось поведение аморфного 8Ю2 в форме параллелепипеда с размерами 10 х 10 х 5 нм. Боковой слой толщиной 1 нм термостатировался методом Берендсена [6]. При расчете применялись периодические граничные условия. Передача энергии атомам осуществлялась в соответствии с результатами предшествующего расчета (рис. 3). Поскольку основной вклад в передачу энергии обусловлен процессами рекомбинации электронов и практически иммобилизованных дырок, пространственное распределение передаваемой энергии определялось расположением дырок после прохождения частицы.
Результаты проведенного расчета в целом аналогичны представленным в работе [6]. Характерный поперечный размер области структурной модификации материала составил около 8 нм при ЛПЭ 57.3 МэВ см2/мг.
Оценку влияния характера пространственного распределения структурных дефектов на величину тока утечки необходимо проводить с учетом вида модели переноса заряда через диэлектрик с
участием ловушек. В настоящее время предложен ряд моделей [10—18], учитывающих вклад множественных перколяционных каналов в пленке окисла. Значительная часть предложенных моделей описывает туннельную проводимость окисла с участием одной или нескольких ловушек. С другой стороны, значительная часть моделей, описывающих туннелирование с участием многих дефектов окисла, опирается на подход, изложенный в работе [17]. Предлагаемый в [11] подход основан на предположении, что ток утечки пропорционален минимальной вероятности перехода с дефекта на дефект при перколяционном прохождении электронов через диэлектрик. В результате может быть рассчитан ток в обособленных перколяционных каналах. В реальных условиях, однако, каналы перколяции не всегда могут считаться обособленными и протекание тока необходимо анализировать с учетом их взаимного влияния. Кроме того, заселенность// уровня /-го дефекта в работе [11] принималась термодинамически равновесной. Данный подход представляется некорректным, поскольку при наличии электрического поля величина // в диэлектрике зависит от локализации дефекта, варьируется на несколько десятков порядков величины и определяется процессами перехода с ближайших дефектов.
В данной работе проводилось моделирование токов утечки в приборах энергонезависимой памяти при воздействии тяжелых заряженных частиц [11]. Для решения данной задачи рассматривалась система дифференциальных уравнений для вероятности заполнения дефектов:
Рис. 4. Схема электронных переходов.
Ь-МЫ-Ж-М-Ш.................^ (!)
] =1 ] =1
где п — число дефектов, созданных ионом в окис- /-м и у-м дефектами. В системе (1) катод и анод ле, /() — вероятность заполнения /-го дефекта имеют номера п + 1 и п + 2. Темпы переходов Щ электроном, Щ — темп перехода электрона между определялись в соответствии с моделями [19-22]:
Щ « РиРт(ЕиЕ]), Рип = ехр \--:-ЦЕ) !>,
Й
-.3/2
еТ - [[ - дЕ^ -х 2 дЕ(1, - z)¡Ei 3'
где Р1ип — вероятность туннельного перехода с центра / на центр у, Е, Еу — глубины залегания уровней центров относительно дна зоны проводимости, Гу — расстояние между центрами, — координаты центров, РТ — фактор, учитывающий
вероятность многофононных переходов между центрами / и у [18—20].
На рис. 4 показана схема переходов между уровнями дефектов. Для темпа переходов между дефектами 1 и 2:
В соответствии с моделью [18—20],
Я,2 -
42
|Рп(Е, й, Е)№ (Е - Ел + Бя) V (2 + Ей - Е - Бм>)Е
т=\ —
да , .
У Рп(Еа + я(т - Б), й, Е№т(т)Фт (2 - Ел + Ей - ж),
Е+Рй 1 \ я Ы (2)
Я ]
Жт(ж) - ехр( + Жя)1(ж,с), V(Е) - 1 Wж (|, 2кТ) я ш
с1 - Б/Ш (я/кТ), с - Б/(((2кТ)),
МОДЕЛЬ ФОРМИРОВАНИЯ ТОКОВ УТЕЧКИ ДИЭЛЕКТРИКОВ
69
101 100 10-1 10-2
о
10-3 10-4 10-5
3 4 5 6 7 8
Ут, В
Рис. 5. Распределения пороговых напряжений при воздействии ионов йода; фигуры — эксперимент [11]; линии — расчет.
: 4 часа
, 481 час
! -Модель /о jfa
- ^^ ъУ J-C0° - У JF"
; у ¿f
vf сР . J o°J 1 1 1 1
где w ~ 0.06 эВ, S — фактор Хуанга—Риса.
Для глубоких (>1 эВ) уровней выражение (2) можно упростить:
Ru » Pu(Et, d, F)—exp(-20! + I(a,2c). (3)
w \ 2kTl
Для темпа переходов с ловушки на катод
I w
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.