Студопедия КАТЕГОРИИ: АвтоАвтоматизацияАрхитектураАстрономияАудитБиологияБухгалтерияВоенное делоГенетикаГеографияГеологияГосударствоДомЖурналистика и СМИИзобретательствоИностранные языкиИнформатикаИскусствоИсторияКомпьютерыКулинарияКультураЛексикологияЛитератураЛогикаМаркетингМатематикаМашиностроениеМедицинаМенеджментМеталлы и СваркаМеханикаМузыкаНаселениеОбразованиеОхрана безопасности жизниОхрана ТрудаПедагогикаПолитикаПравоПриборостроениеПрограммированиеПроизводствоПромышленностьПсихологияРадиоРегилияСвязьСоциологияСпортСтандартизацияСтроительствоТехнологииТорговляТуризмФизикаФизиологияФилософияФинансыХимияХозяйствоЦеннообразованиеЧерчениеЭкологияЭконометрикаЭкономикаЭлектроникаЮриспунденкция |
Исследуемая гетероструктура
Интересующие нас структуры ранее исследовались не только теоретически, но и экспериментально в работах [1,2,new2/ new1, new3]. Они были выращены в Институте Фердинанда Брауна (Берлин) методом металлоорганической парофазной эпитаксии на ориентированных в плоскости (100) подложках из легированного кремнием GaAs (Рис.2.1). Сначала на подложку наносились слои AlxGaAs1-x n-типа с плавно меняющейся концентрацией алюминия и различным уровнем легирования. Далее для создания квантовой ямы размещался слой GaAsyP1-y, причем из-за различия постоянных решетки слоев GaAsyP1-y и слои AlxGa1-xAs квантовая яма получалась биаксиально растянутой на 0,58%. Затем размещались слои AlxGa1-xAs только уже p-типа. Состав каждого слоя (представлен на Рис.11), а также их толщины и уровень легирования (на Рис.11 по просьбе изготовителя не представлены) использованы в качестве вводимых параметров при расчёте. P-GaAsp-Al0.7-0Ga0.3-1Asp-Al0.7Ga0.3Asp-Al0.45Ga0.55Asp-Al0.45Ga0.55Asp-Al0.45Ga0.55Asp-Al0.45Ga0.55Asp-Al0.3-0.45Ga0.7-0.55AsGaAs0.84P0.16n- Al0.3-0.45Ga0.7-0.55Asn-Al0.45Ga0.55Asn-Al0.45Ga0.55Asn-Al0.45Ga0.55Asn-Al0.45Ga0.55Asn-Al0.45-0.7Ga0.55-0.3Asn-Al0.7Ga0.3Asn-Al0.7GaAs0.3n-Al0.7-0Ga0.3-1Asn-GaAs
Рис.11 Схематическое изображение структуры
Особенности используемой программы.
Расчеты были выполнены с помощью программы “Heterostructure Design Studio 2.1”, написанной бывшим аспирантом нашей кафедры Колоколовым К.И.. В основе программы заложен конечно-разностный метод, позволяющий самосогласованно решать уравнение Шредингера с к-р гамильтонианом в представлении Латтинжера-Кона с учётом членов, описывающих деформацию, и уравнение Пуассона для электростатического потенциала. Программа позволяет проводить численный расчет зонной структуры, электронных и оптических свойств квантовых ям различной формы и состава, а также определять положение уровней размерного квантования, непосредственно в самой яме, как при нормальных условиях, так и при одноосном сжатии вдоль различных кристаллографических направлений. Помимо этого предусмотрен учет влияния температуры и электрического поля. •Согласно Биру и Пикусу [22] изменение энергии dЕ состояний в зоне проводимости в окрестности точки Г в GaAs, AlAs, GaP, при деформации кристалла, описываемой тензором eij, определяется одним деформационным потенциалом ac = Dxx: dЕ = Dxx(exx +eyy +ezz) , (1) где x, y, z отвечают, соответственно, кристаллографическим направлениям [100], [010] и [001] (в программе “Heterostructure Design Studio 2.1” считается, что две первых координаты отвечают плоскости структуры, а третья “z” – направлению роста структуры). Т. е. деформация приводит к сдвигу энергетического спектра электронов. Воздействие деформации на валентную зону значительно сложнее и требует привлечения трех деформационных потенциалов av, b, d. В представлении Латтинжера Кона [23] гамильтониан 4×4, описывающий состояния легких и тяжелых дырок в окрестности точки Г, при наличии деформации, может быть представлен в виде [24]: , (2) где: , , , kx, ky, kz – компоненты квазиимпульса, отсчитываемого от точки Г, m – масса свободного электрона, gi – параметры Латтинжера. Если одноосная деформация прикладывается вдоль направлений [110], [100] или [001] (к последней сводятся плоско-напряжённые состояния, обусловленные несоответствием решёток), матрица (1) может быть унитарным преобразованием [25] приведена к виду: (3), где: (4), (5). •Для обеспечения самосогласованности расчета необходимо совместное решение уравнения Пуассона, которое позволяет рассчитать пространственный потенциал, связанный с неоднородным распределением заряда, и уравнение Шредингера с гамильтонианом, включающим как этот пространственный потенциал, так и матрицу Н' (3). Для решения такой системы интегро-дифференциальных уравнений обычно используется самосогласованный метод: сначала вычисляются собственные значения и собственные функции с гамильтонианом H’, т.е. без пространственного потенциала. После этого из уравнения Пуассона определяется потенциал уже с учетом вклада от появившегося неоднородного распределения заряда, который снова подставляется в гамильтониан. Эту операцию необходимо проводить до тех пор, пока каждая последующая итерация уже не будет приводить к изменению решения. В общем случае такая схема не всегда работает и процесс самосогласования необязательно сходится. Однако для системы с уравнением Шредингера с гамильтонианом, учитывающим наличие квантовой ямы, такой метод имеет хорошую сходимость и 20 итераций достаточно для устойчивости седьмого значащего разряда в значении энергии, получаемой из нахождения собственных значений гамильтониана. • Расчет отвечающего за учет деформации члена гамильтониана дырок H’ производился по формулам (3), (4) и (5) с использованием литературных данных о значениях деформационных потенциалов и вычисляемых самой программой значений тензора деформаций. Значения тензора деформации при одноосном сжатии вдоль направления [110] находились из следующих соображений. Считалось, что компоненты εijтензора определены в системе координат xyz. Тогдавсистеме координат x’y’z’ , повернутой вокруг оси z относительно xyz на 45°, тензор напряжений имеет вид: (6) где Р-величина одноосного давления. При переходе от одной системы координат к другой, компоненты тензора изменяются по закону: (7) где суммирование идет по дважды повторяющимся индексам, а матрица Аil представляет собой матрицу перехода из одной системы координат в другую. В рассматриваемом случае она имеет вид: . (8) После преобразования тензора, в системе координат xyz получается: . (9) После использования закона Гука, связывающего тензор напряжений и тензор деформации, получаются следующие выражения для компонент тензора εijпри одноосной деформации вдоль направления [110]: , (10) где S11, S12, S44 – модули упругости. В случае одноосной деформации вдоль направления [100] тензор напряжений имеет вид (6) в исходной системе координат xyz. Применяя закон Гука в этом случае, получаем следующие выражения для ненулевых компонент тензора деформации: . (11) Наконец, в случае одноосной деформации вдоль направления [001] тензор напряжений в системе координат xyz имеет вид: . (12) Соответственно, ненулевыми в этом случае являются следующие компоненты тензора деформации: . (13) • При расчетах учитывалось, что температура влияет на ширину запрещенной зоны и вызывает изменение постоянной решетки. Для описания температурной зависимости ширины запрещенной зоны использовалось эмпирическое выражение [27,28]:
где α и β - подгоночные параметры, а eg(0) - ширина запрещенной зоны при нулевой температуре. Изменение постоянной решетки с температурой учитывалось в простейшем линейном приближении: alc(T)=alc(0)+aT(T-300) где aT- линейный коэффициент расширения, . alc(0) – постоянная решетки при комнатной температуре. Все необходимые значения параметров вводились на основании литературных данных [27,28].
Порядок расчета.
Для расчета уровней размерного квантования, квазиуровней Ферми и волновых функций электронов и дырок в условиях одноосного сжатия при использовании программы необходимо было проделать следующие процедуры: (1) ввод необходимых параметров соединений, из которых состоит изучаемая структура; (2) моделирование гетероструктуры; (3) далее указываются внешние параметры (давление, температура); (4) расчёт зонной диаграммы гетероструктуры; (5) наконец после выбора модели расчета из 5-ти возможных (модели: Латтинжера-Кона 4x4+conduction band, Латтинжера-Кона 6x6+conduction band, Латтинжера-Кона 8x8, Латтинжера-Кона 4x4, Латтинжера-Кона 6x6), программа рассчитывает значения уровней размерного квантования, волновые функции, уровни Ферми. •Расчёт начинался с того, что на основании литературных данных [27,28] для GaAs, AlAs, AlxGa1-xAs, GaP вводились значения необходимых для проведения расчётов параметров: постоянная кристаллической решетки и величина энергетической щели, константы упругости, значения деформационных потенциалов и диэлектрической проницаемости, параметры Латтинжера и т.д. (Рис.12). Поскольку для GaAsyP1-y в литературе данные о параметрах отсутствовали, то необходимые величины определялись предусмотренной в программе процедурой линейной экстраполяции между значениями, отвечающими бинарным сплавам GaAs и GaP. Эта же процедура применялась для определения значений некоторых параметров в AlxGa1-xAs по литературным данным для GaAs и AlAs. Чтобы убедиться в разумности результатов экстраполяционной процедуры тем же методом был проведён расчёт значений таких параметров, как энергетическая щель и постоянной решетки, для хорошо представленных в литературе [27] четверных сплавов InGaAsP на основании литературных данных [28] о параметрах InGa и рассчитанных нами значений для GaAsyP1-y.
Рис.14. Таблица вводимых параметров.
• При моделировании структуры было необходимо ввести для подложки и каждого слоя состав, толщину, концентрацию и тип легирующих примесей (в дальнейшем при расчётах примесь считалась полностью ионизованной, что обеспечивало моделирование встроенного поля p-n-перехода (Рис. 13)).
Рис.13. Зонная диаграмма центральной части гетероструктуры
При расчёте конкретной структуры программой возникла следующая проблема: в программе не предусмотрено описание варизонных структур. Поэтому, прилегающие к квантовой яме из GaAs0.84P0.16 слои AlxGa1-xAs с плавно изменяющимся содержанием алюминия от х = 0.3 до x = 0.45 (Рис.11) были разбиты на три независимых слоя одинаковой толщины и имеющие состав x=0.3, x=0.4 x=0.45 (Рис.15).
• Для моделирования реальных условий эксперимента вводится необходимое значение температуры, а также предусмотрен учёт приложенного одноосного давления с указанием его направления. Учёт влияния внешних воздействий описан выше.
• При расчёте зонной диаграммы структуры учитывалось, что из-за значительной разницы постоянной решетки слой GaAs0.84P0.16 находится во внутренне напряженном, а именно биаксиально растянутом состоянии. Полагая, что тонкий слой GaAsP повторяет структуру слоёв AlxGa1-xAs, рассчитывались значения деформации в нём: , (14) где a0 и а - постоянные решетки AlGaAs и GaAsP соответственно (при этом a0 > а), C11 и C12 ,- упругие константы в GaAsP. Затем, используя необходимое из выражений (10), (11) или (13), рассчитывался вклад в деформацию, вызванный внешним воздействием, и, наконец, используя значения деформационных потенциалов определялись значения энергетических щелей между дном зоны проводимости и потолком лёгких и тяжёлых дырок в GaAsP. Также считалось, что скачок потенциала на гетерогранице составляет 60% и 40% от величины изменения энергетической щели для зоны проводимости и валентной зоны соответственно [new8].
• Далее проводился расчёт волновых функций электронов и дырок и соответствующих значений энергии пространственного квантования. Чтобы иметь возможность сравнить результаты с данными работ [1,2,new1,new2,new4], считалось, как и в перечисленных работах, что концентрация неравновесных носителей в КЯ n = p = 2×1012 см-2. Возможности использованной вычислительной техники позволяли провести расчёт за разумное (порядка часа) время на сетке не более чем из 200 точек только при использовании для описания спектра дырок наиболее простой модели - гамильтониана Латтинжера 4x4. В результате при равномерном распределении точек по всей структуре реальная точность воспроизведения вида волновой функции, которая, очевидно, сосредоточена прежде всего в квантовой яме, оказывается недостаточной, чтобы корректно рассчитать пространственный потенциал из уравнения Пуассона, а, значит, гарантировать достаточную точность расчёта энергии пространственного квантования. Поэтому расчёт проводился в два этапа. На первом использовалось равномерное распределение точек, и на результате полученных волновых функций выяснялась реальная область локализации, которая фактически ограничивалась квантовой ямой и прилегающими к ней ближайшими слоями (огибающая волновых функций затухала в 100 раз на расстоянии 18 нм от границы КЯ). На следующем этапе расчёт проводился для сетки, имеющей большую концентрацию точек в значимой части структуры: в квантовой яме GaAs0.84P0.16 и четырёх ближайших к ней слоях AlxGa1-xAs. Следует сказать, что в использованной программе волновые функции электронов и дырок вычисляются в виде разложения по базисным функциям в представлении Латтинжера-Кона с полным угловым моментом J=3/2 и его проекциями mj = ±1/2 и mj = ±3/2 для лёгких и тяжёлых дырок соответственно. При одноосном сжатии в плоскости структуры симметрия падает, и это приводит к перемешиванию базисных функций, описывающих лёгкие и тяжёлые дырки в недеформированном кристалле. Вследствие этого при такой деформации определить связь уровня с лёгкими или тяжёлыми дырками становится непросто. Однако интерфейс используемой программы предусматривает вывод огибающих волновых функций по базисным для каждого из уровней в КЯ. В результате, следуя подходу работы [32], если проинтегрировать квадрат модуля огибающей в значимой области её существования, то можно установить вклад в волновую функцию от каждой из базисных функций и связать природу уровня с лёгкими или тяжёлыми дырками. После расчета спектра и волновых функций использованная программа позволяла ввести схему возможных оптических переходов между системой уровней и рассчитать значения матричных элементов оператора электрон фотонного взаимодействия для этих переходов между состояниями i и j: Aepji, (15) где A = A0e- векторный потенциал электромагнитной волны, e – единичный вектор, задающий поляризацию фотонов, – оператор импульса, e – заряд и m0 – масса электрона, а затем по известной методике [26, new9] и коэффициент оптического усиления g для энергии Eph = çepvc÷2. (16) Здесь суммирование идет по волновому вектору k, E(c) и E(v) – энергии соответствующих состояний в зоне проводимости и в валентной зоне, nr – коэффициент преломления, fc и fv – Фермиевские функции распределения электронов и дырок, квазиуровни Ферми которых определяются концентрацией носителей: ò fc(v).dE = n(p). Коэффициент оптического усиления показывает, как нарастает интенсивность волны I по мере прохождения среды на расстояние d за счет рекомбинации электронов и дырок. По сути это отрицательный коэффициент поглощения a и рассчитывается по той же формуле во втором порядке теории возмущений, о чем говорит квадрат матричного элемента в выражении (16). Только положительный вклад, как видно, дают теперь переходы между занятыми состояниями электронов в зоне проводимости c и свободными состояниями в валентной зоне v, т.е. процессы рекомбинации электронов и дырок с энергиями E(с) - E(v) = |
||
Последнее изменение этой страницы: 2018-05-29; просмотров: 209. stydopedya.ru не претендует на авторское право материалов, которые вылажены, но предоставляет бесплатный доступ к ним. В случае нарушения авторского права или персональных данных напишите сюда... |