В последнее десятилетие развитие микроэлектроники затормозилось, поскольку уже практически достигнуты ограничения по быстродействию стандартных полупроводниковых устройств. Все большее число исследований посвящается разработке альтернативных полупроводниковой электронике областей - это спинтроника, микроэлектроника со сверхпроводящими элементами, фотоника и некоторые другие.
Новый принцип передачи и обработки информации с помощью светового, а не электрического сигнала может ускорить наступление нового этапа информационного века.
От простых кристаллов к фотонным
Основой электронных устройств будущего могут стать фотонные кристаллы - это синтетические упорядоченные материалы, в которых диэлектрическая проницаемость периодически меняется внутри структуры. В кристаллической решетке традиционного полупроводника регулярность, периодичность расположения атомов приводит к образованию так называемой зонной энергетической структуры - с разрешенными и запрещенными зонами. Электрон, энергия которого попадает в разрешенную зону, может передвигаться по кристаллу, а электрон с энергией в запрещенной зоне оказывается «запертым».
По аналогии с обычным кристаллом возникла идея кристалла фотонного. В нем периодичность диэлектрической проницаемости обуславливает возникновение фотонных зон, в частности, запрещенной, в пределах которой распространение света с определенной длиной волны подавлено. То есть, будучи прозрачными для широкого спектра электромагнитного излучения, фотонные кристаллы не пропускают свет с выделенной длиной волны (равной удвоенному периоду структуры по длине оптического пути).
Фотонные кристаллы могут иметь различную размерность. Одномерные (1D) кристаллы представляют собой многослойную структуру из чередующихся слоев с разными показателями преломления. Двумерные фотонные кристаллы (2D) можно представить в виде периодической структуры из стержней с разной диэлектрической проницаемостью. Первые синтетические прообразы фотонных кристаллов были трехмерными и созданы еще в начале 1990-х годов сотрудниками исследовательского центра Bell Labs (США). Для получения периодической решетки в диэлектрическом материале американские ученые высверливали цилиндрические отверстия таким образом, чтобы получить трехмерную сеть пустот. Для того, чтобы материал стал фотонным кристаллом, его диэлектрическая проницаемость была модулирована с периодом в 1 сантиметр во всех трех измерениях.
Природными аналогами фотонных кристаллов являются перламутровые покрытия раковин (1D), усики морской мыши, многощетинкового червя (2D), крылья африканской бабочки парусника и полудрагоценные камни, например, опал (3D).
Но и сегодня, даже с помощью самых современных и дорогостоящих методов электронной литографии и анизотропного ионного травления, с трудом удается изготовить бездефектные трехмерные фотонные кристаллы с толщиной более 10 структурных ячеек.
Фотонные кристаллы должны найти широкое применение в фотонных интегральных технологиях, которые в перспективе заменят электрические интегральные схемы в компьютерах. При передаче информации с использованием фотонов вместо электронов резко сократится энергопотребление, увеличатся тактовые частоты и скорость передачи информации.
Фотонный кристалл из оксида титана
Оксид титана TiO 2 обладает набором уникальных характеристик, таких как высокий показатель преломления, химическая стабильность и низкая токсичность, что делает его наиболее перспективным материалом для создания одномерных фотонных кристаллов. Если рассматривать фотонные кристаллы для солнечных батарей, то здесь оксид титана выигрывает из-за своих полупроводниковых свойств. Ранее было продемонстрировано увеличение КПД солнечных элементов при использовании слоя полупроводника с периодической структурой фотонного кристалла, в том числе фотонных кристаллов из оксида титана.
Но пока применение фотонных кристаллов на основе диоксида титана ограничивается отсутствием воспроизводимой и недорогой технологии их создания.
Сотрудники химического факультета и факультета наук о материалах МГУ - Нина Саполетова, Сергей Кушнир и Кирилл Напольский - усовершенствовали синтез одномерных фотонных кристаллов на основе пористых пленок оксида титана.
«Анодирование (электрохимическое окисление) вентильных металлов, в том числе алюминия и титана, является эффективным методом получения пористых оксидных пленок с каналами нанометрового размера», - пояснил руководитель группы электрохимического наноструктурирования, кандидат химических наук Кирилл Напольский.
Анодирование обычно проводят в двухэлектродной электрохимической ячейке. В раствор электролита опускают две металлические пластины - катод и анод, и подают электрическое напряжение. На катоде выделяется водород, на аноде происходит электрохимическое окисление металла. Если периодически менять прикладываемое к ячейке напряжение, то на аноде формируется пористая пленка с заданной по толщине пористостью.
Эффективный показатель преломления будет модулироваться, если диаметр пор будет периодически меняться внутри структуры. Разработанные ранее методики анодирования титана не позволяли получать материалы с высокой степенью периодичности структуры. Химики из МГУ разработали новый способ анодирования металла с модуляцией напряжения в зависимости от заряда анодирования, который позволяет с высокой точностью создавать пористые анодные оксиды металлов. Возможности новой методики химики продемонстрировали на примере одномерных фотонных кристаллов из анодного оксида титана.
В результате изменения напряжения анодирования по синусоидальному закону в диапазоне 40–60 Вольт ученые получили нанотрубки анодного оксида титана с постоянным внешним диаметром и периодически изменяющимся внутренним диаметром (см. рисунок).
«Применяемые ранее методики анодирования не позволяли получать материалы с высокой степенью периодичности структуры. Мы разработали новую методику, ключевым составляющим которой является in situ (непосредственно во время синтеза) измерение заряда анодирования, что позволяет с высокой точность контролировать толщину слоев с различной пористостью в формируемой оксидной пленке», - пояснил один из авторов работы, кандидат химических наук Сергей Кушнир.
Разработанная методика упростит создание новых материалов с модулированной структурой на основе анодных оксидов металлов. «Если в качестве практического использования методики рассматривать применение в солнечных батареях фотонных кристаллов из анодного оксида титана, то еще предстоит провести систематическое исследование влияния структурных параметров таких фотонных кристаллов на эффективность преобразования света в солнечных батареях», - уточнил Сергей Кушнир.
Показано, что в зависимости от полярности включения фотодиодов в состав резонатора происходит частотный сдвиг отклика вверх или вниз по частоте при увеличении освещенности. Предложено использовать систему связанных кольцевых резонаторов для увеличения чувствительности исследуемых резонаторов к величине освещенности. Продемонстрировано, что для фиксированного расстояния между связанными резонаторами происходит частотное расщепление отклика системы на четную (яркую) и нечетную (темную) моды при помощи света. Мы уверены, что предложенный метод создания перестраиваемых кольцевых резонаторов позволит создать новый класс метаматериалов, управляемых светом.
Работа поддержана Министерством образования Российской Федерации (соглашения № 14.В37.21.1176 и № 14.В37.21.1283), Фондом «Династия», Фондом РФФИ (проект № 13-02-00411), стипендией Президента Российской Федерации молодым ученым и аспирантам 2012.
Литература
1. Linden S., Enkrich C., Wegener M., Zhou J., Koschny T., Soukoulis C.M. Magnetic Response of Metamaterials at 100 Terahertz // Science. - 2004. - V. 306. - P. 1351-1353.
2. Shelby R., Smith D.R. and Schultz S. Experimental Verification of a Negative Index of Refraction // Science. - 2001. - V. 292. - P. 77-79.
3. Gansel J.K., Thiel M., Rill M.S., Decker M., Bade K., Saile V., von Freymann G., Linden S., Wegener M. Gold Helix Photonic Metamaterial as Broadband Circular Polarizer // Science. - 2009. - V. 325. - P. 15131515.
4. Belov P.A., Hao Y. Subwavelength imaging at optical frequencies using a transmission device formed by a periodic layered metal-dielectric structure operating in the canalization regime // Physical Review B. - 2006. - V. 73. - P. 113110.
5. Leonhardti U. Optical conformal mapping // Science. - 2006. - V. 312. - P. 1777-1780.
6. Кившарь Ю.С., Орлов А.А. Перестраиваемые и нелинейные метаматериалы // Научно-технический вестник информационных технологий, механики и оптики. - 2012. - № 3 (79). - C. 1-10.
7. Shadrivov I.V., Morrison S.K. and Kivshar Yu.S. Tunable split-ring resonators for nonlinear negative-index metamaterials // Opt. Express. - 2006. - V. 14. - P. 9344-9349.
8. Kapitanova P.V., Maslovski S.I., Shadrivov I.V., Voroshilov P.M., Filonov D.S., Belov P.A. and Kivshar Y.S. Controlling split-ring resonators with light // Applied Physics Letters. - V. 99. - P. 251914 (1-3).
9. Marques R., Martin F. and Sorolla M. Metamaterials with Negative Parameters: Theory, Design and Microwave Applications. - NJ: Wiley&Sons, Inc., Hoboken, 2008. - 315 p.
Капитонова Полина Вячеславовна - Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, кандидат технических наук, научный сотрудник, [email protected], [email protected]
Белов Павел Александрович - Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, доктор физ.-мат. наук, главный научный сотрудник, [email protected]
АНАЛИЗ ЗОННОЙ СТРУКТУРЫ ФОТОННОГО КРИСТАЛЛА С КРАТНЫМИ ОПТИЧЕСКИМИ ДЛИНАМИ СЛОЕВ ДЛЯ ТЕРАГЕРЦОВОГО ДИАПАЗОНА
А.Х. Денисултанов, М.К. Ходзицкий
Из дисперсионного уравнения для бесконечного фотонного кристалла выведены формулы для точного расчета границ запрещенных зон, ширины запрещенных зон и точного положения центров запрещенных зон фотонных кристаллов с кратными оптическими длинами слоев в двухслойной ячейке для терагерцового диапазона частот от 0,1 до 1 ТГц. Формулы проверены при численном моделировании фотонных кристаллов методом матриц передачи и методом конечных разностей временной области для первой, второй и третьей кратностей оптических длин в двухслойной ячейке фотонного кристалла. Формулы для второй кратности подтверждены экспериментально. Ключевые слова: фотонный кристалл, запрещенная зона, граничные частоты, кратные оптические длины, матрица передачи, метаматериал.
Введение
В последние годы исследование искусственных сред с необычными свойствами («метаматериа-лов») привлекает интерес достаточно большого круга ученых и инженеров, что обусловливается перспективным использованием этих сред в промышленной и военной индустрии при разработке новых типов фильтров, фазосдвигателей, суперлинз, маскирующих покрытий и т.д. . Одним из видов мета-материалов является фотонный кристалл, который представляет собой слоистую структуру с периодиче-
ски изменяющимся показателем преломления . Фотонные кристаллы (ФК) активно используются в лазерных технологиях, средствах коммуникации, фильтрации, благодаря таким уникальным свойствам, как наличие зонной структуры в спектре, сверхразрешение, эффект суперпризмы и т.д. . Особый интерес проявляется к исследованию фотонных кристаллов в терагерцовом (ТГц) диапазоне для спектроскопических, томографических исследований новых типов материалов и биообъектов . Исследователями уже разработаны двумерные и трехмерные ФК для ТГц диапазона частот и изучены их характеристики , но, к сожалению, на данный момент нет точных формул для расчета характеристик зонной структуры фотонного кристалла, таких как ширина запрещенной зоны, центр запрещенной зоны, границы запрещенной зоны . Целью настоящей работы является получение формул для расчета характеристик одномерного фотонного кристалла для первой, второй и третьей кратностей оптических длин в двухслойной ячейке ФК и проверка этих формул с помощью численного моделирования методом матриц передачи и методом конечных разностей во временной области, а также эксперимента в ТГц диапазоне частот.
Аналитическое и численное моделирование
Рассмотрим бесконечный фотонный кристалл с показателями преломления слоев в двухслойной ячейке п1 и п2 и толщинами слоев й1 и й2 соответственно. Данная структура возбуждается линейно-поляризованной поперечной электрической волной (ТЕ-волной). Волновой вектор к направлен перпендикулярно слоям ФК (рис. 1). Дисперсионное уравнение для такого ФК, полученное с использованием теоремы Флоке и условия непрерывности тангенциальных компонент поля на границе слоев, имеет следующий вид :
С08[кв(йх + й2)] = со8[кг й^]х со$[к2 й2]-0,5)
с бт[кг ё1] х бт[кг й2
где кв - блоховское волновое число; к^ =
ли преломления; й1, й2 - толщины слоев.
2 л х / х п1
; / - частота; пг, п2 - показате-
Рис. 1. Рассматриваемая слоисто-периодическая структура
Л. и Л 1 ! I х. ] л!/ л Пил! л «
и " и | Г ¡4 1 ! 1) 1 1 N V и | 1 У " 11
0,1 0,2 0,3 0,4 0,5 0,6
Частота/ ТГц
Рис. 2. Частотная дисперсия комплексного блоховского волнового числа
Дисперсия комплексного блоховского волнового числа, полученная с использованием уравнения (1), показана на рис. 2. Как видно из рис. 2, на границах запрещенных зон аргумент косинуса кв (й1 + й2) будет принимать значения либо 0, либо п . Следовательно, исходя из этого условия, можно рассчи-
тать значения граничных частот, ширины запрещенных зон и центры запрещенных зон фотонного кристалла. Однако для фотонного кристалла с некратными оптическими длинами слоев внутри двухслойной ячейки данные формулы могут быть получены только в неявном виде. Для получения формул в явном виде нужно использовать кратные оптические длины: пхёх = п2ё2; пхёх = 2хп2ё2; пхёх = 3хп2ё2... . В работе были рассмотрены формулы для 1-й, 2-й и 3-й кратности.
Для фотонного кристалла первой кратности (пхёх = п2ё2) формулы граничных частот, ширины
запрещенной зоны и центра запрещенной зоны имеют следующий вид:
(/п 1 Л (/п «и 1 Л
0,256-1,5 . „ агссо81---I + 2лт
а/ = /1 -/2; /33 = /+/2-; /рз =
/ 2а; /2 = я(т +1)
0,256-1, 5 . „, 1Ч -агссо81 ----- | + 2л(т +1)
где /1 и /2 - низкочастотная и высокочастотная границы запрещенной зоны соответственно; А/ - ширина запрещенной зоны; /зз - центр запрещенной зоны; с - скорость света; / - центр разрешенной
о пх п2 зоны 6 = - +-;
Для ФК с параметрами слоев пх = 2,9; п2 = 1,445; ёх = 540 мкм; ё2 = 1084 мкм для второй запрещенной зоны в диапазоне 0,1-1 ТГц имеют место следующие параметры зонной структуры: /1 = 0,1332 ТГц; /2 = 0,1541 ТГц; А/ = 0,0209 ТГц; /зз = 0,1437 ТГц.
Для ФК, оптические длины слоев которого связаны равенством пхёх = 2п2ё2, получены следующие формулы для параметров зонной структуры:
4+в+У в2-4 6 + 3в-4в2 -4
4 + в-V в2 - 4 6 + 3в + ^в2 - 4
2 + в -V в2 - 4
2ят х с агссоБ
В-#^4 2 + в + 4 в2 - 4
В-#^4 2 + в + л/в2 - 4
4 + в-Vв2 -4 6 + 3в + 4в2 - 4
4 + в + Ув2 - 4 6 + 3в-4в2 -4
где (/1 и /11), (/2 и /21), (/3 и /31), (/4 и /41) - низкочастотная и высокочастотная границы запрещен-
ных зон с номерами (4т+1), (4т+2), (4т+3), (4т+4) соответственно; с - скорость света; Р= - + -;
т = 0,1,2,.... Ширина запрещенной зоны рассчитывается как А/ = /-/х; центр запрещенной зоны
, / + /х. й /зз = ^ ; /рз - центр разрешенной зоны.
Для ФК с параметрами пх = 2,9; п2 = 1,445; ёх = 540 мкм; ё2 = 541,87 мкм для второй запрещенной зоны в диапазоне 0,1-1 ТГц имеем
/2 = 0,116 ТГц; /2х = 0,14 ТГц; А/ = 0,024 ТГц; /зз = 0,128 ТГц.
Для фотонного кристалла, оптические длины которого связаны равенством пхёх = 3п2ё2, получены следующие формулы для параметров зонной структуры:
1 -0,5ß + ^/2,25ß2 -ß-7 3 + 2,5ß-^/ 2,25ß2-ß-7
1 -0,5ß-^2,25ß2 -ß-7 3 + 2,5ß + V 2,25ß2-ß-7
1 -0,5ß-J2,25ß2 -ß-7 3 + 2,5ß + yl2,25ß2 - ß - 7
1 - 0,5ß + 72,25ß2 - ß - 7 3 + 2,5ß-sj2,25ß2 -ß-7
где (/1 и /11), (/2 и /2), (/3 и /) - низкочастотная и высокочастотная границы запрещенных зон с
номерами (3т+1), (3т+2), (3т+3) соответственно; с - скорость света; р = - + -; т = 0,1,2,.... Ширина
запрещенной зоны рассчитывается как Д/ = / - /1; центр запрещенной зоны /зз =
разрешенной зоны.
Для ФК с параметрами п1 = 2,9; п2 = 1,445; = 540 мкм; й2 = 361,24 мкм для второй запрещенной зоны в диапазоне 0,1-1 ТГц имеем
/2 = 0,1283 ТГц; = 0,1591 ТГц; Д/ = 0,0308 ТГц; /зз = 0,1437 ТГц.
Для моделирования ФК конечной длины нужно использовать метод матриц передачи , который позволяет рассчитать значение электромагнитного поля волны, проходящей через фотонный кристалл, в произвольной точке 2 слоя. Матрица передачи для одного слоя имеет следующий вид:
cos(k0 x n x p x sin(k0
: z x cos 0) x n x z x cos 0)
(-i / p) x sin(k0 x n x z x cos 0)
где k0 = -; p = - cos 0 ; n = ; z - координата на оси Oz; 0 - угол падения волны на первый слой.
Используя метод матриц передачи, в математическом пакете MATLAB была построена зонная структура фотонного кристалла для оптических длин слоев в двухслойной ячейке 1-й, 2-й и 3-й кратно-стей), в ТГц диапазоне частот (для 0=0) с 10 элементарными ячейками с параметрами слоев, указанными выше (рис. 3).
Как видно из рис. 3, в спектре пропускания ФК 1-й, 2-й и 3-й кратности выпадают запрещенные зоны, кратные двум, трем, четырем соответственно, по сравнению с зонной структурой ФК с некратными оптическими длинами слоев внутри элементарной ячейки. Для всех трех случаев кратности относительная погрешность вычислений параметров зонной структуры конечного ФК не превышает 1% по сравнению с формулами для бесконечного ФК (ширина запрещенной зоны рассчитывалась на уровне 0,5 коэффициента пропускания для конечного ФК).
Также структура одномерного ФК была рассчитана методом конечных разностей во временной области с помощью программного пакета трехмерного моделирования CST Microwave Studio (рис. 4). Видно такое же поведение зонной структуры конечного ФК, что и для спектров пропускания, полученного методом матриц передачи. Относительная погрешность вычислений параметров зонной структуры конечного ФК в данном пакете моделирования не превышает 3% по сравнению с формулами для бесконечного ФК.
Цж.М"."ш ЩШШ Ш Щ"ДЦ Щ
пШшиЩШ) щщм
пхёх=3п2ё2 Частота / ТГц
Рис. 3. Зонная структура фотонного кристалла для трех кратностей, оптических длин слоев в двухслойной ячейке в ТГц диапазоне частот (цифры указывают номер запрещенной зоны, стрелки - выпадающие
запрещенные зоны)
Я -е -е т о
пхёх=2п2ё2 -ДА/ ут1
пхёх=3п2ё2 Частота, ТГц
Рис. 4. Трехмерная модель ФК в ОЭТ (а) и коэффициент пропускания ФК для трех кратностей (б)
Экспериментальная часть
Случай 2-й кратности был проверен экспериментально методом непрерывной ТГц спектроскопии в диапазоне 0,1-1 ТГц . Был использован метод смешения частот инфракрасного излучения на фото-проводящей (ФП) антенне для генерации ТГц излучения. Вторая ФП антенна была использована в качестве приемника. Между излучающей и принимающей ФП антенной устанавливался собранный ФК (рис. 5).
Исследованный фотонный кристалл имеет следующие параметры: количество бислойных ячеек -3; показатели преломления слоев - пх = 2,9 и п2 = 1,445 ; толщины слоев - ёх = 540 мкм и ё2 = 520 мкм (ё2 на 21 мкм меньше, чем для случая идеальной 2-й кратности). На рис. 5 показано сравнение экспериментального и теоретического спектра для 4 и 5 запрещенных зон. Как видно из экспериментального графика, так же как и для моделирования, наблюдается выпадение запрещенной зоны, кратной трем, по сравнению с зонной структурой ФК с некратными оптическими длинами слоев внутри элементарной ячейки. Небольшое несоответствие положения центров запрещенных зон в экспериментальном и теоре-
тическом спектре связано с отличием толщины слоев тефлона в эксперименте от идеальной 2-й кратности.
1,0 0,9 0,8 0,7 0,6 0,5 0,4 0,3
0,3 0,35 0,4 0,45 0,5 Частота, ТГц
Эксперимент
Моделирование
Рис. 5. Фотография установки, фотография макета фотонного кристалла (а) и сравнительный график экспериментального и теоретического коэффициента пропускания ФК с тремя элементарными
ячейками (б)
Заключение
Таким образом, были получены точные формулы для расчета параметров зонной структуры (ширина запрещенной зоны, границы запрещенной зоны и центр запрещенной зоны) одномерных фотонных кристаллов с кратными оптическими длинами слоев внутри двухслойной элементарной ячейки для случая TE-волны с волновым вектором, перпендикулярным плоскостям слоев фотонного кристалла. Было продемонстрировано для фотонного кристалла 1-й, 2-й и 3-й кратности исчезновение запрещенных зон, кратным двум, трем, четырем соответственно, по сравнению с зонной структурой фотонных кристаллов с некратными оптическими длинами слоев внутри элементарной ячейки. Формулы для 1-й, 2-й и 3-й кратностей были проверены с помощью метода матриц передачи и трехмерного численного моделирования методом конечных разностей во временной области. Случай 2-й кратности был проверен в эксперименте в ТГц диапазоне частот от 0,1 до 1 ТГц. Полученные формулы могут быть использованы для разработки широкополосных фильтров на основе фотонных кристаллов для промышленного, военного и медицинского применения без необходимости моделирования зонной структуры фотонного кристалла в различных математических пакетах.
Работа была частично поддержана грантом № 14.132.21.1421 в рамках реализации Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2013 гг.
Литература
1. Вендик И.Б., Вендик О.Г. Метаматериалы и их применение в технике сверхвысоких частот (Обзор) // Журнал технической физики. - СПбГЭТУ «ЛЭТИ». - 2013. - Т. 83. - Вып. 1. - С. 3-26.
2. Возианова А.В., Ходзицкий М.К. Маскирующее покрытие на основе спиральных резонаторов // Научно-технический вестник информационных технологий, механики и оптики. - 2012. - № 4 (80). -С. 28-34.
3. Терехов Ю.Е., Ходзицкий М.К., Белокопытов Г.В. Характеристики метапленок для терагерцового диапазона частот при масштабировании геометрических параметров // Научно-технический вестник информационных технологий, механики и оптики. - 2013. - № 1 (83). - С. 55-60.
4. Yablonovitch E. Inhibited Spontaneous Emission in Solid-State Physics and Electronics // Physical Review Letters. - 1987. - V. 58. - № 20. - P. 2059-2062.
5. Figotin A., Kuchment P. Band-Gap Structure of Spectra of Periodic Dielectric and Acoustic Media. II. Two-Dimensional Photonic Crystals // SIAM Journal on Applied Mathematics. - 1996. - V. 56. - № 6. - P. 1561-1620.
6. Smolyaninov Igor I., Davis Christopher C. Super-resolution optical microscopy based on photonic crystal materials // Physical review B. - 2005. - V. 72. - P. 085442.
7. Kosaka Hideo, Kawashima Takayuki, Tomita Akihisa. Superprism phenomena in photonic crystals // Physical review B. - 1998. - V. 58. - № 16. - P. 10096-10099.
8. Kurt Hamza, Erim Muhammed Necip, Erim Nur. Various photonic crystal bio-sensor configurations based on optical surface modes // Department of Electrical and Electronics Engineering. - 2012. - V. 165. - № 1. - P. 68-75.
9. Ozbay E., Michel E., Tuttle G., Biswas R., Sigalas M., and Ho K.M. Micromachined millimeter-wave photonic band-gap crystals // Appl. Phys. Lett. - 1994. - V. 64. - № 16. - P. 2059-2061.
10. Jin C., Cheng B., Li Z., Zhang D., Li L.M., Zhang Z.Q. Two dimensional metallic photonic crystal in the THz range // Opt. Commun. - 1999. - V. 166. - № 9. - P. 9-13.
11. Nusinsky Inna and Hardy Amos A. Band-gap analysis of one-dimensional photonic crystals and conditions for gap closing // Physical review B. - 2006. - V. 73. - P. 125104.
12. Басс Ф.Г., Булгаков А.А., Тетервов А.П. Высокочастотные свойства полупроводников со сверхрешетками. - М.: Наука. Гл. ред. физ.-мат. лит., 1989. - 288 с.
13. Борн М., Вольф Э. Основы оптики. - М.: Наука. Гл. ред. физ.-мат. лит., 1973. - 733 c.
14. Gregory I.S., Tribe W.R., Baker C. Continuous-wave terahertz system with a 60 dB dynamic range // Applied Phisics Letters. - 2005. - V. 86. - P. 204104.
Денисултанов Алауди Хожбаудиевич
Ходзицкий Михаил Константинович
Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, студент, [email protected]
Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кандидат физ.-мат. наук, ассистент, [email protected]
Фотонные кристаллы по характеру изменения коэффициента преломления можно разделить на три основных класса:
1. Одномерные, в которых коэффициент преломления периодически изменяется в одном пространственном направлении как показано на рисунке 2. На этом рисунке символом Л обозначен период изменения коэффициента преломления, и - показатели преломления двух материалов (но в общем случае может присутствовать любое число материалов). Такие фотонные кристаллы состоят из параллельных друг другу слоев различных материалов с разными коэффициентами преломления и могут проявлять свои свойства в одном пространственном направлении, перпендикулярном слоям.
Рисунок 1 - Схематическое представление одномерного фотонного кристалла
2. Двухмерные, в которых коэффициент преломления периодически изменяется в двух пространственных направлениях как показано на рисунке 2. На этом рисунке фотонный кристалл создан прямоугольными областями с коэффициентом преломления, которые находятся в среде с коэффициентом преломления. При этом, области с коэффициентом преломления упорядочены в двумерной кубической решетке. Такие фотонные кристаллы могут проявлять свои свойства в двух пространственных направлениях, и форма областей с коэффициентом преломления не ограничивается прямоугольниками, как на рисунке, а может быть любой (окружности, эллипсы, произвольная и т. д.). Кристаллическая решётка, в которой упорядочены эти области, также может быть другой, а не только кубической, как на приведённом рисунке.
Рисунок - 2 Схематическое представление двумерного фотонного кристалла
3. Трёхмерные, в которых коэффициент преломления периодически изменяется в трёх пространственных направлениях. Такие фотонные кристаллы могут проявлять свои свойства в трёх пространственных направлениях, и можно их представить как массив объёмных областей (сфер, кубов и т. д.), упорядоченных в трёхмерной кристаллической решётке.
Как и электрические среды в зависимости от ширины запрещённых и разрешённых зон, фотонные кристаллы можно разделить на проводники - способные проводить свет на большие расстояния с малыми потерями, диэлектрики - практически идеальные зеркала, полупроводники - вещества способные, например, выборочно отражать фотоны определённой длины волны и сверхпроводники, в которых благодаря коллективным явлениям фотоны способны распространяться практически на неограниченные расстояния.
Также различают резонансные и нерезонансные фотонные кристаллы. Резонансные фотонные кристаллы отличаются от нерезонансных тем, что в них используются материалы, у которых диэлектрическая проницаемость (или коэффициент преломления) как функция частоты имеет полюс на некоторой резонансной частоте.
Любая неоднородность в фотонном кристалле называются дефектом фотонного кристалла. В таких областях часто сосредотачивается электромагнитное поле, что используется в микрорезонаторах и волноводах, построенных на основе фотонных кристаллов.
Как и электрические среды в зависимости от ширины запрещённых и разрешённых зон, фотонные кристаллы можно разделить на проводники - способные проводить свет на большие расстояния с малыми потерями, диэлектрики - практически идеальные зеркала, полупроводники - вещества способные, например, выборочно отражать фотоны определённой длины волны и сверхпроводники, в которых благодаря коллективным явлениям фотоны способны распространяться практически на неограниченные расстояния. Также различают резонансные и нерезонансные фотонные кристаллы. Резонансные фотонные кристаллы отличаются от нерезонансных тем, что в них используются материалы, у которых диэлектрическая проницаемость (или коэффициент преломления) как функция частоты имеет полюс на некоторой резонансной частоте.
Любая неоднородность в фотонном кристалле называются дефектом фотонного кристалла. В таких областях часто сосредотачивается электромагнитное поле, что используется в микрорезонаторах и волноводах, построенных на основе фотонных кристаллов. Существует ряд аналогий при описании распространения электромагнитных волн в фотонных кристаллах и электронных свойств кристаллов. Приведем некоторые из них.
1. Состояние электрона внутри кристалла (закон движения) задается решением уравнения Шрлдингера, распространение света в фотонном кристалле подчиняется волновому уравнению, являющемуся следствием уравнений Максвелла:
- 2. Состояние электрона описывается скалярной волновой функцией ш(r,t), состояние электромагнитной волны описывается векторными полями - напряженностью магнитной или электрической компонент, H (r,t) или E(r,t).
- 3. Волновая функция электрона ш(r,t) может быть разложена в ряд по собственным состояниям шE(r), каждому из которых соответствует собственная энергия E. Напряженность электромагнитного поля H(r,t) может быть представлена суперпозицией монохроматических компонент (мод) электромагнитного поля Hщ(r), каждой из которой соответствует собственное значение - частота моды щ:
4. Атомный потенциал U(r) и диэлектрическая проницаемость е(r), фигурирующие в уравнениях Шрлдингера и Максвелла, представляют собой периодические функции с периодами, равными любымвекторам R решетки кристалла и фотонного кристалла, соответственно:
U(r) = U(r + R), (3)
5. Для волновой функции электрона и напряженности электромагнитного поля выполняется теорема Блоха с периодическими функциями u k и u k .
- 6. Возможные значения волновых векторов k заполняют зону Бриллюэна кристаллической решетки или элементарной ячейки фотонного кристалла, задаваемую в пространстве обратных векторов.
- 7. Энергия электрона E, являющаяся собственным значением уравнения Шрлдингера, и собственное значение волнового уравнения (следствия уравнений Максвелла) - частота моды щ - связаны со значениями волновых векторов k блоховских функций (4) законом дисперсии E(k) и щ(k).
- 8. Примесный атом, нарушающий трансляционную симметрию атомного потенциала, является дефектом кристалла и может создавать примесное электронное состояние, локализованное в окрестности дефекта. Изменения диэлектрической проницаемости в определенной области фотонного кристалла нарушают трансляционную симметрию е(r) и приводит к появлению разрешенной моды внутри фотонной запрещенной зоны, локализованной в ее пространственной окрестности.
) — материал, структура которого характеризуется периодическим изменением показателя преломления в 1, 2 или 3 пространственных направлениях.
Описание
Отличительная особенность фотонных кристаллов (ФК) - наличие пространственно периодического изменения показателя преломления. В зависимости от числа пространственных направлений, вдоль которых показатель преломления периодически изменяется, фотонные кристаллы называются одномерными, двумерными и трехмерными, или сокращенно 1D ФК, 2D ФК и 3D ФК (D - от английского dimension) соответственно. Условно структура 2D ФК и 3D ФК показана на рис.
Наиболее яркой чертой фотонных кристаллов является существование в 3D ФК с достаточно большим контрастом показателей преломления компонентов определенных областей спектра, получивших название полных фотонных запрещенных зон (ФЗЗ): существование излучения с энергией фотонов, принадлежащей ФЗЗ в таких кристаллах, невозможно. В частности, излучение, спектр которого принадлежит ФЗЗ, извне в ФК не проникает, существовать в нем не может и полностью отражается от границы. Запрет нарушается только при наличии дефектов структуры или при ограниченных размерах ФК. При этом целенаправленно созданные линейные дефекты являются с малыми изгибными потерями (до микронных радиусов кривизны), точечные дефекты - миниатюрными резонаторами. Практическая реализация потенциальных возможностей 3D ФК, основанных на широких возможностях управления характеристиками световых (фотонных) пучков только начинается. Она затруднена отсутствием эффективных методов создания 3D ФК высокого качества, способов целенаправленного формирования в них локальных неоднородностей, линейных и точечных дефектов, а также методов сопряжения с другими фотонными и электронными устройствами.
Существенно больший прогресс достигнут на пути практического применения 2D ФК, которые используются, как правило, в виде планарных (пленочных) фотонных кристаллов или в виде (ФКВ) (см. подробнее в соответствующих статьях).
ФКВ представляют собой двумерную структуру с дефектом в центральной части, вытянутую в перпендикулярном направлении. Являясь принципиально новым типом оптических волокон, ФКВ предоставляют недоступные другим типам возможности по транспортировке световых волн и управлению световыми сигналами.
Одномерные ФК (1D ФК) представляют собой многослойную структуру из чередующихся слоев с разными показателями преломления. В классической оптике задолго до появления термина «фотонный кристалл» было хорошо известно, что в таких периодических структурах характер распространения световых волн существенно изменяется из-за явлений интерференции и дифракции. Например, многослойные отражающие покрытия давно и широко используются для изготовления зеркал и пленочных интерференционных фильтров, а объемные брэгговские решетки в качестве спектральных селекторов и фильтров. После того, как стал широко употребляться термин ФК, такие слоистые среды, в которых показатель преломления периодически изменяется вдоль одного направления, стали относить к классу одномерных фотонных кристаллов. При перпендикулярном падении света спектральная зависимость коэффициента отражения от многослойных покрытий представляет собой так называемый «брэгговский столик» - на определенных длинах волн коэффициент отражения быстро приближается к единице при увеличении числа слоев. Световые волны, попадающие в спектральный диапазон, показанный на рис. б стрелкой, практически полностью отражаются от периодической структуры. По терминологии ФК эта область длин волн и соответствующая ей область значений энергий фотона (или энергетическая зона) является запрещенной для световых волн, распространяющихся перпендикулярно слоям.
Потенциал практических применений ФК огромен благодаря уникальным возможностям управления фотонами и еще не до конца раскрыт. Нет сомнения, что в ближайшие годы будут предложены новые устройства и конструктивные элементы, возможно принципиально отличающиеся от тех, которые используются или разрабатываются сегодня.
Огромные перспективы применения ФК в фотонике были осознаны после выхода статьи Э. Яблоновича, в которой было предложено использовать ФК с полными ФЗЗ для управления спектром спонтанного излучения.
Среди фотонных устройств, появление которых можно ожидать в ближайшем будущем, следующие:
- низкопороговые ФК лазеры сверхмалых размеров;
- сверхяркие ФК с управляемым спектром излучения;
- сверхминиатюрные ФК волноводы с микронным радиусом изгиба;
- фотонные интегральные схемы с высокой степенью интеграции на основе планарных ФК;
- миниатюрные ФК спектральные фильтры, в том числе перестраиваемые;
- ФК устройства оперативной оптической памяти;
- ФК устройства обработки оптических сигналов;
- средства доставки мощного лазерного излучения на основе ФКВ с полой сердцевиной.
Наиболее заманчивое, но и наиболее трудное в реализации применение трехмерных ФК - создание сверхбольших объемно интегрированных комплексов фотонных и электронных устройств для обработки информации.
Среди других возможных применений трехмерных фотонных кристаллов - изготовление ювелирных украшений на основе искусственных опалов.
Фотонные кристаллы встречаются и в природе, придавая дополнительные оттенки цветовой окраске окружающего нас мира. Так, перламутровое покрытие раковин моллюсков, таких, как галиотисы, имеет структуру 1D ФК, усики морской мыши и щетинки многощетинкового червя представляют собой 2D ФК, а природные полудрагоценные камни опалы и крылья африканских бабочек-парусников (Papilio ulysses) являются природными трехмерными фотонными кристаллами.
Иллюстрации
|
а – структура двумерного (сверху) и трехмерного (снизу) ФК; б – запрещенная зона одномерного ФК, образованного четвертьволновыми слоями GaAs/AlxOy (величина запрещенной зоны показана стрелкой); в – инвертированный ФК никеля, полученный сотрудниками ФНМ МГУ им. М.В. Ломоносова Н.А. Саполотовой, К.С. Напольским и А.А. Елисеевым |
