CC BY
47
1
ФОТОНИКА И ОПТОИНФОРМАТИКА PHOTONICS AND OPTOINFORMATICS
УДК 535.3
ЭФФЕКТ ПЁРСЕЛЛА В ПРЕДЕЛЬНО АНИЗОТРОПНЫХ ЭЛЛИПТИЧЕСКИХ
МЕ ТАМАТЕРИАЛАХ А.В. Чебыкин", А.А. Орлов", Ф. Хайслерь, К.В. Барышникова", П.А. Белов"
a Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected] ь Йенский университет имени Фридриха Шиллера, Йена, 07737, Германия
Аннотация. Теоретически продемонстрирован эффект Пёрселла в предельно анизотропных метаматериалах с эллиптической изочастотной поверхностью. В отличие от гиперболических метаматериалов, данный эффект не связан с расходящейся плотностью состояний. Показано, что большой фактор Пёрселла можно наблюдать и без возбуждения мод с большими волновыми векторами в одном из направлений, при этом нормальная слоям компонента волнового вектора будет меньше k0. Это дает возможность получить в данных материалах увеличение не только мощности, излучаемой в среду, но и мощности, передаваемой в свободное пространство через границу среды, расположенную поперечно слоям материала. Методами анализа являлись построение изочастотных контуров, а также расчет аналитической зависимости фактора Пёрселла от частоты для бесконечной структуры слоистого метаматериала. В диапазоне видимого света сильная пространственная дисперсия не позволяет получить усиление спонтанного излучения в метаматериале с двухслойной элементарной ячейкой. Эффект может быть реализован в периодических слоистых металлодиэлектрических наноструктурах с элементарной ячейкой, содержащей два металлических и два диэлектрических слоя. Анализ полученных зависимостей фактора Пёрселла от частоты показывает, что спонтанное излучение усиливается на порядок и более только для случая ориентации возбуждающего диполя вдоль слоев метаматериала, а для случая поперечной ориентации излучение максимально может усиливаться лишь в 2-3 раза. Результаты работы могут быть использованы для создания нового типа метаматериалов с эллиптическими изочастотными контурами, обеспечивающих более эффективное излучение света в дальнее поле.
Ключевые слова: метаматериалы, эффект Пёрселла, спонтанное излучение, поверхностные плазмоны. Благодарности. Работа выполнена при государственной финансовой поддержке ведущих университетов Российской Федерации (субсидия 074-U01), фонда РФФИ (проект 14-02-31720), а также стипендии Президента Российской Федерации (грант СП-2154.2012.1). Авторы приносят отдельную благодарность анонимному рецензенту за ценные замечания и помощь в подготовке статьи.
PURCELL EFFECT IN EXTREMELY ANISOTROPIC ELLIPTIC METAMATERIALS A.V. Chebykin", A.A. Orlov", F. Heislerb, K.V. Baryshnikova", P.A. Belov"
a ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected] b Friedrich-Schiller-Universitat Jena, Jena, 07737, Germany
Abstract. The paper deals with theoretical demonstration of Purcell effect in extremely anisotropic metamaterials with elliptical isofrequency surface. This effect is free from association with divergence in density of states unlike the case of hyperbolic metamaterials. It is shown that a large Purcell factor can be observed without excitation of modes with large wave vectors in one direction, and the component of the wave vector normal to the layers is less than k0. For these materials the possibility is given for increasing of the power radiated in the medium, as well as the power radiated from material into free space across the medium border situated transversely to the layers. We have investigated isofrequency contours and the dependence of Purcell factor from the frequency for infinite layered metamaterial structure. In the visible light range strong spatial dispersion gives no possibility to obtain enhancement of spontaneous emission in metamaterial with unit cell which consists of two layers. This effect can be achieved in periodic metal-dielectric layered nanostructures with a unit cell containing two different metallic layers and two dielectric ones. Analysis of the dependences for Purcell factor from the frequency shows that the spontaneous emission is enhanced by a factor of ten or more only for dipole orientation along metamaterial layers, but in the case of the transverse orientation radiation can be enhanced only 2-3 times at most. The results can be used to create a new type of metamaterials with elliptical isofrequency contours, providing a more efficient light emission in the far field.
Keywords: metamaterials, Purcell effect, spontaneous radiation, surface plasmons.
Acknowledgements: The work is partially financially supported by the Government of the Russian Federation (grant 074-U01),the Russian Foundation for Basic Research (project 14-02-31720),and the Russian Federation President scholarship (grant СП-2154.2012.1). The authors express their special thanks to an anonymous reviewer for valuable remarks and rendering assistance in paper preparation.
Введение
В настоящее время гиперболические метаматериалы рассматриваются в качестве наиболее подходящих структур для увеличения силы взаимодействия света с веществом в широком диапазоне [1-4]. Такие среды характеризуются тем, что имеют эффективные диэлектрические проницаемости разных знаков, т.е. е^ = еуу< 0 и ezz > 0 (где еxx, еуу и ezz - диагональные компоненты тензора диэлектрической проницаемости) [5]. Плотность состояний фотонов в таких материалах расходится благодаря гиперболическому
виду изочастотных поверхностей, что обеспечивает сильный эффект Пёрселла [6]. Обычно эффект Пёр-селла связан и со скоростью излучательного затухания, и с мощностью, излучаемой в окружающее пространство [7]. Тем не менее, переход от эллиптического режима к гиперболическому сопровождается сокращением времени жизни в возбужденном состоянии и уменьшением наблюдаемой интенсивности излучения в свободное пространство [8, 9]. Причина заключается в том, что спонтанное излучение в гиперболических метаматериалах возникает в основном из-за возбуждения мод с такими волновыми векторами, у которых обе компоненты (вдоль и поперек оптической оси) превышают волновое число свободного пространства. Такие пространственные гармоники при преломлении через границу материала в свободное пространство становятся эванесцентными в обоих практических случаях расположения границы -как вдоль оптической оси, так и ортогонально ей. Эванесцентные волны ввиду экспоненциального затухания при удалении от границы не могут быть детектированы в дальнем поле [10].
В настоящей работе мы теоретически рассматриваем возможность реализации сильного эффекта Пёрселла в предельно анизотропных эллиптических метаматериалах, где выполняется соотношение
0 < ехх = £уу << егг . (1)
В таких средах изочастотная поверхность имеет вид сильно сплющенного эллипсоида (рис. 1). Ус -ловие (1) выполняется по одну сторону от топологического перехода между гиперболическим и эллиптическим режимами одного и того же материала (под эллиптическим и гиперболическим режимами здесь понимается форма изочастотных контуров материала). Данный переход происходит на определенной частоте. Далее мы покажем, что анизотропный эллиптический метаматериал может обладать большими значениями фактора Пёрселла, несмотря на конечную величину плотности фотонных состояний. Более того, в этом случае по одному из направлений распространения (по оси 2) волны с большими волновыми векторами не вносят вклад в спонтанное излучение, поэтому излученная волна может покинуть структуру. В результате этого увеличивается как скорость излучательного затухания (что есть и в гиперболических метаматериалах), так и излучаемая через границу по плоскости Х2 в свободное пространство мощность (чего в гиперболических материалах с плоской границей при ее обычном расположении относительно оптической оси не наблюдается ни экспериментально, ни теоретически).
Изочастотная xj
{ — У
Рис. 1. Изочастотная поверхность материала в предельно анизотропном эллиптическом режиме
Физика эффекта Пёрселла в эллиптических метаматериалах
Причину усиления спонтанного излучения в среде, удовлетворяющей условию (1), можно наиболее просто описать при помощи золотого правила Ферми:
1 2п v-ч
- = y •Ek
5(fc>k- Йш0).
(2)
Здесь т - скорость излучательного затухания; ю0 - частота излучения двухуровневой системы; ё -дипольный матричный элемент; Екс - оператор амплитуды электрического поля, соответствующий единичному кванту излучения; Й - постоянная Планка, 5 - дельта-функция. Суммирование осуществляется по волновому вектору к излучаемой волны. Символом с обозначена поляризация (ТЕ либо ТМ). Электрическое поле может быть записано в следующем виде:
Ek
2пйю
'k-TM (cosek фsinek z):
' тоёе,ТМ,к £гг
где Уто<1е - это эффективный объем моды, 0 и ф - сферические координаты волнового вектора к. ¥тоёе определяется из условия квантования для плоской волны УЕк (еЕк) = 2яЙюкТМ . Здесь V - нормированный объем. Эффективный объем моды может быть выражен через эффективный показатель преломления пТМ:
^тоёе,ТМ,к = ^хс / Пк,ТМ(0к ) , (3)
птм, в свою очередь, определяется дисперсией ТМ-мод: юк,тм = ск / птм (0к), где с - скорость света в вакууме.
Выражение для эффективного показателя преломления может быть записано в виде
,(0,)=
( sin2 0, + cos2 0 v"2
(4)
Подставив уравнения (3) и (4) в уравнение (2) и осуществив интегрирование по к и ф, получим вклад ТМ-мод в скорость излучательного затухания для поляризации по оси X:
1 = ¿1 (Т Г ^к 0к ^ 0кПТм(0к ) = ¿1 (^ (5)
тх,тм 2Й I с ) Г е^ 3ЙI с ) л/ёХТ '
Уравнение (5) расходится при ехх = 0. Эта расходимость имеет место на границе гиперболического (£хх < 0) и эллиптического (0 < ехх << е11) режима, поэтому природа возникновения усиления излучения в данном случае совершенно отличается от гиперболического случая. В частности, в предельно анизотропном эллиптическом режиме (ехх << е11) в спонтанном излучении преобладают волны, распространяющиеся в направлении плоскости симметрии, 0к и п/2:
1 = ¿1 {? ¥2 % (6)
Тх^"2й I с ) ^ V + вхх/егг|5/2 , (6)
где у = 0к - п/2. Резкий максимум в уравнении (6) при у ~ ^ехх / егг << 1 обусловлен расхождением эффективного показателя преломления в уравнении (4), что соответствует стремлению эффективного объема моды в уравнении (3) к 0. Малый эффективный объем моды приводит к сильному эффекту Пёрселла аналогично данному эффекту в резонаторе.
Учтем теперь вклад от ТЕ-мод (он определяется эффективным показателем преломления пТЕ = ) и нормируем результат на скорость излучательного затухания в свободном пространстве
1/ т = 4й2ю3 /(3Йс3). Тогда выражения для фактора Пёрселла в случае излучателей, поляризованных по осямXи 2, можно записать в следующей форме: е 3. е
[ = [ (7)
ригс,х л ригс, у I г у V /
! =-ч/ е . (8)
и ригс ,1 V хх ^ '
Первый член в уравнении (7) представляет собой вклад от ТМ-волн в скорость спонтанного излучения, и он может быть существенно больше в случае предельно анизотропного эллиптического режима. Стоит также заметить, что относительно большой рост фактора Пёрселла по сравнению с единицей будет наблюдаться только для излучателей, поляризованных в плоскости Х7.
Эффект Пёрселла в предельно анизотропном металлодиэлектрическом
метаматериале
Далее перейдем к рассмотрению возможного варианта реализации данного эффекта. Наиболее естественным подходом является применение плазмонного многослойного металлодиэлектрического мета-материала [11], состоящего из периодически чередующихся слоев металла и диэлектрика. Обычно такие структуры описываются в эффективной модели следующими диэлектрическими проницаемостями [12]:
, (ef) _ £(eff) _ = Me Me ьdiel"diel
_ £ _ N ' d + d
Me diel
£(ef ) _ /£-1 (z)\~1 _ dMe ! £Me + ddiel ! £diel
zz * ' d + d
V UMe^ diel ,,
где угловые скобки обозначают пространственное усреднение; ddieh dMe - толщины диэлектрического и металлического слоев; £diel, £Me - их диэлектрические проницаемости. Анализ данных уравнений показал, что условие в уравнении (1) удовлетворяется на частоте, немного превышающей частоту перехода от эллиптического к гиперболическому режиму о ,
£ Me (о*) _- ^т £ « (о*),
dMe
когда обращается в 0. Из условия £f )(о ) > 0 следует, что для выполнения соотношения dMe < ddiei
металлические слои должны быть тоньше диэлектрических. В частности, это значит, что абсолютное
* *
значение £Me больше £diel на частоте о , а также то, что частота о ниже частоты поверхностного плазмо-
n
на, определяемой условием ЕМе = - е^л Это приводит к тому, что в данном случае изочастотный контур содержит не только эллиптическую, но также и дополнительную гиперболическую составляющую из-за эффектов сильной пространственной дисперсии [13]. Следовательно, в данном случае в оптическом диапазоне невозможно получить эффект усиления спонтанного излучения при эллиптической изочастоте метаматериала.
Для получения изолированного эллиптического изочастотного контура мы предлагаем использовать плазмонные многослойные структуры со сложной элементарной ячейкой, обладающей бипериодич-ностью [14]. В частности, мы рассмотрели систему с четырьмя слоями в элементарной ячейке - двумя разными металлическими слоями и двумя одинаковыми диэлектриками, как показано на рис. 2.
Рис. 2. Структура рассматриваемого металлодиэлектрического метаматериала с элементарной ячейкой, состоящей из двух чередующихся слоев диэлектрика с диэлектрической проницаемостью 4,6 и двух слоев
из металлов с разными плазменными длинами волн
Диэлектрическая проницаемость диэлектриков £d была выбрана равной 4,6. Диэлектрическая проницаемость металлических слоев соответствовала модели Друде с разными плазменными длинами волн 1р1 = 250 нм и Xp2 = 200 нм. Толщины всех слоев равнялись 30 нм.
В настоящей работе мы аналитически построили спектральную зависимость фактора Пёрселла для излучателя, помещенного в центр диэлектрического слоя в бесконечной периодической структуре, описанной выше. В расчетах использовался стандартный математический аппарат функций Грина для слоистых систем [15]. Результаты представлены на рис. 3. Сплошная синяя кривая соответствует случаю, когда излучатель расположен вдоль слоев, красная пунктирная кривая - случаю перпендикулярной ориентации излучателя относительно слоев. На рис. 3, б-г, показаны изочастотные контуры на трех фиксированных частотах. Частота 570 ТГц соответствует главному результату: как видно из рисунка, на этой частоте изочастотный контур имеет сильно анизотропный эллиптический вид, а фактор Пёрселла достигает значения 24 (синяя линия). Для частот ниже 570 ТГц большая часть изочастотного контура находится в диапазоне волн, которые при преломлении через вертикальную границу среды (ориентированную поперечно слоям) оказываются распространяющимися в свободном пространстве: kz < k0 = œ/c. Энергия этих волн может эффективно излучаться в дальнее поле через границу материала, а пространственный спектр таких волн в предельно анизотропном эллиптическом режиме довольно широк, потому что горизонтальная полуось эллипса kymax >> k0, причем большая нормальная к границе компонента волнового вектора не препятствует преломлению (эффект полного внутреннего отражения полностью связан с kz). Однако здесь требуется уточнение. Энергия пространственных гармоник, создаваемых оптически коротким диполем, распределена по изочастотному контуру неравномерно. Так, у дипольного излучателя, ориентированного нормально к слоям и создающего поэтому в основном вертикальную поляризацию электрического поля, превалируют пространственные гармоники с малыми значениями ky, потому что именно у гармоник ky << kz поле поляризовано практически вдоль z. Для такого излучателя вытянутый вдоль оси у изочастотный контур особого значения не имеет. Большие значения фактора Пёрселла по сравнению с единицей в эллиптическом режиме наблюдаются только для излучателя, ориентированного параллельно слоям. В полном соответствии с физикой в эллиптическом режиме, который наблюдается в диапазоне 550-570 ТГц, фактор Пёрселла для поперечно ориентированного диполя сравнительно мал (красный пунктир), а для тангенциального диполя весьма велик (15-20). Такая зависимость фактора Пёрселла от поляризации, полученная в результате точных симуляций, согласуется с предсказаниями эффективной среды в уравнениях (7), (8). На рис. 3, а, зеленым штрих-пунктиром показан результат расчета дипольного излучения в модели эффективной среды, полученный с использованием значений £xx и £zz, извлеченных из изочастотных контуров. Функция в уравнении (7) умножена на 1,4, что может быть объяснено поправ-
кой локального поля [16]. Данное полуаналитическое выражение хорошо описывает полученную численно частотную зависимость для фактора Пёрселла (см. синюю и зеленую кривую на рис. 3, а).
а б
Рис. 3. Зависимость фактора Пёрселла от частоты. Синяя сплошная кривая соответствует случаю ориентации диполя вдоль слоев, красный пунктир - случаю поперечной ориентации. Зеленый штрих-пунктир - результат эффективной модели (7), скорректированный коэффициентом 1,4 (а). Изочастотные контуры структуры на частотах 570, 595, 605 ТГц соответственно (б)-(г).
В диапазоне 570-595 ТГц изочастотный контур превращается из эллипса в «гантелю», а на частотах свыше 595 ТГц распадается на две гиперболические ветви. В данном режиме скорость спонтанного излучения квантового излучателя велика, но в нем доминируют волны с большими значениями вертикальной компоненты волнового вектора (kz > k0), которые не могут покинуть структуру и распространяться в дальнем поле.
Заключение
В работе продемонстрирован сильный (порядка 20) эффект Пёрселла в физически реализуемых предельно анизотропных эллиптических метаматериалах. Такие метаматериалы могут быть реализованы в виде слоистых металло-диэлектрических с элементарной ячейкой, состоящей из двух разных металлических слоев (например, золото и серебро) и двух одинаковых диэлектриков. Результаты работы могут быть полезны для получения абсолютно нового типа метаматериалов, которые заслуживают дальнейших исследований, а также для усиления спонтанной эмиссии квантовых излучателей за счет эффекта Пёр-селла.
Литература
1. Jacob Z., Shalaev V.M. Plasmonics goes quantum // Science. 2011. V. 334. N 6055. P. 463-464.
2. Cortes C.L., Newman W., Molesky S., Jacob Z. Quantum nanophotonics using hyperbolic metamaterials // Journal of Optics. 2012. V. 14. N 6. Art. 063001.
3. Drachev V.P., Podolskiy V.A., Kildishev A.V. Hyperbolic metamaterials: new physics behind a classical problem // Optics Express. 2013. V. 21. N 12. P. 15048-15064.
4. Poddubny A., Iorsh I., Belov P., Kivshar Y. Hyperbolic metamaterials // Nature Photonics. 2013. V. 7. N 12. P. 958-967.
5. Felsen L., Marcuvitz N. Radiation and Scattering of Waves. NY: Wiley, 2003. 464 p.
6. Jacob Z., Smolyaninov I.I., Narimanov E.E. Broadband Purcell effect: radiative decay engineering with metamaterials // Applied Physics Letters. 2009. V. 100. N 18. Art. 181105.
7. Kavokin A., Baumberg J.J., Malpuech G., Laussy F.P. Microcavities. Oxford University Press, 2007. 430 p.
8. Tumkur T., Zhu G., Black P., Barnakov Y.A., Bonner C.E., Noginov M.A. Control of spontaneous emission in a volume of functionalized hyperbolic metamaterial // Applied Physics Letters. 2011. V. 99. N 15. Art. 151115.
9. Kim J., Drachev V.P., Jacob Z., Naik G.V., Boltasseva A., Narimanov E.E., Shalaev V.M. Improving the radiative decay rate for dye molecules with hyperbolic metamaterials // Optics Express. 2012. V. 20. N 7. P. 8100-8116.
10. Lu D., Kan J.J., Fullerton E.E., Liu Z. Enhancing spontaneous emission rates of molecules using nanopatterned multilayer hyperbolic metamaterials // Nature Nanotechnology. 2014. V. 9. N 1. P. 48-53.
11. Orlov A.A., Zhukovsky S.V., Iorsh I.V., Belov P.A. Controlling light with plasmonic multilayers // Photonics and Nanostructures - Fundamentals and Applications. 2014. V. 12. N 3. P. 213-230.
12. Agranovich V.M., Kravtsov V.E. Notes on crystal optics of superlattices // Solid State Communications. 1985. V. 55. N 1. P. 85-90.
13. Orlov A.A., Voroshilov P.M., Belov P.A., Kivshar Y.S. Engineered optical nonlocality in nanostructured metamaterials // Physical Review B - Condensed Matter and Materials Physics. 2011. V. 84. N 4. Art. 045424.
14. Orlov A.A., Krylova A.K., Zhukovsky S.V., Babicheva V.E., Belov P.A. Multiperiodicity in plasmonic multilayers: general description and diversity of topologies // Physical Review A - Atomic, Molecular, and Optical Physics. 2014. V. 90. N 1. Art. 013812.
15. Tomas M.S., Lenac Z. Spontaneous-emission in an absorbing Fabry-Perot cavity // Physical Review A -Atomic, Molecular, and Optical Physics. 1999. V. 60. N 3. P. 2431-2437.
16. Poddubny A.N., Belov P. A., Ginzburg P., Zayats A.V., Kivshar Y.S. Microscopic model of Purcell enhancement in hyperbolic metamaterials // Physical Review B - Condensed Matter and Materials Physics. 2012. V. 86. N 3. Art. 035148.
Чебыкин Александр Васильевич
Орлов Алексей Анатольевич Хайслер Фабиан
Барышникова Ксения Владимировна Белов Павел Александрович
Alexander V Chebykin Alexei A. Orlov Fabian Heisler
Ksenia V. Baryshnikova Pavel A. Belov
Принято к печати 09.09.14 Accepted 09.09.14
— аспирант, инженер-исследователь, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected]
— аспирант, младший научный сотрудник, Санкт-Петербург, 197101, Российская Федерация, [email protected]
— студент, Йенский Университет имени Фридриха Шиллера, Йена, 07737, Германия; стажер-исследователь, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected]
— аспирант, инженер-исследователь, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected]
— доктор физико-математических наук, главный научный сотрудник, заведующий лабораторией, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected]
— postgraduate, research engineer, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
— postgraduate, junior scientific researcher, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
— student, Friedrich-Schiller-Universität Jena, Jena, 07737, Germany; trainee researcher, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
— postgraduate, research engineer, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
— D.Sc., chief research fellow, head of laboratory, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
