УДК 535.42
ЗОННАЯ СТРУКТУРА ФОТОННОГО КРИСТАЛЛА С СИММЕТРИЕЙ РЕШЕТКИ КЛАТРАТА SI34
© 2008 П. Н. Дьяченко1, Н. Д. Кундикова2, Ю. В. Микляев2, В.С. Павельев1
1Институт систем обработки изображений РАН 2Вузовско-академический отдел нелинейной оптики ИЭФ УрО РАН и ЮУрГУ, Челябинск
Исследована зонная структура фотонного кристалла с симметрией решетки клатрата Si34, содержащего 34 узла в примитивной ячейке. Обнаружено, что решетка имеет большую изотропную запрещенную зону для широкого диапазона диэлектрической постоянной. Получена зависимость изотропности фотонных зон от диэлектрической постоянной.
Фотонные кристаллы
Фотонными кристаллами называют периодические структуры, состоящие из двух и более материалов с различными диэлектрическими постоянными, обладающие полной фотонной запрещенной зоной (ФЗЗ) [1]. При наличии полной ФЗЗ распространение света подавлено в любом направлении в спектральном диапазоне, совпадающем с запрещенной зоной. Получение трехмерных фотонных кристаллов макроскопических размеров, обладающих ФЗЗ в ближней инфракрасной и видимой областях спектра, позволяет наблюдать целый ряд ранее недоступных для наблюдения эффектов, среди которых наиболее значимым является возможность управления спонтанным излучением [2].
Решетка фотонных кристаллов имеет субмикронные размеры, поэтому для их изготовления требуется разработка новых методов. Кроме того, на данный момент наименьший показатель преломления материала решетки, которая может обладать полной ФЗЗ, равен п = 1,9 для симметрии алмазной решетки [3]. Это накладывает серьезные ограничения на выбор материала, поскольку, например, показатель преломления полимеров и стекол в видимой области обычно ниже этого значения. Одной из возможностей решения этой проблемы является использование трехмерных фотонных квазикристаллов и периодических структур с большим числом узлов («атомов») в примитивной ячейке. Можно ожидать, что большая изотропия позволит снизить порог существова-
ния ФЗЗ по показателю преломления. В микроволновом диапазоне было экспериментально обнаружено [4], что икосаэдрический фотонный квазикристалл, изготовленный методом стереолитографии, имеет большие стоп-зоны. В работе [4] были экспериментально измерены коэффициенты пропускания микроволнового излучения икосаэдри-ческого фотонного квазикристалла, но теоретический анализ не был произведен. В работе [5] был теоретически исследован фотонный кристалл с симметрией решетки пирохлора, содержащий 4 «атома» в примитивной ячейке. Было показано, что такой фотонный кристалл имеет изотропные ФЗЗ, сравнимые по размеру с ФЗЗ алмазной решетки. Решетки с числом «атомов» в примитивной ячейке, большим четырех, до настоящего момента не исследовались.
Рассмотрим в качестве фотонного кристалла с большим числом узлов в примитивной ячейке, обладающем изотропией свойств, кристалл с симметрией решетки клатрата Б134 [6,7]. Данная решетка относится к классу гранецентрированных кубических (ГЦК) и содержит 34 узла в примитивной ячейке (рис. 1). В Б134 каждый «атом» связан с соседними четырьмя «атомами» искаженными тетраэдрическими связями. Как известно, решетка алмаза, имеющая тетраэдрические связи «атомов», позволяет получать самую большую ФЗЗ из всех фотонных кристаллов. Икосаэдриче-ская же структура не имеет тетраэдрических связей [4].
Рис. 1. Четыре примитивных ячейки решетки Б134
Кремниевые клатраты впервые были синтезированы в 1965 году [8]. В последние годы они стали объектом интенсивных экспериментальных и теоретических исследований, прежде всего с активным поиском новых полупроводниковых материалов. В решетке клатрата Б134 расположение атомов наиболее изотропно, то есть форма зоны Бриллюэна наиболее близка к сфере, что благоприятно для открытия полной ФЗЗ. Решетка клатрата Б134 получается при упаковке пентагонального додекаэдра и тетра-кисдекаэдра. Пентагональный додекаэдр -это правильный платонов полиэдр, то есть такое тело, в котором все вершины эквивалентны, а все грани однотипны. Он имеет несколько осей симметрии пятого порядка, поэтому, согласно законам кристаллографии, невозможно заполнение пространства только этими полиэдрами, необходима их комбинация с полиэдрами другого типа, которые содержат еще и шестиугольные грани. Тетра-кисдекаэдр имеет две шестиугольные грани. При этом, как уже отмечалось, каждый атом
связан с соседними четырьмя атомами, искаженными тетраэдрическими связями.
Все вышеизложенное дает основание полагать, что исследование решеток клатра-тов представляет большой интерес для создания на их основе фотонных кристаллов.
Цель настоящей работы - исследование возможности существования полной ФЗЗ в фотонных кристаллах с симметрией Б134.
При проведении исследований использовалась следующая модель. Диэлектрический материал полагался непоглощающим, немагнитным и изотропным, то есть показатель преломления среды п = у[ё, где е - диэлектрическая проницаемость диэлектрика в оптическом диапазоне. Для нахождения собственных мод уравнений Максвелла с периодическими граничными условиями использовали метод разложения поля по плоским волнам [9]. Рассматривались три случая: 1) узлами решетки являются диэлектрические шары, окруженные воздухом; 2) узлами решетки являются сферические полости в диэлектрике; 3) соседние узлы решетки соединены диэлектрическими стержнями.
скорость света в вакууме. По оси х отмечены точки высокой симметрии зоны Бриллюэна ГЦК решетки [10]. Полная ФЗЗ находится между 34-ой и 35-ой зонами и имеет размер Лю / ют =15,6% (при факторе заполнения
f =22%). Под размером полной ФЗЗ (/) - (/' +1) понимается выражение
Л^ = 2 • т1"(Ю'+-) - Иах(Ю-)-100%, где
ют тт(юг+1) + тах(ю)
тт(ю;+1) и тах(юг.) - минимальная и максимальная частота для зон (/ +1) и (/) соответственно. Таким образом, фотонный кристалл с решеткой симметрии клатрата, содержащий более 4 «атомов» в примитивной ячейке, и обладает полной ФЗЗ. Данный кристалл содержит в своей решетке полиэдр, имеющий оси симметрии пятого порядка.
2
1.8
1.6
1.4
О 1.2
с\| 1
го
3 0.8 0.6 0.4 0.2 0
В результате проведенных расчетов оказалось, что фотонный кристалл на основе решетки БІ34, состоящий из диэлектрических шаров в вакууме, не имеет больших запрещенных зон (порядка 5% при диэлектрической постоянной е = 12), а фотонный кристалл, состоящий из воздушных шаров в диэлектрике, не имеет запрещенных зон. Поэтому более подробно исследовался фотонный кристалл, состоящий из диэлектрических стрежней в вакууме, которые соединяют соседние узлы решетки, образовывая тем самым тетраэдрические связи.
На рис. 2 приведена зонная структура фотонного кристалла при е = 12 и оптимальном факторе заполнения f =22% (отношения объема диэлектрика к полному объему ячейки). На графиках частота указана в безразмерных единицах юа /2 га, где ю - циклическая частота, а - размер примитивной ячейки, с -
Рис. 2. Зонная структура Бі34: диэлектрическая постоянная е =12; фактор заполнения / =22%; размер запрещенной зоны Ают / ю =15,6%
Для определения минимального значения диэлектрической постоянной ей, при
котором появляется полная ФЗЗ, т.е. для определения порога возникновения полной ФЗЗ рассчитана зависимость ширины запрещенной зоны от диэлектрической проницаемости. Для расчета зависимости при каждом значении диэлектрической проницаемости е
определялся фактор заполнения диэлектриком f, для которого полная ФЗЗ имеет минимальный размер. Полученные результаты приведены на рис. 3. Как видно из рис. 3, порог возникновения полной ФЗЗ по диэлектрической постоянной равен е& =5. Значение порога оказалось большим, чем у фотонного
кристалла с симметрией алмазной решетки, для которого порог ет =4,0, но меньшим, чем у фотонного кристалла с симметрией инвертированного опала [11]. Это связано с тем, что ФЗЗ структуры с симметрией Б134 более
изотропная, то есть частота граничных зон (зон, ограничивающих ФЗЗ) слабо зависит от направления распространения электромагнитной волны.
О4
0 123456789 10 11 8
Рис. 3. Зависимость размера запрещенной зоны Лю/ Ют от диэлектрической постоянной е
Проведено исследование изотропности запрещенной зоны. Параметр изотропности Б фотонной зоны (/) определялся следующим образом [12]:
„ тах(щ.) - тт(щ.)
Г = 2-------------------.
тах(ю.) + тт(юг.)
Полученная зависимость параметра изотропности Г для зон 34 («нижняя» зона) и 35 («верхняя» зона) от диэлектрической постоянной е при факторе заполнения f =27% представлена на рис. 4. Из рис. 4 видно, что параметр г монотонно уменьшается при увеличении диэлектрической постоянной. Это обусловлено тем, что локализация электромагнитного поля усиливается с увеличением диэлектрической постоянной. Из ранее известных структур фотонных кристаллов наиболее изотропной ФЗЗ обладал фотонный кристалл с симметрией гранецен-трированной кубической решетки, состоящей из воздушных шаров в диэлектрике.
При f =27% и е =12 такая решетка имеет Г =0,08 для «нижней» зоны и Г =0,066 для «верхней» зоны. При тех же условиях фотонный кристалл с решеткой с симметрией Б134 имеет параметр изотропности Г =0,022 для «нижней» зоны и Г =0,056 для «верхней» зоны. Таким образом, ФЗЗ решетки клатрата Б134 является самой изотропной из когда-либо ранее исследованных решеток.
Известно, что групповая скорость распространения электромагнитных волн может обращаться в ноль на границе ФЗЗ. Физическое значение этого явления заключается в том, что на границе ФЗЗ рассеянные волны образуют стоячую волну, и за счет этого возможна генерация когерентного излучения [13]. В обычных фотонных кристаллах это реализуется только в некоторых направлениях, т. к. положение границы ФЗЗ зависит от направления распространения электромагнитного излучения. Для получения когерентного излучения, вне зависимости от направления распространения электромагнит-
ного излучения, необходимо использовать фотонные кристаллы с изотропной ФЗЗ. Изотропность фотонных зон также можно использовать для получения всенаправленной негативной рефракции. Это явление было исследовано для двухмерных фотонных кристаллов [14] и квазикристаллов [15].
Таким образом, впервые обнаружена полная ФЗЗ в фотонном кристалле симметрией решетки клатрата, порог запрещенной
0,08
0,07
0,06
0,05
0,04
0,03 0,02
7 8 9 10 11 є
Рис. 4. Зависимость параметра изотропности Е для зон 34 (сплошная линия) и 35 (пунктирная линия) от диэлектрической постоянной е. Фактор заполнения /=27%.
зоны по диэлектрической постоянной равен е& =5,0. Обнаруженная полная ФЗЗ является самой изотропной из всех известных.
Факторы изотропности оказались равными Г =0,022 для «нижней» зоны и Г =0,056 для «верхней» зоны при факторе заполнения f =27% и е =12.
Библиографический список
1. Joannopoulos J.D., Meade R.D., Winn J.N., Photonic Crystals: Molding the Flow of Light. Princeton University Press, Singapore, 1999.
2. Быков В.П. // ЖЭТФ. - 1972. -Т. 62. - C. 505.
3. Sharp D.N., Turberfield A.J., Denning R.G. // Phys. Rev. B. - 2003. - V. 68. - P. 205102.
4. Man W., Megens M., Steinhardt P.J., Chaikin P.M. // Nature. - 2005. - V. 436. - P. 993.
5. Garcia-Adeva A.J. // Phys. Rev. B . -2006. - V. 73. - P. 073107.
6. Adams G.B., O'Keeffe M., Demkov A.A., Sankey O.F., Huang Y.M. // Phys. Rev. B.
- 1994. - V. 49. - P. 8048.
7. Blase X. // Phys. Rev. B. - 2003. - V. 67. - P. 035211.
8. Kasper J.S., Hagenmuller P. // C. Cros. Science. - 1965. - V. 150. - P. 1713.
9. Johnson S. G., Joannopoulos J.D. // Optics Express. - 2001. - V. 8. - P. 173.
10. Chutinan A., Noda S. // Phys. Rev. B.
- 1998. - V. 57. - P. 2006.
11. Bush K., John S. // Phys. Rev. E. -1998. - V. 58. - P. 3896.
12. Takeda H., Takashima T., Yoshino K. // J. Phys.: Condens. Matter 2004. - V. 16. - P. 6317
13. Meier M., Mekis A., Dodabalapur A., Timko A., Slusher R.E., Joannopoulos J.D. // Appl. Phys. Lett. - 1999. - V. 74. - P. 7.
14. Gajic R., Meisels R., Kuchar F., Hingerl K. // Phys. Rev. B. - 2006. - V. 73. - P. 165310.
15. Feng Z., Zhang X., Wang Y., Li Z.Y., Cheng B., Zhang D.Z. // Phys. Rev. Lett. -
2005. - V. 94. - P. 247402.
References
1. Joannopoulos J.D., Meade R.D., Winn J.N., Photonic Crystals: Molding the Flow of Light. Princeton University Press, Singapore, 1999.
2. Bykov V.P. // “GETF”. - 1972. -Т. б2. - C. 505.
3. Sharp D.N., Turberfield A.J., Denning R.G. // Phys. Rev. B. - 2003. - V. б8. - P. 205102.
4. Man W., Megens M., Steinhardt P.J., Chaikin P.M. // Nature. - 2005. -V. 43б. - P. 993.
5. Garcia-Adeva A.J. // Phys. Rev. B. -
200б. -V. 73. - P. 073107.
6. Adams G.B., O'Keeffe M., Demkov A.A., Sankey O.F., Huang Y.M. // Phys. Rev. B.
- 1994. - V. 49. - P. 8048.
7. Blase X. // Phys. Rev. B. - 2003. -V. б7. - P. 035211.
8. Kasper J.S., Hagenmuller P. // C. Cros. Science. - 19б5. - V. 150. - P. 1713.
9. Johnson S. G., Joannopoulos J.D. // Optics Express. - 2001. - V. 8. - P. 173.
10. Chutinan A., Noda S. // Phys. Rev. B.
- 1998. - V. 57. - P. 200б.
11. Bush K., John S. // Phys. Rev. E. -1998. - V. 58. - P. 389б.
12. Takeda H., Takashima T., Yoshino K. // J. Phys.: Condens. Matter. - 2004. - V. 1б. - P. б317
13. Meier M., Mekis A., Dodabalapur A., Timko A., Slusher R.E., Joannopoulos J.D. // Appl. Phys. Lett. - 1999. - V. 74. - P. 7.
14. Gajic R., Meisels R., Kuchar F., Hingerl K. // Phys. Rev. B. - 200б. - V. 73. - P. 1б5310.
15. Feng Z., Zhang X., Wang Y., Li Z.Y., Cheng B., Zhang D.Z. // Phys. Rev. Lett., 2005.
- V. 94. - P. 247402.
BAND STRUCTURE OF A PHOTONIC CRYSTAL WITH THE CLATHRATE SI-34
CRYSTAL LATTICE
О 2008 P.N. Dyachenko1, N.D.Kundikova2, Yu.V. Miklyayev2, V.S. Pavelyev1
1Image Processing Systems Institute of the RAS, 443001, Samara, Russia College-academic Non-linear Optics department, Experimental Physics Institute, the Ural Branch of the RAS and South-Ural State University, 454080, Chelyabinsk, Russia
The band structure of a photonic crystal with the clathrate Si-34 lattice comprising 34 lattice sites in the unit cell is studied. The lattice is found to have a large isotropic band gap for a wide range of permittivity. The relationship between the photonic band isotropy and the permittivity is deduced.
Photon crystals
Сведения об авторах
Дьяченко Павел Николаевич, Учреждение Российской академии наук Институт систем обработки изображений РАН, Самара, Россия, стажер-исследователь, e-mail: dyachenko@ssau.ru. Область научных интересов: нелинейная оптика
Кундикова Наталия Дмитриевна, Декан физического факультета Южно-Уральского государственного университета. Область научных интересов: нелинейная оптика.
Микляев Юрий Владимирович, Южно-Уральский государственный университет, Челябинск, Россия, доцент, e-mail: miklyaev@mail.ru. Область научных интересов:
Павельев Владимир Сергеевич, Учреждение Российской академии наук Институт систем обработки изображений РАН, Самара, Россия, главный научный сотрудник, e-mail: pavelyev@smr.ru. Область научных интересов: фотонные кристаллы, нелинейная оптика
Dyachenko Pavel Nikolaevich, Establishment of the Russian Academy of Sciences Institute of systems of processing of images of the Russian Academy of Sciences, Samara, Russia, the trainee-researcher, e-mail: dyachenko@ssau.ru. Area of scientific interests: nonlinear optics.
Kundikova Nataliya Dmitrievna, the Ural Branch of the RAS and South-Ural State University, 454080, Chelyabinsk, Russia. Physical faculty, the dean. Area of scientific interests: nonlinear optics.
Miklyayev Yurii Vladimirovich the Ural Branch of the RAS and South-Ural State University, 454080, Chelyabinsk, Russia, The senior lecturer, e-mail: miklyaev@mail.ru. Area of scientific interests: nonlinear optics.
Pavelyev Vladimir Sergeevich, Establishment of the Russian Academy of Sciences Institute of systems of processing of images of the Russian Academy of Sciences, Samara, Russia, the main scientific employee, e-mail: pavelyev@smr.ru. Area of scientific interests: photon crystals, nonlinear optics.