Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
392
Добавлен:
11.03.2016
Размер:
11.85 Mб
Скачать

односторонних отражателей и оптических диодов – устройств, пропускающих свет в одном направлении и блокирующих его в обратном направлении (невзаимное пропускание).

Рассмотрим особенности распространения света через 1D КФК, в котором изменение диэлектрической проницаемости описывается законом

(x)

i

1 a(x)sin2

2 x

.

(7.3.1)

 

(x)

 

 

 

 

 

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

Будем предполагать, что кристаллы немагнитные, т. е. 1. На

рис. 7.4 приведена зависимость коэффициента отражения R от длины волны , рассчитанная в [14] при различных углах падения света

с p-(кривая 1)

и s-(кривая 2)

поляризациями для случая a(x) const ,

(x)

max

min / d x

min (т.е. для линейного профиля измене-

ния шага модуляции). При этом шаг модуляции вдоль направления

распространения света увеличивался от min до max .

 

 

На рис. 7.4, а штриховая кривая соответствует случаю

(x)

const ,

т.е. идеально периодической структуре. В этом случае

(x)

можно

разложить в ряд Фурье:

 

 

 

 

2

 

 

(x)

k eigkx ,

 

 

k

2

 

 

где

g

2

,

0

( i )(1 a / 2)

,

1

1

0

и

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

a(

 

) / 4 . Тогда согласно [15] ширина ФЗЗ равна

 

 

 

 

 

 

 

 

 

 

 

/

 

0 ,

 

 

 

 

 

 

 

 

 

 

 

 

0

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

0

 

 

0 – центральная длина волны ФЗЗ,

 

( max

 

min ) / 2 .

 

 

 

 

 

 

 

291

 

 

 

 

 

 

 

 

Rs ) / 2
I (R p

Согласно рис. 7.4, а в нашем случае для идеально периодической

среды ФЗЗ простирается от

min

0.668 до

max

0.739 мкм. Отме-

 

 

 

тим, что хотя в 1D ФК с чередующимися слоями двух различных веществ каждый отдельный слой изотропен, структура в целом ведет себя как анизотропная среда. S- и p-волны распространяются с разными эффективными фазовыми скоростями и периодическая среда является двулучепреломляющей за счет формы.

Из рис. 7.4 также видно, что при наличии градиента шага модуляции ФЗЗ также появляется. Причем наличие градиента шага модуляции приводит к смещению ФЗЗ и ее существенному расширению. Как и в случае хиральных ФК с градиентным шагом спирали, ширина ФЗЗ зависит от значений min и max , поскольку свет с различными длина-

ми волн претерпевает дифракционное отражение от различных глубин образца, обусловливая расширение ФЗЗ [16, 17]. Однако, как показали численные расчеты [14], в случае наличия градиентного шага модуляции градиент шага также вносит вклад в изменение ширины ФЗЗ.

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

С практической точки зрения важное значение имеет также исследование взаимодействия с рассматриваемой системой неполяризованного света. На рис. 7.5 представлена зависимость интенсивности отра-

жения от длины волны для естественно-

поляризованного падающего света, рассчитанная в [14] при различных углах падения. Видно, что и в этом случае имеется область (от a до б), где имеет место всенаправленное отражение.

Как

показали

расчеты

[14], в

случае если

(x) const , но

a(x)

amax amin

/ d x

amin (т. е.

при изменении амплитуды мо-

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

292

а

б

в

г

д

Рис. 7.4. Зависимость коэффициента отражения R от длины волны , рассчитанная при различных углах падения в случае падения на систему света с p-(кривая 1) и s-(кривая 2)

поляризациями при d

44 мкм; 0 2.25 ; 0

0.000001;

min

0.38

мкм;

max

0.46 мкм; a(x) const

0.5 [14]

 

 

 

 

293

 

Рис. 7.5. Зависимость интенсивности отражения I от длины волны

,

 

рассчитанная при

различных углах падения. Падающий свет

 

 

 

имеет естественную поляризацию:

 

 

 

1

= 0 , 2

= 20 , 3

= 40 , 4

= 60 , 5 – = 80 . Осталь-

 

 

 

ные параметры те же, что и на рис. 6.4 [14]

 

 

Сопоставим теперь спектры отражения КФК двух типов: пер -

вый

– с

(x)

 

max

min / d x

min , а второй – с

(x)

 

min

max / d

x

 

max .

Будем

полагать при

этом,

что

a(x)

const . В первом случае вдоль направления распространения све-

та шаг модуляции линейно увеличивается от значения

min до

max , а

во втором – линейно уменьшается от

max до min . Заметим, что вто-

рой тип совпадает с первым при падении света на данный КФК с обратной стороны.

На рис. 7.6 приведена зависимость коэффициента отражения R от длины волны , рассчитанная в [14] при нормальном падении света для двух рассматриваемых случаев (тип 1 – кривая 1, тип 2 – кривая 2). Видно, что ФЗЗ для этих двух случаев смещены относительно друг друга. В результате имеются области длин волн (от a до b и от c до d, см. рис. 7.6), где имеет место преимущественно одностороннее отражение, и система может работать как широкополосный оптический диод, полностью пропуская свет при его падении на систему с одной стороны и блокируя его при обратном направлении падения света.

Похожий эффект

может иметь место и в случае, когда

(x) const , но a(x)

const . На рис. 7.7 приведены спектры отраже-

ния, рассчитанные в [14] для КФК следующих двух типов: первый –

294

с (x) amax amin / d x amin и второй – с a(x)

amin amax /d x

amax . Видно, что в этом случае ФЗЗ также смещается относительно

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

Рис. 7.6. Зависимость коэффициента отражения R от длины волны , рассчитанная при различных законах изменения (x) :

1 – (x)

[( max

min ) / d]x

min ; 2 – (x) [( min

max ) / d]x

max ;

a(x)

const . Остальные параметры те же, что и на рис. 7.4 [14]

 

Рис. 7.7. Зависимость коэффициента отражения R от длины волны , рассчитанная при различных законах изменения a(x) :

1

a(x) [(amax amin ) / d]x

amin ;

2 a(x) [(amin amax ) / d]x amax ;

(x)

const 0.42 мкм; amin

0.25 ;

amax 0.75 . Остальные параметры

 

те же, что и на рис. 7.4 [14]

 

 

295

Отметим, что так как при нормальном падении света Rs R p , то подобные системы могут работать как оптические диоды и для естест- венно-поляризованного света.

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

7.4. Особенности плотности состояний одномерного разупорядоченного фотонного кристалла

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

Исследования распространения света в экспериментальных образцах ФК (которые всегда в той или иной степени разупорядочены) показали, что в области частот, соответствующих ФЗЗ, в спектрах пропускания света действительно имеется провал, и коэффициент пропускания экспоненциально уменьшается при увеличении толщины образца. Однако длина затухания больше, чем должно быть для идеальной структуры [18–21].

Теоретическое исследование собственных оптических состояний в одномерных разупорядоченных фотонных кристаллах проводилось, как правило, с помощью численных расчетов. При этом было установлено [11], что существует предельно допустимый уровень разупорядочения, до достижения которого вероятность появления собственного состояния в центре ФЗЗ исчезающе мала. Порог достигается, когда относительная флуктуация оптических длин периодов структуры соответствует корню квадратному из одной трети относительной ширины ФЗЗ. На основании численных расчетов в [22] даже была получена

296

0 )ex

эмпирическая формула плотности оптических состояний одномерных фотонных кристаллов ( ) , учитывающая влияние беспорядка. Ана-

литическая теория плотности состояний в одномерных разупорядоченных ФК была построена в [23].

Ограничимся рассмотрением нормального распространения линейно поляризованного света E(r,t) E(z)cos( t в одномерном

ФК, представляющем собой последовательность пар слоев A и B одинаковой толщины d / 2 , где d – период ФК (рис. 7.8). Показатели преломления слоев A и B запишем в виде

 

nA,B (z) n0 (1 )(1

n P(s)) ,

(7.4.1)

где

nA (s) nB (s) / nA (s) nB (s)

– диэлектрический

контраст,

описывающий модуляцию показателя преломления, не зависящий от

номера периода s ;

n – параметр беспорядка, описывающий флуктуа-

ции показателя преломления; n0

nA (s) nB (s)

/ 2

– среднее значе-

ние показателя преломления. Случайная величина

P(s) – с нулевым

средним P(s) 0

предполагается некоррелированной с единичной

дисперсией P(s)P(s )

ss . Отметим,

что функция P(s) остается

постоянной внутри каждого из периодов ФК.

 

 

Согласно (7.1.1а) для немагнитных слоев (

1) распространение

электромагнитной волны частоты

вдоль оси структуры z описыва-

ется волновым уравнением

 

 

 

 

 

 

 

d 2 E

 

n2 (z)

2

E ,

 

(7.4.2)

 

 

dz2

c2

 

 

 

 

 

 

 

 

здесь n(z) – профиль показателя преломления вдоль оси структуры

(рис. 7.8). Учитывая, что показатель преломления внутри каждого слоя остается постоянными, следуя [23], решение уравнения (7.4.2) представим в виде

E(zs

)

E(zs ) cos(kA )

E (zs )sin(kA

) / kA ,

 

 

 

 

(7.4.3)

E (zs

)

E(zs )kA sin(kA

) E (zs ) cos(kA )

 

 

297

 

 

для слоев А и

E(zs

d / 2

)

E(zs

d / 2) cos(kB )

 

 

 

 

 

 

E (zs

d / 2)sin(kB

) / kB ,

 

 

 

 

 

 

(7.4.4)

E (zs

d / 2

)

E(zs

d / 2)kB sin(kB

)

 

 

 

 

 

E (zs

d / 2) cos(kB

)

для слоев типа

B.

Здесь

zs sd , 0

d / 2 , kA,B

nA,B (s) / c , а

E dE / dz – производная электрического поля.

d

 

 

 

d/2

 

d/2

Рис. 7.8. Модель одномерного разупорядоченного фотонного кристалла [23]

Используя (7.4.3) и (7.4.4), можно получить матрицу переноса Ms , связывающую электрическое поле и его производную в точках zs

и zs 1

 

 

 

 

E

 

 

 

E

 

 

 

 

 

 

 

 

 

c

E

Ms

c

E

.

 

 

 

 

(7.4.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s 1

 

 

 

 

s

 

 

 

 

 

 

В нашем случае матрица переноса имеет вид [23]

 

 

 

 

 

CACB

nA (s)

SASB

1

 

CASB

1

CB SA

 

 

 

nB (s)

 

 

 

 

nB

 

 

 

M

s

 

 

 

nA (s)

 

(s)

 

,

(7.4.6)

 

 

 

 

 

 

 

 

 

 

nB (s)

 

 

 

nA (s)CASB nB (s)CB SA

 

CACB

SASB

 

 

 

 

nA (s)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

298

 

 

 

 

 

 

 

 

 

где

CA,B

cos

nA,B (s)d

, SA,B

sin

nA,B

(s)d

 

 

 

.

2c

 

 

 

 

 

 

2c

Используя (7.4.5) и (7.4.6), можно получить дисперсионное выражение для идеальной периодической структуры в виде (7.2.1):

cos(Kd )

cos

nAd

 

cos

nB d

 

 

 

 

2c

 

2c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

nA

 

nB

sin

 

nAd

 

sin

 

nB d

,

(7.4.7)

2

 

nB

nA

 

2c

 

 

2c

 

 

 

 

 

 

 

 

 

 

здесь K – блоховский волновой вектор. Согласно (7.4.7) в спектре такой структуры имеется ФЗЗ, частота центра которой (рис. 7.9).

0

c

.

(7.4.8)

n d

0

 

 

При

этом относительная

ширина

ФЗЗ

составляет

/

0

4 , а

 

 

 

 

затухание света на одном периоде

Im(K0d)

d / 0 ln nA / nB 2 (где K0 ,

0 – соответственно мнимая часть

блоховского волнового вектора

и

длина затухания на частоте

0 )

при 1 на частоте центра ФЗЗ

прямо пропорционально диэлектрическому контрасту [11].

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

На рис. 7.10 представлены плотности состояний, а также час-

Рис. 7.9. Зависимость блоховского волнового вектора K от частоты для идеальной структуры с g =

= n0 = 0.025 [11]

299

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

 

nA,B 2 g

2 P ,

где величина

P принимает хаотическое значение в интервале от –0.5

до 0.5; g

модуляция показателя

преломления; – относительная

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

A

1

1

,

(7.4.9)

n f

M nl

где n f и nl – показатели преломления полубесконечных сред, ограни-

чивающих структуру, M – матрица переноса через структуру, A – некоторая константа.

Для сравнения на рис 7.10, а приведены результаты расчетов с использованием (7.4.9) для периодической структуры конечных размеров. Отметим, что для такой структуры спектр собственных мод будет дискретным, а вследствие утечек света через границы структуры время

жизни собственных состояний

1/ Im i будет конечным, т. е. собст-

венные частоты

i

будут иметь ненулевую мнимую часть.

 

 

 

Расчеты показали, что в спектре пропускания света периодической структуры с модуляцией показателя преломления g = 0.025 и толщиной L = 200d имеется провал, а коэффициент пропускания в минимуме равен 2 10–4.

Рис. 7.10. Собственные частоты идеальной и разупорядоченной структур:

а – слева:

плотность состояний; справа: частоты

Re i и времена жизни

1/ Im i

для идеальной структуры ( 0) . бд – слева: плотность состояний,

усредненная по ансамблю 3000 структур с различными

; справа: частоты Re i

и времена жизни

1/ Im i собственных мод, полученные путем решения (7.4.9)

для 1000 разупорядоченных структур с теми же . Толщина структуры L = 200d, модуляция показателя преломления g = 0.025 [11]

300

Соседние файлы в папке Наноэлектроника лит-ра