WWW.KNIGA.LIB-I.RU
БЕСПЛАТНАЯ  ИНТЕРНЕТ  БИБЛИОТЕКА - Онлайн материалы
 

«Изв. вузов «ПНД», т. 17, № 6, 2009 УДК 621.372; 537.86; 537.87 ВОЛНОВОД С ТОНКИМИ НЕЛИНЕЙНЫМИ СТЕНКАМИ А.Б. Маненков Исследованы ...»

Изв. вузов «ПНД», т. 17, № 6, 2009 УДК 621.372; 537.86; 537.87

ВОЛНОВОД С ТОНКИМИ НЕЛИНЕЙНЫМИ СТЕНКАМИ

А.Б. Маненков

Исследованы характеристики мод планарного диэлектрического волновода с тонкими стенками, изготовленными из нелинейных материалов. Рассмотрен эффект перехода

моды из одного множества в другое при варьировании параметров волновода или передаваемой по нему мощности. Показано, что с ростом амплитуды полей возможно «исчезновение» направляемой или вытекающей мод, то есть при высоком уровне мощности, вводимой в волновод, модовое распространение может быть подавлено.

Ключевые слова: Волноводы с нелинейными средами, вытекающие и антиповерхностные моды, преобразование мод Введение Материалы с нелинейными характеристиками представляют значительный интерес для применения в схемах интегральной и волоконной оптики, а также в системах СВЧ [1, 2]. К настоящему времени исследованы нелинейные диэлектрики с большим разнообразием параметров, которые можно использовать для подобных структур [3, 4]. Среди структур, в которых применяются нелинейные материалы, важную роль играют слоистые волноводы [1, 5]. Ранее было показано, что слоистые волноводы с линейными средами, в которых могут распространяться вытекающие моды (ВМ), обладают рядом уникальных свойств [6–9]. Можно предположить, что, используя нелинейные материалы, удастся существенно расширить области применения подобных систем; поэтому анализ ВМ в таких структурах представляет как теоретический, так и практический интерес.



В работе [10] описан метод расчета характеристик вытекающих мод в планарных волноводах, изготовленных из нелинейных материалов. Был рассмотрен общий случай многослойных волноводов с достаточно произвольными параметрами; анализ проводился численно. В настоящей работе влияние нелинейностей сред изучено более подробно на одном простом примере – на примере волновода с тонкими стенками. Такая структура является двумерной моделью так называемых диэлектрических трубок (брэгговских волокон) [6–8], а также упрощенной моделью волокон на основе фотонных кристаллов [9] и некоторых плазменных структур. Несмотря на простоту, такая модель очень удобна для анализа многих свойств ВМ [11]. Как будет показано ниже, с ее помощью можно рассмотреть различные явления, происходящие в волноводах с нелинейными средами. В частности, на этом примере удается получить аналитические результаты, которые позволяют детально изучить области параметров, где нелинейности проявляются наиболее сильно, а также исследовать вопрос об изменении структуры и типа мод при варьировании параметров волновода и передаваемой по нему мощности.

Отметим, что рассматриваемая задача тесно связана с задачей об отражении плоской волны от нелинейного слоя [12]. В частности, используя результаты указанной задачи, можно определить характеристики вытекающих мод в приближении геометрической оптики.

–  –  –

В рассматриваемой системе могут распространяться моды различных типов.

Ниже будут рассматриваться моды, волновые числа которых являются корнями комплексного дисперсионного уравнения, например, моды дискретного спектра. Заметим, что в волноводах с нелинейными средами, по-видимому, нельзя представить поля в виде разложений по собственным модам, которые были бы аналогичны разложениям в линейных задачах [14]. Поэтому классификация мод в указанных задачах, вообще говоря, может различаться. Однако при малых нелинейностях можно определить тип моды с помощью предельного перехода, полагая 0 и определяя тип моды в соответствующей линейной системе (см. ниже).





Для направляемых мод (НМ) в случае отсутствия диэлектрических потерь решение уравнения (1) имеет вид бегущих волн Ex = As s (y) exp [i(s z t)], (2) где s – постоянная распространения НМ, а As – вещественный амплитудный множитель; для простоты начальную фазу поля положили равной нулю. Индексом s отмечаем, что мода является направляемой (поверхностной). Отметим, что для НМ имеем Im s = 0 и s kng. Для симметричных мод функция s (y) удовлетворяет уравнению d2 s + k 2 (y, A2 |s |2 ) 2 s = 0, s (0) = 1, (3) s s dy 2 и дополнительному условию ds /dy = 0 при y = 0. В приведенном уравнении проницаемость (y, A2 |s |2 ) зависит от поперечной координаты, но не зависит от s координаты z. Это свойство объясняется тем, что для НМ постоянная распространения s вещественна, а значит, модуль поля не зависит от координаты z, то есть в этом случае |Ex | = |As s (y)| зависит только от y. Таким образом, для направляемой моды в уравнении (1) переменные разделяются, что и позволило искать решение задачи в виде (2).

Для другого класса мод, а именно для вытекающих (квазисобственных, англ. – quasi-modes), это свойство не имеет места; для них переменные не разделяются, так что решение задачи имеет более сложную форму. Рассмотрим характеристики вытекающей моды низшего номера ТЕq0. Здесь и ниже добавочным нижним индексом q отмечаем, что мода вытекающая. Ищем приближенное решение задачи, предпоn2 ). Поле Eq0 вытекающей лагая, что нелинейные эффекты малы (то есть |Ex | 1 моды низшего номера в нелинейных слоях в первом приближении удовлетворяет уравнению [10] 2 Eq0 2 Eq0 + [k 2 (n2 + |Eq0 |2 )]Eq0 = 0.

+ (4) y 2 z 2 Из приведенного уравнения следует, что для полей вытекающих мод диэлектрическая проницаемость зависит как от поперечной координаты y, так и от z, поскольку для ВМ приближенно |Eq0 | exp (z Im q0 ), причем Im q0 = 0 (здесь и ниже q0 – комплексный коэффициент распространения ВМ типа ТЕq0 ). Для определенности исследуем прямые моды, у которых Im q0 0. Далее предположим, что радиационные потери мод малы, то есть Im q0 kng. При указанных выше условиях решение уравнения (1) для вытекающей моды можно записать в следующей приближенной форме:

Eq0 = Aq0 (z)q0 (y, Aq0 (z)) exp {i[q0 (z) t]}, (5) где Aq0 (z) – вещественная амплитудная функция, q0 (z) – вещественная фаза (фазовый набег). Функция q0 (y, Aq0 (z)) описывает распределение поля ВМ в поперечной плоскости при фиксированном значении z; она является решением уравнения d2 q0 + {k 2 (y, A2 |q0 |2 ) 2 }q0 = 0, q0 (0) = 1. (6) q0 q0 dy 2 В этом уравнении амплитуда Aq0 = Aq0 (z) является параметром, зависящим от z.

Таким образом, для вытекающей моды проницаемость (y, |A2 q0 |2 ) зависит от q0 двух координат, y и z. В силу медленности спадания модуля поля |Eq0 | можно разбить рассматриваемый интервал оси z на малые интервалы длиной z, на которых амплитуду ВМ можно считать практически постоянной (при условии z Im q0 1).

В этом случае на каждом малом интервале решение (1) можно искать в виде выражения (5), в котором набег фазы приближенно равен q0 (z) Re q0 (z)z. Далее в окрестности любой плоскости z = const можно определить локальные характеристики вытекающей моды, в частности, фазовый коэффициент Re q0 и коэффициент затухания Im q0, которые для нелинейной задачи зависят от амплитуды Aq0 (z), а значит, и от координаты z. Приведенные выше качественные соображения можно уточнить, применяя метод возмущения, неполный метод Галеркина или методики, похожие на метод поперечных сечений [14–17].

Отметим, что рассматриваемая нелинейная задача похожа на линейную задачу распространения волн в нерегулярном волноводе с заполнением, параметры которого зависят от обеих координат. Однако существует ряд принципиальных различий между указанными двумя классами задач. В линейных задачах изменение показателя преломления вдоль оси z может происходить только за счет «внешнего» изменения параметров диэлектриков. В нелинейной задаче изменение эффективного значения диэлектрической проницаемости происходит за счет самого поля (самовоздействие поля), поскольку оно спадает вдоль оси z. Заметим также, что в нелинейной задаче, для того чтобы при анализе получить достаточно простые уравнения, приходится делать ряд дополнительных предположений, в частности, считать, что нелинейные эффекты малы, так что диэлектрическая проницаемость в нелинейных слоях волновода слабо зависит от координаты z, а амплитуды возбуждающихся высших мод малы. Для линейных задач уравнения связанных мод, которые выводятся, например, методом поперечных сечений, вообще говоря, верны при произвольном соотношении между амплитудами всех мод. Для этих задач требование малости амплитуд часто используется лишь на конечном этапе, если дополнительно применяется метод возмущений. Указанное различие связано, в частности, с тем, что в нелинейных структурах не выполнены условия ортогональности мод разных типов в той форме, в которой они получаются для линейных систем [14, 18].

Для фазовой функции и амплитуды основной вытекающей моды типа ТЕq0 из (4), (5) можно получить приближенные уравнения

–  –  –

которые в малоамплитудном приближении решаются итерациями. Заметим, что система (7)–(8) должна решаться совместно с уравнением (6), которое определяет нелинейную зависимость комплексного коэффициента распространения моды q0 от амплитуды поля на оси волновода; поэтому в общем случае правые части системы нелинейно зависят от амплитуды Aq0. Заметим также, что в следующем приближении из уравнения (1) получится система уравнений, похожая на систему уравнений для связанных мод [14, 17]. Как и в линейном случае волновода с переменными параметрами, она описывает эффекты преобразования мод различных типов (как в прямые, так и во встречные). Отметим, что при анализе преобразования в открытом волноводе возникают дополнительные трудности, связанные с эффектом возбуждения пространственной волны.

В дальнейшем не будем рассматривать зависимость полей от продольной координаты, а сосредоточимся только на анализе локальных свойств вытекающей моды.

Заметим, что для случая, когда длина волновода невелика, эти характеристики играют основную роль при решении многих задач, встречающихся на практике.

–  –  –

где q0 – поперечное волновое число. При выводе этих формул учтена симметрия моды ТЕq0, условия излучения и нормировки (q0 (0) = 1). В частности, из условий излучения следует, что Im q0 0. В силу симметрии системы ниже рассматриваем только верхнюю часть волновода (при y 0).

Считаем, что толщина стенок волновода мала [19], то есть

–  –  –

Вывод дисперсионного уравнения для основной ВМ, распространяющейся в такой структуре, проводится, так же как в линейном случае [11], с помощью формул (9) и (12). Опуская промежуточные выкладки, приведем окончательное соотношение

–  –  –

Заметим, что величина |0 | может быть мала, например, при очень малой толщине стенки (d1 0) и фиксированной разности |n1 ng |. Как видно из приведенных соотношений, при 0 + 0 безразмерное волновое число d стремится к (/2 i).

«Уход» волнового числа в бесконечность объясняется тем, что при 0 + 0 волновод фактически перестает существовать: стенки становятся практически прозрачными для падающих на них парциальных лучей и, естественно, радиационные потери вытекающей моды становятся большими. Выражение, подобное (18), можно получить и для нелинейного случая. Следует, однако, учесть, что для нелинейной задачи соответствующая формула дает только качественную оценку величины d, поскольку при больших отрицательных Im(d) моду ТЕq0 уже нельзя считать слабо вытекающей, как предполагали выше. Другой случай, когда 0 0, то есть когда парциальные лучи туннелируют [12] сквозь стенку волновода, более сложен; он будет рассмотрен в следующем разделе.

Кроме корней, которые являются поперечными волновыми числами вытекающей моды, уравнение (14) имеет еще два множества корней, соответствующих модам других типов. Для основной направляемой моды поперечное волновое число можно определить соотношением s0 d = is, где вспомогательный вещественный параметр s 0 является корнем уравнения

s /[1 + exp (2s )] = 0 + |As0 ch (s )|2 (k 2 dd1 /2). (19)

В этом уравнении через As0 обозначено поле на оси. Поля направляемой моды экспоненциально убывают при y ±, поскольку s0 = i|s0 |. Заметим, что при kd 1 характеристики этой моды легко получить в аналитической форме, если рассматривать данную структуру как два связанных пленочных волновода, расположенных на расстоянии 2d друг от друга (см. рис. 1).

Рассматриваемое дисперсионное уравнение имеет также корни, соответствующие несобственным антиповерхностным модам (АПМ) [21], для которых поперечные волновые числа являются чисто мнимыми, но в отличие от чисел направляемых мод они имеют другой знак мнимых частей.

Для моды низшего номера число a0 (индекс а обозначает АПМ) может быть определено из вещественного уравнения a /[1 + exp (2a )] = 0 + |Aa0 ch (a )|2 (k 2 dd1 /2), a0 d = ia, (20) где параметр a 0, а Aa0 – поле на оси для этой моды. В отличие от НМ, поля АПМ растут на бесконечности при y ±, причем очень быстро, поэтому АПМ не может возбуждаться реальными источниками с конечной мощностью. Несмотря на ее «нефизичность», эта мода важна для интерпретации результатов исследования движения корней дисперсионного уравнения на комплексной плоскости и взаимного преобразования мод различных типов при изменении параметров волноводов1.

Как показано ниже, при изменении параметров, а также амплитуды полей, может изменяться тип моды, например, вытекающая мода может преобразоваться в антиповерхностную. Заметим, что этот эффект существует и в линейных структурах [21];

там он происходит, например, при изменении частоты, а не за счет изменения передаваемой мощности.

3. Численные результаты

В этом разделе рассмотрены характеристики различных мод при варьировании параметров стенки и амплитуды полей. Для изучаемой структуры дисперсионное уравнение решали численно с помощью программ из пакета MINPACK [23]. В области || 1 (как при положительных, так и при отрицательных ) для начального Некоторые свойства АПМ похожи на свойства антисвязанных состояний, которые рассмотрены в работах, посвященных квантовой теории рассеяния [22].

приближения при поиске корней использовали значения d, вычисленные с помощью двух или трех итераций из асимптотической формулы (17). В области || 1 для задания начального приближения применяли процедуру непрерывного продолжения параметра. А именно, параметр 0 менялся с небольшим шагом, начиная со значений |0 | 1, и на каждом последующем шаге в качестве начального приближения использовали значение корня на предыдущем шаге. Заметим, что в нелинейном случае выражение для () содержит член | cos(d)|2, который быстро растет при больших Im(d), так что при «неудачном» выборе начального приближения возможно аварийное переполнение, которое, как правило, не позволяет проводить вычисления дальше. Волновые числа НМ и АПМ находили из более простых вещественных уравнений (19) и (20). Отметим, что точность расчетов корней была порядка 107 или выше.

Для данной задачи поведение мнимых частей корней дисперсионного уравнения является более информативным, поскольку для части рассмотренных выше мод имеем Re(d) 0; поэтому начнем анализ результатов с исследования зависимостей Im(d) от параметров задачи. На рис. 2 изображены зависимости безразмерных величин Im(d) от безразмерного параметра 0 для линейной и нелинейной задач.

Буквами q, s и a отмечены кривые, рассчитанные для ВМ, НМ и АПМ, соответственно. Сплошные кривые построены для случая, когда = 0 (линейная задача). Пунктирные кривые построены для волновода, у которого параметр = 0.1 (нелинейная задача). Предполагали, что амплитуды полей всех мод на оси волновода равны единице (то есть Aq0 = As0 = Aa0 = 1) и, кроме того, параметр k 2 dd1 /2 = 1.

Рассмотрим сначала линейную задачу. В этом случае при 0 0 (то есть при (0) 1 g ) ВМ и НМ существуют независимо во всей указанной области параметра 0, и каких-либо изменений типа мод при варьировании адмитанса не происходит.

Поведение зависимостей Im(d) от параметра 0 можно достаточно легко объяснить из физических соображений (в том числе и из соображений геометрической оптики). В частности, для НМ модуль волнового числа s0 монотонно растет с увеличением 0, поскольку в этом случае, как можно показать, поле моды «втягивается»

в стенку. Для ВМ величина Im(d) убывает при увеличении 0 в соответствии с оценкой (17). Уменьшение Im(d) связано с тем, что с ростом 0 углы скольжения парциальных лучей ВМ уменьшаются и, соответственно, уменьшается (по закону Френеля) доля мощности, проходящей сквозь стенку. Качественное поведение дисперсионных кривых для ВМ в области малых положительных 0 было описано выше (см. формулу (18)) и его также можно объяснить из законов геометрической оптики (предполагая kd 1).

(0) Для случая, когда 0 0 (то есть при 1 g ), поведение кривых иное.

Во-первых, в этой области параметров для линейной задачи (при = 0) направляемые моды не существуют. Во-вторых, в области небольших отрицательных значений 0 возникает ситуация, когда вместо вытекающих мод образуются две ветви антиповерхностных (см. левую часть рис. 2). В рассматриваемом случае при

0 0.139 кривая с меткой q раздваивается и корень дисперсионного уравнения переходит из одного множества в другое. Условно эту точку можно назвать точкой отсечки для ВМ. Соответствующее значение поперечного волнового числа в этой точке q0 d i(0.64).

Перейдем теперь к анализу результатов, полученных для волновода с нелинейной стенкой. Результаты расчетов показаны на рис. 2 пунктирными кривыми.

Качественно поведение кривых можно понять, если использовать соотношение между и 0 (15). Согласно этому соотношению, значение функции (), которое и определяет положение корня дисперсионного уравнения на комплексной плоскости чисел, получается из величины 0 добавлением нелинейного члена, пропорционального |Aq0 cos(d)|2. В той области параметров, где мнимая часть Im(d) не очень велика, а величина | cos(d)|2 остается порядка единицы и при изменении меняется не очень сильно, кривые для нелинейной задачи получаются из кривых для линейной горизонтальным смещением (см. рис. 2), причем величина смещения пропорциональна параметру. Поскольку в данном примере 0, то кривые смещены влево. Эти рассуждения неверны в случае, когда | Im(d)| 1. При таких параметрах из-за экспоненциального роста функции | cos(d)|2 при увеличении | Im(d)| кроме смещения происходит также постепенное искривление кривых, которое видно внизу на рис. 2. Заметим, что в области параметров, где Im(d) 1, радиационные потери ВМ велики и ее уже нельзя считать слабо вытекающей, как это предполагали выше. Кроме того, при условии | cos (d)|2 1 поле в стенке волновода становится 2 (0), поэтому для диэлектрической пронибольшим и нелинейный член |Ex | 1 цаемости нелинейной среды нельзя использовать простое выражение (13). По этим причинам указанные области не анализировали.

Опишем теперь зависимости Re(d) от 0 для тех же случаев. Результаты расчетов этой величины представлены на рис. 3. Параметры задачи и обозначения кривых здесь те же, что и на рис. 2. В соответствии с формулами, приведенными выше, для направляемой и антиповерхностной мод имеем Re(d) 0. Соответствующие значения величин Re(d) для этих мод лежат на лучах, идущих вдоль оси абсцисс; для упрощения рисунка они не показаны. Для вытекающей моды (метка q) при 0 ± имеем Re(d) /2, при 0 0 получаем Re(d) 0, а при 0 +0 имеем Re(d) /2. Эти предельные значения определяют области значений вещественной части Re(d) для ВМ в линейном случае. Приведенные выше рассуждения, которые объясняли смещение кривых при изменении параметра нелинейности, применимы и для объяснения поведения зависимостей Re(d) от 0 для нелинейной задачи. Исключение составляет только часть пунктирной кривой для ВМ, которая начинается при положительных 0 и оборвана в левой части рис. 3.

В этой области параметров при 0 0, как указывалось выше, происходит сильный рост | Im(d)|, так что ВМ уже нельзя считать слабо вытекающей и, кроме того, Рис. 2. Зависимость безразмерной величины Рис. 3. Зависимость безразмерной величины Im(d) от параметра 0 для случая, когда = 0.1 Re(d) от параметра 0 для случая, когда = 0.1

–  –  –

ных для линейного случая, но это смещение происходит вправо. Как и в предыдущем примере, такое описание характеристик мод справедливо только до тех пор, пока нелинейные эффекты остаются достаточно малыми, то есть при условии |[A0 cos (d)]2 | 1 (здесь A0 – поле на оси).

Рассмотрим теперь зависимости коэффициентов2 распространения и затухания вытекающей моды низшего номера ТЕq0 от амплитуды поля на оси волновода.

На рис. 5 изображены зависимости Re и Im, где = q0 /(kng ) – безразмерный коэффициент распространения ВМ, от Aq0. Расчеты проводили для волновода с параметрами = 0.1, k 2 dd1 /2 = 1 и 0 = 0.6. Приведенные кривые типичны для широкой области параметров волновода. Как следует из этого рисунка, при увеличении амплитуды поля Aq0 радиационные потери ВМ растут. Такое увеличение потерь можно объяснить, учитывая, что при указанных выше параметрах, когда 0, с ростом поля эффективное значение диэлектрической проницаемости стенки 1 за счет нелинейного члена увеличивается, а значит, уменьшается модуль разности |g 1 |;

за счет сближения проницаемостей всех сред коэффициент прохождения парциальных лучей через стенку растет, что и приводит к увеличению Im(q0 ). Заметим, что при изменении знака по тем же причинам будет наблюдаться уменьшение потерь ВМ с ростом амплитуды поля.

Поведение коэффициента распространения ВМ от Aq0 существенно меняется в узкой области параметров, где возникает корень дисперсионного уравнения, соответствующий АПМ. Подобный случай иллюстрирует рис. 6, на котором представлены действительная и мнимая части величины = /(kng ). Кривые для этого рисунка построены при = 0.1, k 2 dd1 /2 = 1 и 0 = 0.3. В такой области параметров мода ТЕq0 существует не при всех Aq0. При увеличении Aq0 возникают две АПМ, у которых Im(a0 ) 0, то есть эти моды не имеют радиационных потерь.

Таким образом, при увеличении Aq0 величина Im(q0 ) будет уменьшаться и в точке изменения типа моды q0 станет чисто действительной величиной. Точка отсечки, в которой величина Im(q0 ) обратится в ноль, показана на рис. 6 темным кружком.

Следует учесть, что поля АПМ растут при y ± и их нельзя возбудить источниРис. 5. Зависимости безразмерных коэффициен- Рис. 6. Зависимости безразмерных коэффициентов распространения и затухания вытекающей тов распространения и затухания вытекающей моды низшего номера ТЕq0 от амплитуды поля моды низшего номера ТЕq0 от амплитуды поля на оси волновода Aq0 при 0 = 0.6 на оси волновода Aq0 при 0 = 0.3 Для нелинейной задачи величина q0 зависит от z, поэтому эту величину нельзя называть постоянной распространения, как это принято для линейных систем.

ком с конечной мощностью. Поэтому можно предположить, что в указанной выше точке, в которой корень дисперсионного уравнения переходит из одного множества (класса) в другое, модовое распространение излучения (на симметричных модах) будет практически невозможно. Делая этот вывод, следует, конечно, учесть, что в такой ситуации есть возможность возбуждения ВМ высших номеров, которые, как правило, имеют большие потери, а также возможность возбуждения пространственной волны. Детальный анализ этого эффекта должен включать исследование задачи возбуждения и распространения мод вдоль оси волновода, что выходит за рамки данного исследования. Отметим также, что в рассматриваемой системе существуют другие точки, в которых возможен переход моды из одного класса в другой. Например, в окрестности точки 0 = 0 направляемая мода может трансформироваться в антиповерхностную.

Заключение

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

Расчеты показали, что в такой системе при варьировании ее параметров возможны эффекты преобразования типов мод. Под таким преобразованием понимается эффект изменения основных характеристик мод, по которым они классифицируются, например, изменение знака мнимой части поперечного волнового числа. Можно сказать, что в этом случае мода переходит из одного класса (множества), к которому она вначале принадлежала, в другой класс. Подобный переход может происходить также при изменении передаваемой по волноводу мощности.

При этом вытекающая мода может «превращаться» в антиповерхностную. Поскольку антиповерхностные моды не могут возбуждаться реальными источниками, то в области параметров, где корень дисперсионного уравнения соответствует именно такой моде, возможна ситуация, когда в данной структуре модовое распространение волн будет сильно подавлено. В линейной задаче при таких параметрах основная доля излучения уйдет в непрерывный спектр (пространственную волну) и, возможно, в незначительной доле в какие-то моды высших типов. На этом эффекте могут быть сконструированы датчики или элементы схем интегральной оптики.

В работе рассмотрены характеристики мод в волноводе с тонкими диэлектрическими стенками. Описанный подход может быть обобщен на случай, когда стенка волновода образована металлическими решетками (с нелинейной проводимостью) или на случай, когда волновод представляет собой слоистую плазменную структуру.

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

Автор признателен И.А. Молоткову и А.Г. Рожневу за обсуждения вопросов, связанных с данной работой.

Работа выполнена при поддержке РФФИ (гранты 08-02-00621 и 08-02-90002).

Библиографический список

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., Winn J.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.

22. Ньютон Р. Теория рассеяния волн и частиц. М.: Мир, 1969.

23. Каханер Д., Моулер К., Нэш С. Численные методы и программное обеспечение.

М.: Мир, 1998.

–  –  –

Characteristics of modes of the planar waveguide with thin walls, which are made from nonlinear media, are studied. The eect of the conversion of a mode from one set to another is considered. This eect may occur if parameters of the waveguide or transmitted power are varied. It is shown that the eect of the guided or leaky modes «disappearance» can arise when the eld 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.

Похожие работы:

«Грипп (фр. grippe, от нем. grippen — "схватить", "резко сжать") — острое инфекционное заболевание дыхательных путей, вызываемое вирусом гриппа. Входит в группу острых респираторных вирусных инфекций (ОРВИ). Периодически распространяется в виде эпидемий и пандемий. В настоящее время выявлено более 2000 вар...»

«ТАДЖИКИСТАН Пытки и жестокое обращение, безнаказанность, ущемление свободы мнения, насилие над женщинами и сохранение смертной казни Представление Amnesty International в рамках процедуры универсального периодического обзора ООН. Октябрь 2011 года СОДЕРЖАНИЕ ВВЕДЕНИЕ НОРМАТИВНАЯ БАЗА И ВЕДОМСТВЕННЫЕ СТРУКТУРЫ ГОСУДАРСТВА....»

«ПАРАЗИТОЛОГИЯ, 21, 3, 1987 УДК 576.895.421 : 595.421 (574 + 575) НОВЫЙ ВИД ИКСОДОВОГО КЛЕЩА DERMACENTOR USHAKOVAE 1 SP. N. (IXODOIDEA, IXODIDAE) ИЗ КАЗАХСТАНА И СРЕДНЕЙ АЗИИ Н. А. Филиппова, И. В. Панова Новый вид описывается по самке и самц...»

«ПРЕДИСЛОВИЯ К ПЕРВОМУ РУССКОМУ ПЕРЕВОДУ ЮЛИАНЫ НОРИДЖСКОЙ КАЛЛИСТ (УЭР), митрополит Диоклийский Представить читателю первый русский перевод "Откровений Божественной Любви" Юлианы Нориджской – большая радость для меня. Книга эта исключительна по своей мудрости и духовной проницательности. Юлиана исследуе...»

«УДК 621.77 КИРИЕНКО А. А., аспирант, КТУ, г. Кривой Рог РАСПРЕДЕЛЕНИЕ ЭНЕРГИИ ПРИ ПРОКАТКЕ Стаття присвячена розподілу енергії при прокатці. А саме енергії зміни форми і об’єму. У с...»

«Исследования жизни и творчества писателя за 1999 – 2003 годы из фонда Музея "М. А. Шолохов" Год издания 1999 Военная проза : [Антология : В 3 т.] / Ин-т “Открытое о-во”; [Сост. и коммент. Л. И. Лазарева]. – М. : Слово, 1999. – (Пушкинская библиотека / Редкол.: М. Л. Гаспаров и др.).– На авантит.: Мегапроект “Пушкин. б-ка”. [Т. 1] : Виктор...»

«Annotation Юная принцесса Майрам волей богов и ради хрупкого мира для своего народа отдана в жены жестокому варвару. Но быть второй женой слишком унизительно. Не об этом мечтала принцесса, выросшая в почете и уважении. Страшась своей судьбы и нежела...»

«A/HRC/WG.6/23/NRU/3 Организация Объединенных Наций Генеральная Ассамблея Distr.: General 10 August 2015 Russian Original: English Совет по правам человека Рабочая группа по универсальному периодическому обзору Двадцать третья сессия 2–13 ноября 2015 года Резюме, подготовленное Управлени...»

«"Crede Experto: транспорт, общество, образование, язык" — международный информационно-аналитический журнал №2 (09). Сентябрь 2014 (http://ce.if-mstuca.ru/) УДК 808.2(077) ББК Ш141.12-916.3 С544...»








 
2017 www.kniga.lib-i.ru - «Бесплатная электронная библиотека - онлайн материалы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.