УДК 681.7.064.454
О. Н. Гадомский, К. К. Алтунин, Е. Г. Зубков РАДИАЦИОННАЯ ТЕОРИЯ МЕТАЛЛИЧЕСКОГО КЛАСТЕРА
Аннотация. Представлена теория металлического сферического кластера, в которой учтены электрические дипольные квантовые переходы валентных электронов с помощью эффективной поляризуемости. Показано, что эффективная поляризуемость валентных электронов металлического сферического кластера обладает отрицательной дисперсией в широком диапазоне длин волн от 450 до 1000 нм, а это приводит к отрицательной экстинкции металлической сферы.
Ключевые слова: сферический металлический кластер, эффективная поляризуемость валентных электронов, отрицательная дисперсия эффективной поляризуемости, локальные поля.
Abstract. The article introduces a theory of spherical metal cluster, which considers electric dipole quantum transition of valence electrons by means of effective polar-izability. It is shown that effective polarizability of the valence electrons of spherical metal cluster possesses negative dispersion in a broad wave band from 450 to 1000 nm, that leads to negative extinction of a metal sphere.
Key words: spherical metal cluster, effective polarizability of the valence electrons, negative dispersion of the effective polarizability, local fields.
Введение
Оптические свойства металлических сфер описаны в теории Ми [1], в которой на основе уравнений макроскопической электродинамики вычисляются электромагнитные поля внутри и вне сферы.
В данной статье, в отличие от традиционного теоретического подхода [1], будут использованы полевые и материальные уравнения микроскопической электродинамики, что позволит обнаружить новые оптические свойства металлических сфер.
В работе [2] представлены экспериментальные спектры сечения рассеяния изолированных сфер серебра. Мы используем эти данные, чтобы вычислить параметры радиационной теории кластера. Будут получены интерполяционные формулы для частоты квантового перехода, дипольного момента перехода и ширины резонанса валентных электронов в металлической сфере.
Будет получена формула для эффективной поляризуемости валентных электронов в изолированной сфере, которая содержит область отрицательной дисперсии при некоторых радиусах металлических сфер.
Будет показано, что с помощью эффективной поляризуемости валентных электронов может быть вычислен показатель преломления оптической среды из металлических сфер, что позволяет обнаруживать новые свойства композитных материалов, не содержащиеся в теории Максвелла - Гарнетта [3].
1. Уравнения микроскопической электродинамики для металлической сферы в вакууме
Уравнение распространения оптической волны в рассматриваемой граничной задаче имеет следующий вид [1, 4]:
E(г, і) = EI (г, і) +1 ю ю N р(г , ^ dV',
(1)
V
где Е(г, ^) - напряженность электрического поля в точках наблюдения г внутри и вне металлической сферы; Е1 - напряженность электрического поля внешней волны; N - концентрация валентных электронов внутри сферы;
в вакууме; р - индуцированный дипольный момент валентного электрона; V - объем сферической частицы. Дифференцирование в (1) проводится по координатам точки наблюдения.
Уравнение (1) было выведено на основе квантовоэлектродинамического и классического подходов [4]. Будем рассматривать металлическую сферу как непрерывную частицу, пренебрегая пока эффектом ближнего поля, связанного с учетом дискретного распределения атомов [4].
Квантовомеханические средние индуцированных дипольных моментов валентных электронов в кластере представим как
где ю - частота внешнего поля; величина Х удовлетворяет уравнениям для связанных квантовых диполей [5, 6]:
где А = ю0 - ю, ю0 - частота перехода валентного электрона из основного в возбужденное состояние; d0 - дипольный момент перехода; Т2, Ті - времена фазовой и энергетической релаксации; w - инверсия, представляющая собой разность вероятностей обнаружения валентных электронов в основном и возбужденном состояниях; w0 - равновесное значение инверсии, равное - 1, Ео - локальное поле без множителя ехр(—/ю • і).
Найдем стационарное решение уравнений (3) при выполнении следующих условий:
Более того, предположим, что оптическое поле, действующее на валентные электроны в сфере, невелико. Условие малоинтенсивного поля может быть выражено как w ~ w0, т.е. инверсия практически не меняется при облучении сферической частицы. Тогда из уравнения (3а) получим следующее решение:
Я = |г — г |, г - точка внутри сферы или на ее поверхности; с - скорость света
р = X ехр(—ІШ),
(2)
(За)
(36)
w = 0, X = 0.
(4)
р = аЕ0 ехр(—/юі) = а^ Е01 ехр(—іюі),
(5)
где квантовая поляризуемость равна
а =
2d02
h i ю0 -ю-
(6)
T2
aejf представляет собой эффективную поляризуемость валентных электронов внутри сферической частицы.
Для определения явного вида эффективной поляризуемости необходимо решить уравнение (1) для точек наблюдения внутри частицы.
2. Эффективная поляризуемость валентных электронов сферического металлического кластера
Вычислим локальное поле внутри сферической частицы, например, в ее центре.
С помощью математической леммы [1] вынесем оператор rot rot за знак интеграла в уравнении (1). Тогда получим следующее равенство:
J rot rot Nr d V' = exp(-irot) • aTNX =
= exp(-irot)(n2 -1)k0
f* 8n
rotrot J QGdV' --3- Q
V
(7)
ю
где ko = —, Q - некоторая функция координат, определяемая как
Np = (n2 - 1)k2Q, V2Q + k02n2Q = 0 ;
(8)
G(Я) = 6 0—) - функция Грина; ат - внутренний геометрический фак-
Я
тор. Величина щ в (7), (8) представляет собой показатель преломления вещества из структурных элементов с поляризуемостью а.
Учитывая свойства функций Q и G, преобразуем с помощью теоремы Грина объемный интеграл в (7) в поверхностные интегралы по внешней поверхности Е и внутренней поверхности о, окружающей точку наблюдения. Тогда получим
GdV =
1
(П2 - 1)ko2
J( Q !G - G IQ У - J( Q§
-G— dS' dR
(9)
d
где символ —- означает дифференцирование вдоль внешней нормали к по-Эv
верхности Е.
Подставим функцию Q, удовлетворяющую волновому уравнению (8), следующим образом:
Q = Qo
sin(ko nR R
(10)
1
c
№ 3 (23), 2012 Физико-математические науки. Физика
где Qo - постоянный вектор. Тогда поверхностный интеграл по внешней поверхности примет следующий вид:
^ Q 17 - ° ^=т вто, (°)
где
B = iko sin(konR eik(>R - kon cos(ko«1R) eik(>R; (12)
R - расстояние от поверхности E до точки наблюдения r внутри сферы.
Вычисление поверхностного интеграла по сферической поверхности о исчезающе малого радиуса приводит к следующему выражению:
q1Tr - gTr )dS'=4nB(R)• (l3)
О
где R' - расстояние от поверхности о до точки наблюдения r .
Подставим (13), (11) в выражение (9), а затем в (7). После вычисления действия оператора rot rot по координатам точки наблюдения r получим следующую формулу:
4п 4п ikoa
aT =------------------e o
( k sin^ko a)
cos(«1ko a)------------
n^o a
. 4П eik0a sin(«1koa) , (14)
3 «1
где а - радиус сферического кластера.
При выполнении условия к0а << 1 имеем
аТ =~у(1 + ко а). (15)
В электростатическом приближении (с ^ да) внутренний геометрический фактор ат = —4“, что соответствует деполяризующему фактору для сферы [7].
Теперь, определив явный вид внутреннего геометрического фактора, подставим (7) в уравнение (1). Тогда в точке наблюдения г = 0 в центре сферы получим следующее равенство:
Ео = Ео I + атМ аЕо, (16)
откуда с помощью соотношения (5) получим формулу для эффективной поляризуемости валентных электронов сферического кластера:
а*# = 1 а М • (17)
м 1 — атМ а
Как следует из (16), локальное поле внутри кластера может сильно отличаться от внешнего поля. При этом дисперсионные свойства эффективной поляризуемости сильно отличаются от дисперсионных свойств квантовой поляризуемости а = Яе а + / 1т а . Например, 1т ае^ может принимать отрицательные значения при малых радиусах сферы, что соответствует отрицатель-
ной дисперсии эффективной поляризуемости, валентных электронов металлического кластера.
3. Оптическое поле вне металлического кластера
Для вычисления локального поля вне металлического кластера поместим в уравнение (1) точку наблюдения вне сферы и проведем вычисления по изложенному выше методу. Отличие заключается в том, что оператор rot rot может быть вынесен за знак интеграла без применения математической леммы. После вычислений получим следующее выражение для внешнего геометрического фактора aR:
v 4пА
7 =----------
R щкоП -1) IR3
2ik
z = -
R
4пА
niko(n2 -1)
R2
ik0
,2
0 I eikoR
(18)
R R
k0
R
2 I
ik0R
где R - расстояние от центра сферы до точки наблюдения вне сферы,
А = со$,(к0а)$,т(п1коа) - п$т(коа)со$(п1коа) • (19)
С помощью внешнего геометрического фактора cІR можно вычислить локальные поля Е0 в любой точке наблюдения вне сферы в волновой и ближней зонах.
В волновой зоне, когда R >> X, сечение рассеяния QS вычислим с помощью эффективной поляризуемости (17). В результате для сечения рассеяния металлической сферы получим следующую формулу:
Qs =-
4лкХ2^2
(20)
T2
где Го = — а3, По = N 2 0 , П = ПоаТ .
3 Й
С помощью формулы (2о) и экспериментальных спектров рассеяния изолированных сфер серебра [2] найдем следующие интерполяционные формулы для параметров электрических дипольных квантовых переходов:
■л-12
2пс
ю0
= 268,6 • 10-7 + 4,412a +
2,904-10"
■845172,1a2 + 5,669 • 1010a3;
ldp| =
eaB
0,404+
57782,5 • a - 0,298
|845069,4 • a - 4,3502|
0,758
1 = 2|d0\2 ю3
T2
:&1
3 he3
(21)
где gl = N4- а3 - число валентных электронов в сферической частице серебра, N = 5,8 • 1022 см -3, аВ - боровский радиус, аВ = 0,529 • 10 -8 см. В экспери-
ментальных спектрах рассеяния [2] исследовались сферические частицы серебра с радиусами от 20 до 80 нм. С помощью интерполяционных формул (21) мы имеем возможность вычислять параметры электрических дипольных квантовых переходов в радиационной теории металлического кластера при различных радиусах сферических частиц серебра.
Экстинкцию сферы в вакууме вычислим по следующей формуле [1]:
волны (е = I); п0 - единичный вектор вдоль направления распространения линейно поляризованной внешней волны, а(п0) = Nае^е . Очевидно, что Q = Qs + Qa , где Qa - сечение поглощения изолированной сферы [1].
На рис. 1, 2 представлены дисперсионные зависимости эффективной поляризуемости валентных электронов в наночастицах серебра. Видно, что при малых радиусах наночастиц в широком диапазоне длин волн от 450 до 1000 нм образуется область отрицательной дисперсии. Как показывают численные исследования эффективной поляризуемости (17), при радиусах наночастиц а > 12,5 нм знак дисперсии изменяется на противоположный.
(22)
где X - длина волны внешнего излучения, X =
2п-с
; е - амплитуда внешней
^ 9,2
а. нм
X. 10~7 м
9,2
Рис. 1. Реальная часть эффективной поляризуемости валентных электронов в сферических наночастицах серебра как функции длины волны и радиуса наночастиц
а. нм
Рис. 2. Мнимая часть эффективной поляризуемости валентных электронов в сферических наночастицах серебра как функции длины волны и радиуса наночастиц
В области отрицательной дисперсии Im aef < 0 экстинкция становится
также отрицательной согласно формуле (22). Это означает, что при Q < 0 и Qs > 0 имеем отрицательное поглощение Qa < 0.
Заключение
Итак, в данной статье представлена радиационная теория металлического кластера на основе уравнений (1), (3). Решена граничная задача для изолированной сферы, найдены локальные поля внутри и вне сферы, получена формула для эффективной поляризуемости валентных электронов. На основе экспериментальных спектров рассеяния изолированных частиц серебра выведены интерполяционные формулы для параметров электрических ди-польных квантовых переходов, в которых участвуют валентные электроны металлического кластера.
Показано, что при малых радиусах наночастиц серебра в этих частицах формируется область отрицательной дисперсии эффективной поляризуемости валентных электронов, что приводит к образованию в этих частицах отрицательной экстинкции и отрицательного поглощения.
В работе [8] нами сообщалось о синтезе новых композитных материалов (PMMA + Ag) по разработанной в нашей лаборатории технологии. Измерение оптических спектров структур (PMMA + Ag) / glass с толстыми пленками толщиной 80 мкм из этих материалов показало наличие интерференционных максимумов. Наличие этих интерференционных максимумов объясняется, с нашей точки зрения, тем, что в этих материалах формируются квази-нулевые показатели преломления и амплитуды отражения и пропускания представляются как когерентные суперпозиции амплитуд отражения и пропускания с показателями преломления в окрестности его нулевого значения. Для моделирования показателя преломления композитного материала (PMMA + Ag) в [8] применяется соответствующая формула смешения [9] и
показано, что в системе наночастиц серебра формируется отрицательный показатель преломления, обусловленный отрицательной дисперсией эффективной поляризуемости валентных электронов в наночастицах серебра.
Список литературы
1. Born, M. Principle of Optics / M. Bom and E. Wolf. - Cambridge : Cambridge University Press, 1999.
2. Tamaru, H. Resonant light scattering from individual Ag nanoparticles and particle pairs / H. Tamaru, H. Kuwata, H. T. Miyazaki, K. Miyano // Applied Physics Letters. -2002. - V. 80, issue 10. - P. 1826-1828.
3. Garnett, J. C. M. Colours in Metal Glasses and in Metallic Films / J. C. M. Garnett // Philos. Trans. R. Soc. London, Ser. A. - 1904. - V. 203. - P. 385-420.
4. Гадомский, О. Н. Проблема двух электронов и нелокальные уравнения электродинамики / О. Н. Гадомский // Успехи физических наук. - 2000. - Т. 170, № 11. - С. 1145-1179.
5. Гадомский, О. Н. Эффект ближнего поля в квантовом компьютере / О. Н. Гадомский, Ю. Ю. Воронов // Письма в Журнал экспериментальной и теоретической физики. - 1999. - Т. 69, № 10. - С. 750-754.
6. Гадомский, О. Н. Оптические размерные резонансы в наноструктурах / О. Н. Гадомский, Т. Т. Идиатуллов // Журнал экспериментальной и теоретической физики. - 2001. - Т. 119, № 6. - С. 1222-1234.
7. Киттель Ч. Введение в физику твердого тела : пер. с англ. / Ч. Киттель. -2-е изд., перераб. - М. : Физматлит, 1962. - 696 с.
8. Gadomsky, O. N. High-negative effective refractive index of silver nanoparticles system in nanocomposite films / O. N. Gadomsky, K. K. Altunin // Optics Communications. - 2012. - V. 285, № 5. - P. 816-820.
9. Voshchinnikov, N. V. Effective medium theories for irregular fluffy structures: aggregation of small particles / N. V. Voshchinnikov, G. Videen, T. Henning // Applied Optics. - 2007. - V. 46, № 19. - P. 4065-4072.
Гадомский Олег Николаевич
доктор физико-математических наук, профессор, кафедра радиофизики и электроники, Ульяновский государственный университет
E-mail: [email protected]
Алтунин Константин Константинович
кандидат физико-математических наук, доцент, кафедра физики, Ульяновский государственный педагогический университет
E-mail: [email protected]
Зубков Евгений Геннадьевич
аспирант, Ульяновский государственный университет
E-mail: [email protected]
Gadomsky Oleg Nikolaevich Doctor of physical and mathematical sciences, professor, sub-department of radio physics and electronics, Ulyanovsk State University
Altunin Konstantin Konstantinovich Candidate of physical and mathematical sciences, associate professor, sub-department of physics, Ulyanovsk State Pedagogical University
Zubkov Evgeny Gennadyevich Postgraduate student, Ulyanovsk State University
УДК 681.7.064.454 Гадомский, О. Н.
Радиационная теория металлического кластера / О. Н. Гадомский, К. К. Алтунин, Е. Г. Зубков // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2012. - № 3 (23). - С. 144-152.