Изв. вузов «ПНД», т. 17, № 6, 2009
УДК 621.372; 537.86; 537.87
ВОЛНОВОД С ТОНКИМИ нелинейными стенками
А.Б. Маненков
Исследованы характеристики мод планарного диэлектрического волновода с тонкими стенками, изготовленными из нелинейных материалов. Рассмотрен эффект перехода моды из одного множества в другое при варьировании параметров волновода или передаваемой по нему мощности. Показано, что с ростом амплитуды полей возможно «исчезновение» направляемой или вытекающей мод, то есть при высоком уровне мощности, вводимой в волновод, модовое распространение может быть подавлено.
Ключевые слова: Волноводы с нелинейными средами, вытекающие и антиповерхностные моды, преобразование мод
Введение
Материалы с нелинейными характеристиками представляют значительный интерес для применения в схемах интегральной и волоконной оптики, а также в системах СВЧ [1,2]. К настоящему времени исследованы нелинейные диэлектрики с большим разнообразием параметров, которые можно использовать для подобных структур [3,4]. Среди структур, в которых применяются нелинейные материалы, важную роль играют слоистые волноводы [1,5]. Ранее было показано, что слоистые волноводы с линейными средами, в которых могут распространяться вытекающие моды (ВМ), обладают рядом уникальных свойств [6-9]. Можно предположить, что, используя нелинейные материалы, удастся существенно расширить области применения подобных систем; поэтому анализ ВМ в таких структурах представляет как теоретический, так и практический интерес.
В работе [10] описан метод расчета характеристик вытекающих мод в планарных волноводах, изготовленных из нелинейных материалов. Был рассмотрен общий случай многослойных волноводов с достаточно произвольными параметрами; анализ проводился численно. В настоящей работе влияние нелинейностей сред изучено более подробно на одном простом примере - на примере волновода с тонкими стенками. Такая структура является двумерной моделью так называемых диэлектрических трубок (брэгговских волокон) [6-8], а также упрощенной моделью волокон на основе фотонных кристаллов [9] и некоторых плазменных структур. Несмотря на простоту, такая модель очень удобна для анализа многих свойств ВМ [11]. Как будет
3
показано ниже, с ее помощью можно рассмотреть различные явления, происходящие в волноводах с нелинейными средами. В частности, на этом примере удается получить аналитические результаты, которые позволяют детально изучить области параметров, где нелинейности проявляются наиболее сильно, а также исследовать вопрос об изменении структуры и типа мод при варьировании параметров волновода и передаваемой по нему мощности.
Отметим, что рассматриваемая задача тесно связана с задачей об отражении плоской волны от нелинейного слоя [12]. В частности, используя результаты указанной задачи, можно определить характеристики вытекающих мод в приближении геометрической оптики.
1. Исходные уравнения
Рассматриваемый волновод (рис. 1) представляет собой планарную структуру: центральный волноведущий слой, окруженный сверху и снизу нелинейными диэлектрическими слоями малой толщины. Поперечный размер центрального слоя обозначим через 2d, а его постоянную проницаемость - через eg = n2g (ng - показатель преломления). Предполагаем, что проницаемость окружающей среды (выше и ниже стенок) также равна eg. Рассматриваем частный случай симметричного волновода, когда характеристики волновода при у > 0 и у < 0 совпадают. Толщину стенки и ее диэлектрическую проницаемость обозначим через di и ei. В дальнейшем предполагаем, что стенки изготовлены из нелинейного материала с кубической (керровской) нелинейностью [5,13], причем считаем, что амплитуды полей не очень велики; более точные ограничения, накладываемые на параметры задачи, описаны ниже.
Рассмотрим ТЕ-случай, когда электрическое поле имеет только одну компоненту Ex. Исследуем гармонические процессы, при которых временная зависимость имеет вид exp (—гюк), где ю - частота, а t - время. В дальнейшем считаем, что в рассматриваемых средах высшие временные гармоники подавляются. Это приближение оправданно в случае, когда затухание волн на таких гармониках велико.
Поля мод в такой системе должны удовлетворять уравнению [5,13]
dgEx
dyg
+
dgEx
dzg
+ k‘2e(У, \Ex\2)Ex
0,
(1)
где k - волновое число; e - диэлектрическая проницаемость, которая зависит от поперечной координаты у и от квадрата модуля поля \Ex\g. В стенках (покрывающих слоях), изготовленных из нелинейных диэлектриков, проницаемости равны ei =
= ei0 + a\Ex(y)\g, где постоянная ei0
Рис. 1. Геометрия планарного волновода с тонкими стенками
ng - линейная часть (n1 - показатель преломления). Для определенности считаем, что e(0) > 0. Постоянный вещественный коэффициент а может быть как положительным, так и отрицательным. Предполагаем также, что диэлектрические потери на частоте ю отсутствуют, то есть во всех точках пространства Im e = 0. Кроме того, считаем, что магнитная проницаемость всех сред равна единице.
4
В рассматриваемой системе могут распространяться моды различных типов. Ниже будут рассматриваться моды, волновые числа которых являются корнями комплексного дисперсионного уравнения, например, моды дискретного спектра. Заметим, что в волноводах с нелинейными средами, по-видимому, нельзя представить поля в виде разложений по собственным модам, которые были бы аналогичны разложениям в линейных задачах [14]. Поэтому классификация мод в указанных задачах, вообще говоря, может различаться. Однако при малых нелинейностях можно определить тип моды с помощью предельного перехода, полагая а ^ 0 и определяя тип моды в соответствующей линейной системе (см. ниже).
Для направляемых мод (НМ) в случае отсутствия диэлектрических потерь решение уравнения (1) имеет вид бегущих волн
Ex = As<&s(y) exp [*(Ps2 - wt)}, (2)
где |3s - постоянная распространения НМ, а As - вещественный амплитудный множитель; для простоты начальную фазу поля положили равной нулю. Индексом s отмечаем, что мода является направляемой (поверхностной). Отметим, что для НМ имеем Im |3s = 0 и |3s > kng. Для симметричных мод функция Ф^у) удовлетворяет уравнению
^2фв
dy2
+ [k2e(y, As>s|2) - ps2] Фs = 0, Фs(0) = 1,
(3)
и дополнительному условию dФs/dy = 0 при у = 0. В приведенном уравнении проницаемость ^(yjA2^|2) зависит от поперечной координаты, но не зависит от координаты z. Это свойство объясняется тем, что для НМ постоянная распространения |3s вещественна, а значит, модуль поля не зависит от координаты z, то есть в этом случае lExl = |A^s(y)| зависит только от y. Таким образом, для направляемой моды в уравнении (1) переменные разделяются, что и позволило искать решение задачи в виде (2).
Для другого класса мод, а именно для вытекающих (квазисобственных, англ. -quasi-modes), это свойство не имеет места; для них переменные не разделяются, так что решение задачи имеет более сложную форму. Рассмотрим характеристики вытекающей моды низшего номера TEqo. Здесь и ниже добавочным нижним индексом q отмечаем, что мода вытекающая. Ищем приближенное решение задачи, предполагая, что нелинейные эффекты малы (то есть |аЕХ| ^ n\). Поле Eq0 вытекающей моды низшего номера в нелинейных слоях в первом приближении удовлетворяет уравнению [10]
<92EqQ + д 2Eq0 dy2 dz2
+ [k2 (n2 + a|EqQ|2)}EqQ
0.
(4)
Из приведенного уравнения следует, что для полей вытекающих мод диэлектрическая проницаемость зависит как от поперечной координаты y, так и от z, поскольку для ВМ приближенно |Eq0| ~ exp(-zIm|3q0), причем Im|3q0 = 0 (здесь и ниже |3q0 - комплексный коэффициент распространения ВМ типа TEq0). Для определенности исследуем прямые моды, у которых Im |3q0 > 0. Далее предположим, что радиационные потери мод малы, то есть Im |3q0 ^ kng. При указанных выше условиях решение уравнения (1) для вытекающей моды можно записать в следующей приближенной форме:
Eq0 = Aq0(z^q0(y, Aq0(z)) exp {i^q0(z) - wt}}, (5)
5
где Aq0(z) - вещественная амплитудная функция, фчо(<г) - вещественная фаза (фазовый набег). Функция Фч0(у, Aq0(z)) описывает распределение поля ВМ в поперечной плоскости при фиксированном значении z; она является решением уравнения
Л2Фд0
dy2
+ {Л2ф,Адо|Фдо|2) - Рдо}Фд0
0, Фд0(0)
1.
(6)
В этом уравнении амплитуда Ад0 = Ад0(г) является параметром, зависящим от z. Таким образом, для вытекающей моды проницаемость е(у, |Ад0Фд0|2) зависит от двух координат, у и z. В силу медленности спадания модуля поля |Е^д0| можно разбить рассматриваемый интервал оси z на малые интервалы длиной Az, на которых амплитуду ВМ можно считать практически постоянной (при условии Az Im |Зд0^1). В этом случае на каждом малом интервале решение (1) можно искать в виде выражения (5), в котором набег фазы приближенно равен фд0^) w Re ^0(z)z. Далее в окрестности любой плоскости z = const можно определить локальные характеристики вытекающей моды, в частности, фазовый коэффициент Re |Зд0 и коэффициент затухания Im |Зд0, которые для нелинейной задачи зависят от амплитуды A^(z), а значит, и от координаты z. Приведенные выше качественные соображения можно уточнить, применяя метод возмущения, неполный метод Галеркина или методики, похожие на метод поперечных сечений [14-17].
Отметим, что рассматриваемая нелинейная задача похожа на линейную задачу распространения волн в нерегулярном волноводе с заполнением, параметры которого зависят от обеих координат. Однако существует ряд принципиальных различий между указанными двумя классами задач. В линейных задачах изменение показателя преломления вдоль оси z может происходить только за счет «внешнего» изменения параметров диэлектриков. В нелинейной задаче изменение эффективного значения диэлектрической проницаемости происходит за счет самого поля (самовоздействие поля), поскольку оно спадает вдоль оси z. Заметим также, что в нелинейной задаче, для того чтобы при анализе получить достаточно простые уравнения, приходится делать ряд дополнительных предположений, в частности, считать, что нелинейные эффекты малы, так что диэлектрическая проницаемость в нелинейных слоях волновода слабо зависит от координаты z, а амплитуды возбуждающихся высших мод малы. Для линейных задач уравнения связанных мод, которые выводятся, например, методом поперечных сечений, вообще говоря, верны при произвольном соотношении между амплитудами всех мод. Для этих задач требование малости амплитуд часто используется лишь на конечном этапе, если дополнительно применяется метод возмущений. Указанное различие связано, в частности, с тем, что в нелинейных структурах не выполнены условия ортогональности мод разных типов в той форме, в которой они получаются для линейных систем [14,18].
Для фазовой функции и амплитуды основной вытекающей моды типа ТЕд0 из (4), (5) можно получить приближенные уравнения
dA^0
dz
dФq0
dz
Re Рд0,
1 d(Re |3q0)
2 Re |Зд0 dz
+ Im вд0
Ад0,
(7)
(8)
6
которые в малоамплитудном приближении решаются итерациями. Заметим, что система (7)-(8) должна решаться совместно с уравнением (6), которое определяет нелинейную зависимость комплексного коэффициента распространения моды |3qo от амплитуды поля на оси волновода; поэтому в общем случае правые части системы нелинейно зависят от амплитуды Aqo. Заметим также, что в следующем приближении из уравнения (1) получится система уравнений, похожая на систему уравнений для связанных мод [14, 17]. Как и в линейном случае волновода с переменными параметрами, она описывает эффекты преобразования мод различных типов (как в прямые, так и во встречные). Отметим, что при анализе преобразования в открытом волноводе возникают дополнительные трудности, связанные с эффектом возбуждения пространственной волны.
В дальнейшем не будем рассматривать зависимость полей от продольной координаты, а сосредоточимся только на анализе локальных свойств вытекающей моды. Заметим, что для случая, когда длина волновода невелика, эти характеристики играют основную роль при решении многих задач, встречающихся на практике.
2. Локальные характеристики вытекающей моды
Для основной вытекающей моды TEqo решение (6), то есть функцию Ф<^о внутри и вне волновода, можно представить в виде
I cos(Kqoy), при \y\ <d;
Фq0 = < (9)
(Be exp [*Kqo(\y\ - d)}, при \y\ > d,
где Be - амплитуда вытекающей волны в окружающей среде. Волновые числа моды TEqo связаны соотношениями
pqo = k2ng - 4o, (10)
где Kqo - поперечное волновое число. При выводе этих формул учтена симметрия моды TEqo, условия излучения и нормировки (Ф^О) = 1). В частности, из условий излучения следует, что Im Kqo < 0. В силу симметрии системы ниже рассматриваем только верхнюю часть волновода (при y > О).
Считаем, что толщина стенок волновода мала [19], то есть
kd\\n\\ < 1. (11)
Такую тонкую стенку можно заменить импедансной плоскостью, на двух сторонах которой выполнены граничные условия импедансного типа [11,19,20]
E+ = E-, H+ - H+ = -iYwE-, (12)
где E-, H-, E+, H+ - предельные значения полей на нижней и верхней частях импедансной плоскости (при y = d — 0 и y = d + 0), а двухсторонний адмитанс равен
Yw = kd\(e\ - ед), ei = ei° + a\E-\2. (13)
При выводе (12) формально считали, что d1 ^ 0, но так, что параметр Yw остается постоянным. Из (12) следует, что в формуле (9) коэффициент Be = cos (Kd).
7
Вывод дисперсионного уравнения для основной ВМ, распространяющейся в такой структуре, проводится, так же как в линейном случае [11], с помощью формул (9) и (12). Опуская промежуточные выкладки, приведем окончательное соотношение
(v + iKd) exp(-2iKd) = —v, v = v(k) = kdYw/2. (14)
Для упрощения формул опускаем нижние индексы у волновых чисел Kqo и |Зчо. Выражение для функции v(k) можно записать в виде
v(k) = v0 + k2dd1a\Aq0 cos(Kd)|2/2, v0 = k2dd\(ef ° — eg)/2, (15)
где vo - линейная часть выражения для v(k) (при а = 0). Для дальнейшего анализа удобно переписать дисперсионное уравнение в виде
Kd
n
2
1 , / iKd
- ln 1 + ——
2 \ v(k)
(16)
Последнее уравнение в некоторых случаях можно решать итерациями. Например, для линейной задачи (а = 0, v = v0) при условии \v0\ ^ 1 из (16), разлагая логарифм в ряд по степеням 1/v0, получим
Kd
n / 1 1 in
2 V + 2v0 + 4v0 8v0 + " '
(17)
Это соотношение можно обобщить на нелинейный случай, когда \v\ ^ 1. Заметим, что условие \v0 \ ^ 1 (или в общем случае условие \v\ ^ 1) не противоречит неравенству (11); например, величина \v0\ может быть велика, если kd ^ 1 (в области «высоких частот»), или при большой разнице диэлектрических проницаемостей. В таких случаях поле на пленке мало, поскольку Ф<^^) ~ 1/\v0\ ^ 1, а значит, Eq0(d) ~ Aq0/v0; поэтому для таких \v0\ влияние нелинейности материала стенки на характеристики ВМ, как правило, не очень велико.
Из (16) следует, что нелинейность диэлектрика может сильно проявляться, если \v\ < 1. При дополнительном условии \Kd\ < 1 функция Ф<^^) ~ 1, а значит, Eq0(d) ~ Aq0, то есть поле на стенке имеет тот же порядок величины, что и поле на оси. Рассмотрим эту область параметров подробнее, предполагая дополнительно, что v0 > 0 (е1 > eg), то есть считая, что парциальные лучи могут распространяться в стенке. Для линейной задачи оценку величины Kd можно также получить из (16) итерациями. В результате, в первом приближении при v0 > 0 получаем
п — - ln 1 + — ln —
2 ln(1/v0) 2 to о о
0 < v0 < 1.
(18)
Заметим, что величина |v0| может быть мала, например, при очень малой толщине стенки (di ^ 0) и фиксированной разности \ni — ng\. Как видно из приведенных соотношений, при v0 ^ + 0 безразмерное волновое число Kd стремится к (п/2 — iro). «Уход» волнового числа в бесконечность объясняется тем, что при v0 ^ + 0 волновод фактически перестает существовать: стенки становятся практически прозрачными для падающих на них парциальных лучей и, естественно, радиационные потери вытекающей моды становятся большими. Выражение, подобное (18), можно
8
получить и для нелинейного случая. Следует, однако, учесть, что для нелинейной задачи соответствующая формула дает только качественную оценку величины Kd, поскольку при больших отрицательных Im(Kd) моду TEqo уже нельзя считать слабо вытекающей, как предполагали выше. Другой случай, когда v0 < 0, то есть когда парциальные лучи туннелируют [12] сквозь стенку волновода, более сложен; он будет рассмотрен в следующем разделе.
Кроме корней, которые являются поперечными волновыми числами вытекающей моды, уравнение (14) имеет еще два множества корней, соответствующих модам других типов. Для основной направляемой моды поперечное волновое число можно определить соотношением Ksod = ixs, где вспомогательный вещественный параметр Ts > 0 является корнем уравнения
Ts/[1 + exp(-2ts)j = vo + а|Д,о ch (xs)|2(k2ddi/2). (19)
В этом уравнении через Aso обозначено поле на оси. Поля направляемой моды экспоненциально убывают при у ^ ±го, поскольку Ks0 = i|Ks0|. Заметим, что при kd ^ 1 характеристики этой моды легко получить в аналитической форме, если рассматривать данную структуру как два связанных пленочных волновода, расположенных на расстоянии 2d друг от друга (см. рис. 1).
Рассматриваемое дисперсионное уравнение имеет также корни, соответствующие несобственным антиповерхностным модам (АПМ) [21], для которых поперечные волновые числа являются чисто мнимыми, но в отличие от чисел направляемых мод они имеют другой знак мнимых частей. Для моды низшего номера число као (индекс а обозначает АПМ) может быть определено из вещественного уравнения
Kaod = -йа, Та/[1+ exp(2ta)j = vo + а|Д*о ch(Ta)|2(k2 ddi/2), (20)
где параметр та > 0, а Aao - поле на оси для этой моды. В отличие от НМ, поля АПМ растут на бесконечности при у ^ ±го, причем очень быстро, поэтому АПМ не может возбуждаться реальными источниками с конечной мощностью. Несмотря на ее «нефизичность», эта мода важна для интерпретации результатов исследования движения корней дисперсионного уравнения на комплексной плоскости и взаимного преобразования мод различных типов при изменении параметров волноводов1. Как показано ниже, при изменении параметров, а также амплитуды полей, может изменяться тип моды, например, вытекающая мода может преобразоваться в антиповерхностную. Заметим, что этот эффект существует и в линейных структурах [21]; там он происходит, например, при изменении частоты, а не за счет изменения передаваемой мощности.
3. Численные результаты
В этом разделе рассмотрены характеристики различных мод при варьировании параметров стенки и амплитуды полей. Для изучаемой структуры дисперсионное уравнение решали численно с помощью программ из пакета MINPACK [23]. В области |v| ^ 1 (как при положительных, так и при отрицательных v) для начального
1Некоторые свойства АПМ похожи на свойства антисвязанных состояний, которые рассмотрены в работах, посвященных квантовой теории рассеяния [22].
9
приближения при поиске корней использовали значения Kd, вычисленные с помощью двух или трех итераций из асимптотической формулы (17). В области |v| < 1 для задания начального приближения применяли процедуру непрерывного продолжения параметра. А именно, параметр vo менялся с небольшим шагом, начиная со значений |vo | ^ 1, и на каждом последующем шаге в качестве начального приближения использовали значение корня на предыдущем шаге. Заметим, что в нелинейном случае выражение для v(k) содержит член | cos(Kd)|2, который быстро растет при больших Im(Kd), так что при «неудачном» выборе начального приближения возможно аварийное переполнение, которое, как правило, не позволяет проводить вычисления дальше. Волновые числа НМ и АПМ находили из более простых вещественных уравнений (19) и (20). Отметим, что точность расчетов корней была порядка 10_7 или выше.
Для данной задачи поведение мнимых частей корней дисперсионного уравнения является более информативным, поскольку для части рассмотренных выше мод имеем Re(Kd) = 0; поэтому начнем анализ результатов с исследования зависимостей Im(Kd) от параметров задачи. На рис. 2 изображены зависимости безразмерных величин Im(Kd) от безразмерного параметра vo для линейной и нелинейной задач. Буквами q, s и a отмечены кривые, рассчитанные для ВМ, НМ и АПМ, соответственно. Сплошные кривые построены для случая, когда а = 0 (линейная задача). Пунктирные кривые построены для волновода, у которого параметр а = 0.1 (нелинейная задача). Предполагали, что амплитуды полей всех мод на оси волновода равны единице (то есть Aq0 = As0 = Aa0 = 1) и, кроме того, параметр k2ddf/2 = 1.
Рассмотрим сначала линейную задачу. В этом случае при v0 > 0 (то есть при ef0) > eg) ВМ и НМ существуют независимо во всей указанной области параметра v0, и каких-либо изменений типа мод при варьировании адмитанса не происходит. Поведение зависимостей Im(Kd) от параметра v0 можно достаточно легко объяснить из физических соображений (в том числе и из соображений геометрической оптики). В частности, для НМ модуль волнового числа Ks0 монотонно растет с увеличением v0, поскольку в этом случае, как можно показать, поле моды «втягивается» в стенку. Для ВМ величина Im(Kd) убывает при увеличении v0 в соответствии с оценкой (17). Уменьшение Im(Kd) связано с тем, что с ростом v0 углы скольжения парциальных лучей ВМ уменьшаются и, соответственно, уменьшается (по закону Френеля) доля мощности, проходящей сквозь стенку. Качественное поведение дисперсионных кривых для ВМ в области малых положительных v0 было описано выше (см. формулу (18)) и его также можно объяснить из законов геометрической оптики (предполагая kd » 1).
Для случая, когда v0 < 0 (то есть при ef0) < eg), поведение кривых иное. Во-первых, в этой области параметров для линейной задачи (при а = 0) направляемые моды не существуют. Во-вторых, в области небольших отрицательных значений v0 возникает ситуация, когда вместо вытекающих мод образуются две ветви антиповерхностных (см. левую часть рис. 2). В рассматриваемом случае при v0 ~ —0.139 кривая с меткой q раздваивается и корень дисперсионного уравнения переходит из одного множества в другое. Условно эту точку можно назвать точкой отсечки для ВМ. Соответствующее значение поперечного волнового числа в этой точке Kq0d w i(—0.64).
Перейдем теперь к анализу результатов, полученных для волновода с нелинейной стенкой. Результаты расчетов показаны на рис. 2 пунктирными кривыми.
10
Качественно поведение кривых можно понять, если использовать соотношение между v и vo (15). Согласно этому соотношению, значение функции v(k), которое и определяет положение корня дисперсионного уравнения на комплексной плоскости чисел к, получается из величины vo добавлением нелинейного члена, пропорционального a|Aq0 cos(Kd)|2. В той области параметров, где мнимая часть 1ш(к^) не очень велика, а величина | cos(Kd)|2 остается порядка единицы и при изменении a меняется не очень сильно, кривые для нелинейной задачи получаются из кривых для линейной горизонтальным смещением (см. рис. 2), причем величина смещения пропорциональна параметру a. Поскольку в данном примере a > 0, то кривые смещены влево. Эти рассуждения неверны в случае, когда | 1ш(к^)| » 1. При таких параметрах из-за экспоненциального роста функции | cos(Kd)|2 при увеличении | 1ш(к^)| кроме смещения происходит также постепенное искривление кривых, которое видно внизу на рис. 2. Заметим, что в области параметров, где 1ш(к^) < —1, радиационные потери ВМ велики и ее уже нельзя считать слабо вытекающей, как это предполагали выше. Кроме того, при условии | cos (к^)|2 ^ 1 поле в стенке волновода становится большим и нелинейный член a|£^|2 ~ е^, поэтому для диэлектрической проницаемости нелинейной среды нельзя использовать простое выражение (13). По этим причинам указанные области не анализировали.
Опишем теперь зависимости Re(Kd) от vo для тех же случаев. Результаты расчетов этой величины представлены на рис. 3. Параметры задачи и обозначения кривых здесь те же, что и на рис. 2. В соответствии с формулами, приведенными выше, для направляемой и антиповерхностной мод имеем Re(Kd) = 0. Соответствующие значения величин Re(Kd) для этих мод лежат на лучах, идущих вдоль оси абсцисс; для упрощения рисунка они не показаны. Для вытекающей моды (метка q) при v0 ^ ±(Х) имеем Re(Kd) ^ п/2, при v0 ^ —0 получаем Re(Kd) ^ 0, а при
vo ^ +0 имеем Re(Kd) ^ п/2. Эти предельные значения определяют области значений вещественной части Re(Kd) для ВМ в линейном случае. Приведенные выше рассуждения, которые объясняли смещение кривых при изменении параметра нелинейности a, применимы и для объяснения поведения зависимостей Re(Kd) от vo для нелинейной задачи. Исключение составляет только часть пунктирной кривой для ВМ, которая начинается при положительных vo и оборвана в левой части рис. 3. В этой области параметров при vo < 0, как указывалось выше, происходит сильный рост | Тш(к^)|, так что ВМ уже нельзя считать слабо вытекающей и, кроме того,
Рис. 2. Зависимость Тш(к^) от параметра v
безразмерной величины для случая, когда a = 0.1
Рис. 3. Зависимость Re(Kd) от параметра v
безразмерной величины для случая, когда a = 0.1
11
малоамплитудное приближение, которое использовали для расчетов, здесь не работает. Поэтому анализ в этой области не проводили. Заметим, что для направляемой и антиповерхностной мод величины Re(Kd) остаются равными нулю, но смещаются области существования этих мод.
Приведем теперь типичные распределения полей направляемой и вытекающих мод в такой системе. На рис. 4 показаны зависимости \Ex\ от поперечной координаты у для рассмотренных выше мод. Расчеты проводили для волновода с нелинейной стенкой, у которого а = 0.1, амплитуды полей на оси равнялись единице, а параметр k2dd\/2 = 1. Безразмерный параметр vo = 2 при расчете распределения НМ и vo = 10 при расчете ВМ. Значения vo брали разными, чтобы яснее продемонстрировать характерные особенности распределений полей этих мод. Вертикальной штриховой линией показано положение тонкой диэлектрической стенки. Как видно из рисунка, поля НМ убывают при удалении от стенки, поскольку парциальные лучи этой моды распространяются внутри стенок и испытывают полное внутреннее отражение на границах раздела сред. ВМ образованы лучами, которые преимущественно распространяются между стенками, попеременно отражаясь от них. Поля этих мод достаточно велики внутри волновода (при \у\ < d), с наружной стороны стенки они малы, но затем начинают медленно (для слабо вытекающих мод) возрастать при у ^ ±го.
Рис. 4. Распределения полей направляемой (а) и вытекающей (б) мод в поперечной плоскости волновода
Распределение полей антиповерхностных мод, которое здесь не показано, внутри волновода похоже на распределение полей направляемых мод, а вне волновода поля не убывают, как у НМ, а экспоненциально растут (при у ^ ±го) из-за другого знака мнимой части к. Заметим, что как вытекающие моды, так и антиповерхностные являются несобственными, их поля экспоненциально растут в поперечной плоскости, так что формально переносимая ими мощность равна бесконечности. Тем не менее свойства указанных мод существенно различаются [21]. В частности, в отличие от АПМ поля ВМ могут быть выделены из общего поля, возбуждаемого сторонними источниками, и полное поле в большой области пространства с малой погрешностью можно аппроксимировать суммой полей основных ВМ.
Выше был рассмотрен случай, когда параметр нелинейности а > 0. В другом случае, когда а < 0, зависимости Re(Kd) и Im(Kd) от v0 также получаются смещением кривых, рассчитан-
12
ных для линейного случая, но это смещение происходит вправо. Как и в предыдущем примере, такое описание характеристик мод справедливо только до тех пор, пока нелинейные эффекты остаются достаточно малыми, то есть при условии |a[Ao cos (Kd)]2| < 1 (здесь Ao - поле на оси).
Рассмотрим теперь зависимости коэффициентов2 распространения и затухания вытекающей моды низшего номера TEqo от амплитуды поля на оси волновода. На рис. 5 изображены зависимости Re | и Im |, где | = eq0/(kng) - безразмерный коэффициент распространения ВМ, от Aq0. Расчеты проводили для волновода с параметрами а = 0.1, k2ddi/2 = 1 и vo = —0.6. Приведенные кривые типичны для широкой области параметров волновода. Как следует из этого рисунка, при увеличении амплитуды поля Aqo радиационные потери ВМ растут. Такое увеличение потерь можно объяснить, учитывая, что при указанных выше параметрах, когда а > 0, с ростом поля эффективное значение диэлектрической проницаемости стенки ei за счет нелинейного члена увеличивается, а значит, уменьшается модуль разности leg — ei|; за счет сближения проницаемостей всех сред коэффициент прохождения парциальных лучей через стенку растет, что и приводит к увеличению Im(Pqo). Заметим, что при изменении знака а по тем же причинам будет наблюдаться уменьшение потерь ВМ с ростом амплитуды поля.
Поведение коэффициента распространения ВМ от Aqo существенно меняется в узкой области параметров, где возникает корень дисперсионного уравнения, соответствующий АПМ. Подобный случай иллюстрирует рис. 6, на котором представлены действительная и мнимая части величины ^ = в/(кпд). Кривые для этого рисунка построены при а = 0.1, k2dd1/2 = 1 и v0 = —0.3. В такой области параметров мода TEqo существует не при всех Aqo. При увеличении Aqo возникают две АПМ, у которых Im(pao) = 0, то есть эти моды не имеют радиационных потерь. Таким образом, при увеличении Aqo величина Im(Pqo) будет уменьшаться и в точке изменения типа моды pqo станет чисто действительной величиной. Точка отсечки, в которой величина Im(pqo) обратится в ноль, показана на рис. 6 темным кружком. Следует учесть, что поля АПМ растут при у ^ и их нельзя возбудить источни-
Рис. 5. Зависимости безразмерных коэффициентов распространения и затухания вытекающей моды низшего номера TEqo от амплитуды поля на оси волновода Aq0 при v0 = —0.6
Рис. 6. Зависимости безразмерных коэффициентов распространения и затухания вытекающей моды низшего номера TEq0 от амплитуды поля на оси волновода Aq0 при v0 = —0.3
2Для нелинейной задачи величина pq0 зависит от z, поэтому эту величину нельзя называть постоянной распространения, как это принято для линейных систем.
13
ком с конечной мощностью. Поэтому можно предположить, что в указанной выше точке, в которой корень дисперсионного уравнения переходит из одного множества (класса) в другое, модовое распространение излучения (на симметричных модах) будет практически невозможно. Делая этот вывод, следует, конечно, учесть, что в такой ситуации есть возможность возбуждения ВМ высших номеров, которые, как правило, имеют большие потери, а также возможность возбуждения пространственной волны. Детальный анализ этого эффекта должен включать исследование задачи возбуждения и распространения мод вдоль оси волновода, что выходит за рамки данного исследования. Отметим также, что в рассматриваемой системе существуют другие точки, в которых возможен переход моды из одного класса в другой. Например, в окрестности точки Vo = 0 направляемая мода может трансформироваться в антиповерхностную.
Заключение
В работе рассмотрены локальные характеристики симметричных мод, распространяющихся вдоль волновода с тонкими стенками, изготовленными из материала с нелинейной диэлектрической проницаемостью. Задача решалась в приближении, близком к хорошо известному малоамплитудному приближению.
Расчеты показали, что в такой системе при варьировании ее параметров возможны эффекты преобразования типов мод. Под таким преобразованием понимается эффект изменения основных характеристик мод, по которым они классифицируются, например, изменение знака мнимой части поперечного волнового числа. Можно сказать, что в этом случае мода переходит из одного класса (множества), к которому она вначале принадлежала, в другой класс. Подобный переход может происходить также при изменении передаваемой по волноводу мощности. При этом вытекающая мода может «превращаться» в антиповерхностную. Поскольку антиповерхностные моды не могут возбуждаться реальными источниками, то в области параметров, где корень дисперсионного уравнения соответствует именно такой моде, возможна ситуация, когда в данной структуре модовое распространение волн будет сильно подавлено. В линейной задаче при таких параметрах основная доля излучения уйдет в непрерывный спектр (пространственную волну) и, возможно, в незначительной доле в какие-то моды высших типов. На этом эффекте могут быть сконструированы датчики или элементы схем интегральной оптики.
В работе рассмотрены характеристики мод в волноводе с тонкими диэлектрическими стенками. Описанный подход может быть обобщен на случай, когда стенка волновода образована металлическими решетками (с нелинейной проводимостью) или на случай, когда волновод представляет собой слоистую плазменную структуру. Во всех подобного рода задачах дополнительно нелинейной может быть также среда центрального волноведущего слоя.
Автор признателен И.А. Молоткову и А.Г. Рожневу за обсуждения вопросов, связанных с данной работой.
Работа выполнена при поддержке РФФИ (гранты 08-02-00621 и 08-02-90002).
14
Библиографический список
1. Ахмедиев Н.Н., Анкевич А. Солитоны. М.: Физматлит, 2003.
2. Кившарь Ю.С., Агравал Г.П. Оптические солитоны. М.: Физматлит, 2005.
3. Ryasnyansky A.I., Palpant B., et al. Nonlinear optical properties of copper nanoparticles synthesized in indium tin oxide matrix by ion implantation // J. Opt. Soc. Am. B. 2006. Vol. 23, № 7. P. 1348.
4. Сотский А.Б., Хомченко А.В., Сотская Л.И. Измерение параметров кубичнонелинейных волноводов // Письма в ЖТФ. 1994. Т. 20, № 6. С. 49.
5. Ogusu K. Analysis of non-linear multilayer waveguides with Kerr-like permittivities // Opt. Quantum Electron. 1989. Vol. 21, № 2. P. 109.
6. Мелехин В.Н., Маненков А.Б. Диэлектрические трубы как волноводы с малыми потерями // ЖТФ. 1968. Т. 38, № 12. C. 2113.
7. Маненков А.Б. Затухание быстрых волн в диэлектрических трубах // Радиотехника и электроника. 1977. Т. 22, № 10. С. 2043.
8. Li J., Chiang K.S. Disappearance of modes in planar Bragg waveguides // Optics Lett. 2007. Vol. 32, № 16. P. 2369.
9. Joannopoulos J.D, Johnson S.G., WinnJ.N., Meade R.D. Photonic Crystals: Molding the Flow of Light. Princeton: Princeton University Press, 2008 (2-nd edition).
10. Маненков А.Б. Вытекающие моды в многослойном волноводе с нелинейными диэлектриками // Изв. вузов. Прикладная нелинейная динамика. 2008. Т. 16, № 4. С. 20.
11. Маненков А.Б. Возбуждение быстрых волн в открытом волноводе с диэлектрической стенкой // Изв. вузов. Радиофизика. 1975. Т. 18, № 7. С. 1025.
12. Молотков И.А., Маненков А. Б. О нелинейных туннельных эффектах // Радиотехника и электроника. 2007. Т. 52, № 7. С. 799.
13. Ландау Л.Д., Лифшиц Е.М.Электродинамика сплошных сред. М.: Наука, 1992.
14. Маненков А.Б. Возбуждение открытых однородных волноводов // Изв. вузов. Радиофизика. 1970. Т. 13, № 5. С. 739.
15. Никольский В.В. Вариационные методы для внутренних задач электродинамики. М.: Наука, 1967.
16. Вайнштейн Л.А. Теория дифракции. Электроника СВЧ. М.: Радио и связь, 1995.
17. Каценеленбаум Б.З. Теория волноводов с медленно меняющимися параметрами. М.: АН СССР, 1961.
18. Маненков А.Б. Условия ортогональности вытекающих мод // Изв. вузов. Радиофизика. 2005. Т. 48, № 5. С. 388.
19. Вайнштейн Л.А. Теория дифракции и метод факторизации. М.: Радио и связь, 1966.
20. Глущенко А.Г. Теория волноведущих структур СВЧ с нелинейными пленками // Изв. вузов. Радиофизика. 1988. Т. 31, № 9. С. 1098.
21. Вайнштейн Л.А. Электромагнитные волны. М.: Радио и связь, 1982.
15
22. Ньютон Р. Теория рассеяния волн и частиц. М.: Мир, 1969.
23. Каханер Д., Моулер К., Нэш С. Численные методы и программное обеспечение. М.: Мир, 1998.
Поступила в редакцию 11.03.2009 После доработки 11.04.2009
WAVEGUIDE WITH THIN NONLINEAR WALLS
A.B. Manenkov
Characteristics of modes of the planar waveguide with thin walls, which are made from nonlinear media, are studied. The effect of the conversion of a mode from one set to another is considered. This effect may occur if parameters of the waveguide or transmitted power are varied. It is shown that the effect of the guided or leaky modes «disappearance» can arise when the field amplitudes increase, i.e., high powers launched into such a structure can suppress the modal propagation.
Keywords: Waveguides with nonlinear media, leaky and antisurface modes, mode conversion.
Маненков Александр Бенционович - родился в Москве (1943). Окончил Московский физико-технический институт (1965), к.ф.-м.н. (1973), д.ф.-м.н. (1996). С 1965 года работает в Институте физических проблем РАН. В физической лаборатории ИФП под руководством П.Л. Капицы экспериментально и теоретически занимался разработкой микроволновой аппаратуры для плазменных исследований. Совместно с Л.А. Вайнштейном построил общую теорию возбуждения открытых волноводов. В настоящее время - ведущий научный сотрудник ИФП РАН. Область научных интересов - анализ задач возбуждения, распространения и дифракции волн в открытых электродинамических структурах, численное исследование различных задач радиофизики (включая оптику) и электроники.
119334 Москва, ул. Косыгина, 2
Институт физических проблем им. П.Л. Капицы РАН E-mail: manenkov@kapitza.ras.ru
16