CC BY
55
УДК 535.016
НАНОМАТЕРИАЛ ПОВЫШЕННОЙ ПРОЗРАЧНОСТИ
А.С. Шалин1,2, А.С. Кадочкин2
1 Ульяновский филиал Института радиотехники и электроники им. В.А. Котельникова РАН
2 Ульяновский государственный университет
Поступила в редакцию 20.11.2012
Исследованы оптические свойства пространственно упорядоченных нанокомпозитов, состоящих из сферических диэлектрических наночастиц. Представлено аналитическое решение, позволяющее использовать формулы Эйри (метод «мнимой границы») для описания взаимодействия света с упорядоченным нанокомпозитом без усреднения оптических постоянных по объёму. Показано хорошее согласие полученных результатов с точным численным расчетом. Ключевые слова: метаматериал, наночастица, формулы Эйри, мнимая граница.
ВВЕДЕНИЕ
Большое внимание привлекает в последнее время исследование оптических свойств искусственных материалов, состоящих из различного рода наноразмерных объектов, внедренных в некую матрицу, и получивших в литературе название «метаматериалы» [1-4]. Указанный интерес обусловлен возможностью достижения ряда необычных эффектов, обусловленных композиционной структурой среды, а также спецификой рассеяния света наночастицами. Действительно, за счет варьирования материальных и геометрических параметров системы достигается, например, возможность получения сред с гигантским, сверхмалым либо отрицательным показателем преломления [1,5,6], а также показателем преломления, реальная часть которого близка к единице [7].
Для теоретического описания указанных структур в работах [8,9] применяется обобщенная теория Ми. Хотя решение, предложенное Ми [10], получено для дифракции на одной сфере, оно может быть распространено и на системы из любого количества взаимодействующих частиц, путем введения формализма Bj-символов и символов Клебша-Гордана, позволяющих учесть многократность рассеяния света на нанокластерах [8]. Однако, использование данного подхода не всегда целесообразно. Действительно, взаимодействие частиц является весьма дальнодействующим и требует учета весьма большого количества взаимно влияющих элементов. Вследствие этого, точное решение задачи указанным методом возможно лишь численно и требует больших затрат времени. Подобные же затруднения имеет место в случае применения иных методов, использующих Шалин Александр Сергеевич, кандидат физико-математических наук, старший научный сотрудник Ульяновского отделения ИРЭ РАН, старший научный сотрудник НИТИ. E-mail: shalin_a@rambler.ru
Кадочкин Алексей Сергеевич, кандидат физико-математических наук, доцент кафедры радиофизики и электроники. E-mail: askadochkin@sv.ulsu.ru
прямое решение уравнений Максвелла, как то: метод конечных элементов (FEM) [11], метод конечных временных разностей (FDTD) [12,13], метод связанных диполей (CDA) и пр.
В данной работе нами предлагается аналитический подход к описанию оптических свойств упорядоченных нанокомпозитов, позволяющий представить агрегат из наночастиц произвольной толщины в виде математической бесконечно тонкой границы раздела, обладающей нефренелевс-кими коэффициентами отражения и пропускания. При этом, для описания оптических свойств системы, становится возможным использовать простые соотношения Эйри для нескольких границ раздела, одна из которых является «мнимой» и характеризует «сжатый» наноагрегат.
1. ОСНОВНЫЕ УРАВНЕНИЯ
Рассмотрим структуру, показанную на рис. 1 и воспользуемся методом интегральных уравнений [14-16]. Поле, создаваемое композитной средой в каждой точке пространства, может быть записано в следующем общем виде:
, P(r, t-RI c) „ E(r,t)=Ei (r, t) + j rotrot—^—--'-dV+
R
3 J e.(r')-e E.
. 3 ^f . . J\ J m je
+— > rotrot-4.
4=ir
^J ( )+2em
r,t -Rjc) , (1) --jL-kVl Г '
где Ej(r,t) = E0j exp(ik0r -itö)- внешняя волна в точке наблюдения r, k0 - волновой вектор. Первый интеграл в правой части определяет поле, создаваемое средой-матрицей с поляризацией P, пропорциональной падающему на поверхность среды полю, и диэлектрической проницаемостью 8m . При этом R = |r — r'| - расстояние от точки интегрирования r , расположенной внутри среды, до точки наблюдения, V -объем среды, с - скорость света в вакууме, аргумент (t — R / c) характеризует запаздывание со-
R
к /
/
У
/
ды является линейной функцией напряженности поля. Помещая начало координат на поверхность среды-матрицы, а точку наблюдения в центр г-й частицы, имеющий радиус-вектор ri (см. рис. 1), запишем эффективное поле:
j^fW0
!©©©©©© © © © © © ©
_© _© _© _© 1
а®©®©®©®©®©® © © © © ©
Medium © © © © © ©
©©©©©©
Рис. 1. Геометрия системы. Композит состоит из S + 1 слоёв с расстоянием Д между ними. Расстояние между поверхностью среды и плоскостью, проходящей через центры наночастиц нулевого слоя - 5
ответствующей величины. Третий член в (1) определяет поле, создаваемое J взаимодействующими наночастицами c комплексной диэлектрической проницаемостью £j (rj) и объемом Vj. При этом Rj = r — r. , г. - точка интегрирования внутри/-й наночастицы относительно начала координат. Поле Ejef , фигурирующее в выражении (1), представляет собой волну, действующую на каждую точку внутри/-й наночастицы с учетом полей, переизлученных атомами всех наночастиц. Назовем поле Eэффективным. При этом Ejef имеет две составляющие: внешнюю (или действующую со стороны окружения) и внутреннюю, определяющую взаимодействие атомов внутри самой наночастицы и отвечающую за формирование диэлектрической проницаемости среды [17]. Учет внутреннего поля приводит к разделению уравнения (1) на локальное и нелокальное; при этом первое дает известную формулу Лорентц - Лоренца. Поставленная граничная задача сводится, таким образом, к поиску эффективных полей, действующих на каждую наночастицу со стороны ее окружения (обозначим их E..^ ).
Рассмотрим композит из одинаковых однородных нанокластеров ( £ j (r) = 8 ), бесконечно протяженный в плоскости xy, внедренный в полубесконечную среду. Воспользуемся приближением длинных волн [18]:
k0а, k0Re(n)а, к0 Im(n)а 1 << 1, (2)
где n = V8, означающими, что напряженности Е/ и Ejef однородны во всем объеме сферической частицы радиуса а. При этом ограничимся случаем, когда вектор поляризации частиц и сре-
(r,/)=-
8 —£m
4яв +26
■Xrotrot j
j (( —|r —1j(c/nm.
+[ад, (°i t)+TaE, (o, t)\()_ ■
(3)
где пт = , первый член в правой части представляет собой суперпозицию полей, создаваемых наночастицами ансамбля в центре г-й. Второе слагаемое определяет полное поле (с учетом поля, излученного самим слоем и отраженного от границы раздела среда-вакуум), падающее на г-ю частицу извне; Т, Я - тензоры френелевских коэффициентов пропускания и отражения границы раздела, а порядок индексов показывает направление падения исходной волны («12» - из вакуума в среду, «21» - из среды в вакуум); Ет -падающее на границу эффективное поле, аргу-
(
мент
t —
(k о!)
а Л
ко c
характеризует запаздыва-
У
ние на время прохождения волной расстояния от начала координат до центра г-й наночастицы. Выражение для Ет примет следующий вид:
Em (0,t)
= £ — £m
4лё+2ё„
■Xrotrot j
j ((t — «j
c / nm
d
. (4)
Решая в общем виде систему из уравнений (3)-(4), получим значения эффективных полей в центре каждой из наночастиц, а также эффективного поля, падающего из среды на границу раздела. Поле, отраженное от композита, примет вид согласно (1)-(4):
Ег (К,Г) = ЯпЕ1 (0,Г)ехр(к>К) + Г21 [Ет (0,Г) (к), (5)
I V )
где R - радиус-вектор точки наблюдения, первый член в правой части характеризует френе-левское отражение от чистой (в отсутствие нано-частиц) среды, а второй - вклад нанокомпозита в отраженную волну. При этом учтено, что Ет имеет сложную структуру, отличную от плоской волны и, соответственно, запаздывание на время (к0К)/ к0с будет характеризоваться изменением и фазы, и амплитуды поля.
2. МЕТОД «МНИМОЙ ГРАНИЦЫ»
Рассмотрим поле, создаваемое^-й наночасти-цей, расположенной в вакууме, в некоторой точке наблюдения R вне ее объема. Интеграл, отвечающий соответствующей напряженности в (1),
может быть легко вычислен [14-16]. В результате получено следующее соотношение:
Ег) = *р/, (Я), (г,,г) (6)
где введены обозначения:
(Я) = ехр (,к0Я )1-2г - 2,'к°
Я3 Я2
//( Я ) = ехр (гк о Я)
(
,2 Л
++ко
ч Я3 Я2 Яу
Е-1
а = а р Е + 2
(7)
Очевидно, что если окружающая частицу среда не является вакуумом, величины Е и п в (6) -относительные, а не абсолютные, а также необходима замена к0 ^ к0Пт .
Преобразовывая интегральные слагаемые в (3), (4) в соответствии с (6)-(7), получим следующую систему уравнений для полей:
Ее (г"') = арХ I (|г,-г,|) (.г) +
j=l, з* г
Я21Ет (0, г) + Т12ЕI (0, г)"
Ет (аг) = ар£I (1 ). (г.,г).
,=1
(8)
(9)
Е ;=а р Е; Ар л + £ а р Ек С * ((;- к )А)
+
к = 0, к *;
[Я?21Ет (0,г) + Т12ЕI (0,г)}г-к^
, (10)
Ет (0,г) = ар£Ек С+ (-5-кА),
(11)
где индексы к, ; соответствуют номеру слоя. Первому слою присвоен номер 0, а также учтено, что, так как все кластеры в пределах одного слоя находятся в равных условиях, выполняется ра
венство: Е . и Е .
Е'
|=| Е
з;
=| Е ; |. Разность фаз
и учтем в соответствии с принципом трансляционной симметрии [20, 21], согласно которому напряженность поля, падающего на слой из наночастиц, удовлетворяет условию:
Е,пс (г, ) = Етс (0)ехр (,'яг;),
где ] = 1, 2... го, а вектор ц = (,qy,0), где дх = к08т01 ео8ф , qy = к08т018Шф , - угол падения, ф - угол между координатной осью х и плоскостью падения.
В выражениях (10)—(11) введены следующие обозначения:
А р; = £ I (
з =1, з * г
з;
- 5 -;А |)ехр (яг,.;) (12)
Прямое вычисление Е'е/ из системы уравнений (6)-(9) не представляется возможным, поскольку, ввиду наличия в (7) членов, пропорциональных 1/ Я, взаимное поляризующее влияние частиц является дальнодействующим [18,19]. В связи с этим, предлагается использовать метод послойного расчета полей в нано-композите, заключающийся в разделении конечного по толщине агрегата на 5 + 1 взаимодействующих монослоев [15] (см. рис. 1).
Расположим систему координат таким образом, чтобы центр рассматриваемой г - й наночас-тицы имел координаты 5 + ;А = (0,0, -8-;Д), где 5 - расстояние от поверхности до плоскости, проходящей через центры наночастиц первого слоя, Д - расстояние между слоями (рис. 1). При этом система уравнений (8)-(9) примет вид:
- решеточная сумма, описывающая поле, создаваемое наночастицами ; -го слоя в центре принадлежащей ему г-й частицы; - радиус-вектор )-й частицы ; -го слоя,
С;((;-к)А) = £/ (5 + ;А)ехр).(13)
з =1
- решеточная сумма, описывающая поле, создаваемое к-м слоем в г-й частице ; -го слоя. При этом С + имеет место при к > ; (поле излучается в положительном направлении оси 2, аргумент Ср ((; - к ) А ) - положительный вектор), а Ср
- при к < ; (поле излучается в отрицательном направлении оси 2, аргумент С р ((; - к )А) -отрицательный вектор). Вычисление решеточных сумм вида (12) и (13) проделано нами ранее [14-16] и здесь приводиться не будет. Отметим, что учет ненулевых затухающих гармоник в сумме С р необходим лишь в случае больших значений постоянных решетки, малых радиусов частиц либо больших длин волн [14-16], поэтому можно записать следующее выражение для Ср :
С Лг )= СП ехр (+,'к 0 гПт )
С = А.
2Я(к0Пт )2
|Е/
а х а
;|('к0 А )
(14)
где к0г - 2-компонента волнового вектора, а1, а2 - направляющие векторы решетки слоя. С учетом (14) выражение (5) легко привести к виду:
к =0
г =
Ег (0,г) = Д2 + Япс ехр{2шт (к05)} Е, (0,г) 1 -к21кпс ехр{ (коб)}'(15)
что в точности соответствует виду коэффициента отражения Эйри для пленки, расположенной на поверхности подстилающей среды, причем
о
Я = а V С ехр (¡к 0 (М )п )
пс р ^^ п 0 \ / т }
Е',
к = 0
Е, (0,г)
- тензор нефренелевских коэффициентов отражения ансамбля наночастиц, равный отношению напряженности поля, создаваемого им на своей верхней границе, к напряженности поля падающей волны. Следует отметить, что Ек пс - напряженность эффективного поля в &-м слое в отсутствие границы раздела, которая определяется решением укороченной системы уравнений (10) (11):
Е' = Ег (0,г) + а Е' А
' пс жи! \ ) ^ р^' пс р'
V а р Е'кпсС * ((л-к )А).
(16)
к=0, к
Аналогично записывается коэффициент пропускания системы, описывающий полное поле, прошедшее в среду-матрицу:
I =
Е, (б + SA, г) Т^ ехр{¡пп (к0б)}
Е (0,г) 1 -ККс ехр{т (к0б)}'
(17)
где тензор нефренелевских коэффициентов пропускания наноагрегата
Тпс = ехр {¡Б (к0Л)} + Яп
равен отношению напряженности поля, созданного агрегатом наночастиц на своей нижней границе, к напряженности падающего поля.
Таким образом, упорядоченный наноагрегат, внедренный в среду-носитель, может быть трактован как некая мнимая граница раздела двух сред, бесконечно тонкая и обладающая сложными нефренелевскими коэффициентами отражения и пропускания. Другими словами, рассматриваемая система может быть представлена в виде пленки, толщиной 5 (расстояние от поверхности среды-матрицы до мнимой границы), расположенной на поверхности полубесконечной подстилающей среды, хотя в действительности и пленка, и подложка состоят из одного и того же материала и, соответственно, обладают одинаковым показателем преломления. Следует отметить, что ранее [22] предпринимались попыт-
ки трактовать слой наночастиц (либо упорядоченный агрегат из нескольких слоев) как пленку конечной толщины с неким эффективным показателем преломления. Однако, ввиду отсутствия четких границ слоя, которые можно было бы считать границами данной пленки, толщина последней определяется в различных источниках по-разному. Данное затруднение легко понять, если представить монослой не из сферических частиц, а из неких мицеллярных нанообразований, представляющих собой ядро, покрытое иглообразными наростами. При этом неочевидно, что следует понимать под толщиной данной гетерогенной пленки: следует ли ее определить как диаметр ядра либо необходимо включить в пленку также и мицеллярные «иглы». Поскольку от того, как выбраны границы композитной пленки, зависит фактор заполнения ее наночастицами, эффективный показатель преломления является функцией толщины, что влияет на результирующие расчетные оптические характеристики системы.
В рамках рассматриваемого формализма указанные затруднения отсутствуют, поскольку понятия толщины гетерогенной пленки, представляющей собой слой наночастиц, а также ее эффективного показателя преломления теряют смысл. Выражения (15), (17) весьма удобны, поскольку позволяют, единожды рассчитав эффективные поля в наноагрегате из простой системы алгебраических уравнений (16), использовать хорошо проработанный в классической оптике многослойных покрытий формализм Эйри для решения более сложных задач, касающихся исследования электродинамических свойств трехмерных многоуровневых композитных образований как в полубесконечных, так и в слоистых средах. Отметим, что если композит внедрен в пленку конечной толщины, выражения (15), (17) необходимо преобразовать в формулы Эйри для трех границ раздела.
3. СРАВНЕНИЕ С ТОЧНЫМ ЧИСЛЕННЫМ РАСЧЕТОМ
Предложенный подход апробируем сравнением с точным численным решением уравнений Максвелла методом конечных элементов. На рис. 2 показаны расчетные спектральные зависимости коэффициента отражения двух стеклянных пленок, в каждую из которых внедрены три упорядоченных монослоя наночастиц. Очевидно, что полученные результаты находятся в хорошем согласии со строгим решением, что подтверждает применимость и точность предложенного метода моделирования оптических свойств наноком-позитных пленок. Более того, поскольку в рамках метода конечных элементов нановключения
%%
Рис. 2. Спектр отражения пленки с тремя слоями наночастиц. Точками показаны результаты
численного решения, сплошной линией -расчет в рамках предлагаемого метода. Параметры системы: 1 - толщина пленки
h = 100 нм , а = 10 нм, |a 12| = 3а , д = 4a ; 2 -
h = 180 нм , а = 15 нм, |a 12| = 4а , д = 5а .
Среда-матрица - стекло: nm = 1.5 , показатель преломления наночастиц n = 2.5
характеризуются не микроскопическими величинами, как это имеет место в уравнении (1), а макроскопической диэлектрической проницаемостью материала, из которого они состоят, согласие расчетных данных также подтверждает верность соответствующего преобразования (6), (7), позволяющего перейти к макроскопическому описанию как диэлектрических, так и проводящих наночастиц.
Следует отметить, что межчастичное расстояние в рассмотренных упорядоченных агрегатах таково (см. подпись к рис. 2), что в разложении рассеянного частицами поля (7) существенную роль играют члены, пропорциональные ik0 / R2 и 1/ R3, а также запаздывающее электродинамическое взаимодействие, которые, как следует представленных на рис. 2 результатов, адекватно учитываются в рамках предлагаемой модели.
ЗАКЛЮЧЕНИЕ
Таким образом, в данной работе были исследованы оптические свойства упорядоченных на-ноагрегатов из сферических диэлектрических наночастиц. Предложен аналитический метод «мнимой границы» для вычисления оптический параметров нанокомпозитных систем. Показано, что подобные структуры могут быть описаны формулами Эйри, применяемыми обычно для расчёта характеристик многослойных систем и покрытий. Рассматриваемый агрегат из наночастиц играет роль мнимой границы бесконечно малой
толщины со специфическими нефренелевскими коэффициентами отражения и пропускания, в которых учитывается взаимодействие частиц друг с другом и с материалом матрицы. Рассматриваемый метод «мнимой границы» не требует усреднения оптических констант по объёму наноагре-гата (гомогенизации) и является микроскопическим. Показано хорошее согласие расчетных оптических характеристик системы с результатами применения точного численного метода.
Исследование выполнено при финансовой поддержке РФФИ в рамках научных проектов № 12-02-97036-р_поволжье_а; № 12-02-31423-мол_а, а также при поддержке Министерства образования и науки РФ в рамках ФЦП «Кадры» 20122013 г.
СПИСОК ЛИТЕРАТУРЫ
1. Shalaev V.M., Cai W., Chettiar U.K., Yuan H. -K., Sarychev
A.K., Drachev V.P., andKildishev A.V. Negative index of refraction in optical metamaterials. Opt. Lett. 2005. 30. 3356-3358.
2. Smith D.R. and Pendry J.B. Homogenization of Metamaterials by Field Averaging // J. Opt. Soc. Am.
B. 2006. 23. 391-403.
3. Agranovich V.M., Shen Y.R., Baughman R.H., and Zakhidov A.A. Linear and nonlinear wave propagation in negative refraction metamaterials // Phys. Rev. B. 2004. 69. 165112.
4. O.N. Gadomskii and A. S. Shalin // J. Exp. Theor. Phys. 105, 761 (2007).
5. Grigorenko N., Geim A. K., Gleeson H. F., and Zhang Y. Nanofabricated media with negative permeability at visible frequencies // Nature. 2005. 438. 335-338.
6. Experimental demonstration of Near-Infrared NegativeIndex Metamaterials / S. Zhang, W. Fan, N.C. Panoiu, K.J. Malloy, R.M. Osgood, and S.R.J. Brueck // Phys. Rev. Lett. 2005. 95. 137404.
7. Very low-refractive-index optical thin films consisting of an array of SiO2 nanorods / J.-Q. Xi, J. K. Kim, E. F. Schubert, D. Ye, T. M. Lu, S.-Y. Lin, and J. S. Juneja. // Opt. Lett. 2006. 31. 601-603.
8. Mishchenko M.I., Travis L.D., Lacis A.A. Scattering, absorption and emission of light by small particles. Cambridge: Cambridge University Press, 2002. 457 p.
9. Haarmans M.T., Bedeaux D. The polarizability and the optical properties of lattices and random distributions of small metal spheres on a substrate. Thin Solid Films. 1993. 224. 117-131.
10. Mie G. Beiträge zur Optik trüber Medien, speziell kolloidaler Metallösungen // Ann. Phys. 1908. 25. 377.
11. Arfken G.B., Weber H.J. Mathematical Methods for Physicists. New York: Academic Press, 1995. 1195 p.
12. TafloveA., HagnessS.C. Computational Electrodynamics: The Finite-difference Time-Domain Method. - Boston: Artech House, 2000. 886 p.
13. Prather D.W., Shi S. Formulation and application of the finite-difference time-domain method for the analysis of axially symmetric diffractive optical elements // J. Opt. Soc. Am. A. 1999. 16. 1131-1142.
14. Shalin A.S., Moiseev S.G. // Opt. Spectrosc. 2009. V.106. P. 916.
15. Shalin A.S. // J. Commun. Technol. Electron. 2009. V. 54. P. 699.
16. Shalin A.S. // Progress in Electromagnetic Research B. 2011. V. 31. P. 45.
17. Борн М. Вольф Э. Основы оптики. Москва: Наука, 1973. 720 с.
18. Борен К., Хаффмен Д. Поглощение и рассеяние света малыми частицами. Москва: Мир, 1986. 664 с.
19. Mishchenko M.I., Travis L.D., and Lacis A.A. Scattering, Absorption and Emission of Light by Small Particles.
Cambridge: Cambridge University Press, 2002. 351 p.
20. Wijers M.J. and Poppe G.P.M. Microscopic treatment of the angular dependence of surface induced optical anisotropy // Phys. Rev. B. 1992. 46. 7605-7620.
21. Poppe G.P.M., Wijers C.M.J. and Silfhout A. IR spectroscopy of CO physisorbed on NaCl(100) // Microscopic treatment. Phys. Rev. B. 1991. 44. 79177929.
22. Milton G.W. The Theory of Composites. Cambridge: Cambridge University Press, 2004. 420 p.
HIGH TRASPARENCY NANOMATERIAL
© 2012 A.S. Shalin1, 2, A.S. Kadochkin2
1 Ul'yanovsk Branch of Institute of Radio Engineering and Electronics named after V.A. Kotel'nikov of RAS 2 Ul'yanovsk State University
Optical properties of spatially ordered nanocomposites consisting of spherical dielectric nanoparticles have been investigated. We proposed an analytical solution allowing us to use simple Airy formulas ("imaginary boundary" method) to describe interaction of a light with an ordered nanostructure without averaging of optical constants over the volume. It has been shown a good agreement between the results obtained with the exact numerical calculations Keywords: metamaterial, Airy formulae, imaginary boundary.
Alexandr Shalin, Candidate of Physics and Mathematics, Senior Research Fellow at the Ul'yanovsk Branch of IRE RAS, Senior Research Fellow at the Technological Institute. E mail: shalin_a@rambler.ru
Alexey Kadochkin, Candidate of Physics and Mathematics, Senior Research Fellow at the Tech-nological Institute. E mail: askadochkin@sv.ulsu.ru
