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

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

А.А. Дроздов, С.А. Козлов

УДК 535.1

ФАЗОВАЯ САМОМОДУЛЯЦИЯ ОДНОПЕРИОДНЫХ ОПТИЧЕСКИХ ВОЛН
А.А. Дроздов, С.А. Козлов

Получены аналитические и численные решения уравнения динамики поля светового импульса в изотропной диэлектрической среде с нормальной групповой дисперсией и безынерционной кубической нелинейностью для граничного условия вида однопериодной оптической волны. Показано, что уширение спектра таких предельно коротких по числу колебаний импульсов в нелинейной среде происходит как в высоко-, так и в низкочастотную область. Максимум спектральной плотности смещается в низкочастотную область. Возникновения привычного в среде с кубичной нелинейностью излучения утроенных частот не происходит, а генерируется волна с центральной частотой в примерно четыре с половиной раза большей исходной центральной частоты волны. Ключевые слова: однопериодные импульсы, фазовая самомодуляция, нелинейные диэлектрические среды.

Введение

Фазовая самомодуляция интенсивных сверхкоротких световых импульсов в прозрачных оптических средах теоретически изучалась ранее в большом числе работ [1, 2]. В последнее десятилетие особо

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 2 (72)

99

ФАЗОВАЯ САМОМОДУЛЯЦИЯ ОДНОПЕРИОДНЫХ ОПТИЧЕСКИХ ВОЛН
большой интерес вызывали исследования особенностей этого нелинейного явления для импульсов, содержащих лишь несколько колебаний светового поля [3–5]. В работе [6], по-видимому, впервые явление фазовой самомодуляции экспериментально наблюдалось для предельно коротких по числу колебаний однопериодных импульсов (single cycle pulses), которые генерируются в терагерцовом спектральном диапазоне [7]. В настоящей работе особенности фазовой самомодуляции однопериодных волн в изотропных диэлектрических средах с нормальной групповой дисперсией и безынерционной кубической нелинейностью рассмотрены теоретически.

Уравнение динамики поля излучения из малого числа колебаний в нелинейной диэлектрической среде

Уравнение, описывающее эволюцию линейно поляризованного излучения, спектр которого нахо-

дится в области прозрачности нелинейной диэлектрической среды и в той ее части, где групповая дис-

персия нормальная, можно записать в виде [8]

E z



N0 c

E t

a

3E t3



gE 2

E t



0

,

(1)

где Е – электрическое поле светового импульса; z – пространственная координата, вдоль которой распро-

страняется импульс; t – время; c – скорость света в вакууме; N0, a – параметры, характеризующие зависимость линейного показателя преломления n0 среды от частоты  излучения:

n0 ()  N0  ac2 ,

(2)

параметр g характеризует безынерционную кубическую по полю нелинейность поляризационного откли-

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

g  2n2 c [9].

Для получения аналитических решений уравнения (1) целесообразно провести следующую его нормировку. Введем новые безразмерные переменные E  E E0 , z  z λ 0 , t  t T0 , где E0 – макси-

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

период колебаний поля входного излучения, λ0  cT0 N0 – центральная длина волны. В этих перемен-

ных уравнение (1) принимает следующий вид:

E z



E t



A

3E t3

G

E 3 t



0,

(3)

где безразмерный коэффициент A  ac N0T02 характеризует дисперсию линейного показателя прелом-

ления диэлектрической среды; G  cgE02 3N0 описывает вклад безынерционной кубической по полю

нелинейности. При дальнейшем анализе уравнения (3) знак «» будем опускать.

Аналитическое решение нелинейного уравнения динамики поля

Для получения аналитического решения уравнения (3) будем использовать метод последователь-

ных приближений Пикара [10], в котором малым будем предполагать параметр G. В этом разделе работы

ограничимся анализом решений нормированного уравнения (3) для случая, когда дисперсией линейного

показателя преломления среды можно пренебречь, т.е. когда выполняется условие G  A . Например,

это неравенство выполняется для кристалла ниобата лития ( N0  4, 7 , a  3, 41038 с3/см,

n2  5, 4 1012 см2/Вт [6]) в поле терагерцового излучения с интенсивностью I  1, 0 108 Вт/см2 и цен-

тральным периодом T0  4, 7 1012 с. Этим условиям, при которых в работе [6] экспериментально на-

блюдалась фазовая самомодуляция, соответствуют G  7, 7 105 , A  9,8106 . Приведенные параметры

дисперсии ниобата лития N0 и a были рассчитаны авторами методом наименьших квадратов при ап-

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

MgO:LiNbO3 [11]. В соответствии с выбранным методом решение уравнения (3) будем искать в виде

E  E0  GE1  G2E2  ... ,

(4)

в котором ограничимся только двумя первыми слагаемыми. Тогда в представлении (4) нормированное уравнение (3) будет иметь следующий вид:

100

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 2 (72)

А.А. Дроздов, С.А. Козлов

 E 0  z



E 0  t



0

 

E

1

 z



E 1 t



E 0 3 t

 0.

Решение первого уравнения системы (5) имеет вид [12]:

(5)

E0 z,t  E0 t  z .

(6)

Второе уравнение системы (5) также несложно решить в квадратурах, переписав его в новых пе-

ременных z  z , τ  t  z . Тогда его решение примет следующий вид [13]:

E1



z,

τ







z

z0

 τ



E



0



τ



3

dz

.

(7)

Временную эволюцию поля световой волны на входе в среду (при z0  0 ) будем рассматривать в

виде однопериодного импульса

E0 (0, )



E0

 p

 exp 




2 2p

  

,

(8)

где E0 – его амплитуда,  p – длительность. Волна (8) хорошо аппроксимирует, например, терагерцовое

излучение фотопроводящих поверхностей полупроводников при их облучении фемтосекундными им-

пульсами [14]. Соответственно спектр такого излучения имеет вид

G  




E0 0,τ eiωτdτ  



 2



E0

2pie

 

 p 2

2

 

.

(9)

Из (9) следует, что максимум спектра приходится на центральную частоту 0  2  p .

Решение (7) с учетом (8) может быть приведено к виду

E1



z,





3E0

  

2 3p

  1   

2

2 2p

  

   

exp

  



32 2p

  



 

z



z0



,

(10)

спектр которого имеет вид:

G1  

3 36

E032p


2e

1 3

p 2

2

 

 1  

1 2

  

 p 3

2  

   

.

(11)

Из решения (11) следует, что спектральная плотность излучения, генерируемого в нелинейной

среде, обращается в ноль на утроенной центральной частоте   3 2  p  30 . Этот вывод чрезвычайно

любопытен, так как генерация излучения на утроенных частотах – это классическое нелинейное оптиче-

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

пульса (8) для центральной частоты спектра исходного излучения это не так. Оценим теперь частоту,

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

тотах. Взяв производную от выражения (11) по частоте и приравняв ее к нулю, находим, что

m1ax  42  p  4, 60 .

(12)

Таким образом, итерационное решение (4) уравнения динамики светового поля (1) для однопери-

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

E



z,





E 0 ()



GE1 ( z,

)



E0

 p

 exp 




2 2p

  

 1 

 3G 


 2p

  1   

2

2 2p

   

 exp 




22 2p

  



 





z



z0


 .


(13)

Иллюстрации аналитического решения нелинейного уравнения динамики поля для граничного условия вида однопериодной волны

На рис. 1 приведена иллюстрация решения (13) при G(z  z0 )  0,3 . На рис. 1, а, изображена рас-
считанная в первой итерации поправка (10) к электрическому полю волны, а на рис. 1, б, – поправка (11) к модулю спектра. Пунктиром показаны временная и спектральная структуры поля однопериодной волны на входе в среду. На рис. 1, в, приведено общее решение (13), которое описывает результат самовоздействия однопериодной оптической волны в нелинейной среде без учета ее дисперсии.

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 2 (72)

101

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

Рис. 1. Изменение временной и спектральной структур однопериодной оптической волны из-за ее
самовоздействия в нелинейной диэлектрической среде: итерационная добавка к полю E1 (а) и модулю
спектра G1 (б); результирующее поле с совместным действием эффектов фазовой самомодуляции
и генерации утроенных частот (в) и с учетом только эффекта фазовой самомодуляции (г). Пунктиром показаны структуры поля и спектра начального импульса
Рис. 1, б, иллюстрирует обращение в ноль на утроенной центральной частоте модуля спектра G1 поправки (10) к полю волны. Видно, что генерируемая в нелинейной среде добавка к основным частотам излучения простирается от трех до семи центральных частот. Ее максимум находится на частоте в примерно 4,5 раза большей исходной центральной частоты волны. Из рис. 1, в, видно, что максимум и минимум однопериодной волны из-за самовоздействия излучения в нелинейной среде начинают запаздывать во времени. Под самовоздействием здесь будем понимать эффекты фазовой самомодуляции и генерации в поле начального излучения утроенных частот, которые в бездисперсионной среде не разделимы. Видно, что при совместном влиянии этих явлений происходит и искажение временного профиля импульса. Для сравнения на рис. 2, г, приведено изменение временной структуры волны, которое происходило бы без учета эффекта генерации утроенных частот. Оно, естественно, полностью соответствует хорошо известному характеру фазовой самомодуляции для импульсов из большого числа колебаний.

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

Для численного моделирования динамики поля излучения в нелинейной оптической среде удобна другая нормировка уравнения (1). Введем новые переменные E  E E0 , z  a30z , t  0t , где E0 –
максимальное значение напряженности электрического поля на входе в нелинейную среду, 0 – цен-
тральная частота излучения, и перейдем к «запаздывающему» времени   t  N0 z c . В новых переменных уравнение (1) приводится к виду [3, 8]

E z



3E 3

GE2

E 

0,

(14)

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

G  cg E02 3N0  4n2 I ac02 . Например, рассмотренному выше кристаллу ниобата лития в поле терагер-

цового излучения с центральной частотой 1,0 ТГц и интенсивностью I  1, 0 108 Вт/см2 соответствует

G  5, 4 102 . В качестве модели однопериодного импульса на входе среду (при z  0 ) будем по-

прежнему брать зависимость (8), которая в новых переменных принимает вид

E0 (0,

t

)



t 0 p

 exp 




t2 022p

 

,



(15)

На рис. 2 приведены результаты численного моделирования по уравнению (14) эволюции времен-

ной структуры поля однопериодной волны (15) в среде с параметром G  5, 4 102 .

102

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 2 (72)

А.А. Дроздов, С.А. Козлов
Рис. 2. Динамика электрического поля Е, нормированного на максимальное значение входного поля Е0,
в среде с параметром G  5, 4 102 Из рис. 2 видно, что уже на расстоянии z  1,0 ( z  1,2 мм) происходит увеличение длительности импульса, измеряемой по уровню половины максимальной амплитуды поля, за счет дисперсии линейного показателя преломления среды примерно в 4 раза. Расчеты показывают, что уширения спектра импульса при таком параметре G практически не происходит, т.е. среда для излучения является линейной. На рис. 3, 4 приведены аналогичные результаты численного моделирования по уравнению (14) эволюции временной и спектральной структур исходно однопериодной волны (15), распространяющейся в среде с параметром G  5, 4 . Из рис. 3 видно, что по мере генерации новых колебаний электрического поля на заднем фронте волны уже с расстояния z  0,4 ( z  0,5 мм) от основной волны отделяется импульс, максимальная амплитуда поля которого составляет величину порядка 2% амплитуды входного поля.

Рис. 3. Динамика электрического поля Е, нормированного на максимальное значение входного поля Е0, в среде с параметром G  5, 4 . Во вставке изображено электрическое поле образующейся новой времен-
ной структуры
На рис. 4 приведена иллюстрация эволюции спектральной плотности G однопериодной волны в среде с параметром G  5, 4 , нормированной на максимальную входную спектральную плотность G0. Из
рисунка видно, что при распространении исходно однопериодной волны в среде с таким параметром G происходит значительное смещение спектральной плотности как в высокочастотную область с характерным минимумом на частоте в примерно 3,5 раза большей начальной центральной частоты, так и в низкочастотную область с приблизительно 20% сдвигом максимума ее плотности. Видно, что происходит генерация волны с центральной частотой примерно в 4 раза большей центральной, спектральная плотность которой составляет величину порядка 3% от максимальной спектральной плотности.

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 2 (72)

103

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

Рис. 4. Динамика модуля спектральной плотности G, нормированной на максимальную входную
спектральную плотность G0, в среде с параметром G  5, 4 (а–е). Смещение максимума спектральной
плотности волны на выходе из среды в низкочастотную область (ж)
Заключение
В работе рассмотрены особенности фазовой самомодуляции световых волн, содержащих на входе в нелинейную среду лишь одно полное колебание светового поля. Показано, что уширение спектра однопериодных импульсов в изотропной диэлектрической среде с нормальной дисперсией и безынерционной кубической нелинейностью происходит как в высоко-, так и в низкочастотную область. Возникновения излучения утроенных частот не происходит, а генерируется волна с центральной частотой в примерно 4,5 раза большей исходной центральной частоты волны.
Работа поддержана грантами НШ-5707.2010.2 и РНП 2.1.1/4923.
Литература
1. Ахманов С.А., Выслоух В.А., Чиркин А.С. Оптика фемтосекундных лазерных импульсов. – М.: Наука, 1988. – 312 с.
2. Агравал Г. Нелинейная волоконная оптика. – М.: Мир, 1996. – 324 с. 3. Козлов С.А., Самарцев В.В. Основы фемтосекундной оптики. – М.: Физматлит, 2009. – 292 с. 4. Беспалов В.Г., Козлов С.А., Сутягин А.Н., Шполянский Ю.А. Сверхуширение спектра интенсивных
фемтосекундных лазерных импульсов и их временное сжатие до одного колебания светового поля // Оптический журнал. – 1998. – Т. 65. – № 10. – С. 85–88. 5. Белашенков Н.Р., Дроздов А.А., Козлов С.А., Шполянский Ю.А., Цыпкин А.Н. Фазовая самомодуляция фемтосекундных световых импульсов, спектры которых сверхуширены в диэлектриках с нормальной групповой дисперсией // Оптический журнал. – 2008. – Т. 75. – № 10. – С. 3–8. 6. Hebling J., Hoffmann M.C., Ka-Lo Yeh, Tóth G. and Keith A.N. Nonlinear lattice response observed through terahertz SPM // Springer Series in Chemical Physics. – 2009. – V. 92. – Part 8. – P. 651–653. 7. Крюков П.Г. Фемтосекундные импульсы. – М.: Физматлит, 2008. – 208 с. 8. Козлов С.А., Сазонов С.В. Нелинейное распространение импульсов длительностью в несколько колебаний светового поля в диэлектрических средах // ЖЭТФ. – 1997. – Т. 111. – № 2. – С. 404–418. 9. Беспалов В.Г., Козлов С.А., Шполянский Ю.А. Метод анализа динамики распространения фемтосекундных импульсов с континуумным спектром в прозрачных оптических средах // Оптический журнал. – 2000. – Т. 67. – № 4. – С. 5–11. 10. Корн Г., Корн Т. Справочник по математике. – М.: Наука, 1973. – 832 с. 11. Palfalvi L., Hebling J., Kuhl J., Peter A. and Polgar K. Temperature dependence of the absorption and refraction of Mg-doped congruent and stoichiometric LiNbO3 in the THz range // J. Appl. Phys. – 2005. – V. 97. – 123505. 12. Фарлоу С. Уравнения с частными производными для научных работников и инженеров. – М.: Мир, 1985. – 384 с.

104

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 2 (72)

А.В. Альфимов, Е.М. Арысланова, Д.Н. Вавулин и др.

13. Буяновская Е.М., Козлов С.А. Закономерности взаимодействия встречных световых волн из малого числа колебаний в нелинейных диэлектрических средах // Научно-технический вестник СПбГУ ИТМО. – 2010. – Т. 66. – № 2. – С. 23–29.
14. Беспалов В.Г., Крылов В.Н., Путилин С.Э., Стаселько Д.И. Генерация излучения в дальнем ИК диапазоне спектра при фемтосекундном оптическом возбуждении полупроводника InAs в магнитном поле // Оптика и спектроскопия. – 2002. – Т. 93. – № 1. – С. 158–162.

Дроздов Аркадий Анатоьевич Козлов Сергей Аркадьевич

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

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 2 (72)

105