Оптические и электрические параметры гетероструктур в физике полупроводниковых материалов: моделирование и эксперимент
12 1 1 А.В. Благин , Л.В. Благина , И.Г. Попова , В.В. Подольцев
1 Донской государственный технический университет, Ростов-на-Дону 2Институт технологий ДГТУ, г. Волгодонск
Аннотация: В работе содержится анализ результатов моделирования электрофизических параметров светоизлучательных твердых растворов GaInAsSb/GaSb. В работе построены графики зависимостей мощности излучения от силы тока и вольт-амперная характеристика (ВАХ) светодиодного устройства. Массив носителей описан статистикой Ферми. В явном виде раскрыты индикаторы моделирования, основанные на особенностях зонной структуры полупроводниковых систем. Обсуждены результаты экспериментов, показавшие удовлетворительное согласие с данными, полученными на основе расчетов. Ключевые слова: твердые растворы, Sim Windows 1,5, излучающие структуры, токи, ограниченные пространственным зарядом, электрофизические параметры.
Введение
Многокомпонентные гетероструктуры, созданные на основе соединений III-V, II-VI и IV-VI групп Периодической системы вот уже более трёх десятилетий широко используются в электронной технике. Растёт номенклатура полупроводниковых материалов, имеющих прикладное значение.
Это обусловлено с одной стороны, прогрессом полупроводниковых технологий (возникновение и развитие ряда новых методов и модификация известных [1-3]). С другой стороны - возможностью комбинировать различные элементы указанных групп Периодической системы.
Следует отметить, что в последние 20 лет возникли возможности получения многокомпонентных гетероструктур термодинамически нестабильных в условиях кристаллизации [4]. В этой ситуации особую роль играет предварительное моделирование и имитация функций приборного материала, что можно эффективно осуществлять с помощью некоммерческой программы Sim Windows 1,5 (freeware software, автор-разработчик - Д.
Уинстон, Центр оптоэлектронных систем при Колорадском университете, Баулдер, США) [5].
В настоящей работе обсуждаются результаты моделирования и их сопоставление с экспериментальными значениями оптических и электрофизических параметров материалов электронной техники.
Прогнозирование характеристик многокомпонентных гетероструктур с использованием программной среды Sim Windows
Sim Windows - оболочка, позволяющая моделировать свойства гетероструктур в одномерном приближении. Уровни энергии в среде рассчитываются относительно дна зоны проводимости; учитываются дрейф носителей заряда, ток термоэлектронной эмиссии и туннельный ток через поверхность раздела. Процессы рекомбинации в модели включают в себя механизм спонтанной безызлучательной рекомбинации. То есть, Sim Windows рассматривает традиционные модели, что позволяет достичь большей точности.
Ядром программы является симулятор полупроводниковых устройств, который генерирует решение уравнений, статистически описывающих перенос электронов в гетероструктуре как в равновесных, так и в неравновесных условиях.
Базовые уравнения имеют вид [5]:
где и т^ - потенциалы Планка для электронов и дырок; Ес и ¿^ - уровни, соответствующие дну зоны проводимости и потолку валентной зоны
J
соответственно; и Тр - температуры электронного и дырочного газа; ЕРк и 5 - уровни Ферми для электронов и дырок.
Общее уравнение эффективной плотности состояний запишем в виде:
кОг} = ВД^ ; = Ш ■
где Р^р - интеграл Ферми; --
И р(х) = q- nf» -Ь Mü + (х) - NA*(х),
где N№ - эффективная плотность состояний; NDA - концентрация примесей
(доноров и акцепторов); Е0А - энергия соответствующих примесных
уровней; - параметр вырождения примесей; и - концентрация
электронов; р - концентрация дырок; вь„р - эффективная масса носителей.
Неравновесное состояние задаётся внешним полем. В модели, используемой в Sim Windows, это внешнее поле задаётся уравнением Пуассона:
где D - модуль электростатической индукции (смещения), и уравнениями потенциалов Планка для электронов:
J
и дырок
К
рек
где /„ и Л, плотность электронного и дырочного тока;
i.',.^. - скорость
рекомбинации электронов и дырок.
Ранее отмечалось, что в Sim Windows возможен учёт трёх типов тока -дрейфа, термоэлектронной эмиссии и туннелирования через потенциальные барьеры.
Уравнения модели обсуждаются в [6]. Носители с малой энергией описываются в рамках классической статистики, с высокой - распределением Ферми-Дирака.
Модель оболочки оперирует массивами внешних и внутренних фотонов (генерация которых осуществляется внешними устройствами и внутри устройства, например, светоизлучающего диода). Устройство моделируется дискретными элементами в виде двумерных квадратных матриц. Схема дискретизации приведена на рисунке 1.
Граница устройства
Область i
Область i+1
Граница устройства
Падающая во лна
Отражённая волна
7* :l +
ЕГ
Рис. 1. - Дискретизация оптоэлектронного устройства
Уравнение, описывающее волну в ¿-ой области, запишем как
Магнитная компонента имеет вид: с+ р -
н
■eL
где Щ = [й] — i ; Zt = —; к - волновое число; Ei
электрическая
составляющая поля в /-й области; ^ - относительный комплексный
импеданс; •l}i - показатель преломления /-й области; а - показатель
поглощения; [т;] - действительная часть ^
Матрица распространения поля в области (/+7) задаётся параметрами /-й области:
После ряда преобразований выражения - вычислители параметров волны в устройстве принимают вид:
ji+i
-г+1
2Z
(X + Zt+Jet1**^ <:,rf+i tX - Zi+J) ff" 1
Sim Windows допускает иерархическую структуру материалов устройства, так, в роли материалов могут выступать слои многокомпонентных твёрдых растворов. Их характеристики (такие, как
коэффициент преломления, ширина запрещённой зоны) заносятся в файл параметров материала с возможностью изменения и дополнения.
Особый интерес как с практической точки зрения, так и в отношении анализа интересных эффектов, отслеживаемых в процессе моделирования, представляют излучающие структуры с вертикальным резонатором, состоящим из брэгговских отражателей. Их разделяет рабочая область толщиной порядка 50 нм. Схема лазерного светодиода с вертикальным резонатором представлена на рисунке 2.
//////////////////////////////////////////////////////////
В
С "
О Г> Г> 5
Е
а \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
А
Б
Б
Рис. 2. - Схема лазерного светодиода с вертикальным резонатором А - контактный слой ааБЬ(р+); В, Е - брэгговские отражатели; С - слой аа1п8ЬЛБ с мезоструктурой; Б - рабочая область (аа1пБЬ<В1>); Б -контактный слой ааБЬ(п-); а - подложка ааБЬ.
Результаты моделирования и их сравнение с экспериментом
Нами была рассчитана гетероструктура GaInAsSb/GaSb в среде Sim Windows. Построены графики зависимостей мощности излучения от силы тока и вольт-амперная характеристика (ВАХ) светодиодного устройства. Массив носителей описывается статистикой Ферми. На рисунке 3, а) приведена ВАХ устройства и несколько экспериментальных точек. На рисунке 3, б) дана зависимость мощности излучения от силы тока.
A
n
<D
О
<D
o
n
0
1
Я
i di
cd &
10
0,
2,
Voltage, V
(a)
20
10
15
Current, mA (б)
20
25
Рис. 3. - Зависимость электрических параметров лазерного светодиода на основе гетероструктуры GaInAsSb/GaSb
1
а) - вольт-амперная характеристика; б) - мощность излучения; Сплошные кривые - результаты моделирования Sim Windows; Значки ■ , • -экспериментальные данные.
На ВАХ получено более удовлетворительное соответствие модели эксперименту. Экспериментально измеренная мощность излучения (тестер F0D1203A/B, Россия) была ниже спрогнозированной, причём, AP Ф const, что требует проведения дополнительных исследований. Возможно, реальный уровень инжекции ниже расчётного. Однако расхождение в 3-5 мВт не является существенным.
Гетероструктуры GaInAsSb/GaSb получены методом градиентной жидкофазной эпитаксии [3]. Физические свойства базовых твердых растворов, образующих данный материал, описаны в работах [7,8]. Технологические аспекты приведены в [9]. Был проведен также анализ характеристик р-и-переходов. Наблюдалось существенное снижение квантового выхода люминесценции. Нами были проведены исследования причин, вызывающих деградацию люминесцентных свойств, что становится все более актуальной задачей [10].
Вольт-амперные характеристики р-и-переходов хорошо
аппроксимируются выражением:
/
I = Is exp
С тт Л
eU
.вт.
В диапазоне смещений 1,35.. .1,72 В значение в =2. Это свидетельствует о том, что основным процессом является рекомбинация в области пространственного заряда. При более высоких напряжениях, если вычесть из напряжения падение потенциала на последовательном сопротивлении р-п-перехода, наблюдается участок с в=1, что указывает на процесс рекомбинации в нейтральной области. У термообработанных диодов фактор
неидеальности значительно больше двух. Вольт-амперные характеристики этих диодов, перестроенные в двойном логарифмическом масштабе (рис. 4, а), имеют участки, характерные для протекания токов, ограниченных пространственным зарядом (ТОПЗ).
В эксперименте (рис.4, б) ток обратно пропорционален кубу толщины пленки, что согласуется с теорией и свидетельствует о возможности использования ТОПЗ для интерпретации ВАХ. Таким образом, в процессе термообработки вблизи активной зоны возникает скомпенсированная область с высоким сопротивлением.
Рис. 4. - Электрофизические характеристики светодиодных пленок
GaInAsSb/GaSb
а) - ВАХ пленок, 1 - без отжига; 2 - с отжигом при Т= 713 К; 3 - 753 К; 4 -890 К; 5 - 923 К; б) - зависимость тока от толщины высокоомной области.
Оценка содержания рекомбинационных центров в высокоомной области была выполнена двумя способами. На вольт-амперных характеристиках (рис. 1, а) наблюдаются участки резкого возрастания тока, связанные с заполнением ловушек.
:
Потенциал, при котором происходит заполнение ловушек, позволяет оценить их концентрацию:
Uv 8 8^
N =
ed2
где Ux - потенциал компенсации ловушек; d - ширина высокоомной зоны. Концентрация ловушек, определенная по этой формуле, приведена в таблице.
Таблица 1
Концентрация ловушек в р - «-переходах термообработанных эпитаксиальных пленок ОаЛзР<В1>
Температура обработки, К Концентрация, определенная по ВАХ, см-3 Концентрация, определенная из емкостных измерений, см-3 Относительная яркость
Не подвергались 8,5-1013 1,5-1014 1,0
термообработке
713 3,6-1014 2,2-1014 0,85
753 5,4-1014 4,3-1014 0,08
890 7,9-1014 5,5-1014 0,06
923 9,4-1014 5,7-1014 0,03
Из экспериментов следует, что при термообработке пленок определяющую роль в электрофизических параметрах играют токи, ограниченные пространственным зарядом. Концентрация ловушек существенно зависит от структурного совершенства пленок.
Метод моделирования можно считать перспективным для прогнозирования электрических и оптических характеристик новых материалов.
Заключение
В статье обсуждается методика компьютерного моделирования сложных полупроводниковых систем в среде Sim Windows. Методика отработана на четырёхкомпонентных светоизлучающих гетероструктурах GaInAsSb/GaSb.
В явном виде раскрыты индикаторы моделирования, основанные на особенностях зонной структуры полупроводниковых систем. Обсуждены результаты экспериментов, показавшие удовлетворительное согласие с данными, полученными на основе расчетов.
Литература
1. Алфимова Д.Л., Благин А.В., Гусев Д.А., Лунина М.А., Сысоев И.А. Формирование массивов квантовых точек Оа^щ.^АБуР^у/ОаАз в
процессе ионно-лучевого осаждения // Неорганические Материалы, 2014. Т.50. №3. С. 237-243.
2. Patrick McCray. MBE deserves a place in the history books // Nature, Vol.2. 2007. рр. 259.
3. Лунина Л.С., Благин А.В., Алфимова Д.Л., Попов А.И., Разумовский П.И. Физика градиентной эпитаксии полупроводниковых гетероструктур. Ростов-на-Дону. Изд-во СКНЦ ВШ, 2008. 212 с.
4. Середин П. В. Спинодальный распад в эпитаксиальных твёрдых растворах гетероструктур Al^ga^^As/Cи Gct^frti.^P/Getjils (100). Самарского центра РАН. Т.11. №3, 2009. С. 46-52.
5. Winston D.W. Faculty of the Graduate School of the University of Colorado, 1997. Physical Simulation of Optoelectronic Semiconductor Devices. URL: wemif.pwr.wroc.pe/zpp/laboratoria/labo_opto/labo_opto.php.
6. Yang K., East J.R., Haddad G.I. Numerical Modeling of Abrupt Heterojunctions Using a Thermoionic Field Emission Bundary Condition// Solid State Electron. Vol. 36. 1994. рр.321-330.
7. Благин А.В., Благина Л.В., Нефедова Н.А., Попова И.Г., Политов О.А. Физические свойства светоизлучательных твердых растворов, изопериодных антимониду галлия // Инженерный вестник Дона, 2017. №4. URL: ivdon.ru/ru/magazine/archive/n4y2017/4613.
8. Благин А.В., Благина Л.В., Кодин В.В., Нефедова Н.А., Попова И.Г. Физические свойства многокомпонентных узкозонных твердых растворов с мезоструктурой // Инженерный вестник Дона, 2018. №2. URL: ivdon.ru/ru/magazine/archive/N2y2018/4830
9. Благин А.В., Благина Л.В., Кодин В.В., Нефедова Н.А., Попова И.Г. Структура эпитаксиальных слоев узкозонных твердых растворов и компенсация дефектов // Инженерный вестник Дона, 2018. №2. URL: ivdon.ru/ru/magazine/archive/N2y2018/4820
10. Никифоров С. Проблемы, теория и реальность светодиодов для современных систем отображения информации высокого качества // Компоненты и технологии, 2005. №5. URL: kit-e.ru/articles/led/2005_5_48.php.
References
1. Alfimova D.L., Blagin A.V., Gusev D.A., Lunina M.A., Sysoev I.A. /Neorganicheskie Materialy, 2014. Vol.50. №3. рр. 237-243.
2. Patrick McCray. Nature. Vol.2.2007. Р. 259.
3. Lunina L.S., Blagin A.V., Alfimova D.L., Popov A.I., Razumovskij P.I. Fizika gradientnoj epitaksii poluprovodnikovykh geterostruktur. [Gradient epitaxy physics of semiconductor heterostructures]. Rostov-na-Donu. Izd-vo SKNC VSH, 2008. 212 p.
4. Seredin P.V. Samarskogo centra RAN, Vol.11. №3. 2009. рр. 46-52.
5. Winston D.W. Faculty of the Graduate School of the University of Colorado, 1997. Physical Simulation of Optoelectronic Semiconductor Devices. URL: wemif.pwr.wroc.pe/zpp/laboratoria/labo_opto/labo_opto.php.
6. Yang К., East 1.Я., Haddad а1. Solid State Electron. Уо1. 36. 1994. рр. 321-330.
7. В^т А.У., В^т Ь.У., Nefedova КА., Popova Ш., Po1itov О.А. Inzenernyj vestnik Dona (Rus), 2017. №4. ЦКЬ: ivdon.ru/ru/magazine/archive/n4y2017/4613.
8. В^т А.У., B1agina Ь.У., КоШд У.У., Nefedova N.A., Popova Ш. Inzenernyj vestnik Dona ^ш), 2018. №2. URL: ivdon.ru/ru/magazine/archive/N2y2018/4830
9. В^т А.У., B1agina Ь.У., Kodin У.У., Nefedova N.A., Popova Ш. Inzenernyj vestnik Dona (Rus), 2018. № 2. URL: ivdon.ru/ru/magazine/archive/N2y2018/4820.
10. Nikiforov S. Komponenty i tekhno1ogii, 2005. № 5. URL: кк-e.ru/artic1es/1ed/2005_5_48.php.