CC BY
74
РАСЧЕТ СОБСТВЕННЫХ ВОЛН ПЛАНАРНОГО АНИЗОТРОПНОГО ВОЛНОВОДА ДЛЯ РАЗЛИЧНЫХ ПОЛОЖЕНИЙ ОПТИЧЕСКОЙ ОСИ
Моисеева Н.М.
Волгоградский государственный университет
Аннотация
Для однородного планарного анизотропного волновода выполнено решение уравнений Максвелла. Получены фундаментальные матрицы решения для ТЕ- и ТМ-волн в анизотропном плоском однородном слое. Показано, что фазовый сдвиг ТМ-волны при отражении от границы «анизотропный слой - изотропная среда» зависит от угла наклона оптической оси. Выполнено численное решение дисперсионного уравнения ТМ-волны в анизотропном планарном волноводе при различных ориентациях оптической оси в плоскости распространения волны; найдена асимптотика дисперсионных кривых.
Ключевые слова: планарный волновод, анизотропия, уравнения Максвелла, дисперсионное уравнение, полное внутреннее отражение, фазовый сдвиг при отражении света, асимптотика дисперсионных кривых.
Введение
В настоящее время происходит интенсивное развитие технологий оптоэлектроники и интегральной оптики, значительно вырос интерес исследователей к волноведущим оптическим структурам [1]. Это можно объяснить развитием оптических технологий обработки и передачи информации, разработкой фотонных интегральных схем, оптической памяти [17]. Преимуществами оптических волноведущих структур являются высокая скорость передачи информации, низкая стоимость сырья, компактные размеры. Планарные диэлектрические волноводы являются важнейшими элементами новых высокоскоростных технологий. Толщина волновода и его оптические параметры определяют особенности распространения и отражения в нём электромагнитных волн, то есть его собственные моды. Известно, что в анизотропных средах электромагнитные волны распространяются иначе, чем в изотропных средах [2]. В анизотропной среде оптические свойства зависят от направления распространения света, и поэтому ориентация оптической оси волновода должна влиять на условия распространения собственных волн. Новые оптические решения требуют создания сред со специальными оптическими свойствами. Как известно, интегрально-оптические компоненты изготавливаются в очень сложных технологических процессах, например, в процессе ионной имплантации. Чтобы создаваемые устройства функционировали так, как это планирует разработчик, необходим детальный анализ волноводных характеристик распространения, а также разработка простого набора расчётных инструментов для производства. Для этого требуется адекватная математическая модель взаимодействия электромагнитных волн с веществом, построенная на уравнениях Максвелла.
За последние десятилетия разработано множество аналитических и численных методов расчёта волно-водных мод [3 - 4]. Метод распространения луча применён в работе [7] для анизотропного тонкого волновода. Рассматривалась трёхслойная система. В случае произвольной ориентации оптической оси анизотропного волновода в каждом слое будет происходить
кросс-поляризация, взаимодействие волн. Для описания этого явления в работе [8] Юхом была предложена новая 4^4 матричная алгебра, которая комбинирует метод матриц 2*2 и метод матриц Джонса. Она была применена для исследования распространения плоской волны в произвольной анизотропной среде. Для каждого слоя записывалась матрица распространения; было выяснено, что вид матрицы зависит от ориентации оптического тензора относительно осей волновода. Процедура, основанная на формализме матриц Юха 4*4, получила развитие в работе [9]. Было найдено асимптотическое решение для набора очень тонких слоёв, толщиной значительно меньше длины волны. Метод подбора нулевых элементов характеристической матрицы был предложен в статье [10]. Его применяли для многослойного планарного оптического волновода с заданным профилем показателя преломления. Постоянные распространения волновода были получены из условий равенства нулю элементов передаточной матрицы. Метод применялся к волноводам с потерями и анизотропией. В работе [11] был предложен метод передаточной матрицы. Это стандартный матричный метод 2*2, используемый для оптики тонких плёнок. Его применили к планарному многослойному волноводу. Получены выражения для расчёта профилей полей и энергии в каждом слое. Метод, применённый к поглощающему многослойному волноводу и к отражению плоской волны, был предложен в [12]; он является адаптированной для планарного волновода модификацией метода характеристических матриц [2]. В монографии [1] предложен детальный теоретический анализ различных видов оптических волноводов. Детально исследованы условия существования волн, получены дисперсионные уравнения для расчёта зависимости показателя преломления для первых двух мод от относительной толщины волновода, указан интервал для верхнего предела «изменения показателя преломления».
Для расчёта и анализа волноводных мод часто применяют модификации метода конечных разностей для временной и частотной области, а также метод конечных элементов. Оригинальный алгоритм был предложен Йи [13]. Новая полно-векторная конечно-
разностная дискретизация для волноводов с поперечной анизотропией была предложена в [14]. В отличие от предшествующих решений метод позволяет решать задачи для произвольной ориентации оптической оси относительно распространения волны. Метод конечных элементов применялся для расчёта мод при различных углах распространения в работах японских учёных [15] и [16]. Численные методы позволяют выполнить решения для собственных мод для волноводов с различными оптическими свойствами и геометрией. Однако их общими недостатками являются появление нефизических зависимостей, а также затруднённость в анализе решения.
1. Постановка задачи
В настоящей работе, исходя из уравнений классической электродинамики, выполнена сшивка решений на границах трёхслойного однородного волновода с внутренней анизотропной средой с помощью матрицы Коши. Выполнен расчёт решений дисперсионных уравнений для волн ТЕ и ТМ, распространяющихся в анизотропном волноводе, получена зависимость асимптотики дисперсионных кривых от угла наклона оптической оси волновода.
Рассмотрим распространение электромагнитной волны в анизотропном плоском слое. Он представлен на рис. 1. Покровная среда «1» и подложка «3» имеют меньшие значения диэлектрической проницаемости, чем значения еО, ее второй среды, и обеспечивают распространение электромагнитной волны в пределах слоя «2» в плоскости Х02 за счёт явления полного внутреннего отражения. Плоскость распространения волн Х02, показанную на рисунке, будем называть плоскостью падения. Пусть оптическая ось анизотропной среды волновода ориентирована в плоскости падения под некоторым углом ф к оси ОХ.
1 к 81
*Л
O'd 3 А ез
Рис. 1. Однородный анизотропный планарный волновод. 00' - оптическая ось анизотропной среды «2» расположена в плоскости XOZ под углом ф к оси OX
При указанном положении оптической оси диэлектрическая проницаемость одноосной анизотропной среды «2» имеет вид:
(
ee cos2 ф + eo sin2 ф 0 0 e
e -e ■ „, -o sin2f
2
e -e e — sin2f
2
0 eo cos2 ф + ee sin2 ф
(1)
Рассмотрим распространение в среде «2» волн двух типов. Первая волна будет обыкновенной волной в кристалле «2»; это ТЕ-волна. Векторы напря-
женностеи электрического и магнитного полей для неё имеют вид:
= (о
E,
HTE = (Hy
Ey о), (2а)
0 Hz). (2б)
Вторая волна в кристалле «2» является необыкновенной, это ТМ-волна. Вектор напряжённости электрического поля для данной волны лежит в плоскости падения и имеет компоненты ЕХ, Е2, а вектор Н перпендикулярен ей, поэтому имеет только компоненту Иу :
Е2). (3а)
0) . (3б)
ETM - (Ex 0
HT
-(о
H
Подставим векторы (2) и (3) в уравнения Максвелла и учтём, что для всех компонент полей Е, выполняются условия:
^ Е = 0, (4а)
- ПсоЩ .
(4б)
Параметр а зависит от угла падения волны на границы раздела х = 0 и х = ё и от оптических свойств волновода. Значения углов падения 81, 83 в средах «1» и «3», а также значения углов 8ТЕ и 8тм в слое «2» связаны с оптическими свойствами материалов законом Снеллиуса:
а = .у/ё11!5т 8! = ^/е3ц3 83 =
■4
(5)
--y/eyymsin8te -n sin8tm - const.
Здесь ne - cos2 (8 - ф) + ee sin2 (8 - ф) - показатель преломления для волны TM в слое «2» [2, с. 627].
2. Решение уравнений Максвелла Подстановка векторов (2) и (3) в уравнения Максвелла с учётом связей (4) и (5) позволяет получить две системы, состоящие из двух обыкновенных дифференциальных уравнений. Первая система ОДУ описывает распространение в среде «2» TE-волны: d
—E - ikomHz dx
~rHz - iko dx
e yy -
2A
(6)
Ey.
Вторая система из двух ОДУ отвечает за распространение в среде «2» волны ТМ:
d е
—H --iko a—H - iko
dx y e y
e „ --
E
~TEz --iko dx
' a2 ^
m—
e
(7)
H -ikoa^Ez e
Для решения систем ОДУ (6) и (7) были найдены фундаментальные матрицы решения, а затем получены матрицы Коши для ТЕ- и ТМ-волн в анизотропном слое. Данный метод применялся в работах
e
o
e2 -
5, 6, 18]. В работах [19, 20] были получены матричные решения для однородного анизотропного слоя. «Сшивка» решений на границах «1 - 2», когда х = 0, и «2 - 3», когда х = й, для компонент полей ТЕ-волны имеет вид:
Г Еу (й) ^ ( > пТ>V Еу г Л
VН (0) ,
Коэффициенты матрицы Коши в уравнении (8) имеют вид:
) —(ТЕ) - со8(к0^еууц-а2),
Нг (d))
n(TE ) n(TE ) V 21 "22 )
(8)
= n
22
(te ) im sin(^Q dj e ^ m-a2),
e yy m-a
JTE) = '21 _
^e yym-a2
m
(8а) (8б)
sln(k0d^eyym-a2). (8в)
Матричная «сшивка» решения для ТМ-волны на границах слоя «2» записывается в виде:
r Hy (d) v Ez (d)
Г n (TM)
nii
n (TM)
V n21
ATM) ^
12 n (22)
Hy (0) Ez (0)
(9)
Матрица Коши Nр (й,0) в формуле (8) для волны ТМ имеет вид:
e (
-ik0ad
(TM) (TM) 0 е.,
nn — n22 — e cos
k0 d
^A(e jm-a2)
nTM) =
'12
A
e ..m-a
v
k0 d
^A(e ..m-a2)
-ik ad
n2[M) = ie 0 e"
ex.m-a
A
k0 d
•x/A(exIm-a2)
(9а)
(9б)
(9в)
Величина А - еххе** - ее^ - е0ее не зависит от ориентации оптической оси в плоскости Х02. Формулы (8) и (9) позволяют рассчитывать локализацию полей электромагнитной волны в волноводе. Из формул (9а) - (9в) видно, что собственные значения системы ОДУ (5) принимают значения:
^1,2 = ±lk0
^A(e x,m-a2)
- гак0-^. (10)
е е
хх хх
Второе слагаемое в формуле (9) объясняет появление фазового множителя в коэффициентах матрицы (9). Сшивки (8) и (9) получены путём решения систем ОДУ при помощи матрицанта. Как правило, при расчётах применяют метод передаточной матрицы [11], в котором фазовый множитель опускают. Исходя из уравнений Максвелла, связь между компонентами напряжённости электрического и магнитного полей ТМ-волны записывается в виде системы:
| к2Ну - к0 (еххЕх + еI (11)
\кхНу — к0 (е хЕх +е ггЕг ) .
Из формул (10) и (11) следует, что связь между компонентами Ну и Е имеет вид:
Hy =±
A
(12)
je x,m-a2
Угол падения электромагнитной волны TM в слое «2» равен 8TM. Если ETM - вектор напряжённости электрического поля, то ETM cos 0TM - его проекция на ось 0Z. Теперь нам известны тангенциальные составляющие вектора H для падающей и отражённой волн:
Hy = ±
A
2 ETM cos 0TM ■
(13)
,е ххЦ-а2
«Сшивкам» граничных условий на плоскости х = 0 для волны, распространяющейся из среды «2» в среду «1» и отражённой в слой «2», позволяет получить амплитудный коэффициент отражения для ТМ-волны:
rp =
■y/e1 (e ..m-a2) - VXcos01
(14)
>/е1 (е ххЦ-а2) + л/А^! В слое «2» будет происходить полное внутреннее отражение при углах е™, больших угла е , определяемого по формуле:
(
0 = arcsin
^ sln2(0 -f) + ee cos2(0 -f)
(15)
Исходя из формулы (14), по методу, показанному в работе [2, с. 65], получим выражение для сдвига фазы при отражении от границы «2 - 1»:
f
S21 = -2arctg
Va e,
2
a2 -e
e..m-a
(16)
Аналогичная формула справедлива для границы
«2 - 3».
(
d23 = -2arctg
Va
e 3
2
a2 -e
Л
e „m-a
(17)
3. Решение дисперсионных уравнений
Величина А при любом значении угла ф остаётся неизменной, а значение ехх является функцией этого угла.
Для того чтобы волна могла распространяться в волноводе, необходимо, чтобы при её распространении от границы «1 - 2» к границе «2 - 3», с учётом сдвига фаз при отражении, от покровной среды и от подложки изменение фазы было кратно я. Это условие определяет вид дисперсионного уравнения для мод ТМ планарного волновода:
7А(е„Ц-а2) + Гл/а I а2-е --агйя--„I-1
е,
k0d
(
-arctg
a -e,
I ejm-a
exxm -a
= mp, m e N.
(18)
n
e
e
e
e
e
Уравнение, аналогичное (18), было получено в работе [12], однако в этой работе задача решалась для диагонального тензора €2. Величина ехх несёт
информацию об угле ф. Следовательно, несмотря на универсальный характер уравнения (18), его решения для анизотропного волновода должны определяться положением оптической оси.
Дисперсионное уравнение для ТЕ-волны будет совпадать с дисперсионным уравнением ТЕ-волн в изотропной среде с диэлектрической проницаемостью е2 = еуу :
(19)
Между параметром а и постоянной распространения Р существует связь: Р=£0а. На рис. 2 представлены дисперсионные кривые для волны ТЕ для мод т = 0, 1, 2, 3, 4, являющиеся решениями уравнения (19). Случай распространения ТЕ-волны в анизотропной среде при заданной выше ориентации оптической оси аналогичен распространению волны ТЕ типа в изотропной среде, как это и следует из уравнения (19).
' р/ы
0 10 20 30
Рис. 2. Зависимость постоянной распространения в от фазы к0ё волны ТЕ в анизотропном планарном волноводе для мод ТЕ0 - ТЕ4
Решение уравнения (18) для ТМ-волн выполнялось численно для искусственной среды с параметрами тензора диэлектрической проницаемости е0 = 3, ее = 4. Значения угла ф между оптической осью кристалла и нормалью, проведённой к границам раздела сред (рис. 1), изменялись в пределах [0°; 90°].
Расчёт показал, что для ТМ-волны положение дисперсионных кривых зависит от угла ф. Горизонтальная асимптота у = А ограничивает графики, составленные из точек, являющихся решениями дисперсионного уравнения, сверху. Значение А зависит от угла ф и связано с показателем преломления необыкновенной волны. Если угол ф = 0, высота горизонтальной линии
будет равна величине: А = д/ё7. Графики, рассчитанные при значении угла ф = 0, когда оптическая ось совпадает с осью ОХ, представлены на рис. 3а.
б)
Рис. 3. Зависимость постоянной Р от фазы к04 в анизотропном планарном волноводе для мод ТМ0 - ТМ4: угол ф = 0° (а), угол ф = 90° (б)
Параметр А в этом случае будет максимально отличаться от асимптоты графиков ТЕ-волны, приведённых на рис. 2. Если угол ф=90°, положение асимптоты
будет таким же, как и в изотропном случае: А = ^/ёО
(рис. 3б). Это не противоречит результатам, представленным в монографии [1]. Из них следует, что верхний предел изменения параметра пв, равного а, «заключён в интервале между Щ и п3», иными словами, между
¡ёо и Vее
Для произвольного значения угла ф параметр А принимает значение, равное показателю преломления необыкновенной волны для максимально возможного угла (90°- ф), между волновым вектором и оптической осью:
А = ^
2 2 = д/еО бш ф + ее со8 ф .
(20)
Если угол наклона оси ОО' равен ф, а волна распространяется вдоль оси 02, что соответствует максимальному значению а, то угол между волновым вектором и оптической осью будет равен (90°- ф).
На рис. 4а, б представлены дисперсионные кривые для случаев, когда угол ф составляет 30° и 60°. Нижняя штриховая линия на рис. 2 - 4 соответствует наименьшему значению параметра а при угле полного внутреннего отражения 8 : аш1п = ^тахСе,, е3) .
При изменении в формуле (1) для тензора диэлектрической проницаемости среды «2» угла ф на «-ф» характер зависимостей не изменяется.
волны TM в анизотропном волноводе для мод TM0 - TM4 от фазы k0d при различных значениях угла ф: ф = 30° (а), ф = 60° (б)
Заключение
В данной работе была получена математическая модель для расчёта собственных мод планарного анизотропного волновода для произвольного наклона оптической оси в плоскости падения. Исследована асимптотика дисперсионных кривых для TM-волны в зависимости от угла ориентации оптической оси в плоскости распространения волны.
Литература
1. Гончаренко, А.М. Основы теории оптических волноводов / А.М. Гончаренко, В.А. Карпенко. - М.: Едито-риал УРСС, 2004. - 240 с.
2. Борн, М. Основы оптики / М. Борн, Э.Вольф - М.: Наука, 1973. - 720 с.
3. Котляр, В.В. Моды планарного градиентного гиперболического секансного волновода / В.В. Котляр, А.А. Ковалёв, Я.Р. Триандафилов, А.Г. Налимов // Компьютерная Оптика. - 2010.- Т. 34, № 2. - С. 146155. - ISSN 0134-2452.
4. Мусакаев, М.Р. Математические модели поляризационной модовой дисперсии высших порядков для кварцевого анизотропного оптического волновода / М.Р. Мусакаев, А.Х. Султанов // Компьютерная оптика - 2012. - Т. 36, № 2. - С. 165-171.
5. Моисеева, Н.М. Расчёт дисперсионных характеристик поверхностных поляритонов на границе неоднородных сред / Н.М. Моисеева, В.В. Яцышен // Физика волновых процессов и радиотехнические системы. -2004. - Т. 7, № 2. - С. 30-32.
6. Моисеева, Н.М. Расчёт дисперсионных характеристик поверхностных плазмонов на границе неоднородной анизотропной среды / Н.М. Моисеева, В.В. Яцышен //
Физика волновых процессов и радиотехнические системы. - 2005. - Т. 8, № 1. - С. 77-81. - ISSN 1810-3189.
7. Gia Russo, D.P. Wave propagation in anisotropic thin film optical waveguides / D.P. Gia Russo, J.H. Harris // J. Opt. Soc. Am. - 1973. - Vol. 63, N 2. - P. 138-145.
8. Yeh, P. Optics of anisotropic layered media: a new 4x4 matrix algebra / P. Yeh // Surface Science. - 1980. -Vol. 96, N 1-3. - P. 41-53.
9. Visnovsky, S. Magneto-optic effects in ultrathin structures at longitudinal and polar magnetizations / S. Visnovsky // Czech. J. Phys. - 1998. - Vol. 48, N 9 - P. 1083-1104. -(Online).
10. Walpita, L.M. Solutions for planar optical waveguide equations by selecting zero elements in a characteristic matrix / L.M. Walpita //J. Opt. Soc. Am. A. - 1985 -Vol. 2, N 4. - P. 595-602. - (Online).
11. Chilwell, J. Thin-films field-transfer matrix theory of planar multilayer waveguides and reflection from prism-loaded waveguide / J. Chilwell, I. Hodgkinson // J. Opt. Soc. Am. A. - 1984. - Vol. 1. - P. 742-753. - (Online).
12. Liao, W. Explicit analysis of anisotropic planar waveguides by the analytical transfer-matrix method / W. Liao, X. Chen ,Y. Chen, Y. Xia, Y. Chen // J. Opt. Soc. Am. A.
- 2004. - Vol. 21, N 11. - P. 2196-2204. - (Online).
13. Yee, K.S. Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media / K.S. Yee // IEEE Trans. Antennas and Propag. -1966. - Vol. 14, N 3. - P. 302-307.
14. Fallahkhair, A.B. Vector Finite Difference Modesolver for Anisotropic Dielectric Waveguide / A.B. Fallahkhair, K.S. Li, Th.E. Murphy // Journal of Lightwave Technol. -2008. - Vol. 26, N 11. - P. 1423-1431.
15. Koshiba, M. Finite-Element Solution of Anisotropic Waveguides with Arbitrary Tensor Permittivity / M. Koshiba, K. Hayata, M. Suzuki // Lightw. Technol. - 1986. -Vol. 4, N 2. - P. 121-126.
16. Tsuji, Y. Finite element beam propagation method for anisotropic optical waveguides / Y. Tsuji, M. Koshiba, N. Takimoto // J. Lightw. Technol. - 1999. - Vol. 17, N 4.
- P. 723-728.
17. Burmeister, E.F. Photonic integrated circuits optical buffer for packet-switchen networks / E.F. Burmaister, J.P. Mack, H.N. Poulsen, M.L. Masanovic, B. Stamenik, D.J. Blumental // Optics Express - 2009 - Vol. 17, N 8. - P. 66296635. ISSN: 1094-4087.
18. Моисеева, Н.М. Особенности дисперсионных характеристик поверхностных плазмонов на границе неоднородной среды / Н.М. Моисеева, В.В. Яцышен // Вестник Волгоградского государственного университета. Серия 9, выпуск 3. - 2003-2004. - Ч. 2. - С. 37-41.
19. Моисеева, Н.М. Зависимость главного и псевдо-брюстеровского углов от ориентации оптической оси плоской одноосной анизотропной среды / Н. М. Моисеева // Лазеры. Измерения. Информация - 2011. - С.-Пб., Издательство Политехнического университета, 2011. - Т. 3. - С. 430-438.
20. Моисеева, Н.М. Кросс-поляризация электромагнитной волны при отражении на границе анизотропной среды / Н. М. Моисеева // Лазеры. Измерения. Информация - 2011. - С.-Пб. - Издательство Политехнического университета, 2011. - Т. 2: - C. 48-61.
References
1. Goncharenko, A.M. Fundamentals of the theory of optical waveguides / A.M. Goncarenko, V.A. Karpenko -Moscow: Editorial URSS, 2004. 240 p. - (In Russian).
2. Born, M. Principles of optics / M. Born, E. Wolf. - Moscow: "Nauka" Publisher, 1973. - 720 p. - (In Russian).
3. Kotlyar, V.V. Modes of planar gradient-index hyperbolic secant waveguide / V.V. Kotlyar, A.A. Kovalev, Ya.R. Tri-andafilov, A.G. Nalimov // Computer Optics. - 2010. -V. 34, N 2. - P. 146-155. - (In Russian).
4. Musakaev, M.R. Higher-order polarization mode dispersion mathematical models for silica anisotropic optical waveguide / M.R. Musakaev, A.Kh. Sultanov// Computer Optics. - 2012. - V. 36, N 2. - P. 165-171. - (In Russian).
5. Moiseeva, N.M. The calculation of the dispersion characteristics of surface polaritons on the boundary of non-homogeneous film / N.M. Moiseeva, V.V. Yatshishen // Physics of wave propagation and radio systems. - 2004. -V. 7, N 2. - P. 30-32. - (In Russian).
6. Moiseeva, N.M. Calculation of dispersion characteristics of surface plasmons on the boundary a non-homogeus ani-sotropic medium / N.M. Moiseeva, V.V. Yatshishen // Physics of wave propagation and radio systems. - 2005. -V. 8, N 1. - P. 77-81. - (In Russian).
7. Gia Russo, D.P. Wave propagation in anisotropic thin film optical waveguides / D.P. Gia Russo, J.H. Harris // J. Opt. Soc. Am. - 1973. - Vol. 63, N 2. - P. 138-145.
8. Yeh, P. Optics of anisotropic layered media: a new 4* 4 matrix algebra / P. Yeh // Surface Science. - 1980. -Vol. 96, N 1-3. - P. 41-53.
9. Visnovsky, S. Magneto-optic effects in ultrathin structures at longitudinal and polar magnetizations / S. Visnovsky // Czech. J. Phys. - 1998. - Vol. 48, N 9 - P. 1083-1104. ISSN: 1572-9486 - (Online).
10. Walpita, L.M. Solutions for planar optical waveguide equations by selecting zero elements in a characteristic matrix / L.M. Walpita //J. Opt. Soc. Am. A. - 1985 -Vol. 2, N 4. - P. 595-602. - ISSN: 1520-8532 - (Online).
11. Chilwell, J. Thin-films field-transfer matrix theory of planar multilayer waveguides and reflection from prism-loaded waveguide / J. Chilwell, I. Hodgkinson // J. Opt. Soc. Am. A. - 1984. - Vol. 1. - P. 742-753. - (Online).
12. Liao, W. Explicit analysis of anisotropic planar waveguides by the analytical transfer-matrix method / W. Liao, X. Chen ,Y. Chen, Y. Xia, Y. Chen // J. Opt. Soc. Am. A. - 2004. - Vol. 21, N 11. - P. 2196-2204. - (Online).
13. Yee, K.S. Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media / K.S. Yee // IEEE Trans. Antennas and Propag. -1966. - Vol. 14, N 3. - P. 302-307.
14. Fallahkhair, A.B. Vector Finite Difference Modesolver for Anisotropic Dielectric Waveguide / A.B. Fallahkhair, K.S. Li, Th.E. Murphy // Journal of Lightwave Technol. -2008. - Vol. 26, N 11. - P. 1423-1431.
15. Koshiba, M. Finite-Element Solution of Anisotropic Waveguides with Arbitrary Tensor Permittivity / M. Koshiba, K. Hayata, M. Suzuki // Lightw. Technol. - 1986. -Vol. 4, N 2. - P. 121-126.
16. Tsuji, Y. Finite element beam propagation method for anisotropic optical waveguides / Y. Tsuji, M. Koshiba, N. Takimoto // J. Lightw. Technol. - 1999. - Vol. 17, N 4.
- P. 723-728.
17. Burmeister, E.F. Photonic integrated circuits optical buffer for packet-switchen networks / E.F. Burmaister, J.P. Mack, H.N. Poulsen, M.L. Masanovic, B. Stamenik, D.J. Blumental // Optics Express - 2009 - Vol. 17, N 8. - P. 6629-6635.
18. Moiseeva, N.M. Features of the dispersion characteristics of surface plasmons at the boundary of the inhomogeneous medium / N.M. Moiseeva, V.V. Yatshishen // Journal of Volgograd State University. Series 9, Issue 3. - 2003-2004.
- Part 2. - P. 37-41. - (In Russian).
19. Moiseeva, N.M. The dependence of the primary and pseudo-Brewster angles of orientation of the optical axis of the plane uniaxial anisotropic medium / N.M. Moiseeva // Lasers. Measurement. Information - 2011. - St.-Petersburg, Polytechnical University Publisher, 2011. - V. 3. - P. 430438. - (In Russian).
20. Moiseeva, N.M. Cross-polarization of electromagnetic waves reflected at the boundary of an anisotropic medium / N.M. Moiseeva // Lasers. Measurement. Information - 2011.
- St.-Petersburg: Polytechnical University Publisher, 2011. -V. 2. - P. 48-61. - (In Russian).
THE CALCULATION OF EIGENVALUES MODES OF THE PLANAR ANISOTROPIC WAVEGUIDES FOR VARIOUS ANGLES THE OPTICAL AXIS
N.M. Moiseeva Volgograd State University
Abstract
For a homogeneous anisotropic planar waveguide made solution of the Maxwell. We obtained fundamental matrices decision for the TE and TM waves in an anisotropic plane homogeneous layer. Was obtained the dependence of the phase shift of TM waves reflected from the boundaries of "anisotropic layer - isotropic medium" on the inclination angle the optical axis. The numerical solution of the dispersion equation TM waves in an anisotropic planar waveguide with different orientation of the optical axis in the plane of the wave. The asymptotic behavior of the dispersion curves was decided.
Key words: planar waveguide, anisotropy, the Maxwell equations, the dispersion, total internal reflection, the phase shift of the reflection of light.
Сведения об авторе
Моисеева Наталья Михайловна, 1966 года рождения, в 1988 году окончила Волгоградский государственный университет по специальности «Физика». Кандидат физико-математических наук (2005г.). Доцент кафедры судебной экспертизы и физического материаловедения Волгоградского государственного университета. Область научных интересов: оптика, электродинамика, дифференциальные уравнения, математическое моделирование, вычислительная математика.
E-mail: natalja2005@narod. ru .
Natalja Michailovna Moiseeva (b. 1966), was graduated from Volgograd State University in 1988, majoring in physics, candidate of science in physics and the mathematics (2005). She is the assistant professor of forensic science and materials science, Volgograd State University. Research interests: optics, electrodynamics, differential equations, mathematical modeling and computational mathematics.
Поступила в редакцию 26 сентября 2012 г.
