Аннотация результатов, полученных в 2020 году
1. В ходе проекта была разработана модель усилителя с учетом тепловых эффектов и реального распространения гауссова пучка. Модель основа на решении системы уравнений для населенности верхнего уровня, изменения интенсивности излучения накачки и сигнальной волны. Для решения этих уравнений поочередно были учтены следующие моменты: 1. Учет гауссового распределения интенсивности накачки и сигнальной волны 2. Учет распределения температуры в кристалле при использовании металлических оправ с теплопроводящим слоем 3. Учет зависимости времени жизни на верхнем рабочем уровне от температуры 4. Учет изменения показателя преломления от температуры по всему объему кристалла 5. Учет каустики распространения излучения накачки и сигнальной волны внутри кристалла 6. Учет дисперсии материала для излучения накачки и сигнальной волны по всему объему кристалла 2. Найдены предельные характеристики усилителя с имеющимся в лаборатории оборудованием и получены результаты исследования поведения кристалла Cr:CdSe в процессе усиления на основе разработанной модели. Были получены предельные характеристики усилителя для длины волны 2,65 мкм при использовании накачки с помощью тулиевого волоконного лазера с максимальной мощностью 20 Вт. За выходные параметры лазера, излучение которого будет усиливаться, были взять данные из [Tarabrin, Mikhail K., et al. ”High-efficiency continuous-wave single-mode room-temperature operation of Cr2+:CdSe single-crystal laser with output power of 2.3 W.” Opt. Express 27, 12090-12099 (2019)], полученные руководителем настоящего проекта, как наилучший результат в мире. Из полученной зависимости выходной мощности излучения на длине волны 2,65 мкм от мощности накачки предельной выходной мощностью данной системы является величина 4,9 Вт, которая достигается при накачке 14 Вт. После данной величины происходит загиб, который обусловлен доминированием процесса уменьшения времени жизни и ухудшения перекрытия каустик накачки и сигнала над самим процессом усиления. При этом среда нагревается до максимальной температуры 412 К, а время жизни в свою очередь падает до 0,09 мкс, то есть фактически часть кристалла вблизи передней грани не вносит вклад в усиление. 3. Найдено предельное значение эффективности системы, ширина спектра усиления и абсолютных значений мощности на выходе усилителя на основе экспериментально полученных данных по разработке мастер осциллятора. Предельная эффективность системы будет наблюдаться при использовании излучения на длине волны 2,65 мкм, так как она соответствует максимуму усиления среды. При этом ее величина на основе моделирования составила 44,7% по поглощенной мощности. Для получения ширины спектра усиления было проведено моделирования для каждой длины волны, при этом изменялись параметры сечения усиления, а также другие оптические параметры, которые зависят от длины волны. Мощность излучения на каждой длине волны в рамках данного моделирования была равной 1 Вт, но в будущем данные будут взяты с перестраиваемого лазера, излучение которого планируется усилить во время второго года осуществления проекта. Усилитель позволяет осуществлять работу во всем спектре люминесценции активной среды Cr:CdSe. Полученная зависимость коэффициента усиления от длины волны показывает возможность усиливать лазерное излучение от 2 до 3,4 мкм с помощью активной среды Cr:CdSe. На основе данных по расчету коэффициента усиления активной среды Cr:CdSe были получены значения абсолютной мощности на выходе усилителя при тех же условиях. При подаче 1 Вт сигнальной волны и 3 Вт накачки происходит увеличение мощности в максимуме усиления в 1,5 раза, на других длинах волн усиление имеет меньшую величину. 4. В ходе выполнения проекта были разработаны конструкции двух лабораторных установок, которые являются системами мастер осциллятор – усилитель мощности (MOPA). Первая установка представляет собой лазер с неселективным резонатором для генерации на 2,65 мкм для проверки концепции лазерного усилителя на базе кристалла Cr:CdSe. В ходе проекта были подобраны оптимальные параметры кристалла для мастер осциллятора и конфигурация его резонатора. Для корректного расчета системы фокусировки части излучения накачки и генерации мастер осциллятора для оптимального перекрытия их каустик в кристалле усилителя необходимо было измерить геометрические параметры излучения на выходе мастер осциллятора и накачки и их качество пучка. По результатам проведенных измерений качество пучка мастер осциллятора оказалось М2=1,045, а для излучения накачки М2=1,027. Для реализация широкополосного усиления в рамках проекта был разработан второй мастер осциллятор с плавной перестройкой по длине волны. Выходные характеристики данного осциллятора превосходят существующие в мире лазеры на основе кристалла Cr:CdSe. 5. Дополнительно были выполнены работы по исследованию параметров высокоскоростных приемников среднего ИК-диапазона на длине волны 1,9 мкм с помощью разработанного в лаборатории лазера ультракоротких импульсов на основе тулиевого активного волокна и гибридной синхронизации мод. Дополнительно была разработана методика просветления поликристаллических волокон на основе AgClBr для передачи излучения среднего ИК-диапазона на выходе разрабатываемого в ходе проекта усилителя до места его использования для медицинских исследований.

Аннотация результатов, полученных в 2021 году
В ходе второго года выполнения проекта была проведена оптимизация разработанной математической модели лазерного усилителя на предмет быстродействия и апробированы основные инструменты увеличения скорости вычисления, такие как анализ эффективности встроенных решателей дифференциальных уравнений в MATLAB, создание и использование структуры параметров, предварительное выделение памяти для массивов, использование параллельных вычислений. Расчеты по исходной версии математической модели при точности вычислений с шагом 4 мкм по длине кристалла и 0,75 мкм по радиусу занимали 106,5 минут, что ограничивало количество возможных итераций в течении суток. По результатам профилирования был проведен анализ скорости вычислений каждого типа решателя ОДУ и с учетом влияния погрешности вычислений выбран решатель, основанный на методе трапеций с интерполяцией. Данное решение позволило ускорить расчеты на 37,48 минут. Использование структур параметров не показало явного ускорения - увеличилось время вызова функции. Данный метод стоит использовать в случае работы с матрицами, содержащими разные типы данных. Предварительное выделение памяти также увеличило общее время выполнения расчётов, но уменьшило процессорное время, затраченного на выполнение основных функций. Самым эффективным способом оптимизации показал себя метод параллельного расчета и использования многоядерных процессоров. Данный способ ускорил время вычисления программы 35,6 раза, что позволило значительно увеличить эффективность работы математической модели. Для верификации математической модели были разработаны и созданы две конфигурации экспериментальных макетов - работающий на длине волны максимума усиления среды и широкополосный усилитель, работающий в диапазоне 2,56-2,92 мкм. Макет, работающий на максимуме усиления среды, состоит из двух основных блоков – блок лазера, являющимся источником для последующего усиления сигнала и блок однопроходного усилителя. В разработанной схеме в качестве источника излучения использовался лазер с активной средой в виде кристалла Cr2+:CdSe, работающий на длине волны 2650 нм. Усилитель лазерного излучения также был разработан на основе кристалла Cr2+:CdSe и выполнен в конфигурации торцевой накачки для лучшего пространственного перекрытия мод. В качестве источников накачки использовались волоконные тулиевые лазеры TLM-20 и TLM-10 на длине волны 1908 нм. Также были проведены предварительные расчеты геометрии распространения пучка и параметров оптических элементов с помощью программного обеспечения ReZonator. Экспериментальный макет широкополосного усилителя лазерного излучения среднего ИК-диапазона спроектирован также на кристалле Cr2+:CdSe. В качестве источника излучения использовался перестраиваемый по длине волны лазер работающий на кристалле Cr2+:ZnSe в диапазоне 2,56-2,92 мкм. Перестройка лазера осуществляется с помощью интерференционно-поляризационного фильтра Лио. Преимущество данной схемы заключается в неизменности направления выходного излучения лазера, что позволяет создать полностью автоматизированный широкополосный усилитель с торцевой накачкой. В качестве источников накачки также использовались тулиевые волоконные лазеры TLM-20 и TLM-10 на длине волны 1908 нм. С помощью разработанных макетов были получены значения M^2 параметра по методу ножа Фуко. В соответствии с полученными результатами были рассчитаны значения параметра M^2: 2650 нм – 2,3, 2700 нм – 4,2, 2800 нм – 2,1, 2920 нм – 1,95. Из полученных данных видно, что пучок имеет высокое значение M^2 параметра. Это может быть объяснено отсутствием просветляющих покрытий на гранях кристалла, что приводит к возникновению переотражений и искажает профиль пучка, а также сильными тепловыми эффектами, наводящимися в кристалле. Улучшить качество пучка поможет просветление поверхностей кристалла, улучшение параметров охлаждения и расчёт оптимальных характеристик такие как габариты кристалла и степень легирования. Кроме этого, были получены значения максимальной выходной мощности и проведена оценка эффективности разработанного усилителя лазерного излучения на различных длинах волн. Для этого были использованы ранее полученные данные зависимости мощности генерации от падающей мощности накачки. Получены следующие значения эффективности для различных длин волн: 2650 нм – 55%, 2700 нм – 40%, 2800 нм – 29%, 2920 нм – 6%. Полученные данные хорошо согласуются с математическим расчетами что подтверждает универсальность разработанной модели для любых типов четырехуровневых сред и возможность рассчитать оптимальные параметры системы.