Дифракции волноводных мод на стыке планарных волноводов Текст научной статьи по специальности «Физика»

УДК 501

Кочанова Мария Андреевна

ООО «ТЭМ консалтинг» Россия, Москва1 Консультант E-Mail: [email protected]

Дифракции волноводных мод на стыке планарных волноводов

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

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

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

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

На начальном этапе из рассмотрения исключаются излучательные моды.

Ключевые слова: математическое моделирование; численное моделирование;

волноводы; направляемые волноводные моды; дифракция волноводных мод; волновая оптика; оптика наноструктур; проектирование оптических устройств; волновые процессы; электродинамика.

Идентификационный номер статьи в журнале 52TVN314

1 117923 г. Москва, ул. Орджоникидзе, д. 3

Введение

Электромагнитные волны оптического диапазона издавна изучались в рамках геометрической и волновой оптики некогерентного и когерентного излучения. Сегодня оптические каналы передачи, приема и обработки информации широко применяются в современных системах связи и телекоммуникаций, включая всемирную сеть Интернет. Для этих целей используют оптические волноводы различного типа, среди которых преобладают волоконные и планарные конструкции как основные волноведущие структуры. Обработка информации в оптических волноводах осуществляется на основе материальных сред со специфическими свойствами, позволяющими эффективно управлять характеристиками распространения электромагнитных волн. К подобным средам, прежде всего, относятся материалы, обладающие акустооптическими, электрооптическими и магнитооптическими свойствами. Широкое применение также находят периодические волноводы с гофрированной поверхностью и связанные диэлектрические волноводы.

Электродинамические проблемы волновых процессов в разнообразных устройствах волоконной и интегральной оптики составляют сегодня предмет волноводной оптики. Теоретическое изучение электромагнитных волн в любой волноведущей структуре всегда базируется на совместном решении уравнений Максвелла с учетом специфики применяемой среды в форме материальных уравнений и электродинамических граничных условий. Анализ любой волноведущей структуры имеет конечной целью выяснение характера волновых взаимодействий в форме обмена мощностью между модами этой структуры при наличии ее возмущения.

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

Для того чтобы проанализировать взаимодействие направляемых мод дискретного спектра, а впоследствии и излучательных мод непрерывного спектра, в качестве базового волновода выбрана планарная трехслойная диэлектрическая структура. В этом случае можно строго вывести явные выражения для ТЕ- и ТМ-поляризованных мод и последовательно

выполнить анализ скоростей и структуры полей для направляемых (и излучательных)

мод.

Актуальность выбранной темы для исследований обусловлена тем фактом, что в последнее время возрос интерес к изучению и проектированию волноведущих систем, включающих открытые волноводы. Для проектирования необходим адекватный учет амплитудно-фазовых согласований (рассогласований) на границах раздела сред и устройств. Одним из таких сочленений интегрально-оптических волноводов является их жесткий стык. А недостаточно изученными являются процессы отражения и преломления волноводного излучения в терминах базисных систем волноводных мод: волноводные законы Снеллиуса и волноводные формулы Френеля. Их изучению и выводу посвящена данная работа.

Целью работы является получение системы линейных алгебраических уравнений для определения вклада каждой направляемой моды в описании дифракции волноводных мод на стыке планарных волноводов. На начальном этапе из рассмотрения исключаются излучательные моды.

Кратко охарактеризуем задачу, которая решается в данной работе.

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

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

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

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

Математическая модель

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

Ограничим класс рассматриваемых электромагнитных волн плоскими, однородными, монохроматическими волнами, зависимость полевых векторов которых от пространственных

координат г и времени / определяется формулами:

Е{г,1) = £0ехр[-/(ю/-£г)]; (1)

Н(г,Х) = Н0ещ>[Ч(Ш-кг)\; (2)

/)(г,0 = /)0ехр[-/(ю/-£г)]; (3)

В{г,1) = 50 ехр[-/(ю/ -кг)]; (4)

где со - циклическая частота, г - радиус-вектор точки наблюдения, к - волновой вектор, Е - вектор напряженности электрического поля, О - вектор индукции электрического ПОЛЯ, Н -вектор напряженности магнитного поля, В - вектор индукции магнитного поля.

Теперь опишем процесс распространения электромагнитной волны в среде из плоскопараллельных диэлектрических подобластей. В дальнейшем под термином «волновод»

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

Процесс распространения электромагнитной волны в такой среде описывается уравнениями Максвелла. Запишем их в системе единиц Гаусса:

гоШ = ——; (5)

С дt

го{е = -——\ (6)

С дt

сИл’Б = 0; (7)

(Ну В = 0. 8)

Уравнение (5) показывает, что магнитное поле порождается переменным электрическим полем, уравнение (6) представляет собой математическую формулировку закона электромагнитной индукции, уравнение (7) выражает факт отсутствия статического

электрического поля, уравнение (8) постулирует отсутствие в природе магнитных зарядов.

Для совместности системы (5)-(8), дополним уравнения электромагнитного поля материальными уравнениями, в которых учитывается соотношение между векторами £> и Е, В и Н . Для изотропной среды связь между индукцией и напряженностью электромагнитных волн определяется скалярами в и ^. Таким образом, материальные уравнения приобретают

В = еЕ\ (9)

В = /иН. 10)

Так как мы рассматриваем только не обладающие магнитными свойствами среды, связь между векторами напряженности и индукции магнитного поля определяется следующим образом:

В = Н. (11)

Следует также сформулировать граничные условия для уравнений электромагнитного поля, которые, в случае нашей задачи, сводятся к равенству тангенциальных компонент

векторов Н и Е на плоской границе раздела двух планарных волноводов:

к=к> <12)

Также добавим асимптотические граничные условия на бесконечности, так как в

рассматриваемой системе присутствуют полубесконечные среды:

| Ё\\ —» 0; || Я || -» 0. (13)

Итак, система уравнений, в которую входят уравнения электромагнитного поля (5)-(8), материальные соотношения (9)-(10) и граничные условия (12)-(13), полностью описывает процесс распространения электромагнитного поля в выбранной волноведущей структуре.

Физическая постановка задачи

Трехслойный планарный волновод, представленный на рисунке ниже, состоит из трех слоев, но образован четырьмя подобластями:

• полубесконечный покровный слой (воздух) с показателем преломления пс;

• волноводный слой (пленка), состоящий из двух частей с показателями преломления п^ и п^ , плотно пристыкованных друг к другу перпендикулярно к границам раздела слоев;

• полубесконечная подложка с показателем преломления пв;

На рисунке <Л = а2 — ах - толщина волноводного слоя.

Все слои однородны в направлениях х и г.

Рисунок - Схематическое изображение стыка двух полубесконечныхрегулярных

вдоль оси Oz планарных волноводов

Опишем явление дифракции волноводных мод на стыке планарных волноводов волноведущей структуры.

Показатели преломления слоев пу, пу, пу; у = 1,2 вещественны. Толщина ё допускает существование N ТЕ- и ТМ-поляризованных волноводных мод в I (левой) подобласти и Ы2 ТЕ- и ТМ-поляризованных волноводных мод во II (правой) подобласти, задаваемых компонентами Еу и Ну соответственно: Е\ м ;Н\ м и Н11!1^г.

Слева на границу раздела набегает выделенная мода, например ТЕ\ , описываемая стоячей вдоль оси Ох волной ЛЕ!2 (х) . После дифракции на стыке волноводов влево отразятся моды Я^е!(х) и ЩМН!(х); у = 1,...Ж1 . Во второй подобласти (в правой части волновода)

будут распространяться моды Т^Е1](х) и Т^И1](х); ] = 1,...,Щ . Здесь А - амплитуда

7 7 4 7 7 7

ТЕ грТ V ’ ~ V ’ ^ 7 , І

падающей волны, ЯТ’Е, Я™, ТТЕ, Т™ -амплитуды отраженных и прошедших волн

соответственно.

Таким образом, поле в подобласти I - это совокупность полей падающих и отраженных мод. Поле в подобласти 1 будет иметь вид:

Е‘ = ЕП (х,2, і) + Е? (х, 2, і) (14)

Слагаемое Еус (х, z) описывает падающие (набегающие на границу раздела подобластей) волны, слагаемое Еуе (х, г) описывает волны, отраженные от границы раздела подобластей.

С учетом граничных условий на бесконечности поля в подобласти I (слева) примут вид:

Щі

Е!у = АЕ[ (х) ехр{/(юі - к0 рТ1Е2)} + Ъ ЯЕІ (х) ехр{/(юі + к0 рт}Е г)};

7=1

щ (15)

К = Ъ ЯИ (х)ЄХр{і(®і + к„р™ г)}.

7=1

Здесь 3 - коэффициент фазового замедления.

Поле в подобласти II будет представлять собой совокупность прошедших мод:

е; = Еуг (х, 2, і) (16)

С учетом граничных условий на бесконечности поля в подобласти II (справа) примут

вид:

ЕУ = 2 ТЕ! ( х) ехр{^—ко^.Гг)};

у=1

N (17>

Н? = 2 Т.Н- ( х) ехрМ®( — ко,р™ г )}.

У =1

На границе раздела подобластей I и II должно выполняться равенство тангенциальных компонент.

Для ТЕ-мод:

ЕП ( х, г, г) + е;Г ( х, г, г) = Е“” ( х, г, г) (18)

Для ТМ-мод:

Н1? ( х, г, г) = Н’~(х,г,г) (19)

На границе раздела подобластей тангенциальными компонентами также являются компоненты Нх для ТЕ-мод и Ех для ТМ-мод, которые выражаются через Е и Н с помощью уравнений Максвелла.

Выразим их через компоненты Е и Н . Используя уравнения Максвелла, получим:

£ дг

(20)

Согласно граничным условиям, справедливы следующие равенства:

ИП (х, I) + И? (х, г, I) = И'Г (х, г, I)

Б? (х, г, I) = Е'Г (х, г, I)

(21)

(22)

На границе раздела подобластей 2=0, на оси Ох тангенциальные компоненты электромагнитного поля Б и Их для ТЕ-мод и И и Ех для ТЕ-мод совпадают. В силу того,

что 2=0, множители ехр(±/£0Дг) = 1. Равенство тангенциальных компонент выполняется в

любой момент времени, а значит и при I = 0, тогда множитель ехр(/®1) = 1.

Итоговая система уравнений, согласно граничным условиям, примет вид:

Наша задача - найти неопределенные амплитудные коэффициенты этой системы уравнений и таким образом определить вклад каждой отраженной или прошедшей моды

Решение задачи

Чтобы найти неопределенные амплитудные коэффициенты этой системы уравнений и, таким образом, определить вклад каждой отраженной или прошедшей моды, проделаем следующую процедуру. Спроектируем отрезки функциональных рядов на подпространства биортогонального базиса. Более детально - умножим каждое уравнение системы на комплексно-сопряженные функции к функциям Е1 (х); Е^ (х); Н1 (х); И^ (х) , а затем

проинтегрируем полученные выражения по всей оси Ох. Затем воспользуемся условием ортогональности:

(23)

]=1

]=1

Умножив систему уравнений сначала на Е*1 (х) и И *! (х) (комплексно-сопряженные к

функциям в левой части системы) и проинтегрировав по х, получим следующую систему (для удобства запишем ее в матричном виде):

^АЫ1х2Ы1 '^2^x1 — '^т,х1 _ х1 (25)

где А и В - блочные матрицы с числовыми элементами, Я и Т - , векторы неопределенных коэффициентов, С - вектор свободных членов.

Умножив систему на Е*11 (х) и И*п (х) (комплексно-сопряженные к функциям в левой части системы) и проинтегрировав по х, получим:

/V —► /V

х2Щ ' ^2Иг х1 — ^4Ы2 x2N2 ' х1 _ ^ЧЛ?2 х1 ’ (26)

где I и С - блочные матрицы с числовыми элементами.

Итоговая система в матричном виде :

^4(Л?1+Л?2)х2Л?1 '^2^x1 — С?4(Л?1+Л?2)х2№ '^2Л?2х1 _ ^~'4(.Ы1+Ы2 )х1 >

(27)

где Р и Q - блочные матрицы, составленные из А и I, В и С соответственно. Учитывая (24), примем следующие обозначения:

да

(Е‘, Е‘‘) =| Е‘ (х)Е]п (х)Л;

—да

(Е1, Ек ) = 81к; (28)

Е, ] = 8,.

да

И, И‘‘) =| И (х)И'/ (х)Л;

—да

И, И[) = 5к; (29)

И, ] = 8,.

Используя (28) и (29), представим матрицы, фигурирующие в системе, более подробно.

Матрица Р в общем виде:

ґ ' Е1 Е1

Еі, еі

Р

~(Е1, ЕІ) ц' '

Е1 Е1)

ЕЫ , Е1 /

И

Е1, ЕІ

0

0

Р

0

(н{ ,Н) ... (НІ.НІ)

Р / тт 1 тт11 \Н N ,Н ы2 є \ 1 2 , у

(30)

Матрица Q в общем виде:

е=

Е£, Е1

.Л/ е11 Е

\ЕМ, , Е№,

И '

0

0

0

Н1, Н) ... (нНN ,Н/)

0

0

Р/и11 и11 \ Р / и11 и11 \

\И1 , - \HN2, HN2)

£ £

Для примера покажем подробнее получение одного из элементов матриц, например,

да

{е;,Е‘2')=1 Е/(х)Е-"(х) :

—да

да а; а2

Е, Е11) = | Е[ (х)Е*2И (х) = | Е;1 (х)Е;ш (х^ +1Е? (х)Е\П1 (х^х +

—да —да a^

да а; да

+| Е1(х) Е2с11 (х)^х = 144 ехр((^ + ^2) х№ +144* ехр(—(гс + г2) х¥х +

а

а2

+| {Аг + ехр(/^?х) + Аг — ехр(—/*?х)} {Аг *+ ехр(/^?х) + А г *— ехр(—/^х)} dx

а1

Остальные элементы матриц Р и Q вычисляются аналогично.

(31)

Выводы

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

ЛИТЕРАТУРА

1. Барыбин A.A. Электродинамика волноведущих структур. - М.: ФИЗМAТЛИТ, 2007.

2. Кочанова МА., Севастьянов ЛА. Дифракции волноводных мод на стыке планарных волноводов // Материалы всероссийской конференции с международным участием «Информационно-телекоммуникационные технологии и математическое моделирование». - Москва, РУДН, 22-25 апреля 2014. - С. 248-251.

3. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. - М.: Наука, 1982. - Т. 8.

4. Калитеевский Н.И. Волновая оптика. - М.: Высшая школа, 1978.

5. Лхманов СА., Никитин С.Ю. Физическая оптика. - М.: Наука, 2004.

6. Хансперджер Р. Интегральная оптика. - М.: Мир, 1985.

7. Терлецкий ЯП., Рыбаков Ю.П. Электродинамика. - М.: Высшая школа, 1990.

S. Никоноров Н.В., Шандаров С.М. Волноводная фотоника. - СПб.: СПбГУ ИТМО,

200S.

9. Ayrjan E.A., Egorov A.A., Michuk E.N., Sevastyanov A.L., Sevastianov L.A., Stavtsev A.V. Representations of guided modes of integrated-optical multilayer thin-film waveguides // Preprint JINR E11-2011-31, Dubna, 2011, 52 P.

10. Oulton R. F., Pile D. F. P., Liu Y., Zhang X. Scattering of surface plasmon polaritons at abrupt surface interfaces:Implications for nanoscale cavities // Phys. Rev. B 76, 03540S (2007).

11. Mahmoud S. F., Beal J. C.. Scattering of Surface Waves at a Dielectric Discontinuity on a Planar Waveguide // IEEE Trans. Microwave Theory Tech. MTT-23, 193 1975

Рецензент: Севастьянов ЛА., д.ф-м. н., профессор ФГБОУ ВПО РУДН.

Maria Kochanova

TEM consulting Russia, Moscow E-Mail: [email protected]

Waveguide modes diffraction at a dielectric discontinuity on a planar waveguide

Abstract: This work deals with the problem of waveguide mode diffraction at a dielectric discontinuity of two half-infinite planar waveguides that are regular on each side of the waveguide system given. Physically a wave guide is a channel that allows electromagnetic wave spread along an axial line with relative low fading effect and localizes it within the aria close to the axial line.

The ability of artificial waveguides to transmit energy or information (signals) of different kind is based on the total reflection effect.

There can be two types of waveguide modes in open planar waveguides: guided modes, which form discontinuous spectrum, and radiation modes, which form continuous spectrum. When waveguide modes diffract at a dielectric discontinuity of a planar waveguide, the effect is similar to refraction and reflection of plane body waves at a bedding interface.

The aim of this work is derivation of simultaneous linear algebraic equations in order to evaluate the contribution of each mode to the resulting electromagnetic field while considering waveguide mode diffraction at a dielectric discontinuity of half-infinite planar waveguides. The main research methods are: analytical reduction of Maxwell equations considering boundary conditions for electromagnetic field and the application of boundary conditions of Maxwell equations diffraction at a discontinuity of planar waveguides for expansion in orthogonal modes in two parts of the waveguide system using biorthogonal basis technique.

At the beginning radiation modes are excluded from the consideration.

Keywords: Mathematical modeling; numerical modeling; waveguides; guided waveguide modes; waveguide mode diffraction; wave optics; nanostructural optics; optical devices design; wave process; electrodynamics.

Identification number of article 52TVN314

REFERENCES

1.

2.

3.

4.

5.

6.

7.

8.

9.

10. 11.

Barybin A.A. Jelektrodinamika volnovedushhih struktur. - M.: FIZMATLIT, 2007.

Kochanova M.A., Sevast'janov L.A. Difrakcii volnovodnyh mod na styke planarnyh volnovodov // Materialy vserossijskoj konferencii s mezhdunarodnym uchastiem «Informacionno-telekommunikacionnye tehnologii i matematicheskoe modelirovanie». - Moskva, RUDN, 22-25 aprelja 2014. - S. 248-251.

Landau L.D., Lifshic E.M. Teoreticheskaja fizika. - M.: Nauka, 1982. - T. 8.

Kaliteevskij N.I. Volnovaja optika. - M.: Vysshaja shkola, 1978.

Ahmanov S.A., Nikitin S.Ju. Fizicheskaja optika. - M.: Nauka, 2004.

Hansperdzher R. Integral'naja optika. - M.: Mir, 1985.

Terleckij Ja.P., Rybakov Ju.P. Jelektrodinamika. - M.: Vysshaja shkola, 1990.

Nikonorov N.V., Shandarov S.M. Volnovodnaja fotonika. - SPb.: SPbGU ITMO, 2008.

Ayrjan E.A., Egorov A.A., Michuk E.N., Sevastyanov A.L., Sevastianov L.A., Stavtsev A.V. Representations of guided modes of integrated-optical multilayer thin-film waveguides // Preprint JINR E11-2011-31, Dubna, 2011, 52 P.

Oulton R. F., Pile D. F. P., Liu Y., Zhang X. Scattering of surface plasmon polaritons at abrupt surface interfaces:Implications for nanoscale cavities // Phys. Rev. B 76, 035408 (2007).

Mahmoud S. F., Beal J. C.. Scattering of Surface Waves at a Dielectric Discontinuity on a Planar Waveguide // IEEE Trans. Microwave Theory Tech. MTT-23, 193 1975.