Например, Бобцов

МОДЕЛИРОВАНИЕ НЕУСТОЙЧИВОГО РЕЗОНАТОРА МОНОИМПУЛЬСНОГО Nd:YAG-ЛАЗЕРА С ГРАДИЕНТНЫМ ВЫХОДНЫМ ЗЕРКАЛОМ

И.Н. Дубинкин, В.В. Назаров, Л.В. Хлопонин, В.Ю. Храмов

УДК 621.373.826
МОДЕЛИРОВАНИЕ НЕУСТОЙЧИВОГО РЕЗОНАТОРА МОНОИМПУЛЬСНОГО Nd:YAG-ЛАЗЕРА
С ГРАДИЕНТНЫМ ВЫХОДНЫМ ЗЕРКАЛОМ
И.Н. Дубинкин, В.В. Назаров, Л.В. Хлопонин, В.Ю. Храмов
Представлен метод моделирования генерации лазерного излучения в неустойчивом радиально-симметричном резонаторе с выходным градиентным зеркалом. Приведены результаты исследований влияния параметров градиентного зеркала на пространственные и энергетические характеристики лазерного излучения. Проведено моделирование лазерной генерации в режимах активной и пассивной модуляции добротности. Проведен сравнительный анализ пространственных и энергетических характеристик излучения для указанных режимов. Ключевые слова: моделирование, неустойчивый резонатор, градиентное зеркало, модуляция добротности.
Введение
Среди разрабатываемых в настоящее время лазеров твердотельные лазеры, в частности Nd:YAGлазеры, занимают одно из лидирующих мест по объему производства. Мощные импульсные Nd:YAGлазеры могут использоваться в различных технологических процессах благодаря компактности и надежности конструкции, стабильности параметров и возможности применения оптических волокон для доставки излучения до объекта воздействия. Они нашли широкое применение для обработки материалов, в медицине, научных исследованиях, лазерных системах с преобразованием частоты излучения.
Проблеме распространения лазерного излучения в открытых неустойчивых резонаторах, начиная с работы Сигмана [1], посвящено большое количество теоретических и экспериментальных работ, основные результаты которых подытожены в известных монографиях [2, 3]. Математические модели и основные численные методы решения указанной проблемы изложены в [4–8]. Однако предыдущие исследования [9, 10], имеющие непосредственное отношение к тематике настоящей работы, как правило, не учитывали такие факторы, как динамика поля лазерного излучения внутри резонатора, дифракционные эффекты, влияние усиления излучения в активной среде.
Математическая модель моноимпульсного лазера
Одним из наиболее эффективных методов решения задачи о лазерной генерации в случае радиально симметричного резонатора является метод, основанный на быстром преобразовании Ганкеля [1].

Научно-технический вестник информационных технологий, механики и оптики, 2013, № 3 (85)

103

МОДЕЛИРОВАНИЕ НЕУСТОЙЧИВОГО РЕЗОНАТОРА МОНОИМПУЛЬСНОГО …

Кратко рассмотрим разработанный авторами алгоритм, реализующий математическую модель, основанную на указанном методе.
Первый шаг алгоритма заключался в том, что на основе быстрого преобразования Ганкеля был реализован метод Фокса и Ли для пустого резонатора, содержащего тонкую линзу, эквивалентную тепловой линзе активного элемента, а также градиентное выходное зеркало. В результате было получено распределение комплексной амплитуды поля, соответствующее суперпозиции собственных мод резонатора, что позволило определить потери для исследуемого резонатора. Указанное распределение задает стартовое «шумовое» поле, которое будет усиливаться в активной среде в процессе лазерной генерации.
На втором шаге был реализован метод Фокса и Ли для резонатора, содержащего активную среду и насыщающийся поглотитель. Начальным распределением поля здесь служило распределение, полученное на предыдущем шаге для пустого резонатора. Для моделирования процесса взаимодействия излучения внутри резонатора с активной средой и насыщающимся поглотителем были использованы уравнения Франца–Нодвика [6]. Стоит отметить, что на данном шаге критерием прекращения вычислений являлось снижение интенсивности излучения на несколько порядков относительно пикового значения.
В результате указанных процедур с шагом дискретизации по времени, равным времени обхода резонатора, определялось распределение интенсивности поля на выходном зеркале резонатора в ходе генерации. Полученные данные использовались для определения энергетических и пространственных характеристик импульса излучения лазерной генерации с учетом характеристик выходного зеркала резонатора.
Основными выражениями при реализации метода Фокса и Ли были уравнения, полученные в работах [1, 4–6] и преобразованные в соответствии с рассматриваемой задачей. Эти уравнения можно записать в следующем виде:

  

fn



r0

exp n  E1 r0

exp n  exp 

i B

Ar02

exp2n ;

 

gn



m







2 B

r00

exp n  m  J0

 2  B

r00

exp n



m



 

;

(1)

 



IFFT

 IFFT



fn





FFT



gnm



;

  

E2



0



exp



m







i



0 exp
0 exp 

m m







exp

 



i B

D02

exp



2m



,

где n, m – целочисленные индексы, соответствующие узлам расчетной сетки; A, B, D – элементы лучевой

матрицы оптической системы; E1 – начальное распределение комплексной амплитуды поля на опорной плоскости; E2 – искомое распределение комплексной амплитуды поля; λ – длина волны излучения; J0 – функция Бесселя нулевого порядка; r0, ρ0,  – числовые параметры, определяемые в соответствии с [4]; FFT, IFFT – операторы прямого и обратного быстрого преобразования Фурье соответственно.

Основными выражениями при моделировании процесса взаимодействия лазерного излучения с

активной средой и насыщающимся поглотителем были уравнения, полученные в [7] и модифицирован-

ные в соответствии с решаемой задачей:



U 


i

AE

,

SH



 



exp

2

AE

U


i

,SH 

AE,SH   exp  AE,SH

c n AE,SH 

I

i



 



exp



U


i

AE

,

SH



L

AE

, SH





AE

, SH

U


i

AE

,

SH

  L

AE

, SH







1

;

 


Ii





 1

Ii

1



exp





AE

,

SH

U


i

AE

, SH



L

AE

,

SH







exp

 

2

AE

, SH



c n AE,SH 

Ii 

,

(2)

где i – целочисленные индексы, соответствующие узлам расчетной сетки; U(AE,SH) – разность населенно-

стей рабочих уровней в активной среде или насыщающемся поглотителе; σ(AE,SH) – сечение рабочего пе-

рехода в активной среде или насыщающемся поглотителе; L(AE,SH) – характерный размер активной среды

или насыщающегося поглотителя; n(AE,SH) – показатель преломления активной среды или насыщающегося

поглотителя; τ – длительность импульса излучения постоянной интенсивности на входе в среду; Ii – интенсивность излучения в i-ой точке расчетной сетки; с – скорость света в вакууме.

Заметим, что уравнения (2) являются частным решением задачи о взаимодействии импульса излу-

чения с двухуровневой средой и применимы только в случае, если со средой взаимодействует импульс

постоянной во времени интенсивности. По этой причине в рамках данной модели усиливаемый импульс

представлялся последовательностью постоянных во времени элементарных импульсов.

Для определения пространственных характеристик импульса излучения лазерной генерации был

проведен анализ распределения плотности энергии в фокальной плоскости модельной тонкой линзы,

размещенной на выходе лазерного резонатора. С помощью выражений (1) по уровню 90% энергии опре-

104

Научно-технический вестник информационных технологий, механики и оптики, 2013, № 3 (85)

И.Н. Дубинкин, В.В. Назаров, Л.В. Хлопонин, В.Ю. Храмов

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

Результаты исследования характеристик лазерного излучения

На основе разработанного алгоритма было проведено численное исследование генерации излучения в неустойчивом резонаторе, схема которого показана на рис. 1. С целью упрощения задачи оптимизации геометрические параметры резонатора были зафиксированы.

R1 = 5000

R2 = –1000

М1 А

F АЭ

М2 Д

8 8

150 100 150 Б

Рис. 1. Геометрия исследуемого резонатора: АЭ – активный элемент; M1, M2 – зеркала резонатора; А – блок пассивной модуляции добротности; Б – блок активной модуляции добротности; Д – диафрагма; F – тепловая линза. Оптическая сила тепловой линзы 0,5 м–1. Все линейные величины приведены в мм

Для данной геометрии коэффициент увеличения резонатора составляет 1,6, а эквивалентное число Френеля равно 7,5. В качестве насыщающегося поглотителя в режиме пассивной модуляции добротности был использован кристалл GSGG:Cr:Mg.
Профиль коэффициента отражения выходного градиентного зеркала определяется выражением

R



R0

exp

  

2

rn nm

  

,

где R0 – максимальный коэффициент отражения; r – радиальная координата, ωm – радиус (характерный размер) градиентного зеркала; n – показатель степени супергауссова профиля коэффициента отражения

выходного зеркала.

E, мДж

1000

900 1
800 2
700

600 1 2 3 4 5 6 7 8 9 10 n1,1отн.12ед.
Рис. 2. Зависимость энергии импульса излучения от показателя степени супергауссова профиля коэффициента отражения выходного зеркала: 1 – активная модуляция добротности;
2 – пассивная модуляция добротности с использованием насыщающегося поглотителя с начальным пропусканием 60%. Параметры градиентного зеркала: максимальный коэффициент отражения – 0,25;
радиус градиентного зеркала – 1,35 мм

Научно-технический вестник информационных технологий, механики и оптики, 2013, № 3 (85)

105

МОДЕЛИРОВАНИЕ НЕУСТОЙЧИВОГО РЕЗОНАТОРА МОНОИМПУЛЬСНОГО …

Параметры градиентного зеркала Максимальный коэффициент отражения Радиус (характерный размер) градиентного зеркала, мм Показатель степени супергауссова профиля коэффициента отражения выходного зеркала
Параметры насыщающегося поглотителя Начальный коэффициент поглощения, см–1 Характерный линейный размер, см Концентрация активатора, 1018 см–3

0,25–0,3 1,2–1,5
6–8
1,0–1,5 0,1–0,3 5,5–8,5

Таблица 1. Характеристики насыщающегося поглотителя и градиентного зеркала

6I,0МВт/см2 I, МВт/см2
5500 4400
3300

11 2
33

2200

1100

0
00,00

00,,0055 00,,1100 00,,1155

00,2,200 00,,2255

00,,3300 00,,3355 r, сr,м0с,4м0

Рис. 3. Поперечное распределение интенсивности в лазерном пучке: 1 – активная модуляция добротности; 2 – пассивная модуляция добротности с использованием насыщающегося поглотителя
с начальным пропусканием 80%; 3 – пассивная модуляция добротности с использованием насыщающегося поглотителя с начальным пропусканием 60%. Параметры градиентного зеркала: минимальный коэффициент отражения – 0,25; радиус градиентного зеркала – 1,35 мм; показатель степени
супергауссова профиля коэффициента отражения выходного зеркала – 6

Результаты оптимизации параметров градиентного зеркала и насыщающегося поглотителя с це-

лью достижения более однородного поперечного распределения интенсивности в импульсе лазерного

излучения и повышения энергии импульса излучения в режиме пассивной модуляции добротности при-

ведены в табл. 1. Зависимость энергии импульса излучения от показателя степени супергауссова профи-

ля коэффициента отражения выходного зеркала показана на рис. 2. Дальнейшие результаты приведены

для указанных в табл. 1 значений параметров поглотителя и зеркала.

Поперечное распределение плотности мощности в импульсе лазерного излучения для случаев ак-

тивной модуляции добротности и пассивной модуляции добротности с использованием насыщающихся

поглотителей с различным начальным пропусканием показано на рис. 3. Результаты приведены для мо-

мента достижения пиковой мощности в импульсе генерации. Из приведенных на рис. 3 данных очевидно,

что в режиме пассивной модуляции добротности повышается однородность лазерного пучка.

Оценка неоднородности поперечного распределения интенсивности лазерного пучка  проводи-

лась при помощи выражения





Imax  I0 I max

,

где I0 – минимальное значение интенсивности в приосевой области; Imax – максимальное значение интенсивности лазерного пучка.

Основные характеристики излучения для различных режимов модуляции приведены в табл. 2.

Следует отметить, что для оценки значения коэффициента усиления слабого сигнала в активной среде,

исходя из энергии накачки, были использованы выражения из работы [10].

106

Научно-технический вестник информационных технологий, механики и оптики, 2013, № 3 (85)

И.Н. Дубинкин, В.В. Назаров, Л.В. Хлопонин, В.Ю. Храмов

Параметр
Энергия импульса генерации при энергии импульса накачки 50 Дж, мДж Длительность импульса излучения по полувысоте, нс Угловая расходимость излучения, мрад Угловая расходимость излучения, мрад Неоднородность по поперечному сечению, %

Значение Активная модуляция Пассивная модуляция

1026

814

30 50 1,8 1,5

18 11

Таблица 2. Характеристики излучения лазерной генерации

Заключение

В работе приведено описание разработанной математической модели, позволяющей исследовать генерацию лазерного излучения в неустойчивом резонаторе с выходным градиентным зеркалом. Была проведена оптимизация параметров градиентного зеркала с целью достижения более однородного поперечного распределения плотности мощности излучения и максимальной энергии импульса генерации. Представлены результаты численного моделирования параметров лазерного излучения для режимов активной и пассивной модуляции добротности. Сравнительный анализ указанных режимов показал, что использование метода пассивной модуляции добротности позволяет повысить однородность поперечного распределения интенсивности в импульсе лазерного излучения и снизить угловую расходимость излучения. Проведен анализ влияния параметров градиентного зеркала и насыщающегося поглотителя на характеристики лазерного излучения.

Литература

1. Siegman A.E., Miller H.Y. Unstable optical resonator loss calculations using the prony method // Applied Optics. – 1970. – V. 9. – № 12. – P. 2729–2763.
2. Ананьев Ю.А. Оптические резонаторы и проблемы расходимости лазерного излучения. – М.: Наука, 1979. – 328 с.
3. Ананьев Ю.А. Оптические резонаторы и лазерные пучки. – М.: Наука, 1990. – 264 с. 4. Murphy W.D., Bernabe M.L. Numerical procedures for solving nonsymmetrical eigenvalue problems associ-
ated with optical resonators // Applied Optics. – 1978. – V. 17. – № 15. – P. 2358–2365. 5. Карамзин Ю.Н., Конев Ю.Б. Численные методы исследования работы неустойчивых телескопических
резонаторов с учетом дифракции и эффекта насыщения в активной среде. – Препринт ИПМ им. М.В. Келдыша АН СССР, 1974. – № 57. – 14 с. 6. Елкин Н.Н., Напартович А.П. Прикладная оптика лазеров. – М.: ЦНИИатоминформ, 1989. – 183 с. 7. Lee M. Frantz, John S. Nodvik. Theory of Pulse Propagation in a Laser Amplifier // J. Appl. Phys. – 1963. – V. 34. – P. 2346–2349. 8. Аладов А.В., Беззубик В.В., Белашенков Н.Р., Карасев В.Б., Назаров В.В. и др. Применение зеркал с переменным коэффициентом отражения в компактных твердотельных лазерных системах // Изв. вузов. Приборостроение. – 1998. – Т. 41. – № 3. – С. 53–56. 9. Morin M. Graded reflectivity mirror unstable laser resonators // Optical and Quantum Electronics. – 1997. – V. 29. – № 8. – P. 819–866. 10. Зверев Г.М., Голяев Ю.Д., Шалаев Е.А., Шокин А.А. Лазеры на алюмоиттриевом гранате с неодимом. – М.: Радио и связь, 1985. – 144 с.

Дубинкин Илья Николаевич Назаров Вячеслав Валерьевич Хлопонин Леонид Викторович Храмов Валерий Юрьевич

– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, студент, ilya.dubinkin@mail.ru
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кандидат технических наук, ст. научный сотрудник, lab255@grv.ifmo.ru
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кандидат технических наук, ст. научный сотрудник, l_khloponin@yahoo.com
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, доктор технических наук, профессор, зав. кафедрой, khramov@grv.ifmo.ru

Научно-технический вестник информационных технологий, механики и оптики, 2013, № 3 (85)

107