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

МОДЕЛЬ МНОГОЧАСТОТНОЙ 3 МКМ-ГЕНЕРАЦИИ ИЗЛУЧЕНИЯ ЛАЗЕРОВ НА ЭРБИЕВЫХ КРИСТАЛЛАХ С ДИОДНОЙ НАКАЧКОЙ

25
УДК 535.374:621.375.8

М. В. ИНОЧКИН, В. В. НАЗАРОВ, Д. Ю. САЧКОВ, Л. В. ХЛОПОНИН, В. Ю. ХРАМОВ
МОДЕЛЬ МНОГОЧАСТОТНОЙ 3 МКМ-ГЕНЕРАЦИИ ИЗЛУЧЕНИЯ ЛАЗЕРОВ НА ЭРБИЕВЫХ КРИСТАЛЛАХ С ДИОДНОЙ НАКАЧКОЙ

Теоретически исследованы параметры свободной генерации излучения в области длин волн 3 мкм Er:YLF-, Er:YAG- и Er:YSGG-лазеров при селективном возбуждении верхнего лазерного уровня перехода 4I11/2—4I13/2 одиночными импульсами накачки.

Ключевые слова: эрбиевый кристалл, селективная накачка, диодная накачка, многочастотная генерация.

Введение. Задачи по созданию источников лазерного излучения 3 мкм-диапазона не теряют

своей актуальности в связи с максимальным коэффициентом поглощения биологических тканей в

данной области спектра. Особый интерес представляют лазеры на эрбиевых кристаллах, накачи-

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

рацию в различных режимах. В рамках настоящей статьи теоретически исследованы параметры

излучения наиболее распространенных эрбиевых Er:YLF-, Er:YAG- и Er:YSGG-лазеров в режиме

свободной генерации при одиночных импульсах диодной накачки. Теоретическая модель. Для описания процессов 3 мкм-генерации (переход 4I11/2—4I13/2)
излучения лазера на эрбиевом кристалле с селективной накачкой была использована матема-

тическая модель, основанная на следующей системе балансных уравнений:

∑ ( )dN2
dt

=



i

Si

σe (λi ) N2 − σa (λi ) N1

+ R − ( A20 + A21 + W21 ) N2 + γ1N12 − γ2 N22;⎫⎪


∑ ( )dN1
dt

=

i

Si

σe (λi ) N2 − σa (λi ) N1

− A10N1 + ( A21 + W21 ) N2 − γ1N12;

⎪ ⎪ ⎪⎪

dS1 dt

=

VcµS1

(σe

( λ1 )

N2



σa

( λ1 )

N1



δ ( λ1 ) )

+

N2q;

⎬ ⎪ ⎪

.....................................................................................



dSi dt

= VcµSi

(σe (λi ) N2

− σa

(λi ) N1 − δ(λi )) +

N2q,

⎪ ⎪ ⎭⎪

(1)

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 9

26 М. В. Иночкин, В. В. Назаров, Д. Ю. Сачков и др.

где N2 — населенность верхнего лазерного уровня 4I11/2; N1 — населенность нижнего лазерного уровня 4I13/2; Amn — вероятность спонтанного излучательного перехода между уровнями m и n; W21 — вероятность безызлучательного перехода между уровнями 4I11/2 и 4I13/2; Si — плотность потока фотонов в резонаторе на длине волны λi; σe(λi) — сечение вынужденного излучения на длине волны λi; σa(λi) — сечение поглощения излучения из возбужденного состояния на длине волны λi; R — скорость накачки на уровень 4I11/2; γk — коэффициенты, характеризующие скорости процессов апконверсии; δ(λi) — показатель полных потерь излучения в резонаторе на длине волны λi; q — коэффициент, учитывающий „вклад“ спонтанного излучения в плотность потока фотонов в резонаторе; Vc=с/nar — скорость света в активном элементе, здесь nar — показатель преломления активного элемента; µ — коэффициент заполнения резонатора активной средой.
В рассматриваемой модели учитываются только два процесса апконверсии, возникающие вследствие взаимодействия между собой ионов эрбия Er3+, находящихся в состояниях 4I13/2 (γ1) и 4I11/2 (γ2). В ходе взаимодействия ионов, возбужденных в состояние 4I13/2, один из них переходит в основное состояние 4I15/2, а второй — через короткоживущее состояние 4I9/2 — на верхний лазерный уровень 4I11/2. Полное описание процесса миграции энергии в ходе апконверсии с верхнего лазерного уровня 4I11/2 требует дополнительно принимать во внимание населенности нескольких возбужденных состояний, лежащих выше верхнего лазерного уровня, а также коэффициенты ветвления переходов с этих уровней. В рамках данной модели предлагается принять, что в результате взаимодействия оба иона переходят из состояния 4I11/2 в основное состояние 4I15/2.
При расчете коэффициента заполнения резонатора µ учитывается, что активный эле-
мент заполняет не весь его объем, таким образом, µ = narlar / (L + (nar −1)lar ) , где lar — длина
активного элемента, L — длина резонатора. В расчетах принимается, что в течение импульса накачки ее скорость R постоянна по
всему объему активного элемента; значение R вычисляется по формуле
R = (λ p Ppη) / (sarlar hc),

где Pp — мощность излучения накачки, η — эффективность осветительной системы, sar — площадь поперечного сечения активного элемента, h — постоянная Планка, с — скорость

света.

В рамках данного приближения спектральные параметры излучения накачки определя-

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

излучения накачки в активном элементе.

Потери излучения в резонаторе рассчитывались по формуле

δ(λ)

=

δ0

(λ)



ln

(ρoc (λ))
2lar

,

(2)

где ρoc — коэффициент отражения выходного зеркала. Для режима свободной генерации одиночных импульсов, плотности поглощенной им-
пульсной мощности накачки до 15 кВт/см3 и длительности импульсов накачки до 5 мс (что

соответствует возбуждению не более 5 % атомов активной среды) расчеты в соответствии с

системой уравнений (1) дают расхождение не более 1 % с системой балансных уравнений и

аналитической моделью, использованными ранее [1]. Система уравнений (1) решалась чис-

ленными методами в программном пакете MathCAD.

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

чения Er:YLF-, Er:YAG- и Er:YSGG-лазеров. Рассматривались как случай многочастотной

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

волн, соответствующих основным максимумам сечения вынужденного излучения перехода

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 9

Модель многочастотной 3 мкм-генерации излучения лазеров

27

4I11/2—4I13/2 в данном кристалле, так и случай одночастотной генерации, когда предполагается наличие излучения только на одной определенной длине волны. Последний случай соответствует методам спектральной селекции излучения.
Для одночастотной генерации принималось, что она осуществляется на наиболее сильной линии перехода 4I11/2—4I13/2. В качестве таких линий выбраны известные из экспериментов: λ = 2,81 мкм для Er:YLF, λ = 2,94 мкм для Er:YAG и λ = 2,79 мкм для Er:YSGG.
При рассмотрении многочастотной генерации использованы зависимости σe(λ), полученные в работах [2, 3, 4] для кристаллов Er:YLF, Er:YAG и Er:YSGG соответственно. При этом
значения сечений поглощения σa(λ) рассчитаны исходя из зависимости σe(λ) по формуле

σa

(λ)

=

σe

(λ)

1

− F (λ) F (λ)

,

где

F (λ)

=

σe (λ) σe(λ) + σa (λ)

=

1+

Zu Zl

exp

⎛ ⎜⎝⎜

1
hc kT

⎛ ⎜ ⎝

1 λ



1 λ0

⎞⎞

⎟ ⎠

⎟⎠⎟

,

где Zu и Zl — числа заполнения для верхнего (4I11/2) и нижнего (4I13/2) лазерных уровней [5]; k — постоянная Больцмана; T — температура активной среды; λ0 — длина волны перехода между нижними подуровнями штарковской структуры уровней 4I11/2 и 4I13/2; значения F(λ) рассчитывались по данным об энергиях подуровней, приведенным в работе [6] для Er:YLF и Er:YAG и в работе [7] для Er:YSGG.
Отметим, что для кристалла Er:YSGG не удалось найти работ, в которых зависимость
σe(λ) была бы приведена в абсолютных величинах. Поэтому спектр σe(λ), рассмотренный в работе [4], был оцифрован и нормирован к сечению σe=6,5⋅10–21 см2 для длины волны 2,79 мкм.
Параметры активных сред, использованные в расчетах, приведены в табл. 1 и 2.

Параметр
A10, с–1 A21, с–1 W21, с–1 Время жизни верхнего лазерного уровня, мс
λ, мкм σe⋅10–20, см2 σa⋅10–20, см2 γ1⋅10–17, см3/с γ2⋅10–17, см3/с Концентрация ионов эрбия, ат. %

Er:YAG 200 17 [6] 9855
0,1 2,94 0,45 0,1 1,5 1,0
15

Er:YSGG 294 20 [6] 136
1,3 2,79 0,65 0,37 0,7 2,3
15

Таблица 1 Er:YLF
113 17 [8]
90
4 2,81 0,93 0,53 2,0 0,1
15

Таблица 2

Параметр

Er:YAG

Er:YSGG

Er:YLF

λ, мкм 2,69 2,79 2,83 2,87 2,94 2,79 2,83 2,87 2,93 2,66 2,71 2,81 2,84 2,85

σe(λ)

0,98 0,54 0,54 0,39 0,45 0,65 0,30 0,15 0,08 2,33 1,51 0,93 0,40 0,31

σa(λ)

1,04 0,31 0,26 0,17 0,1 0,37 0,15 0,06 0,02 3,49 1,62 0,53 0,19 0,14

Во всех расчетах принималось, что активный элемент имеет форму цилиндра диаметром 2 мм и длиной 40 мм, длина резонатора 60 мм, неактивные потери δ0=0,005 см–1, излучение как накачки, так и генерации заполняет весь объем активного элемента.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 9

28 М. В. Иночкин, В. В. Назаров, Д. Ю. Сачков и др.

Результаты расчетов и их обсуждение. Рассчитанные зависимости энергии выходного

излучения Eg рассматриваемых лазеров от поглощенной энергии импульса накачки Eр для

режима одиночных импульсов длительностью τр, равной 150 и 1000 мкс, приведены на

рис. 1, а, б соответственно. Для плотности поглощенной импульсной мощности накачки

Qp ≈ 15 кВт/см3 и τр = 150 мкс рассчитанные дифференциальные КПД генерации излучения

составили 8,5, 8 и 4 % для Er:YLF-, Er:YSGG- и Er:YAG-лазеров соответственно. Как видно,

наибольшими дифференциальными КПД обладают Er:YSGG- и Er:YLF-лазеры. Сравнивая

абсолютные КПД, отметим близкие значения данного параметра для Er:YSGG- и Er:YLF-

лазеров при существенно меньшей эффективности Er:YAG-лазера, что, очевидно, является

следствием малого времени жизни верхнего лазерного уровня в данной среде.

а) τр = 150 мкс

б) τр = 1000 мкс

Eg, Дж 0,06

YLF YAG YSGG

Eg, Дж 0,8

YLF YAG YSGG

0,6
0,04 0,4

0,02 0,2

00

0,1 0,2 0,3 0,4 0,5 Eр, Дж

1 2 3 Eр, Дж

Рис. 1

При исследовании процессов многочастотной генерации расчеты проводились при по-

стоянных по спектру потерях излучения в резонаторе (см. формулу (2)) и высоких значениях поглощенной плотности мощности накачки (Qp = 10…15 кВт/см3 и соответственно R=4,9…7,4⋅1022 с–1).

Генерация излучения Er:YLF-лазера происходит последовательно (в течение импульса

накачки) на длинах волн 2,66, 2,71, 2,81, 2,84 и 2,85 мкм (что в целом совпадает с получен-

ными ранее теоретическими и экспериментальными результатами [1]). Следует отметить, что

изменение параметров апконверсии γ1 и γ2 оказывает заметно большее влияние на параметры генерации (задержка начала генерации, энергия генерируемого импульса) на длинах волн 2,84

и 2,85 мкм, чем на длинах волн 2,66, 2,71 и 2,81 мкм. В табл. 3 приведены рассчитанные зна-

чения задержек (Тg) начала генерации (отсчет ведется от начала импульса накачки) для Qp = =10 кВт/см3 и τр = 3 мс при двух вариантах выбора параметров апконверсии: 1) γ1=4⋅10–17 см3/с, γ2= 1⋅10–17 см3/с; 2) γ1= 2⋅10–17 см3/с, γ2= 1⋅10–18 см3/с.
Таблица 3

λ, мкс

Тg, мкс, — 1-й вариант

Тg, мкс, — 2-й вариант

2,66 13

13

2,71 51

51

2,81 66

66

2,84 764

705

2,85 Нет генерации

2146

Указанная особенность генерации, по-видимому, связана с тем, что на длинах волн 2,66,
2,71 и 2,81 мкм она осуществляется при относительно малой населенности энергетических уровней перехода 4I11/2—4I13/2, при которой интенсивность апконверсионных переходов невелика. Так, при генерации на длине волны 2,81 мкм средняя скорость апконверсионных пере-

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 9

Модель многочастотной 3 мкм-генерации излучения лазеров

29

ходов с нижнего лазерного уровня (определяется как γ1N12 ) составляет 0,5⋅1022 см–3⋅с–1, тогда как при генерации на длине волны 2,85 мкм скорость данного процесса равна 3,3⋅1022 см–3⋅с–1.
При Qp < 12 кВт/см3 и τр < 3 мс в среде Er:YSGG излучение генерируется на одной дли-

не волны (2,79 мкм). При более высоких значениях Qp и τр как в Er:YLF-, так и в Er:YSGGлазере происходит изменение длины волны излучения до значений 2,83, 2,87 мкм.
Для Er:YAG-лазера в случае постоянного для всех длин волн перехода 4I11/2—4I13/2 зна-

чения потерь излучения в резонаторе δ(λ) генерация излучения происходит на длинах волн

2,69 и 2,94 мкм, что иллюстрируется представленным на рис. 2, а графиком зависимости

плотности потока фотонов от времени S(t). При этом относительно небольшие (в пределах 1 %)

изменения δ(λ) по спектру, например за счет изменения коэффициента отражения выходного

зеркала лазера, приводят к возникновению пиков генерации на других длинах волн внутри

перехода 4I11/2—4I13/2 (например, при λ=2,83 мкм — рис. 2, б). Отметим, что длина волны 2,94 мкм является наиболее сильной в спектре вынужденного излучения Er:YAG-лазера и

проявляется при различных сочетаниях параметров потерь в резонаторе. Данный факт хоро-

шо согласуется с известными экспериментальными данными.

а)

S⋅1018, см–2⋅с–1 20

λ=2,69 мкм λ=2,94 мкм

10

0 б) S⋅1018, см–2⋅с–1
20
10

37,5 75 112,5 150 187,5 225 262,5 t, мкс
λ=2,83 мкм λ=2,71 мкм λ=2,94 мкм

0 37,5 75 112,5 150 187,5 225 262,5 t, мкс Рис. 2
Заключение. Представленные теоретические исследования параметров свободной генерации излучения Er:YLF-, Er:YAG- и Er:YSGG-лазеров при одиночных импульсах накачки показали, что при селективной накачке Er:YLF- и Er:YSGG-лазеры имеют близкие КПД генерации, превышающие таковые для Er:YAG-лазера. Обеспечение эффективной генерации излучения Er:YAG-лазера требует использования излучения накачки с существенно большей импульсной мощностью. При постоянных по спектру потерях излучения в резонаторе лазера генерация излучения лазеров на эрбиевых кристаллах происходит на нескольких длинах волн.
Статья подготовлена по результатам работы, выполненной при поддержке Российского фонда фундаментальных исследований, грант №12-02-31203.
СПИСОК ЛИТЕРАТУРЫ
1. Иночкин М. В., Назаров В. В., Сачков Д. Ю. и др. Динамика спектра генерации трехмикронного Er:YLFлазера при полупроводниковой накачке // Оптич. журн. 2009. Т. 76, № 11. С. 62—67.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 9

30 М. В. Иночкин, В. В. Назаров, Д. Ю. Сачков и др.

2. Labb´e C., Doualan J.-L., Girard S. et al. Absolute excited state absorption cross section measurements in Er3+:LiYF4 for laser applications around 2.8 µm and 551 nm // J. Physics: Condensed Materials. 2000. Vol. 12. P. 6943—6957.

3. Koetke D., Huber G. Infrared excited-state absorption and stimulated-emission cross sections of Er3+-doped crystals // Applied Physics B. 1995. Vol. 61. P. 151—158.

4. Smirnov V. A., Shcherbakov I. A. Rare-Earth scandium chromium garnets as active media for solid-state lasers // IEEE J. of Quantum Electronics. 1988. Vol. 24, N 6. P. 949—959.

5. Payne S. A., Chase L. L., Smith L. K. et al. Infrared cross-section measurements for crystals doped with Er3+, Tm3+, and Ho3+ // IEEE J. of Quantum Electronics. 1992. Vol. 28, N 11. P. 2619—2630.

6. Sardar D. K., Bradley W. M., Perez J. J. et al. Judd—Ofelt analysis of the Er3+(4f11) absorption intensities in Er3+doped garnets // J. of Appl. Phys. 2003. Vol. 93, N 5. P. 2602—2607.

7. Gruber J. B., Quagliano J. R., Reid M. F. et al. Energy levels and correlation crystal-field effects in Er3+-doped garnets// Physical Rev. B. 1993. Vol. 48, N 21. P. 15561—15573.

8. Ткачук А. М., Разумова И. К., Мирзаева А. А. и др. Up-конверсия и заселение возбужденных уровней иона эрбия в кристаллах LiY1-xErxF4 (x=0,003-1) при непрерывной накачке излучением InGaAs-лазерных диодов // Оптика и спектроскопия. 2002. Т. 92, № 1. С. 73—88.

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

Сведения об авторах — канд. физ.-мат. наук; Санкт-Петербургский национальный исследова-
тельский университет информационных технологий, механики и оптики, кафедра лазерной техники и биомедицинской оптики; E-mail: m_inochkin@mail.ru — канд. техн. наук; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра лазерной техники и биомедицинской оптики; E-mail: lab255@grv.ifmo.ru — канд. техн. наук, доцент; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра лазерной техники и биомедицинской оптики; E-mail: dsachkov@gmail.com — канд. техн. наук; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра лазерной техники и биомедицинской оптики; E-mail: l_khloponin@yahoo.com — д-р техн. наук, профессор; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра лазерной техники и биомедицинской оптики; заведующий кафедрой; E-mail: khramov@grv.ifmo.ru

Рекомендована кафедрой лазерной техники и биомедицинской оптики

Поступила в редакцию 26.04.13 г.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 9