Прикладные задачи
^^^^^^^^^^»нелинейной теории колебаний и вслн
УДК 621.372; 537.86; 537.87
РАСПРОСТРАНЕНИЕ ВОЛН ВДОЛЬ ДИФФУЗНОЙ ГРАНИЦЫ НЕЛИНЕЙНОГО МЕТАМАТЕРИАЛА
А. Б. Маненков
Институт физических проблем им. П.Л. Капицы РАН
Исследованы характеристики поверхностных ТМ мод, направляемых диффузной границей раздела двух нелинейных сред (метаматериала и обычного диэлектрика). Рассчитаны распределения полей и зависимости постоянных распространения от параметров переходного слоя. Изучены эффекты, возникающие из-за нелинейностей сред.
Ключевые слова: Поверхностные электромагнитные моды, нелинейные метаматериалы, диффузные границы раздела сред.
Введение
Поверхностные электромагнитные моды (ПЭМ)1 достаточно давно являются объектом многочисленных исследований [1,2]. Эти моды представляют собой мощный инструмент для различных физических экспериментов, например, при изучении свойств поверхностей материалов на СВЧ и в оптическом диапазоне [3-5]. В последние годы интерес к ПЭМ существенно возрос в связи с предложениями использовать их в устройствах интегральной оптики, в системах ближнепольной микроскопии, при изучении различных нелинейных эффектов, а также свойств искусственных сред (метаматериалов).
В данной работе рассмотрены характеристики ПЭМ, которые могут распространяться вдоль плоской нерезкой границы, которая разделяет среду из метаматериала и обычную среду. Исследован общий случай, когда среды могут быть нелинейными. К настоящему времени опубликовано несколько работ, посвященных изучению свойств ПЭМ, распространяющихся вдоль структуры со скачком параметров на границе раздела сред (см., например, [6,7]). Представляет интерес рассмотрение альтернативной модели, когда граница раздела является «размытой». В такой геометрии проницаемости непрерывно изменяются в тонком слое, который расположен между двумя однородными полубесконечными средами (структура с диффузной границей).
1Такие моды называют также поверхностными поляритонами, а в случае, когда одна из сред является металлом, - поверхностными плазмон-поляритонами.
Важно отметить, что в такой системе проницаемости могут быть близки к нулю, и это свойство может сильно изменять структуру полей. Кроме чисто теоретического интереса, изучение таких структур представляет также практический интерес, в первую очередь потому, что многие реальные структуры могут быть ближе именно к этой модели. Отметим, что, строго говоря, подобную систему следует рассматривать как слоистый волновод. Однако в случае, когда толщина переходного слоя намного меньше длины волны, в первом приближении моду, распространяющуюся вдоль такой структуры, можно считать поверхностной, а не объемной.
В указанных выше системах может быть существенна роль эффектов, вызванных нелинейностями сред из-за сильного локального увеличения полей (см. ниже). Этот эффект близок к эффекту, возникающему в известной задаче о наклонном падении волн на плазменные слои [8].
Заметим, что нелинейности в метаматериалах могут быть существенно разнообразнее, чем в обычных средах. В последние годы наметился чрезвычайно быстрый прогресс в технологии создания различных композитов, особенно в диапазоне СВЧ. В частности, недавно созданы нелинейные метаматериалы, в которых в диэлектрическую среду внедрены микрочастицы (мета-атомы) с нелинейными характеристиками [9], а именно, частицы, нагруженные микродиодами. Нелинейные свойства таких сред можно достаточно легко менять за счет внешних воздействий. Подобные нелинейные материалы могут представлять интерес при создании различных устройств на основе структур с ПЭМ [10].
1. Линейная задача
Рассмотрим сначала систему с линейными средами. Как известно, нелинейности обычно невелики, поэтому анализ линейной задачи, как правило, может служить отправной точкой для исследования нелинейной задачи. Геометрия изучаемой структуры схематично показана на рис. 1. Предполагаем, что метаматериал расположен в нижней части системы, а обычная среда - в верхней. Действительные части проницаемостей метаматериала отрицательны, а обычных сред - положительны. Предполагаем, что в области \у\ < й расположен промежуточный (переходный) слой, в котором проницаемости е(у) и ^(у) зависят от поперечной координаты у. Ниже основное внимание уделено случаям, когда проницаемости - всюду непрерывные функции у. Вне переходного слоя проницаемости постоянны; при у < —й
диэлектрическую и магнитную проницаемости обозначим через 81 и Ц1, а при у > й - через 82 и ц2.
В работе рассматриваем двумерную задачу, причем исследуем характеристики ТМ мод. Для мод этого типа магнитное поле имеет только одну компоненту Нх. Предполагаем, что процесс монохроматический; поля изучаемых мод зависят от координат у и г, Рис. 1. Геометрия задачи и система координат а также от времени Ь. Как указывалось
выше, граничный слой предполагаем тонким, то есть к! ^ 1. Считаем, что все среды имеют малые диэлектрические и магнитные потери, то есть во всех точках пространства 1т е > 0 и 1т ^ > 0, причем почти везде |1т е/е\ ^ 1 и |1т ^ 1; эти условия будут уточнены ниже. Наличие диссипации в средах в ряде случаев является необходимым условием для существования решения рассматриваемой задачи. Метод расчета ПЭМ в данной структуре во многом повторяет методику, описанную в [11], поэтому опишем его кратко, останавливаясь только на основных моментах.
Поля ПЭМ имеют вид волн, бегущих вдоль оси г; все компоненты полей пропорциональны множителю2 ехр [г($г — ю£)], где в - комплексная константа распространения, ю - частота. Такое представление решений возможно, поскольку в уравнениях Максвелла для линейных сред переменные разделяются. Для решения уравнений Максвелла
Ух Е = гк^Н, Ух Н = —гкеЕ (1)
(к = ю/с - волновое число, с - скорость света в вакууме) удобно ввести две вспомогательные функции
и (у) = Нх, V (у) = г Ех, (2)
которые удовлетворяют системе обыкновенных дифференциальных уравнений
и' = kеV, V' = [в2/(ке) — кц]и. (3)
Здесь и далее штрихами обозначены производные по координате у. Поперечная компонента электрического поля выражается через магнитное поле по формуле
Еу = —Нхв/(ке). (4)
При анализе ПЭМ уравнения (3) должны быть дополнены стандартными условиями спадания полей при у
Дисперсионное уравнение (ДУ) для рассматриваемой моды получали по следующей схеме. С учетом условия на бесконечности магнитное поле вне промежуточного слоя \у\ > ! можно представить в виде
[Бх ехр [рх(у + !)], у < —!, Нх = и = I (5)
В ехр[—р2(у — !)], у >
где рх и р2 - внешние поперечные волновые числа (в общем случае - комплексные), Б\ и В2 - амплитуды поля на границах переходного слоя. Все волновые числа связаны соотношениями
в2 = к2ех + р1 = к2е2^2 + (6)
Для вывода ДУ интегрируем3 систему (3) от точки у = —! до точки у = используя начальные значения, которые следуют из (5),
и (—!) = Бх , V (—!) = и'(—!)/(ке{), (7)
2Как правило, указанный множитель опускаем.
3При численном решении (3) следует учесть, что эта система может быть жесткой.
где и'(—й) = —р1 Вь В конечной точке у = й полученные значения функций и(й) и V (й) приравнивали к тем значениям, которые следуют из (5). В результате приходим кДУ
k82V (й) + Р2и (й)=0. (8)
Напомним, что для ПЭМ поперечные волновые числа должны удовлетворять следующим условиям: И,е Р1 > 0 и И,е р2 > 0 [12].
Опишем очень кратко свойства ПЭМ для системы с резкой границей (й = 0); детальный анализ этого случая можно найти в [6,7]. В такой системе ДУ имеет вид
Р1/81 = —Р2/82. (9)
Для простоты дополнительно предположим, что среды не диссипативные. Из формул (6) и (9) нетрудно получить явные выражения для решений ДУ. В частности, выражения для поперечных волновых чисел имеют вид
р? = к2Д21/[1 — (82/81)2],
(10)
Р2 = — Р1(82/81), Д21 = 82^2 — 81Ц1.
Для рассматриваемого случая, когда среды бездиссипативные и граница раздела резкая, ПЭМ существуют, если выполнены условия
Д21 < 0, 82 > |81|, (11)
Д21 > 0, 82 < \81 \. (12)
При выполнении условий (11) мода будет обратной [1,13], а при выполнении (12) -прямой4.
Рассмотрим теперь случай, когда граница раздела двух сред размыта (то есть й > 0), причем функции 8(у) и ц(у) - непрерывные. В такой системе в поперечной плоскости действительные части проницаемостей изменяются от отрицательных значений в подложке до положительных значений в покрытии, поэтому в одной или нескольких точках (например, в точке у = уо) функции И,е 8(у) и И,е ц(у) должны обращаться в нуль. Диэлектрическая проницаемость входит в знаменатель коэффициента в правой части второго уравнения (3), поэтому, если потери в средах отсутствуют, при у = уо этот коэффициент обращается в бесконечность, и тогда уравнения (3) могут иметь сингулярные решения [14].
Если среды имеют потери, то поля уже не будут сингулярными. Для простоты предположим, что обе функции И,е 8(у) и И,е ц(у), которые меняют знаки на интервале (—й, й), обращаются в нуль в одной и той же точке; для определенности считаем, что координата этой точки уо = 0. В малой окрестности точки уо можно считать, что функции И,е 8(у) и Ие ц(у) изменяются линейно. В этом случае приближенное решение уравнений (3) в окрестности уо можно записать в виде [8,11]
Нх » со + у2[С1 + (со|372)1п(у)], ( )
Еу(у) - —(Р/к)Нх(0)/[у(Ие 8)' + 11т 8)],
4Напомним, что для прямых мод фазовые и групповые скорости направлены в одну сторону, а для обратных - в противоположные.
где Со и С1 - некоторые константы. Из (13) следует, что магнитное поле в окрестности уо конечно. В системе с потерями параметр 1т е = 0, а значит, электрическое поле будет также конечно. Следует учесть, что при малом поглощении в среде амплитуда электрического поля в точке у = 0 будет велика и при увеличении \у\ она будет быстро убывать.
2. Результаты расчетов
Исследуем характеристики ПЭМ, распространяющихся вдоль рассматриваемых структур. Прежде всего опишем структуру полей в системе с резкой границей (при ( = 0); соответствующий профиль И,е е(у) приведен на рис. 2 (сплошная кривая 1). На этом рисунке вертикальные штриховые линии показывают границы переходного слоя и его середину. Приведенные на рисунке для сравнения кривые 2 и 3 получены для структур с нерезкими границами (см. ниже).
На рис. 3 представлены распределения модулей полей \ЕУ(у)\ (кривая 1) и \НХ(у)\ (кривая 2) для системы с резкой границей. Поля нормированы так, что
Их(0) = 1.
Параметры волновода: е1
И-1 = -2 + ¿10-2, и-2 = 1 + ¿Ю-2 и (/X = 0 (X - длина волны в вакууме). В данном примере, как и в рассмотренных ниже, предполагали, что среды дис-сипативные. Поля ПЭМ экспоненциально убывают при удалении от границы раздела сред, а на границе (у = 0) они имеют изломы.
Заметим, что на рис. 3 и на других рисунках данные расчетов представлены в безразмерной форме; кроме того, параметры материалов в рассматриваемых структурах были подобраны так, чтобы полученные результаты можно было легко использовать для оценки характеристик ПЭМ как в оптическом диапазоне, так и в диапазоне СВЧ. В частности, проницаемость покрытия была близка к проницаемости арсенида галлия на оптических частотах, а также к проницаемости керамик на основе окиси алюминия в микроволновом диапазоне. Отметим, что проницаемости подложки, изготовленной из мета-материала, могут изменяться в достаточно широких пределах за счет изменения плотности мета-атомов (внедренных частиц) или их геометрии. Для
= -6.48 + ¿10-2 Яее 10.0 5.0 0
-5.0 -10.0
е2 = 11.664 + ¿10
2
1
2/7 уЗ .... --------
.... ....
-2.0
-1.0
0
1.0
у/й
Рис. 2. Распределения Ие е(у) вблизи границ раздела сред. Вертикальные линии показывают границы слоя
\Еу\ 0.8
0.6
0.4
0.2
1 2
| 1—
—Н
\
\НХ\ 0.8
0.6
0.4
0.2
-2.0 -1.0
0
1.0 у/Х
Рис. 3. Распределения электрического (кривая 1) и магнитного (кривая 2) полей ПЭМ в структуре с кусочно-постоянными профилями проницае-мостей. Вертикальная штриховая линия показывает границу раздела сред
определенности в работе рассмотрены случаи, когда модули показателей преломления подложки и покрытия имели близкие значения. Пока остается не ясным вопрос о том, как на практике реализовать плавные профили проницаемостей. По-видимому, в области низких частот эту проблему решить не очень сложно, но в общем случае вопрос требует специального изучения.
Рассмотрим теперь структуру с переходным слоем (при й > 0). Предполагаем, что внутри слоя при —й < у < й функции 8(у) и ц(у) изменяются линейно, а на концах интервала их значения совпадают со значениями в нижней и верхней средах, соответственно. Для такой структуры при \у\ < й действительная часть е(у) имеет вид
Ие е(у)= + Ыу/й), (14)
где обозначили = Ие (81 + 82)/2, = Ие (е2 — 8^/2. На рис. 2 показано распределение Ие 8(у) (кривая 2) для рассматриваемой структуры. Аналогичные соотношения верны для ц(у). Для простоты предполагали, что мнимые части проницаемостей постоянны на всем бесконечном интервале изменения координаты у. Нули функций Ие 8(у) и Ие ц(у) лежат на интервале (—й, й), но они могут быть смещены относительно друг друга. На рис. 4 приведено распределение электрического поля \ЕУ(у)\ (кривая 1) для структуры с размытой границей и линейными профилями вида (14). На этом рисунке вертикальные штриховые линии показывают границы переходного слоя и его середину. Систему (3) и уравнение (8) решали численно [11]. При расчетах использована следующая нормировка поля: Нх(—й) = 1. Параметры подложки и покрытия были такими же, как и в предыдущей задаче, а толщина тонкого переходного слоя составляла й/Х = 0.0002.
Во всех представленных примерах потери намеренно взяты не очень малыми5 для того, чтобы яснее продемонстрировать структуру полей. Видно, что распределение электрического поля \ЕУ \ имеет очень узкий и высокий пик в окрестности точки у0, где Ие 8(у) равно нулю. Высота этого пика приближенно равна \|ЗНх(0)\/(Мт 8), а полная ширина Дур & 2й1т 8/\^1\. Вне интервала (—й,й) поле \ЕУ(у)\ спадает по экспоненциальному закону. В этой области функция \ЕУ(у)\ мала по сравнению
с ее максимальным значением, поэтому при выбранном нами масштабе кривая 1 на рис. 4 идет практически вдоль оси абсцисс.
Пунктирная кривая 2 на рис. 4 показывает распределение магнитного поля \Нх(у)\ в переходной области. Видно, что на интервале (—й, й) модуль магнитного поля почти постоянен. Вне переходной области поле \Нх(у)\ очень близко к полю в системе с резким изменением проницаемости (см. рис. 3); при |у| > й оно спадает экспоненциально в соответствии с формулами (5). Таким
5 Отметим, что многие метаматериалы имеют достаточно большой тангенс угла потерь.
Рис. 4. Распределения электрического и магнитного полей ПЭМ для структуры с линейными профилями проницаемостей в области переходного слоя. Вертикальные линии показывают границы слоя
образом, наличие переходной области мало изменяет магнитное поле; сглаживание профиля £(y) приводит лишь к сглаживанию излома в распределении \Hx(y)\.
В рассматриваемой системе формируется своеобразная мода, у которой доминирует электрическое поле. Заметим, что в обычных диэлектрических волноводах, как правило, амплитуды полей E и H одного порядка. В данной системе электрическое поле втягивается в область, где диэлектрическая проницаемость мала, за счет чего и происходит его заметный рост. Отметим, что за счет диэлектрических потерь рассматриваемая мода может приобретать некоторые черты волны Ценнека [1].
При варьировании толщины переходного слоя изменяются дисперсионные характеристики исследуемой структуры. На рис. 5 показаны зависимости действительных и мнимых частей безразмерного продольного волнового числа ß/к от безразмерной толщины переходного слоя 2d/L Параметры подложки и покрытия: £1 = -6.48 + i10-2, £2 = 11.664 + i10-2, щ = -2 + i10-2, и-2 = 1 + i10-2. В переходном слое проницаемости изменялись линейно в соответствии с формулой (14). В рассматриваемой области изменения толщины слоя исследуемая ПЭМ является обратной, поэтому значения Re ß < 0.
При d ^ 0 постоянная распространения ß моды в структуре с диффузной границей стремится к значению для структуры с резкой границей. Тем не менее предельный переход при d ^ 0 имеет ряд особенностей, которые отличают структуры, содержащие метаматериалы, от систем из обычных диэлектриков. Предположим, что внутри тонкого переходного слоя проницаемости изменяются линейно в соответствии с (14), и при этом толщина слоя стремится к нулю (d ^ 0). В такой системе амплитуда пика электрического поля слабо зависит от размера слоя; как указывалось ранее, величина максимума имеет порядок max\Ey\ ~ \ßHx(0)\/(kIm £). При уменьшении толщины слоя ширина указанного пика Ayp уменьшается (см. выше); она пропорциональна 2d. Таким образом, при d ^ 0 пик распределения сужается, оставаясь почти постоянным по высоте. Описанный предельный переход напоминает процесс сходимости рядов Фурье разрывных функций (явление Гиббса) [15].
3. Нелинейная задача
Как показано выше, в линейной задаче при плавном изменении проницаемо-стей возможно очень сильное возрастание электрического поля, которое, в свою очередь, может привести к заметным нелинейным эффектам. Рассмотрим такие эффекты, предполагая, что нелинейность керровская [16], то есть диэлектрическая проницаемость £ определяется соотношением
£(y,z)= £(0)(y)+ ae(y)\E (y,z)\2, (15)
где £(0) - линейная часть. Коэффициент нелинейности ae может зависеть от y.
Для линейных задач переменные в уравнениях Максвелла разделяются, поэтому выше их решения искали в виде функций, которые пропорциональны exp [i(ßz — rnt)j. В общем же случае, когда присутствуют одновременно и потери и нелинейности сред, переменные не разделяются и решение становится более сложным [17]. В такой системе за счет потерь амплитуды полей ПЭМ будут убывать
при увеличении г, а следовательно, проницаемость е будет зависеть не только от поперечной координаты у, но и от г (за счет зависящего от г второго слагаемого в формуле (15)). Если потери малы, то есть 1т в ^ к, то изменение проницаемо-стей с ростом г будет небольшим и, самое главное, медленным; в этом случае в первом приближении все характеристики моды можно считать медленно меняющимися функциями осевой координаты г [17]. Процесс распространения мод в такой системе будет напоминать процесс распространения мод в волноводе с линейными средами, у которого медленно меняются его параметры (см., например, [18,19]).
Учитывая сказанное, решение рассматриваемой нелинейной задачи можно искать приближенно в виде бегущих волн с медленно меняющимися вдоль оси г амплитудой, фазовой скоростью и распределением полей. В каждом сечении волновода в первом приближении функции и и V удовлетворяют той же системе уравнений (3); единственное отличие заключается в том, что диэлектрическая проницаемость удовлетворяет соотношению (15), в котором координата г играет роль параметра. Для расчетов, как и выше, удобно в формуле (15) заменить компоненты электрического поля на функции и и V
е(у, г) = е(0)(у) + ае(у)[\в(г)и(у, г)/(ке(у, г))\2 + V(у, г)\2], (16)
где координата г является параметром.
Как следует из приведенной формулы, зависимость е от введенных выше функций \и\ и V\ оказывается неявной, поскольку проницаемость входит в знаменатель второго члена в правой части формулы (16). Следует, однако, учесть, что, как правило, представление диэлектрической проницаемости в таком виде является приближенным. В исходной формуле (15) отброшены малые члены, содержащие \Е(у)\4 и члены с более высокими степенями. Поэтому, оставаясь в рамках малоамплитудного приближения, можно упростить расчеты, подставляя в правую часть (в знаменатель) линейное значение е(у, г) = е(0)(у).
При расчетах считали, что среда нелинейна только в некотором слое ограниченного размера, а именно при \у\ < (п1. Это ограничение не является слишком жестким при анализе ПЭМ, если считать, что (п1 > (, так как в рассматриваемой структуре электрическое поле этой моды локализовано преимущественно в переходном слое, а вне его экспоненциально спадает (см. выше). Далее предполагаем, что коэффициент нелинейности при \у\ < (п1 определяется соотношением
ае(у) = ао[1 - (у/(п1)2], (17)
где ао - значение коэффициента нелинейности на оси.
Исследуем теперь дисперсионные характеристики ПЭМ для нелинейной задачи. Предполагаем, что коэффициент нелинейности очень мал, в том числе даже там, где мала проницаемость е. В этом случае для расчета в можно использовать метод возмущений. Все величины, относящиеся к линейной задаче, отмечаем верхним индексом (0), а величины для нелинейной задачи оставляем без дополнительного индекса. Используя лемму Лоренца, можно получить выражение для сдвига посто-
янной распространения (в — |(0)):
те
(| — |(0))/к = I (г — г(0))[ЕуЕ(°) — ЕхЕ(°)] йу/Рс, (18)
—те
где
те
Рс = I [Еу+ Е^ИХ] йу. (19)
—те
Формула (18) является точной, если в нее подставлять точные значения всех величин. Приближенное выражение для сдвига (| — |(0)) получается, если в эту формулу подставить вместо решения нелинейной задачи соответствующее решение линейной [19] и при этом считать, что (г — г(0)) и ае\Е(0) |2. При вычислениях удобно выразить все компоненты полей через функции и(0) и V (0) . Считая толщину промежуточного слоя малой, после некоторых преобразований получаем явное выражение для сдвига постоянной распространения за счет нелинейности среды
—1
а0\иХ°)(0)\2, (20)
предполагая, что йа\ > 2й. Поскольку магнитное поле в переходном слое меняется очень мало, то для упрощения записи часто будем подставлять вместо иХ°^(у0) величину иХ0)(0).
Из полученного соотношения видно, что для применимости метода возмущений необходимо, чтобы было выполнено условие
Ы|нХ°)(у00))|2 = о(1)(1ш г)3/\г1г2^1^2Р1Р2й2\1/2, (21)
где у00) - координата точки, в которой И,е г(0) = 0. Здесь, как обычно, через о(1) обозначили бесконечно малую величину. В рассматриваемой системе за счет нелинейности сред изменение коэффициента распространения ПЭМ мало. Оценим это смещение для структуры, у которой параметры среды были те же самые, что и для примера, данные которого показаны на рис. 5. При 2й/\ = 0.01 модуль относительной величины сдвига постоянной распространения будет порядка 2.4 • 105\а0\\И0)(0)\2, причем для линейной задачи |(0) и (—3.4 + 0.3г)к.
В рассматриваемой системе за счет нелинейности может происходить несколько интересных явлений, например, смещение максимума электрического поля. Величину эффекта смещения при увеличении передаваемой мощности ПЭМ можно оценить аналитически.
в — |(0) я\р(0)\2№)
|(0)
4к2(1ш г)3
. 2пй
— г-
„(0)Р „(0)Р " Р2 г2 Р1 г1
1
1
Рис. 5. Зависимости действительных и мнимых частей безразмерных волновых чисел от нормированной толщины переходного слоя для линейных профилей г и
Нетрудно показать, что координата максимума поля с высокой степенью точности совпадает со значением корня уравнения И,е е(уо) = 0, где е - нелинейная проницаемость. Учитывая (16), получаем, что за счет нелинейности смещенное значение координаты уо максимума поля равно
(уо - уоо))/( « ^1ао\нХ0)(0)\2, ^ = -\в\2/[?1(к1т е)2],
где уоо) - значение корня уравнения И,е е(о) =0 в линейном случае (когда ао = 0 и е = е(о)), а уо - в нелинейном; параметр определен в (14).
При указанных условиях сдвиг пика существенно меньше его ширины: \ уо -уоо)\ ^ Аур. При рассматриваемом здесь линейном профиле е(у) и при условиях (21) смещение максимума поля при увеличении передаваемой мощности очень мало. Например, для системы с параметрами, которые были использованы при построении рис. 5, и при значении \аоНХ(0)2\ = 10-7 величина смещения имеет порядок 10-3(. Следует, однако, заметить, что для других профилей е(у) смещение может быть намного больше (см. ниже).
4. Слой с малой проницаемостью
Выше исследована структура, у которой в переходном слое функция И,е е(у) была линейной и имела простой нуль. Рассмотрим кратко структуру с другой формой профилей, а именно, когда профили уплощенные6 в области малых И,е е. Начнем анализ опять с линейной задачи. Предположим, что действительная часть проницаемости И,е е имеет вид
{(у/()2И,е е1, -(<у< 0,
(23)
(у/()2И,е е2, 0 < у < (.
Профиль действительной части диэлектрической проницаемости для такой структуры показан на рис. 2 (кривая 3). Для этого случая функция И,е е(у) в точке уо = 0 обращается в нуль вместе со своей первой производной, но следующая производная не равна нулю. Заметим, что разрыв второй производной функции е(у) при у = 0 не оказывает значимого влияния на полученные результаты. Магнитная проницаемость подчиняется сходным соотношениям (с заменой в уравнении вида (23) величины И,е е1 на И,е и т.д.). Как и выше, при расчетах предполагали, что величины 1т е и 1т (1 малы и постоянны. Считали также, что промежуточный слой тонкий, то есть к( < 1.
На рис. 6 (кривая 1) приведено распределение электрического поля в переходном слое и его окрестности для данного примера. Параметры подложки и покрытия были следующие: е1 = -6.48 + ¿10-2, е2 = 11.664 + ¿10-2, = -2 + ¿10-2, = 1 + ¿10-2. В промежуточном слое, толщина которого составляла 2(/Х = 0.0002, проницаемости изменялись непрерывно в соответствии с формулой (23). В отличие
6В последние годы было предложено несколько необычных волноводных структур, в которых применяются диэлектрические слои с малой проницаемостью («zero index metamaterials») [20].
от системы, рассмотренной ранее (см. рис. 4), в данном примере пик распределения \ЕУ(у)\ оказывается достаточно широким. На этом же рисунке показано распределение магнитного поля \Нх(у)\ внутри промежуточного слоя (пунктирная кривая 2). Распределение магнитного поля в данной структуре напоминает поле ПЭМ в предыдущем примере (см. раздел 2).
На рис. 7 показаны зависимости безразмерных действительных и мнимых частей продольного волнового числа в/к от безразмерной толщины переходного слоя 2!/Х. Горизонтальная штриховая линия показывает нулевое значение 1т в/к. Рисунок построен для структуры с описанными выше распределениями проницаемостей (23). Параметры сред были такие же, как и при расчете полей (см. рис. 6). Отметим, что при малых к! рассматриваемая мода является обратной. Дисперсионные зависимости для рассматриваемой структуры отличаются от тех, которые были получены для системы с линейными профилями проницаемостей. В частности, оказывается, что при увеличении толщины слоя 2! очень быстро мода из поверхностной (собственной, направляемой) становится вытекающей (квазисобственной) [12,21]; для вытекающей моды значения И,е р\ < 0 и ее поля растут на бесконечности. Темным треугольником на рис. 7 показана точка, которая определяет границу, правее которой мода становится вытекающей. Заметим также, что при изменении знака 1т в мода из обратной становится прямой.
Исследуем характеристики мод в данной структуре при наличии нелинейно-стей сред. Как и выше, предполагаем, что нелинейность керровская, то есть справедливо соотношение (15). В этой структуре, как и в рассмотренной ранее, за счет нелинейности происходит изменение формы дисперсионных кривых и смещение пика в распределении поля \ЕУ \. Рассмотрим второй эффект. Из (16) следует, что при росте амплитуды полей положение смещенного максимума поля определяется формулой
(уо - у00))/! « ^ао\нХО)(0)\2, ^ = -\eAkVRlm £)\. (24)
При выводе последней формулы считали, что а0 > 0.
11е(Р/1)_
-3.4
-3.5
-3.6
-3.7
2 ..........,. -
1 / 4 -
.........
Лтф/1)
0.10
0.05
-0.05
Рис. 6. Распределения электрического и магнитного полей ПЭМ для структуры с параболическими профилями проницаемостей в области переходного слоя. Вертикальные линии показывают границы слоя
0 0.0002 0.0004 0.0006 0.0008 ЪНХ
Рис. 7. Зависимости действительных и мнимых частей безразмерных волновых чисел от нормированной толщины переходного слоя для параболических профилей £ и
Сравнение соотношений (22) и (24) показывает, что сдвиг максимума поля в системе с уплощенным профилем И,е г больше, чем в случае, когда профиль И,е г линейный; для первого примера сдвиг пропорционален \а0И;2(0)\, тогда как для второго - л/\а0И,2(0)\. Нетрудно показать, что сдвиг увеличится, если профиль сделать еще более уплощенным, например, если в распределении (23) вместо показателя степени «2» взять показатель «4».
Обратим внимание на еще одно свойство рассматриваемой системы, у которой в некоторой точке обращаются в нуль как сама функция И,е г(у), так и ее производная. Как отмечалось выше, при увеличении толщины переходного слоя мода из направляемой может превратиться в вытекающую. За счет нелинейности точка перехода (см. рис. 7) может смещаться, если увеличивать мощность моды. Анализ влияния нелинейностей на характеристики вытекающих мод сложнее, чем для собственных мод; некоторые вопросы этого анализа до сих пор строго не решены. Укажем на одну из проблем. Напомним, что при изучении ПЭМ параметр йп\ можно устремить к бесконечности, то есть считать, что нелинейность присутствует во всем пространстве. Поля же вытекающих мод растут на бесконечности (при \у\ ^ то), поэтому, если размер йп\ достаточно велик, то методы возмущений перестают работать. Учет эффектов насыщения нелинейности не спасает положение. Заметим, что слабовытекающие моды в весьма грубом приближении можно исследовать, пренебрегая радиационными эффектами. Однако такой способ не дает полную картину нелинейных явлений. Для строго анализа, по-видимому, следует решать проблему возбуждения открытого нелинейного волновода, то есть рассчитывать совместное распространение вытекающей моды и мод непрерывного спектра. К сожалению, на этом пути возникает много дополнительных проблем. В частности, характеристики решения задачи возбуждения, а значит, и результаты расчетов нелинейных эффектов могут достаточно сильно зависеть от выбора спектра радиационных мод. Поэтому возникает задача оптимального выбора этого спектра, когда можно было бы в основном получить характеристики вытекающей моды.
Заключение
В работе исследованы волноводные характеристики полубесконечного слоя метаматериала с нерезкой границей. Показано, что в такой структуре электрическое поле ПЭМ типа ТМ локализовано в очень тонком слое и, как правило, принимает гигантские значения. Это свойство может представлять интерес для практических приложений (например, при создании различных микроволновых и оптических датчиков). Локализация поля в таких системах обычно существенно выше, чем в структуре с резкой границей раздела сред.
При наличии нелинейностей распределение электрического поля в переходном слое зависит от передаваемой мощности, что также может быть использовано на практике. В частности, при изменении передаваемой мощности может происходить смещение максимума электрического поля. Отметим, что нелинейные эффекты зависят от формы распределения линейной части диэлектрической проницаемости.
В работе рассмотрены достаточно простые структуры, в частности, предполагали, что мнимые части проницаемостей не зависят от поперечной координаты.
Однако, если профили проницаемостей формируются за счет изменения концентрации одинаковых мета-атомов, то переменным будет также профиль Im е, что, в свою очередь, может заметно изменить линейные и нелинейные характеристики мод. Также при анализе предполагали, что нелинейны только действительные части диэлектрических проницаемостей, то есть коэффициент ae в формуле (15) действительный. Дополнительные возможности возникают, если он является комплексным; в этом случае из-за нелинейности может изменяться поглощение в средах. В подобных структурах при изменении передаваемой мощности могут варьироваться характерные параметры распределения полей (например, ширина пика и его высота), а также соотношение между электрическим и магнитным полями. Отметим также, что нелинейности часто приводят к анизотропии проницаемостей [16]; эти эффекты также могут найти практические применения. Оценки указанных эффектов можно сделать методами, подобными тем, что были использованы выше.
Автор признателен А.Г. Рожневу и А.А. Башарину за обсуждение ряда вопросов, связанных с данной работой.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (гранты 13-02-00732 и 14-02-00976).
Библиографический список
1. Вайнштейн Л.А. Электромагнитные волны. М.: Радио и связь, 1982.
2. Вайнштейн Л.А. Теория дифракции и метод факторизации. М.: Радио и связь,
1966.
3. Поверхностные поляритоны / Ред. В.М. Агранович, Д.Л. Миллс. М.: Наука, 1985.
4. Либенсон М.ДПоверхностные электромагнитные волны оптического диапазона // Соросовский образоват. журнал. 1996. № 10. С. 92.
5. Князев Б.А., Кузьмин А.В. Поверхностные электромагнитные волны: От видимого диапазона до микроволн // Вестник НГУ. Серия: Физика. 2007. Т. 2, № 1. С. 108.
6. Ruppin R. Surface polaritons of a left-handed medium // Phys. Lett. A. 2000. Vol. 277, № 1. P. 61.
7. Shadrivov I.V., Sukhorukov A.A., Kivshar Yu.S., et al. Nonlinear surface waves in left-handed materials // Phys. Rev. E. 2004. Vol. 69, № 1. P. 016617.
8. Гинзбург В.Л.Распространение электромагнитных волн в плазме. М.: Наука,
1967.
9. Kapitanova P. V., Slobozhanyuk A.P., Shadrivov I.V., et al. Competing nonlinearities with metamaterials // Appl. Physics Lett. 2012. Vol. 101. P. 231904.
10. Lapine M., Shadrivov I.V., Kivshar Yu.S. Colloquium: Nonlinear metamaterials // Reviews Modern Physics. 2014. Vol. 86, № 3. P. 1093.
11. Маненков А.Б. Дисперсионные характеристики мод волновода из метаматери-ала // Радиотехника и электроника. 2012. Т. 57, № 9. С. 983.
12. Вайнштейн Л.А., Маненков А.Б. // Сб. Лекции по электронике СВЧ и радиофизике. Саратов: Изд-во Сарат. ун-та, 1986. Кн. 1. С. 141.
13. Шевченко В.В. Прямые и обратные волны: Три определения, их взаимосвязь и условия применимости // Успехи физ. наук. 2007. Т. 177, № 3. С. 301.
14. Трикоми Ф. Дифференциальные уравнения. М.: Изд-во иностр. лит., 1962.
15. Арнольд В.И. Лекции об уравнениях с частными производными. М.: Фазис, 1979.
16. Ogusu K. TM waves guided by nonlinear planar waveguides // IEEE Trans. Microwave Theory Tech. 1989. Vol. 37, № 6. P. 941.
17. Маненков А.Б. Характеристики волновода из нелинейного метаматериала // Изв. вузов. Прикладная нелинейная динамика. 2010. Т. 18, № 6. С. 93.
18. Маненков А.Б. Распространение волн в открытых волноводах с анизотропным диэлектриком // Изв. вузов. Радиофизика. 1981 Т. 24, № 1. С. 84.
19. Никольский В.В. Вариационные методы для внутренних задач электродинамики. М.: Наука, 1967.
20. Basharin A.A., Mavidis C., Kafesaki M. et al. Epsilon near zero based phenomena in metamaterials // Phys. Rev. B. 2013. V. 87, № 15. P. 155130.
21. Маненков А.Б. Условия ортогональности вытекающих мод // Изв. вузов. Радиофизика. 2005. Т. 48. № 5. С. 388.
Поступила в редакцию 10.11.2014 После доработки 15.12.2014
PROPAGATION OF WAVES ALONG DIFFUSE BOUNDARY OF NONLINEAR METAMATERIALS
A. B. Manenkov P. Kapitza Institute for Physical Problems RAS
The characteristics of surface TM modes, guided by diffuse interface between two nonlinear media (metamaterial and conventional dielectric), are investigated. The field distributions and dependencies of propagation constants on parameters of the transition layer are calculated. Effects, caused by media nonlinearity, are studied.
Keywords: Surface electromagnetic modes, nonlinear metamaterials, diffuse boundaries of media interfaces.
Маненков Александр Бенционович - родился в Москве (1943). Окончил Московский физико-технический институт (1965), к.ф.-м.н. (1973), д.ф.-м.н. (1996). С 1965 года работает в Институте физических проблем РАН. В физической лаборатории ИФП под руководством П.Л. Капицы экспериментально и теоретически занимался разработкой микроволновой аппаратуры для плазменных исследований. Совместно с Л.А. Вайнштейном построил общую теорию возбуждения открытых волноводов. Провел (совместно с В.Н. Мелехиным) пионерские исследования открытых волноводов брэгговского типа. В настоящее время - ведущий научный сотрудник ИФП РАН. Область научных интересов -исследование задач возбуждения, распространения и дифракции волн в открытых электродинамических структурах, численный анализ различных задач радиофизики (включая оптику) и электроники.
119334, Москва, ул. Косыгина, 2
Институт физических проблем им. П.Л. Капицы РАН
E-mail: manenkov@kapitza.ras.ru