Генератор излучения на основе нелинейного метаматериала с отрицательным показателем преломления
Автор: Коробко Дмитрий Александрович, Остаточников Владимир Александрович, Золотовский Игорь Олегович, Абрамов Алексей Сергеевич, Минвалиев Рамиль Наильвич
Журнал: Известия Самарского научного центра Российской академии наук @izvestiya-ssc
Рубрика: Физика и электроника
Статья в выпуске: 4-5 т.15, 2013 года.
Бесплатный доступ
В работе рассмотрено трехволновое параметрическое взаимодействие в планарном волноводе из нелинейного метаматериала, обладающего в определенном частотном диапазоне свойствами “левой” среды. Получены условия параметрической генерации излучения в квадратично нелинейном метаматериале, в том числе и с учетом влияния керровской нелинейности. Предложенная схема параметрического генератора может быть использована для построения модели микроволнового или ТГц лазерного источника на основе нелинейного метаматериала.
Генератор излучения, параметрическое взаимодействие, метаматериал, левая среда
Короткий адрес: https://sciup.org/148202436
IDR: 148202436
Текст научной статьи Генератор излучения на основе нелинейного метаматериала с отрицательным показателем преломления
В последние годы проводится активное исследование оптических свойств метаматериалов – искусственных композитных материалов, которые имеют особые электромагнитные свойства, не встречающиеся в природе. Одним из наиболее интересных типов метаматериалов являются структуры, у которых в определенной частотной области действительная часть показателя преломления имеет отрицательные значения. Это происходит только в том случае, если в этой области отрицательными одновременно являются действительные части диэлектрической и магнитной проницаемостей ε и μ . Отличительным свойством таких материалов является то, что вектор Пойнтинга распространяющейся в нем электромагнитной волны противоположен по направлению волновому вектору и образует с векторами электрического и магнитного полей левую тройку, в силу чего такие среды получили название «левых» [1-3]. В настоящее время известны метаматериалы с отрицательными значениями показателя преломления в микроволновой, инфракрасной, и даже в оптической областях спектра [4-6].
Из-за уникального свойства волн, распространяющихся в «левой» среде в частотном диапа-
зоне с показателем преломления n ( to ) < 0 , переносить энергию в направлении противоположном своему волновому вектору, их принято называть «обратными» волнами. Многие интересные свойства «левых» сред проявляются при наличии связи между «обратной» и обычной «прямой» волной, переносящей энергию по направлению своего волнового вектора. Реализация такой связи возможна, например, в туннельно-связанных волноводах изготовленных из различных материалов [7]. В простейшем случае оптические свойства «левой» среды рассматривают в линейном приближении, и связанные волны имеют одну частоту. Показано, что при этом распространение волн описывается теми же уравнениями, что и в лазере с распределенной обратной связью [8, 9] и при выполнении соответствующих условий в подобной структуре также можно возбудить генерацию излучения [10].
Одним из новых направлений оптики является исследование нелинейных свойств «левых» сред. При этом следует отметить, что в материале подобного типа, представляющем сложную композитную структуру, оптическая нелинейность определяется комбинацией нелинейностей отдельных компонент и может принимать крайне высокие значения, например, при включении элементов с высоконелинейными характеристиками [11, 12] или вблизи плазмонного резонанса [13]. Учет нелинейных факторов, в частности, позволяет рассмотреть в подобном метаматериале параметрические процессы и исследовать связь волн, распространяющихся на разных частотах, при этом одна из волн может являться обратной, т.е. ее частота лежит в диапазоне n(to) < 0. В работах [14-16] рассмотрены модели «левой» среды с квадратичной нелинейностью и описаны генерация второй гармоники и параметрическое усиление при трехволновом взаимодействии в такой среде. Трехволновое параметрическое усиление и генерация второй гармоники в метаматериале, обладающем в некотором диапазоне показателем преломления n < 0 , наблюдались экспериментально [17, 18]. В работе [14], так же как и в линейном случае, показана возможность резонансной связи между прямой и обратной волнами, что дает возможность компенсировать потери в метаматериале и, в принципе, осуществить параметрическую генерацию. В работах [16, 19] рассмотрены также четырехволновой параметрический процесс в метаматериале с кубической нелинейностью и генерация третьей гармоники. В отличие от этих результатов в настоящей работе мы обращаем внимание на то, что нелинейный отклик метаматериала имеет комплексный характер, и при рассмотрении параметрических процессов в квадратично-нелинейной среде обязательно необходимо учитывать и кубически-нелинейный (керровский) отклик метаматериала. В работе получены условия параметрической генерации излучения в квадратично-нелинейном метаматериале, обладающими свойствами «левой» среды, причем в отличие от указанных выше работ, рассмотрен случай не только высокочастотной накачки, но и накачки с частотой меньшей чем частота «обратной» сигнальной волны. Обсуждается влияние на квадратичную параметрическую генерацию керровской нелинейности метаматериала, показано, что при корректном подборе параметров среды и волны накачки керровская нелинейность облегчает выполнение условий параметрической генерации в квадратично-нелинейном метаматериале.
ОСНОВНЫЕ УРАВНЕНИЯ
В настоящей работе будет рассмотрена модель планарного волновода из нелинейного композитного материала, обладающего значительным коэффициентом квадратичной нелинейности. В этом случае при параметрическом взаимодействии входных волн – сигнальной волны с частотой ωs и интенсивной волны накачки на частоте ωp будет происходить генерация «холостой» волны с частотой toq = top — tos . Рассмотрим простую геометрию взаимодействия, соответствующую сонаправленным волновым векторам. Наиболее эффективным усиление сигнальной и холостой волн будет при выполнении условия фазового синхронизма k (to ) = k (toq ) + k (to ) , (1) здесь k (to ) = ton (to )| С — коллинеарные волновые вектора взаимодействующих волн. Вследствие дисперсии показателя преломления n(to)
условие фазового синхронизма точно практически никогда не выполняется. При этом величина отстройки от фазового синхронизма "k = k ( to p J — k ( tos ) — k ( to q ) может быть достаточно большой, что является часто существенным препятствием для достижения эффективного параметрического преобразования частот.
Параметрическое взаимодействие волн будем рассматривать в монохроматическом приближении, которое вполне оправдано при описании импульсов излучения с длительностью т > 10 - 8 c , при этом уравнения трехволнового взаимодействия в среде данного типа могут быть записаны следующим образом [20, 21]:
' = — bA- * e " k — a A ,
∂z spq ss
∂A
-q = — zb q - p - * e z " kz — a A , ∂ z
∂A
"A = —zbPAAqe— ap-p. ∂z
Здесь Ai – комплексные амплитуды соответствующих монохроматических волн, bi и αi - коэффициенты квадратичной нелинейности и показатели поглощения волновода на этих частотах.
Далее перейдем к описанию параметрического взаимодействия в волноводе из материала, в определенном диапазоне частот проявляющем свойства «левой» среды. Предполагаем, что в волновод подается волна накачки с частотой ωp , лежащей в области n (to) > 0 , при этом одна из взаимодействующих волн (для определенности положим, что это сигнальная волна) попадает в диапазон, характеризующийся отрицательным показателем преломления, переносит, таким образом, энергию в навстречу волне накачки и становится «обратной» волной. Холостая волна, для которой действительная часть показателя преломления n (to) > 0 , распространяется в прямом направлении. В дальнейшем необходимо рассмотреть два варианта, имеющих принципиальное отличие. В первом из них интенсивное излучение накачки имеет максимальную частоту, и параметрический процесс при этом характеризуется распадом кванта излучения с частотой накачки ωp на два фотона с частотами tos и toq: top = tos + toq. В частотную область с n(to) < 0 попадает волна с частотой ωs . Этот случай назовем случаем высокочастотной накачки. Во втором случае интенсивная накачка на частоте ωp приводит к усилению сигнальной компоненты с частотой более высокой, чем частота накачки tos > top , причем эта частота tos лежит в диапазоне n(to) < 0. Таким образом, сигнальная волна вновь является обратной, и параметрический процесс сопровождается излучением фотона с частотой toq = tos — top. Отстройка от фазового синхронизма тогда может быть записана как Ak = k (tos) — k (top) — k (toq). Это случай низкочастотной накачки. Его уравнения могут быть получены из системы (2) перестановкой индексов s О p . Схемы параметрическо-
Aq (0) = A,, As (L) = 0. (4)
Запишем решения системы (3 а) для высокочастотной накачки с учетом этих граничных условий:
го взаимодействия волн для двух случаев приведены на рис. 1. Воспользуемся тем, что интенсивности сигнальной и холостой волн малы
| a s —a q + iA k 1 2d cos d ( L — z ) + iA sin d ( L — z )
A,(z) = A, expl------;------z I-------------------------------- q I 2 J 2d cos dL+zAsindL
a —a — iAk A * ( z ) = A ,expl -----
2 io * sin d ( L—z ) 2 d cos dL + iA sin dL
(5а)
по сравнению с интенсивностью волны накачки, и запишем систему (2) в приближении неистощимой накачки Ap = const :
—
^ = ia* Ae—iAkz az s q
A* — as As ,
aa
-= = — ia A* eAkz — a A . dz q s q q
(3а)
d A s
az
= — ia Ae-iAkz — aA,
sq ss
5A q az
=— i^ AA^—a A.
qs qq
(3б)
Здесь учтено то, что в выбранной нами геометрии прямая волна переносит энергию по направлению оси z , а обратная – в противоположном направлении. Полученные системы уравнений для связанных волн описывают параметрический процесс в волноводе, проявляющем в определенном диапазоне свойства «левой» среды, в случае высокочастотной (3 а) и низкочастотной (3 б) накачки. Введенные коэффициенты связи волн определяются через коэффициенты квадратичной нелинейности волновода и амплитуду накачки a s = bsAp , a q = bqAp . Для волновода длины L рассмотрим следующие граничные условия, соответствующие вводу прямой волны в волновод слева и отсутствию обратной волны на противоположном конце волновода:
В случае низкочастотной накачки решение для обратной волны записывается для амплитуды, а не для ее сопряжения:
As ( z ) = A n exp

— iAk
2io * sin d ( L — z ) 2 d cos dL + id sin dL
, (5б)
при этом введены следующие обозначения: условно называемая комплексной отстройкой величина 5 = A k — i ( a s + aq ) и параметр «периодичности» d = ( 5 74 ± g 2 ) 1/2, g 2 = a * a q = b s b q|Ap |2 . Знак «+» соответствует высокочастотной накачке, «-» – низкочастотной. В данных условиях коэффициенты отражения и пропускания волновода, т.е. отношения мощностей выходящей обратной и прошедшей прямой волны к входной мощности прямой волны, выглядят как
R = PM P q (0)
T P q ^ LL
P q (0)
2 a s sin dL
2 d cos dL + i5 sin dL
2 d
2 d cos dL + iA sin dL
2 (6)
exp ( ( a s — a q ) L ) .
Графики коэффициентов R и T для случаев высоко- и низкочастотной накачки приведены на рис. 2 и 3.
В случае высокочастотной накачки (рис. 2) можно отметить периодическую зависимость коэффициентов R и T от длины волновода L . В отличие от прохождения излучения в «обычной»

Рис. 1. Схемы трехволнового параметрического усиления в «левом» метаматериале для высокочастотной (а) и низкочастотной (б) накачки: k i – волновые вектора; S i – вектора Пойнтинга взаимодействующих волн;
A k - отстройка от фазового синхронизма

Рис. 2. Коэффициент отражения (а) и пропускания (б) планарного волновода длины L в случае высокочастотной накачки при параметрах g =100м - 1 , a =5м \ a s =0.5м1 , отстройкам от фазового синхронизма А к = (0,25,50) м соответствуют сплошная, штриховая и пунктирная линии
поглощающей среде, характеризующимся экспоненциальным падением интенсивности, здесь при определенных длинах волновода можно видеть рост коэффициентов пропускания и отражения. Такое «просветление» волновода связано с наличием резонансной связи между прямой и обратной волнами. При выполнении условия фазового синхронизма А к = 0 коэффициенты R и T могут обращаться в бесконечность, что говорит о принципиальной возможности параметрической генерации излучения сразу на двух частотах ωs и ωq при накачке волновода на частоте ωp . В отсутствии поглощения a s = a q = 0 условие обращения в ноль знаменателей в выражениях (6) и, соответственно, бесконечно больших коэффициентов отражения и пропускания имеет простой вид L = (2 n + 1) п /2 g . При наличии поглощения период по L увеличивается. Варьируя мощность накачки, можно изменять значения коэффициента связи g прямой и обратной волн и добиться совпадения длины волновода с характерной длиной возникновения генерации. Таким образом, рассматриваемая модель параметрического генератора аналогична лазеру с распределенной обратной связью [22].

Рис. 3. Коэффициент отражения (а) и пропускания (б) планарного волновода длины L в случае низкочастотной накачки при параметре g = 100 м - 1 и параметрах a s = a q = А к = 0 - сплошная линия, А к = 250 м " 1, a s = 0.5 м-1, a q = -0.5 м-1 -штриховая линия,
А к = 250 м " 1, a s = 0.5 м " 1, a q = -10 м " 1 -пунктир
При низкочастотной накачке (рис. 3) также можно отметить перекачку энергии из прямой в обратную волну, однако в этом случае фотон накачки имеет частоту top = tos — toq , и накачка не способна усилить излучение одновременно на двух частотах ωs и ωq . Однако, из рис. 3 можно видеть, что при наличии дополнительного усиления на одной из частот и достаточной отстройке А к > 2 g + as + aq, также возможно получить сингулярные значения коэффициентов отражения и пропускания рассматриваемого волновода. Таким образом, здесь также можно говорить о параметрической генерации излучения. Важно отметить то, что в отличие от предыдущего случая при низкочастотной накачке не требуется выполнения условия фазового синхронизма А к = 0, наоборот отстройка в этом случае должна иметь определенное достаточно высокое значение, при котором действительная часть пара- метра «периодичности» d = (52/4 — g2) пре
-
взойдет его мнимую часть. Так как условие фазового синхронизма практически трудно реализуемо, это является преимуществом способа низкочастотной накачки. Необходимые для возникновения генерации значения отстроек и отвечающие им необходимые значения усиления могут быть приближенно вычислены. Для этого рассмотрим условие возникновения генерации излучения в рассматриваемом волноводе (справедливое как для низко- так и для высокочастотной накачки), которое вытекает из условия равенства нулю знаменателей в (6):
2d tg dL = i —.
Решения этого уравнения в общем случае могут быть получены только численными методами. Получим его приближенные решения в случае низкочастотной накачки и отстроек, достаточно высоких по сравнению с коэффициентом параметрической генерации и суммарным коэф- фициентом усиления Ak > 2g, |az | = |aq + as . Так как в данном случае генерация без дополнительного усиления невозможна, параметр суммарного поглощения полагается отрицательным a^ < 0. Приближенные выражения для действительной и мнимой частей коэффициента «периодичности» d могут быть записаны как

d i = -Aka
^ Ak 2 -a ^
I 4
A - 1/2
Уравнение (7) преобразуем к виду exp( idL) + exp( - idL) = ± —. (8)
g
Пользуясь тем, что мнимая часть di существенно больше от нуля, в левой части этих выражений можно пренебречь вкладом одного из слагаемых по сравнению с другим. В данном случае при усилении a^ < 0 пренебрегаем слагаемым ^ exp(idL) . Приравнивая модули и фазы в левых и правых частях равенств (8), получаем условия генерации diL = ln
( 4 A k 2 + a | / g )
d r L = nm + arctg( a s/A k ),
m = 0, ± 1, ± 2...
Ограничиваясь случаем сильной связи g >> | a z| и пренебрегая слагаемыми более чем первого порядка малости по ( g A k ) - можно дать оценку для необходимой отстройки A k и отвечающего ей усиления | a z |, обеспечивающих генерацию излучения в волноводе длины L :
Ak « 2л](лт /L)2 + g2, m = 1,2,...,
I a z|”
In ( A k / g )
2 L
Из полученных выражений видно, что при высоких коэффициентах связи между прямой и обратной волной, возможно получать параметрическую генерацию волн с существенной отстройкой от фазового синхронизма. Изменяя мощность накачки, можно менять «критическое» значение отстройки и в определенном диапазоне осуществлять перестройку генерации по частотам to s и toq . Так как обычно область частот, связанная с отрицательными значениями показателя преломления метаматериала весьма невелика, перестройке в основном будет подвержена частота прямой волны toq . Важным следствием полученных соотношений является также то, что при генерации волн с большими отстройками от синхронизма необходимо прикладывать существенное усиление. Реализовать дополнительное усиление в волноводе возможно разными способами. Некоторые из них могут быть основаны на использовании нелинейных свойств композитного материала, образующего волновод. В частности при достаточно высоком значении коэффициента кубической (керровской) нелинейности это усиление может осуществляться за счет модуляционной неустойчивости мощной волны накачки или за счет ее вынужденного комбинационного рассеяния (ВКР) [23]. Подробнее это будет рассмотрено ниже.
ПАРАМЕТРИЧЕСКАЯ ГЕНЕРАЦИЯ В МЕТАМАТЕРИАЛЕ С КВАДРАТИЧНОКУБИЧНОЙ НЕЛИНЕЙНОСТЬЮ
Как уже отмечалось метаматериалы, обладающие свойствами «левой» среды, часто характеризуются сложным нелинейным откликом. Это связано с тем, что эти материалы представляют собой гетероструктуры с пространственными неоднородностями размером менее длины волны, в качестве примера можно привести двумерную композитную структуру с квадратной решеткой периодических ячеек из микроскопических проводящих проволочек и разомкнутых кольцевых резонаторов [2, 3, 17, 24]. В работе [24] экспериментально продемонстрирована одновременная генерация второй и третьей гармоники в метаматериале на основе разомкнутых кольцевых резонаторов, при этом квадратичные и кубичные (керровские) нелинейные коэффициенты преобразования имели сопоставимые значения. Это показывает, что при описании нелинейных свойств подобной среды необходимо учитывать, по меньшей мере, два порядка нелинейности – квадратичную и кубичную. Другим примером могут быть металлические частицы малых размеров в диэлектрической матрице [13]. Подобные композиционные материалы проявляют уни- кальные нелинейно-оптические свойства, обусловленные, в том числе и тем, что их нелинейный отклик определяется, как минимум, двумя факторами: объемной диэлектрической проницаемостью матрицы и проницаемостью включенной в нее решетки металлических включений. Если в качестве диэлектрической матрицы выбран материал с сильной квадратичной нелинейностью, то вследствие того, что металлические микрочастицы обладают керровским нелинейным откликом, общий нелинейный отклик композита будет иметь сложный квадратично-кубичный характер, причем его величина будет значительно выше, чем в стандартном диэлектрике (например, кварцевом стекле). На сегодняшний день получены нелинейные композитные материалы с диэлектрической восприимчивостью третьего порядка X(3) около 10-8 ед СГСЭ [25], что более чем на пять порядков превышает значения керровской нелинейности для кварцевых световодов и соответствует параметру нелинейности у(3) ~ 102 Вт-1м■’ . Таким образом, из вышеизложенного следует вывод о том, что при рассмотрении квадратично-нелинейных эффектов в композитных материалах должна учитываться и кубическая составляющая нелинейного отклика.
При рассмотрении трехволнового взаимодействия в метаматериале, обладающем свойствами «левой» среды, наиболее важными эффектами, связанными с проявлениями керровской нелинейной поляризации являются кросс-модуляционное взаимодействие волн, вынужденное комбинационное рассеяние (ВКР) и модуляционная неустойчивость мощной волны накачки [23]. Влияние фазовой самомодуляции накачки на параметрическое взаимодействие сказывается лишь опосредованно через усиление волн, попадающих в диапазон ее модуляционной неустойчивости. Фазовой самомодуляцией сигнальной и холостой волн, вследствие их малой мощности, можно пренебречь. Учет кубичнонелинейных членов сохраняет общий вид уравнений для связанных волн (3 а, б), но при этом показатели поглощения сигнальной и холостой волн заменяются эффективными параметрами, зависящими от мощности накачки:
аs = а Пр —2iY(3)Pp, аq = aq -ПР —2iY(3)Pp.(11) Здесь через у(3) обозначен коэффициент керровской нелинейности волновода, а последние слагаемые в правой части соотношений (11) описывают кросс-модуляционное воздействие между накачкой с одной стороны и сигнальной или холостой волной с другой [23]. Через nsPp, nqPp обозначены коэффициенты усиления сигнальной и холостой волн, пропорциональные мощности накачки. Как уже отмечалось выше, их происхождение может быть связано с двумя видами проявления кубично-нелинейного отклика среды. Если частота накачки ωp находится в области аномальной дисперсии волновода, то монохроматическая волна малой амплитуды с частотой ω, находящейся в полосе |to- top < 2(у(3)Pp]|в21) , ( в2 — дисперсия групповых скоростей волновода) испытывает усиление, обусловленное модуляционной неустойчивостью волны накачки. На частотах to = top ±(2у(3)Pp /|в21) коэффициент модуляционного усиления достигает максимального значения
П (mod) p p = у (3) p p
ВКР накачки, в свою очередь, представляет собой еще один эффект, определяемый нелинейностью керровского типа. ВКР не зависит от дисперсии среды и связано с задержкой нелинейного отклика, приводящей к перекачке энергии накачки в область более низких частот (в стоксову компоненту) [20, 23]. В спектральном представлении ВКР может быть представлено как параметрический процесс рассеяния волны на колебаниях частиц среды – фононах, который приводит к появлению сопутствующих стоксовой и антистоксовой компонент с частотами to p — Ato и to p + Ato ( A to — частота фонона). Интенсивность антистокосовой компоненты значительно ниже, чем у стоксовой, поэтому в большинстве случаев ей можно пренебречь. Величина частотного сдвига Δ ω определяется параметрами среды. Процесс ВКР можно охарактеризовать коэффициентом рамановского усиления на частоте ω , пропорционального мощности накачки – П ( Ram ) ( to ) P p . На рис. 4 схематично приведены возможные варианты частотной зависимости коэффициента усиления, вызываемого ВКР и модуляционной неустойчивостью накачки. В большинстве случаев на частоте пробной монохроматической волны превалирует усиление, вызванное одним из механизмов, связанных с кубично-нелинейным откликом среды. Таким образом, волна накачки в рассматриваемой задаче может одновременно играть две роли: обеспечивать квадратично-нелинейный параметрический процесс ( to p = tos ± toq ) и усиливать, по крайней мере, одну из двух других участвующих в этом процессе волн за счет кубично-нелинейного эффекта (ВКР или модуляционной неустойчивости). Усиление сигнальной или холостой монохроматических волн, обусловленное этим эффектом, будем описывать при помощи эффективных коэффициентов ηsPp и ηqPp . Частоты параметри-

Рис. 4. Схематическое изображение частотной зависимости усиления, вызываемого ВКР (сплошные линии) и модуляционной неустойчивостью (штриховые линии) накачки: (а), (б) – частота накачки ωp находится в области аномальной дисперсии волновода; (в), (г) – нормальной дисперсии;
P p 1 - мощность накачки, соответствующая рис. (а) и (в), P p 2 отвечает рис. (б) и (г), P p 1 < P p 2 .
пропускания и отражения для случаев высоко- и низкочастотной накачки. На рис. 5 представлены линии уровня коэффициента отражения волновода в случае высокочастотной накачки при различных уровнях ее мощности. Также как и в рассмотренном выше случае строго квадратичной нелинейности (рис. 2) можно отметить точки с координатами (Lt, Aк t), в которых значение коэффициента стремится к бесконечности. Физически это означает, что в таком волноводе длины Li при выполнении условия на отстройку Δk ≈ Δki из шумовых компонент (например, спонтанной люминесценции) на частотах ωq и ωs может быть возбуждена параметрическая генерация. Следует отметить, что для высокочастотной накачки учет керровской нелинейности позволяет исключить требование полного фазового синхронизма из условий генерации и тем самым облегчить их выполнение. Также из рис. 5 можно видеть, что ВКР-усиление снижает критическую длину волновода Li и существенно расширяет (по Δk ) области высоких значений коэффициентов отражения и пропускания волновода, что дополнительно способствует чески взаимодействующих волн можно подобрать так, что как минимум один из коэффициентов усиления будет значительным, что в свою очередь существенно облегчит достижение параметрической генерации. Очень важным является также учет фазовой кросс-модуляции. Ее действие приводит к снижению эффективной отстройки от фазового синхронизма Ak = А к - 2/(3) Pp и при корректном подборе параметров позволяет добиваться генерации в значительно более широком диапазоне отстроек. Важным следствием действия кросс-модуляции является дополнительная возможность регулировки «критической» отстройки от синхронизма, и, соответственно, возможность подстройки частоты генерируемых параметрически волн.
Рассмотрим численный пример, соответствующий модели планарного волновода из метаматериала с размером в несколько см, в который вводится излучение накачки мощностью несколько Вт. Нелинейные свойства волновода описываются параметрами квадратичной b s = b q = gjP^ = 100 м -1 Вт " 1/2 и керровской Y (3) = 5 м-1Вт-1 нелинейности. Положим, что накачка обеспечивает полосу ВКР-усиления, в которую попадает прямая волна с частотой ωq , при этом ее коэффициент усиления составляет ηqPp . На частотах ωq и ωs коэффициент поглощения метаматериала равен αq и αs соответственно. Рассмотрим резонансные характеристики волновода, которые представлены коэффициентами

Рис. 5. Линии уровня коэффициента отражения в случае высокочастотной накачки для волновода с параметрами b s = b q = 100 м-1Вт-1/2, Y (3) = 5 м-1Вт-1, n q = 5 м -1 Вт’’, n s = 0 , a q = a s = 0.5 м-1 . P p = 1 Вт (а), Pp = 3 Вт (б)
выполнению условий параметрической генерации. Обращаем внимание на то, что область частот, в которой действительная часть показателя преломления метаматериала отрицательна, не очень широка: в микроволновом диапазоне при to s ~ 5 ГГц, его ширина не более нескольких сотен МГц, при ω s ~ ГГц не более 5 ГГц, в ТГц диапазоне при ω s ~ 1 ТГц не более 200 ГГц и т.п [6, 26]. При этом следует помнить о том, что в силу резонанс-ности процесса прохождения в этой полосе присутствует зона максимума поглощения, которую приходится исключить. Таким образом, в небольшом диапазоне частот отрицательный показатель n ( to ) изменяется довольно значительно (от минимально возможного значения до нуля), что позволяет получить необходимое значение отстройки. Изменяя частоту накачки ωp , можно осуществлять перестройку частот параметрического генератора, при этом обратная сигнальная волна остается в узком диапазоне частот с n ( to ) < 0 , а частота прямой холостой волны ωq может изменяться довольно широко. Пределами ее изменений могут быть частотные границы области высоких значений квадратичной нелинейности метаматериала. При постоянной частоте накачки тонкую регулировку параметрически генерируемых частот можно осуществлять при помощи изменения мощности накачки.
На рис. 6 показаны линии уровня коэффициента отражения волновода в случае низкочастотной накачки при различных ее мощностях. Так же как и в предыдущем случае можно отметить существование точек, в которых значение коэффициента стремится к бесконечности, т.е. может возбуждаться параметрическая генерация. При изменении мощности накачки изменяется значение критической отстройки от синхронизма и коэффициент ВКР-усиления, что в итоге приводит к изменению частот параметрической генерации. В отличие от предыдущего в данном случае при низкочастотной накачке возрастает роль ВКР-усиления. Если при высокочастотной накачке оно играло вспомогательную роль, то при низкочастотной для обеспечения генерации усиление обязательно, причем значение требуемого усиления увеличивается с ростом необходимой отстройки Δk (см. (10)). Таким образом, важно, чтобы частота усиливаемой волны совпадала с максимумом ВКР-усиления накачки. В результате получаемый параметрический генератор может быть охарактеризован как ВКР-лазер, обеспечивающий генерацию обратной волны с частотой to = to + to , где to — частота сигна-spq q ла ВКР рассеяния в рассматриваемом волноводе от накачки частоты ωp . Так как частота обратной волны может изменяться в небольших пределах в области n(to) < 0, возможность сильной

Рис. 6. Линии уровня коэффициента отражения в случае низкочастотной накачки для волновода с параметрами b s = b q = 100 м -1 Вт - 1/2, Y (3) = 5 м "’Вт "’ , n q = 5 м "’ Вт "’ , n s = 0 , a q = a s = 0.5 м "’ . P p = 1 Вт (а), P p = 3 Вт (б)
перестройки частоты генерации в этом случае ограничена. Тонкую подстройку частот также как и в предыдущем случае можно осуществлять, изменяя мощность накачки.
ЗАКЛЮЧЕНИЕ
В работе рассмотрено трехволновое параметрическое взаимодействие монохроматических волн в планарном волноводе из метаматериала, обладающего в определенном частотном диапазоне свойствами «левой» среды и характеризующегося высокой квадратичной нелинейностью. В приближении «неистощимой» накачки выделены два случая параметрического усиления, соответствующие высоко- и низкочастотной накачке, при этом одна из усиливаемых волн является «обратной». Для обоих случаев получены выражения для коэффициентов отражения и пропускания рассматриваемого волновода (6). Показано, что данная структура аналогична лазеру с распределенной обратной связью и получены условия параметрической генерации излучения.
Метаматериалы, обладающие свойствами «левой» среды представляют собой композитную структуру, их нелинейный отклик имеет сложный комплексный вид, при этом, если метаматериал обладает высокой квадратичной нелинейностью, его кубическая (керровская) нелинейность также должна быть учтена, что подтверждается и экспериментально [24]. Описанный выше процесс параметрического усиления прямой и обратной волн в метаматериале в настоящей работе рассмотрен с учетом влияния основных эффектов, основанных на керровской нелинейности – фазовой кросс-модуляции волн, модуляционной неустойчивости и вынужденного комбинационного рассеяния мощной волны накачки. Показано, что при соответствии параметров накачки и волновода керровская нелинейность метаматериала облегчает выполнение условий генерации. Для высокочастотной накачки это соответствие заключается в том, чтобы частоты трех волн (обратной, прямой и накачки) находились в частотном диапазоне высокой нелинейности метаматериала. В этом случае при достаточном уровне накачки в волноводе определенной длины будут выполнены условия параметрической генерации. Изменяя частоту и мощность накачки, можно управлять частотами генерации. В случае низкочастотной накачки в первую очередь необходимо обеспечить требуемый уровень усиления. Для этого частота прямой волны должна быть совмещена с частотой сигнала ВКР накачки. При достаточно высоком уровне накачки, обеспечивающем необходимое усиление, и подборе частот, при котором частота обратной волны to s = to p + to q ( ωq – частота сигнала ВКР рассеяния, ωp – частота накачки) также выполняются условия параметрической генерации. Следует отметить, что кросс-модуляционный керровский эффект в обоих случаях снижает эффективную отстройку от синхронизма, т.е. способствует генерации при более высоких отстройках Δ k и может быть использован для тонкой регулировки генерируемых частот. Предложенная схема параметрического генератора может быть использована для построения модели лазера (в первую очередь микроволнового и ТГц диапазона) на основе нелинейного метаматериала.
Работа выполнена при поддержке Министерства образования и науки России.
Список литературы Генератор излучения на основе нелинейного метаматериала с отрицательным показателем преломления
- Веселаго В.Г.//УФН, 1967. Т. 92. Вып. 3. С. 517 -526.
- Smith D.R., Padilla W.J., Vier D.C. et al.//Phys. Rev. Lett. 2000. Vol. 84. P. 4184 -4187.
- Shelby R.A., Smith D.R., Schultz S.//Science. 2001. Т. 292. ¹ 5514. P. 77 -79.
- Shalaev V.M., Cai W., Chettiar U.K. et al.//Optics Letters. 2005. Vol. 30. ¹ 24. P. 3356 -3358
- Вендик И.Б., Вендик О.Г.//ЖТФ, 2013, T. 83. Вып. 1. С. 3 -28.
- Cai W., Shalaev V.M.//Optical Metamaterials: Fundamentals and Applications, New York: Springer, 2010. 212 p.
- Барыкина Е.И., Золотовский И.О., Семенцов Д.И.//Опт. и спектр. 2012. Т. 112. Вып. 1. С. 117 -124.
- Alu A., Engheta N. Negative-refraction Metamaterials: Fundamental Principle and Applications/eds. Elefteriades G.V., Baltman K.G. New York: Wiley. 2005. 440 p.
- Маймистов А.И., Казанцева В.Е.//Опт. и cпектр. 2012. Т. 112. Вып. 2. С. 291 -298.
- Золотовский И.О., Коробко Д.А., Семенцов Д.И., Фотиади А.А.//Опт. и спектр. 2013. Т. 114. Вып. 2. С. 331 -337.
- Zharov A.A., Shadrivov I.V, and. Kivshar Y.S. //Nonlinear properties of left-handed metamaterials // Phys. Rev. Lett., Vol. 91, 037401, 2003.
- Lapine M.M., Gorkunov, and Ringhofer K.H. Nonlinearity of a metamaterial arising from diode insertions into resonant conductive elements//Phys. Rev. E, Vol. 67, 065601, 2003.
- Ахмеджанов И.М., Тищенко А.В., Щербаков А.А. Моделирование рассеяния света на наночастицах сложной формы методом обобщенных источников//Опт. и спектр. 2008. Т. 105. № 6. С. 1033 -1038.
- Popov A.K. and Shalaev V.M. Negative-index metamaterials: Second-harmonic generation, Manley-Rowe relations and parametric amplification//Appl. Phys. B, 2006. Vol. 84. P. 131 -137.
- Shadrivov I.V., Zharov A.A., and Kivshar Y.S. Second-harmonic generation in nonlinear left-handed metamaterials.//J. Opt. Soc. Am. B, 2006. Vol. 23. P. 529 -534.
- Maimistov A.I., Gabitov I.R., Eur. Phys. J. Special Topics, 2007. Vol. 147. P. 265 -286.
- Klein M.W., Enkrich C., Wegener M., Linden S.//Science. 2006. Vol. 313. P. 502 -504.
- Kozyrev A.V., Kim H., van der Weide W.R.//Appl. Phys. Lett., 2006. Vol. 88. 264101-3
- Popov A.K., Myslivets S.A., Shalaev V.M.//J. Opt. A: Pure Appl. Opt., 2009. Vol. 11. 114028-41.
- Ахманов С.А., Выслоух В.А., Чиркин А.С. Оптика фемтосекундных лазерных импульсов. М.: Наука. 1988.
- Сухоруков А.П. Нелинейные волновые взаимодействия в оптике и радиофизике. М.: Наука. 1988.
- Ярив А. Введение в оптическую электронику. M.: Мир: 1983. 398 с.
- Агравал Г. Нелинейная волоконная оптика. М.: Мир: 1996. 323 с.
- Klein M.W., Wegener M., Feth N., and Linden S. Experiments on second-and third-harmonic generation from magnetic metamaterials//Opt. Express, 2007. Vol. 15. P 5238 -5247.
- Olivares J., Requejo-Isidro J. del Coso R. de Nalda R, Solis J., Afonso C.N., Stepanov A.L., Hole, Townsend P.D., Naudon A. Large enhancement of the third order optical susceptibility in Cu-silica composites produced by low energy high current ion implantation//J.Appl. Phys. 2001. Vol. 90. P. 1064.
- Padilla W.J., Basov D.N., Smith D.R., Materials Today, 2006. Vol. 9. P. 28 -35.