Известия РАН. Серия физическая, 2022, T. 86, № 7, стр. 1008-1012

Конкретизация срока службы мощного ИК лазерного диода с использованием модели автономного канала генерации

В. В. Близнюк 1*, В. А. Паршин 1, А. Г. Ржанов 2, А. Е. Тарасов 1

1 Федеральное государственное бюджетное образовательное учреждение высшего образования “Национальный исследовательский университет “МЭИ”
Москва, Россия

2 Федеральное государственное бюджетное образовательное учреждение высшего образования “Московский государственный университет имени М.В. Ломоносова”, физический факультет
Москва, Россия

* E-mail: 40595141@mail.ru

Поступила в редакцию 14.02.2022
После доработки 28.02.2022
Принята к публикации 23.03.2022

Полный текст (PDF)

Аннотация

Предлагается методика оценочного расчета срока службы серийных мощных лазерных диодов с широким контактом по 50-часовому натурному исследованию выборки из большой партии этих приборов и спектру конкретного прибора в начале его эксплуатации. Для обоснования методики анализируется распад излучения в таких лазерах на отдельные нефазированные каналы генерации. При этом лазерные диоды представляются как набор одномодовых излучателей.

ВВЕДЕНИЕ

Мощные полупроводниковые лазерные диоды (ЛД) с квантово-размерными гетероструктурами широко используются в разных областях науки, техники и медицины. К таким лазерам принято относить ЛД с мощностью непрерывной генерации от сотен милливатт до десятков ватт [15]. При серийном производстве ЛД на первый план выходит их технико-экономический параметр – срок службы [6]. Известные методики прогнозирования срока службы ЛД основаны на использовании дорогостоящей контрольно-измерительной аппаратуры [1, 5]. В настоящее время для определения срока службы ЛД используются испытания по их ускоренному старению при повышенных значениях температуры окружающей среды и тока накачки. Это во многом обусловлено тем, что первые признаки деградации проявляются не ранее, чем через две-три тысячи часов наработки.

Современные методики прогнозирования срока службы мощных ЛД основаны на анализе временнóй зависимости медленно изменяющихся параметров лазерного излучения: мощности, диаграммы направленности и степени линейной поляризации (контраста). Данные методики позволяют определить срок службы некоторого числа экземпляров из партии ЛД, а потом распространить полученные результаты на всю партию [6]. Однако их реализация связана со значительным расходованием ресурса лазеров и использованием трудоемких статистических методов обработки больших массивов чисел. Основным же недостатком подобных методик является невозможность перейти от обобщенного прогноза срока службы партии лазеров к определению срока службы конкретного экземпляра ЛД из серии. На важность определения качества отдельно взятого прибора электронной техники обращено внимание в [711].

ПРЕДЛАГАЕМЫЙ МЕТОД ПРОГНОЗА СРОКА СЛУЖБЫ ЛАЗЕРНЫХ ДИОДОВ С ШИРОКИМ КОНТАКТОМ

Одним из возможных решений данной проблемы является использование методики [10, 11]. Она базируется на численном анализе формы линии ${{f}_{{{\text{эксп}}}}}\left( {\frac{{\nu - {{\nu }_{0}}}}{{{{\Delta }}\nu }}} \right),$ огибающей спектр излучения одномодовых лазерных диодов (ОЛД), и результатов измерений токовых зависимостей спектральных характеристик, полученных в ходе ускоренных испытаний контрольной партии лазерных диодов через 200 ч наработки в реальном масштабе времени, то есть на начальном этапе их эксплуатации. Применение метода иллюстрирует рис. 1. Анализ проводится путем сравнения $~{{f}_{{{\text{эксп}}}}}\left( {\frac{{\nu - {{\nu }_{0}}}}{{{{\Delta }}\nu }}} \right)$ с нормированной гауссовой функцией. Вводится интегральный параметр Aнач, характеризующий степень отклонения функции ${{f}_{{{\text{эксп}}}}}\left( {\frac{{\nu - {{\nu }_{0}}}}{{{{\Delta }}\nu }}} \right)$ от функции Гаусса ${{f}_{G}}\left( {\frac{{\nu - {{\nu }_{0}}}}{{{{\Delta }}\nu }}} \right)$ в пределах ширины линии, огибающей спектр излучения, в начальной стадии эксплуатации лазера при разных токах накачки, как это изображено на рис. 1а и 1б). С математической точки зрения интегральный параметр A представляет собой норму близости нормированной экспериментальной функции огибающей спектра ${{f}_{{{\text{эксп}}}}}\left( {\frac{{\nu - {{\nu }_{0}}}}{{{{\Delta }}\nu }}} \right)$ от идеализированной спектральной функции, соответствующей генерации на фундаментальной моде ${{f}_{G}}\left( {\frac{{\nu - {{\nu }_{0}}}}{{{{\Delta }}\nu }}} \right){\text{:}}$

(1)
$A = 1 - \int\limits_{{{\nu }_{{\frac{1}{2}min}}}}^{{{\nu }_{{\frac{1}{2}max}}}} {\frac{{\left| {\frac{{{{f}_{{{\text{эксп}}}}}\left( {\frac{{\nu - {{\nu }_{0}}}}{{\Delta \nu }}} \right)}}{{{{f}_{G}}\left( {\frac{{\nu - {{\nu }_{0}}}}{{\Delta \nu }}} \right)}} - 1} \right|d\nu }}{{\Delta \nu }}} ,$
где ∆ν – полуширина линии спектра фундаментальной моды, ${{\nu }_{{\frac{1}{2}max}}}$ и ${{\nu }_{{\frac{1}{2}min}}}$ – частоты, соответствующие полуширине спектральной линии моды.

Рис. 1.

Кривая ${{f}_{i}}\left( {\frac{{{{\nu }} - {{{{\nu }}}_{{0i}}}}}{{{{\Delta }}{{{{\nu }}}_{i}}}}} \right)$ (сплошная линия), характеризующая спектральное распределение интенсивности длинноволновой составляющей спектра излучения мощного лазерного диода АТС-С200-100-980 с наработкой 200 ч при токе накачки: 940 (а); 980 мА (б). Пунктирная линия – нормированная функция Гаусса.

Разработанная методика прогнозирования срока службы ОЛД по его спектральной характеристике может быть использована только в том случае, если имеет место режим генерации на фундаментальной моде.

Нами предлагается распространить методику прогнозирования полезного времени наработки ОЛД [711] для анализа мощных ЛД с широким контактом – ЛДШК [8]. Измерение формы линии, огибающей спектр излучения ЛДШК в течение первых часов работы лазера, позволило нам зафиксировать два-три ярко выраженных спектральных пика (линии) генерации. Каждый пик соответствовал автономному пространственному каналу генерации в активной области ЛДШК. В [11] была разработана методика, которая позволяет прогнозировать срок службы серийного изготовленного ЛДШК за время, не превышающее пятидесяти часов, в число которых входит время ускоренных испытаний на старение, соответствующее двумстам часам работы в реальном масштабе времени. Через двести часов испытаний в реальном масштабе времени, количество каналов генерации и линий в спектре излучения лазера ${{f}_{{emp}}}\left( {{\nu }} \right),$ возрастает от 2–3 до 4–5 (рис. 2) [11].

Рис. 2.

Спектр излучения мощного лазерного диода АТС-С200-100-980 с наработкой 200 ч. ${{f}_{{{\text{эксп}}}}}\left( {{\nu }} \right),$ и его составляющие ${{f}_{i}}\left( {\frac{{{{\nu }} - {{{{\nu }}}_{{0i}}}}}{{{{\Delta }}{{{{\nu }}}_{i}}}}} \right)$ при двух значениях тока накачки: 940 (а); 980 мА (б).

Мощные ЛДШК, в отличие от ОЛД, характеризуются сложной формой спектральной характеристики. что исключает простую трактовку токовой зависимости функции ${{f}_{{emp}}}\left( {{\nu }} \right)$ при разработке алгоритма прогнозирования срока службы конкретного экземпляра лазера. Именно поэтому в [10, 11] было предложено проводить анализ спектра ЛДШК путем разложения его на N составляющих с профилями линии ${{f}_{i}}\left( {\frac{{\nu - {{\nu }_{{0i}}}}}{{{{\Delta }}{{\nu }_{i}}}}} \right),$ характерными для излучения одномодового лазера с центральной частотой ${{{{\nu }}}_{{0i}}}.$

Нами экспериментально установлено, что с ростом тока накачки спектральная характеристика ЛД ${{f}_{{emp}}}\left( {{\nu }} \right)$ претерпевает значительные изменения, которые проявляются в изменении соотношения амплитуд спектральных составляющих и ширин линии излучения ${{f}_{i}}\left( {\frac{{{{\nu }} - {{{{\nu }}}_{{0i}}}}}{{{{\Delta }}{{{{\nu }}}_{i}}}}} \right),$ а также числа таких составляющих. Так, на рис. 2 видно, что число составляющих спектра ЛДШК ${{f}_{{emp}}}\left( {{\nu }} \right)$ увеличивается от четырех до пяти. Численный анализ каждой из четырех кривых ${{f}_{i}}\left( {\frac{{{{\nu }} - {{{{\nu }}}_{{0i}}}}}{{{{\Delta }}{{{{\nu }}}_{i}}}}} \right)$ на рис. 2а, показывает, что соответствующие им значения Ai лежат в интервале от 0.93 до 0.95. Из анализа пяти контуров линий ${{f}_{i}}\left( {\frac{{{{\nu }} - {{{{\nu }}}_{{0i}}}}}{{{{\Delta }}{{{{\nu }}}_{i}}}}} \right)$ на рис. 2б следует, что им соответствуют значения параметра Ai, изменяющиеся от 0.82 до 0.86. Таким образом, увеличение тока накачки приводит в целом к уменьшению параметров Ai, а значит, согласно [11], и к уменьшению срока службы ЛДШК. Однако из-за большого числа параметров Ai возникает вопрос, по которому из них следует прогнозировать конкретный срок службы лазера. Принимая во внимание, что меньшему значению параметра Ai соответствует меньший срок службы τ, для оценки срока службы $\tau $ конкретного ЛДШК при фиксированном токе накачки предлагается использовать наименьшее из всех значений параметра Ai при каждом значении тока накачки и подставлять его в эмпирическую формулу, подобную той, которая использовалась в [10, 11]:

(2)
$\frac{\tau }{{{{\tau }_{{max}}}}} = {{\left( {{{A}_{{min}}}} \right)}^{m}},$
где значение параметра m находится по той же методике, что в [10, 11]. В проведенных ранее экспериментах с ОЛД средней мощности вычисленное по статистическим данным значение показателя m оказалось равным 1/3. Однако, значение этого параметра, т.е. вид конкретной функциональной зависимости (2), в общем случае не носит универсального характера, и может отличаться для разных промышленных партий и структур ЛД. Данный аспект требует дополнительных исследований.

Для прогнозирования срока службы ЛДШК при рассмотрении партии из 100 приборов при браковочном уровне качества 10% был использован ГОСТ 16493–70 “Межгосударственный стандарт. Качество продукции. Статистический приемочный контроль по альтернативному признаку. Случай недопустимости дефектных изделий в выборке”, согласно которому, достаточно обеспечивать репрезентативную выборку из 25 экземпляров лазеров.

По результатам проведенных авторами исследований 25 ЛДШК АТС-С200-100-980 при токе накачки 940 мА из контрольной партии установлены зависимости τ/τmax (Анач), которые представлены на рис. 3. Стоит отметить, что из контрольной партии было отбраковано 3 экземпляра ЛДШК, у которых мощность излучения после 200 ч наработки в реальном масштабе времени уменьшилась более чем на 10% от начальной мощности. Время наработки ЛДШК варьировалось от 2230 до 2400 ч. Было определено, что погрешность определения интегрального параметра Анач составила 4.8%.

Рис. 3.

Зависимости τ/τmax(Aнач) ЛДШК из контрольной партии (${{{{\tau }}}_{{max}}}$ = 2400 ч) при значении тока накачки 940 мА: сплошная линия – расчетная функция $T = \sqrt[3]{{{{A}_{{{\text{нач}}}}}}};$ точки – результаты эксперимента.

Анализ распада излучения в ЛДШК на отдельные нефазированные каналы генерации был проведен в работах [1215]. В этих работах были подробно описаны причины такого распада и предложены способы аналитических и численных способов оценки количества каналов и ширины спектральных составляющих этих приборов. Причиной разделения излучения на независимые каналы генерации оказывается воздействие нескольких факторов. Можно объяснить это явление следующим образом: для того, чтобы поле излучения в определенной части волновода, например, на краю активной области, могло быть связано по фазе с полем на другом краю этой области, волновому лучу необходимо несколько раз пробежать между зеркалами. Если длина когерентности меньше оптического пути такого пробега, то излучение разделяется на нефазированные каналы генерации. В ином случае распадения излучения на каналы не будет. Это необходимое условие такого распада. Таким образом, это явление обусловлено, в первую очередь, малой длиной когерентности ЛД, которая составляет 2–6 см [1617], а также превышением порога генерации по всей площади контакта (W = 100 мкм), а также в пределах всего контура усиления квантовой ямы (около 8 мэВ). Кроме того, на формирование канала влияет нелинейный эффект пространственной самофокусировки за счет зависимости рефракции от концентрации неравновесных носителей, и трансформации латерального профиля инверсии, что также было проанализировано в работах [1214].

В работе [18] была получена связь между длиной продольной когерентности излучения ${{L}_{{{\text{ког}}}}}$ и минимальной шириной когерентного канала ${{w}_{0}}$ в латеральной плоскости ЛДШК. Наиболее близкий к экспериментальным данным [4, 7, 8, 11] результат получается при использовании модели резонатора ЛД как планарного волновода, ограниченного плоскими зеркалами [18]:

(3)
${{w}_{0}} = \sqrt {\frac{{\lambda {{L}_{{{\text{ког}}}}}}}{{2\pi {{n}_{{{\text{эф}}0}}}}}} ,$
где $\lambda $ и ${{L}_{{{\text{ког}}}}}$ – соответственно, длина волны и длина продольной когерентности излучения, ${{n}_{{{\text{эф0}}}}}$ – эффективный показатель преломления фундаментальной латеральной моды. Учет самофокусировки излучения посредством нелинейной рефракции в этом случае позволяет рассмотреть в качестве модельного приближения локальный волновод с квадратичным профилем, и также оценить ширину установившегося канала генерации через длину когерентности.

Превышение порога генерации ЛД хотя бы на 10–15% означает, что усредненное по времени поле излучения должно присутствовать на всей ширине активной области W. Изменения интенсивности излучения по латеральной координате могут быть заметны только в пределах длины диффузии неравновесных носителей (примерно 5–15 мкм). Это обосновывает появление каналов, разделенных в пространстве при наличии ограниченной когерентности излучения. Также подразумевается сплошное заполнение спектральной линии усиления ЛД по всей ее ширине, где есть превышение порога генерации (см. рис. 2). Это объясняет разнесение каналов по спектру (отсутствие вырождения по частоте), что наблюдалось в эксперименте [7, 15], и определяет связь между заполнением каналами генерации как в пространственной, так и в спектральной областях. В пределах ограничений, накладываемых пространственной диффузией и внутризонной термолизацией носителей, излучение должно заполнять соответствующие области относительно равномерно.

ЗАКЛЮЧЕНИЕ

Таким образом показано, что для конкретизации срока службы ЛДШК его следует рассматривать как набор ОЛД с изменяющимся в процессе наработки числом излучающих элементов. При этом для каждого элемента можно использовать методику быстрого прогнозирования срока службы ОЛД по токовой зависимости интегрального параметра Aнач в начальный период эксплуатации лазера. В этом случае срок службы ЛДШК определяется по наименьшему из всей совокупности Ai значению этого параметра. Максимальный срок службы, равный 2400 часам, был определен в ходе ускоренных испытаний ЛДШК с наибольшим значением Aнач = 0.95 из контрольной партии.

Список литературы

  1. Богатов А.П., Дракин А.Е., Стратонников А.А., Коняев В.П. // Квант. электрон. 2000. Т. 30. № 5. С. 401.

  2. Слипченко С.О., Бондарев А.Д., Винокуров Д.А. и др. // ФТП. 2009. Т. 43. № 1. С. 119.

  3. Слипченко С.О., Винокуров Д.А., Пихтин Н.А. и др. // ФТП. 2004. Т. 38. № 12. С. 1477.

  4. Лютецкий А.В., Пихтин Н.А., Фетисова Н.В. и др. // ФТП. 2009. Т. 43. № 12. С. 1646.

  5. Жуков А.Е. Физика и технология полупроводниковых наноструктур. СПб.: Элмор. 2007. 304 с.

  6. Воробьев Л., Софронов А., Фирсов Д. и др. // Фотоника. 2012. Т. 31. № 1. С. 20.

  7. Близнюк В.В., Брит М.А., Гадаев И.С. и др. // Изв. РАН. Сер. физ. 2015. Т. 79. № 12. С. 1666; Bliznyuk V.V., Brit M.A., Gadaev I.S. et al. // Bull. Russ. Acad. Sci. Phys. 2015. V. 79. No. 12. P. 1458.

  8. Близнюк В.В., Коваль О.И., Паршин В.А. и др. // Учен. зап. физ. фак-та МГУ. 2018. № 6. С. 1860501.

  9. Близнюк В.В., Коваль О.И., Паршин В.А. и др. // Изв. РАН. Сер. физ. 2018. Т. 82. № 11. С. 1504; Bliznyuk V.V., Koval O.I., Parshin V.A. et al. // Bull. Russ. Acad. Sci. Phys. 2018. V. 82. No. 11. P. 1367.

  10. Близнюк В.В., Коваль О.И., Паршин В.А. и др. // Изв. РАН. Сер. физ. 2020. Т. 84. № 2. С. 225; Bliznyuk V.V., Koval O.I., Parshin V.A. et al. // Bull. Russ. Acad. Sci. Phys. 2020. V. 84. No. 2. P. 173.

  11. Близнюк В.В., Паршин В.А., Ржанов А.Г., Тарасов А.Е. // Изв. РАН. Сер. физ. 2021. Т. 85. № 2. С. 255; Bliznyuk V.V., Parshin V.A., Rzhanov A.G., Tarasov A.E. et al. // Bull. Russ. Acad. Sci. Phys. 2021. V. 85. No. 2. P. 184.

  12. Ржанов А.Г. // Изв. РАН. Сер. физ. 2021. Т. 85. № 2. С. 250; Rzhanov A.G. // Bull. Russ. Acad. Sci. Phys. 2021. V. 85. No. 2. P. 180.

  13. Ржанов А.Г. // Изв. РАН. Сер. физ. 2020. Т. 84. № 2. С. 220; Rzhanov A.G. // Bull. Russ. Acad. Sci. Phys. 2020. V. 84. No. 2. P. 169.

  14. Rzhanov A.G. // EPJ Web Conf. 2019. V. 220. Art. No. 02013.

  15. Koval O.I., Rzhanov A.G., Solovyev G.A. // Phys. Wave Phenom. 2013. V. 21. No. 4. P. 287.

  16. http://www.holography.ru/files/holmich.htm#top.

  17. Адамов А.А., Баранов М.С., Храмов В.Н. // Научн.-техн. вестн. ИТМО. 2018. Т. 18. № 3. С. 356.

  18. Ржанов А.Г. // Изв. РАН. Сер. физ. 2018. Т. 82. № 11. С. 1508; Rzhanov A.G. // Bull. Russ. Acad. Sci. Phys. 2018. V. 82. No. 11. P. 1371.

Дополнительные материалы отсутствуют.