Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

rQ3PkAPUVO

.pdf
Скачиваний:
1
Добавлен:
13.02.2021
Размер:
744.91 Кб
Скачать

Исключив из этих уравнений отношение A/B, получим характеристическое уравнение для ТЕ-мод:

tg2aq = ( p + r )q . q2 pr

С точностью до Nπ это уравнение совпадает с характеристическим уравнением, полученным методами геометрической оптики.

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

C exp(-rx) ,

x ³0;

 

0³ x ³- 2a;

H y = C cos qx + D sin qx,

(C cos 2aq - Dsin 2aq)exp p(x + 2a),

- 2a ³ x;

 

 

 

 

 

(-rC / n2 ) exp(-rx) ,

 

 

 

3

E

 

 

z

= (q / n2 )(-C sin qx + D cos qx) ,

 

 

1

 

 

( p / n2 )(C cos 2aq - D sin 2aq) exp p(x + 2a) ,

 

 

 

2

 

 

 

 

Характеристическое уравнение для ТМ-мод:

tg2aq = (n32 p + n22r) n12q . n22n32q2 n14 pr

x ³0;

0 ³ x ³ - 2a ;

- 2a ³ x.

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

Для обозначения моды, в распределении поля которой имеется N узлов, используется обозначение TEN или TMN. Для низшей (основной) моды N = 0. На рис. 2.2 приведены распределения электрического поля для ТЕ0- и ТЕ1-мод.

Определим число направляемых

мод, используя условие отсечки β =

Рис. 2.2

= n2k, которое на языке модового опи-

11

(2a +1/ p +1/ r) = w

сания означает нарушение локализации поля излучения и расплывание этого поля в область с n2. Общее число направляемых ТЕ и ТМ-мод

 

1

 

 

n

2

n2

1 / 2

 

 

 

 

 

 

 

M =

 

 

2ν − arctg G

 

2

3

 

 

,

 

 

2

2

 

π

 

 

 

 

 

 

n1

n2

 

 

 

 

 

 

 

 

 

 

 

 

int

где G = 1 для ТЕ-мод и G = (n1/n3)2 для ТМ-мод;

ν = ak(n12 n22)1/2 – нор-

мированная частота.

 

 

 

 

 

 

 

 

 

 

Усредненное по времени значение потока энергии Р в волноводе равно интегралу по поперечному сечению волновода от компоненты вектора Пойтинга Sz в направлении оси z:

1

 

 

 

 

P = Sxdx =

Re(E × H * )z dx.

2

−∞

−∞

 

 

 

Для ТМ-мод компонента Sz задается выражением

Sz = − 1 Ey H x* =

β

 

E y .

 

 

 

 

 

 

 

 

2

2

 

 

 

2ωμ0

 

 

 

 

 

 

 

 

 

Интеграл в выражении для Р следует разложить на три слагаемых, соответствующих трем слоям волновода. Тогда полная мощность задается выражением

 

β A2

q2

+ r 2

 

 

1

 

1

 

 

P = P + P + P =

 

 

 

 

 

 

 

 

2a +

 

+

 

 

,

 

 

 

 

 

2

 

 

1 2 3

2ωμ0

2

 

q

 

 

 

 

 

 

 

 

 

 

 

p r

 

где – эффективная ширина волновода (ТЕ-мода), кото-

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

 

n2n

2 (q2 + p2 )

 

n2n2 (q2 + r 2 )

1

2

 

 

1

3

 

 

w = 2a +

 

 

 

+

 

 

 

.

p(q2n

4 + p

2n4 )

r(q

2n4

+ r 2n4 )

 

 

 

2

1

 

 

3

1

 

Отношение мощности в сердцевине Pcore к полной мощности Р задается

выражением Pcore / P = (2a + J )/ w,

где

 

 

 

 

 

 

12

 

 

 

 

p

 

 

+

 

 

r

 

 

для ТЕ-моды;

 

 

 

 

 

 

 

 

 

 

 

 

 

q2 + p2

q2 + r 2

 

 

 

 

 

J =

n2n

2 p

 

 

 

 

 

n2n2r

 

 

1

 

2

 

 

+

 

1

3

 

для ТМ-моды.

 

q

2n4

+ p2n

4

q

2n4

+ r

2n4

 

2

 

1

 

 

 

 

3

 

 

1

 

2.2.Задания

Всоответствии со своим вариантом (табл. 2.1) произвести расчет планарного трехслойного оптического волновода и проанализировать результаты.

 

 

 

 

 

 

 

 

Таблица 2.1

 

Исходные данные для различных вариантов задания

 

 

 

 

 

 

 

 

 

 

 

 

Длина вол-

Толщина

 

 

 

 

 

 

Вариант

волновода,

n1

n2

n3

 

 

N

 

ны, мкм

мкм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,625

20

1,6

1,55

1

 

0…5

 

 

 

 

 

 

 

 

 

 

2

1,3

52

1,32

1,2

1,18

 

2…7

 

 

 

 

 

 

 

 

 

 

3

0,85

100

1,55

1

1

 

0…3;

10; 22

 

 

 

 

 

 

 

 

 

4

0,68

35

1,78

1,3

1,2

 

3…8

 

 

 

 

 

 

 

 

 

5

1,8

75

1,56

1,4

1,3

 

1; 3; 5; 7; 9; 11

 

 

 

 

 

 

 

 

 

Для этого выполнить следующие задания:

1. Построить распределение электрического поля в поперечном сечении трехслойного планарного волновода для ТЕN- и ТМN-направляемых мод для заданных N.

2.Определить количество направляемых ТЕ- и ТМ-мод в волноводе.

3.Определить эффективную ширину для каждой моды.

4.Определить долю мощности от общей мощности излучения, распространяющегося в сердцевине волновода.

5.Построить зависимости эффективной ширины и доли мощности в сердцевине волновода от номера моды.

2.3.Содержание отчета

Вотчете по практической работе необходимо представить: 1. Цель работы.

2. Краткое описание методики расчета волновода.

3. Графики рассчитанных зависимостей.

4. Определенное количество направляемых мод.

13

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

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

7.Выводы.

3. РАСЧЕТ СВЯЗАННЫХ ВОЛНОВОДОВ

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

Предмет исследования – изучение влияния параметров волноводов и межволноводной среды на коэффициент связи волноводов.

3.1. Общие сведения

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

Рассмотрим некоторые свойства системы из двух связанных планарных волноводов, расположенных на расстоянии d друг от друга (рис. 3.1). Для описания системы двух связанных волноводов могут быть использованы методы теории возмущений. Второй волновод

Рис. 3.1 рассматривается как возмущающий фактор, приводящий к изменению показате-

14

ля преломления первого волновода [4], [5]. Теория возмущений позволяет получить систему из двух уравнений, описывающих амплитуды поля А(z) в двух связанных волноводах:

A1,2 = k1,2 A2,1 exp( i (β2 − β1)z).

z

Здесь индексы 1 и 2 относятся к первому и второму волноводам соответственно; K – коэффициент связи между волноводами; β – константа распро-

странения. При этом K равно

 

 

 

k

 

K1,2 = − i

0

∫∫ R1R2δn2dxdy .

 

1,2

 

В этом выражении k0 – волновое число; R = R(x, y) – поперечные компоненты поля мод; δn – возмущающая добавка к показателю преломления первого волновода.

Практический интерес представляет случай, когда связь между волноводами максимальна. Такая связь возникает при выполнении условия фазового синхронизма – равенства фазовых скоростей мод в двух волноводах. При этом

∆β = β1 – β2 = 0, k1 = k2 = k.

В данном случае система уравнений, описывающая компоненты поля A(z) связанных мод, примет простой вид:

 

A1 = kA ;

 

z

2

 

 

 

A

 

 

2

= kA .

 

 

z

1

 

 

Отсюда можно получить выражения для компонент поля мод связанных волноводов в условиях фазового синхронизма:

A1 = A0 cos ( k z)exp(iβ z);

A2 = − iA0 sin( k z)exp(iβ z).

Здесь А0 – константа, определяемая из граничных условий; β = β1 = β2. В вы-

ражение для А1 входит функция cos(kz), а в выражение для А2 – функция sin(kz). Это означает, что амплитуды поля мод в первом и во втором волноводах изменяются по периодическому закону вдоль направления распространения электромагнитных волн. Причем при увеличении мощности волны в од-

15

ном волноводе происходит уменьшение мощности волны во втором волноводе (рис. 3.2), т. е. происходит периодическая перекачка мощности из одного волновода в другой. Полная перекачка мощности из первого волновода во второй происходит на расстоянии z = = L (рис. 3.2). Параметр L называется длиной связи и в условиях фазового синхронизма определяется из выра-

жения

L = πk/2.

Для двух одинаковых (n2 = n4)

планарных волноводов с ТЕ-модами

Рис. 3.2

при ∆β = 0 и условии (n2 n3)/n2 << 1

выражение для коэффициента связи примет вид

k ≈ − i

 

 

 

 

n

 

 

 

n

2

 

 

 

 

 

2

 

 

exp

 

d ,

 

 

 

n

 

 

 

 

 

 

 

n

2

 

D n

2

 

 

λ

 

 

 

λ

 

3

 

 

 

+ 2

 

 

 

 

 

 

n2

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

где D – толщина волновода.

 

Из данного выражения следует, что коэффици-

ент связи экспоненциально зависит от расстояния между волноводами d, а также является функцией показателей преломления волноводов n2 и среды между ними n3. Изменяя данные параметры, можно изменять коэффициент связи и тем самым управлять процессом перекачки мощности из одного волновода в другой.

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

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

На рис. 3.3, б показан волноводный фильтр, состоящий из кольцевого волноводного резонатора, связанного с основным волноводом. Кольцевой ре-

16

Рис. 3.3

зонатор представляет собой замкнутый волновод, в котором электромагнитная волна может циркулировать по кругу. Если при полном обходе волновода фаза волны оказывается равной 2πm (m = 0, 1, 2 …), то происходит сложение амплитуд – резонанс. В устройстве, показанном на рисунке, в этом случае мощность оптического сигнала в основном волноводе будет уменьшаться за счет ее перекачки в кольцевой волновод. На других (нерезонансных) длинах волн оптический сигнал после прохода по кольцевому резонатору будет возвращаться в волновод без потерь мощности.

3.2.Задания

Всоответствии со своим вариантом (табл. 3.1) произвести расчет свя-

занных оптических волноводов и проанализировать результаты.

 

 

 

 

 

 

Таблица 3.1

 

Исходные данные для различных вариантов задания

 

 

 

 

 

 

 

 

Длина вол-

Толщина

Расстояние

 

 

Доля перека-

 

волново-

между вол-

 

 

Вариант

n2

n3

ченной мощ-

ны

дов D,

новодами d,

 

λ, мкм

 

 

ности P2

 

мкм

мкм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1,3

100

0,18

1,55

1

0,3

 

 

 

 

 

 

 

2

0,68

85

0,09

1,2

1,18

0,5

 

 

 

 

 

 

 

3

0,85

45

0,3

1,65

1,55

0,1

 

 

 

 

 

 

 

4

1,8

96

0,23

1,3

1,2

0,7

 

 

 

 

 

 

 

5

0,53

35

0,16

1,45

1,3

1

 

 

 

 

 

 

 

Для этого выполнить следующие задания:

1. Для пары планарных связанных оптических волноводов построить за-

висимости нормированных мощностей излучения в каждом плече от длины связи волноводов (Р1(z) и Р2(z)).

17

2. Определить три набора оптимальных параметров (z, d) по набору входных данных (λ, D, n2, n3):

а) при полной перекачке оптической мощности во второй волновод;

б) при заданной доле перекачанной мощности P2.

3. Определить максимальное расстояние между волноводами для задан-

ных деления по мощности P2 и набору параметров (λ, D, n2, n3) при z = 2 мм.

3.3.Содержание отчета

Вотчете по практической работе необходимо представить: 1. Цель работы.

2. Краткое описание явления электромагнитной связи в оптических вол-

новодах.

3.Графики рассчитанных зависимостей.

4.Таблицы наборов оптимальных параметров.

5.Значение максимального расстояния между волноводами и описание методики ее определения.

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

7.Описание применения данного физического явления в интегральной

иволоконной оптике.

8.Выводы.

4. РАСЧЕТ ОПТИЧЕСКИХ СВОЙСТВ КОМПОЗИТА

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

Предмет исследования – расчет коэффициентов поглощения композитной среды.

4.1. Общие сведения

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

18

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

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

Коэффициент пропускания Т слоя среды толщиной d определяется как

Т = Евых/Е0.

Здесь Евых – энергия прошедшего через среду излучения; Е0 – энергия пада-

ющего излучения. Для композитной среды, содержащей частицы одного типа с концентрацией N и с известными сечениями поглощения и рассеяния σa и

σs, в приближении однократного рассеяния энергия излучения прошедшего через слой среды строго по оси z (рис. 4.1, а) выражается как

Eвых = E0exp[– σeNd] = E0exp[– ( σa + σs)Nd].

Следовательно, коэффициент пропускания можно представить в виде

Т = exp[– ( σa + σs)Nd].

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

б). Очевидно, что приближение однократного Рис. 4.1 рассеяния справедливо для тонких слоев среды и/или малой концентрации частиц.

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

 

 

+ σ

 

)N

 

T = exp − (σ

a,i

s,i

d .

 

 

i

 

i

 

 

 

 

 

19

Здесь индекс i соответствует группе частиц с одинаковыми характеристиками.

Оптическими характеристиками частиц являются сечение поглощения σa, сечение рассеяния σs и сечение экстинкции σe. Сечения поглощения и рассеяния определяются следующим образом:

σa = Wa/S, σs = Ws/S,

где Wa – скорость поглощения энергии частицей; Ws – скорость рассеяния энергии частицей; S – количество энергии, падающей в единицу времени на единичную площадку в сечении частицы. Единица измерения сечений поглощения и рассеяния – квадратные сантиметры. Сечение экстинкции представляет собой сумму сечений поглощения и рассеяния:

σe = σa + σs.

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

p = εhδE.

Здесь р – вектор дипольного момента; εh – диэлектрическая проницае-

мость среды, в которой находится диполь; Е – вектор электрического поля электромагнитной волны.

Показатель преломления n и диэлектрическая проницаемость ε среды связаны соотношением ε = n2. Для среды с поглощением показатель преломления является комплексной величиной n* = n – ik , где k – коэффициент экс-

тинкции, связанный с коэффициентом поглощения α соотношением

α = 4π kλ .

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

Сферическая наночастица:

 

2

 

εh

r3

 

128π5ε

2r6

 

 

2

 

 

 

 

 

 

 

 

 

σa =

 

 

 

 

 

Im(δ),

σs =

 

h

 

δ

 

,

(4.1)

 

 

 

 

 

4

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

δ=

ε p − εh

,

 

 

 

 

 

 

(4.2)

 

 

 

ε p + 2εh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]