Анализ конструктивно-технологических ограничений при проектировании лавинных фотодиодов, работающих в режиме счета фотонов
1 2 11 Ю.Ф. Адамов , А.В. Верховцева , И.О. Шурчков , Е.С. Балака
Институт проблем проектирования в микроэлектронике РАН, Москва
Институт радиотехники и электроники им. В. А. Котельникова РАН, Москва
Введение
В перечне быстродействующих фоточувствительных элементов особое место занимают фотодиоды с лавинным умножением носителей в сильном электрическом поле. Структура лавинного фотодиода (ЛФД) отличается от структуры обычного р4-п фотодиода введением тонкого легированного р+ слоя на границе с п+ слоем. В четырехслойной структуре ЛФД между контактными п+, р+ областями присутствуют слаболегированная светопоглощающая область и сильнолегированная область, в которой происходит умножение носителей. Процесс умножения имеет управляемый характер только в случае, когда коэффициент умножения носителей не превышает отношения коэффициентов ионизации электронов и дырок. В кремнии это отношение всегда меньше 10, а для реальных размеров фотодиодов менее 5. Из этого следует, что максимальное усиление фототока в интегральном ЛФД не может значительно превышать величину 5 в стационарном режиме [1].
Перспективным считается импульсный режим включения ЛФД, в котором внешняя цепь прерывает ток нагрузки и снижает напряжение на диоде при развитии неуправляемой лавины с коэффициентом умножения, стремящимся к бесконечности. В этом режиме достигается предельная чувствительность к обучению, когда каждому поглощенному фотону соответствует импульс в нагрузке фотодиода [2].
Разработка фотоприемников на основе работы ЛФД в режиме счета фотонов требует исследования процессов зарождения, прерывания и стационарного состояния лавины с учетом конструктивно-технологических ограничений.
Ударная ионизация и лавинное умножение в полупроводниках и полупроводниковых приборах.
Как известно, в полупроводниках, под действием электрического поля, носители заряда приобретают определенную направленную скорость цЕ, а, следовательно, и дополнительную энергию цЕ2 т, где ц- подвижность, а т- так называемое время релаксации энергии. Эта энергетическая добавка перераспределяется процессами испускания и поглощения акустических фононов и превращается в электронную температуру Те = Т0 + цЕ2 т, превышающую температуру кристаллической решетки, т.е. фононной подсистемы Т0.
При более подробном исследовании процессов электроразогрева носителей, выделяют область малых электрических полей Е < 103 В/см, когда излучаются только акустические фононы, которые «уносят» электронную энергию сравнительно маленькими порциями. При этом температура электронного газа остается меньше энергии оптического фонона, поэтому последние и не излучаются.
В области промежуточных полей 103 В/см < Е < 104 В/см процесс терморелаксации усиливается и усложняется за счет испускания оптических фононов. При этом существенно уменьшается время релаксации импульса, а дрейфовая скорость носителей перестает зависеть от электрического поля (выходит на насыщение).
В области больших полей Е >10 В/см электронная температура увеличивается примерно пропорционально полю Е и достигает достаточно больших значений (порядка нескольких тысяч градусов), что и активизирует ударную ионизацию зона-зона.
Как правило, при интерпретации экспериментальных данных, общепринята довольно простая аппроксимация полевой зависимости коэффициентов ударной ионизации
(l)
Конкретные значения параметров А, Ь и т для основных полупроводников в наше время надежно установлены.
Ионизационные акты могут инициироваться как достаточно разогретым электроном так и не менее горячей дыркой. Именно эта биполярность ударной ионизации обуславливает эффективную внутреннюю положительную обратную связь ионизационного процесса в области сильного поля р-п перехода. Именно эти парциальные обратные связи суммируются уравнением непрерывности тока в известных моделях стационарного лавинного пробоя, в которых напряжение лавинного пробоя определяется, как особенность соответствующего коэффициента умножения:
Поскольку рассматриваемые в диссертационной работе лавинные фотодиодные структуры работают именно в области надпробойных напряжений, а их реакция на облучение имеет ярко выраженный импульсный (нестационарный) характер, при построении адекватного теоретического описания функционирования лавинных фотодиодов, модель стационарного пробоя оказывается практически бесполезной [3,4]. В настоящей работе в качестве исходного момента используется лишь формула (1) с экспериментально установленными значениями коэффициентов:
Рис. 1. Поперечный разрез структуры планарного кремниевого лавинного фотодиода, сформированного на подложке р-типа с удельным сопротивлением р ~ 10 Ом-см (р0 ~
Глубина р-слоя ё определяющим образом влияет на напряжение лавинного пробоя Уп. Отметим, для конкретики, что при ё = 0.7 мкм Уп ~ 20 В, а уже при ё = 1 мкм Уп ~ 30 В. В этом диапазоне напряжений нарастание величины электрического поля, обусловленное
позволяет считать электрическое поле между слоями постоянным E = const.
Во втором разделе, в пренебрежении актами дырочных ионизаций, методами классической статистики исследован процесс формирования электронных сгустков в области сильного поля из единственного электрона, стартующего с катода. Получена функция распределения кратности электронного умножения при однократном пролете области сильного поля:
dx
(2)
(3)
Физическая модель интегрального фотодиода в режиме счета фотонов.
Простейшая структура лавинного фотодиода приведена на рис. 1.
1015 см-3) с помощью двух базовых имплантаций: поверхностной п++ с дозой 10-100 мкКл/см2 и заглубленной р+ с дозой порядка 1 мкКл/см2
концентрацией фоновой примеси (р0 ~ 1015 см'3'), оказывается сравнительно малым, что и
Функция распределения (4), представленная на рис.2 при трех различных значениях параметра аё, максимальна при п = 1 (безионизационый пролет), демонстрирует значительную протяженность в сторону больших п. Ее максимальное значение уменьшается с ростом аё, но ее протяженность в область больших п увеличивается более значительно, так, что среднее число электронов в сгустке у анода существенно растет.
ад
(п) = 2 пРп = ехр (а). (5)
п-1
р.,
0.12
0.10
0.0$
0.06
0.04
0.02
0.00
Рис. 2. Распределение вероятности кратности электронного умножения при различных значениях параметра аё: 1 - аё = 2; 2- аё = 3; 3 - аё = 4
Следует обратить внимание также на аномально широкую дисперсию этого распределения при больших средних (п)
((п - (П) ^ = ехр (а/) • (ехр (а/) -1) = (п) 2 - (П• (6)
Адекватность полученных классически особенностей функции распределения кратности электронного умножения подтверждена математическим моделированием по методу Монте-Карло.
Вводится понятие коэффициента усиления петли обратной связи лавинного процесса, для замыкания которой необходим учет дырочных ионизаций. Предполагая выполненным условие р/а << 1 (в силу асимметрии коэффициентов ударной ионизации), в рассматриваемой статистике можно ограничиться учетом лишь двух альтернативных событий: 1) пролет дыркой интервала а-1 без ионизации с вероятностью (1- Ра-1), и 2) пролет с единичной ионизацией, т.е. генерацией нового электрона с вероятностью ра-1 << 1.
Использование такой бинарной функции распределения вероятности для дырок совместно с функцией распределения (4) позволяет записать вероятности различных реализаций первой стадии развития лавинного процесса, заканчивающейся генерацией дырками некоторого количества новых электронов в прикатодном слое пролетной области протяженностью а-1. Соответствующие аналитические выражения имеют вид:
1 -РТ-1,
а )
Q(1) = 2(п- 1)Рехр(-а/М1 -ехр(-а/))”-1^1 -Р] , (7)
п=1 а V а)
6(2) =]Г(п 1Хп 2) ^ ехр (-а/ )^(1 - ехр (-а ))И-1[^1 -р •
6(0) = 2 ехр(-а) •(1 - ехр(- а ))п-11
п =1
После сравнительно несложных преобразований формулы (7), при условии exp (ad)>> 1 >> —, трансформируются в довольно простые выражения
a
Q(0) ={l + q Г, Qa) = q{l + q)
(8)
(0) - q {1 + q) > Q(3)
Q(0)- q2 {l + q\\ Q(3)- q3 {l + q)
— „ad
где обозначено q — — e
a
Обобщая соотношения (8), получим искомую формулу для вероятности того или иного значения к - парциального коэффициента усиления единичной петли обратной связи:
1
к
Qk) — -^- .к — °А2... (9)
{) q +111 + q)
Эта функция распределения парциальных значений коэффициента усиления петли обратной связи (9) нормирована на единицу. Как и функция (4), выражение (9) максимально при минимальном к (к = 0), последние означает, что наивероятнейшим результатом регенерации единичного электрона является негативное событие, т.е. выключение лавины уже на первой стадии ее развития. Функция распределения (9) дает правильную величину среднего коэффициента регенерации:
к—q-——E) exP {d a{E) (10)
a(E)
Итак, предположим, что вследствие неких причин в момент времени ti в окрестности катодного края (i-p переход, х = 0) диода появилось ni(0) электронов. В свою очередь для каждого из этих электронов имеется вероятность совершить da(E) > 1 актов ударной ионизации на пути к анодной границе (i-n переход, х = d). При этом в среднем за время пролета tnp= d/vs (vs - скорость насыщения) в n+- область синфазно уйдут n(0) электронных сгустков со средним числом exp(d a(E)) электронов в каждом, то есть
и;+i{d) — П (о)-exp {d •a{E)) (11)
В результате тех же актов электронной ударной ионизации в течение временного интервала At = 2tnp окрестность катодной границы пересечет такое же количество дырок. Эти дырки, несмотря на сравнительно небольшую величину ионизационного интеграла, создадут в слое протяженностью а-1 у катодной границы (х = 0) в момент времени ti+1= ti + At новые электроны в количестве
—(e) _
<Е )
которые и «запустят» следующую стадию формирования лавинного процесса.
В четвертом разделе второй главы построен эффективный алгоритм расчетной процедуры определения формы токовых импульсов и приведены результаты численного моделирования.
Чтобы проследить развитие процесса разрядки емкости диода при превышении напряжения на диоде над пороговым, необходимо, наряду с формулой для среднего коэффициента усиления (10), также учитывать изменение напряжения на лавинном фотодиоде, определяющего численное значение коэффициента (10). При этом диод с одной стороны, разряжается внутренним током электронных сгустков и дырочных «хвостов», а с другой - дозаряжается током от источника питания Vu через нагрузочное сопротивление RH. Упрощенная эквивалентная схема, иллюстрирующая этот процесс, показана на рис 3(а).
n+1{о)—n{о) • exP{d -a {e )) • —^ — q • П {о), (12)
Ин
б)
Рис. 3. а)Упрощенная эквивалентная схема однофотонного лавинного фотодиода. б) Характерные временные зависимости напряжения на диоде (результаты численного
моделирования ё = 0,7 мкм, С = 250 фФ)
Приведенной схеме включения отвечает парциальное изменение напряжения на диоде за интервал At = 2^р дискретизации текущего времени:
V, = V -
’г +1 ’г
1 (с1 )-
(V - V К
я
■ с 1,
(13)
где С = еа 0---емкость диода, е - элементарный заряд.
(I
Система уравнений (11, 12, 13) представляет собой простую итерационную модель развития лавинного процесса в фотодиоде. Для корректной реализации модели использовались только целочисленные значения щ что особенно важно для числа инициирующих электронов. При этом, как только п(0) становится равным нулю в предложенной итерационной схеме, лавинный пробой автоматически «выключается», что позволяет току во внешней цепи через Ян восстановить напряжение на диоде до исходного напряжения источника питания. При моделировании использовались следующие начальные значения параметров: У0=Уи - напряжение источника питания, п00(ф = 0 - число электронов у анодной границы в начальный момент времени t0, п0(0) = 1 - первичный фотоэлектрон у катодной границы в начальный момент времени 1;0. Остальные необходимые для расчета технологические параметры лавинного фотодиода и эквивалентной схемы: емкость диода С, длину пролетной области ё, сопротивление нагрузки Ян - легко задавать и варьировать.
На рис. 3(б) приведены зависимости напряжения на лавинном фотодиоде, включенном одиночным электроном, от времени для напряжения источника Уи = 26 В при нескольких характерных значениях сопротивления нагрузки Ян = 1, 2, 4 кОм. Оказалось, что значению Ян = 1 кОм отвечает своеобразное апериодическое «включение» стационарного пробоя с напряжением Уп = 23,8 В и током, определяемым падением избыточного напряжения (перенапряжения) (Уи - Уп) на нагрузочном сопротивлении. А при Ян > 2 кОм происходит восстановление напряжения питания на фотодиоде в результате полного выключения лавинного тока во временном промежутке, когда напряжение на диоде было ниже порогового.
Следует отметить практически одинаковое поведение всех этих кривых на первом этапе релаксационного процесса, который распространяется до напряжения существенно меньшего, чем Уп, что обусловлено специфической инерционностью лавинного процесса. Анализируя результаты построенной простой теоретической модели работы кремниевого лавинного фотодиода, можно сформулировать несколько основных выводов:
— в режиме лавинного пробоя в области умножения кремниевого фотодиода подавляющее число ионизаций инициируется электронами (а > в);
с
— преобладание электронных ионизаций компенсирует недостаток дырочных, так что результирующий коэффициент положительной обратной связи не превышает разумных значений ~ 2,5;
— форму однофотонного импульса отличает быстрая стадия лавинной разрядки емкости диода с характерной длительностью около 0,1 - 0,3 нс, и относительно медленная стадия (восстановления) дозарядки, длительность которой определяют электрические параметры внешней цепи;
— при слишком быстрой дозарядке диод переходит в состояние стационарного лавинного пробоя с напряжением V = Уп и стационарным током (Уи - Уп)/Ян.
Динамика токовых импульсов с учетом латеральной неоднородности.
Рассматривается сложная пространственная структура задачи, а именно, высокая локализация исследуемого нестационарного процесса лавинной перезарядки по площади диодной структуры. Поскольку главная фаза нестационарного процесса лавинной разрядки емкости диода, инициированная единичным электроном в какой-то случайной точке его площади, столь быстротечна (несколько десятых наносекунд), естественно предположить, что в это малое время перезарядиться до глубоко подпороговых значений сможет относительно малая часть площади диода. Поэтому в конкурирующем процессе дозарядки емкости диода от источника напряжения появится дополнительная временная задержка на растекание локального инжектированного заряда по площади катодного коллекторного р- слоя структуры (рис. 1), что может кардинально проявиться в кинетике дозарядки, а, следовательно, повлиять (количественно) на соответствующий критерий самовыключения лавинного процесса.
Рассматривается электрическая реакция лавинного фотодиода (рис. 4(а)) на «прибытие» в центральную (для простоты) точку р+ - катода диода очередной порции положительного дырочного заряда, в характерный промежуток времени Лt = 21пр ~ 2*10-11 с.
а)
б)
Рис. 4. а) Условная картина растекания электрического тока при локальном пробое лавинного фотодиода. Уи - напряжение источника питания, Ян - сопротивление нагрузки, Ус - сигнальное напряжение. б) Результаты моделирования возможных режимов ЛФД с емкостью 90 фФ и сопротивлением 10 кОм для Ян = 30 кОм (кривая 1) и Ян =12 кОм (кривая 2). Для сравнения приведена кривая релаксации (3) ЛФД с малым (0.1 кОм) сопротивлением растекания
Вводится определение локального потенциала в катодном р+-слое диодной структуры ф(г,(), зависящего от времени и внутренней латеральной координаты г, тогда ток во внешней цепи задается интегральным соотношением:
где р8
ъ.
р*К пАК К ■К
10 Ом-см удельное сопротивление подложки, ~ 10 мкм
(14)
толщина ее слоя между
р - обкладкой и заземляющей шиной, Я - площадь диода (р -слоя). Соответственно полный емкостной заряд определяется выражением:
( К
Я = к. С - К + 1 С|ЁКр(г,г) = е0-ЩЗ(г), (15)
V К
где удельная емкость диода С0 = ££0Ш ~ 10" Ф/см . Комбинируя выражения (14), (15), получим
(К + к )с а
Отсюда следует интересная особенность «сконструированной» для упрощения задачи внешней цепи (протяженный электрод заземления), заключающаяся в том, что любой нестационарный заряд 0 релаксирует с одной и той же постоянной времени С(Яп+Ян) вне зависимости от формы его распределения по площади р+ - обкладки конденсатора (диода). Тем не менее, это обстоятельство необходимо учитывать для расчета текущих значений коэффициентов умножения а и в, поскольку для этого необходимы текущие значения потенциала ф(0,0 = Фь при г = 0, т.е. в «точке» пробоя.
Для этого рассматривается известное решение двумерного телеграфного уравнения [9], описывающего центрально-симметричное расплывание заряда из окрестности точки г = 0
Р(-. г) = Ро -. ехр - , (17)
мг
г ч ^
где Я ~ 10 кОм - сопротивление р+ - слоя, ф0Х0 - постоянная интегрирования. В качестве ^0, естественно, выбирается характерное время задачи т ~ 2^пр ~ 2*10-11 с. Интегрирование по площади потенциала (17), умноженного на С0, дает соответствующий заряд
q = | акСоР(г. г ) = (Ро (18)
Здесь коэффициент 4пт№=АС, тем больший, чем меньше сопротивление растекания Я, отражает характерное расплывание отдельной порции локального дырочного заряда за время его «введения» в систему. Он играет роль дискреты эффективной емкости растекания, которая, увеличиваясь со временем, достигает номинальной емкости диода С за время I = пт, где п = С/АС- число временных интервалов т, за которое введенный заряд успевает распространиться по всей площади диода. При этом, в соответствии с формулой (18), интересующий потенциал центральной точки Ф1 уменьшается обратно пропорционально дискретному времени гт, достигая стационарного значения при г = п. Сама же формула (18) позволяет по величине парциального заряда Ц\, локально инжектированного на соответствующей стадии лавинного процесса в р+ - обкладку исследуемой диодной структуры, рассчитать амплитуду ф(0,г) - компоненту потенциала, отвечающего компоненте заряда
Предлагается модифицированная итерационная схема моделирования с учетом динамики растекания в катодном слое. Обозначив число электронов, инициирующих г-ю стадию процесса, генерированных в прикатодном слое протяженностью а"1, на предыдущей (г-1) стадии как N1, можно записать сначала парциальный заряд электронного сгустка, пришедшего на анод в г интервал дискретного времени:
ql = еЫ^еяр &(Г-Фг-1 )•а). (19)
Точно такой же величине равен и парциальный заряд дырок, пришедших на катод за временной интервал 2^пр, сгенерировав при этом N новых электронов для инициации следующей (г+1) стадии процесса
N = N1 -1 ехр {а(¥-Ф1-1). а \Р\ ~Ф 1 ’ (20)
а(к -ф1 -1)
Полный дополнительный заряд на обкладках диода увеличивается за счет поступления очередной порции лавинного тока и уменьшается, «вытекая» во внешнюю цепь:
(21)
где тп = (Яп+Ян)С.
Соответственно, ток во внешней цепи равен
1
1=у 2 qJ ехР \-—(/ -1 )г,
—
(22)
а напряжение на диоде Уи - 1Ян.
Чуть более сложный вид имеет формула для локального потенциала Ф в точке пробоя
Ф1=2
q,
ехр
тт Я/ -1 +1);
С
АС
(23)
Следует обратить внимание на то, что числитель в формуле (23) отвечает «утеканию» инжектированных зарядов во внешнюю цепь, а знаменатель - их «растеканию» из «точки» пробоя по площади диода.
На рис. 4(б) приведены кривые релаксации для значений Ян = 30 и 12 кОм. И в том, и в другом случае зависимости демонстрируют быструю лавинную стадию разрядки емкости диода до значений существенно ниже напряжения пробоя (горизонтальный участок линии 2).
Тем не менее, на стадии восстановления явно различаются два этапа. В быстром этапе
+
превалирует восстановление локального потенциала за счет растекания по площади р -катода дырочного заряда, инжектированного в точку пробоя. Этот этап достаточно близко подводит локальный потенциал (снизу) к потенциалу Уп стационарного режима лавинного процесса (д = 1), облегчая при этом задачу второго, более медленного этапа обычной ЯнС дозарядки (уход оставшегося неравновесного дырочного заряда во внешнюю цепь через подложку). Задачей второго этапа является смещение общего уровня потенциала так, чтобы его локальное (в точке пробоя) значение превзошло критический потенциал Уп. В этом случае лавина сможет развиться снова, если к тому времени в ней останется сколько -нибудь значимое количество инициирующих электронов. Последнее обстоятельство объясняет тот важный аспект, что для выключения лавинного процесса необходимо соответствующее замедление ЯнС дозарядки и увеличение критического сопротивления нагрузки с 1 до порядка 20 кОм.
Проделав значительную серию численных экспериментов для широкой области значений внешних параметров С и Я, были получены зависимости критического сопротивления нагрузки кк от сопротивления имплантированного слоя Я (рис. 5)
Рис. 5. Результаты серии численных экспериментов по определению значений критического сопротивления нагрузки: а) в линейном масштабе; б) в логарифмическом
п
масштабе. Кривые на графиках соответствуют следующим значениям емкости С, фФ:
1 - 10, 2 - 30 , 3 - 90, 4 - 270 Результаты, представленные на рис. 4(б) и рис. 5(а,б), убедительно доказывают важность рассмотренного физического явления (растекания тока) в формировании импульсной реакции лавинного фотодиода на однофотонное возбуждение. Поскольку именно учет своеобразной динамики растекания заряда из точки пробоя переводит расчетные значения требуемых нагрузочных сопротивлений в область реальных величин, соответствующих практике эксперимента и проектирования сенсорных устройств на основе лавинных фотодиодов. Не менее важно, что данная модифицированная модель позволяет определить необходимую дозу имплантации р+ - слоя, которую обычно, в целях минимизации дефектности активной области приборных структур, стремятся уменьшить.
Нестабильность начальной стадии лавинного процесса.
В начальной стадии развития лавинного процесса со сравнительно малым количеством носителей в области сильного поля в наибольшей степени проявляется случайность процесса регенерации тока [5].
Анализируя элементарный акт регенерации, получены функции распределения его результатов для одного, двух и трех электронов на «входе» единичной петли обратной связи. Показано, что с ростом числа носителей на входе петли обратной связи наблюдается тенденция к уменьшению относительных флуктуаций числа частиц в лавине уже на первых этапах развития процесса. Вероятность нулевой результативности единичной регенерации существенно подавляется, что позволяет сформулировать критерий эффективности, как по меньшей мере удвоение инициирующего электрона уже на первой стадии процесса размножения:
ш ш 1 ( q У ( q I2
* = 2 ^> = 21-1 гЬ)=( гЪ). (24)
к=2 к=2 1 + q V1 + q У V1 + У
Учитывая случайность координаты фотогенерации первичного электрона в области сильного поля, получена спектральная зависимость Кф(Х) усреднением аналога формулы (24) по функции поглощения световых квантов Р(х) = ц(Х). е , где 1](Л) -
коэффициент поглощения в кремнии:
к=КхЧ<т+&)К+^' (25)
Соответствующие зависимости Кф(Х) приведены на рис. 6(а)
0 ----------1--------1--------*-------г-------->--------г-------*--------1--------1-------1 1-----------I---------I----------,----------.---------.----------1----------,----------,---------.
2* 29 30 31 32 33 34 35 36 3? 38 062 0Я2 1.08 1.5 2 2.68 2.65 5 7 4
У,В
а) б)
Рис. 6. а) Зависимость эффективности фотоэлектрического преобразования Кф(Х) от напряжения при возбуждении монохроматическим излучением (1 - X = 0,5 мкм; 2 - X = 0,6 мкм; 3 - X = 0,8 мкм; 4 - К по формуле (24); б) Зависимости коэффициента эффективности
одноэлектронного возбуждения лавины от напряжения, рассчитанные: 1 - теоретически (по результативности одной петли обратной связи); 2 - упрощенная Монте-Карло модель (раздел 4.3); 3 - обобщенная Монте-Карло модель с высокой степенью дискретизацией
Графики 1,2,3 рис. 6(б) демонстрируют рост эффективности для более длинноволновых фотонов, поскольку последние генерируют первичный фотоэлектрон в прикатодном а"1 слое с большей вероятностью.
Чтобы преодолеть сомнения в адекватности оценки эффективности по результатам усиления первой петли обратной связи, была разработана специализированная Монте-Карло модель эволюции лавинного процесса на более длительном промежутке времени, включающим несколько (вплоть до десяти) петель обратной связи. Критерием неэффективности в модели служило отсутствие электронов на выходе из последней петли. Полученные численные результаты (кривая 2 на рис. 6(б)) демонстрируют ожидаемое существенное падение эффективности в окрестности напряжения пробоя и достаточную близость к аналитической зависимости (24) (кривая 1 на рис. 6(б)) при сравнительно больших надпробойных напряжениях. Полученный результат позволил оценить характерную длительность импульсов квазистационарного лавинного тока (микроплазменный режим), сформулировав критерий самовыключения лавинного процесса. Соответствующие вычисления показывают, что при сопротивлениях нагрузки порядка нескольких сотен КОм, длительность микроплазменных импульсов находится в диапазоне 0,1 -10 мкс, что согласуется с известными экспериментальными данными. Построена обобщенная Монте-Карло модель лавинного процесса с высокой степенью дискретизации пространственно-временных координат. В принципе, эта модель, свободная от использованных ранее упрощений (а<<в, Аt = 2tпр и т.п.), дает наглядную картину статистически разупорядоченных начальных стадий развития лавины.
Полученная в рамках этой модели, эффективность преобразования фотоэлектрона в сигнальный импульс (кривая 3 на рис. 6(б)), практически полностью совпадает с предыдущей зависимостью (кривая 2 на рис. 6(б)), что подтверждает правомочность разработанной модели включения, равно как и адекватность использованных упрощений при анализе импульсных переходных процессов в ЛФД.
Результаты работы
Новое направление быстродействующих детекторов слабого излучения основано на использовании ЛФД в режиме счета фотонов. Каждый пиксель включает матрицу ЛФД, а простая цифровая схема регистрирует суммарное число лавинных импульсов за интервал экспозиции. Этот метод обеспечивает предельно высокую чувствительность фотоприемника и высокое быстродействие. До настоящего времени не существовало методов и средств проектирования и оптимизации интегральных фотоприемников на основе ЛФД.
• Разработана новая модель перезарядки емкости лавинного фотодиода, учитывающая ограничения конструкции и профили легирования слоев и позволяющая оптимизировать конструкцию интегрального лавинного фотодиода для счета фотонов.
• Установлено, что в режиме счета фотонов емкость лавинного фотодиода разряжается до напряжений существенно меньше пробивного. Режим счета фотонов реализуется при задержке сигнала во внешней цепи не менее 1 нс, что на порядок больше времени нарастания тока лавинного импульса.
• Нестабильность фронта сигнального импульса определяется случайным характером регенерации лавинного пробоя и составляет 3 -5 времен пролета носителей.
• Теоретические расчеты и эксперименты показывают, что регистрация световых импульсов в режиме счета фотонов может осуществляться с разрешением 5 нс. Литература:
[1] Верховцева А.В., Ванюшин И.В., Гергель В.А., Зеленый А.П., Зимогляд В.А., Тишин Ю.И. Статистические особенности зарождения процесса лавинного умножения в
полупроводниках с различающимися коэффициентами ударной ионизации // Радиотехника и электроника. 2009. Т. 54. № 11. С. 1403-1414.
[2] Адамов Ю.Ф., Сибагатулин А.Г., Сомов О.А. Тенденции развития сенсорных систем и интеллектуальных датчиков // Датчики и системы. 2011. № 5. С. 58-59.
[3] Верховцева А.В. Моделирование динамики развития локального микропробоя в кремниевых лавинных фотодиодных структурах // Нелинейный мир. 2009. Т. 7. № 3. С. 212.
[4] Верховцева А.В. Исследование динамики микропробоя в кремниевых лавинных фотодиодных структурах // Сб. тезисов докладов 16-й Всероссийской межвузовской научно-технической конференции студентов и аспирантов «Микроэлектроника и информатика - 2009». М.: МИЭТ, 22-24 апреля 2009. С. 35.
[5] Ruida Yun, Valencia Joyner A Monolithically Integrated Phase-sensitive Optical Sensor for Frequency-Domain NIR Spectroscopy // IEEE Sensors Journal. July 2010. V. 10. № 7. pp. 12341242.