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

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

О.А. Столповская, С.Э. Путилин, В.Г. Беспалов, С.А. Козлов

УДК 535.1
ИНТЕРФЕРЕНЦИЯ ИЗЛУЧЕНИЯ НА УТРОЕННЫХ И КОМБИНАЦИОННЫХ ЧАСТОТАХ ПРИ ОДНОВРЕМЕННОМ НЕЛИНЕЙНОМ ОТРАЖЕНИИ ОТ ДИЭЛЕКТРИКА ДВУХ ФЕМТОСЕКУНДНЫХ ИМПУЛЬСОВ НА ОСНОВНОЙ И УДВОЕННОЙ ЧАСТОТАХ
О.А. Столповская, С.Э. Путилин, В.Г. Беспалов, С.А. Козлов

Представлены теоретические и экспериментальные результаты одновременного нелинейного отражения от диэлектрика двух фемтосекундных импульсов на основной ω1 и удвоенной ω2 = 2ω1 частотах при падении света на границу раздела сред, близком к нормальному. Показано, что при отражении такого двухкомпонентного излучения происходит генерация импульсов на утроенной и комбинационных частотах. При этом спектр излучения, отраженного на утроенных частотах, представляет собой наложение двух компонент: утроенной 3ω1 и комбинационной 2ω2 – ω1. Интерференция этих двух компонент приводит к модуляции спектра излучения, отраженного на утроенных частотах, при изменении временной задержки между падающими импульсами. Ключевые слова: нелинейное отражение, утроенные и комбинационные частоты, фемтосекундный импульс.

Введение

Распространение фемтосекундного излучения в прозрачных твердых телах без оптического пробоя вещества оказывается возможным при значительно больших интенсивностях, чем для более длинных оптических импульсов [1]. Это позволяет при взаимодействии в нелинейных оптических средах двух или более фемтосекундных импульсов разных спектральных составов эффективно генерировать излучение на кратных и комбинационных частотах. Обзор работ по взаимодействию в нелинейных средах сонаправленных и встречных фемтосекундных импульсов можно найти в [2], а также, например, в [3, 4]. В настоящей работе рассмотрено взаимодействие двух разночастотных фемтосекундных импульсов при нелинейном отражении.
В работах [5–7] были получены соотношения, связывающие параметры поля отраженного излучения с характеристиками падающего на нелинейный диэлектрик света, спектр которого может занимать значительную часть диапазона прозрачности диэлектрической среды, для случая нормального падения плоской поперечно однородной волны. Особенностью этих работ было использование непосредственно полевых уравнений при выводе зависимостей. В работе [8] также на основе полевых уравнений были получены аналоги формул Френеля для случая отражения от нелинейного диэлектрика поперечно слабонеоднородной (параксиальной) волны со сверхшироким временным спектром. В этой работе была выведена зависимость пространственно-временного спектра светового импульса, отраженного от диэлектрической среды с нерезонансной дисперсией и нелинейностью, от параметров пространственно-временного спектра падающего параксиального пучка, временной спектр которого в общем случае рассматривается сверхшироким, вида

 Gref

(, kx , ky

)



n1 n1

 

n2 n2

 1 

c2 n1n2

k

2 x



k

2 y

2

  Ginc



3

gcS n1  n2

,

(1)

где Ginc и Gref – пространственно-временные спектры падающего и отраженного излучения соответственно; n1(ω) и n2(ω) – линейные показатели преломления граничащих сред, для которых зависимость от частоты имеет вид

ni ()  Ni  ni (),

ni

()



cai 2



c

bi 2

,

где Ni, ai, bi – константы, характеризующие нерезонансную дисперсию линейного показателя преломления среды; g  2n2 / c – описывает безынерционную нелинейность поляризационного отклика среды, n2

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



  S(, kx , ky ) 

E3 inc

(t,

x,

y)e dti(t kx xky  y)

dx

dy





      1 (2)6



Ginc

(



,

kx



k


x

,

ky



k


y

)

Ginc

(



,

k


x



kx,

k


y



k y) 



 Ginc (, kx, ky) ddkxdky ddkxdky,

где Еinc – напряженность электрического поля падающей волны. В работах [8, 9] были рассмотрены сценарии нелинейного отражения одиночных импульсов из

малого числа колебаний, а также экспериментально исследовано отражение высокоинтенсивного фемто-

секундного излучения от поверхности диэлектрика.

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

33

ИНТЕРФЕРЕНЦИЯ ИЗЛУЧЕНИЯ НА УТРОЕННЫХ И КОМБИНАЦИОННЫХ …

Теоретический анализ одновременного нелинейного отражения от диэлектрической среды двух фемтосекундных импульсов разных частот и генерации при этом излучения кратных и комбинационных частот в настоящей работе проведем на основе формулы (1).

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

Применим (1) для анализа нелинейного отражения нормированной суперпозиции двух фемтосекундных гауссовых импульсов:

  Einc (t, x,

y)



e 2

t 1

2 

e  2

x 1

2 

e  2

y 1

2 

 sin

1

t  1  t 2



   E  e  e  e sin2

 2

t

 t 2

2 

 2

x 2

2 

 2

y 2

2 

E1

2

t  t  2  t 2

,

(2)

где Е1 и Е2 – амплитуды первого и второго импульсов соответственно; τ1 и τ1 – их временные длительности; ρ1 и ρ2 – поперечные размеры; 1 = 2с / 1 и 2 = 2с / 2 – центральные частоты импульсов; 1 и 2

– их центральные длины волн; t – временная задержка второго импульса относительно первого;  и

 – коэффициенты фазовой модуляции.

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

модуляции (     0 ) с амплитудами Е1 = 2Е2, центральными длинами волн λ1 = 800 нм и λ2 = 400 нм,
пространственно-временными параметрами ρ1 = ρ2 = 10 ∙ λ1 и τ1 = τ2 = 26 фс и пиковой интенсивностью I = 5 ∙ 1013 Вт/см2, которые падают из воздуха (N1 = 1, a1 = b1 = 0) на кварцевое стекло (N2 = 1,4508, a2 = 2,7401 ∙ 10–44 с3/см, b2 = 3,9437 ∙ 1017 1/(с ∙ см), n2 = 2,9 ∙ 10–16 см2/Вт [10]) с временной задержкой
Δt = 20 фс. На рис. 1 представлены нормированные поле и спектр падающего излучения. На рис. 2 изображены нормированные поле и спектр отраженного излучения. Можно показать, что
дифракционно-дисперсионным изменением поля и спектра при указанных параметрах излучения и вещества можно пренебречь, и основными эффектами, как видно из рисунка, являются нелинейные, заключающиеся в генерации импульсов на кратных и комбинационных частотах 3ω1 (совпадающей с 2ω2 – ω1), 2ω1 + ω2, 2ω2 + ω1 и 3ω2. На утроенных частотах отраженного излучения наблюдается провал спектральной плотности, который объясняется наложением комбинационной частоты 2ω2 – ω1 на утроенную частоту 3ω1.
На рис. 3 проиллюстрирован волновой пакет, отраженный на утроенной частоте. Из рисунка видно, что этот импульс имеет меньшую длительность по сравнению с импульсами, отраженными на основ-

ных частотах ω1 и ω2 (в 3 раз). Введем в падающее излучение (2) фазовую самомодуляцию, которая может возникнуть в оптиче-
ском тракте лазерной системы. Положим  = 6 ∙ 10–4ω1,  = 2 ∙ 10–4ω1. Тогда излучение, отраженное на
утроенных частотах, будет иметь вид, который представлен на рис. 4. Из рисунка видно, что спектр третьей гармоники уширяется, провал в спектральной плотности увеличивается и смещается.
1 10–3 2

1,5
0 1

Einc, отн. ед. Ginc, отн. ед.

0,5 –1

–1/2 0 1/2 t, фс
а

0 12 3 45 6 /1
б

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

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

О.А. Столповская, С.Э. Путилин, В.Г. Беспалов, С.А. Козлов

1

10–4

10–6

43

Gref, отн. ед.

Eref, отн. ед.

32 0
21

03 4

56

1 /1

–1

–1/2 0 1/2 t, фс
а

0 12 3 4 56 /1 б

Рис. 2. Нормированные поле (а) и спектр (б) отраженного излучения, включающего импульсы на кратных и комбинационных частотах 3ω1 (совпадающей с 2ω2 – ω1), 2ω1 + ω2, 2ω2 + ω1 и 3ω2. На вставке изображен спектр кратных и комбинационных частот отраженного излучения (увеличенный)

10–3 1

10–6 4

2

Gref, отн. ед.

Eref, отн. ед.

0

0 2,5

3 /1

3,5

–1

–1/2

0 1/2 t, фс

Рис. 3. Нормированное поле импульса, отраженного на утроенной частоте 3ω1. На вставке показан его спектр

10–3 1

10–6 4

Gref, отн. ед.

2

0

0 2,5

3 3,5

/1

Eref, отн. ед.

–1

–1/2

0 t, фс

1/2

Рис. 4. Нормированное поле импульса, отраженного на утроенной частоте 3ω1, при фазовой модуляции падающего излучения. На вставке показан его спектр

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

35

ИНТЕРФЕРЕНЦИЯ ИЗЛУЧЕНИЯ НА УТРОЕННЫХ И КОМБИНАЦИОННЫХ …
Рассмотрим теперь излучение, отраженное на утроенных частотах, при изменении временной задержки между падающими импульсами. Пусть Δt изменяется от 15 фс до 25 фс с шагом 0,5 фс. На рис. 5 продемонстрировано изменение спектра третьей гармоники при увеличении временной задержки Δt. Из рисунка видно, что провал в спектральной плотности смещается в зависимости от временной задержки между импульсами. При этом изменяется ширина и амплитуда спектральной плотности.
410–6

G, отн. ед.

2,5 3 3,5

/1

Рис. 5. Спектр импульса, отраженного на утроенной частоте, при изменении задержки между падающими импульсами на 0,5 фс

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

В экспериментах нами была использована фемтосекундная лазерная система на кристаллах титансапфира с усилением чирпованного импульса c пиковой мощностью 30 ГВт [11]. Система имела следующие параметры: длительность импульса 40 фс, энергия одиночного импульса до 1,5 мДж, частота повторения импульсов 50 Гц.
Излучение от фемтосекундной системы, пройдя линзу L (рис. 6) с фокусным расстоянием 51 см, попадало на кристалл второй гармоники SHG, установленный на подвижку, с помощью которой можно было изменять положение кристалла вдоль оси распространения пучка. Далее излучение основной частоты и второй гармоники падало на переднюю грань образца кварцевого стекла С. Чтобы избежать влияния излучения, отраженного от задней грани, образец был наклонен относительно падающего луча (угол падения составлял 5,5°), а область взаимодействия была максимально приближена к матовой боковой грани образца. Фокальная плоскость линзы была расположена в 6 см за передней гранью кристалла. Плотность мощности была близка к энергии пробоя, но не превышала ее и составляла около 4 ТВт/см2.
Отраженное от образца излучение проходило через ультрафиолетовые стекла (УФС-5) F суммарной толщиной 10 мм. За фильтрами на линейный транслятор с шаговым двигателем MTS был установлен спектрограф ASP100 SD. Работа шагового двигателя линейного транслятора и спектрографа управлялись при помощи компьютера PC. На рис. 7 представлен зарегистрированный на спектрографе спектр отраженного излучения без внесения кристалла второй гармоники в схему. Из рисунка видно, что спектр третьей гармоники имеет плавный контур.

Рис. 6. Схема экспериментальной установки
36 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 3 (79)

О.А. Столповская, С.Э. Путилин, В.Г. Беспалов, С.А. Козлов

G, отн. ед. G, отн. ед.

1400 1200 1000 800 600 400

1500 0
1000
500
0 255

260 265 270 , нм

200

0 200

300 400

500 600 700 800 900 , нм

Рис. 7. Спектр отраженного излучения, прошедшего через УФС-фильтры. На вставке изображен спектр излучения, отраженного на утроенных частотах, имеющий плавный контур

G, отн. ед.

2000 1500 1000

500 0 250 260 270

, нм

Рис. 8. Спектр импульса, отраженного на утроенной частоте, при изменении положения кристалла второй гармоники с шагом 2 мм
Установим кристалл второй гармоники так, как показано на схеме (рис. 6), на расстоянии 7,5 см от линзы L, и с помощью подвижки будем изменять его положение по оси распространения пучка с шагом 2 мм в сторону линзы. Изменяя положение кристалла второй гармоники, мы можем управлять временной задержкой между импульсами основной частоты и второй гармоники вследствие дисперсии воздуха.
На рис. 8 проиллюстрировано изменение спектра третьей гармоники в зависимости от положения кристалла второй гармоники. Как видно из рисунка, в спектре отраженного на утроенной частоте излучения присутствуют провалы, вызванные наложением двух спектральных компонент – утроенной и комбинационной, что доказывает генерацию более высоких комбинационных частот в спектре отраженного излучения. Также из рисунка видно, что провалы в спектральной плотности излучения, отраженного на утроенной частоте, смещаются в зависимости от временной задержки между импульсами основной частоты и второй гармоники, что подтверждает вышеприведенные теоретические результаты.
Заключение
Рассмотренный случай одновременного нелинейного отражения от диэлектрика двух фемтосекундных импульсов на основной и удвоенной частотах расширяет теорию нелинейного отражения, полученную в работе [8], а также экспериментально демонстрирует возможность эффективной генерации утроенных и комбинационных частот. Интерференция между сгенерированными спектральными компонентами на утроенной частоте позволяет качественно изменять отраженное излучение. В частности, изменяя положение кристалла второй гармоники по оси распространения падающего пучка, можно управлять спектральной шириной импульса, отраженного на утроенной частоте.
Работа поддержана грантами Министерства науки и образования Российской Федерации 16.740.11.0459 и РФФИ 11-02-01346а.
Литература

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

37

ИНТЕРФЕРЕНЦИЯ ИЗЛУЧЕНИЯ НА УТРОЕННЫХ И КОМБИНАЦИОННЫХ …

1. Von der Linde D., Schuler H. Breakdown threshold and plasma formation in femtosecond laser-solid interaction // J. Opt. Soc. Am. B. – 1999. – V. 13. – № 1. – Р. 216–222.
2. Козлов С.А., Самарцев В.В. Основы фемтосекундной оптики. – М.: Физматлит, 2009. – 292 с. 3. Сазонов С.В., Устинов Н.В. Резонансно-параметрический механизм оптического выпрямления и ге-
нерации высших гармоник // Оптика и спектроскопия. – 2009. – Т. 106. – № 3. – С. 473–480. 4. Буяновская Е.М., Козлов С.А. Взаимодействие встречных световых импульсов из малого числа коле-
баний в нелинейных диэлектрических средах и генерация излучения на комбинационных частотах в этом процессе // Оптика и спектроскопия. – 2011. – Т. 111. – № 2. – С. 325–332.
5. Bespalov V.G., Kozlov S.A., Oukrainski A.O., Sazonov S.V., Shpolyanskiy Yu.A. Self-action of femtosecond pulses with continuum spectrum // Proc. of SPIE. – 1999. – V. 3735. – Р. 43–54.
6. Розанов Н.Н. Отражение сверхкоротких импульсов от границы среды Друде–Лоренца // Оптика и спектроскопия. – 2003. – Т. 94. – № 3. – С. 449–452.
7. Ястребова Н.В., Шполянский Ю.А., Козлов С.А. Нелинейное отражение импульсов из малого числа колебаний светового поля от просветленной границы раздела сред // Оптический журнал. – 2004. – Т. 71. – № 6. – С. 78–83.
8. Мохнатова О.А., Козлов С.А. Нелинейное отражение фемтосекундного лазерного суперконтинуума // ЖЭТФ. – 2008. – Т. 13. – Вып. 2. – С. 260–270.
9. Мохнатова О.А., Путилин С.Э., Беспалов В.Г., Козлов С.А. Генерация третьей гармоники при отражении от кристалла фемтосекундного излучения высокой мощности // Научно-технический вестник СПбГУ ИТМО. – 2010. – № 4 (68). – С. 11–15.
10. Berkovsky A.N., Kozlov S.A., Shpolyansky Yu.A. Self-focusing of few-cycle light pulses in dielectric media // Phys. Rev. A 72. – 2005. – V. 72. – P. 043821.
11. Беспалов В.Г., Киселев В.М., Кисляков И.М. и др. Антистоксов самосдвиг и уширение спектра излучения фемтосекундного лазера в сильно поглощающей среде // Оптика и спектроскопия. – 2009. – Т. 106. – № 4. – С. 670–679.

Столповская Ольга Александровна Путилин Сергей Эдуардович Беспалов Виктор Георгиевич Козлов Сергей Аркадьевич

– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, аспирант,
o.stolpovskaya@gmail.com – Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, кандидат физ.-мат. наук, доцент, seputilin@yandex.ru – Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, доктор физ.-мат. наук, профессор, victorbespaloff@gmail.com – Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, доктор физ.-мат. наук, профессор, декан, kozlov@mail.ifmo.ru

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