CC BY
46
РАССЕЯНИЕ СВЕТА НА ДИЭЛЕКТРИЧЕСКОМ ЦИЛИНДРЕ, ВКЛЮЧАЮЩЕМ ДВУМЕРНУЮ РЕШЕТКУ МЕТАЛЛИЧЕСКИХ НАНОСТЕРЖНЕЙ
Д.В. Нестеренко, В.В. Котляр Институт систем обработки изображений РАН Самарский государственный аэрокосмический университет им. С.П. Королева
Аннотация
Рассмотрена задача распространения ТМ и ТЕ поляризованного света в двумерных массивах серебряных наностержней различного радиуса и концентрации в диэлектрической структуре. Было проведено сравнение точных численных результатов с расчетами на основе нелокальной модели однородной среды для периодичных массивов стержней с отрицательной диэлектрической постоянной. Показано сильное влияние мультипольных вкладов на рассеяние, что не может быть учтено рассматриваемой моделью.
Введение
Металлические частицы с размерами от 5 нм до нескольких сотен нанометров, распределенные периодично или случайно на диэлектрической подложке, или слои таких частиц имеют интересные оптические свойства и изучались долгое время [1]. В последнее время расширились возможности по изготовлению таких систем с периодичным распределением одинаковых сфер. В зависимости от размера и распределения сфер они могут быть очень хорошими поглотителями света [2, 3] или фотонными кристаллами действующие как непоглощающие зеркала в определенном частотном диапазоне электромагнитного спектра [4].
Такие структуры являются композитными средами, включающими два и более различных материалов с характерными размерами включений, много меньшими длины волны падающего света. Для анализа подобных структур применяют теорию эффективных сред (ТЭС) в виде соотношений Мак-свелла-Гарнета [5], когда один из компонентов может рассматриваться как основной, в который встроены включения другого материала. Они предполагают точное вычисление поля в однородном основном материале, произведенного одиночным сферическим или эллиптическим включением, и модификация его с учетом взаимодействия между включениями. Теория Максвелла-Гарнета пригодна для низких концентраций (менее 30 %) относительно маленьких наночастиц (до нескольких десятков нанометров).
В работе [6] предложен метод расчета эффективной диэлектрической проницаемости композитной среды на основе детального учета запаздывания радиационного взаимодействия между составляющими среду сферическими частицами, расстояние между которыми сопоставимо с длиной волны излучения. Представленный метод позволяет рассчитывать оптические свойства композитных сред, в которых расстояние между частицами сопоставимо с длиной волны внешнего излучения. Исследования сферических серебряных включений в диэлектрическую среду показали влияние радиуса сфер и расстояний между ними на отклонение результатов теории
Максвела-Гарнета от решения, учитывающего зависимость оптических свойств металлических нано-шаров от их размеров. В частности, для среды с на-ношарами радиусом 10 нм и концентрацией 10 % отличие в результатах численного расчета показателя преломления композитного материала несущественно. Однако уже при радиусе 20 нм разница достигает нескольких десятков процентов.
В работе [7] было исследовано прохождение ТМ и ТЕ поляризованных волн в двумерных массивах серебряных наностержней различного диаметра в желатиновой среде. Авторы провели сравнение метода рекурсивных функций Грина с теорией Мак-свелла-Гарнета на примере изучения характеристик пропускания периодических и непериодических массивов наностержней. Было показано сильное влияние взаимодействия между наностержнями, вкладов мультиполей и формирования фотонных запрещенных зон на оптические спектры, что не может быть учтено в рамках классической ТЭС. Авторы отметили частичную адекватность теории Максвела-Гарнета для ТМ поляризованного света и ее непригодность для случая ТЕ поляризации.
В [8, 9] была предложена нелокальная модель однородной среды для периодичных массивов стержней с отрицательной диэлектрической постоянной на основе определения поляризуемости одиночного стержня и применения методов интегральных уравнений. Исследование модели было проведено для низких концентраций стрежней (менее 0,8 %) и малого радиуса стержней (до 12 нм для видимого диапазона). Сравнение результатов моделирования задачи отражения на конечной пластине со стержнями, ориентированными параллельно поверхности пластины, полученных с использованием нелокальной модели и электромагнитного подхода, показало адекватность подхода для ТЕ и ТМ поляризаций.
В настоящей статье проводится моделирование микроструктур, включающих массивы наностерж-ней, объединенным методом «конечных элементов Галеркина - граничных элементов» (МКЭГ-ГЭ) [10], основанным на точной электромагнитной теории. Полученные результаты сравниваются с моделированием данных структур с усредненной диэлек-
трическои проницаемостью, рассчитанной в нелокальной модели однородной среды для периодичных массивов стержней с отрицательной диэлектрической постоянной. Прототипы подобных массивов были недавно изготовлены экспериментально [11, 12]. Эти массивы представляют собой упорядоченные длинные стержни из диэлектрического или металлического материала диаметром нескольких десятков нанометров.
Модель эффективной диэлектрической проницаемости
В работе [9] рассматривается матрица из диэлектрического материала с двумерной решеткой нанос-тержней с отрицательной диэлектрической постоянной. Радиус стержней Я, они расположены в квадратной сетке, расстояние между стержнями (постоянная сетки) а. Диэлектрическая постоянная стержней ет. Диэлектрическая постоянная материала е^. В данной работе рассматривается распространение монохроматических полей с временной зависимостью вида в+ш'.Структура бесконечна по оси г, поэтому решение уравнений Максвелла разделяется на ТЕ (вектор электрического поля параллелен г) и ТМ (вектор магнитного поля параллелен г) случаи.
В [9] было показано, что для длин волн много больших радиуса стержней поле в композитной периодической структуре может быть описано с использованием нелокальной диэлектрической проницаемости:
2е,
= е +----, (1)
1 ет +е-
/ ет-е-
еТЕ К ) = е- +-
-1
Р2 - К
(2)
(ет-е- )/ в
где еТМ ,еТЕ - диэлектрические проницаемости для
ТМ и ТЕ поляризаций, / = пЯ2/а2 - концентрация стержней, Кг - г-компонента волнового вектора,
в = соу!е-^0 - волновое число в диэлектрической
структуре, вр - плазменное волновое число для идеально проводящих стержней приводится в [8]:
(вра) -. (3)
1п
а - | + 0,5275
2пЯ )
Заметим, что в [8, 9] ур. (1) для еш приводится только для стационарного предела (к = 0 и в = 0), тогда как ур. (2) для еТЕ выводится для динамического случая. Воспользовавшись вычислениями авторов в [8], запишем выражение для еш для динамического случая:
етм (а) = е- +-
1 ет + е- в2
(4)
Из ур. (2) и (4) следует, что среда обладает пространственной дисперсией. Так, диэлектрическая проницаемость зависит не только от частоты падающего света в, но и от компонента волнового вектора, параллельного стержням. Предложенная модель также действует для диэлектриков с положительной вещественной частью диэлектрической проницаемости.
Массивы наностержней в диэлектрическом цилиндре Для проверки предложенной модели эффективной диэлектрической проницаемости мы рассмотрим рассеяние плоской волны с длиной волны 0,5 мкм на диэлектрическом цилиндре с серебряными стержнями. Как показано на рис. 1, стержни находятся в цилиндре круглого сечения с диэлектрической проницаемостью е- = 2,25. Радиус цилиндра во всех экспериментах был принят равным 0,5 мкм. Материал стержней - серебро (ет = - 9,49 + 1,4831). Для простоты мы ограничимся случаем падения волны с Кг = 0.
&
2/ ет -е- Р1
Рис. 1. Массивы наностержней в диэлектрическом цилиндре, бесконечно простирающемся по оси г
Далее, для демонстрации точности теоретических результатов, проведена серия численных экспериментов моделирования дифракции света на структуре с представленной геометрией с вариацией значений концентрации стержней / от 0 до 1 и вариацией значений радиуса стержней от 5 нм до 30 нм. Содержание стержней 100 % соответствует полностью металлическому цилиндру, содержание 0 % -полностью диэлектрическому. Выбор значений радиуса стержней был продиктован значительным отличием поляризационных характеристик маленьких и больших наночастиц [13]. Если наночастица достаточно мала (Я << Х0, где Х0 - длина волны падающего света), то, в соответствии с теорией Ми, вклад в спектр поглощения дают только дипольные плазмон-ные колебания. Тогда как, для больших частиц вклад в спектр дают также и резонансы высокого порядка. Когда наноразмерная частица освещается светом, электрическая компонента электромагнитного поля возбуждает коллективные колебания электронной плазмы внутри частиц - плазмоны. Если частицы расположены в цепочки эти плазмонные колебания принимают резонансный характер, что облегчает прохождение волн резонансных частот. Такие цепочки исследовались в литературе [9, 14] как потенци-
альные кандидаты для субволнового волноводного распространения света.
Эффективная диэлектрическая проницаемость рассчитывалась по формулам (2) и (4). Объединенным методом МКЭГ-ГЭ проводилось моделирование дифракции плоской волны на цилиндре с вычисленной эффективной диэлектрической проницаемостью. С другой стороны, этим же методом также моделировалась дифракция на структуре с детальной геометрией стержней. Для оценки дифракционного рассеяния использована диаграмма направленности рассеяния, зависящая от угловой координаты р и определенная в бесконечно удаленных точках (l ^ да) как а(р) = lim 2nl \ usc \2/ \u'n \2, где usc и um - рассеянное и падающее поля. В частности, рассеяние определялось в прямом (р= 0°), в обратном (р= 180°) и в поперечном направлениях (р = 90°).
ТЕ поляризация
Рассмотрим прохождение ТЕ поляризованного света через структуру с геометрией, представленной на рис. 1. На рис. 2 представлены зависимости от концентрации f отношения значений рассеяния TE волн а к значениям рассеяния а0 диэлектрического цилиндра для расчетных углов р для стержней радиусом 5, 10, 20, 30 нм.
ст/ст „ 4 3 2 1
а/а,
ТЕ, ср=0°
— 5 нм 10 нм
м ■ 20 нм ■ 30 нм
L b......... ™ ■--
\ р
ТЕ, ф=180°
Л ------ 5 нм - 10 нм ............ 20 нм 30 нм
! \ \ — ■—
£ Yt
% Ж
0
20
40
60
80
Рис. 2. Зависимость рассеяния ТЕ волны на структуре точной геометрии от концентрации стержней/
Рассматривая зависимости рассеяния в прямом направлении можно отметить незначительное отли-
чие распределении для стержней различных радиусов. В диапазоне концентрации от 6 % до 20 % с увеличением радиуса стержней наблюдается увеличение рассеяния вперед в несколько раз. Рассеяние в поперечном направлении для стержней радиусом 5, 10 и 20 нм и концентрацией от 6 % до 12 % существенно снижается, для стержней радиусом 30 нм в данном диапазоне концентрации наблюдается обратная зависимость. Рассеяние в обратном направлении для стержней радиусом 5 и 10 нм имеет минимумы на концентрациях 3 и 11 %. Минимумы смещаются на концентрации 3 и 7 % соответственно для радиуса стержней 20 нм и 3 и 6 % для радиуса 30 нм.
На рис. 3 представлены зависимости от концентрации / отношения отклонений Аа значений рассеяния ТЕ волн а от значений а,$, рассчитанных на основе эффективной диэлектрической проницаемости, к значениям рассеяния а0 диэлектрического цилиндра для расчетных углов ф. По представленным зависимостям можно отметить рост значений отклонений с ростом радиуса стержней на всем диапазоне концентраций. Неравномерность роста можно объяснить резонансными эффектами, связанными с взаимодействием плазмонов в стержнях и зависящих от расстояния между ними.
Аа/а0 1,5
1,0 0,5
0
Аа/а„ 3,0
2,5
2,0
1,5
1,0
0,5
0
Аа/а0 3,0 2,5 2,0 1,5 1,0 0,5
0
ТЕ, ф=0°
1
5 нм 10 нм ■ 20 нм 30 нм
¡\ \
—.— -----
20
40
60
80
ТЕ, ф =90°
------ 5 НМ - 10 НМ ............ 20 нм ■ 30 нм
Л '
1 1
____
| V/ Ufr.—¿-п. ....."-SSs:
20
40
60
80
А ТЕ, ф=180°
------ 5 нм - 10 нм
—4 I V ............ 20 нм 30 нм
1 \
---------------- ■ ■—— ^.....
20
40
60
80
Рис. 3. Зависимость от концентрации стержней/ отклонений рассеяния ТЕ волны на структуре точной геометрии от отклонений рассеяния на структуре с эффективной диэлектрической проницаемостью
На рис. 4 представлены зависимости от концентрации ф относительных отклонений Ао значений рассеяния ТЕ волн а от значений Оф рассчитанных на основе эффективной диэлектрической проницаемости.
Аст/а0
0,7 0,6 0,5 0,4 0,3 0,2 0,1
О
А а/а0 1,6 1,2 0,8 0,4
О
А а/а0 1,6
1,2
0,8
0,4
О
ТЕ, ф=0°
------ 5 нм
Н ! 10 нм
1 ............ 2 и нм ---их А
1 *
—___
|\ \ __
20
40
60
80
1-7 "Е, ф=90°
1
1 ------ 5 нм - 10 нм ............ 20 нм ---30 нм
■—-
1Ш1 "—___
^—
20
40
60
80
11В Е, ф=180°
нв и.
! 1 М \ \ \ 1
11М
II/ \
1 = у
¡у , .__
|| 1 \ _____
ш ¡5
20
40
60
80
Рис. 4. Зависимость от концентрации стержнейф относительных отклонений рассеяния ТЕ волны на структуре точной геометрии от отклонений рассеяния на структуре с эффективной диэлектрической проницаемостью
Среднее по концентрации относительное отклонение представлено в таблице:
Радиус стержней Направление рассеяния
(р= 0° р= 90° р= 180°
5 нм 1 % 8 % 5 %
10 нм 2 % 16 % 13 %
20 нм 11 % 19 % 41 %
30 нм 15 % 60 % 579 %
Таким образом, из сводной таблицы следует, что для случая прохождения ТЕ волны через массив металлических стержней в диэлектрической структуре рассматриваемая нелокальная модель однородной среды адекватно описывает среду со стержнями радиусом Я < Х0/50.
ТМ поляризация Рассмотрим прохождение ТМ поляризованного света через рассматриваемую структуру. На рис. 5
представлены зависимости от концентрации ф отношения значений рассеяния ТМ волн а к значениям рассеяния а0 диэлектрического цилиндра для расчетных углов (р .
Аст/ст0 4
3
2
1
_ ТМ, Ф=0°
/ \
/ \
1 '
—
II 4 ■ 5 нм
/// г-чГч)«") III
Щ/
20
40
60
80
Аа/а0 5
4
3
2
1
О
А а/а0 200
150
100
50
О
ТМ, ф =90°
1 | ------ 5 нм
1 | Л) нм ............ 20 нм ■ 30 нм
» 1 И |
"" "ч.
тргч
20
40
60
80
ТМ, <р=180°
■■■ 5 нм ............ 20 нм - 10 нм ---30 нм
У—
1 I ¿х/ /
V
20
40
60
80
Рис. 5. Зависимость рассеяния ТМ волны на структуре точной геометрии от концентрации стержней ф
Зависимости рассеяния в прямом направлении слабо отличаются для радиусов 5, 20 и 30 нм. Рассеяние в поперечном направлении для радиуса стержней 5 нм не имеет резонансных экстремумов. Для Я > 20 нм появляются ярко выраженные минимумы на порядок меньшие рассеяния диэлектрической структуры. Рассеяние назад для металлической структуры на два порядка больше рассеяния назад диэлектрической структуры. При увеличении радиуса стержней до 30 нм в зависимости рассеяния назад от концентрации появляются минимумы одного порядка с рассеянием на диэлектрической структуре.
На рис. 6 представлены зависимости от концентрации f отношения отклонений Ао значений рассеяния ТМ волн о от значений ое(Ь рассчитанных на основе эффективной диэлектрической проницаемости, к значениям рассеяния а0 диэлектрического цилиндра для расчетных углов р. По представленным зависимостям не отмечается явная зависимость значений отклонений от радиуса стержней.
Аа/а0 2,0
1,5
1,0
0,5
0
Аа/а0 4
3
2
1
Л ТМ, ср=0°
А ■■■ 5 нм - 10 нм ............ 20 нм ---30 нм
¡Л
и/ \ /
20
40
60
80
1 ТМ, ф=90°
1 1 | ------ 5 нм - 10 нм ............ 20 нм ---30 нм
\ , 1
(и
20
40
60
80
[\
\/ \
V,
ш^
ТМ, ч=180°
5 нм 10 нм 20 нм 30 нм
20
40
60
80
Рис. 6. Зависимость от концентрации стержней/ отклонений рассеяния ТМ волны на структуре точной геометрии от отклонений рассеяния на структуре с эффективной диэлектрической проницаемостью
На рис. 7 представлены зависимости от концентрации / относительных отклонений Аа значений рассеяния ТМ волн а от значений ст^, рассчитанных на основе эффективной диэлектрической проницаемости.
Аст/ст0 1,5
1,0
0,5
0
Дс/а„ 1,5
1,0
0,5
0
ТМ. т=0°
1 ■■■ 5 нм - 10 нм ............ 20 нм ---30 нм
1 Лгч
Ил
1 / ч л 1 / V ---
г .... ./
20
40
60
80
ТМ. т=90°
\ 1
—И ■■■ 5 нм - 10 нм ............ 20 нм
\ \
\ -- 30 нм
\ А- —
./ "ч„ 'V.
\ / \Ч "Чц.
Аа/а0 1,5
1,0
0,5
1 ТМ (п=1Я0°
\ 1
1 ■■■ 5 нм - 10 нм ............ 20 нм ---30 нм
д \ 1
1 /•
\ \ / V / 'Ч.
V
И\ .- у ----
0
20
40
60
80
Рис. 7. Зависимость от концентрации стержней/ относительных отклонений рассеяния ТМ волны на структуре точной геометрии от отклонений рассеяния на структуре с эффективной диэлектрической проницаемостью
Среднее по концентрации относительное отклонение представлено в таблице:
Радиус стержней Направление рассеяния
< = 0° <= 90° <= 180°
5 нм 29 % 27 % 50 %
10 нм 28 % 162 % 23 %
20 нм 33 % 483 % 21 %
30 нм 34 % 169 % 106 %
20
40
60
80
Среднее относительное отклонение рассеяния для случая прохождения ТМ волны заметно больше среднего относительного отклонения в случае прохождения ТЕ волны. Можно отметить, что для концентраций / < 1 % среднее отклонение не превышает 10 %, если радиус стержней меньше 20 нм.
Таким образом, рассматриваемая нелокальная модель однородной среды не может адекватно описать резонансные процессы высокого порядка, происходящие в стержнях большого радиуса. Также стоит отметить, что учет взаимодействия между на-ностержнями, осуществленный в модели, не полностью описывает реальные процессы. В случае малых концентраций стержней и, соответственно, большого расстояния между ними расчеты показывают низкое отклонение. С уменьшением расстояния между стержнями, погрешность увеличивается.
Заключение
В работе рассмотрена задача распространения ТМ и ТЕ поляризованного света в двумерных массивах серебряных наностержней различного радиуса и концентрации в диэлектрической структуре. Сравнение точных численных результатов с расчетами на основе нелокальной модели однородной среды для периодичных массивов стержней с отрицательной диэлектрической постоянной было проведено объединенным методом МКЭГ-ГЭ. На основе сравнения были установлены границы применимости нелокальной модели.
Работа выполнена при поддержке российско-американской программы «Фундаментальные исследования и высшее образование» (грант СКОБ ЯиХ0-014-8Л-06), и грантов РФФИ №№05-08-50298, 07-07-97600 и 08-07 99007.
Литература
1. Abeles F., Optical properties of discontinuous thin films and rough surfaces of silver // Abeles F., Borensztein Y., Lopez-Rios T., -Advances in solid state physics, Braunschweig: Vieweg, 1984. Vol. 24. P. 93-117.
2. Taleb A., Collective optical properties of silver nanoparti-cles organized in two-dimensional superlattices // Taleb A., Russier V., Courty A., Pileni M.P.,- Phys. Rev. B,
1999. Vol. 59(20). P. 13350-13358.
3. Yannopapas V., Scattering and absorption of light by periodic and nearly periodic metallodielectric structures // Yannopapas V., Modinos A., Stefanou N.,- Opt. Q. Electr., 2002. Vol. 34. N.1-3. P. 227-234.
4. Zhang W.Y., Robust photonic band gap from tunable scatterers // Zhang W.Y., Lei X.Y., Wang Z.L., Zheng D.G., Tam W.Y., Chan C.T., Sheng P.,- Phys. Rev. Lett.,
2000. Vol. 84(13). P. 2853-2856.
5. Maxwell-Garnett J. C., Colours in metal glasses and in metallic films // Maxwell-Garnett J. C.,- Philos. Trans. R. Soc. London Ser. A, 1904. Vol. 203. P. 385-420.
6. Сухов С.В. Нанокомпозитный материал с единичным показателем преломления // Квантовая электроника, 2005. В. 35. № 8.
7. Rahachou A.I., Light propagation in nanorod arrays // Rahachou A.I., Zozoulenko I.V.,- J. Opt. A: Pure Appl. Opt., 2007. Vol. 9. P. 265-270.
8. Silveirinha M.G. Nonlocal homogenization model for a periodic array of g-negative rods // Phys. Rev. E, 2006. Vol. 73. P. 046612.
9. Silveirinha M.G., Subwavelength imaging at infrared frequencies using an array of metallic nanorods // Silveirinha M.G., Belov P.A., Simovski C.R.,- Phys. Rev. B, 2007. Vol. 75. P. 035108.
10. Нестеренко Д.В., Анализ дифракции света на элементах цилиндрической микрооптики объединенным методом конечных элементов Галеркина и граничных элементов // Нестеренко Д.В., Котляр В.В.,- Компьютерная оптика, 2007. В. 31. №2. С. 9-15.
11. Ahsanulhaq Q, Controlled selective growth of ZnO nano-rod arrays and their field emission properties // Ahsanul-haq Q, Kim J.H., Hahn Y.B.,- Nanotechnology, 2007. Vol. 18. N. 40. P. 485307.
12. Yang J., Controlled growth of aluminium nitride nanorod arrays via chemical vapour deposition // Yang J., Liu T.W., Hsu C.W., Chen L.C., Chen K.H., Chen C.C.,-Nanotechnology, 2006. Vol. 17. N. 11. P. 321-326.
13. Kelly K., The optical properties of metal nanoparticles: the influence of size, shape, and dielectric environment // Kelly K., Coronado E., Zhao L., Schatz G.,- J. Phys. Chem. B, 2003. Vol. 107. P. 668-677.
14. Quidant R., Frustrated energy transport through micro-waveguides decorated by gold nanoparticle chains // Quidant R., Leveque G., Weeber J.-C., Dereux A., Girard C., Weiner J.,- Europhys. Lett., 2004. Vol. 66 (6). P. 785-791.
