УДК 535.374:621.375.8
ГЕТЕРОСТРУКТУРА С РАСШИРЕННЫМ ОПТИЧЕСКИМ ВОЛНОВОДОМ ДЛЯ ДИОДНЫХ ЛАЗЕРОВ
Т. И. Гущик, Н. В. Дьячков, А. П. Богатов, А.Е. Дракин
Осуществлена оптимизация толщин и состава слоев диодного гетеролазера по параметру отношения, оптической ширины пучка d к фактору Г - оптического ограничения. Проверена устойчивость распределения, вертикального поля, в ближней зоне и наиболее важных волно-водних параметров к статическим изменениям, а, также к динамическим вариациям показателей преломления, слоев гетероструктури в рабочем, режиме лазера.
Ключевые слова: полупроводниковые лазеры, гетероструктуры. расширение ближнего вертикального поля, вытекание излучения в эмиттерньте слои.
Основным фактором, ограничивающим яркость многих типов одномодовьтх лазеров, является проблема катастрофического оптического разрушения или деградации (КОД) [1 4]. Эта проблема вызвана тем. что размер пятна излучения типичных одномодовьтх лазеров, как правило, ограничен конструкцией лазера в горизонтальном направлении 3 6 мкм. и параметрами слоев структуры в вертикальном направлении типично 0.5 мкм. и таким образом, при достижении на выходном зеркале лазера критической плотности мощности, называемого порогом КОД. грань разрушается. Очевидно, что увеличение выходной мощности при неизменной плотности оптического пучка (ниже критического значения) возможно за счет увеличения ширины пучка в вертикальном направлении (перпендикулярно слоям структуры). Отметим, что увеличение размера пучка в горизонтальном направлении сопряжено со срывом поперечно-одномодовой генерации [5] и явлением филаментации [6]. В свою очередь, расширение вертикального поля, при сохранении толщины активной области, уменьшает тепловыделение на поверхности зеркала, которое пропорционально плотности оптического потока в активной области и является затравочным механизмом КОД [3, 4]. Исходя из этих соображений, нами был выбран критерий оптимизации структуры по отношению d/Г - оптической ширины пучка d к фактору Г - оптического ограничения.
ФИАН, Россия 119991, Москва, Ленинский пр-т, 53; e-mail: [email protected].
х, мкм Угол, град
Рис. 1: Распределение ближнего (а) и дальнего (б) поля фундаментальной моды и параметры качества для волновода подобранной структуры. На рис. 1 (а) приведен также профиль показателей преломления слоев подобранной структуры, d - ширина распределения ближнего поля основной моды по полувысоте; р - ширина диаграммы направленности по полувысоте; d/Г - отношение ширины ближнего поля к фактору оптического ограничения; Ад - наименьший дефицит усиления ненулевых мод относительно усиления нулевой моды на пороге генерации; а& - относительное уменьшение модового усиления фундаментальной моды за счет туннелирования (вытекания) излучения из волновода в эмиттерные слои.
Основной подход к проблеме расширения оптического пучка связан с ослаблением волновода (уменьшение Ап слоев). Главными трудностями на этом пути являются срыв поперечно-одномодовой генерации и повышение чувствительности конфигурации и ширины ближнего поля к отклонениям толщин либо показателей преломления слоев от заданных [7]. Такие изменения слоев могут иметь место за счет технологического разброса параметров слоев гетероструктуры при ее изготовлении, а также за счет изменения внешних факторов: температуры слоев при работе лазера, пространственно-неоднородного выедания носителей в активной области при больших мощностях и т.д.
С нашей точки зрения перспективным способом расширения ближнего поля является использование структур с возможностью вытекания части излучения из волновода в эмиттерные слои за счет туннелирования (вытекания) оптической волны через обкла-дочный слой. Благодаря этому происходит увеличение дефицитов усиления ненулевых мод относительно порога генерации фундаментальной моды [8]. В таких структурах эффективные показатели преломления мод волновода меньше показателей преломления эмиттерных слоев и не более чем на 0.05 больше, чем у обкладочных слоев волново-
да. Малое отличие эффективных показателей преломления мод и обкладочньтх слоев приводит к вытеканию излучения из волновода.
Целью нашей работы был подбор описанной в предыдущем абзаце структуры, оптимизированной по параметру отношения d/Г, а также по устойчивости этого и других наиболее важных волноводньтх параметров к вариациям составов слоев, допустимых при производстве и при работе лазера. Это делалось путем численных расчетов критичных для нас параметров качества волноводов для структур, получающихся из исходной путем варьирования в широких пределах толщин и показателей преломления некоторых ее слоев.
Методика. Расчеты конфигурации ближнего и дальнего полей и интегральных параметров мод подбираемых волноводов осуществлялись с помощью методик [10. 11]. реализованных в программном продукте Wave. На рис. 1 представлены примеры рассчитанных такими методиками распределений ближнего и дальнего полей основной моды, а также значений волноводных параметров d, d/Г, Ag, ad для соответствующего пространственного профиля показателя преломления слоев гетероструктурьт и с пороговым значением модового усиления, взятого равным во всех расчетах 10 см-1.
Наша методика подбора слоев структуры состояла в варьировании толщин и показателей преломления волноводных и обкладочньтх слоев структуры, параметры которой приведены в табл. 1. В подбираемой структуре было предложено использовать 1 квантовую яму из-за критерия оптимизации по параметру d/Г. Длина волны излучения, соответствующая выбранному составу квантовой ямьт, составляла 980 нм. Расчеты начинались с толщин волноводных слоев в 1 мкм^ и в дальнейшем толщина волноводных слоев изменялась на 0.1 мкм. При этом для уменьшения числа вариантов рассчитываемых структур рассматривались случаи с суммарной толщиной волноводньтх слоев свьтттте 1.5 мкм, так как нас интересовало широкое ближнее вертикальное поле нулевой модьт, лежащее внутри волноводньтх слоев. Толщина и показатель преломления обкладочньтх слоев согласовывались с каждым из рассчитываемых вариантов слоев волновода для селективного вытекания ненулевых мод в подложку. В дальнейшем в структуру добавлялись еще волноводные слои толщинои от 0.1 мкм и отличием в показателе преломления от начальных волноводньтх слоев на 0.05 и 0.1. Потом толщина этих слоев также варьировалась с шагом 0.1 мкм.
Для изготавливаемых структур крайне важно знать, насколько чувствительно поперечное распределение поля и сохраняется ли поперечно-одномодовая генерация при изменениях показателей преломления волноводньтх и обкладочньтх слоев, возможных
при производстве структуры. Поэтому при подборе структуры мы учитывали возможность возникновения отклонений состава и толщин слоев. Особое внимание следует уделять отклонениям показателя преломления обкладочньтх и волноводньтх слоев, так как показатель преломления подложки нам известен заранее, до выращивания на ней гетероструктурьт. а окружающие квантовую яму слои растятся поатомно с высокой
ТОЧНОСТЬЮ.
Таблица 1
Названия и общая, формула составов и толщин слоев подбираемой гетероструктуры
Л"2 слоя Назначение слоя Состав слоя Толщина слоя
1 Металлический контакт Аи, Ъ\\
2 Полупроводниковый контактный (р-эмиттерный) слой СаАэ 100 150 нм
3 Обкладочньтй слой А1ЖС а1-жАэ ¿1
4 Волноводньтй слой Al0.11Ga0.s9As ¿2
5 Барьерный слой Alo.3Gao.7As 40 нм
6 Переходной слой СаАэ 40 нм
7 Активный квантово-размерньтй слой ¡пСаАэ 8 нм
8 Переходной слой СаАэ 40 нм
9 Барьерный слой А10.^а0.7Аэ 40 нм
10 Волноводньтй слой Al0.11Ga0.89As ¿3
11 Обкладочньтй слой А1ув&1-уАб ¿4
12 Подложка (и-эмиттерный слой) ваАэ 100 мкм
В качестве критерия стабильности структуры были взяты условия дефицита усиления мод ненулевого порядка относительно усиления нулевой моды на пороге генерации не менее 3 см-1 и изменения отношения ¿/Г основной моды не более 30% при вариациях показателей преломления люоых волноводньтх и обкладочньтх слоев в пределах ±0.01, а также при отклонениях иа до ±0.1, моделирующих отклонения показателя преломления активной области иа из-за неоднородности концентрации носителей по длине и ширине лазера [5]. Также важно, чтобы величина материального усиления в активной области, подбираемая при численных расчетах по условию преодоления фундаментальной модой порога генерации, не превышала 4000 см-1. Такое наибольшее значение усиления дает эксперимент.
Для проверки выполнения наших критериев необходимо было проверять большое количество вариантов изменения слоев. Однако были выявлены 4 ключевых граничных для структур варианта изменения показателей преломления слоев, которые показаны на рис. 2.
Варианты, в которых уменьшается показатель преломления обкладочньтх слоев и одновременно обеспечивается сохранность нулевой моды при максимальном смещении пика интенсивности от активной области (рис. 2(а), (б)), важны для проверки сохранения условий одномодовости генерации. Варианты, в которых увеличивается показатель преломления обкладочньтх слоев и также одновременно обеспечивается максимальное смещение пика интенсивности нулевой модьт от активной области (рис. 2(в). (г)), важны для проверки работоспособности структуры при максимальном вытекании в поглощающие слои. Также эти описанные 4 случая выявляли наибольшее отклонение параметра ¿/Г основной моды.
Таблица 2
Названия, составы, и толщины слоев подобранной гетероструктуры
№ Назначение Состав Толщина слоя. Показатель
слоя слоя слоя мкм преломления
1 Металлический контакт Ди 0.24
2 Полупроводниковый контактный (р-эмиттерный) слой СаАэ 0.16 3.53
3 Обкладочньтй слой Al0.32Ga0.68As 0.3 3.32
4 Волноводньтй слой Al0.11Ga0.s9As 1.1 3.45
5 Барьерный слой Al0.29Ga0.71As 0.05 3.34
6 Переходной слой СаАэ 0.04 3.53
7 Активный квантово~ размерный слой ¡пСаАэ 0.008 3.48
8 Переходной слой СаАэ 0.04 3.53
9 Барьерный слой А10.^а0.7Аэ 0.04 3.33
10 Волноводньтй слой Al0.11Ga0.89As 1.0 3.45
11 Обкладочньтй слой Al0.17Ga0.83As 1.2 3.41
12 Подложка (и-эмиттерный слой) СаАэ 3.528
Рис. 2: Вертикальные распределения ближнего поля фундаментальной моды и параметры качества для волноводов, соответствующих структурам, отличающимся от подобранной структуры уменьшением показателя преломления активной области па на 0.1 и: (а), (в) - увеличением показателя преломления всех слоев левого волновода на 0.01 при его уменьшении для всех слоев правого волновода на 0.01 и уменьшении (а) или
0.01 0.01
0.01
0.01
соответствующих структурах.
Результаты. На рис. 1 представлено распределение вертикального профиля ближнего и дальнего поля, наборы волноводньтх параметров и профиль показателя преломления слоев структуры, подобранной с помощью предложенного нами метода. В табл. 2 приведены параметры слоев этой гетероструктурьт. На рис. 2 представлены искажения вертикального профиля ближнего поля и наборы волноводных параметров для подо~ бранной нами структуры в случае отклонений параметров слоев, введенных по нашей методике проверки стабильности излучательньтх характеристик.
Обсуждение. Как показано на рис. 1(6) доля излучения нулевой моды, туннелирую-
тцего из волновода в сторону подложки, составила 3.3%. На диаграмме направленности
дальнего поля виден боковой пик, соответствующий этому вытеканию. Ад составил при -1
Из данных, полученных нами при расчетах искажений конфигураций ближнего поля и изменений волноводных параметров, следует, что подобранная структура устойчива к отклонениям показателей преломления слоев и удовлетворяет введенному нами критерию стабильности. Граничные варианты просчитывавшихся отклонений слоев представлены на рис. 2. Видно, что сохраняется оптимизация структуры по отношению ¿/Г,
туннелирование излучения из волновода нулевой моды в эмиттерньте слои не превыттта-
Ад -1
одномодовой генерации.
Попытки дальнейшего расширения ближнего вертикального поля за счет добавления слоев волновода ограничивались уменьшением дефицита ненулевых мод относительно нулевой моды или ростом вытекания излучения основной моды из волновода в эмиттерньте слои. Добавление слоев, смещающих максимум интенсивности первых и вторых мод из активной области, решает проблему недостаточных, дефицитов усиления мод 1 и 2 порядка, но при этом возникает вероятность генерации мод более высоких порядков при определенном сочетании отклонений параметров слоев волновода от заданных. Увеличение вытекания излучения в поглощающие слои для решения проблемы возбуждения ненулевых мод также неэффективно, т.к. дальнейшее ослабление волновода сопровождалось значительным увеличением потерь для фундаментальной модьт. При данном выборе критерия стабильности, по-видимому, существует ограничение расширения вертикального профиля ближнего поля на уровне 1.6 1.7 мкм, которое и достигнуто в данной работе. Оно обусловлено "силой" волновода, и не зависит от конфигурации (числа, толщин и показателей преломления слоев) волновода, эту "силу" обеспечивающего.
Выводы. Предложена методика подбора структуры с расширенным волноводом, обладающей устойчивостью волноводньтх параметров к вариациям составов слоев, по принципу максимизации селективности туннелирования излучения из волновода в эмиттерньте слои для ненулевых мод. По этой методике была подобрана структура, оптимизированная по параметру отношения ¿/Г и отвечающая указанным критериям. Подобранная структура имеет ширину области, занимаемой излучением, Н ~ 1.6 мкм, отвечает критерию устойчивости, а также отличается малым вытеканием излучения в подложку и имеет повышенную устойчивость распределения поля волноводной моды к отклонениям составов слоев.
Работа была выполнена в рамках темы Л"2 01201156501 "Исследование свойств опто-электронных материалов и структур и возможностей их использования в лазерной технике, информатике и медицине" при частичной финансовой поддержке гранта РФФИ М()Л_А_2012 № 12-02-31345 и УНК ФИАН.
ЛИТЕРАТУРА
[1] И. В. Акимова, А. П. Богатов, А. Е. Дракин, В. П. Коняев, Квантовая электроника 38(11), 647 (1998).
[2] Д. Р. Мифтахутдинов, И. В. Акимова, А. П. Б огатов и др.? Квантовая электроника 25(7), 993 (2008).
[3] Д. Р. Мифтахутдинов, А. П. Богатов, А. Е. Дракин, Квантовая электроника 40(7), 583 (2010).
[4] Д. Р. Мифтахутдинов, А. П. Богатов, А. Е. Дракин, Квантовая электроника 40(7), 589 (2010).
[5] С. А. Плисюк, Д. В. Батрак, А. П. Богатов, А. Е. Дракин, Квантовая электроника 36(11), 1058 (2006).
[6] А. П. Богатов, А. Е. Дракин, А. А. Стратоников, В. П. Коняев, Квантовая электроника 30(5), 401 (2000).
[7] А. П. Богатов, Т. И. Гутцик, А. Е. Дракин и / ш.. Квантовая электроника 38(10), 935 (2008).
[8] G. Erbert, F. Bugge, J. Fricke, et al., IEEE J. Sel. Topics Quantum Electron. 11(5), 1217 (2005).
[9] S. Ishizaka, K. Muro, Т. Fujimoto, Y. Yamada, Патент US5764668 (1998).
[10] A. A. Stratoimikov, А. P. Bogatov, А. Е. Drakin, F. F. Kamenets. Appl. Opt. 4. 535 (2002).
[11] В. В. Поповичев, Е. И. Давыдова. А. А. Мармалюк и / Квантовая электроника 32(12), 1099 (2002).
Поступила в редакцию 1 апреля 2013 г.