CC BY
52
ЗАТУХАНИЕ ВОЛНОВОДНЫХ МОД ТОНКИХ ПЛЕНОК: РАЗДЕЛЕНИЕ ВКЛАДОВ ПОГЛОЩЕНИЯ И ПОВЕРХНОСТНОГО РАССЕЯНИЯ СВЕТА
Сотский А.Б., Чудаковский П.Я.
Могилёвский государственный университет им. А.А. Кулешова
Аннотация
В рамках допущений о хорошей локализации энергии волноводных мод в тонкой плёнке и о малости по сравнению с длиной волны излучения среднеквадратичных отклонений и интервалов автокорреляции шероховатостей границ плёнки получено аналитическое решение задачи о разделении вкладов в коэффициенты затухания мод ТЕ-поляризации эффектов материального поглощения и поверхностного рассеяния света. На основе этого решения сформулирован алгоритм восстановления мнимой части комплексного показателя преломления плёнки, использующий вещественные части постоянных распространения и коэффициенты затухания двух ТЕ-мод плёнки, измеряемые методом волноводной спектроскопии. Эффективность подхода подтверждена расчётом дифракционных полей, возникающих при распространении волноводных мод вдоль плёнок с шероховатыми границами.
Ключевые слова: волноводная спектроскопия, волноводная мода, коэффициент затухания волноводных мод, тонкие плёнки.
Введение
Оптический контроль шероховатости поверхности и коэффициентов поглощения тонких плёнок играет важную роль в микро-, наноэлектронике и интегральной оптике [1 - 5]. Для соответствующих измерений перспективен метод волноводной спектроскопии, основанный на призменном возбуждении в плёнке волноводных мод [2, 3]. Его высокая чувствительность объясняется значительной длиной взаимодействия зондирующего излучения, распространяющегося вдоль плёнки в виде волноводной моды, с материалом плёнки и её границами. В настоящей работе основное внимание уделено проблеме восстановления указанным методом мнимой части комплексного показателя преломления плёнки kf. Эта проблема возникает при спектроскопии материалов тонких плёнок.
В методе волноводной спектроскопии осуществляется обработка углового распределения интенсивности отражённого светового пучка в призменном устройстве связи, дающая вещественную часть постоянной распространения h и коэффициент затухания аt возбуждаемой моды [2, 3]. Значение at определяется материальным поглощением света и его рассеянием на микроскопических дефектах волноведу-щей структуры. Данное рассеяние имеет место на внутренних неоднородностях материала плёнки и на её границах. В настоящей работе рассматривается последний из указанных механизмов рассеяния, который является, как правило, превалирующим [6].
Если предположить, что потери мощности моды плёнки, вызванные поглощением и рассеянием света, мультипликативны, а процесс рассеяния является стационарным и эргодическим, то можно принять приближение [6, 7]
а, = -21m h + a
(1)
где 1т И - мнимая часть постоянной распространения моды плёнки с идеально гладкими границами, связанная с материальным поглощением света, а^ -вклад в коэффициент затухания моды поверхностного рассеяния света.
Существующие алгоритмы решения обратных задач волноводной спектроскопии используют комплексные постоянные распространения мод плёнки И = Re И + i 1т И и сводят дело к определению толщины d и комплексного показателя преломления плёнки п- ikf из известных дисперсионных уравнений [2, 8, 9]. При этом значения d и п могут быть выражены через Re И , измеренные как минимум для двух мод, а значение kf - через 1т И одной из мод. Следовательно, определение 1т И, а значит, и kf по измеряемым значениям а^ предполагает такую обработку экспериментальных данных, которая позволила бы разделить вклады в at материального поглощения и поверхностного рассеяния света.
Ныне в методе волноводной спектроскопии используется приближение at = -21т И, которое представляет собой простейший вариант указанного разделения [2, 3, 6]. Но такой подход может быть оправдан только при исследовании сильно поглощающих плёнок (затухание мод >10 дБ/см) с почти гладкими поверхностями [2, 3, 6].
Попытка селективной оценки величин 1т И и а^ была предпринята в работе [10]. Соответствующая техника, основанная на анализе углового спектра энергетического коэффициента отражения светового пучка при призменном возбуждении мод плёнки, предполагает проведение измерений при различных толщинах воздушного зазора между призмой связи и исследуемой плёнкой и при этом довольно трудоёмка. Кроме того, известно, что на указанный спектр оказывает влияние эффект выхода излучения возбуждаемой моды за пределы призмы связи [2, 11, 12]. Поскольку в работе [10] названный эффект игнорируется, можно ожидать, что предлагаемый в ней подход будет эффективен только при условии достаточно сильного (>25 дБ/см [12]) затухания мод, когда данный эффект не существен. Ограничение применимости своего метода случаем сильно затухающих мод качественно оговорено и авторами работы [10].
Ниже представлен иной способ определения коэффициентов Imh, а; и kf. В нём указанные коэффициенты аналитически выражаются через значения Reh и at для двух мод ТЕ-поляризации, направляемых плёнкой (экспериментальное определение последних величин методом волноводной спектроскопии не вызывает принципиальных затруднений [2, 11, 12]). Эффективность предлагаемого метода подтверждена численными расчётами затухания ТЕ-мод различного порядка, выполненными с помощью спектрального разложения дифракционного поля по локальным модам нерегулярной волноводной структуры.
1. Разделение вкладов поглощения и поверхностного рассеяния в затухание мод Пусть мода ТЕ-поляризации, поле которой зависит от времени как exp(/rot), распространяется в z -направлении вдоль диэлектрической плёнки с двумя шероховатыми границами, отделяющими её от одной среды (например, воздуха) и другой среды (подложки). Эти границы описываются независимыми случайными функциями y = y+ (х, z) и y = y_ (x, z) соответственно. В большинстве представляющих практический интерес ситуаций можно считать, что средние значения указанных функций равны:
y+ (х, z) = 0 , y_ (х, z) = -d , и что выполнены неравенства
х-1^« i, ^ << i, ^ << i, (2)
y[Äyl и y[Äy-
где у/ Ду+ и у! Ду- - среднеквадратичные отклонения функций у+ (х, г) и у_ (х, z) ; }, В1+}, В(х}, В^ - интервалы автокорреляции этих функций в х -
и г-направлениях, Х0 - длина волны излучения в вакууме. Тогда потерю мощности моды dPs на отрезке вызванную рассеянием излучения на границах раздела сред, можно приближённо представить в виде суммы [7]
dP; = dP;(+) + dP;
(3)
где <Р/+} и } - мощности излучения, рассеянного верхней и нижней поверхностями плёнки. Согласно [7] данные величины пропорциональны интенсивности моды в окрестности соответствующих поверхностей:
dP;+-1 =-W\A(х,z) Y(0) dz,
dpy> =-W21 A(x, z)\"\Y(—d)|"dz,
(4)
(5)
где W\ и W2 - некоторые коэффициенты, определяемые статистическими свойствами функций у+ (х, г), у_ (х, г), показателями преломления сред и длиной волны излучения [7], А(х, г) - амплитуда моды, У (у) - поперечное распределение х -компоненты электрического поля моды регулярного плёночного волновода толщины d.
В аддитивном приближении [6] полные потери мощности моды на отрезке представляются в виде
<Р = <Ра + , (6)
где <Ра - потери мощности, вызванные поглощением излучения в средах. Как известно [13],
йРа = 21т hP(х, , (7)
где
Im h =-
k02 J Im [E(y)] |Y(y)2dy
—да
да -
2Reh J |Y(y)2dy
P(x,z) =| A(x, z)2 J |Y(y)2dy,
(8)
(9)
к0 = 2ЛХ-1 - волновое число света в вакууме, е( у) -относительная диэлектрическая проницаемость регулярной волноведущей структуры, в которой границы плёнки описываются уравнениями у=0 и у=- d. Предположим, что выполнено неравенство к/■ << п, характерное для большинства представляющих практический интерес ситуаций, и что энергия мод хорошо локализована в плёнке. Тогда в (8) можно приближённо положить
Re h = к0п г, 1т[е(у)] = -2п/к/ , (10)
в результате чего получим:
1т h = -к0 к г .
В соответствии с (3) - (7), (9)
Р = (21т к - Ш),
(11)
(12)
где
K =-
Y (2)"
J Y(y)2dy
w = w +
W2 |Y(—d)|'
l l2 Y (2)
или после учёта известного дисперсионного уравнения для ТЕ-мод плёночного волновода [13]
K =
2k[n) — (Re h)2 k22] ö
(ko d + Y—1 + Y—1)(n) — na2)' W = Wj + W2(n) — n2a)(n) — n2)—
(13)
где
Уа = ^(Яе к)2к- - п2а , у, = ^(Яе к)2к- - п2 ,
па и - показатели преломления сред, примыкающих к плёнке.
Заметим, что из (12) следует выражение (1), где
а, = Ш , (14)
а 1т к имеет вид (11).
Предположим теперь, что методом волноводной спектроскопии измерены вещественные части по-
—да
стоянных распространения и коэффициенты затухания двух ТЕ-мод плёнки, равные Re И,, Re hj и а('-1,
а^-1, где ' и . - номера мод. Если величины па и П известны (что обычно имеет место), то из дисперсионного уравнения для ТЕ-мод плёнки несложно рассчитать величины и d [2]. В результате будут
определены коэффициенты К1 и К. вида (13) для каждой из мод. После этого из (1), (13), (14) найдём:
(15)
а)=к^-.
kf = (2К0)-1 (а® -а?>), Ж = (а('> -а(Л)(К, - К,. )-1,
(16) (17)
где К, и а;; - коэффициенты К и а;, отнесённые к ' -й моде.
Выражения (11), (13), (15) - (17) определяют аналитическое решение обратной оптической задачи о разделении вкладов поглощения и поверхностного рассеяния поля волноводной моды и об определении коэффициента поглощения тонкой плёнки через измеряемые величины.
Полученное решение основано на ряде приближений, поэтому его корректность нуждается в проверке. Соответствующие оценки можно получить, численно рассчитав постоянные распространения и коэффициенты затухания ТЕ-мод структуры с заданными значениями d , па, ka, пг, kf, и , К
(па - ika, - iks - комплексные показатели преломления окружающих плёнку сред) и сравнив результат решения обратной задачи для kf с исходным значением kf. Подобные расчёты представлены в следующем разделе.
2. Численное моделирование затухания мод
Проблема расчёта поверхностного рассеяния полей мод оптических волноводов исследовалась в ряде работ (см. [7, 14, 15] и цитированную там литературу). При этом в основном рассматривались различные варианты теории возмущений для дифракционных полей. Подобные решения не имеют точно очерченных границ применимости, поэтому их использование для тестирования полученного выше решения обратной задачи проблематично.
Более точные вычислительные схемы для исследования дифракционных полей в продольно-нерегулярных волноводах были предложены в [11, 12, 16 - 19].
В [16] применительно к закрытым планарным волноводам рассмотрен метод Галёркина, в котором используется линейная сплайн-интерполяция компонент поля в направлении, ортогональном оси волновода. Расчёт амплитуд сплайнов на продольно-нерегулярном участке волновода сведён к решению системы обыкновенных дифференциальных уравнений. В [16] обоснована принципиальная сходимость предложенной расчётной схемы, но примеры её практического применения не представлены.
Открытые планарные волноводы с продольно-нерегулярными границами исследовались в [17, 18]. В данных работах использована ступенчатая интерполяция границ плёнки на нерегулярных участках. Такая модель позволяет рассматривать нерегулярный волновод как последовательность отрезков регулярных волноводов. Компоненты поля в каждом из этих волноводов представляются в виде рядов Фурье по координате, ортогональной оси волновода, причём границы интервалов периодичности гармоник Фурье располагаются внутри полностью согласованных слоёв (ПСС). Амплитуды гармоник Фурье выражаются через матрицы рассеяния, которые строятся из условий согласования компонент электромагнитного поля на границах регулярных волноводов. В указанных работах рассчитаны коэффициенты прохождения мод через нерегулярные участки тонких плёнок с шероховатыми границами [17] и с периодическими дифракционными решётками, вытравленными на поверхности плёнок [18]. В [17] исследованы только симметричные волноводы
( Па = п, , К = К ).
В [11, 12, 19] выполнен расчёт дифракционных полей, возникающих при возбуждении волноводных мод плёнок с шероховатыми границами посредством ограниченной призмы связи. Использованный при этом метод представляет собой модификацию известного метода спектральных разложений дифракционных полей в открытых планарных структурах по модам таких структур [14, 20]. В своей классической формулировке этот метод приводит к интегро-дифференциальным уравнениям относительно амплитуд мод непрерывного спектра, строгое решение которых представляет сложную проблему [14, 20]. В [11, 12, 19] вместо открытых плёночных волноводов рассматриваются закрытые планарные структуры с дискретным спектром мод. Отражение излучения от внешних границ данных структур минимизируется за счёт использования ПСС, которые моделируют открытое пространство. В результате, решение дифракционной задачи сводится к численному интегрированию систем либо обыкновенных дифференциальных, либо интегральных уравнений относительно амплитуд мод и не вызывает принципиальных затруднений.
Представленная ниже расчётная схема представляет собой развитие результатов работ [11, 12, 19]. Мы рассматриваем затухание мод в общем случае асимметричных плёночных волноводов (па ф ,
ка ф ). Дифракционное поле, связанное с рассеянием излучения на шероховатых границах плёнки, при каждом значении г представляется в виде разложения по локальным модам закрытой волноведу-щей структуры с ПСС. Данные моды строятся с учётом шероховатости границ плёнки, что обеспечивает достаточно быструю сходимость результатов (см. ниже). Кроме того, использование локальных мод, в отличие от рядов Фурье, позволяет естественным
образом выделить в каждом сечении волновода амплитуду той моды, затухание которой анализируется.
Как и авторы работ [14-18], мы ограничиваемся рассмотрением двумерной дифракционной задачи (д / д х = 0, у+ - (х, г) = у+ - (г)). Основанием для этого служит то, что двумерная модель шероховатости границ оптических волноводов приводит к корректному описанию реальных экспериментов по измерению затухания мод [14]. Кроме того, решение обратной задачи (11), (13), (15) - (17), как это очевидно из его вывода, идентично при рассмотрении как трёхмерной, так и двумерной моделей шероховатости границ плёнки. Поэтому для тестирования данного решения естественно исследовать более простую двумерную дифракционную задачу.
Рассматриваемая структура представлена на рис. 1. Здесь плёнка со средней толщиной d окружена средами с конечными средними толщинами ds и da. Плёнка имеет шероховатые границы на конечном участке 0 < г < гг оси 0г . Вне этого участка структура является регулярной и плоскопараллельной. На вход данного участка из области г<0 падает некоторая мода плёнки. При у=А и у=В структура ограничена идеально проводящими экранами. К экранам примыкают плоскопараллельные ПСС толщины Д. П'/////////////////////////////////////////у/л ПСС
AV//////^////////////////////////////////^//.
О гг
Рис. 1. Рассматриваемая структура
Тензоры относительной магнитной и диэлектрической проницаемости ПСС удовлетворяют соотношениям [21]
hi/ = h ji , hxy = hxz = h yz = 0
е.. = е.., е =е =е = 0,
ij jl > xy xz yz >
-in
hxx = hzz = ho =a eXP| —
hyy = he =
'1 ^
е =е = е =е.ц. ,
xx zz o jt^o >
е =е =
yy e
ho
(18)
(19)
где а > 0 - некоторая вещественная константа, еj -
диэлектрическая проницаемость среды, контактирующей с ПСС (j принимает значения 1 либо 3;
е1 = (п* - 'К )2, е3 = (П - 'ка )2).
Нетрудно показать, что при выполнении (18), (19) коэффициенты отражения плоских волн ТЕ-и ТМ-поляризации от границы раздела «диэлектрик-ПСС» равны
R = - exp(-2i ho^ jk0 Д)
(20)
где стj = cos 6 , 6 - угол падения волны. В соответствии с (20) при достаточно большом значении Д отражение излучения от ПСС становится пренебрежимо малым (|R ^ 0) независимо от его поляризации и угла падения. Таким образом, ПСС моделируют открытое пространство в закрытой волновод-ной структуре.
Согласно уравнениям Максвелла, x - компонента напряжённости электрического поля ТЕ-волн в области A < y < B (рис. 1) - описывается уравнением
д
д
д2
he f- h-1 — Ex + -т Ex + ^2е( y, z)Ex = 0 : д y д y д z
(21)
решаемым с нулевыми граничными условиями при у = А , у = В . Здесь в областях ПСС функция
е(у, г) принимает определённые выше значения е}, а и определяются по формулам (19). Вне этих областей = = 1, а функция е(у, г) описывает диэлектрическую проницаемость нерегулярной пла-нарной структуры.
Рассмотрим задачу Штурма-Лиувилля
д , д 2 2 he^T h- Т" + k^е(У, z) - h v( z) -у -У
9v (У, z) = 0,
(22)
(А, г) = 0, (В, г) = 0,
где г является фиксированным параметром. Эта задача имеет дискретный набор нетривиальных решений (мод), которые формируют полный набор функций на отрезке А < у < В [22]. При этом к„(г) имеет смысл постоянной распространения, а Ф„ (у, г) - поперечного распределения поля V -й локальной моды, эквивалентной моде продольно-регулярной планарной структуры с е(у) = е(у, г). В соответствии с (22), в I -м слое структуры
Фv(y, г) = Фvl (у - уI)] + +фVl (у - у1)]/ ,
(23)
где Фv1 = Фv (у1) , фVl = ^ / дА у=у1 , у! = у!(г) -
нижняя граница слоя (рис. 1), ст1 =^к02е1 - кV (г) , е1 -
диэлектрическая проницаемость слоя. Согласно (22) и (23), имеют место соотношения
Фv3 = Фv2 С0^СТ2^ ) + Ф^' ^П(СТ2^ ) / СТ2 "I
фVз = ^2СТ2 ^П(СТ2d2) + фV2 С0^СТ2d2 )I ' К (1 - Ч)
(24)
9V2 = К , ^v2 =Т
1СТ1 (1 + q)
h
o
е
q = ехр^/стД^ + ц„ А)], где = (г) - локальная толщина 1 -го слоя, ^ = ^ (г) - множитель, определяемый из условия нормировки полей мод
{ц-Ч(У, г) фц (У, г)(у =
(25)
(26)
где 5 - символ Кронекера (ортогональность функций с различными номерами V следует из (22)).
Из (20), (22), (23) получаем дисперсионное уравнение относительно И Дг):
ст3 ^3)ехр[-/СТ3^3 + ц„А)] +
+ст3 (iстзФvз +фVз)ехр^^ + ц„А)] = 0.
В соответствии с (24), левая часть (26) является целой трансцендентной функцией комплексной переменной И. В результате число корней уравнения (26) на полной комплексной плоскости неограниченно велико [22]. Значения этих корней в любой замкнутой области комплексной плоскости могут быть рассчитаны методом контурного интегрирования [23]. Последующее построение функций фv (у, г) на основании (23) - (25) не вызывает затруднений.
В силу полноты набора функций ф1, (у, г) на отрезке А < у < В решение уравнения (21) может быть представлено в виде равномерно сходящихся рядов
Ех (у, г) = УЦ (г) фv (у, г),
д
— Ех(у, г) = ^К (г)Фv (у, г), д г v
(27)
где Uv (г) и Vv (г) - некоторые функции г . Подставляя (27) в (21) и учитывая (25), приходим к обыкновенным дифференциальным уравнениям с переменными коэффициентами
А и = V -У D и ,
(г " " ^ уц ц
АV = -и И2-У D V
(г v v ^ vц ц
(28)
где е2 = (п/ - ikf )2.
Ограничимся в бесконечных рядах (27) К слагаемыми. В этом случае (28) представляет собой замкнутую систему 2К уравнений. Поскольку вне отрезка г0 < г < гг рассматриваемая волноведущая структура регулярна, в соответствующих областях выполняются соотношения (у± / (г = 0, Dvц (г) = 0 . Тогда согласно (28)
Uv = av ехр(-/И ^ г) + Ь ехр(/И ^ г), 1
Vv = -/И V0)[av ехр(-/И V0)г) - ^ ехр(/И V0)г)] при г < г0 и
Uv= с ехр(-/И V0'г), 1 ^=-/И V0) Cv ехр(-/И V0) г)1
(30)
(31)
при г > гг. Здесь И ^ - постоянная распространения V -й моды регулярных планарных структур, расположенных в областях г < г0 и г > гг; av, bv, cv -константы, имеющие смысл амплитуд падающих, отражённых и прошедших мод указанных регулярных структур. Значения av задаются внешними источниками излучения, а амплитуды bv и cv должны определяться в результате решения дифракционной задачи. Поскольку нас интересует затухание конкретной моды плёнки с номером ц, в (30) следует положить
av = a5vц , (32)
где а - некоторая постоянная.
Для численного решения задачи (28) - (32) заменим её эквивалентными интегральными уравнениями. Эти уравнения получаются из уравнений, приведённых в [11, 19], путём исключения членов, описывающих поле светового пучка, падающего на волноведущую структуру, и имеют вид
ЦV (г) = ^ ехр (-/'И ^ г) +
1 г 1г; (33)
+- { Ф^ (г, г') & — { ФV+) (г, г') дг\ 22
2
где Dцц= 0,
Яц= Яц (г) = К2[И2(г) - И 2(г)]-1 х
г де
х{ ц-Ч (y, г)фц (y, г) дг (у ^цХ
А+А дг
В последнее выражение входит производная ступенчатой функции де / дг , которая выражается через дельта-функции Дирака. Интегрирование этих функций даёт аналитический результат:
Яц (г) =
(е2 -^"Г^Кфц) +
г1-т V н-/у= у+ (г)
(г
у=у+ (г )
+(е1 -^^^фц)
аг
' у=у- (г )
(29)
Иц2(г) - И V(г)
Vv (г) = -/И V0) av ехр (-/И V0) г) -
-1 /И(0) 2 v
{ ФV-) (г, г ' + { ФV+) (г, г') (г'
где
(34)
Ф^ (г, г') = ехр[±/И ^(г - г')] [^ (г') ± Gv (г')],
Fv (г) = -у Дц (г)Ц ц (г),
2
Gv (г) = (/И V0))-1 [(И ^ - И 2(г)]ЦV (г) --(/И V0))-1 У Ац (г V (г).
Целесообразность перехода от задачи (28) - (32) к уравнениям (33), (34) объясняется тем, что они
А
V
2
г
г
В-А
0
г
2
К
0
содержат малые параметры (И V)) - И Д г), Dvц (г) и
допускают решение методом итераций (см. ниже).
Следует отметить, что объём вычислений в рассматриваемой расчётной схеме возрастает при увеличении порядка редукции рядов (27). Значение К при заданной точности решения определяется спектром используемых мод фv (у, г), который, в свою очередь, зависит от толщины ПСС А. В частности, при А ^ ж данный спектр становится квазинепрерывным, поэтому К ^ ж . С другой стороны, при А ^ 0 спектр мод является разреженным (К можно выбрать конечным), но получаемое решение будет искажено вследствие паразитного отражения излучения от вспомогательных идеально проводящих экранов. Как показано в [19], данное отражение становится пренебрежимо малым при условиях а = 1, А > 1,1Х0. Влияние выбора К на получаемое решение при конкретном значении А из указанного диапазона обсуждается ниже.
С целью исследования затухания мод мы использовали для границ плёнки (рис. 1) модель у2 = у- (г) = ^ + Ау- (г) (0 < г < гг), у2 = у- (г) = ^ (г < 0, г > гг), уз = у+ (г) = Ау+ (г) (0 < г < гг). уз = у+ (г) = 0 (г < 0, г > гг). Здесь Ау- (г) и Ау+ (г) - независимые случайные функции. В эквидистантных узлах гj = ¡Аг
(- = 0, N, Аг = гг / N) отрезка 0 < г < гг они принимают значения
Ау±°).(ю = 0 , Ау±) = Ь(г±(л - 0,5) (N > - > 0), где Ау± -1 = Ау± (гj); г+( ¡) и г-¡) - независимые случайные числа, равномерно распределённые на отрезке (0, 1) (при расчётах г+(¡) и г-¡) задавались равномерным генератором псевдослучайных чисел); 0,5Ь - максимальная амплитуда колебаний функций Ау_ (г) и Ау+ (г) . Между узлами (г,- < г < г-+1) данные функции линейны:
Ау± (г) = Ау±¡) + (г - г, )(Ау±-+1) - Ау±5)/ Аг .
В рамках принятой модели автокорреляционные функции р± (ст) = Ау± (г)Ау± (г + ст) равны
р±(ст)=12
- -3с+1С2 -^(1 -о3 6 2 2 2
(С< 1),
Р± (ст) = Ь (272 (1 <С< 2), р± (ст) = 0 (С> 2),
где С = ст / Аг. Согласно последним формулам, для интервалов корреляции можно использовать
оценку В^т -1 — А г.
При расчётах мы полагали Х0 = 0,6328 мкм, а = 1, А = 1 мкм , Ь = 0,05 мкм, Аг = 0,1 мкм, п = 1,457 ,
К = 10 , ds = 3 мкм , = 1,5147, К, = 1,1326-10-5, па = 1, Ка = 0 , da = 3 мкм. Здесь ds и da - толщины окружающих плёнку сред в диапазонах г < 0 , г > гг (рис. 1), а значения показателей преломления соответствуют плёнке из стекла К8, нанесённой на кварцевую подложку и контактирующей с воздухом. Вычисления выполнены для плёнок двух толщин: d = 3 мкм и d = 1,5 мкм , которые направляют, соответственно, 4 и 2 собственных моды. Опишем подробно результаты расчётов, относящиеся к плёнке с d = 3 мкм (рис. 2-4).
Ключевую роль при решении уравнений (33), (34) играет построение спектра комплексных постоянных распространения мод И, (г). Фрагмент такого спектра, рассчитанный методом контурного интегрирования, приведён на рис. 2. Он относится к регулярной планарной структуре (г < 0, г > гг). Четыре крайние правые точки спектра (вставка на рис. 2) соответствуют собственным модам плёнки с номерами 0, 1, 2, 3. Постоянные распространения этих мод равны:
И0(0) = К0 (1,51172454 - / -1,131 -10-5), И1(0) = К0(1,50281694 - /-1,123-10-5), И20) = К0 (1,48808230 - /-1,103-10-5), И3(0) = К0 (1,46806775 - / -1,032 -10-5).
1тИ/ко-0
-0,5 -1,0 -1,5 -2,0
(35)
_1_
. • • 1 -
0 — • ■ ■ м Л
•
-0,01 • • • 1
V и 1,5 )
1 ЯеЫко
0 0,4 0,8 1,2 1,6
Рис. 2. Расположение постоянных распространения мод на комплексной плоскости
В качестве базиса разложений (27) мы использовали К мод, постоянные распространения которых И<^) наиболее близки по модулю к постоянной распространения моды, падающей на вход нерегулярного участка плёнки. Данный выбор объясняется тем, что коэффициенты связи мод (28) обратно пропорциональны разностям Иц (г) - И ^(г).
Численное решение уравнений (33), (34) было выполнено методом итераций. В качестве нулевого приближения выбирались функции V,(г) = -/И(V0)Uv(г),
и V (г) = а5 ехр(-/И ^ г).
Практически стопроцентная сходимость решений при всех использованных значениях ^ указанных ниже, наблюдалась при 6 - 7 итерациях.
На рис. 3 серыми быстро колеблющимися кривыми представлены расчётные зависимости |иц (г)/«|2 на нерегулярном участке плёнки, соответствующие поступающим на вход данного участка модам с номерами ц = 0 и ц = 1. Данные кривые рассчитаны
при k = 70 и конкретной реализации случайных функций Ду- (z) и Ду+ (z) (фрагмент графика толщины плёнки приведён на вставке рис. 3).
\и0!а\2 \иМ2
600
г,гг,мкм
Рис. 3. Затухание основной (0) и первой высшей (1) мод на нерегулярном участке плёнки.
На вставке - фрагмент зависимости d(z)
Известно, что увеличение порядка (или номера) моды сопровождается усилением рассеяния её излучения на шероховатых границах плёнки [7]. Этим объясняется то, что колебания зависимости |и1^)/а|2 на рис. 3 являются более выраженными, чем в случае зависимости |и0 ^)/а|2. Интересно, что для обеих рассмотренных мод амплитуда указанных колебаний максимальна в начале нерегулярного участка, причём |и0д ^)/а|2 Ф 1. Эта особенность объясняется отражением мод от нерегулярного участка плёнки.
Заметим, что зависимости, аналогичные приведённым на рис. 3, измеряют в экспериментах по исследованию затухания волноводных мод методом фотомет-рирования их треков, при этом коэффициенты затухания мод определяют путём аппроксимации данных зависимостей экспонентами вида ехр(-а^) [24]. Результаты подобной аппроксимации зависимостей 0 и 1 на рис. 3 представлены чёрными линиями. Этим линиям соответствуют значения а1 = 2,387 -10-4мкм-1 (ц = 0) и а( = 2,791 •10-4мкм-1 (ц = 1), которые определены методом наименьших квадратов.
Строго говоря, рассмотренная оценка а1 является приближённой, так как в функцию иц (¿) при z < zr дают вклад моды, распространяющиеся как в положительном, так и в отрицательном направлениях оси 0z. Более корректно коэффициент затухания ц -й моды может быть найден путём экспоненциальной интерполяции зависимости |иц ^г)/а|2, поскольку
при z > zr моды, распространяющиеся в отрицательном направлении оси 0z, отсутствуют (см. (31)) [10]. Экспериментально функция |иц ^г)/а|2 может быть оценена известным методом двух призм [25], в котором с помощью подвижной призмы связи осуществляется полный вывод излучения моды, прошедшей отрезок волновода контролируемой длины zr. Погрешность восстановления а этим методом составляет около 0,02 дБ/см [25].
Результаты расчёта функций | и01 ^г)/а|2 представлены на рис. 3 крестиками. Как видно, эти данные практически не отклоняются от графиков описанных выше экспонент, что подтверждает корректность подхода, использованного в [24]. Аналогичная картина имеет место и при рассмотрении мод плёнки с номерами 2 и 3 (эти графики не приведены). Представленные ниже значения а найдены обоими описанными способами, которые привели к идентичным (в приведённых ниже десятичных разрядах) результатам. Расчёты выполнены в диапазоне zr <500 мкм.
Рис. 4 иллюстрирует внутреннюю сходимость принятой вычислительной схемы.
1оЕ1о{а(1ко)
-5 -
-6 -
-7 ~
-8
А *
О
Д
А * * О
т
т
А *
О
®
А д А А А А А А А
А * * * * * * * *
О О О О О О О О О о
® <§> ® <§> ® ® ® ® ® ®
®
т ц=0 О
* р.=2
А ц=3
_1_
_1_
_1_
0 20 40 60 80
Рис. 4. Зависимости расчётных коэффициентов затухания мод плёнки от порядка редукции рядов (29)
Как видно из рис. 4, использование в качестве базиса для представления дифракционных полей локальных мод нерегулярной волноведущей структуры обеспечивает сравнительно быструю сходимость решения (практически стопроцентная стабилизация результатов для а наблюдается при k > 60). Для сравнения, в [17, 18] устойчивые данные для коэффициентов затухания мод получались при использовании более 200 гармоник Фурье.
Аналогичные результаты мы получили и при исследовании плёнки с d = 1,5 мкм, которая направляет две собственные моды с постоянными распространения
й0(0) = ^(1,50490664-/1,113-10-5),1 = ^(1,47660575 - /1,009 • 10-5).
Обратимся теперь к оценке корректности решения обратной задачи (11), (13), (15) - (17). Исходные данные для такой оценки, полученные описанным методом, представлены в табл. 1.
(36)
Таблица 1. Модовые характеристики плёнок
В табл. 2 приведены результаты расчёта коэффициента Ж, отвечающего за поверхностное рассеяние света и итоговое решение обратной задачи о восстановлении мнимой части показателя преломления плёнок. Значения К и Ж рассчитаны на основании выражений (16) и (17) при использовании пар мод с номерами, указанными во втором столбце табл. 2 (при вычислении К на основании (16) конкретный выбор одного из двух указанных номеров мод не существен). Использование основной моды (с номером 0) во всех парах мод в табл. 2 объясняется тем, что для данной моды приближение (10), (11) имеет максимальную точность [6].
Таблица 2. Восстановление мнимой части показателя преломления плёнок
Плёнка Моды 103W 105kf
d = 3мкм 0, 1 3,2823 1,1331
0, 2 3,2788 1,1332
0, 3 3,1500 1,1359
d = 1,5 мкм 0, 1 3,5945 1,1311
Согласно табл. 1, величины а /(2К0), относящиеся к модам различного порядка, существенно отличаются как от соответствующих значений -1т И/К) (см. (35), (36)), так и друг от друга. Это указывает на грубый характер упомянутого во введении приближения а=-2 1т И. Из табл. 1 следует также резкое возрастание затухания моды при увеличении её номера. В соответствии с (1), (35), (36), этот эффект можно связать только с поверхностным рассеянием света. В принятой модели затухания мод данный эффект описывается выражениями (14), (17).
Количественно судить о корректности приближения (14) позволяет табл. 2. Из неё видно, что, несмотря на существенный разброс значений -1 и К, для различных мод (см. табл. 1), коэффициент Ж, характеризующий рассеивающие свойства поверхностей плёнки, остаётся почти постоянным. Некоторый рост Ж при уменьшении ё объясняется, по-видимому, усилением корреляции излучения, рассеянного двумя границами плёнки, которая при записи (3) игнорировалась.
Об эффективности полученного решения обратной оптической задачи по определению К свидетельствует последний столбец табл. 2. Представленные в нём данные позволяют сделать вывод о целесообразности использования пары мод низшего порядка. В этом случае модули отклонения восстанов-
ленных значений kf от исходного значения kf = 1,1326405 минимальны и равны 5^10-9 (четы-рёхмодовая плёнка) и 15-10-9 (двухмодовая плёнка). При таких уровнях погрешности предлагаемого метода основным источником ошибок при восстановлении kf будут погрешности измерения at. Например, при расчёте kf на основании (16) экспериментальной погрешности at, равной 0,02 дБ/см, будет соответствовать погрешность kf, равная 40-10~9 .
Заключение При условиях малых по сравнению с длиной волны излучения среднеквадратичных отклонений и интервалов автокорреляции шероховатостей границ тонкой плёнки (2) и хорошей локализации энергии волновод-ных мод в плёнке (10), (11) получено аналитическое решение задачи о разделении вкладов в затухание вол-новодных мод эффектов поглощения и поверхностного рассеяния света. В рамках данного решения мнимая часть показателя преломления плёнки может быть рассчитана по формуле (16) после измерения вещественных частей постоянных распространения и коэффициентов затухания двух ТЕ-мод плёнки. Корректность предложенного подхода подтверждена численными расчётами дифракционных полей, возникающих при распространении волноводных мод вдоль плёнок с шероховатыми границами.
Литература
1. Глянько, М.С. Программное обеспечение для устройства контроля чистоты и шероховатости оптических подложек / М.С. Глянько, П.Ю. Изотов // Компьютерная оптика. - 2012. - Т. 36, № 2. - С. 242-248.
2. Khomchenko, A. Waveguide spectroscopy of thin films / A. Khomchenko. - Elsevier, 2005. - 236 p.
3. Cardin, J. Determination of refractive index, thickness, and the optical losses of thin films from prism-film coupling measurements/ J. Cardin, D. Leduc // Appl. Opt. -2008. - Vol. 47, N 7. - P. 894-900.
4. Биленко, Д.И. Определение оптических свойств и толщины нанослоёв по угловым зависимостям коэффициента отражения / Д.И. Биленко, А.А. Сагайдач-ный, В.В. Галушка, В.П. Полянская // ЖТФ. - 2010. -Т. 80, № 10. - С. 89-94.
5. Тарасов, И.А. Эллипсометрическая экспресс-методика определения толщины и профилей оптических постоянных в процессе роста наноструктур Fe/Si02/Si(100) / И.А. Тарасов, Н.Н. Косырев, С.Н. Варнаков, С.Г. Овчинников, С.М. Жарков, В.А. Швец, С.Г. Бондаренко, О.Е. Терещенко // ЖТФ. - 2012. - Т. 82, № 9. - С. 44-48.
6. Хансперджер, Р. Интегральная оптика: Теория и технология / Р. Хансперджер. - М.: Мир, 1985. - 384 с.
7. Унгер, Х.Г. Планарные и волоконные оптические волноводы / Х.Г. Унгер. - М.: Мир, 1980. - 656 с.
8. Сотский, А.Б. Анализ распределения интенсивности отражённого пучка в схеме призменного возбуждения диэлектрических волноводов / А.Б. Сотский, А.А. Ро-маненко, А.В. Хомченко, И.У. Примак // Радиотехника и электроника. - 1999. - Т. 44, № 6. - С. 687-695.
9. Хомченко, А.В. Волноводный метод измерения параметров тонких плёнок / А.В. Хомченко, А.Б. Сотский, А.А. Романенко, Е.В. Глазунов, А.В. Шульга // Журнал технической физики. - 2005. - Т. 75, № 6. - С. 98-106.
Плёнка Мода 105 at / (2k0 ) 103 K / k0
d = 3мкм 0 1,2020 0,419709
1 1,4054 1,65903
2 1,7286 3,63191
3 2,0684 5,92077
d = 1,5 мкм 0 1,5765 2,47799
1 2,7138 8,80641
10. Zhang, X.-J. Simple and efficient technique for evaluating the optical losses from surface scattering and volume attenuation in a thin film / Xi-Jing Zhang, Xi-Zhi Fan, Hui-Tian Wang, Jing-Liang He, N.B. Ming // Opt. Express. - 2002. - Vol. 10, N 25. - P. 1485-1490.
11. Сотский, А.Б. Отражение светового пучка от призмы связи. Ч.1 / А.Б. Сотский, П.Я. Чудаковский, И.У. Примак, Л.И. Сотская // Вестник Могилёвского государственного университета им. А.А. Кулешова. - 2012. -Т. 39, № 1. - С. 44-59.
12. Сотский, А.Б. Отражение светового пучка от призмы связи. Ч.2 / А.Б. Сотский, П.Я. Чудаковский, И.У. Примак, Л.И. Сотская // Вестник Могилёвского государственного университета им. А.А. Кулешова. - 2012. -Т. 40, № 2. - С. 45-59.
13. Тамир, Т. Интегральная оптика / Т. Тамир. - М.: Мир, 1978. - 344 с.
14. Маркузе, Д. Оптические волноводы / Д. Маркузе. -М.: Мир, 1974. - 576 с.
15. Лахад, Л.С.А. Рассеяние лазерного излучения в градиентных волноводах с шероховатой границей / Л.С.А. Лахад, А.Н. Осовицкий // Вестник РУДН. серия: Физика. - 2005. - Т. 13, № 1. - С. 60-67.
16. Ильинский, А.С. Математические модели в электродинамике / А.С. Ильинский, В.В. Кравцов, А.Г. Свешников. - М.: Высшая школа, 1991. - 224 с.
17. Elson, J.M. Propagation in planar waveguides and the effects of wall roughness / J.M. Elson // Opt. Express. -2001. - Vol. 9, N 9. - P. 461-475.
18. Silberstein, E. Use of grating theories in integrated optics / E. Silberstein, P. Lallane, J.P. Hugonin, Q. Cao // J. Opt. Soc. Am. A. - 2001. - Vol. 18, N 11. - P. 2865-2875.
19. Сотский, А.Б. Дифракция света в планарных структурах с полностью согласованными слоями / А.Б. Сотский, П.Я. Чудаковский, И.У. Примак // Вестник Моги-лёвского государственного университета им. А.А. Кулешова. - 2010. - Т. 29, № 1. - С. 45-50.
20. Шевченко, В.В. Плавные переходы в открытых волноводах / В.В. Шевченко. - М.: Наука, 1969. - 192 с.
21. Sacks, Z.S. A perfectly matched anisotropic absorber for use as an absorbing boundary condition / Z.S. Sacks, D.M. Kingsland, R. Lee, J. Lee // IEEE Trans. Antennas Propagat. - 1995. -Vol. 43, N 12. - P. 1460-1463.
22. Шевченко, В.В. О спектральном разложении по собственным и присоединённым функциям одной несамосопряженной задачи типа Штурмa-Лиувилля на всей оси / В.В. Шевченко // Дифференциальные уравнения. - 1979. - Т. 15, № 11. - С. 2004-2020.
23. Сотский, А.Б. Теория оптических волноводных элементов / А.Б. Сотский. - Могилёв: МГУ им. А.А. Кулешова, 2011. - 456 с.
24. Strasser, T.A. Optical loss measurement of low-loss thin-film waveguides by photographic analysis / T.A. Strasser, M.C. Gupta // Appl. Opt. - 1992. - Vol. 31, N 12. -P. 2041-2046.
25. Weber, H.P. Loss measurement in thin-film optical waveguides / H.P. Weber, F.A. Dunn, W.N. Leibolt // Appl. Opt. - 1973. - Vol. 12, N 4. - P. 755-757.
References
1. Glyanko, M.S. Software for the device for detection of cleanliness and roughness of optical substrates / M.S. Glyanko, P.Yu. Izotov // Computer Optics. - 2012. - Vol. 36, N 2. - P. 242-248. - (In Russian).
2. Khomchenko, A. Waveguide spectroscopy of thin films / A. Khomchenko. - Elsevier, 2005. - 236 p.
3. Cardin, J. Determination of refractive index, thickness, and the optical losses of thin films from prism-film coupling measurements/ J. Cardin, D. Leduc // Appl. Opt. -2008. - Vol. 47, N 7. - P. 894-900.
4. Bilenko, D.I Determination of optical properties and thickness of the nanolayers on the angular dependence of the reflection coefficient / D.I. Bilenko, A.A. Sagaidach-ny, V.V. Galuska, V.P. Polyanskaya // The Journal of Technical Physics. - 2010. - Vol. 80, N 10. - P. 89-94. -(In Russian).
5. Tarasov, I.A. Ellipsometric express method for determining the thickness and optical constants of the profiles in the growth process of nanostructures Fe/SiO2/Si (100) / I.A. Tarasov, N.N. Kosyrev, S.N. Varnakov, S.G. Ovchin-nikov, S.M. Zharkov, V.A. Shvets, S.G. Bondarenko,
0.E. Tereshchenko // The Journal of Technical Physics. -2012. - Vol. 82, N 9. - P. 44-48. - (In Russian).
6. Hunsperger, R. Integrated optics: Theory and Technology / R. Hunsperger. - Moscow: "Mir" Publisher, 1985. -384 p. - (In Russian).
7. Unger, H.G. Planar optical waveguides and fibers / H.G. Unger. - Moscow: "Mir" Publisher, 1980. - 656 p. -(In Russian).
8. Sotsky, A.B. Analysis of the reflected beam intensity distribution in the scheme of prism excitation of dielectric waveguides / A.B. Sotsky, A.A. Romanenko, A.V. Khom-chenko, I.U. Primak // Radiotekh. Electron. - 1999. -Vol. 44, N 6. - P. 687-695. - (In Russian).
9. Khomchenko, A.V. The waveguide method of measuring the parameters of thin films / A.V. Khomchenko, A.B. Sot-sky, A.A. Romanenko, E.V. Glazunov, A.V. Shulga // The Journal of Technical Physics. - 2005. - Vol. 75, N 6. -P. 98-106. - (In Russian).
10. Zhang, X.-J. Simple and efficient technique for evaluating the optical losses from surface scattering and volume attenuation in a thin film / Xi-Jing Zhang, Xi-Zhi Fan, Hui-Tian Wang, Jing-Liang He, N.B. Ming // Opt. Express. - 2002. - Vol. 10, N 25. - P. 1485-1490.
11. Sotsky, A.B. Light beam reflection from the prism-film coupler. Part I. Theoretical model / A.B. Sotsky, P.J. Chuda-kouski, I.U. Primak, L.I. Sotskaya // Vestnik Mogilev State University. - 2012. - Vol. 19, N 1. - P. 44-59. - (In Russian).
12. Sotsky, A.B. Light beam reflection from the prism-film coupler. Part II. The solution of the inverse problem for the waveguide / A.B. Sotsky, P.J. Chudakouski,
1.U. Primak, L.I. Sotskaya // Vestnik Mogilev State University. - 2012. - Vol. 40, N 2. - P .45-59. - (In Russian).
13. Tamir, T. Integrated optics / T. Tamir. - Moscow: "Mir" Publisher, 1978. - 344 p. - (In Russian).
14. Marcuse, D. Light transmission optics / D. Marcuse. -Moscow: "Mir" Publisher, 1972. - 576 p. - (In Russian).
15. Lahad, L.S.A. Dispersion of light in the gradient wave guides with rough border / L.S.A. Lahad, A.N. Osovitsky // Vestnik RUDN (Physics). - 2005. - Vol. 13, N 1. -P. 60-67. - (In Russian).
16. Ilyinsky, A.S. Mathematical models of electrodynamics / A.S. Ilyinsky, V.V. Kravtsov, A.G. Sveshnikov. - Moscow: "Vysshaya Shkola" Publisher, 1991. - 224 p. - (In Russian).
17. Elson, J.M. Propagation in planar waveguides and the effects of wall roughness / J.M. Elson // Opt. Express. -2001. - Vol. 9, N 9. - P. 461-475.
18. Silberstein, E. Use of grating theories in integrated optics / E. Silberstein, P. Lallane, J.P. Hugonin, Q. Cao // J. Opt. Soc. Am. A. - 2001. - Vol. 18, N 11. - P. 2865-2875.
19. Sotsky, A.B. Diffraction of light in planar structures with perfectly matched layers / A.B. Sotsky, P.J. Chudakouski,
I.U. Primak // Vestnik Mogilev State University. - 2010. -Vol. 29, N 1. - P. 45-50. - (In Russian).
20. Shevchenko, V.V. Smooth transitions in open waveguides / V.V. Shevchenko. - Moscow: "Nauka" Publisher, 1969. - 192 p. - (In Russian).
21. Sacks, Z.S. A perfectly matched anisotropic absorber for use as an absorbing boundary condition / Z.S. Sacks, D.M. Kingsland, R. Lee, J. Lee // IEEE Trans. Antennas Propagat. - 1995. -Vol. 43, N 12. - P. 1460-1463.
22. Shevchenko, V.V. About a spectral resolution on characteristic and associated functions of one non self-conjugate problem of type of Sturm-Liouville on all axis /
V.V. Shevchenko // Differential Equations. - 1979. -Vol. 15, N 11. - P. 2004-2020. - (In Russian).
23. Sotsky, A.B. Theory of optical waveguide elements / A.B. Sotsky. - Mogilev: MSU Publisher, 2011. - 456 p. -(In Russian).
24. Strasser, T.A. Optical loss measurement of low-loss thin-film waveguides by photographic analysis / T.A. Strasser, M.C. Gupta // Appl. Opt. - 1992. - Vol. 31, N 12. -P. 2041-2046.
25. Weber, ftP. Loss measurement in thin-film optical waveguides / H.P. Weber, F.A. Dunn, W.N. Leibolt // Appl. Opt. - 1973. - Vol. 12, N 4.- P. 755-757.
ATTENUATION OF WAVEGUIDE MODES OF THIN FILMS: SEPARATION OF CONTRIBUTIONS OF ABSORPTION AND SURFACE SCATTERING OF THE LIGHT
A.B. Sotsky, P.Ya. Chudakouski Mogilev State A. Kuleshov University
Abstract
Analytical solution for the problem on separation of contributions of the material effects of absorption and surface scattering of the light in the attenuation coefficients of TE polarization modes is obtained under the assumption of a good localization of the energy of waveguide modes in a thin film and mean square deviations and the autocorrelation intervals of rough boundaries of film are small than the wavelength of the radiation. Based on this solve the algorithm of reconstruction of the imaginary part of the complex refractive index for the film is formulated. This algorithm uses of the real parts of the propagation constants and attenuation coefficients for the two TE modes of film, that are measured by means the waveguide spectroscopy. Effectiveness of the approach is confirmed by calculation of the diffraction fields arising at the propagation of waveguide modes along the film with rough boundaries.
Key words: waveguide spectroscopy, waveguide mode, attenuation coefficient of waveguide modes, thin film.
\ '
Сведения об авторах
Сотский Александр Борисович, 1953 года рождения. В 1976 году окончил Белорусский государственный университет. Доктор физико-математических наук (2002 г.), профессор, работает профессором кафедры экспериментальной и теоретической физики Мо-гилёвского государственного университета им. А.А. Кулешова. Область научных интересов: дифракционная оптика, интегральная и волоконная оптика, обратные оптические задачи.
Alexandr Borisovich Sotsky (b. 1953). In 1976 graduated from Belarusian State University. Doctor of Physical and Mathematical Sciences (2002), professor, a professor of the Department of Experimental and Theoretical Physics in Mogilev State A. Kuleshov University. Research interests: diffractive optics, integrated and fiber optics, optical inverse problems.
Чудаковский Павел Яковлевич, 1985 года рождения. В 2008 году окончил Моги-лёвский государственный университет имени А.А. Кулешова. Магистр физико-математических наук (2009), аспирант кафедры экспериментальной и теоретической физики Мо-гилёвского государственного университета имени А.А. Кулешова. Область научных интересов: волноводная оптика.
E-mail: pawel206265@mail.ru.
Pavel Yakovlevich Chudakouski, (b. 1985). In 2008 graduated from Mogilev State A. Kuleshov University. Master of Mathematical and Physical Sciences (2009), post-graduate student of the Department of Experimental and Theoretical Physics in Mogilev State A. Kuleshov University. Research interests: waveguide optics.
Поступила в редакцию 9 октября 2012 г.
