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

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

ГЕНЕРАЦИЯ ИЗЛУЧЕНИЯ КРАТНЫХ ЧАСТОТ ПРИ ВОЗБУЖДЕНИИ …
УДК 535.3, 535:530.182
ГЕНЕРАЦИЯ ИЗЛУЧЕНИЯ КРАТНЫХ ЧАСТОТ ПРИ ВОЗБУЖДЕНИИ В ГАЗЕ ПЛАЗМЫ СИЛЬНЫМ ПОЛЕМ ИМПУЛЬСОВ ИЗ МАЛОГО ЧИСЛА
КОЛЕБАНИЙ
С.А. Штумпф, А.А. Королев, С.А. Козлов
Рассмотрена генерация высокочастотного излучения при оптическом пробое неона фемтосекундными импульсами из малого числа колебаний светового поля. Показано, что при уменьшении числа колебаний в импульсе, вплоть до всего одного, происходит сильное смещение кратных гармоник в коротковолновую область спектра. Плазменная нелинейность диэлектриков приводит к усилению перераспределения энергии излучения в коротковолновую область. Ключевые слова: предельно короткие импульсы, плазменная нелинейность, генерация кратных гармоник.
Введение
Развитие нелинейной оптики импульсов из малого числа колебаний светового поля, экспериментально получаемых сегодня уже во многих научных лабораториях [1, 2], ставит новые задачи по созданию теории взаимодействия таких предельно коротких импульсов (ПКИ) с веществом. Одной из таких задач является построение математической модели распространения в диэлектрических средах настолько сильных полей ПКИ, что оптические электроны вещества становятся квазисвободными, но разрушение среды из-за скоротечности взаимодействия и малой энергии ПКИ не происходит.
К настоящему времени разработан ряд моделей генерации плазмы в диэлектрических средах в поле фемтосекундного излучения, содержащего большое число колебаний оптического поля. Эти модели, основанные на рассмотрении динамики огибающих световых волн, позволили изучить процессы зарождения филаментов, локализации плазменного канала и т.п. в поле длинных импульсов (например, [3, 4]). Однако применение метода медленно меняющейся огибающей импульса для ПКИ, например, содержащего всего одно полное колебание поля, становится неплодотворным, поскольку понятие огибающей для ПКИ теряет свое физическое содержание (например, [5–7] и обзоры в них).
В работе [8] авторами было выведено уравнение динамики сильного поля ПКИ в диэлектрической среде, которое учитывает электронную нелинейность, обусловленную изменением населенностей высоковозбужденных энергетических состояний и движением электронов в квазисвободном состоянии. В работе [9] на основе анализа решений этого уравнения изучены условия доминирования различных физических факторов при самовоздействии ПКИ и приведены иллюстрации диапазонов интенсивности, длительности и спектрального состава излучения, когда плазменная нелинейность среды становится преобладающей. Показано, что модель корректно описывает возрастание коэффициента нелинейного показателя преломления оптической среды в фиолетовой части видимого спектра излучения и ближнем УФ диапазоне. В работах [10, 11] с помощью этой модели объяснен ряд особенностей генерации низкочастотного терагерцового излучения при оптическом пробое воздуха двуцветным фемтосекундным излучением, в частности, квазипериодическая зависимость длины плазменного филамента от исходной временной расстройки между фемтосекундными импульсами на основной и удвоенной частотах.
В настоящей работе методами численного моделирования выведенного в [8, 9] полевого уравнения определены особенности генерации излучения высоких частот при оптическом пробое газа сильными полями световых импульсов из малого числа, вплоть до всего одного, колебаний поля.
Уравнение динамики сильного поля импульса из малого числа колебаний в диэлектрической среде
Уравнение динамики поля светового импульса высокой интенсивности из малого числа колебаний в диэлектрической среде было выведено в [8] на основе формализма матрицы плотности в приближении трехзонной энергетической модели среды [12]. При описании плазменной нелинейности третье энергетическое состояние рассматривалось как зона квазисвободного движения электронов [10, 11]. Результатом вывода явилась система волнового и динамических материальных уравнений, учитывающая как инерционную часть электронной кубической нелинейности, так и генерируемую в сильном поле плазменную нелинейность. В приближении однонаправленного распространения светового импульса [5, 6] (при пренебрежении самоотражаемым излучением) эта система может быть представлена в виде
28 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 3 (79)

С.А. Штумпф, А.А. Королев, С.А. Козлов

 E  z



a

3E 3



g

E 3 



 

 

g (1) E

 

E 

2 



g(2) E2

2E 2

 



2 cn0

j0

 

 





 р

 E2

 



j

j

 E3,

 c

(1)

где E – электрическое поле излучения; z – направление его распространения;   t  n0 z – время в соc
провождающей системе координат; c – скорость света в вакууме; n0 и a характеризуют зависимость линейного показателя преломления n диэлектрической среды от частоты  вида

n()  n0  ca2 ,

(2)

которую несложно при необходимости уточнить слагаемыми более высокого порядка по  [4]; коэффи-

циенты g , g(1) , g(2) описывают дисперсию коэффициента нелинейного показателя преломления диэлек-

трика [9] вида

n2 ()  n20  A2 ,

где

n20

 3 cg ; 2

A  1 c(g(1)  3g(2) ) ; коэффициенты 2

  n02 1 , 4T21 21

  e2 me* (

p223 32 )2

характеризуют инерци-

онную плазменную нелинейность диэлектрика [10], pij и ij – дипольные моменты и частоты, соответ-

ствующие переходам i  j в трехзонной модели; T21 ,  p – времена поперечной релаксации и релакса-

ции населенностей в паре состояний (1), (2), me* – эффективная масса электрона в зоне проводимости

диэлектрика; c – среднее время столкновительной релаксации свободных электронов. Первое уравнение системы (1) описывает самовоздействие поля излучения в среде в условиях ее

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

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

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

жение этих электронов под действием электромагнитных волн.

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

На рис. 1, 2 приводятся результаты численного моделирования изменения спектра ПКИ высокой интенсивности вида
( t )2
E  E0e tи sin(0t) , где E0 – амплитуда; tи – длительность; 0 – центральная частота излучения, при его распространении в неоне. При моделировании использованы значения линейного показателя преломления n0  1,000067
[11] и коэффициента нелинейного показателя преломления неона n2  0,0061022 м2 Вт [12].

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

29

ГЕНЕРАЦИЯ ИЗЛУЧЕНИЯ КРАТНЫХ ЧАСТОТ ПРИ ВОЗБУЖДЕНИИ …

G/G0 1
10–1
10–2

E/E0 1
0
–1 –10 0 10 20 . фс

10–3

10–4 0 2 4 6 8 10 12 14 16 18 /0

Рис. 1. Спектр излучения в нелинейной среде (неоне) после распространения на расстояние z = 5 см (пунктирная линия – при учете в расчете только безынерционной кубичной нелинейности среды,
сплошная линия – при дополнительном учете инерционности электронной нелинейности).
Вставка на графике показывает общий вид импульса на входе в среду

G/G0 1

E/E0 1

0 10–1
–1 10–2 –5 0 5 . фс

10–3 0

1

2

G/G0 1

345

67 8 а

9 10 /0

E/E0

1

10–1 0
–1 10–2 –5 0 5 . фс

10–3 0 1 2 3 4 5 6 7 8 9 10 /0 б
Рис. 2. Спектр излучения в нелинейной среде (неоне) после распространения на расстояние z = 5 см импульса из трех полных колебаний поля (пунктирная линия – спектр импульса на входе в среду) (а); одного колебания светового поля (пунктирная линия – спектр импульса при тех же условиях без учета плазменной нелинейности) (б). Вставка на графике показывает общий вид поля входных импульсов
30 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 3 (79)

С.А. Штумпф, А.А. Королев, С.А. Козлов
На рис. 2 показано, что высокоинтенсивный импульс всего из одного полного колебания светового поля при нелинейном распространении в диэлектрической среде генерирует второй максимум спектра излучения не в области третьей гармоники основной частоты излучения (как для более длинных импульсов), а в области четвертой (что отмечалось также в работе [14]). Аналогично смещаются и прочие компоненты спектра, соответствующие кратным гармоникам более высоких порядков (рис. 2, б). Видно, что процесс носит характер, близкий к пороговому; генерация кратных гармоник для импульса из двух полных колебаний светового поля (спектр приведен на рис. 2, а) носит, напротив, вполне традиционный для изучаемого явления характер (генерируются гармоники нечетных порядков, картина аналогична представленной на рис. 1). Влияние ионизационных компонент модели (1) проявляется в усилении перераспределения спектра излучения в высокочастотную область, что также можно видеть из иллюстрации рис. 2, б: амплитуда излучения в области 6-й гармоники увеличивается при учете плазменных компонент приблизительно в 2 раза, а в области 9-й гармоники – уже примерно в 10 раз.
Рис. 1 демонстрирует сверхуширение спектра импульса длительностью tи  10фс с центральной
длиной волны излучения 0  2с 0  780нм при его входной интенсивности I  11015 Вт/см2 , на расстоянии распространения в неоне z  5см (генерация гармоник в этой среде, вплоть до 103-й, описана,
например, в [13]). Из рисунка видно, что сверхуширение спектра излучения реализуется в виде генерации гармоник нечетных порядков. При учете инерционности электронной нелинейности в указанных параметрах импульса происходит ослабление амплитуды каждой последующей гармоники по отношению к предыдущей в 2–5 раз. Учет только безынерционной части кубичной по полю нелинейности позволил бы реализовать формирование «плато» с 9-й по 17-ю гармонику практически без относительного убывания их интенсивности.
Заключение
В работе продемонстрирован эффект сильного смещения частоты генерируемых кратных гармоник в высокочастотную область при распространении в неоне высокоинтенсивного импульса из всего одного полного колебания светового поля (например, второй максимум спектра излучения наблюдается в области четвертой гармоники, а не третьей). Показано усиление перераспределения энергии в высокочастотную область при генерации квазисвободных электронов в диэлектрической среде сильным полем светового излучения.
Работа поддержана грантом Министерства образования и науки Российской Федерации 16.740.11.0459 и РФФИ 11-02-01346а.
Литература
1. Brabec Th., Krausz F. Intense few-cycle laser fields: Frontiers of nonlinear optics // Rev. Mod. Phys. – 2000. – V. 72. – № 2. – P. 545–591.
2. Paul P.M., Toma E.S. et al. Observation of a train of attosecond pulses from high harmonic generation // Science. – 2001. – V. 292. – Р. 1689–1692.
3. Панов Н.А., Косарева О.Г., Кандидов В.П., Акозбек Н., Скалора М., Чин С.Л. Локализация плазменного канала при множественной филаментации в воздухе // Квантовая электроника. – 2007. – Т. 37. – № 12. – С. 1153–1158.
4. Kosareva O.G., Panov N.A., Urypina D.S., Kirilova M.V., Mazhorova A.V., Savel’ev A.B., Volkov R.V.,
Kandidov V.P., Chin S.L. Optimization of a femtosecond pulse self-compression region along a filament in air // Applied Physics B. – 2008. – V. 91. – P. 35–43. 5. Козлов С.А., Сазонов С.В. Нелинейное распространение импульсов длительностью в несколько колебаний светового поля в диэлектрических средах // ЖЭТФ. – 1997. – Т. 111. – В. 2. – С. 404–418.
6. Bespalov V.G., Kozlov S.A., Shpolyansky Yu.A., Walmsley I.A. Simplified field wave equations for nonlinear propagation of extremely short light pulses // Phys. Rev. A. – 2002. – V. 66. – № 1. – Р. 013811–013821.
7. Шполянский Ю.А. Комплексная огибающая произвольного оптического сигнала для описания взаимодействия импульсов с континуумными спектрами и различными центральными частотами в прозрачных нелинейных средах // Научно-технический вестник СПбГУ ИТМО. – 2009. – № 3. – C. 26–34.
8. Штумпф С.А., Королев А.А., Козлов С.А. О природе и инерционности нелинейной поляризованности диэлектриков в поле световых импульсов из малого числа колебаний // Известия РАН. Серия физическая. – 2006. – Т. 70. – № 1. – С. 118–125.
9. Штумпф С.А., Королев А.А., Козлов С.А. Распространение сильных полей световых импульсов из малого числа колебаний в диэлектрических средах // Известия РАН. Серия физическая. – 2007. – Т. 71. – № 2. – С. 158–161.

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

31

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

10. Андреев А.А., Беспалов В.Г., Козлов С.А., Штумпф С.А. и др. Генерация сверхширокополосного терагерцового излучения при оптическом пробое воздуха двумя разночастотными фемтосекундными импульсами // Оптика и спектроскопия. – 2009. – Т. 107. – № 4. – С. 569–576.
11. Leonard P.J. Refractive Indices, Verdet Constants, and Polarizabilities of the Inert Gases // Atomic Data and Nuclear Data Tables. – 1974. – V. 14. – P. 21–37.
12. Shaw M.J., Hooker C.J. and D.C. Wilson. Measurement of the non- linear refractive index of air and other gases at 248 nm // Opt. Commun. – 1993. – V. 103. – P. 153–160.
13. Miyazaki K., Takada H. High order harmonic generation in the tunneling regime // Phys. Rev. A. – 1995. – V. 52. – P. 3007–3021.
14. Kozlov S.A., Drozdov A.A., Suchorukhov A.A., Kivshar Y.S. On peculiarities of self-action of one-period light waves // Book of Abstracts, Int. Workshop Nonlinear Photonics. – 24–26 Aug. – 2011. – P. 37.

Штумпф Святослав Алексеевич – Санкт-Петербургский национальный исследовательский университет

информационных технологий, механики и оптики, кандидат физ.-мат.

наук, доцент, stumpf@rain.ifmo.ru

Королев Александр Александрович – Санкт-Петербургский национальный исследовательский университет

информационных технологий, механики и оптики, кандидат физ.-мат.

наук, доцент, korolev@phd.ifmo.ru

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

– Санкт-Петербургский национальный исследовательский университет

информационных технологий, механики и оптики, доктор физ.-мат. наук,

профессор, декан, kozlov@mail.ifmo.ru

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