CC BY
57
Прикладные задачи
нелинейной теории колебании и волн
УДК 535.015; 535.14; 535.530; 537.86 https://doi.org/10.18500/0869-6632-2019-27-3-61-72
Влияние запаздывания на динамику оптического параметрического осциллятора с внутрирезонаторной накачкой полупроводниковым дисковым лазером
Л. А. Конкуров1, М.И. Балакин1, П. В. Купцов1,2, Ю.А. Морозов2
Саратовский государственный технический университет им. Гагарина Ю.А. Россия, 410054 Саратов, ул. Политехническая, 77 2 Саратовский филиал Института радиотехники и электроники им. В.А. Котельникова РАН Россия, 410019 Саратов, ул. Зеленая, 38 E-mail: lkochkurov@gmail.com, balakinmaxim@gmail.com, p.kuptsov@rambler.ru, yuri.mor@rambler.ru Автор для переписки Кочкуров Леонид Алексеевич, lkochkurov@gmail.com Поступила в редакцию 28.03.2019; принята к публикации 15.05.2019
Тема. Исследована динамика численной модели нелинейного оптического взаимодействия в резонаторе полупроводникового дискового лазера с учетом временного запаздывания. Проанализированы условия самовозбуждения колебаний, режимы стационарной генерации и их устойчивость. Методы. Для анализа устойчивости стационарной генерации был использован пакет ВБЕ^йоок Анализ динамических режимов более высокой размерности проводился с помощью численного интегрирования, построения фазовых портретов, спектров и расчета показателей Ляпунова. Результаты. Проведено численное моделирование динамики излучения в области неустойчивости стационарного состояния. Показано, что неустойчивость носит квазигармонический характер только вблизи бифуркации Андронова-Хопфа, и быстро принимает квазипериодический характер при вариации управляющих параметров. Изучены переходные процессы в системе. Обсуждение. Полученные результаты могут быть использованы для оптимизации параметров лазерных генераторов в устройствах спектроскопии высокого разрешения.
Ключевые слова: оптический параметрический осциллятор, полупроводниковый дисковый лазер, нелинейно-оптическое взаимодействие, динамическая система с запаздыванием.
Образец цитирования: Кочкуров Л.А., Балакин М.И., Купцов П.В., Морозов Ю.А. Влияние запаздывания на динамику оптического параметрического осциллятора с внутрирезонаторной накачкой полупроводниковым дисковым лазером// Известия вузов. ПНД. 2019. Т. 27, № 3. С. 61-72. https://doi.org/10.18500/0869-6632-2019-27-3-61-72
Финансовая поддержка. Работа выполнена при поддержке грантов РФФИ № 18-32-01028-мол_а и 18-08-00599-а, а также в рамках государственного задания ИРЭ им. В.А. Котельникова РАН.
https://doi.org/10.18500/0869-6632-2019-27-3-61-72
Impact of time delay on the dynamics of optical parametric oscillator with intra-cavity pumping by semiconductor disk laser
L. A. Kochkurov1, M. I. Balakin1, P. V. Kuptsov1'2, Yu. A. Morozov2
1Yuri Gagarin State Technical University of Saratov 77, Politechnicheskaya st., 410054 Saratov, Russia 2 Kotelnikov Institute of RadioEngineering and Electronics of RAS, Saratov Branch 38, Zelenaya st., 410019 Saratov, Russia E-mail: lkochkurov@gmail.com, balakinmaxim@gmail.com, p.kuptsov@rambler.ru, yuri.mor@rambler.ru Correspondence should be addressed to Kochkurov Leonid A., lkochkurov@gmail.com Received 28.03.2019; accepted for publication 15.05.2019
Theme. The dynamics of a nonlinear numerical model of a nonlinear optical interaction in the semiconductor disk laser resonator under influence of the time delay is investigated. The conditions of self-excitation, stationary generation modes and their stability are studied. Methods. The analysis of stationary generation stability was performed with DDE-Biftool package. Analysis of higher dimensional regimes was performed using numerical integration, construction of phase portraits, spectra and calculation of Lyapunov exponents. Results. A numerical simulation of the dynamics in the region of steady state instability shown, that the instability is quasi-harmonic only in the vicinity of Andronov-Hopf bifurcation, and quickly turns into quasi-periodic instability with variation of control parameters. Transient dynamics is studied. Discussion. The results can be used for optimization of laser generator parameters in high resolution spectroscopy devices.
Key words: optical parametric oscillator, semiconductor disk laser, nonlinear optical interaction, time delayed systems.
Reference: Kochkurov L.A., Balakin M.I., Kuptsov P.V., Morozov Yu.A. Impact of time delay on the dynamics of optical parametric oscillator with intra-cavity pumping by semiconductor disk laser. Izvestiya VUZ, Applied Nonlinear Dynamics, 2019, vol. 27, no. 3, pp. 61-72. https://doi.org/10.18500/0869-6632-2019-27-3-61-72
Acknowledgements. This work was supported by the Russian Foundation for Basic Research (projects no. 18-32-01028-мол_а and no. 18-08-00599_a), and was carried out at IREE RAS in the framework of the state task.
Введение
Когерентные источники непрерывного излучения в среднем и дальнем инфракрасном диапазонах находят широкое применение во многих прикладных областях, среди которых одной из самых значимых является спектроскопия высокого разрешения. Вследствие этого, исследование рабочих характеристик генераторов разностной частоты (ГРЧ) и оптических параметрических осцилляторов (ОПО) становится важной задачей [1]. Некоторые реализации таких источников могут быть достаточно просты для изготовления, компактны и удобны в применении [2,3]. Компактность и эффективность этих устройств может быть достигнута за счет использования схемы с внутрирезонаторным нелинейным преобразованием частоты. Пример подобного рода был, в частности, представлен в работе [4]. Однако предложенный в этой работе подход не позволяет реализовать условия, при которых длины волн сигнального и накачивающего излучений были бы достаточно близки друг другу для обеспечения холостого излучения с длиной волны порядка десятков микрометров. В работе [2] было предложено использовать общий резонатор как для волны накачки, так и для сигнальной волны, вследствие чего эта особенность была устранена. Кроме того, как было теоретически показано в работе [3], использование данного подхода позволяет создать как внутрирезонаторный ГРЧ, так и внутрирезонаторный ОПО.
Одним из факторов, определяющих возможность применимости таких устройств для прикладных задач, является устойчивость их стационарного режима генерации. Для анализа режимов генерации и их устойчивости используются математические модели, формализующие поведение динамических систем. В большинстве математических моделей внутрирезонаторных
генераторов с нелинейным частотным преобразованием применяется метод разложения поля по нормальным модам резонатора [3,5-9]. В приближении одной моды (или для несфазированных мод) такой метод не позволяет изучать возможные быстрые (на масштабе времени обхода по резонатору) колебания амплитуд оптических полей.
Для анализа быстрых амплитудных осцилляций (которые могут явиться причиной нестабильности стационарного режима) мы используем модель лазера со связанными резонаторами, один из которых является внешним по отношению к другому [10-13]. Внутренний резонатор образован полупроводниковым брэгговским зеркалом (БЗ) и излучающей поверхностью полупроводникового чипа. Внешний резонатор ограничивается выходным сферическим зеркалом с высоким коэффициентом отражения, с одной стороны, и упомянутым БЗ с другой. Было показано, что такой лазер может рассматриваться как динамическая система с запаздыванием, демонстрирующая сложное поведение [13]. Подобная модель кардинально отличается от широко известной модели Ланга-Кобаяши [14], в которой учитываются лишь слабые отражения от внешнего зеркала.
Недавно в работе [15] нами было выполнено развитие модели [13], позволяющее учесть влияние нелинейно-оптического взаимодействия в резонаторе полупроводникового дискового лазера на характеристики его излучения. В данной работе мы проводим углубленный анализ динамики излучения внутрирезонаторного ОПО.
1. Динамические уравнения оптического параметрического осциллятора
Устройство однорезонаторного ОПО, построенного на основе нелинейно-оптического взаимодействия в резонаторе полупроводникового дискового лазера (ПДЛ), схематически показано на рис. 1 [2,3]. Колебательная система выполнена по трехзеркальной схеме. Внешний резонатор, ограниченный вогнутым зеркалом с высоким коэффициентом отражения, связан с внутренним резонатором, образованным распределенным двухполосным БЗ и выходным торцом полупроводникового чипа. Внутренний резонатор по сути является активным зеркалом. Волна накачки (1.98 мкм) и сигнала (2.25 мкм) являются модами общего внешнего резонатора. Предполагается, что трехволновое нелинейно-оптическое взаимодействие осуществляется в нелинейном кристалле с регулярной доменной структурой на основе ваЛз. Выходное сферическое зеркало является прозрачным для холостого излучения с длиной волны 16.5 мкм. Для обеспечения высокого коэффициента отражения активного зеркала на волнах накачки и сигнала было использовано двух-
QW gain region
Nonlinear crystal
Double-band mirror
Primary pump (980 nm)
Pump (1.98 pm)
Signal (2.25 pm)
Idler radiation (16.5 pm)
Concave mirror
Рис. 1. Схема однорезонаторного оптического параметрического осциллятора Fig. 1. Schematic setup of a single cavity optical parametric oscillator (OPO)
полосное БЗ [16]. Первичная оптическая накачка активной области осуществляется с помощью диодного лазера с длиной волны 980 нм.
Устройство активного зеркала существенным образом не отличается от устройства этого элемента обычного ПДЛ: это зеркало содержит квантовые ямы на основе Gao.74Ino.26Sb, расположенные в узлах стоячей волны накачки для обеспечения максимально возможного усиления [17].
Построение модели внутрирезонаторного ОПО [15] основано на рассмотрении отражений в колебательной системе, состоящей из связанных резонаторов (см. рис. 3 из [15] и последующие пояснения). В соответствии с подходом [12,13], уравнения математической модели формулируются для переменных, относящихся к внутреннему резонатору (область длиной Lin на этом рисунке).
Следуя [13,15], можно получить скоростные (балансные) уравнения, представляющие динамическую систему с запаздыванием, обусловленным задержкой сигнала при обходе по внешнему резонатору:
a 1 = п
a 2 V1
п
(Gi - 1) + 1 - l) - I(a2 + а2т)
-1 + 1 (а^ - V + T (ai + а1т) а2 '
ai,
= о — v1 — G1a1.
(1)
Здесь число фотонов ai и число носителей vi нормализованы на a10 и v1th, соответственно; aïo = vith/n, где n = тг/xph, тг и Tph = as)-1 время жизни носителей и фотонов во внутреннем резонаторе (vg - групповая скорость для излучения накачки и сигнала), v1th - пороговое значение числа носителей, T = 2asLin - потери в резонаторе за проход, причем as - ain (Lin) 1 ln[rDBM(1 - r2)R], ain - внутренние потери вследствие рассеяния и нерезонансного поглощения во внутреннем резонаторе, rDBM, r и R - коэффициенты отражения двухполосного БЗ, внешней поверхности активного зеркала и сферического выходного зеркала, соответственно. Нормированный коэффициент усиления Gi принимается в форме
G1 = 1 + G10 ln V1,
(2)
где G10 = 4mGQw/T, m — число квантовых ям в активной области, Gqw коэффициент усиления, приходящийся на одну квантовую яму. Принимая во внимание эти параметры, пороговое значение числа носителей можно записать в виде v1th = mn^^Nt exp(1/G10), где wpp радиус пучка первичной оптической накачки, Nt плотность числа носителей, соответствующая просветлению активной среды (transparent carrier density). Коэффициент нелинейного взаимодействия определяется в форме 8 = /%т, где Xin — 2Lin/vg — время обхода внутреннего резонатора, ^ = 2у/[л(^2 + w|)j. При выводе выражений предполагалось, что оптические поля накачки и сигнала представлены в форме гауссовых пучков с радиусами перетяжек внутри кристалла wi и W2,
Y =
32Zo
П1П2П3
d2 d14
jLL
X1X2
(3)
Zo = 120п - волновое сопротивление вакуума, ^14 - элемент тензора нелинейной восприимчивости кристалла GaAs, Ьс - длина кристалла, Х1 и Х2 - длина волны накачки и сигнальной волны, П1,2,3 - коэффициенты преломления, Ншз - энергия фотонов холостой волны. В уравнениях (1) о соответствует мощности внешней первичной накачки Р, нормированной на пороговое значение Р^. Точки над переменными означают дифференцирование по величине Ь/хг.
2. Численные результаты
Численный анализ динамики системы выполнялся с применением стандартных методов нелинейной динамики: построение временных реализаций, фазовых портретов, фурье-спектров. Бифуркационный анализ состояния равновесия системы проводился с помощью пакета ЭЭБ-В1Иоо1 [18]. Также была выполнена проверка наличия хаоса, для чего вычислялся старший показатель Ляпунова. Необходимо отметить, что для изучения динамики системы были выбраны значения запаздывания т, нормированной мощности первичной накачки о/ооро и длины кристалла Ьс, которые отвечают возможности практической реализации. Результаты, полученные разными методами, хорошо согласуются между собой.
Ранее полученные результаты [15] показывают, что в случае внутрирезонаторной генерации разностной частоты / для всех практически важных значений управляющих параметров состояние равновесия системы остается устойчивым. Другими словами, наблюдается генерация с постоянной амплитудой. Этот результат представляется важным с точки зрения использования таких генераторов в спектроскопии высокого разрешения.
В случае параметрической генерации ситуация сложнее. На плоскости управляющих параметров имеются области, где устойчивость состояния равновесия нарушается. На рис. 2 в виде Э, Э и С и Э и С и В обозначены такие области. Каждая из областей отвечает определенному значению длины нелинейного кристалла Ьс при прочих фиксированных значениях параметров. При переходе через линии 1АВ,1ВС,1СВ наблюдается потеря устойчивости для Ьс = 7, 6 и 5 мм, соответственно. При построении графиков о нормировалась на значение Ооро = 11.03, соответствующее значению величины первичной накачки на пороге параметрической генерации при Ьс = 5 мм. При переходе через границы данных областей старшая пара комплексно-сопряженных корней характеристического уравнения приобретает положительную действительную часть, наблюдается супекритическая бифуркация Андронова-Хопфа. Однако область устойчивости возникающих периодических колебаний ограничена, при вариации управляющих параметров они быстро сменяются квазипериодическими колебаниями.
Рис. 2. Разбиение плоскости параметров а/аоро — т на области с различными типами устойчивости. Линии суперкритической бифуркации Андронова-Хопфа обозначены lAB,lBC,lCD. Область А отвечает устойчивой генерации для всех рассматриваемых значений длины кристалла. Области D, D U C и D U C U B соответствуют неустойчивости стационарного состояния системы при Lc = 5; Lc = 6; Lc = 7 мм
Fig. 2. Areas on the parameter plain а/аоро — т with different types of stability. The lines of the supercritical Andronov-Hopf bifurcation are marked lAB, lBC, lCD. Region A corresponds to stable generation for all considered values of the crystal length. The areas D, D U C and D U C U B correspond to the instability of the steady state of the system for Lc = 5 mm, Lc = 6 mm, Lc = 7 mm
Для иллюстрации на рис. 3 и 4 приведены характерные колебательные режимы излучения для различных значений нормированной мощности первичной накачки о/ооро при фиксированном времени запаздывания т = 0.1 и длине кристалла Ьс = 5 мм. Для малых значений о, то есть за пределами области Б на рис. 2, после завершения переходного процесса устанавливается режим стационарной генерации (см. рис. 3). Переход в заштрихованную область Б сопровождается возбуждением периодических колебаний амплитуды на частоте/ = 14п/т, что соответствует седьмой гармонике межмодовых биений для рассматриваемого резонатора или 35 ГГц. Необходимо отметить, что область устойчивости периодических колебаний находится в непосредственной близости от бифуркационной линии стационарного состояния. Так, для т = 0.1 бифуркация Андронова-Хопфа происходит при о/ооро ~ 2.250, а потеря устойчивости периодических колебаний - при о/оОРО « 2.257. Бифуркационный анализ показывает, что характер развивающейся неустойчивости меняется на квазипериодический в результате бифуркации Неймарка-Сакера, в фазовом пространстве системы возникает устойчивый двухчастотный тор (см. рис. 4).
Для подтверждения вывода о квазипериодическом характере динамики системы в области неустойчивости состояния равновесия был вычислен старший показатель Ляпунова. Применявшаяся для этого стандартная процедура состоит в том, что сначала для изучаемой системы выводится вариационное уравнение, представляющее собой линеаризацию исходной системы вдоль траектории и описывающее эволюцию вектора инфинитезимального возмущения этой траектории. Так как вариационное уравнение линейно по отношению к вектору возмущения, норма этого вектора растёт или убывает в среднем по закону экспоненты, показатель которой называется показателем Ляпунова. Для вычисления этого показателя нужно решать вариационное уравнение совместно с исходной системой, периодически выполняя перенормировку вектора возмущения. При этом логарифмы накопленных к моменту очередной перенормировки норм следует суммировать. Тогда, пройдя достаточно протяжённый отрезок траектории и усреднив накопленные логарифмы норм по времени, мы получим старший показатель Ляпунова. Отметим, что при вычислении более одного показателя Ляпунова применяют алгоритм, идея которого впервые была предложена независимо Бенеттином с соавторами [19] и Шимадой и Нагашимой [20]. Подробное описание этого алгоритма, включающего наряду с перенормировками также ортогонализации векторов возмущений, можно найти, например, в статьях [21,22].
Рис. 3. Устойчивый режим генерации при а/аоро = 1.5. Показаны (a) временная реализация сигнального излучения и (b) фазовый портрет в координатах (а1; a2)
Fig. 3. Stable generation regime for а/аоро = 1.5. Time series (a) of signal radiation and (b) phase portrait in coordinates (a1, a2) are shown
-12
-16 ----1-
20 25 30 35 40 45 /GHz
e
Рис. 4. Квазипериодические колебания амплитуды при а/аоро = 2.26. Показаны временные реализации (a, b) и фазовые портреты переходного процесса (с) и установившегося режима (d) (черным цветом показано сечение Пуанкаре); e - амплитудный спектр установившегося режима колебаний
Fig. 4. Quasi-periodic oscillations of the amplitude for а/аоро = 2.26. Time series (a, b) and phase portraits the transition process (c) and the steady state regime (d) (black limit cycle denotes the Poincare section); e - the amplitude spectrum of the stationary regime
Уравнение в вариациях для системы уравнений (1) имеет следующий вид:
_ дР(й1,0,2, VI, Й1т, Й2т) / | + дР(а1,а2, У1,а1т,Й2т) / . (4)
д(аг,а2, Vl) д(01т,02т) 10
Здесь Р(а1,а2, v1,а1т,а2т) - правая часть системы (1), а дР(а1,а2, v1, а1т, а2т)/д(а1,а2, v1) и дР(а1,а2, v1,a1т,a2т)/д(а1т,а2т) обозначают матрицы производных (матрицы Якоби), вычисленные по отношению к переменным а1, а2, Vl и запаздывающим переменным а1т, а2т, соответственно. Переменные с тильдой 0ь а2 и V! обозначают возмущения траектории, которые мы вычисляем. С формальной точки зрения система с запаздыванием (1) имеет бесконечно большую размерность фазового пространства, а её вектор возмущения представляет собой объединение двух отрезков траектории [а^ — т), а1 (¿)], [а2(^ — т), а2(¿)] и текущего значения переменной Vl(t), для которой нет запаздывания. При численном счёте дискретизация по времени одновременно задаёт дискретизацию вектора возмущения так, что его размерность равна 2(т/Н + 1) + 1, где Н - шаг численного счёта. Более подробное обсуждение теории и численных методов ля-пуновского анализа систем с запаздыванием можно найти, например, в работах [23,24]. Для системы (1) мы вычислили старший показатель Ляпунова для восьми произвольно выбранных точек внутри заштрихованной области О на рис. 2, а именно для о/ооро = 2.72 и 3.17 для запаздываний т = 0.1, 0.2, 0.3, 0.4. Вычисления проводились с постоянным шагом Н = 10-5. Во всех случаях показатель попал в диапазон от —7 ■ 10-4 до 2 ■ 10-4, то есть оказался достаточно малым по абсолютной величине. Это даёт основания полагать, что фактическое значение старшего показателя Ляпунова во всех восьми рассмотренных случаях равно нулю, а отклонения полученных значений от нуля следует считать численной погрешностью. Таким образом, анализ показателей Ляпунова подтверждает сделанный выше вывод о квазипериодическом характере динамики изучаемой системы при рассматриваемых значениях параметров.
Заключение
В работе представлена модель внутрирезонаторного оптического параметрического генератора с внутрирезонаторной накачкой полупроводниковым дисковым лазером. Исследование динамики системы проведено с использованием средств численного интегрирования, бифуркационного анализа систем с временным запаздыванием и расчета спектра показателей Ляпунова. Показано, что, в отличие от генератора разностной частоты, вариация управляющих параметров в исследуемой системе может приводить к потере устойчивости состояния равновесия, соответствующего стационарной генерации. Развивающаяся неустойчивость носит квазипериодический характер для всех рассмотренных случаев, лежащих в пределах физической реализуемости параметров устройства. Предметом дальнейших исследований будет изучение причин возникновения неустойчивости, а также изучение характерных бифуркационных переходов, в том числе приводящих к мультистабильности.
Библиографический список
1. Sorokina I.T., Vodopyanov K.L., eds. Solid-State Mid-Infrared Laser Sources. Ser. Topics in Applied Physics, Vol. 89. Springer-Verlag, 2003. 558 p.
2. Morozov Y.A., MorozovM.Y., Kozlovsky V.I., Okhotnikov O.G. Compact intracavity singly-resonant optical parametric oscillator pumped by GaSb-based vertical external cavity surface-emitting laser: Concept and the main operational characteristics // IEEE J. of Selected Topics in Quantum Electron. 2015.Vol. 21, no. 1. 1603105-1-1603105-5.
3. Morozov Y.A. Multi-mode dynamics of optical oscilators based on intracavity nonlinear frequency down-conversion // Appl. Phys B. 2018. Vol. 124, no. 1. 12-1-12-7.
4. Stothard D.J.M., Hopkins J.-M., Burns D., Dunn M.H. Stable, continuous-wave, intracavity, optical parametric oscillator pumped by a semiconductor disk laser (VECSEL) // Optics Express. 2009. Vol. 17, no. 13. P. 10648-10658.
5. Siegman A.E. Lasers, University Science Book, 20 Edgehill Road, Mill Valley, California, 1986.
6. Oshman M.K. and Harris S.E. Theory of optical parametric oscillation internal to the laser cavity // IEEE J. Quantum Electron. 1968. Vol. QE-4, no. 8. P. 491-502.
7. Hodges S.E., Munroe M., Cooper J., Raymer M.G. Multimode laser model with coupled cavities and quantum noise // JOSA B. 1997. Vol. 14, no. 1. P. 191-199.
8. Turnbull G.A., Dunn M.H., Ebrahimzadeh M. Continuous-wave, intracavity optical parametric oscillators: an analysis of power characteristics // Appl. Phys B. 1998. Vol. 66. P. 701-710.
9. Morozov Y.A. Transient power characteristics of a compact singly resonant intracavity optical parametric oscillator pumped by a semiconductor disk laser // JOSA B. 2016. Vol. 33, no. 7. P. 1470-1475.
10. Park J.-D., Seo D.-S., McInerney J.Self-pulsations in strongly coupled asymmetric external cavity semiconductor lasers // IEEE J. Quantum Electron. 1990. Vol. 26, no. 8. P. 1353-1362.
11. Hui R.-Q., Tao S.-P. Improved rate equations for external cavity semiconductor lasers // IEEE J. Quantum Electron. 1989. Vol. 25. P. 1580-1584.
12. van Tartwijk G.H. M., Lenstra D. Semiconductor laser with optical injection and feedback // Quantum Semiclass. Opt. 1995. Vol. 7. P. 87-143.
13. Morozov Y.A., Leinonen T., Harkonen A., Pessa M. Simultaneous dual-wavelength emission from vertical external-cavity surface-emitting laser: A numerical modeling // IEEE J. Quantum Electron. 2006. Vol. 42. P. 1055-1061.
14. Lang R., Kobayashi K. External optical feedback effects on semiconductor injection laser properties // IEEE J. Quantum Electron. 1980. Vol. 16. P. 347-355.
15. Морозов Ю.А., Балакин М.И., Кочкуров Л.А., Морозов М.Ю. Генератор разностной частоты и оптический параметрический генератор с внутрирезонаторной накачкой полупроводниковым дисковым лазером: Сопоставительный анализ в модели с запаздыванием // Изв. Сарат. ун-та. Нов. сер. Сер. Физика. 2019. Т. 19, № 1. С. 34-42.
16. Calvez S., Burns D., Dawson M.D. Optimization of an optically pumped 1.3-^m GaInNAs vertical-cavity surface-emitting laser // IEEE Phot. Techn. Lett. 2002. Vol. 14, no. 2. P. 131-133.
17. Tropper A.C., Hoogland S. Extended cavity surface-emitting semiconductor lasers // Prog. Quantum Electronics. 2006. Vol. 30. P. 1-43.
18. Engelborghs K., Luzyanina T., Samaey G. DDE-BIFTOOL. v.2.00 user manual: a Matlab package for bifurcation analysis of delay differential equations. Technical Report TW-330, Katholieke Universiteit Leuven, Leuven, Belgium (2001).
19. Benettin G., GalganiL., GiorgilliA., Strelcyn J.-M. Lyapunov characteristic exponents for smooth dynamical systems and for Hamiltonian systems: A method for computing all of them. Part 1: Theory // Meccanica. 1980. Vol. 15, no. 1. P. 9-20.
20. Shimada I., Nagashima T. A numerical approach to ergodic problem of dissipative dynamical systems // Prog. Theor. Phys. 1979. Vol. 61, no. 6. P. 1605-1616.
21. Kuptsov P., Parlitz U. Theory and computation of covariant Lyapunov vectors // Journal of Nonlinear Science. 2012. Vol. 22, no. 5. P. 727-762. DOI 10.1007/s00332-012-9126-5.
22. Купцов П.В. Вычисление показателей Ляпунова для распределенных систем: Преимуще-
ства и недостатки различных численных методов // Изв. вузов. Прикладная нелинейная динамика. 2010. Т. 18, № 5. С. 93-112.
23. Kuptsov P.V., Kuznetsov S.P. Numerical test for hyperbolicity of chaotic dynamics in time-delay systems // Phys. Rev. E. 2016. Vol. 94. 010201(R).
24. Kuptsov P. V, Kuznetsov S.P. Numerical test for hyperbolicity in chaotic systems with multiple time delays // Communications in Nonlinear Science and Numerical Simulation. 2018. Vol. 56. P. 227-239.
References
1. Sorokina I.T., Vodopyanov K.L., eds. Solid-State Mid-Infrared Laser Sources. Ser. Topics in Applied Physics, Vol. 89, Springer-Verlag, 2003, 558 p.
2. Morozov Y.A., Morozov M.Y., Kozlovsky V.I., Okhotnikov O.G. Compact intracavity singly-resonant optical parametric oscillator pumped by GaSb-based vertical external cavity surface-emitting laser: Concept and the main operational characteristics, IEEE J. of Selected Topics in Quantum Electron, 2015, Vol. 21, no. 1, 1603105-1-1603105-5.
3. Morozov Y.A. Multi-mode dynamics of optical oscilators based on intracavity nonlinear frequency down-conversion. Appl. Phys B, 2018, vol. 124, no. 1, 12-1-12-7.
4. Stothard D.J.M., Hopkins J.-M., Burns D., Dunn M.H. Stable, continuous-wave, intracavity, optical parametric oscillator pumped by a semiconductor disk laser (VECSEL). Optics Express, 2009, vol. 17, no. 13, pp. 10648-10658.
5. Siegman A.E. Lasers, University Science Book, 20 Edgehill Road, Mill Valley, California, 1986.
6. Oshman M.K. and Harris S.E. Theory of optical parametric oscillation internal to the laser cavity. IEEE J. Quantum Electron., 1968, Vol. QE-4, no. 8. pp. 491-502.
7. Hodges S.E., Munroe M., Cooper J., Raymer M.G. Multimode laser model with coupled cavities and quantum noise. JOSA B, 1997, vol. 14, no. 1, pp. 191-199.
8. Turnbull G.A., Dunn M.H., Ebrahimzadeh M. Continuous-wave, intracavity optical parametric oscillators: An analysis of power characteristics. Appl. Phys B, 1998, vol. 66, pp. 701-710.
9. Morozov Y.A. Transient power characteristics of a compact singly resonant intracavity optical parametric oscillator pumped by a semiconductor disk laser. JOSA B, 2016, vol. 33, no. 7, pp. 1470-1475.
10. Park J.-D., Seo D.-S., McInerney J. Self-pulsations in strongly coupled asymmetric external cavity semiconductor lasers. IEEE J, Quantum Electron., 1990, vol. 26, no. 8, v. 1353-1362.
11. Hui R.-Q., Tao S.-P. Improved rate equations for external cavity semiconductor lasers. IEEE J. Quantum Electron., 1989, vol. 25. pp. 1580-1584.
12. van Tartwijk G.H. M., Lenstra D. Semiconductor laser with optical injection and feedback. Quantum Semiclass. Opt., 1995, vol. 7, pp. 87-143.
13. Morozov Y.A., Leinonen T., Harkonen A., Pessa M. Simultaneous dual-wavelength emission from vertical external-cavity surface-emitting laser: A numerical modeling. IEEE J. Quantum Electron., 2006, vol. 42, pp. 1055-1061.
14. Lang R., Kobayashi K. External optical feedback effects on semiconductor injection laser properties. IEEE J. Quantum Electron., 1980, vol. 16, pp. 347-355.
15. Morozov Y.A., Balakin M.I., Kochkurov L.A., Morozov M.Y. Generator difference frequency and optical parametric generator with intracavity-pumped semiconductor disk laser: A comparative analysis of the models with a delay. Izv. of Saratov State University. New Ser. Ser. Phyzika, 2019, vol. 19, no. 1, pp. 34-42 (in Russian).
16. Calvez S., Burns D., Dawson M.D. Optimization of an optically pumped 1.3-^m GalnNAs vertical-cavity surface-emitting laser. IEEE Phot. Techn. Lett. 2002, vol. 14, no. 2, pp. 131-133.
17. Tropper A.C., Hoogland S. Extended cavity surface-emitting semiconductor lasers, Prog. Quantum Electronics, 2006, vol. 30, pp. 1-43.
18. Engelborghs K., Luzyanina T., Samaey G. DDE-BIFTOOL. v.2.00 user manual: a Matlab Package for Bifurcation Analysis of Delay Differential Equations. Technical Report TW-330, Katholieke Universiteit Leuven, Leuven, Belgium (2001).
19. Benettin G., Galgani L., Giorgilli A., Strelcyn J.-M. Lyapunov characteristic exponents for smooth dynamical systems and for Hamiltonian systems: A method for computing all of them. Part 1: Theory. Meccanica, 1980, vol. 15, no. 1, pp. 9-20.
20. Shimada I., Nagashima T. A numerical approach to ergodic problem of dissipative dynamical systems. Prog. Theor. Phys., 1979, vol. 61, no. 6, pp. 1605-1616.
21. Kuptsov P., Parlitz U. Theory and computation of covariant Lyapunov vectors. Journal of Nonlinear Science, 2012, vol. 22, no. 5, pp. 727-762, DOI 10.1007/s00332-012-9126-5.
22. Kuptsov P.V. Computation of Lyapunov exponents for spatially extended systems: Advantages and limitations of various numerical methods. Izvestiya VUZ. Applied Nonlinear Dynamics, 2010, vol. 18, no. 5, pp. 93-112 (in Russian).
23. Kuptsov P.V., Kuznetsov S.P. Numerical test for hyperbolicity of chaotic dynamics in time-delay systems. Phys. Rev. E. 2016. Vol. 94. 010201(R).
24. Kuptsov P.V., Kuznetsov S.P. Numerical test for hyperbolicity in chaotic systems with multiple time delays. Communications in Nonlinear Science and Numerical Simulation, 2018, vol. 56, pp. 227-239.
Кочкуров Леонид Алексеевич. Окончил Саратовский государственный университет (2012). Кандидат физико-математических наук (2017), доцент кафедры «Радиоэлектроника и телекоммуникации» Саратовского государственного технического университета имени Гагарина Ю.А. Научные интересы - лазерная физика, волоконная оптика, нелинейная оптика. Автор более 30 научных публикаций.
Россия, 410054 Саратов, ул. Политехническая, 77
Саратовский государственный технический университет имени Гагарина Ю.А. E-mail: lkochkurov@gmail.com
Балакин Максим Игоревич. Окончил Саратовский государственный университет (2011). Кандидат физико-математических наук (2015), доцент кафедры «Радиоэлектроника и телекоммуникации» Саратовского государственного технического университета имени Гагарина Ю.А. Научные интересы - нелинейная динамика, теория колебаний, системы с временным запаздыванием и мультистабильность. Автор более 30 научных публикаций.
Россия, 410054 Саратов, ул. Политехническая, 77
Саратовский государственный технический университет имени Гагарина Ю.А. E-mail: balakinmaxim@gmail.com
Купцов Павел Владимирович - профессор кафедры «Приборостроение» Саратовского государственного технического университета имени Гагарина Ю.А., а также сотрудник научной группы теоретической нелинейной динамики в Саратовском филиале Института радиотехники и электроники РАН. Область научных интересов - сложные колебания в распределённых системах, хаотическая динамика в системах высокой размерности, численные методы.
Россия, 410054 Саратов, ул. Политехническая, 77
Саратовский государственный технический университет имени Гагарина Ю.А. E-mail: p.kuptsov@rambler.ru
Морозов Юрий Александрович - кандидат физико-математических наук, ведущий научный сотрудник Института радиотехники и электроники им. В.А. Котельникова Российской академии наук.
Россия, 410019 Саратов, ул. Зеленая, 38
Саратовский филиал Института радиотехники и электроники им. В.А. Котельникова РАН E-mail: yuri.mor@rambler.ru
