CC BY
173
К. А. Насыров, В. А. Гриценко, Ю. Н. Новиков, Д. В. Гриценко, Д.-В. Ли, Ч. В. Ким
МЕХАНИЗМ ПЕРЕНОСА ЗАРЯДА В НИТРИДЕ КРЕМНИЯ: ЭФФЕКТ ФРЕНКЕЛЯ, ИЛИ МНОГОФОНОННЫЙ МЕХАНИЗМ ИОНИЗАЦИИ ЛОВУШЕК1
Экспериментально и теоретически изучался механизм переноса заряда в аморфном нитриде кремния 81Ых. Обнаружено, что модель Френкеля, широко применяемая для объяснения ионизации ловушек, описывает экспериментальные результаты, но при этом необходимо использовать нефизично малый частотный фактор и аномально большую эффективную туннельную массу электрона и дырки. Экспериментальные данные количественно описываются многофононным механизмом ионизации ловушек.
По определению, проводимость диэлектрика должна равняться нулю. Тем не менее, Пулл в 1916 году обнаружил экспоненциальное увеличение проводимости слюды в зависимости от электрического поля [1]. В 1938 году Френкель предложил теоретическое объяснение «экспоненциального закона» для диэлектриков с ловушками [2]. В настоящее время модель Френкеля используется для интерпретации переноса заряда в широком классе диэлектриков, в том числе в оксидах и нитриде кремния [3], в диэлектриках с высокой диэлектрической проницаемостью [4], в ферроэлектриках [5]. Ионизация ловушек в полупроводниках описывается эффектом Френкеля только в слабых электрических полях (<103 В^м) [6]. В сильных электрических полях ионизация ловушек в полупроводниках объясняется с помощью многофононного механизма [7, 8].
Общим свойством аморфных полупроводников и диэлектриков является возможность локализовать электроны и дырки [9]. Аморфный нитрид кремния SiNx (х < 4/3) обладает способностью локализовать электроны и дырки c гигантским временем удержания более 10 лет при 400 K [10]. Эффект памяти в SiNx широко используется в энергонезависимых микросхемах памяти. Перенос заряда в Кремний-Оксид-Нитрид-Оксид-Кремний (КОНОК)-устройствах памяти определяется механизмом ионизации глубоких центров в SiNx. В настоящее время принято считать, что механизм ионизации в SiNx обусловлен эффектом Френкеля. Эффект Френкеля основан на понижении кулоновского потенциала в сильном электрическом поле [2]. Целью настоящей работы является экспериментальное изучение переноса заряда в SiNx в широком диапазоне электрических полей и температур, количественное сравнение экспериментальных данных с предсказаниями модели Френкеля и квантовой модели многофононного механизма ионизации в SiNx.
1 Настоящая работа поддержана интеграционным проектом № 116 Сибирского отделения Академии наук и национальной программой Корейского министерства науки и технологии по наноэлектронике терабитных масштабов.
This work was supported by project No116 of Siberian Branch of Russian Academy of Sciences, and by National Program for Tera-Level Nanodevice of the Korea Ministry of Science and Technology as one of the 21st Century Frontier Programs.
147
Образцы. Образцы Металл-Нитрид-Оксид-Полупроводник (МНОП)-структур изготавливались на кремниевой подложке п-типа с сопротивлением 7,5 О-см, выращенной методом Чохральского. Тонкий, туннельно-прозрачный окисел толщиной 1,8 нм выращивался при температуре 750 °С. Нитрид кремния толщиной 53 нм был получен при температуре 760 °С методом осаждения в реакторе пониженного давления. Соотношение SiH2Cl2/NH3 составляло 0,1. Алю-
3 2
миниевые электроды площадью 5-10" см напылялись через маску.
Используемая в расчетах модель. На рис. 1 показана энергетическая диаграмма МНОП-структуры для двухзонной модели, учитывающей перенос электронов и дырок в Si3N4 при разных полярностях потенциала [11-13]. Эта модель учитывает инжекцию электронов из отрицательно смещенного электрода и ин-жекцию дырок из положительно смещенного электрода. Для расчета инжекци-онного тока на контакте использовался модифицированный механизм Фаулера—Нордгейма. Модель учитывает рекомбинацию свободных электронов с локализованными дырками и свободных дырок с локализованными электронами. В работе бралось одинаковое сечение рекомбинации для электронов и дырок
13 2
ог = 5 х 10 ст [11]. Перенос заряда описывался с помощью уравнений Шокли—Рида—Холла и уравнения Пуассона, учитывающего неоднородное распределение электрического поля в нитриде кремния.
дп{х, ,) = ^ д(х’ ,) -а ш(х, ,)( - (х, t)) + п((х, t)Р(п)(х, ,) -огип(х, t) д (х, t); (1)
д, е дх г
дпА£t) = оип( х, t) (, - п (х, t)) - п (х, t) Р(п) (х, t) - огир( х, t )п, (х, t); (2)
- аур(х, t) ( - Р, (х, t)) + Р, (х, t)Р(р) (х, t) - агур(х, ,)п, (х, ,); (3 )
^t) = оир(х, t) ( - р,(х, t))- )(х, t)Р(р)(х, t) -огир(х, t)п,(х, t); (4)
д^ = е (п, (х, 11) - Р, (X 11)) (5)
— = -е------------------. (5)
дх 880
Здесь п и п, — концентрации свободных и захваченных электронов, р и р, — концентрации свободных и захваченных дырок, Ые( и — концентрации электронных и дырочных ловушек соответственно, ^ (х,,) — локальное электрическое поле, е — заряд электрона, и — дрейфовая скорость, а = 7,5 — низкочастотная диэлектрическая проницаемость SiзN4. В работе для электронов и дырок использовалось одинаковое сечение захвата о = 5 х 10" см . Дрейфовая скорость электронов связана с плотностью тока соотношением у = епи, а дырок — соотношением у = -ери. Р^п’рР — вероятности ионизации ловушки в Si3N4 в единицу времени.
Модель Френкеля рассматривалась с учетом термически облегченного туннелирования (ТОТ).
Рис. 1. Энергетические диаграммы МНОП-структуры для двухзонной модели переноса заряда в а — при положительном потенциале на алюминии,
б — при отрицательном потенциале на алюминии
Вероятность ионизации в модели Френкеля определяется как
(6)
Здесь Ш, — энергия ловушки, / — постоянная Френкеля, аш = 4,0 — высокочастотная диэлектрическая проницаемость, V — частотный фактор. Вероятность ионизации для ТОТ определяется как
Р =
V
~кТ
щ -р4¥
Г с1аехр - —-----------------------Гб/хд/2т {е¥(х)- а)
* 1гТ Р) *
(7)
Здесь У(х) — кулоновский потенциал в электрическом поле:
V (х) = Ж, -
■ - Гх, а — энергия, соответствующая возбужденному уров-
4пе1е0 х
ню (рис. 2), т* — эффективная туннельная масса. Классические точки поворота х1, х2 даются выражением
х1,2 =
еГ
(
1 +
1 --
еГ
\
1/2 Л
паап
-а)2
(8)
Туннелирование электрона в этой модели рассматривается в квазикласси-ческом приближении и сводится к вычислению эллиптического интеграла.
Другая модель основывается на многофононном механизме ионизации ловушки. В этой модели предполагается, что ловушка для электрона или дырки
3
е
х
0
е
является нейтральной. В качестве ловушки выступает «осциллятор», или «ядро», встроенное в решетку нитрида. Ловушка может притягивать электрон (электронная ловушка) или дырку (дырочная ловушка). Эта модель основана на квантовой теории Макрам—Эбеи и Ланну [7]. Для описания ловушки вводятся: энергия фонона Жф,, термическая Жт и оптическая энергия (рис. 3). Внеш-
нее электрическое поле может способствовать процессу ионизации. В ре-
Рис. 2. Ионизация кулоновской ловушки в сильных электрических полях. Здесь: прямое туннелирование при низких температурах и высоких полях, термически облегченное туннелирование при умеренных температурах, термическая ионизация (эмиссия Френкеля) в слабых полях и высоких температурах
Рис. 3. Конфигурационная диаграмма нейтральной ловушки для модели многофононной ионизации
зультате ионизации получается пустое «ядро» и свободный носитель, суммарная энергия которых равна начальной энергии заполненной ловушки. Как правило, после ионизации ловушки энергия «ядра» соответствует возбужденному состоянию. Избыток энергии расходуется на колебательную моду решетки. Согласно квантовым расчетам [7], для скорости ионизации получается следующее выражение:
пЖи Жр1г
----^ 5 соШ——
2кТ 2кТ
(
5
\
81пЬ(Жрй /2кт)
Р (Жт + пЖрк)
(9)
Р (Ж) =
вГ
2^2ш * Ж
ехр
3/2 ^
3 НвЕ
5 =
Жор{ - Жт
Ж,
рИ
Здесь 1п — модифицированная функция Бесселя, а Рг(Ж) описывает вероятность туннелирования носителя через треугольный барьер с высотой Ж.
п
п=-го
Сравнение эксперимента с расчетом. Экспериментальная часть работы включала две серии измерений, в одной из которых измерялись температурные зависимости тока для трех напряжений J-T (рис. 4) в другой — вольт-амперные характеристики J-V для трех температур (рис. 5). Вольт-амперные измерения проводились при двух полярностях потенциала на алюминиевом электроде. Температурные зависимости (см. рис. 4) представлены в аррениусовских координатах (log J - T-1).
о < -D
О
С
20
30
40
U(V)
50
60
1000/Т, к
Рис. 4. Зависимости тока от температуры для МНОП-структуры при различных прикладываемых положительных потенциалах на алюминии. Экспериментальные результаты показаны точками, расчет с использованием модели Френкеля показан пунктирной линией. На рисунке приведены значения эффективной туннельной массы для каждого напряжения. Расчет с использованием теории многофо-нонной ионизации показан сплошной линией
Рис. 5. Вольт-амперные характеристики МНОП-структуры для разных температур. Измерения проводились при двух полярностях на алюминиевом электроде. Эксперимент показан точками, а расчет с использованием теории многофононной ионизации — сплошной линией. Знаки (+) и (—) соответствуют положительной и отрицательной полярностям прикладываемого напряжения на алюминии
Вначале для описания экспериментальных данных (log J - T-1) использовалась модель Френкеля с учетом ТОТ. Наилучшее согласие эксперимента с расчетом получено для равных энергий (Wt = 1,2 эВ), эффективных масс и концентраций (N t= 7 х 1019 см3) электронных и дырочных ловушек. Согласия эксперимента с расчетом удалось достичь только с использованием нефизично малого частотного фактора v = 6-106 с-1. В оригинальной работе Френкеля величина частотного фактора оценивалась как v&Wt/h « 1015с-1 [2]. Ранее в работах [12, 13] для согласия с экспериментом также использовался заниженный частотный фактор v= 106 - 109 с-1. В области низких температур доминирует туннельный эффект, и основную роль играет эффективная туннельная масса. Для достижения наилучшего согласия эксперимента с расчетом для разных напряжений не-
обходимо было использовать разные эффективные туннельные массы: для V = 44 В — т* = 4,8 те, для V = 37 В — т* = 3,8 те, и для V = 30 В — т* = 2,0 те. Используемые в расчетах величины эффективных масс на порядок больше величины эффективной туннельной массы 0,3-0,6 те электрона и дырки, определенной в Б1КХ из эксперимента [14-16]. Отсюда мы сделали вывод, что модель Френкеля формально описывает перенос заряда в Б1КХ, но при этом необходимо использовать в расчетах нефизично малый частотный фактор и аномально большую эффективную туннельную массу.
Далее для описания эксперимента была использована теория многофо-нонной ионизации. В этой теории варьировались следующие параметры: ЖорЬ ЖрИ. Эффективная туннельная масса для электрона и дырки везде принималась равной те*=ть*=0,5т0. Наилучшее согласие расчета с экспериментом получено при следующих параметрах: Жор1 = 2,8 эВ, Жт = 1,4 эВ, Жри = 60 мэВ (см. рис. 4).
Для моделирования вольт-амперных характеристик (рис. 5) использовалась только теория многофононной ионизации. Параметры для расчета брались те же, что и в моделировании температурных зависимостей тока. В эксперименте ток при отрицательной полярности потенциала на алюминии в три раза больше тока при положительной полярности. Эта же особенность наблюдалась при моделировании. Такое различие тока при разной полярности на металле объясняется разными условиями инжекции электронов (при положительной полярности на алюминии) и дырок (при отрицательной полярности) в нитрид через границу 81/8Ю2. Эффективный барьер для дырок на границе 81/Б1КХ ниже, чем для электронов [15]. Поэтому инжекционный ток дырок на границе 81/Б1КХ больше инжекционного тока электронов при той же величине поля в туннельнотонком оксиде.
Наилучшее согласие расчета с экспериментом при использовании теории многофононной ионизации было получено при одинаковых параметрах электронных и дырочных ловушек, таких как концентрация, сечение захвата, оптическая и термическая энергия. Раннее предположение об идентичности параметров электронных и дырочных ловушек было сделано в работе [12] для объяснения экспериментальных результатов по стеканию заряда при различных «стягивающих» напряжениях.
Величина оптической энергии в расчетах получилась в два раза больше величины термической энергии Жор(~2Жт. Это означает, что энергия активации при захвате электрона или дырки близка к нулю [8]. Этот факт согласуется со слабой температурной зависимостью сечения захвата для электронов и дырок в Б1КХ, наблюдаемой экспериментально [17].
Энергия фонона Жри =60 мэВ, полученная в расчетах, совпадает с энергией фонона в аморфном кремнии, величина которой определена по Рама-новским спектрам в работах [18, 19]. Это говорит о том, что в качестве электронных и дырочных ловушек могут выступать аморфные кремниевые нанокластеры в нитриде кремния. Такая гипотеза высказывалась в работах [10, 20]. Квантово-размерные эффекты в аморфных кремниевых точках, встроенных в нитрид кремния, наблюдались в экспериментах по оптическому поглощению и фотолюминесценции [21, 22]. Недавно было показано, что
аморфные Si нанокластеры могут захватывать электроны и дырки в нитриде кремния [23].
* * *
Таким образом, J-V и J-T экспериментальные данные для SiNx формально могут быть описаны моделью Френкеля с учетом двухзонного механизма проводимости. Интерпретация экспериментальных результатов с использованием модели Френкеля даёт нефизично малое значение частотного фактора и аномально большую эффективную туннельную массу. С другой стороны, перенос заряда в нитриде кремния количественно описывается с использованием двухзонной модели и многофононного механизма ионизации ловушек. Наилучшее согласие между экспериментом и расчетом получается при одинаковых параметрах для электронных и дырочных ловушек, таких как концентрация, сечение захвата, туннельная масса, оптическая и термическая энергия. Величина энергии фонона в расчете совпадает с энергией фонона в аморфном кремнии. Этот факт свидетельствует в пользу того, что электронными и дырочными ловушками в нитриде кремния являются кремниевые нанокластеры.
БИБЛИОГРАФИЧЕСКИЕ ССЫЛКИ
1. Pool H. H. On the Dielectric Constant and Electrical Conductivity of Mica in Intense Fields // Phil. Mag. 1916. V. 34. № 187. Р. 112-128.
2. Френкель Я. И. К теории электрического пробоя в диэлектриках и электронных полупроводниках // ЖЭТФ. 1938. Т. 8. № 12. С. 1292-1301.
3. Sze S. M. J. Appl. Phys. 1967. V. 38. № 7. Р. 2951-2956.
4. Zhu W. J., Ma T.-P., Tamagawa T., Kim J., Di Y. Current transport in metal/hafnium oxide/silicon structure» // IEEE Electron Device Letteres. 2002. V. 23. P. 97-100.
5. Laha A., Krupanidhi S. B. Leakage current conduction of pulsed excimer laser ablated BaBi2Nb2O9 thin films // Appl. Phys. 2002. V. 92. № 1. Р. 415-420.
6. Ganichev S. D., Ziemann E., Prettl W., Yassievich I. N., Istratov A. A., Weber E. R. Distinction between the Poole-Frenkel and tunneling models of electric-field-stimulated carrier emission from deep levels in semiconductors // Phys. Rev. B. 2000. V. 61. Р. 10361-10365.
7. Makram-Ebeid S. S., Lannoo M. Quantum model for phonon-assisted tunnel ionization of deep levels in a semiconductor // Phys. Rev. B. 1982. V. 25. Р. 6406-6424.
8. Абакумов В. Н., Перель В. И., Ясиевич И. Н. Безызлучательная рекомбинация в полупроводниках. СПб., 1997. С. 375.
9. Мотт Н., Дэвис Э. Электронные процессы в некристаллических веществах. М., 1982. С. 662.
10. Gritsenko V. A. Electronic Structure and Optical Properties of Silicon Nitride, in Silicon Nitride in Electronics. New York, 1986. Р. 236.
11. Гадияк Г. В., Обрехт М. A, Синица С. П. Расчет биполярной инжекции и рекомбинации в МНОП-структурах // Микроэлектроника. 1985. Т. 14. С. 512-516.
12. Гриценко В. A., Меерсон E. E., Травков И. В., Голтвянский Ю. В. Нестационарный перенос электронов и дырок в режиме деполяризации МНОП-структур: эксперимент и численное моделирование // Микроэлектроника. 1987. Т. 16. № 1. С. 42-50.
13. Bachhofer H., Reisinger H., Bertagnolli E., Philipsborn H. Transient conduction in multidielectric silicon-oxide-nitride-oxide semiconductor structures // Appl. Phys. 2001. V. 89. № 5. Р. 2791-2800.
14. Miyazaki S., Ihara Y., Hirose M. Resonant tunneling through amorphous silicon— silicon nitride double-barrier structures // Phys. Rev. Lett. 1987. V. 59. Р. 125-127.
15. Gritsenko V. A., Meerson E. E., Morokov Yu. N. Thermally assisted hole tunneling at the Au-Si3N4 interface and the energy-band diagram of metal-nitride-oxide-semiconductor structures // Phys. Rev. B. 1989. V. 57. Р. R2081-R2083.
16. Yeo Y. C., Lu Q., Lee W. C., King T. -J., Hu C., WangX., Guo X., Ma T. P. Direct tunneling gate leakage current in transistors with ultrathin silicon nitride gate dielectric // IEEE Electron Dev. Lett. 2000. V. 21. Р. 540-542.
17. Hampton F. L., Cricci J. R. Steady-state electron and hole space charge distribution in LPCVD silicon nitride films // Appl. Phys. Lett. 1979. V. 35. Р. 802-804.
18. Smith J. E., Jr., Brodsky M. H., Crowder B. L., Nathan M. I., Pinczuk A. Raman Spectra of Amorphous Si and Related Tetrahedrally Bonded Semiconductors // Phys. Rev. Lett. 1971. V. 26. Р. 642-646.
19. Giorgis F. Optical absorption and photoluminescence properties of a-Si1-xNx:H films deposited by plasma-enhanced CVD // Phys. Rev. B. 2000. V. 61. Р. 4693-4698.
20. Гриценко В. А. Строение и электронная структура аморфных диэлектриков в кремниевых МДП структурах. Новосибирск, 1993. С. 340.
21. Park N.-M., Choi C.-J., Seong T.-Y., Park S.-J. Quantum Confinement in Amorphous Silicon Quantum Dots Embedded in Silicon Nitride // Phys. Rev. Lett. 2001. V. 86. Р. 1355-1357.
22. Park N.-M., Kim T.-S., Park S.-J. Band gap engineering of amorphous silicon quantum dots for light-emitting diodes // Appl. Phys. Lett. 2001. V. 78. Р. 2575-2577.
23. Park N.-M., Choi S.-H., Park S.-J. Electron charging and discharging in amorphous silicon quantum dots embedded in silicon nitride // Appl. Phys. Lett. 2002. V. 81. Р. 1092-1094.
K. Nasyrov, V. Gritsenko, Yu. Novikov, D. Gritsenko, D.-V. Lee, Ch. Kim
THE CHARGE TRANSFER MECHANISM IN SINX
The Charge Transfer Mechanism in the amorphous SiNx is described. The
experimental data are qualitatively described by the mechanism of ionization
of traps.
