Поверхность. Рентгеновские, синхротронные и нейтронные исследования, 2021, № 3, стр. 31-38

Транспорт вторичных ионов – продуктов неупругого взаимодействия протонов с материалами интегральной схемы

Н. В. Новиков a*, Н. Г. Чеченин a, А. А. Широкова a

a Московский государственный университет им. М.В. Ломоносова, НИИЯФ им. Д.В. Скобельцына
119991 Москва, Россия

* E-mail: nvnovikov65@mail.ru

Поступила в редакцию 04.06.2020
После доработки 27.08.2020
Принята к публикации 30.08.2020

Полный текст (PDF)

Аннотация

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

Ключевые слова: сбои электроники вследствие воздействия отдельного иона, линейная передача энергии, тормозная способность материала, количество электронно-дырочных пар.

ВВЕДЕНИЕ

В составе космических лучей протоны составляют более 90% и являются основным источником сбоев электроники в околоземном пространстве [1]. Одной из причин таких сбоев в элементах интегральных схем являются ионизационные столкновения заряженных продуктов неупругого взаимодействия с атомами кремния в чувствительном объеме (SV). Обычно такие события происходят в ячейках памяти или в статических триггерах при попадании в них ионов. Возникший при этом импульс тока переводит ячейку или триггер в противоположное состояние [2, 3]. Если количество электронно-дырочных пар в чувствительном объеме и накопленный избыточный заряд Qdep превышает критическое значение Qc [4], то прогнозируется изменение вольтамперной характеристики этого узла интегральной схемы, которое может привести к пробою диэлектрического слоя и к переключению логического устройства в виде ложного срабатывания [5, 6]. Величина критического заряда Qc зависит от конструкции интегральной схемы – толщины чувствительного слоя dSV, количества и взаимного расположения чувствительных областей, состава и расположения элементов вблизи чувствительной области, их удаленности от трека первичного иона.

Ионизационные столкновения одного первичного протона создают небольшое количество электронно-дырочных пар, и, как правило, накопленный (осажденный) избыточный заряд в этом случае меньше критического значения: Qdep < Qc. Это связано с небольшими неупругими потерями энергии протонов и ионов отдачи Si+ вследствие неупругих столкновений с атомами. Однако при увеличении заряда первичных ионов и их потока могут произойти накопление электронно-дырочных пар, увеличение их концентрации и генерация локального заряда в чувствительной области. Ключевым приближением в этом подходе является предположение о том, что вероятность сбоя интегральной схемы не зависит от вида и спектрального состава первичного излучения, а определяется линейной передачей энергии первичной частицы [7].

Кроме того, накопление в чувствительном объеме избыточного заряда Qdep> Qc возможно и вследствие образования заряженных продуктов неупругого взаимодействия первичного иона с ядром в чувствительном объеме или в материалах, близко расположенных к нему. Неупругие потери энергии заряженных продуктов ядерной реакции могут быть значительно больше, чем первичной частицы и ядер отдачи Si+. Если количество электронно-дырочных пар, образующихся в результате взаимодействия вторичных ионов, достаточно для выполнения условия Qdep> Qc, то может произойти переключение логического устройства вследствие воздействия отдельного иона с энергией E0. Вероятность таких сбоев пропорциональна количеству первичных частиц, попавших в чувствительную область N0, и потерям энергии заряженных фрагментов ядерной реакции. Соотношение вкладов разных источников в накопление Qdep зависит от энергии первичного иона E0, конструкции интегральной схемы (количества и взаимного расположения чувствительных областей), состава материалов элементов вблизи чувствительной области и их удаленности от нее.

Цель настоящей работы – оценка вклада вторичных ионов в материалах, близко расположенных к чувствительному объему, в вероятность сбоя в работе интегральной схемы под действием быстрых протонов. Исследовано влияние материала и расстояния до чувствительной области на накопление избыточного заряда.

СБОЙ ЭЛЕКТРОНИКИ ВСЛЕДСТВИЕ ВОЗДЕЙСТВИЯ ОТДЕЛЬНОГО ИОНА

Оценка воздействия одного иона на накопление избыточного заряда Qdep связана с расчетами количества электронно-дырочных (eh) пар, которые образуются в чувствительном объеме. Эта величина пропорциональна длине трека первичного иона или эффективной хорде d ' [4] внутри чувствительной области толщиной dSV (dSVd '). Транспорт одного иона в среде может привести к накоплению в чувствительном объеме избыточного заряда:

(1)
${{Q}_{{{\text{dep}}}}} = {{q}_{e}}{{F}_{{eh}}}d{\kern 1pt} '~\, \approx {{{{q}_{e}}{{E}_{{de}}}} \mathord{\left/ {\vphantom {{{{q}_{e}}{{E}_{{de}}}} {{{\varepsilon }_{{{\text{Si}}}}},}}} \right. \kern-0em} {{{\varepsilon }_{{{\text{Si}}}}},}}$
где Feh – количество электронно-дырочных пар на единице длины трека, Edep – поглощенная мишенью энергия, εSi = 3.6 эВ – энергия связи электрона в атоме кремния, qe – заряд электрона. Поглощенная энергия Edep пропорциональна длине хорды d ' и для иона c энергией E выражается через его неупругие потери энергии Se(E), [эВ/A] или [кэВ/мкм]:
(2)
${{E}_{{{\text{dep}}}}} = {{S}_{e}}(E)d{\kern 1pt} '.$
В задаче описания сбоев электроники принято использовать величину энергии, переданной ионом материалу мишени LET, которая измеряется в МэВ · см2/мг и получила название линейной передачи энергии (LET). В случае мишени из кремния величина LET связана с Se(E) коэффициентом [8]:
(3)
$LET(E) = 0.043087{{S}_{e}}(E),$
и отличается от тормозной способности материала, в которую кроме электронных (неупругих) потерь энергии Se входят также ядерные (упругие) потери энергии SN. Максимальное количество электронно-дырочных пар на единице длины трека одного иона вычисляется из (1) и (2):
(4)
${{F}_{{eh}}} \approx {{{{S}_{e}}(E)} \mathord{\left/ {\vphantom {{{{S}_{e}}(E)} {{{\varepsilon }_{{{\text{Si}}}}}.}}} \right. \kern-0em} {{{\varepsilon }_{{{\text{Si}}}}}.}}~$
Неупругое столкновение первичного иона с энергией E0 и атома мишени может привести к появлению либо одной электронно-дырочной пары с образованием электрона и иона остатка Si+, либо к взаимодействию с ядром атома и рождению нескольких вторичных частиц. Вероятности этих двух процессов отличаются на несколько порядков. Количество вторичных ионов, образующихся в результате одного неупругого взаимодействия с ядром, определяется суммой:
(5)
${{N}_{{ZA}}}({{E}_{0}}) = \sum\limits_j {{{N}_{j}}({{E}_{j}},{{E}_{0}})} = \sum\limits_i {{{N}_{i}}(LE{{T}_{i}},{{E}_{0}})} ,$
где ${{N}_{j}}({{E}_{j}},{{E}_{0}})$ и ${{N}_{i}}(LE{{T}_{i}},{{E}_{0}})$ – распределения количества вторичных ионов по их энергии Ej и линейной передаче энергии LETi. Отметим, что в сумму (5) входят также вторичные ионы, образующиеся в ходе упругого рассеяния первичного иона на ядре. Энергия таких ионов в столкновениях Si(p, p')Si на порядки превышает энергию ионов отдачи в столкновениях с атомом. Прохождение вторичных ионов через среду приводит к образованию на единице длины ${{F}_{{eh}}}$ электронно-дырочных пар:
(6)
$\begin{gathered} {{F}_{{eh}}}({{E}_{0}}) \approx \frac{1}{{{{\varepsilon }_{{{\text{Si}}}}}}}\sum\limits_j {{{N}_{j}}({{E}_{j}},{{E}_{0}})} {{S}_{e}}({{E}_{j}}) = \\ = \frac{1}{{{{\varepsilon }_{{{\text{Si}}}}}}}\sum\limits_i {{{N}_{i}}(LE{{T}_{i}},{{E}_{0}}){{S}_{e}}(LE{{T}_{i}})} . \\ \end{gathered} $
Для анализа удобно также использовать величины, относящиеся к количеству определенных вторичных ионов:
(7)
$\begin{gathered} {{f}_{{eh}}}({{E}_{j}},{{E}_{0}}) = \frac{1}{{{{\varepsilon }_{{{\text{Si}}}}}}}N({{E}_{j}},{{E}_{0}}){{S}_{e}}({{E}_{j}})\,\,\,\,{\text{и}} \\ {{f}_{{eh}}}(LE{{T}_{i}},{{E}_{0}}) = \frac{1}{{{{\varepsilon }_{{{\text{Si}}}}}}}N(LE{{T}_{i}},{{E}_{0}}){{S}_{e}}(LE{{T}_{i}}). \\ \end{gathered} $
Тогда количество электронно-дырочных пар, образующихся на единице длины трека под воздействием всех вторичных ионов:
(8)
$\begin{gathered} {{F}_{{eh}}}({{E}_{0}}) = \int\limits_0^{{{E}_{{\max }}}} {dE{\kern 1pt} '{{f}_{{eh}}}} (E{\kern 1pt} ',{{E}_{0}}) = \\ = \int\limits_0^{LE{{T}_{{\max }}}} {dLET{\kern 1pt} '{{f}_{{eh}}}(LET{\kern 1pt} ',{{E}_{0}})} . \\ \end{gathered} $
Максимальное количество неупругих столкновений с ядром в чувствительном объеме определяется соотношением:
(9)
${{N}_{{SV}}}({{E}_{0}}) = {{N}_{0}}({{E}_{0}}){{\sigma }_{{{\text{inel}}}}}({{E}_{0}}){{d}_{{SV}}}{{\rho }_{{{\text{Si}}}}},$
где ${{N}_{0}}({{E}_{0}})$ – количество первичных протонов, проходящих через чувствительный объем толщиной dSV, ${{\rho }_{{{\text{Si}}}}}$ – плотность материала мишени, [ат./см3], ${{\sigma }_{{{\text{inel}}}}}({{E}_{0}})$ – сечение неупругого взаимодействия с ядром, в результате которого образуется несколько вторичных частиц, и трек первичной частицы разделяется на несколько новых треков.

Для сбоя в работе интегральной схемы в одном из чувствительных объемов должно произойти, по крайней мере, одно неупругое столкновение протона с ядром ${{N}_{{SV}}}({{E}_{0}}) > 1.$ Вероятность этого события пропорциональна количеству попавших в чувствительный объем протонов ${{N}_{0}}({{E}_{0}}).$ Характерные размеры чувствительной области кремния в современных электронных схемах составляют dSV = 0.01–2 мкм [9], а критическое значение заряда Qc = 0.01–0.03 пКл [4]. В настоящей работе для оценки вероятности сбоя в электронике вследствие воздействия отдельного иона будем использовать значения dSV = 0.5 мкм и Qc = = 0.016 пКл. Тогда из соотношения (1) получаем, что деградация параметров и сбой в работе интегральной схемы может произойти при условии, если сумма электронно-дырочных пар, образующихся под влиянием первичного протона и продуктов его взаимодействия с атомом и ядром кремния в чувствительном объеме, превышает критическое значение (Feh)c = 200 нм–1. Эта величина для кремния соответствует значению LETc = = 3.1 МэВ · см2/мг.

РЕЗУЛЬТАТЫ РАСЧЕТОВ

Первичные протоны и ионы отдачи

При прохождении протона с энергией E0 > 1 МэВ через мишень из кремния образуются ионы отдачи Si+ с энергией ESi < 0.3 кэВ [8]. Столкновения ионов Si+ с атомами мишени также могут приводить к появлению дополнительных электронно-дырочных пар. Большинство электронов в этих столкновениях образуется в процессе ионизации внешних оболочек атома. В общем случае дрейф и диффузия этих заряженных частиц с разной подвижностью, а также их накопление на границе раздела разных материалов зависят от температуры [10] и режима работы интегральной схемы.

Первичные протоны и ионы Si+ при ионизации атома быстрыми протонами создают в чувствительном объеме электронно-дырочные пары, количество которых $F_{{eh}}^{b}({{E}_{0}})$ в работе интегральной схемы (с выбранными в настоящей работе значениями) к сбою от одного первичного протона привести не может, так как $F_{{eh}}^{b}({{E}_{0}})$ < (Feh)c (рис. 1). Но такой сбой может произойти в менее жестких, с точки зрения радиационной стойкости, устройствах. Результаты (рис. 1) показывают, что вероятность сбоя, вызванного ионизационными потерями энергии первичной частицы, достигает максимума в области энергии E0 ≈ 10 МэВ и определяется, прежде всего, ионами отдачи Si+.

Рис. 1.

Количество электронно-дырочных пар, создаваемых первичным протоном и ионами отдачи Si+ при прохождении ими мишени из кремния толщиной 1 нм. Результаты расчетов по программе SRIM для: 1 – протонов; 2 – ионов отдачи Si+; 3 – суммы протонов и ионов отдачи.

Вторичные ионы, образующиеся в результате столкновений с ядром в чувствительном объеме

Столкновение первичного иона с ядром атома – событие редкое, поскольку сечение такого неупругого процесса на несколько порядков меньше сечения ионизации атома. Для прохождения протонов через слой кремния толщиной d  ' = 0.5 мкм условие ${{N}_{{SV}}}({{E}_{0}}) > 1$ приводит к соотношению ${{N}_{0}}({{E}_{0}}) > {{10}^{6}}.$ Однако при неупругом столкновении с ядром рождается нескольких вторичных частиц, потери энергии которых могут значительно превышать фон $F_{{eh}}^{b}({{E}_{0}}),$ создаваемый в ионно-атомных столкновениях. Как правило, для протонов с E0 < 1 ГэВ при неупругом столкновении с ядром атома образуется один тяжелый ион с зарядом ядра Z ≥ 3 и несколько легких заряженных фрагментов (p, d, t, 3He, α, π ±). В этом случае ${{\sigma }_{{{\text{inel}}}}}({{E}_{0}})$ = ${{\sigma }_{{ZA}}}({{E}_{0}}),$ где ${{\sigma }_{{ZA}}}({{E}_{0}})$ – сечение образования одного вторичного иона с Z ≥ 3. Для протонов с энергией E0 > 1 ГэВ соотношение между сечениями ${{\sigma }_{{{\text{inel}}}}}({{E}_{0}})$ и ${{\sigma }_{{ZA}}}({{E}_{0}})$ может измениться, так как становятся возможными процессы рождения нескольких тяжелых ионов в процессе деления ядра, а также его фрагментация, когда ядро распадается на легкие частицы, и тяжелого иона с Z ≥ 3 может не остаться.

Если длина хорды для тяжелого вторичного иона близка к dSV, то возможна ситуация, при которой Feh > (Feh)c, и общее количество электронно-дырочных пар в чувствительном объеме может стать достаточным для накопления избыточного заряда Qdep> Qc (рис. 2). Результаты расчетов по программе GEANT4 [11] (рис. 2) показывают, что вклад тяжелых фрагментов ядерной реакции в общее количество электронно-дырочных пар является доминирующим во всем энергетическом диапазоне из-за небольших потерь энергии легких вторичных ионов с LET ≤ 1.5 MэВ · см2/мг. Однако учет концентраций электронно-дырочных пар, создаваемых легкими вторичными ионами, может стать значимым для выполнения условия Feh > (Feh)c в области энергии E0 > 1 ГэВ из-за увеличения количества легких фрагментов ядерной реакции. Поправка на количество электронно-дырочных пар, создаваемых первичным протоном с энергией E0 > 20 МэВ, и ионами отдачи Si+ (рис. 1) влияет слабо.

Рис. 2.

Количество электронно-дырочных пар, создаваемых в неупругих столкновениях протона в кремнии. Результаты расчетов продуктов ядерной реакции по программе GEANT4: 1 – легкие ионы (Z ≤ 2); 2 – тяжелые ионы (Z ≥ 3); 3 – сумма легких и тяжелых вторичных ионов, образующихся в результате взаимодействия с ядром Si; 4 – результаты расчетов по программе SRIM для суммы первичных протонов и ионов отдачи Si+.

Диапазон энергии в распределении количества тяжелых вторичных ионов Z ≥ 3 при одном неупругом столкновении c ядром составляет от сотен кэВ до десятков МэВ (рис. 3а). Каждому иону с зарядом ядра Z, массой A и энергией E можно поставить в соответствие неупругие потери энергии LET(E) (3) и представить распределение количества вторичных ионов NZA(E) в виде NZA(LET) (рис. 3б). Эти распределения характеризуются параметрами, определяющими диапазон их изменения ($E_{{ZA}}^{{\max }},$ LETmax) и средними значениями $\left\langle {{{E}_{{ZA}}}} \right\rangle ,$ $\left\langle {LET} \right\rangle .$ Соотношение ${{N}_{{SV}}}({{E}_{0}}) > 1$ является необходимым, но не достаточным условием сбоя электроники под действием одного быстрого иона. Не все ионы в распределении NZA(LET) (рис. 3б) могут привести к накоплению достаточного избыточного заряда. При небольшой длине трека вторичных ионов в чувствительном объеме (d ' ≤ dSV) или в случае их небольших потерь энергии накопленного избыточного заряда Qdep может оказаться недостаточно (Qdep < Qc).

Рис. 3.

Зависимость количества тяжелых (Z ≥ 3) вторичных ионов от их энергии EZA (а) и линейной передачи энергии LET (б) в точке, в которой произошло неупругое столкновение с ядром. Результаты расчетов по программе GEANT4 при энергии первичного протона E0: 1 – 0.1; 2 – 1; 3 – 10 ГэВ.

Вероятность того, что при неупругом столкновении рождается тяжелый ион, потери энергии которого при прохождении чувствительного объема будут достаточны для накопления избыточного заряда Qdep > Qc, можно вычислить с помощью отношения:

(10)
$P_{{\text{c}}}^{{\text{0}}}({{E}_{0}}) = {{\int\limits_{LE{{T}_{{\text{c}}}}}^{LE{{T}_{{\max }}}} {dLET{\kern 1pt} 'f_{{eh}}^{0}} (LET{\kern 1pt} ',{{E}_{0}})} \mathord{\left/ {\vphantom {{\int\limits_{LE{{T}_{{\text{c}}}}}^{LE{{T}_{{\max }}}} {dLET{\kern 1pt} 'f_{{eh}}^{0}} (LET{\kern 1pt} ',{{E}_{0}})} {\int\limits_0^{LE{{T}_{{\max }}}} {dLET{\kern 1pt} 'f_{{eh}}^{0}(LET{\kern 1pt} ',{{E}_{0}})} }}} \right. \kern-0em} {\int\limits_0^{LE{{T}_{{\max }}}} {dLET{\kern 1pt} 'f_{{eh}}^{0}(LET{\kern 1pt} ',{{E}_{0}})} }} \leqslant 1.$
При увеличении энергии E0 ширина распределения NZA(LET) (рис. 3б) и параметр $P_{{\text{c}}}^{0}({{E}_{0}})$ уменьшаются (табл. 1). Тогда достаточное условие сбоя в работе интегральной схемы под действием тяжелых вторичных ионов, которые рождаются, можно записать в виде:
(11)
${{N}_{{SV}}}({{E}_{0}})P_{{\text{c}}}^{0}({{E}_{0}}) > 1.$
Таким образом, для сбоя в работе интегральной схемы необходимо не только чтобы в чувствительную область попало достаточно большое количество первичных ионов ${{N}_{0}}({{E}_{0}}),$ но и чтобы в распределении вторичных частиц, которые рождаются в результате взаимодействия с ядром, доля вторичных ионов с большими неупругими потерями энергии была значительной.

Таблица 1.  

Параметры распределения вторичных тяжелых ионов (рис. 3)

E0, ГэВ EZA〉, МэВ LET〉, MэВ · см2/мг LETmax, MэВ · см2/мг Feh〉, нм–1 $P_{{\text{c}}}^{0}({{E}_{0}})$
  0.1 2.3 5.81 13.8 374.8 0.785
  1.0 4.8 5.57 13.6 351.6 0.775
10.0 5.5 5.20 12.3 181.6 0.438

Примечание. Значения ${{P}_{{\text{c}}}}({{E}_{0}})$ вычисляли для LETc = 3.1 МэВ · см2/мг.

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

Пробеги вторичных ионов зависят от энергии E0 и в случае быстрых протонов с E0 > 100 МэВ могут превышать десятки микрометров. Взаимодействие первичного протона может произойти с одним из материалов за пределами чувствительной области. Прохождение вторичных ионов, образующихся в результате такого взаимодействия, через одну или несколько чувствительных областей интегральной схемы также может привести к накоплению в них избыточного заряда и изменению токов утечки. Состав и потери энергии таких вторичных ионов зависят от материла мишени в интегральной схеме (Si, O, Al, Ti, W и так далее) и расстояния R от точки, в которой произошло неупругое столкновение с ядром атома, до чувствительной области.

При удалении от точки, в которой произошло неупругое столкновение с ядром, вторичные ионы с зарядом ядра Z и массой A замедляются, а их средняя энергия EZA и LET изменяются (рис. 4). В результате уменьшается количество вторичных ионов, которые могут привести к накоплению избыточного заряда в чувствительном объеме, достаточному для сбоя интегральной схемы. Распределение вторичных тяжелых ионов на расстоянии R будем характеризовать функцией ${{f}_{{eh}}}(LET{\kern 1pt} ',{{E}_{0}},R).$ Ее можно вычислить для каждого иона из распределения в точке, в которой произошли неупругое столкновение (R = 0) $f_{{eh}}^{0}(LET{\kern 1pt} ',{{E}_{0}})$ и потери энергии Se(EZA) + SN(EZA) [8], после прохождения этим ионом расстояния R. При увеличении расстояния R количество тяжелых вторичных ионов:

(12)
${{N}_{{ZA}}}({{E}_{0}},R) = \int\limits_0^{LE{{T}_{{\max }}}} {dLET{\kern 1pt} '{{f}_{{eh}}}} (LET{\kern 1pt} ',{{E}_{0}},R),$
уменьшается в результате их замедления ${{N}_{{ZA}}}({{E}_{0}},R)$ < $ < {{N}_{{ZA}}}({{E}_{0}}).$ Изменяется также доля ионов, у которых LET > LETc. Вероятность сбоя электроники под действием одного вторичного тяжелого иона (Z ≥ 3), образующегося в результате неупругого столкновения с ядром атома, на расстоянии R от чувствительной области записывается в виде:
(13)
${{P}_{{\text{c}}}}({{E}_{0}},R) = {{\int\limits_{LE{{T}_{{\text{c}}}}}^{LE{{T}_{{\max }}}} {dLET{\kern 1pt} '{{f}_{{eh}}}} (LET{\kern 1pt} ',{{E}_{0}},R)} \mathord{\left/ {\vphantom {{\int\limits_{LE{{T}_{{\text{c}}}}}^{LE{{T}_{{\max }}}} {dLET{\kern 1pt} '{{f}_{{eh}}}} (LET{\kern 1pt} ',{{E}_{0}},R)} {\int\limits_0^{LE{{T}_{{\max }}}} {dLET{\kern 1pt} '{{f}_{{eh}}}(LET{\kern 1pt} ',{{E}_{0}})} }}} \right. \kern-0em} {\int\limits_0^{LE{{T}_{{\max }}}} {dLET{\kern 1pt} '{{f}_{{eh}}}(LET{\kern 1pt} ',{{E}_{0}})} }} \leqslant P_{{\text{c}}}^{0}({{E}_{0}}).$
При прохождении в среде часть вторичных ионов останавливается, что приводит к уменьшению площади под кривой на рис. 4 и уменьшению количества ионов ${{N}_{{{\text{tot}}}}}$ после одного неупругого столкновения (табл. 2). Изменяются также параметры распределений оставшихся вторичных ионов ($E_{{ZA}}^{{\max }},$ LETmax, $\left\langle {{{E}_{{ZA}}}} \right\rangle ,$ $\left\langle {LET} \right\rangle $) и максимальное количество электронно-дырочных пар $\left\langle {{{F}_{{eh}}}} \right\rangle ,$ которое эти ионы могут образовать. Возрастание средней энергии ионов $\left\langle {{{E}_{{ZA}}}} \right\rangle $ при их замедлении объясняется уменьшением количества ионов с небольшой энергией и связано с переходом этой доли продуктов ядерной реакции с максимальными потерями энергии в область, прилегающую к пику Брэгга. Отметим, что даже для расстояния R в десятки микрометров в распределении вторичных ионов остаются ионы с большими потерями энергии LET > LETc, которые могут привести к накоплению избыточного заряда (Qdep > Qc) в одной или нескольких чувствительных областях интегральной схемы. Относительная доля таких ионов ${{P}_{{\text{c}}}}({{E}_{0}},R)$ быстро уменьшается с увеличением расстояния R. Следовательно, для регистрации эффекта воздействия неупругого столкновения с ядром атома в материале вблизи чувствительной области на концентрацию электронно-дырочных пар требуется существенное увеличение количества падающих на мишень ионов ${{N}_{0}}({{E}_{0}}) \gg {{10}^{6}}.$

Рис. 4.

Распределение количества тяжелых (Z ≥ 3) фрагментов ядерной реакции по линейной передаче энергии (LET) после столкновения протона с энергией E0 = 1 ГэВ с ядром кремния и прохождения вторичными ионами слоя кремния толщиной R. Результаты расчетов по программам GEANT4 для R = 0 (1) (рис. 3б) и SRIM для R: 1 (2); 5 (3); 20 мкм (4).

Таблица 2.  

Параметры распределения вторичных тяжелых ионов (рис. 4) для взаимодействия протонов с энергией E0 = 1 ГэВ и ядра атома кремния

R, мкм ${{N}_{{{\text{tot}}}}}$ EZA〉, МэВ LET〉, MэВ см2/мг Feh〉, нм–1 ${{P}_{{\text{c}}}}({{E}_{0}},R)$
0.0 1.000 4.0 6.22 393.1 0.775
1.0 0.688 3.7 5.29 234.7 0.492
5.0 0.106 6.7 5.68 38.8 0.149
20.0 0.009 8.8 4.29 2.5 0.019

Примечание. Значения ${{P}_{{\text{c}}}}({{E}_{0}},R)$ вычисляли для LETc = 3.1 МэВ · см2/мг.

Результаты расчетов (рис. 5) показывают, что максимум в распределении количества электронно-дырочных пар приходится на область энергии E0 > 1 ГэВ. При удалении от чувствительной области на расстояние R > 1 мкм количество электронно-дырочных пар, создаваемых вторичными ионами в столкновении Ti(p, x) при их прохождении чувствительной области толщиной dSV ≤ 0.5 мкм, будет недостаточным для накопления избыточного заряда Qdep> Qc.

Рис. 5.

Количество электронно-дырочных пар в чувствительной области, создаваемыми тяжелыми вторичными ионами при одном неупругом взаимодействии p–Ti. Результаты расчетов по программам GEANT4 и SRIM для расстояния R: 1 – 0; 2 – 0.5; 3 – 1 мкм.

Область, в которой внешний источник фрагментов ядерной реакции может изменить количество электронно-дырочных пар в чувствительном объеме Feh, ограничена. Будем характеризовать размер этой области критическим расстоянием Rc. Способов определения границы области R < Rc может быть несколько. Например, параметр Rc можно определить из расстояния, на котором количество электронно-дырочных пар Feh, создаваемых вторичными ионами, на порядок превышает фон, создаваемый одним первичным протоном [12]. Параметр Rc можно уточнить, принимая во внимание изменение параметров распределения ${{N}_{i}}(LE{{T}_{i}}).$ В случае распределения с $\left\langle {LET} \right\rangle $ = LETc потери энергии приблизительно половины вторичных ионов превышают значение LETc. Кроме того, из распределения на рис. 6 можно найти критическое расстояние Rc, при котором LETmax(Rc) ≤ ≤ LETc. При расстоянии R > Rc в распределении тяжелых ионов NZA(LET) не останется ионов, которые при прохождении чувствительной области толщиной dSV могли бы привести к накоплению заряда, достаточного для сбоя в работе интегральной схемы.

Рис. 6.

Изменение параметров распределения тяжелых (Z ≥ 3) ионов, образующихся в столкновениях протона c ядром кремния, при прохождения слоя кремния толщиной R. Результаты расчетов по программам GEANT4 и SRIM для протонов с энергией E0 = 1 ГэВ.

Увеличение критического расстояния Rc при возрастании E0 (рис. 7) объясняется увеличением средней энергии вторичных ионов $\left\langle {{{E}_{{ZA}}}} \right\rangle $ (табл. 1) и более высокими значениями их пробегов. Величина Rc может достигать сотен микрометров, что на порядок превышает глубину чувствительного слоя. При прохождении протонов с энергией E0 ≥ ≥ 2 ГэВ через материалы интегральной схемы величина Rc возрастает при уменьшении заряда ядра атома мишени. Это связано с увеличением доли более легких фрагментов ядерной реакции, которые характеризуются более высокими пробегами в кремнии.

Рис. 7.

Зависимость критического расстояния Rc от энергии E0 для распределения вторичных тяжелых ионов при взаимодействии протонов с материалами интегральной схемы. Результаты расчетов по программам GEANT4 и SRIM: 1 – в кремнии; 2 – в титане; 3 – в вольфраме.

ЗАКЛЮЧЕНИЕ

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

Результаты расчетов показали, что ионизационные столкновения одного первичного протона с энергией E0 > 20 МэВ в чувствительном объеме создают небольшое количество электронно-дырочных пар, и накопленного избыточного заряда для сбоя в работе интегральной схемы с высокой радиационной стойкостью, как правило, недостаточно.

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

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

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

Список литературы

  1. Модель космоса. Том II. Воздействие космической среды на материалы и оборудование космических аппаратов / Ред. Новиков Л.С. М.: Изд-во НИИЯФ МГУ, 2007. 1145 с.

  2. Chechenin N.G. // J. Nucl. Eng. Technol. 2012. V. 2. № 1. P. 22.

  3. Чуманов В.Я., Кадменский А.Г., Чеченин Н.Г. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2014. № 12. С. 36.

  4. Petersen E.L., Pickel J.C., Adams J.H. et al. // IEEE Trans. Nucl. Sci. 1992. V. 39. P. 1577.

  5. Holbert K. E. // Single Event Effects. 2008. http://holbert.faculty.asu.edu/eee560.

  6. Read R.A., Weller R.A., Akkerman A. et al. // IEEE Trans. Nucl. Sci. 2013. V. 60. P. 1876.

  7. Чумаков А.И. Действие космической радиации на интегральные схемы. М.: Радио и связь, 2004. 320 с.

  8. Ziegler J., Biersack J.P. // SRIM: the Stopping and Range of Ions in Matter (www.srim.org).

  9. Dodd P.E., Schwank J.R., Shaneyfelt M.R. et al. // IEEE Trans. Nucl. Sci. 2007. V. 54. P. 2303.

  10. Таперо К.И., Улимов В.Н., Членов А.М. Радиационные эффекты в кремниевых интегральных схемах космического применения. М.: Изд-во БИНОМ. Лаборатория знаний, 2014. 302 с.

  11. Allison J., Amako K., Apostolakis J. et al. // Nucl. Instrum. Methods Phys. Res. A. 2016. V. 835. P. 186. https://geant4.web.cern.ch

  12. Новиков Н.В., Чеченин Н.Г., Чувильская Т.В. и др. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2020. № 2. С. 72.

Дополнительные материалы отсутствуют.