CC BY
42
СИММЕТРИЯ ИНТЕНСИВНОСТИ И ПОТОКА МОЩНОСТИ СУБВОЛНОВОГО ФОКУСНОГО ПЯТНА
Котляр В.В., Ковалёв А.А., Стафеев С.С.
Институт систем обработки изображений РАН,
Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет)
Аннотация
С помощью разложения светового поля с линейной поляризацией по плоским волнам показано, что эллиптическая форма поперечного сечения интенсивности определяется продольной составляющей электрического вектора, и так как проекция вектора Умова-Пойнтинга на оптическую ось (поток мощности) не зависит от этой продольной составляющей электрического вектора, то сечение потока мощности имеет форму круга. Экспериментально с помощью микроскопа ближнего поля с металлическим кантилевером и малым отверстием показано, что при фокусировке линейно поляризованного гауссова пучка бинарной стеклянной зонной пластинкой с фокусным расстоянием, равным длине волны, формируется фокусное пятно в виде слабого эллипса с размерами по декартовым осям = (0,44 ± 0,02)1 и = (0,52 ± 0,02)1 и глубиной фокуса БОГ = (0,75 ± 0,02)1, 1 -
длина волны. Сравнение экспериментальных данных с результатами моделирования позволило заключить, что микроскоп ближнего поля измеряет поперечную интенсивность (плотность мощности), а не поток мощности и не полную интенсивность. То, что металлический кантилевер с малым отверстием измеряет поперечную интенсивность, следует из теории Бете-Баукампа.
Ключевые слова: субволновый фокус, эллиптичное фокусное пятно, теория Бете-Бау-кампа.
Введение
Субволновая фокусировка света с помощью элементов микрооптики - актуальная задача, так как уменьшение размеров фокусного пятна увеличивает не только разрешение при использовании этого пятна в задачах литографии, микроскопии, оптической памяти, но и плотность мощности излучения, которая важна при манипулировании микрообъектами. В работах последних лет в этом направлении произошло существенное продвижение. Например, для субволновой фокусировки используются планарные плазмонные структуры [1, 2] или плазмонные линзы [3, 4]. Лазерный свет также остро фокусируется вблизи поверхности традиционных оптических элементов: микроаксикона [5], зонной пластинки [6, 7], микролинзы [8], твёрдотельной иммерсионной линзы (81Ь) [9] и обычного высокоапертурного объектива [10 - 13]. Также свет собирается в субволновый фокус на вершине диэлектрических [14, 15] и металлических микроконусов [16]. С помощью элементов микрооптики можно не только фокусировать свет в субволновое пятно, но и изображать микро- и нанообъекты с субволновым разрешением (со сверхразрешением) [17].
Однако в перечисленных работах не рассматриваются особенности регистрации светового излучения в ближнем поле с помощью кантилевера с малым отверстием. До сих пор нет ответов на следующие вопросы. Что регистрирует микроскоп ближнего поля
- плотность энергии (мощности) или поток энергии (мощности)? Почему для линейно-поляризованного света интенсивность (плотность мощности) фокусно-
го пятна имеет вид эллипса, а проекция на оптическую ось вектора Умова-Пойнтинга (поток мощности) в фокусном пятне имеет вид круга? Также в перечисленных работах нет детального сравнения данных эксперимента по измерению субволнового фокуса с результатами строгого моделирования на основе решения уравнений Максвелла.
В данной работе в общем виде с помощью разложения светового поля с линейной поляризацией по плоским волнам показано, что эллиптическая форма поперечного сечения интенсивности определяется продольной составляющей электрического вектора и, так как проекция вектора Умова-Пойнтинга на оптическую ось (поток мощности) не зависит от этой продольной составляющей электрического вектора, сечение потока мощности имеет форму круга. Экспериментально с помощью микроскопа ближнего поля с металлическим кантилевером и малым отверстием показано, что при фокусировке линейно-поляризованного гауссова пучка бинарной стеклянной зонной пластинкой с фокусным расстоянием, равным длине волны, формируется фокусное пятно в виде слабого эллипса с размерами по декартовым осям ГШЫМХ = (0,44 ± 0,02)1 и ГШЫЦ, = (0,52 ± 0,02)1 и глубиной фокуса БОГ = (0,75 ± 0,02)1, 1 - длина волны. Сравнение экспериментальных данных с результатами моделирования ГБТБ-методом позволило однозначно заключить, что микроскоп ближнего поля измеряет поперечную интенсивность (плотность мощности), а не поток мощности и не полную интенсивность. То, что металлический кантилевер с малым отверстием измеряет поперечную интенсивность, следует из теории Бете-Баукампа.
1. Интенсивность и проекция вектора Умова-Пойнтинга для линейно-поляризованного света
Рассмотрим распространение света вдоль оптической оси г от одной поперечной плоскости Р1 (будем называть её начальной плоскостью) до другой поперечной плоскости Р2, параллельной начальной и находящейся от неё на расстоянии г. Введём в этих плоскостях декартовы координаты (х, у) и (и, V) и полярные координаты (г, ф) и (р, 8). Пусть в начальной плоскости сформировано линейно-поляризованное электромагнитное поле с радиальной симметрией:
Ех (Г, ф,0)° Ех (Г ) ,
■Еу (г, ф,0) ° 0, (1)
Ег (^ ф,0)° 0
В плоскости Р2 будем измерять интенсивность
I 12 I 12 I |2 | |2
I = |е| = |Ех | +| Еу | +| Е2| (2)
и поток мощности (составляющую вектора Умова-Пойнтинга, параллельную оптической оси г)
= ^{(ЕхН*)г} = ^Яе^Я* -ЕуЯ*} . (3)
Согласно дифракционным интегралам Рэлея-Зоммерфельда, составляющая Еу и в плоскости Р2 также будет равна нулю. Поэтому (2) и (3) можно переписать в виде:
(4)
і і2 і |2 I = \е + E ,
Sz = ^Re^H*} . (5)
Используя уравнение Максвелла для монохроматического света с частотой w
rot E = -i wm0 m H, (6)
где m - магнитная проницаемость среды, m0 - магнитная постоянная, получим из (5):
S = Re
I-
E
f dE* _dE^ I
dz du
(7)
Разложим составляющую Ех в угловой спектр плоских волн. Введём для этого декартовы (а, Р) и полярные (£, ф) координаты в спектральной плоскости. Тогда
Ex (u, v, z) = JJ A (a, p) exp{ik [au + pv +
к2
+zj 1 -a2-p2 }dadp.
Из третьего уравнения Максвелла dE dEv dE
+
+
du dv dz
=0
(8)
(9)
с учётом того, что Еу = 0, получим выражение для составляющей Ег:
Ez (u,v, z)=-JJ-r
к2 V1
a
-a2-p:
=A (a, p)x
x exp {ik
+C (u, v).
au + pv + z-J 1 -a2 -p2 }dadp
(10)
Константа C (u, v), возникшая при интегрировании по z, означает постоянное поле бесконечной протяжённости вдоль оси z. Из физических соображений приравняем её нулю. Преобразуем также выражение в скобках в (7):
E Ж=_ik ff 1 -b2
5z Эи i^Vl-a2 _b (11)
x exp {-ik
^A* (a, p)>
{-ік аи + Ру + 1 -а2 -р2 }dadp.
Т ак как в начальной плоскости составляющая Ех радиально симметрична, её угловой спектр также радиально симметричный, т.е. А (£, ф) = А (£). Учитывая это, перепишем (8), (10) и (11) в полярных координатах. При этом все интегралы по ф выразятся через функции Бесселя:
Z2dZ
Ex (p,Є,z) = 2pJA(Z)exp(/^д/ЇЧ2) Jo (kpZ)ZdZ
0
Ez (p, Є, z) = -2p cos eJ A(C)exp(ikW1-C2) J1 (kpc^-r—2.
0 V1 -Z
¥ Ґ C- "> \
2pikJ A* (Z) exp (-ikz^/1 -Z2)
dE* - dE*
z u 0
Подставим (12)-(14) в (4) и (7):
1
2
Z2
Jo (kpZ)-^ J2 (kpz)cos(2Є)
ZdZ
I = 4p2
+ 4p2 cos2 Є
= - 2ri Re
JA(Z)exp(kz^V ) Jo (kpZ)ZdZ
JA(Z)exp(ikzV1 -Z2)Jo (kpZ)ZdZ^x
JA(Z)exp(ikzV1-Z2)J1 (kpZ) Z~dZ7
o V1 -Z
wmom
x IJ A* (Z) exp (-ikzVb-Z7) f1 -£. I Jo (k pZ)--Z2 J2 (kpZ) cos (2Є)
ZdZ
(12)
(13)
(14)
(15)
(16)
лА-с2
+
-
2
2
Из (15) и (16) видно, что и интенсивность, и поток мощности симметричны относительно оси х. Т.е. для любых пар точек с полярными координатами (р, 8) и (р, -8) они равны. Из (15) также видно, что для фиксированного р интенсивность в точках (р, 0) и (р, р) максимальна, а в точках (р, р/2) и (р, 3р/2) минимальна. Этим и объясняется формирование фокального пятна в форме эллипса, вытянутого вдоль оси х (плоскость поляризации хг). Заметим, что нарушение радиальной симметрии в выражении для интенсивности определяется вторым слагаемым, а в выражении для потока мощности - вторым множителем. В случае малой числовой апертуры множитель С 2/(1 - С 2)12 во втором интеграле в (15) близок к нулю и не позволяет второму слагаемому давать существенный вклад в интенсивность. Поэтому фокальное пятно имеет круглую форму. Если же фокусировка острая, то в угловом спектре существенную часть занимают волны с С, близким к единице. В этом случае вклад второго слагаемого возрастает, может превысить вклад первого слагаемого в несколько раз, появляется асимметрия, и пятно принимает эллиптическую форму или даже форму «кости» или «гантели».
Для потока мощности зависимость иная. Множитель С /(1 - С 2)12 одинаково влияет и на симметричную, и на асимметричную часть выражения в квадратных скобках в (16). В случае малой числовой апертуры (С << 1) поток мощности определяется радиально симметричным слагаемым (1 -С"2/2) Л (к р С). Поэтому фокальное пятно, измеренное не по интенсивности, а по потоку мощности, имеет круглую форму. В случае большой числовой апертуры, когда С близко к единице, вклад обоих слагаемых в квадратных скобках примерно одинаков, но вклад второго слагаемого не может существенно превысить вклад первого, и, кроме того, вблизи фокуса (р = 0) функция Бесселя нулевого порядка оказывает большее влияние по сравнению с функцией Бесселя второго порядка. Поэтому пятно получается более круглой формы, чем в случае измерения интенсивности.
Для примера рассмотрим бесселев пучок. У него угловой спектр имеет кольцевую форму
А (С) = 8(С-а), (17)
где 8 (х) - дельта-функция Дирака, а - параметр бесселева пучка. В этом случае интегралы в (15) - (16) пропадают и интенсивность с потоком мощности легко вычисляются:
I = [ 2ра/0 (к ар) ]2 + 2р2 к
/0 (кар)
2ра2
л/ъ-
а2
/1 (к ар) соб 0
( а2 ^
1 -—
2
а
/0 (к ар)—— /2 (к ар) соб (20)
а
(18)
(19)
а
2
Для моделирования случая малой числовой апертуры выберем а = 0,3 (рис. 1). Для моделирования острой фокусировки выберем а = 0,8 (рис. 2) и а = 0,9 (рис. 3). Интенсивность и поток мощности для разных значений а показаны на рис. 1 - 3: горизонтальная ось - ось х, вертикальная ось - ось у. Другие параметры моделирования: длина волны 1 = 532 нм, расстояние г = 101, область моделирования -51 £ х £ +51, -51 £ у £ +51.
Рис. 1. Интенсивность (а) и поток мощности ( для бесселева пучка с а= 0,3
Рис. 1 - 3 подтверждают ранее сделанные предположения. При увеличении числовой апертуры фокальное пятно, измеренное по интенсивности, сначала имело круглую форму, затем эллиптическую, затем форму «кости». При этом фокальное пятно,
измеренное по потоку мощности, всё время оставалось почти круглым (на рис. 3б фокусное пятно имеет форму слабого эллипса, вытянутого вдоль оси у, это нетрудно показать с помощью уравнения (19)).
Рис. 2. Интенсивность (а) и поток мощности (б) для бесселева пучка с а= 0,8
Рис. 3. Интенсивность (а) и поток мощности для бесселева пучка с а= 0,9
Из рис. 2 и 3 можно видеть, что размеры фокусных пятен для интенсивности и потока мощности вдоль вертикальной оси близки.
2. Острая фокусировка линейно-поляризованного
света с помощью субволнового бинарного аксикона
В предыдущем параграфе с помощью соотношений (18), (19) моделировалось распространение линейно-поляризованного светового пучка с начальной амплитудой в виде функции Бесселя. В этом параграфе, наоборот, рассматривается формирование бесселева пучка с линейной поляризацией с помощью бинарного микроаксикона РБТБ-методом. Мы покажем, что дифракционные картины в обоих случаях имеют одинаковую структуру.
Моделировалось прохождение линейно-поляризованного гауссова пучка с длиной волны 1 = 532 нм и радиусом перетяжки а = 71 через бинарный микроак-сикон с радиусом Я = 8 мкм, показателем преломления материала п = 1,52 и высотой рельефа И = 532 нм. Моделирование проводилось с помощью программы РиІ^АУБ [18], реализующей ББТБ-метод, шаг сетки по пространству - 1/50, по времени - 1/100с, где с
- скорость света в вакууме. Интересно заметить, что при расчёте не интенсивности (плотность световой мощности), а потока световой мощности (модуль проекции на оптическую ось вектора Умова-Пойнтинга) вместо эллиптического фокусного пятна возникает круглое пятно с диаметром немного большим, чем меньший диаметр эллиптического фокусного пятна интенсивности.
В этом номере журнала в статье этих же авторов «Моделирование фокусировки линейно-поляризованного света с помощью субволнового бинарного аксикона» показана двумерная картина интенсивности (в полутонах серого) в непосредственной близости от бинарного аксикона с периодом Т = 1 = 532 нм, а также картина потока световой мощности в той же плоскости.
Из этих рисунков видно, что поток мощности формирует круглый фокус с диаметром по полуспаду, равным FWHM=0,361, в то время как интенсивность формирует эллиптический фокус с диаметрами по декартовым осям FWHMx=0,751 и FWHMy=0,301. Это объясняется тем, что продольная составляющая вектора электрического поля Ег не даёт вклада в продольную составляющую вектора Умова-Пойнтинга. Составляющая Ег возникает в плоскости поляризации Х2 при схождении лучей в фокус. В перпендикулярной плоскости У2 при схождении лучей в фокус составляющая Ег не возникает. Картины фокусов на этих рисунках качественно совпадают с картинами на рис. 3.
3. Острая фокусировка линейно-поляризованного
света с помощью зонной пластинки
3.1. Моделирование
Рассматривалась фокусировка линейно-поляризованного гауссова пучка с длиной волны 1 = 532 нм и
радиусом w = 71 бинарной зонной пластинкой (ЗП) с фокусным расстоянием, равным одной длине волны фокусируемого света f= 1, радиусом 10,64 мкм (201), показателем преломления материала 1,52. На рис. 4 показан шаблон такой ЗП в поле расчёта. Числовая апертура ЗП равна NA = 0,997. Рассчитанная эффективность фокусировки - 42%.
Contour Map of Transverse Index Profile at Z=1
Y, \ш
1,52
Рис. 4. Шаблон ЗП с фокусным расстоянием, равным длине волны / = 1 (а), и распределение интенсивности в плоскости фокуса (б). Ось у - ось поляризации
Радиусы ЗП (рис. 4а) рассчитывались по известной формуле гт = (ш1/+ т212/4)1/2, где f= 1 = 532 нм -фокусное расстояние, т - номер радиуса зоны. Моделирование проводилось БОЯ-РБТБ методом, при этом проводилось сравнение результатов, полученных при разной дискретности разбиения вычисляемой области по пространству. На рис. 4б показано в псевдоцветах распределение интенсивности в фокусе.
На рис. 5 показаны сечения по декартовым осям: по оси х (ф = 0) и по оси у (ф = п/2) интенсивности (а) и потока мощности (модуль проекции на ось г вектора Умова-Пойнтинга) (б). Из рис. 4б и рис. 5а видно, что фокальное пятно интенсивности эллиптическое. А из рис. 5б также видно, что фокальное пятно потока мощности - круглое. В табл. 1 первые две строчки значений показывают размер фокусного пятна по по-луспаду интенсивности (первая строка) и потока мощности (вторая строка).
б) -2 -1 0 1 г,Х
Рис. 5. Распределения (а) интенсивности и (б) модуля проекции на ось 2 вектора Умова-Пойнтинга в фокусе
при дискретности разбиения Аг = 1/50. Сечение по оси х (ф= 0) и по оси у (ф = п/2) Таблица 1. Значения диаметров фокусного пятна
FWHM* (Ф = 0), 1 FWHMj, (ф = п/2), 1 DOF, 1
Интенсивность 0,42 і 0,01 0,84 і 0,01 0,8б і 0,01
Модуль проекции на ось 2 вектора Умова-Пойнтинга 0,45 і 0,01 0,45 і 0,01 -
Эксперимент на ШОМ 0,44 і 0,02 0,52 і 0,02 0,75 і 0,02
На рис. 6 показаны распределения интенсивности вдоль оптической оси для ЗП, рассчитанные с помощью метода BOR-FDTD в пакете Matlab (кривая 1) и методом FDTD, реализованного в программе FullWAVE (кривая 2).
Значения глубины фокуса для ЗП показаны в последнем столбце табл. 1.
3.2. Эксперимент
Зонная пластинка высокого качества была изготовлена по технологии литографии из резиста ZEP (показатель преломления резиста n = 1,52). На рис. 7 показаны вид сбоку (а) и сверху (б) ЗП в сканирующем микроскопе на атомных силах: глубина рельефа - 510 нм, диаметр - 14 мкм, крайняя зона -
0,51 = 266 нм. ЗП имеет 12 колец и центральный диск.
Экспериментально с помощью сканирующего ближнепольного оптического микроскопа (СБОМ) Ntegra Spectra (рис. 8) изучалось прохождение ли-
нейно поляризованного гауссова пучка с длиной волны 1 = 532 нм через зонную пластинку с фокусным расстоянием /= 1.
\Е2\, а.и.
Рис. 6. Сравнение распределения интенсивности вдоль оптической оси для ЗП при моделировании в БОК-РОТО
Matlab (кривая 1) и FDTD FullWAVE (кривая 2)
j__________і_________і_________і_________і__________і_________і_________і_________L
б) 2 4 6 8 10 12 14 16 х, мкм
Рис. 7. Вид ЗП сбоку (а) и сверху (б), полученные на сканирующем микроскопе на атомных силах
На рис. 9а показано экспериментальное распределение интенсивности вдоль оптической оси для ЗП (кривая, ось слева), а также значения диаметра фокусного пятна (квадраты, ось справа). На рис. 9б показан пример распределения интенсивности в фокусе (изображение, полученное непосредственно с микроскопа).
Усреднённые значения диаметра эллиптического фокусного пятна ЗП приведены в третьей строке значений в табл. 1. Из рис. 9а видно, что максимальное значение интенсивности на оси смещено от плоскости геометрического фокуса /= 1 = 532 нм ближе к ЗП (2 = 400 нм). Заметим, что на этом расстоянии
(2 = 400 нм) меньшие диаметры фокусного пятна (расчётный и экспериментальный) совпадают ^ИМ = 0,421).
Рис. S. СБОМ, используемый в эксперименте
Intensity, а.и. FWHM, X
0,46
а) 0,2 0,4 0,6 0,8 z, мкм
0,40
б) 0 0,4 0,8 1,2 х, мкм
Рис. 9. Экспериментальное распределение интенсивности вдоль оптической оси (а) для ЗП пластинки на рис. 7 (кривая, ось слева) и значения меньших диаметров фокусных пятен (квадраты, ось справа); поперечное сечение фокусного пятна на фокусном расстоянии /= 1 = 532 нм (вертикальная ось совпадает с плоскостью поляризации) (б)
Для сравнения на рис. 10 приведены сечения вдоль оси х, перпендикулярной плоскости поляризации, интенсивности (кривая 1), потока мощности (кривая 3) и экспериментальная кривая, измеренная с помощью К80М (кривая 2). Кривые почти совпадают (см. второй столбец табл. 1): их отличия
меньше ошибки измерения (±0,021). Отличия заметны только в боковых лепестках. Но и по ним нельзя однозначно сказать, что измеряется в эксперименте - интенсивность или поток мощности, так как боковые лепестки расчётного потока мощности немного меньше, чем у экспериментальной кривой, а боковые лепестки кривой расчётной интенсивности немного больше, чем у экспериментальной кривой.
-2 -1 0 1 х, X
Рис. 10. Сравнение экспериментального и расчётного распределения в фокусе вдоль оси х: расчётное распределение интенсивности (кривая 1), экспериментальное распределение интенсивности (кривая 2) и расчётное распределение модуля проекции на ось 2 вектора Умова-Пойнтинга (кривая 3)
На рис. 11 показаны сечения кривых в фокусе вдоль оси у (параллельной плоскости поляризации): расчётное распределение модуля проекции на ось 2 вектора Умова-Пойнтинга (кривая 3), экспериментальное распределение интенсивности (кривая 2) и расчётное распределение интенсивности (кривая 1), вычисленной как сумма всех составляющих (а) и как сумма только поперечных составляющих (б).
Из рис. 11а видно, что в эксперименте не измеряется продольная составляющая интенсивности (см. третий столбец Таблицы 1), так как кривая полной интенсивности (Б’^НМ = 0,841) шире экспериментальной кривой (FWHM = 0,521) на величину, большую ошибки измерения (±0,021). Но, в свою очередь, экспериментальная кривая шире расчётной кривой потока мощности (FWHM = 0,451) также на величину, большую ошибки эксперимента. Но что же измеряется в эксперименте? На рис. 11 б приведены сравнения экспериментальной кривой (кривая 2) с поперечной интенсивностью |Ех|2 + |Еу|2 (кривая 1). Их ширины совпадают (FWHM = 0,521). Поэтому из рис. 11 можно однозначно заключить, что микроскоп ближнего поля К80М с кантилевером в виде металлической пирамиды с отверстием 100 нм (рис. 12) измеряет поперечную интенсивность |Ех|2 + |Еу|2 (плотность мощности), а не поток мощности и не полную интенсивность |Ех|2 + |Еу|2 + |Е2|2. Поэтому эллипс фокусного пятна на рис. 9б менее выражен, чем расчётный эллипс на рис. 4б.
\Е2\, а.и.
\Е\ а.и.
Рис. 11. Сравнение экспериментального и расчётного распределения в фокусе вдоль оси у (параллельной плоскости поляризации): экспериментальное распределение интенсивности (кривая 2), расчётное распределение модуля проекции на ось г вектора Умова— Пойнтинга (кривая 3) и расчётное распределение интенсивности (кривая 1), вычисленной как сумма всех составляющих (а) и как сумма только поперечных составляющих (б)
Рис. 12. Электронное изображение верхней части металлического кантилевера сканирующего ближнеполного микроскопа в виде пирамиды с отверстием размером 100 нм
Проникновение электромагнитного поля в малое отверстие в металлическом экране описывается в рамках теории Бете-Баукампа [19-21]. Согласно этой теории плоская линейно-поляризованная волна, падающая под углом на металлический экран с малым круглым отверстием диаметром а << 1, индуцирует
электрическии диполь, ориентированный перпендикулярно отверстию, и магнитный диполь, расположенный в плоскости отверстия. И поэтому поле в дальней зоне для малого отверстия, освещённого наклонной плоской волной Е, задаётся излучением электрического Р и магнитного М диполей с моментами
4 3 / \
Р= - 3е0а (ЕП )Пz,
38 (20)
М = - 3 а3 [п х[Е х п ]],
где п2 - единичный вектор вдоль оптической оси (перпендикулярной плоскости отверстия). Из (20) видно, что электрический диполь формирует только продольная составляющая электрического поля Е. Но диполь, ориентированный вдоль оптической оси, излучает в поперечном направлении и не излучает вдоль самой оптической оси. Магнитный диполь в (20) формируется, наоборот, только поперечными составляющими электрического поля, так как внутреннее векторное произведение в правой формуле в (20) равно нулю для продольной составляющей электрического поля. То есть продольная составляющая электрического поля не будет регистрироваться фотоприёмником, расположенным на некотором расстоянии от малого отверстия в металле вдоль оптической оси.
Заключение
В данной работе получены следующие результаты.
1. В общем случае с помощью разложения векторного светового поля по плоским волнам показано, что для линейно поляризованного начального светового поля с высокой числовой апертурой (полуугол расходимости светового пучка близок к 90 градусам) поперечное распределение интенсивности (плотности мощности) имеет вид эллипса или «гантели», вытянутых вдоль линии, параллельной плоскости поляризации начального светового поля, а поперечное распределение потока мощности (проекции вектора Умова-Пойнтинга на оптическую ось) имеет вид круга или эллипса, но вытянутого вдоль линии, перпендикулярной плоскости поляризации начального поля.
2. С помощью моделирования ББТБ методом показано, что при субволновой фокусировке линейно поляризованного света бинарным аксиконом с периодом, равным длине волны, в непосредственной близости от поверхности аксикона формируется фокусное пятно интенсивности в виде «гантели», при этом фокусное пятно потока мощности (модуля проекции на оптическую ось вектора Умова-Пойн-тинга) в том же месте имеет вид круга.
3. С помощью моделирования ББТБ методом показано, что при фокусировке линейно поляризованного света бинарной зонной пластинкой из стекла с фокусным расстоянием, равным длине волны, фокусное пятно интенсивности имеет вид эллипса с меньшим диаметром по полуспаду FWHM=(0,42 ± 0,01)1, а фокусное пятно потока мощности имеет форму круга с диаметром по полуспаду FWHM = (0,45 ± 0,01)1.
4. Экспериментально с помощью микроскопа ближнего поля с металлическим кантилевером с отверстием 100 нм показано, что при фокусировке линейно поляризованного гауссова пучка бинарной зонной пластинкой с фокусным расстоянием, равным длине волны, формируется фокусное пятно в виде слабого эллипса с размерами по декартовым осям FWHM* = (0,44 ± 0,02)1 и FWHMy = (0,52 ± 0,02)1 и глубиной фокуса DOF = (0,75 ± 0,02)1. Сравнение экспериментальных данных с результатами моделирования позволило заключить, что микроскоп ближнего поля измеряет поперечную интенсивность (плотность мощности), а не поток мощности и не полную интенсивность. То, что металлический кантилевер с малым отверстием измеряет поперечную интенсивность + |Ey|2, следует из теории Бе-
те-Баукампа.
Благодарности
Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» (госконтракт № 14.740.11.0016), грантов Президента РФ поддержки ведущих научных школ (НШ-4128.2012.9) и молодого кандидата наук (МК-3912.2012.2) и гранта РФФИ (12-07-00269).
Литература
1. Jia, B. Near-field visualization of focal depth modulation by step corrugated plasmonic slits / B. Jia, H. Shi, J. Li, Y. Fu, C. Du, M. Gu // Appl. Phys. Lett. - 2009. -Vol. 94. - P. 151912.
2. Chen, K.R. Beyond-limit light focusing in the intermediate zone / K.R. Chen, W.H. Chu, H.C. Fang, C.P. Liu, C.H. Huang, H.C. Chui, C.H. Chuang, Y.L. Lo, C.Y. Lin, H.H. Hwung, A.Y.-G. Fuh // Opt. Lett. - 2011. - Vol. 36, N 23. - P. 4497-4499.
3. Yu, Y. Effect of lens size on the focusing performance of plasmonic lenses and suggestions for the design / Y. Yu,
H. Zappe // Opt. Express. - 2011. - Vol. 19, N 10. -P. 9434-9444.
4. Liu, Y. Far-field superfocusing with an optical fiber based surface plasmonic lens made of nanoscale concentric annular slits / Y. Liu, H. Xu, F. Stief, N. Zhitenev, M. Yu // Opt. Express. - 2011. - Vol. 19, N 21. - P. 20233-20243.
5. Kotlyar, V.V. Tight focusing with a binary microaxicon / V.V. Kotlyar, S.S. Stafeev, L. O’Faolain, V.A. Soifer // Opt. Lett. - 2011. - Vol. 36, N 16. - P. 3100-3102.
6. Mote, R.G. Experimental demonstration of near-field focusing of a phase micro-Fresnel zone plate (FZP) under linear polarized illumination / R.G. Mote, S.F. Yu, A. Kumar, W. Zhou, X.F. Li // Appl. Phys. B. - 2011. -Vol. 102, N 1. - P. 95-100.
7. Стафеев, C.C Субволновая фокусировка с помощью зонной пластинки Френеля с фокусным расстоянием 532 нм / C.C Стафеев, Л. О’Фаолейн, М.И. Шанина,
В.В. Котляр, В.А. Сойфер // Компьютерная оптика. -2011. - Т. 35, № 4. - с.460-461.
8. Ye, J.-S. Long-focal-depth cylindrical microlens with flat axial intensity distributions / J.-S. Ye, G.-A. Mei, X.-H. Zheng, Y. Zhang // J. Mod. Opt. - 2012. - Vol. 59, N 1. - P. 90-94.
9. Huang, K. Realization of a subwavelength focused spot without a longitudinal field component in a solid immer-
sion lens-based system / K. Huang, Y. Li // Opt. Lett. -2011. - Vol. 36, N 18. - P. 3536-3538.
10. Yuan, G.H. Nondiffracting transversally polarized beam / G.H. Yuan, S.B. Wei, X.-C. Yuan // Opt. Lett. - 2011. -Vol. 36, N 17. - P. 3479-3481.
11. Li, X. Superresolution-focal-volume induced 3.0 Tbytes/disk capacity by focusing a radially polarized beam / X. Li, Y. Cao, M. Gu // Opt. Lett. - 2011. - Vol. 36, N 13.
- P. 2510-2512.
12. Lin, J. Achievement of longitudinally polarized focusing with long focal depth by amplitude modulation / J. Lin, K. Yin, Y. Li, J. Tan // Opt. Lett. - 2011. - Vol. 36, N 7. -P. 1185-1187.
13. Lin, H. Generation of an axially super-resolved quasi-spherical focal spot using an amplitude-modulated radially polarized beam / H. Lin, B. Jia, M. Gu // Opt. Lett. - 2011.
- Vol. 36, N 13. - P. 2471-2473.
14. Kotlyar, V.V. Modeling the sharp focus of a radially polarized laser mode using a conical and a binary microaxi-con / V.V. Kotlyar, S.S. Stafeev // J. Opt. Soc. Am. B. -
2010. - Vol. 27, N 10. - P. 1991-1997.
15. Martin, J. Plain Intense Bessel-like beams arising from pyramid-shaped microtips / J. Martin, J. Proust, D. Gerard, J.-L. Bijeon, J. Plain // Opt. Lett. - 2012. - Vol. 37, N 7. -P. 1274-1276.
16. De Angelis, F. Breaking the diffusion limit with superhydrophobic delivery of molecules to plasmonic nanofocus-ing SERS structures / F. De Angelis, F. Gentile, F. Me-carini, G. Das, M. Moretti, P. Candeloro, M.L. Coluccio,
G. Cojoc, A. Accardo, C. Liberale, R.P. Zaccaria, G. Peroz-ziello, L. Tirinato, A. Toma, G. Cuda, R. Cingolani, E. Di Fabrizio // Nature Photonics. - 2011. - Vol. 5. - P. 682-687.
17. Rogers, E.T.F. A super-oscillatory lens optical micro-
scope for subwavelength imaging / E.T.F. Rogers, J. Lind-berg, T. Roy, S. Savo, J.E. Chad, M.R. Dennis, N.I. Zhe-ludev // Nature Materials Lett. - 2012. -
doi:10.1038/nmat3280.
18. http://www.rsoftdesign.com/products.php?sub=Componen t+Design&itm=F ullWAVE.
19. Новотный, Л. Основы нанооптики / Л. Новотный, Б. Хехт. - М.: Физматлит, 2011. - 484 с.
20. Michalski, K.A. Complex image method analysis of a plane wave-exited subwavelength circular aperture in a planar screen / K.A. Michalski // Prog. Electromag. Res. B. -
2011. - Vol. 27. - P. 253-272.
21. Wu, J.H. Modeling of near-field optical diffraction from a subwavelength aperture in a thin conducting film / J.H. Wu // Opt. Lett. - 2011. - Vol. 36, N 17. - P. 3440-3442.
References
1. Jia, B. Near-field visualization of focal depth modulation by step corrugated plasmonic slits / B. Jia, H. Shi, J. Li, Y. Fu, C. Du, M. Gu // Appl. Phys. Lett. - 2009. -Vol. 94. - P. 151912.
2. Chen, K.R. Beyond-limit light focusing in the intermediate zone / K.R. Chen, W.H. Chu, H.C. Fang, C.P. Liu,
C.H. Huang, H.C. Chui, C.H. Chuang, Y.L. Lo, C.Y. Lin,
H.H. Hwung, A.Y.-G. Fuh // Opt. Lett. - 2011. - Vol. 36, N 23. - P. 4497-4499.
3. Yu, Y. Effect of lens size on the focusing performance of plasmonic lenses and suggestions for the design / Y. Yu, H. Zappe // Opt. Express. - 2011. - Vol. 19, N 10. -P. 9434-9444.
4. Liu, Y. Far-field superfocusing with an optical fiber based surface plasmonic lens made of nanoscale concentric an-
nular slits / Y. Liu, H. Xu, F. Stief, N. Zhitenev, M. Yu // Opt. Express. - 2011. - Vol. 19, N 21. - P. 20233-20243.
5. Kotlyar, V.V. Tight focusing with a binary microaxicon / V.V. Kotlyar, S.S. Stafeev, L. O’Faolain, V.A. Soifer // Opt. Lett. - 2011. - Vol. 36, N 16. - P. 3100-3102.
6. Mote, R.G. Experimental demonstration of near-field focusing of a phase micro-Fresnel zone plate (FZP) under linear polarized illumination / R.G. Mote, S.F. Yu, A. Kumar, W. Zhou, X.F. Li // Appl. Phys. B. - 2011. -Vol. 102, N 1. - P. 95-100.
7. Stafeev, S.S. Subwavelength focusing using Fresnel zone plate with focal length 532 nm / S.S. Stafeev, L. O’Faolain, M.I. Shanina, V.V. Kotlyar, V.A. Soifer // Computer optics. - 2011. - Vol. 35(4). - P. 460-461. - (In Russian).
8. Ye, J.-S. Long-focal-depth cylindrical microlens with flat axial intensity distributions / J.-S. Ye, G.-A. Mei, X.-H. Zheng, Y. Zhang // J. Mod. Opt. - 2012. - Vol. 59, N 1. - P. 90-94.
9. Huang, K. Realization of a subwavelength focused spot without a longitudinal field component in a solid immersion lens-based system / K. Huang, Y. Li // Opt. Lett. -2011. - Vol. 36, N 18. - P. 3536-3538.
10. Yuan, G.H. Nondiffracting transversally polarized beam / G.H. Yuan, S.B. Wei, X.-C. Yuan // Opt. Lett. - 2011. -Vol. 36, N 17. - P. 3479-3481.
11. Li, X. Superresolution-focal-volume induced 3.0 Tbytes/disk capacity by focusing a radially polarized beam / X. Li, Y. Cao, M. Gu // Opt. Lett. - 2011. - Vol. 36, N 13.
- P. 2510-2512.
12. Lin, J. Achievement of longitudinally polarized focusing with long focal depth by amplitude modulation / J. Lin, K. Yin, Y. Li, J. Tan // Opt. Lett. - 2011. - Vol. 36, N 7. -P. 1185-1187.
13. Lin, H. Generation of an axially super-resolved quasi-spherical focal spot using an amplitude-modulated radially
polarized beam / H. Lin, B. Jia, M. Gu // Opt. Lett. - 2011.
- Vol. 36, N 13. - P. 2471-2473.
14. Kotlyar, V.V. Modeling the sharp focus of a radially polarized laser mode using a conical and a binary microaxi-con / V.V. Kotlyar, S.S. Stafeev // J. Opt. Soc. Am. B. -
2010. - Vol. 27, N 10. - P. 1991-1997.
15. Martin, J. Plain Intense Bessel-like beams arising from pyramid-shaped microtips / J. Martin, J. Proust, D. Gerard, J.-L. Bijeon, J. Plain // Opt. Lett. - 2012. - Vol. 37, N 7. -P. 1274-1276.
16. De Angelis, F. Breaking the diffusion limit with superhydrophobic delivery of molecules to plasmonic nanofocus-ing SERS structures / F. De Angelis, F. Gentile, F. Me-carini, G. Das, M. Moretti, P. Candeloro, M.L. Coluccio, G. Cojoc, A. Accardo, C. Liberale, R.P. Zaccaria, G. Peroz-ziello, L. Tirinato, A. Toma, G. Cuda, R. Cingolani, E. Di Fabrizio // Nature Photonics. - 2011. - Vol. 5. - P. 682-687.
17. Rogers, E.T.F. A super-oscillatory lens optical micro-
scope for subwavelength imaging / E.T.F. Rogers, J. Lind-berg, T. Roy, S. Savo, J.E. Chad, M.R. Dennis, N.I. Zhe-ludev // Nature Materials Lett. - 2012. -
doi:10.1038/nmat3280.
18. http://www.rsoftdesign.com/pro-ducts.php?sub=Component+Design&itm=FullWAVE.
19. Novotny, L. Principles of nano-optics / L. Novotny, B. Hecht.
- Cambridge University Press, 2006. - 539 p.
20. Michalski, K.A. Complex image method analysis of a plane wave-exited subwavelength circular aperture in a planar screen / K.A. Michalski // Prog. Electromag. Res. B. -
2011. - Vol. 27. - P. 253-272.
21. Wu, J.H. Modeling of near-field optical diffraction from a subwavelength aperture in a thin conducting film / J.H. Wu // Opt. Lett. - 2011. - Vol. 36, N 17. - P. 3440-3442.
INTENSITY AND POWER FLOW SYMMETRY OF SUBWAVELENGTH FOCAL SPOT
V. V. Kotlyar, A.A. Kovalev, S.S. Stafeev Image Processing Systems Institute of the RAS,
Samara State Aerospace University named after S.P. Korolyov (National Research University)
Abstract
Using plane wave expansion method it was shown that elliptical shape of focal spot is determined by the longitudinal component of electric field. However the cross section of power flow has the shape of circle because it is independent of the longitudinal component of electric field. Using near-field scanning optical microscope with metallic cantilever it was experimentally shown that binary Fresnel zone plate with focal length equals to wavelength illuminated by linearly polarized Gaussian beam forms an elliptical focal spot with diameters FWHMX = (0,44 ± 0,02)1 and FWHMy = (0,52 ± 0,02)1 and depth of focus DOF = (0,75 ± 0,02)1. The comparison between experiment and simulation confirms that near-field scanning optical microscope measures the transverse intensity rather than the power flow and the total intensity. Bethe-Bouwkamp law explains that metallic cantilever with small hole measures only the transverse intensity.
Key words: subwavelength focusing, elliptical focal spot, Bethe-Bouwkamp law.
Сведения об авторе Котляр Виктор Викторович — см. стр. 164 этого номера.
Сведения об авторах Ковалёв Алексей Андреевич, Стафеев Сергей Сергеевич - см. стр. 188-189 этого номера.
Поступила в редакцию 11 марта 2012 г.
