Кучеренко М.Г.
Центр лазерной и информационной биофизики Оренбургского государственного университета
E-mail: rphys@mail.osu.ru
ДИНАМИЧЕСКАЯ ПОЛЯРИЗУЕМОСТЬ НАНОШАРА В СЛУЧАЕ ВЫРОЖДЕННОГО ЭЛЕКТРОННОГО ГАЗА И ЕЕ РОЛЬ В ПЛАЗМОННОМ МЕХАНИЗМЕ ПЕРЕДАЧИ ЭНЕРГИИ
Показано, каким образом вырожденность электронного газа в сферической металлической наночастице сказывается на резонансном поведении ее дипольной поляризуемости в монохроматическом электромагнитном поле на оптических частотах при наличии диссипативных процессов в проводнике (комплексности диэлектрической проницаемости). Частотные зависимости действительной и мнимой частей поляризуемости сферической наночастицы получены в виде графиков для различных значений радиуса частиц, длины томас-фермиевского экранирования и постоянной затухания. Представлены выражения для скорости безызлучательного переноса энергии между молекулами вблизи наночастицы и скорости передачи энергии к ней, включающие в себя динамическую поляризуемость наносферы с вырожденным электронным газом.
Ключевые слова: поляризуемость сферической наночастицы, вырожденный электронный газ, томас-фермиевское экранирование, плазмоны, перенос энергии
Влияние нанотел на радиационные и бе-зызлучательные переходы в молекулах оказалось весьма существенным, о чем неоднократно сообщалось разными авторами в ряде работ [1-15]. В металлооптике рассеяние света сферическими частицами изучалось Ми [16-17], причем в то время электронный газ проводника рассматривался как классический, а потенциалы электромагнитного поля, внутри частицы-глобулы, удовлетворяли соответствующим уравнениям электродинамики.
В духе современных представлений квантовые свойства электронного газа (ферми-газа) типичны для металлов вплоть до температур их плавления и для полупроводников - при низких температурах [18]. Функция распределения электронов по энергиям для такого - вырожденного газа, в равновесных условиях является ферми-дираковской, а распределение электронной плотности п(г) в неоднородных системах - металлических кластерах, может быть найдено из решения уравнения Томаса-Ферми для потенциала ф(г) электрического поля [19]. В этой связи представляет интерес определение степени влияния вырожденности электронного газа наночастиц на безызлучательную передачу энергии электронного возбуждения от молекул-доноров к акцепторным центрам возле металлических наночастиц, включая и случай выполнения этими частицами функций тушителей возбужденных состояний. Безусловно важным, при этом, оказывается установление плазмонных резонансов дипольной поляризуемости наноглобул с учетом диссипативных процессов в металле.
Для нейтрального кластера радиуса Я, с числом атомов N>>1, уравнение Томаса-Ферми для потенциала ф(г) в «модели желе» приобретает в правой части слагаемое, представляющее однородную плотность положительных ионов металла, и при г < Я оно может быть записано в виде[20-21]
1
dr
dr
4e[2meф(r)f// 2 3eN
3nh3
R3
(1)
В большей части ионной сферы этот потенциал постоянен ф = ф0 и определяется числом ионов N и радиусом Я кластера [21]
. _(9^)2 73 Й2
27 / 3 emeR2
(2)
Линеаризация уравнения Томаса-Ферми (1) проведенная для малого монохроматического возмущения 5ф(г,ґ) электростатического потенциала ф —> ф 0 + §ф, 5ф<< ф0 приводит к следующему результату с использованием атомных единиц ( те = е = Н = 1) при г < Я [21]
(3)
У28ф = ±,/2Ф08ф, Фо = (9^
пе 2773 Я
Уравнение (3) известно в математической физике как уравнение Гельмгольца, и оно появляется как базовое в теории рассеяния электромагнитных волн, однако постоянная 4
k2 =-
пе
фигурирующая в (3), носит совершенно иной физический смысл, нежели квадрат волнового числа в классическом уравнении Гельмгольца.
2
r
2
r
Спектр собственных плазмонных колебаний в наночастице
При фиксированной плотности ионов в глобуле потенциал ф0 не зависит от радиуса кластера [21]. В отсутствие внешнего поля диэлектрическая проницаемость е(О), обусловленная наличием свободных электронов в кластере, определяется на частоте О собственных колебаний электронной плотности в наноглобуле как действительная функция аргумента О
е(О) = 1 -(р/о)2 <0. (4)
Здесь юр = д/4ле2пе/ш - плазменная частота; пе - концентрация электронов. Введем обозначение
к2 =-—72ф0 = -
4
к2 = —
пе ' '" пе/2
4т3 / 2е5 /2
> 0
или
пй3е
л/2ф0 > 0,
(5)
где / = (2ф0)—/4 - длина томас-фермиевского экранирования.
Тогда уравнение (3) в сферических координатах при г<Я принимает вид
1 д
д —<
дг
1
д
в дв
этв
дв
= —к
<Ф.(6)
Решение уравнения (6) в форме собственных дипольных колебаний
8ф(г, 0 = и(г)р(соз0)соз Оt. (7)
Здесь Р1 (cos0) - полином Лежандра степени I.
Подставляя (7) в (6) получаем обыкновенное дифференциальное уравнение для функции и(т) (уравнение Бесселя)
1 А
г2 Аг
. йы йг
(8)
Его решение, регулярное в начале координат, имеет вид
и(г) = А(кт)~1/2]з/2(кт),т< Я, (9)
где ]г+1 /2(г) - функция Бесселя, причем для диполя /=1. Константа А в (9) определяет амплитуду собственных колебаний плотности вырожденного электронного газа.
Решение для возмущающей добавки Ьф(т^) потенциала вне ионной сферы т > Я строится на основе уравнения Лапласа
1 " ■ " (10)
1 д Г 2 д х і
72 д Г 57 5ф|+
2 ІІйвїІ^ 0
г яивдві дв
Решение уравнения (10), регулярное на бесконечности можно записать в виде
8ф(г, г) = -^2 Р(соэв)соэ Ш. (11)
Константа В в (11) представляет собой амплитуду осциллирующего дипольного момента р (ґ) нанокластера. Из условий непрерывности потенциала 5ф(г,ґ) на границе ионной сферы г=Я, а также непрерывности нормальной проекции вектора электростатической индукции на основе (9) и (11) находим неявное выражение для собственной частоты Ш диполь-ных колебаний электронного облака кластера в приближении Томаса-Ферми [21]
2(кЯ ^ кЯ—1)
(кЯ)2
(12)
По сути, с учетом (4) и (5), выражение (12) представляет собой трансцендентное уравнение относительно О . Гадомский и Шалин [22] сообщают о наборе корней (12), находящихся при определенных радиусах Я нанокластера в видимой спектральной области. Оказалось, что спектр собственных частот О чувствителен к изменениям радиуса Я как в отношении области локализации дискретных уровней, так и в отношении их взаимного расположения. В пределе очень малого кластера кЯ<<1 из (12) следует ом = юр/43, то есть из всего спектра остается единственная хорошо известная частота Ми поверхностных ди-польных колебаний малого шара.
На рис. 1. представлены вычисленные на основе (12) зависимости собственных частот О различных ветвей плазмонных колебаний наношара от его радиуса Я, выраженного в длинах томас-фермиевского экранирования I = (2ф0)_1774. Нижняя ветвь условно названа нами «ветвью частот Ми». При кЯ ^ 0 для нее характерно значение Ом = юр /л/э. Более высокочастотные ветви появляются лишь при больших значениях Я/1. На рис. 1 показаны две из них. Для собственных частот о всех ветвей характерно монотонное уменьшение их значений с ростом параметра Я /I. Переход к невырожденному электронному газу осуществляется при I ^ ^, и тогда нижняя ветвь вырождается в прямую ( Ом/ Юр)2 ^ 1 / 3, а верхние ветви исчезают.
Из (12) и рис. 1-2 следует, что с ростом радиуса кластера И. частота О нижней ветви убывает по сравнению с частотой Ми ом = юр /-Уэ. Для золота концентрация электронов пе = 5.9 1022 см-3
+
[22], плазменная частота юр = 13.7 1015 рад/с, k2 =-3.17 ■ 1016 е-1 см-2. Тогда для частот оптического диапазона и нанокластера радиусом R > 1 нм выполняется неравенство kR>>1 и явная зависимость от радиуса кластера R частоты Q из (12) существенна.
На рис 2. представлена зависимость квадрата (Q / юр)2 собственной частоты колебаний плотности вырожденного электронного газа в единицах плазменной частоты от радиуса R наночастицы в единицах длины l томас-фермиевс-кого экранирования, а также зависимость функции kR от параметра R/l. Из графика kR видно, что имеет место нелинейная связь томас-ферми-евских параметров k и 1/1 друг с другом. В следующем разделе будет показано, что всем собственным частотам Q локальных плазмонных колебаний отвечают резонансы дипольной динамической поляризуемости нанокластера, причем появление новых ветвей при больших значениях параметра kR>>1 делает спектр поляризуемости все более сложным, в отличие от случая малых kR<< 1, когда электронный газ становится практически невырожденным и в спектре остается одна частота QM = юр / V3 - частота Ми.
Динамическая поляризуемость кластера с учетом вырожденности электронного газа
Рассмотрим теперь вынужденные диполь-ные колебания, вызываемые внешним квазиод-нородным электромагнитным (лазерным) полем, напряженностью E0cos ot. Решение для возмущающей добавки 8ф(r,t) потенциала вне ионной сферы r>R строится в виде
a (о)
<ф(г, 0 = Е0
-r+
P (cos в) cos (Ot, (13)
где a(о) - динамическая поляризуемость кластера на частоте поля. Решение для возмущающей добавки S$(r,t) внутри сферы (r<R) представляем в форме вынужденных дипольных колебаний
8ф(г, t) = C (kr )-1/2 J3/2(kr )P (cos в) cos ot. (14)
Используя обозначение (5) и выражение для диэлектрической проницаемости е(о) электронного газа в форме (4) для функции k2 (о) получаем
4
kZ(O) =------= І 24Ю 2\,2 ,
пе(ю) n(COp -Ю )l)
(15)
где, по-прежнему, l = (2ф0 )_1 / 4 фермиевского экранирования.
длина томас-
Из условий непрерывности потенциала <ф(г,ґ) на границе ионной сферы г=И и непрерывности нормальной проекции вектора электростатической индукции, на основе (13) и (14) для поляризуемости а(а) получаем [21]
a( ю) = R‘
1+-
3(1-kR ctg kR)
. (16)
2( е( ю) -1)(1- kR ctg kR) - е( o)(kR)2 _
При kR<<1 из (16) следует известный классический результат [23]
е(ю) -1
а(ю) ^ — Я (17)
Е(Ю) + 2 • V
Из выражения е( ю)=1-( юр/ю)2 следует е(Юр)=0 и тогда для поляризуемости а(ю) на плазменной частоте получаем а(юр) = -Я / 2.
При совпадении частоты ю с частотой О собственных колебаний электронного газа на-
Радиус наночастиц, Я/1
Рисунок 1. Зависимости квадратов собственных частот о различных ветвей плазмонных колебаний наношара от его радиуса И, выраженного в длинах томас-фермиевского экранирования I = (2ф0 )-1 /4.
Радиуснаночастиц, Я/1
Рисунок 2. Зависимость квадрата собственной частоты О нижней ветви плазмонных колебаний наношара от его радиуса И, выраженного в длинах I=(2ф0 )-1/4 томас-фермиевского экранирования. Зависимость функции кИ от параметра И/1, реализована посредством неявного
О / ю.
2 ( R
соотношения kR = I
1 - (О / Юр)
ступает резонанс и поляризуемость а(о) обращается в бесконечность. Условие равенства нулю знаменателя во втором слагаемом (16) приводит к уравнению (12). Рассмотрим ту часть формулы (16), которая отвечает за резонансный эффект
_______3(1- kR ctg kR)/(kR)2_____=
2(е(о) -1)(1 - kRctg kR)/(kR)2 - е(о) ~
= 3(1 - kR ctg kR) ю2
(кЯ)2 ю2 -О2
Вблизи собственной частоты ю«О резонансный характер динамической поляризуемости отображается зависимостью
2
О
юр
V 1 /
ю2р -О2 . (18)
О2-ю2
Для устранения расходимости а(ю) на резонансной частоте рассмотрим, далее, диссипативный вариант плазменной системы. Тогда диэлектрическая проницаемость е( ю) вместо (4) может быть выбрана в форме Друде-Зоммер-фельда
е( ю) = 1 - 2 р •
Ю2 + iOY
(19)
Параметр у в (19) имеет смысл частоты столкновений электронов в модели свободного газа, и он определяет мощность диссипативных процессов в металле. В связи с комплексным обобщением проницаемости г(ю) соответственно изменяется и томас-фермиевский параметр (15)
4(ю2 + 1юу)
24 k (ю) = —
(20)
же(ю) 0 п(сор -ю2 - г'ю/)2
Таким образом, в выражении (20), величины к(ю) и г(ю) являются комплекснозначными функциями частоты ю . Везде, далее, будем считать диэлектрическую проницаемость е( ю) и динамическую поляризуемость а(ю) комплексными величинами.
Рассмотрим, теперь, более подробно поведение динамической поляризуемости а( Ю) вблизи резонансной частоты ю = о при наличии диссипации. Вместо резонансного выражения (18) в новых условиях получаем
а(ю)=3я3 ■Ю(Ю-2)[-о2) +,г'Ю'у] ,(21)
2 юр (ю2 -О2)2 + Ю2у2
Rea(o^O).IR-О^^.;Ю[р±Г. ,
2 ;2 (;2-О2)2 + ; у
(22)
з 3 о2 - ю2
Ima(; ^ О) . — R ----------2—p
2 ;2
;уО2
(ю2 -О2)2 +ю у
, (23)
или, оставляя частоту ю « О лишь в резонансном знаменателе
Y( y+ ІО )
4( ю-О )2 + Y2
(24)
Отметим, что кроме пренебрежения 1 вблизи резонансной частоты при выводе (24) использовалось выражение (12). Поэтому формулы (21)-(24), отражающие лоренцеву форму резонансной линии, носят лишь оценочный характер. Точное выражение для поляризуемости сферической наночастицы при произвольных частотах ю^ О, включая и большие отстройки от резонанса, можно записывать только в виде (16).
На рис. 3-8 представлены частотные зависимости действительной Ие а( ю) и мнимой 1т а( ю) частей дипольной поляризуемости а(ю) сферической металлической наночастицы, а также ее модуля | а( ю)| , определенные на основе (16) для различных значений длины I томас-фермиевско-го экранирования. Для кривых на рис. 9 и 10 изменялся, кроме того, радиус Я наношара.
В результате анализа приведенных зависимостей можно сделать следующие выводы. При неизменной величине радиуса наночастицы (И=1 ед. длины) с уменьшением параметра Я /1 от 20 до 4, изменение модуля поляризуемости I а(ю)| остается в пределах от 1 до 14. Для Я/I равном 2,5 (рис. 7) и 0,2 (рис. 8) эти изменения составляют 30 и 60, соответственно. При этом спектральные кривые Яе а(ю), 1ш а(ю) и | а(ю)| испытывают существенные изменения. Наиболее богаты резонансами кривые на рис. 3 -для наибольшего использованного значения параметра Я/1 = 20. При этом глубина модуляции кривых Иеа(ю) , 11ша(ю) и | а(ю)| относительно невелика и составляет, приблизительно 0,4 в диапазоне частот, не выходящем за предел в виде плазменной частоты ю< юр. С уменьшением параметра Я/1 до 10 (рис. 4), число резонансов на кривых Иеа(ю), 1та(ю) и |а(ю)| уменьшается, но при этом растет глубина их модуляции, достигая значения в 2 единицы.
Указанная тенденция сохраняется при уменьшении параметра Я/1 до значения 1: 0,15=6,66(6) (рис. 5) и далее до значения в 4 единицы (рис. 6). Глубина модуляции кривых на рис. 5 возрастает до 5, а на рис. 6 - до 14. При этом
2
ю
р
2
Яе а(а) принимает отрицательные амплитудные значения -1 и -6, соответственно. Для рис. 7 и 8 они достигают значений -12 и -29. Низкочастотный резонансный пик на рис. 7 имеет амплитуду, более чем в 6 раз превышающую высоты других резонансов, число которых уменьшается до 2-3-х. Наконец при Я/1 =0,2 (рис. 8) остается один резонанс, приблизительно отвечающий частоте Ми
0,м =юр/л/3 ~ 0,577ар.
С возрастанием радиуса Я наночастицы до 5 при одновременным увеличении длины I то-мас-фермиевского экранирования до 0,4 (рис. 9), а затем и до 5 (рис. 10), наблюдается частичный возврат к ранее полученным спектральным зависимостям рис. 3 и 4. Качественно кривые рис. 9 подобны кривым рис.3, а рис. 10 - кривым рис. 4, однако масштаб этих графиков существенно различается. При возрастании радиуса Я наносферы поляризуемостьа(а) растет как Я3. Поэтому амплитудные значения Яеа(а), 1та(а) и | а(а)| на рис. 9 и 10 превышают аналогичные величины рис. 3 и 4 на два порядка.
Таким образом, число резонансов поляризуемости определяется лишь значением параметра Я/1, причем чем больше Я/1 по сравнению с 1, тем больше резонансов. Амплитудные значения Иеа( а), 11та( а) и | а( а)|, определяются при неизменном отношении Я/1 еще и радиусом И наночастицы. Увеличение диссипационного параметра у/ ар > 0,01 приводит к вполне ожидаемому уширению резонансных линий и уменьшению амплитуд резонансных пиков.
Межмолекулярный безызлучательный перенос энергии вблизи металлической наночастицы
Для результирующей скорости безызлуча-тельного переноса энергии электронного возбуждения между молекулами в присутствии сферической металлической наночастицы в [13-15] нами получено выражение
и (гОА , ГО , ГА , ^ОЛ , 0) = и Г (гОА , ^ол ) +
+ 9с4хЧПол) 1 Г„, _ ^ . Ла
Не а, 1т о-, \а\
\в(—)Рв —)рА (а) — (25)
ПА ГПА J —
в котором функция —) спектрального отклика наночастицы на внешнее электромагнитное возмущение определяется резонансными плаз-монными факторами и конфигурационными параметрами трехчастичного комплекса «до-нор-наноглобула-акцептор» (рис. 11)
Рисунок 3. Частотные зависимости действительной (1) и мнимой (2) частей поляризуемости а( а) сферической наночастицы, а также ее модуля | а( а)| (3). Значения параметров: а„ = 1, у = 0,01. И=1, 1=0,05. е(а) = 1 -ар /(а2 + шу)
Рисунок 4. Частотные зависимости действительной (1) и мнимой (2) частей поляризуемости а( а) сферической наночастицы, а также ее модуля | а( а)| (3). Значения параметров: а„ = 1, у = 0,01. И=1, 1=0,10. е( а) = 1 -ар/( а2 + гау)
Рисунок 5. Частотные зависимости действительной (1) и мнимой (2) частей поляризуемости а( а) сферической наночастицы, а также ее модуля | а( ю)| (3). Значения параметров: ар = 1, у = 0,01. И=1, 1=0,15.
Ке ач 1т ач |а |
Рисунок 6. Частотные зависимости действительной (1) и мнимой (2) частей поляризуемости а(а) сферической наночастицы, а также ее модуля| а(а)| (3). Значения параметров: ар = 1, у = 0,01. И=1, 1=0,25.
Рисунок 7. Частотные зависимости действительной (1) и мнимой (2) частей поляризуемости а( а) сферической наночастицы, а также ее модуля | а( а)| (3). Значения параметров: юр = 1, у = 0,01.И=1, 1=0,4.
Кеа, 1то', |о-1
со
Рисунок 8. Частотные зависимости действительной (1) и мнимой (2) частей поляризуемости а( а) сферической наночастицы, а также ее модуля | а(ю)| (3). Значения параметров: а = 1, у = 0,01.
И=1, 1=5,0.
Рисунок 9. Частотные зависимости действительной (1) и мнимой (2) частей поляризуемостиа( а) сферической наночастицы большого радиуса И=5, а также ее модуля | а(а)| (3). Значения других параметров: ар = 1, у = 0,01.1=0,4.
Рисунок 10. Частотные зависимости действительной (1) и мнимой (2) частей поляризуемости а( а) сферической наночастицы большого радиуса И=5, а также ее модуля | а( а)| (3). Значения других параметров: юр = 1, у = 0,01.1=5,0.
в(—) = 8
ТпТ.
<|2 а
■ (—)'
ТпТа
х2( &мра)
Х2( пш)
-Х^Яеіа'( а)]
(26)
а( ю) = а'( ю)Я3,
Х(&ва) = йп’&р ■яи'бА со5(ф0 -Фа)-2со5ЬвсоздА —
ориентационный фактор для векторов диполь-ных моментов донора pD и акцептора pA, %(0.ЮА) — аналогичный фактор для векторов дипольных моментов наночастицы PNp и акцептора pA, а функция иЕ(тш,пш) представляет собой скорость переноса в отсутствие наночастицы, совпадающую с выражением теории Ферстера
иг(Гол,Пол) = Г Рв(а)Цл(а)-4 . (27)
пл Гол а
Подинтегральные функции в (25) и (27) ¥в(а), цА(а) — спектры испускания и поглощения молекул донора и акцептора [13-15]. Знак Ие[а'(а)] может быть произвольным [14]. Слагаемое, содержащее Ке[а'( а)] в (26) появляется в результате интерференции амплитуд прямого DA— перехода и опосредованного перехода через наночастицу. В случае
2\а'—)\
гпЛ
ГгТ,
а) <Х(^ N
Л
Re[а'(—)] (28)
ХЧ^пл ) Х(&пл )
интеграл в (25) становится отрицательным и суммарная скорость межмолекулярного переноса энергии в присутствии наночастицы становится меньше скорости переноса без нее. Если неравенство в (28) изменяет знак, общая скорость безызлучательного переноса (25) возрастает по сравнению с невозмущенным случаем (27). Таким образом, поведение общей скорости и(тВА,тв,гА,пш,в) межмолекулярного переноса энергии существенно зависит от величины и знака интеграла двойного спектрального перекрытия
/вОО (а)^л (а) юЮ . (29)
Наличие многих плазмонных резонансов в металлической наночастице с вырожденным электронным газом существенным образом влияет на такое перекрытие. Отметим, что в более ранних моделях [1-15] безызлучательного переноса энергии вблизи наночастицы учет вырождения электронного газа не был произве-
ден, и в качестве дипольнои поляризуемости а(а) металлического шара радиуса Я на частоте ю использовалось выражение
a( a)-
е( а) -еm R3
г( a)+2em ’
(30)
где 8т — диэлектрическая проницаемость прозрачной среды, в которую помещен шар. Очевидно, что в общем случае при расчете функции скорости и(тш,тЕ),гА,Qш,0) для поляризуемости а( а) следует использовать более точную формулу (16).
Перенос энергии возбуждения к металлической наночастице
Средняя диссипируемая в проводящей частице мощность 2 электромагнитного поля определяется через вектор его напряженности E и вектор Р дипольного момента частицы известными соотношениями [23]
(31)
, <//' </г.
о 1 г .
dt dt
Или, используя комплексное представление
(J- ' j- ' RcjWlK j - Jnij /Ч: J.
„ (32)
Переходя к динамической поляризуемости a(w) наночастицы получаем
2 = yIm{<PE*}=^®Ima|E|2. (33)
В целях перехода к выражению для скорости ум диссипации энергии в наночастице разделим (33) на энергию кванта hw
<34>
Здесь, в (34), E=E(r, t) - поле молекулярного диполя донора, осциллирующее на частоте w, поэтому из (34) следует (r=rD)
1 Р2
Yм =Ima(w)~ppDQ (1 + 3cos2 Qd) . (35)
2h e2r°
Угол 0D в (35) образован векторами pD и rD. Важно отметить, что здесь, как и в предыдущих разделах, дипольная поляризуемость а( w) сферической частицы определяется формулой (16) и лишь в предельном случае отсутствия вырождения - формулой (17) или (30). Только в этом - последнем случае, который равносилен выполнению сильного неравенства kR << 1, мы вправе записать
a(a) = R3
Ima(a) = R
[е(ю) - 1][е* (a) + 2 ]
[(г'(a) + 2)2 +г"2(ю)] ‘ Зе\a>)
[ (г '(a) + 2)2 + s" 2(w) ]'
(Зб)
Соответственно, для скорости диссипации (35) в отсутствие вырождения электронного газа можем привести итоговое выражение
3е "(о)
Гм ■ - PD R
2h r6 [(г'(a) + 2)2 + г"2 (a)]
(1 + 3cos20D)
(37)
Отметим, что формула (37) фактически совпадает с формулой Говорова-Ли-Котова [6], которая была получена на основе флуктуаци-онно— диссипационной теоремы. Эти авторы исходили из выражения
їм =
= 2ітГ h J
г (а) Pd Згт
4п є^г з г (а) + 2гт
d r,
в результате чего они и пришли к аналогу (37). Однако, в общем случае вырожденности электронного газа наночастицы вместо (37) необходимо использовать следующее выражение для скорости необратимого монохроматического переноса энергии от возбужденного донора к металлической наночастице
1 Р2
Y м = ^ Ef R3(1 + 3 cos2 0d) x 2h r
X Im
3(1 - k( a)R ctg k( a)R)
_ 2( е( о) -1)(1 - Ь( о)Я ctg Н о)Я) - е( о)(£( о)Я)2
(38)
Монохроматичность энергопереноса подразумевает практическое отсутствие спектральной ширины возбужденного уровня донора. При необходимости такой учет уширения легко мо-
Рисунок 11. Конфигурация трехчастичного комплекса «донор-наноглобула-акцептор»
жет быть произведен путем вычисления интеграла перекрытия
I ¥и (о)1та(о)Ло. (39)
Вблизи резонанса о « О используя (23) получаем
3
Yм(r, 0D, а) = ТГ R 4h
Q2 - а'
ayQ
Р2
- 2 2 2 ?к(1 +3 cos2 0D)
(a2 -Q2 )2 + а2у2 є2гпгб D
(40)
Спектр собственных частот Q для его использования в (40) по-прежнему определяется на основе уравнения (12).
Таким образом, в данной работе исследовано резонансное поведение дипольной поляризуемости сферической металлической наночастицы с учетом вырожденности электронного газа и диссипативных процессов в проводнике. Частотные зависимости действительной и мнимой частей поляризуемости сферической наночастицы демонстрируют наличие большого числа резонансных пиков в области больших значений относительного радиуса R/l частиц, выраженного в длинах томас-фермиевского экранирования. Полученные выражения для скорости безызлуча-тельного переноса энергии между молекулами вблизи наночастицы с вырожденным электронным газом и скорости необратимой передачи энергии к ней, включают в себя модуль динамической поляризуемости a( w) наношара, а также действительную и мнимую части a( w) по отдельности. Мультирезонан-сный характер спектра поляризуемости позволяет проще обеспечить частотное перекры-
2
,, Y w p
тие £ (w)=Ime( w)=-— 1 со спектральны-w ( w2 +y 2)
ми полосами молекул, участвующих в безызлучательном переносе энергии.
Приложение
Для приложений важно рассчитать действительные величины |а(о)|2, Яеа(о) и 1та(о) через комплексную функцию к(а). Не аппелируя к модели Друде-Зоммерфельда вместо (19) можем записать
е(о) = е'(о) + ге"(о). (П1)
Тогда
4 _ 4 -е'( о)+ге"( о)
пе(о)12 п12 е'2(о)+е"2(о) Выделяя действительные и мнимые части 4 -е' ( о )
k2( а) = -
(П2)
x = Rek (а) =
nl2 є"2(a)+є"2(a)
4 є" ( a)
nl2 є"2(a)+є"2(a)
для комплексного k(а) получаем
y = Imk2( a) =
k( a) = yj x + iy = +
r+x r-x
\ 2 + WT
' = •/
x2 + y2 =
nl2 є"2(a) + є"2(a) nl21єI
(П3)
(П4)
(П5)
(Пб)
Заметим, что |k( a)| = r.
Для аналитической функции ctg kR можем
записать
ctg kR =
sin(2RekR)-ish(2 ImkR)
RekR = +
ch(2 ImkR.)-cos(2RekR)
12 і V n
(П7)
I є I -є"
ImkR = +.
Г2 yjI є I +є" V п I є 11
I е\1 ’ | е|1
Таким образом, все величины в (19) могут быть определены через (П1)-(П7).
В модели Друде-Лоренца-Зоммерфельда
є" (а) = Re є( а) = 1-
22 а2 + Y2
є" ( а) = Мє( а) =
Y ар .
a ( a2 + Y2)
30.07.2011
Работа поддержана РФФИ (проект № 10-02-96021-р_урал_а) и Минобрнауки РФ (АВЦП «Развитие научного потенциала ВШ» М.1. Проект № 1.3.06), а также ФЦП «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России на 2007-2011 годы» (ГК № 16.513.11. 3015 и ГК № 16.513.11.3042)
Список литературы:
1. Gersten J., Nitzan A. Spectroscopic properties of molecules interacting with small dielectric particles // J. Chem. Phys. 1981. -V. 75. 1139; doi:10.1063 /1.442161 (14 pages).
2. Gersten J. I., Nitzan A. Accelerated energy transfer between molecules near a solid particle // Chem. Phys. Letters. 1984.
V. 104. -Issue 1. -P. 31-37.
X
а
4
4
1
2
а
р
3. Hua X.M., Gersten J. I., Nitzan A. Theory of energy transfer between molecules near solid state particles //J. Chem. Phys. 1985. V. 83. - Issue 7. -P. 3650 -3659.
4. Klimov V. V., Letokhov V. S. Resonance interaction between two atomic dipoles separated by the surface of a dielectric nanosphere // Phys. Rev. A. 1998. -V. 58. N 4. -P. 3235-3247.
5. Novotny L., Hecht B. Principles of Nano-Optics. Cambridge University Press, New York, 2006. -539 p.
6. Govorov A.O., Lee J., Kotov N.A. Theory of plasmon-enhanced Forster energy transfer in optically excited semiconductor and metal nanoparticles // Phys. Rev. B. 2007. - V. 76. -P. 125308.
7. Durach M., Rusina A., Klimov V. I., Stockman M. I. Nanoplasmonic renormalization and enhancement of Coulomb interactions // New J. Phys. 2008. V.10. 105011.
8. Kucherenko M. G., Chmereva T. M., Kislov D. A. Energy Transfer in Molecular Systems at the Surface of Metal Solids and Nanoparticles // High Energy Chemistry. 2009. Vol. 43. -No. 7. -P. 587-591.
9. Marocico C. A., Knoester J. Intermolecular resonance energy transfer in the presence of a dielectric cylinder // Phys. Rev. A 79. 2009. P. 053816 (1-15)
10. Davis T. J., Gomez D. E., Vernon K. C. Interaction of molecules with localized surface plasmons in metallic nanoparticles // Phys. Rev. B 81. 2010. -P. 045432 (1-11)
11. Chung H. Y., Leung P. T., Tsai D. P. Enhanced Intermolecular Energy Transfer in the Vicinity of a Plasmonic Nanorice // Plasmonics. 2010. DOI 10.1007/s11468-010-9151-x. # Springer Science+Business Media, LLC 2010
12. Pustovit V. N., Shahbazyan T.V. Resonance energy transfer near metal nanostructures mediated by surface plasmons // Phys. Rev. B. 2011. V.83. -P. 085427 (1-5).
13. Кучеренко М.Г., Кислов Д.А., Чмерева Т.М. Повышение качества FRET-SNOM изображений посредством плазмонного резонанса в наноантеннах // Сборник материалов международной научной конференции: «Наука и образование: фундаментальные основы, технологии, инновации». Часть 5. Оренбург: ОГУ. 2010. - 369 с. - С. 351-356.
14. Кучеренко М.Г., Чмерева Т.М. Процессы с участием электронно- возбужденных молекул на поверхностях твердых адсорбентов. Оренбург: Оренбургский государственный университет. Монография. 2010. -346 с.
15. Кучеренко М.Г., Чмерева Т.М., Кислов Д.А. Увеличение скорости межмолекулярного безызлучательного переноса энергии электронного возбуждения вблизи плоской границы твердого тела // Вестник ОГУ. 2011. №1. С. 170-181. http://vestnik.osu.ru/2011_1/30.pdf
16. Давыдов А.С. Теория твердого тела. М.,Наука, 1976, 639 с.
17. Ландау Л.Д., Лифшиц. Е.М. Квантовая механика. Нерелятивистская теория. Т. III. М.: Наука. 1974. - 752 с.
18. Born M., Wolf E. Principles of Optics. Electromagnetic theory of propagation, interference and diffraction of light. Forth edition. Pergamon Press. Oxford-London-Edinburgh-New York-Raris-Frankfurt. 1968.
19. Mie G. // Ann. der Physik. 1908. -V. 25. -P. 377-445.
20. Смирнов М. Б., Крайнов В.П. Многократная ионизация томас- фермиевского кластера сильным электромагнитным полем // ЖЭТФ. 1999. -Т. 115. -Вып. 6. -С. 2014-2019.
21. Крайнов В.П., Смирнов М. Б. Эволюция больших кластеров под действием ультракороткого сверхмощного лазерного импульса // Успехи физ. наук. 2000. -Т. 170. -№ 9. -С. 969-990.
22. Гадомский О.Н., Шалин А.С. Электронные состояния в металлических кластерах // ЖЭТФ. 2007. -Т. 131. - Вып. 5. -С. 5-13.
23. Ландау Л.Д., Лифшиц. Е.М. Электродинамика сплошных сред. Т. VIII. М.: ФИЗМАТЛИТ. 2003. - 656 с.
Сведения об авторе: Кучеренко Михаил Геннадьевич, директор Центра лазерной и информационной биофизики Оренбургского государственного университета, заведующий кафедрой радиофизики и электроники, доктор физ.-мат. наук, профессор 460018, Оренбург, пр-т Победы, 13, тел.: (3532) 364653, 372457, е-mail: rphys@mail.osu.ru
UDC 539.23; 539.216.1 Kucherenko M.G.
Orenburg state university, Centre of Laser and Information Biophysics, Orenburg, e-mail: rphys@mail.osu.ru NANOSPHERE DYNAMIC POLARIZABILITY IN THE CASE OF A DEGENERATE ELECTRONIC GAS AND ITS ROLE IN ENERGY TRANSFER PLASMONIC MECHANISM
The author shows how the degeneracy of electronic gas in spherical metal nanoparticles effect on the resonance behavior of the dipole polarizability in monochromatic electromagnetic field on the optical frequencies with dissipative processes in a conductor (complex dielectric permittivity). Frequency dependence of the imaginary and the real parts of polarizability of a spherical nanoparticles are obtained in the form of graphs for different values of particle radius, length of Thomas-Fermievskiyof shielding and damping constant. The author provides expressions for speed radiationless transfer of energy between the molecules in the vicinity of nanoparticles, and rate of energy to it, which include dynamic polarisability on nanoball with a degenerate electronic gas.
Key words: polarisability spherical nanoparticles, a degenerate electronic gas, Thomas-Fermievskiy, shielding, energy transfer, plasmons.