Градиентные элементы микрооптики для достижения сверхразрешения

Котляр Виктор Викторович; Ковалёв Алексей Андреевич; Налимов Антон Геннадьевич

ГРАДИЕНТНЫЕ ЭЛЕМЕНТЫ МИКРООПТИКИ ДЛЯ ДОСТИЖЕНИЯ СВЕРХРАЗРЕШЕНИЯ

Виктор Викторович Котляр1'2 (заведующий лабораторией, e-mail: kotlyar@smr.ru), Алексей Андреевич Ковалев1'2 (научный сотрудник, e-mail: alanko@smr.ru), Антон Геннадьевич Налимов1'2 (научный сотрудник, e-mail: anton@smr.ru) 1 Самарский государственный аэрокосмический университет имени академика С.П. Королева, 2 Учреждение Российской академии наук Институт систем обработки изображений РАН

Аннотация

Получены модовые решения уравнения Гельмгольца для произвольного градиентного планарного волновода. Амплитуда моды представлена как экспонента с показателем в виде ряда Тейлора, коэффициенты которого находятся из рекуррентных соотношений. Показано, что минимальная ширина моды в квадратичном и секонсном градиентных планарных волноводах равна 0,4 от длины волны в вакууме, деленной на значение показателя преломления вдоль оси волновода. С помощью моделирования программой FullWAVE показано, что градиентные микролинзы Микаэляна и «рыбий глаз» Максвелла могут формировать изображения со сверхразрешением. Для 2D микролинзы Микаэляна показано, что точечный источник света изображается вблизи поверхности линзы в виде светового пятна с шириной по полуспаду интенсивности 0,12 от длины волны в вакууме. Это меньше, чем дифракционный предел для кремния с показателем преломления 3,47, который равен 0,144 от длины волны в вакууме. Также показано, что микролинза Микаэляна разрешает по полуспаду два близких точечных источника, разделенных расстоянием 0,3 от длины волны в вакууме.

Ключевые слова: сверхразрешение, дифракционный предел, градиентный волновод, пла-нарный волновод, модовое решение, ширина моды, микролинза Микаэляна, «рыбий глаз» Максвелла.

Введение

Фокусировать когерентный лазерный свет в пятно меньшее дифракционного предела или изображать точечный источник света со сверхразрешением стало возможно в последние годы с развитием микрооптики и нанофотоники. Преодолевать дифракционный предел, равный 0,51Х/П, где X - длина волны света в вакууме, п - показатель преломления вещества в месте фокуса, можно с помощью суперлинз [1]. Суперлинза - это 2Б плоская пластина из метаматериала, который представляет собой чередующиеся субволновые слои металла и диэлектрика. Диэлектрические проницаемости таких слоев подбирают таким образом, чтобы эффективный показатель преломления композитного вещества был равен п = -1. Эксперименты по сверхразрешению с помощью такой линзы были осуществлены в [2, 3]. В эксперименте было достигнуто сверхразрешение, равное 0,4Х [2]. Аналогичный эксперимент был проведен в [4] с помощью субволнового слоя серебра, который также работает как суперлинза. В этой работе были разрешены две линии, разделенные расстоянием 145 нм при освещении их ультрафиолетовым светом с длиной волны 365 нм. То есть в [4] было также достигнуто сверхразрешение, равное 0,4Х.

В [5] была реализована гиперлинза, которая изображает со сверхразрешением в дальней зоне. При этом были разрешены две линии шириной по 35 нм, разделенные расстоянием 150 нм, с помощью света с длиной волны 365 нм. То есть опять было достигнуто сверхразрешение равное 0,4 X. Хотя при моделировании гиперлинзы в виде решетки в [6] было получено сверхразрешение равное 0,05Х на расстоянии 1,5Х от поверхности. По-видимому, при моде-

лировании не учитывается поглощение и рассеяние света метаматериалом, которое имеет место в реальных экспериментах. Это косвенно подтверждают результаты работы [7]. В [7] исследовалась фокусировка лазерного света с помощью зонной пластинки из золотой пленки толщиной 100 нм. Расчет показал, что диаметр фокусного пятна по полуспаду интенсивности должен быть равен 0,35Х, а в эксперименте было получено, что диаметр фокусного пятна по полуспаду интенсивности равен X.

В [8-10] численно анализируются многослойные и анизотропные наноструктуры, которые позволяют получить субволновое разрешение. Так, в [8] подобраны параметры 8-ми слойной Ш структуры на основе А^Ю2, которая имеет толщину 400 нм и изображает источник шириной 0,4Х в световое пятно такой же ширины. В [9] предложена анизотропная 2Б наноструктура, у которой составляющие тензора диэлектрической проницаемости имеют значения ех = 0,01-/0,01 и ег = -100, а толщина структуры вдоль оси г равна 400 нм (длина волны X = 700 нм). Моделирование показало, что такая структура разрешает две линии толщиной по 3 нм, разделенные расстоянием 23 нм. То есть было достигнуто сверхразрешение, равное 0,03Х. В [10] теми же авторами предложена другая 2Б наноструктура, состоящая из двух разных анизотропных слоев, которая разрешает две узкие щели, разделенные расстоянием 50 нм, при освещении их светом с длиной волны 1550 нм.

Еще одним кандидатом на достижение сверхразрешения являются фотонные кристаллы. В [11] с помощью моделирования показано, что 2Б фотон-но-кристаллический слой из диэлектрика с проницаемостью е = 12 и треугольной решеткой круглых отверстий с радиусом г = 0,4а, а - период решетки

отверстий, имеет эффективный показатель преломления п = -1 и является изображающей линзой для длины волны X = а/0,3. При этом точечный источник изображается в виде пятна диаметром по полуспаду интенсивности, равным 0,3Х, а два точечных источника разрешаются по критерию 20%, если расстояние между ними равно 0,5Х.

В недавно проведенных экспериментах по использованию 2Б фотонно-кристаллического слоя в качестве суперлинзы было получено изображение точечного источника шириной 0,4Х [12]. Появились также работы, в которых предлагается усовершенствовать суперлинзу с помощью специальной наноразмерной оболочки [13] или с помощью градиентной границы для слоя вещества с отрицательной рефракцией [14]. В [15] экспериментально продемонстрировано, что слой анизотропного вещества способен усиливать и преобразовывать поверхностные неоднородные волны в распространяющиеся световые волны.

В данной работе численно показано, что хорошо известные в оптике градиентные линзы Микаэляна и «рыбий глаз» Максвелла также могут рассматриваться как кандидаты на получение изображений со сверхразрешением. Для 2Б микролинзы Микаэляна показано, что точечный источник света изображается вблизи поверхности линзы в виде светового пятна с шириной по полуспаду интенсивности БЩИМ = 0,12Х. Это меньше, чем дифракционный предел для кремния (п = 3,47): 0,5Х/п = 0,144 X, и меньше, чем сообщалось в работах [7, 8, 11, 12]. Также показано, что микролинза Микаэляна разрешает по полуспаду два близких точечных источника, разделенных расстоянием 0,3 X, что меньше, чем сообщалось в работах [2, 4, 5, 11].

В [16,17] получены аналитические соотношения для расчета мод в градиентном планарном волноводе. В [16] на основе аппроксимации непрерывной функции показателя преломления градиентного волновода кусочно-постоянной функцией и на основе применения метода передаточных матриц записаны рекуррентные уравнения для поиска неопределенных постоянных локальных мод в каждом слое с постоянным показателем преломления. В [17] на основе модифицированных функций Эйри и метода ВКБ получены уравнения на собственные значения для расчета параметров модовых функций градиентной среды. Общим ограничением в этих работах является предположение о наличии одной [16] или двух [17] точек возврата, которые имеет функция распределения показателя преломления планарного волновода.

В данной работе получены более общие аналитические выражения для модовых решений в пла-нарном волноводе без ограничений на число точек возврата в функции показателя преломления.

1. Решение уравнения Гельмгольца для 2Б градиентного волновода

На рис. 1 показана схема задачи. Рассмотрим градиентную двумерную среду с показателем пре-

ломления п = п(х), в которой распространяется электромагнитная волна с ТЕ- поляризацией, вектор напряженности электрического поля которой направлен вдоль оси у.

XI

Еу(х,г)

Рис. 1. Схема задачи распространения ТЕ-поляризованной волны в 2Б градиентном волноводе

Амплитуда электрического вектора Еу(х,г) удовлетворяет уравнению Гельмгольца [18]:

э7+кп( х)

Еу (х, г) = 0,

(1)

где к - волновое число. Разложим амплитуду электромагнитной волны по поперечным модам градиентной среды:

Еу (х, г) = X Сп (х) ехр(/Дг):

(2)

где вп = кгп - постоянная распространения п-ой моды. Из представления (2) следует, что у светового поля имеется продольный период Т, так что вп = 2пп/Т. Например, для градиентной среды с квадратичной зависимостью показателя преломления модами являются функции Эрмита-Гаусса, составляющие счетный базис [19].

Подставляя (2) в (1), получим уравнение для амплитуды мод градиентной среды:

д 1(Сп(х)

дх

- + Р„ (х)Си (х) = 0,

где

Рп (х) = к 2 «2(х) "А»2. Заменой переменных

Си (х) = Си (0)ехр

} / (Х)дХ

(3)

(4)

(5)

уравнение (3) сводится к нелинейному дифференци альному уравнению Уиттекера [20] для функции/,^:

д/п (х)

Лх

- + /Д х) + Рп (х) = 0.

(6)

Уравнение (6) может быть решено с помощью разложения в ряд Тейлора функций /п(х) и Рп(х):

/п (х) = X С^ хт

Рп (х) = X Р.

<Х> хт

т

(7)

(8)

п= 0

0

т= 0

т=0

где СтП и ртП неизвестные и известные коэффициенты разложения соответствующих функций. Подставляя (7) и (8) в (6), для неизвестных коэффициентов ряда (7) получим рекуррентные соотношения (т > 0):

сг =-;->;-. + X с*ст*,

(9)

П(п)

где п - положительное целое число, С0 - неопределенные константы. Тогда для амплитуд мод градиентной волноводной среды можно записать явное выражение:

(

Сп (х) = Сп (0)ехр

X С

(п)

х

т+1 Л

т +1

(10)

в котором Ст(п) находятся из рекуррентных соотно-

шений (9). Коэффициенты р находятся по формулам:

(п) = атп\х)

Рт т! с1хт

(п)

т '

входящие в (8) и (9),

(11)

(12)

при т > 0 и при т = 0 из соотношений

р0п) = к 2п 2(0)-0

Окончательное выражение для амплитуды светового поля ТЕ-поляризованной волны в градиентной 2Б среде имеет вид:

Еу (х, г) = £ Сп (0) ехр

0+X с;

(п)

т +1

(13)

В (13) для каждой моды имеются две неопределенные константы Сп (0) и С0(п). Моды (10) не являются ортогональными и нормированными, поэтому для разложения поля (13) по этим модам (то есть нахождения коэффициентов Сп (0)) потребуется

усечение обоих рядов в (13) до конечных сумм и решение системы линейных алгебраических уравнений.

Постоянные С0(п) для каждой моды нужно выбирать специальным образом. Для примера рассмотрим два частных случая.

1) Пусть п(х) = по и моды распространяются в однородном пространстве. Из (11) и (12) получим,

что р0п) = к2п2(0)-Ь2 при т = 0 и р^ = 0 при т > 0. Из (9) следует, что

С,(п) =-( р0п) + (С0(п))2)

Пользуясь произволом в выборе С

С„(п) =

(п) , 0

= ф п -0

(14) положим:

(15)

00

Тогда из (14) следует, что С1(п) = 0, а из (9) следует, что и все остальные коэффициенты равны ну лю: Ст(п) = 0 при т > 0. Итак, для однородного про странства п(х) = п0 вместо (13) получим:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Еу (х,г) = £ Сп(0)ехр 0 + ¡ф2-00

(16)

Решение (16) представляет собой хорошо известное решение уравнения Гельмгольца (1) в виде суперпозиции плоских волн.

2) Рассмотрим другой частный случай: п (х) = Щ - а2х - квадратичная среда. В этом случае

, Ч (п) , (и) (и) 2 (п) ,22 а2

рп (х) = р0; + Р1 х + р2;x, где р0; = к п - Рп,

С0(п) = 0 , (п) =

р1(п) = 0 , р2(п) = -к2а2 и рт(п) = 0 при т > 2. Пусть тогда С2(п) = - (С0(п)С1(п)) = 0 , а С1(п) = -р0п) = 00 - к2 Щ. Чтобы все остальные коэффициенты были равны нулю С^ = 0 при т > 2, достаточно положить равным нулю третий коэффициент: С3(п) = -р2п) / 3 - (С1(п) )2 / 3 = 0 . Тогда на параметр среды а получим условие: к2а2 = (к2щ -0,2)2. Положим

0п = кп0142 , тогда получим, что а = кп^ / 2 . Итак, получается, что модой квадратичного волновода с показателем преломления п2(х) = п^ (1 - к2 п^х2 /4) является гауссовая экспонента:

Е (х,г) = С(0)ехр

^ /кп0 г к2 п х2 ^

72

4

(17)

В заключении этого раздела отметим, что так как производная в нуле модовой функции (10) пропорциональна нулевому коэффициенту в ряде (7): СС(х)

Сх

= Сп (0)^

(18)

то выбором С0(п) = 0 (как это было сделано в последнем примере) обеспечивается наличие экстремума у моды на оптической оси.

2. Общее решение для секонс-градиентного волновода

Для градиентного волновода с секонсной зависимостью показателя преломления от поперечной координаты

п(х) = п0 сИ 11 —п°х

(19)

известно частное решение уравнения (1) следующего вида [21]:

„ ( 1кпп2 | , , ( кппх \

Е„<х,„ = ехр[-2-]-п.] . 00)

Уравнение Гельмгольца (1) в этом случае имеет

вид:

Э2 + Э2 2 , -2 (кп0х

эх2+кп"сЬ [^Т

Еу (х, г) = 0.

Будем искать общее решение (21) в виде Е1у (х, г) = А( х)ехр(/уг).

(21)

(22)

п= 0

т-1

1= 0

т=0

х=0

V

у

х

п=0

т=0

х=0

Подставляя (22) в (21), получим уравнение

д2 А( х) ск2

+ я (х)А(х) = 0 ,

где

Я (х) = к 2п1 сЬ

-2 | кп0 х | 72

-Л

(23)

(24)

Известно [20], что уравнение вида (23) имеет общее решение

А(х) = А (х)

С + С2} А"2(%)д%

(25)

где А1(х) - частное решение (23), С1, С2 - неопределенные константы. В нашем случае в качестве частного решения может быть выбрано решение (20), то есть

А( х)=1 кп2х

7 =

кп„

72'

(26)

Тогда общее решение для моды секонс-градиентного волновода будет иметь вид:

Ех (х, г) = ехр

1кп0 г

7Т

сЬ "

11 х | /

IX

X <! С, +

С

1 кппу[2

— $Ь(кп0

кп0 х

Ж

(27)

где

С, = К (х = 0, г)

С2 =

СЕ. (х = 0, г)

сх

(28)

Из (20) можно получить выражение для ширины моды по полуспаду интенсивности в секонс-градиентном волноводе:

1п(3 + 2^2) 0,41

^ГЯМ = -

лп,

л/2

(29)

"0

где X - длина волны света в вакууме, пс - как и прежде, показатель преломления на оси волновода.

3. Частное решение для квадратичного волновода В первом разделе данной работы показано, что для квадратичной среды с определенными параметрами модой является гауссовая экспонента (17). В этом разделе покажем, что это остается верным и для квадратичной среды с произвольными параметрами:

'(х) = Я! (1" х2 ):

(30)

где V - произвольная постоянная. Тогда уравнение (1) будет иметь вид:

-\2 г\2

Э , г 2 2 л 2 2ч

э7+э7+кп0(1" "х)

Еу (х, г) = 0,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

а его решение будем искать в виде: Е2у (х,г) = Е0ехр (р - q2х2) .

(31)

(32)

Подставляя (32) в (31), получим уравнение:

(^2 -Р2 + к2п02 + 4q4х2 - V2к2п02х2)X

хЕ0 ехр(р - q2х2) = 0.

Из (33) определим неизвестные параметры:

wknn

(33)

2

q =

2

Р=кп-'1 - кп:

Таким образом, частное модовое решение уравнения (31) имеет вид:

( I- > А

Е2 у (х, г) = Е0 ехр

¿кп0 2. 11 -

wkn,

кп0 2

0 х2

(34)

Заметим, что при V = кпс/2 решение (34) совпадает с решением (17). При V = кпс/2 из (34) следует, что ширина (диаметр) гауссовой моды по полуспаду интенсивности равна:

Т1П4я 0.381

FWHM = -

лпп

Сравнивая (35) с (29), видим, что ширина обеих мод (секонсной и гауссовой) почти одинаковая. Эффективная ширина квадратичного волновода находится из условия п(хс) = 1 и равна

2^

2 х =-

0

/п2 -1

-1.

(36)

При пс = 1,5 из (36) получим: 2хс ~ 0,48Х. То есть эффективная ширина стеклянного планарного волновода с квадратичной зависимостью показателя преломления, в котором может распространяться только гауссовая мода (34), почти равна половине длины волны света в вакууме. Заметим, что такие «полуволновые» волноводы в настоящее время активно используются в приложениях [22, 23].

4. Моделирование прохождения света через градиентные элементы микрооптики: сверхразрешение Несколько лет назад были проведены эксперименты по сверхразрешению в оптическом диапазоне с помощью суперлинз [3,5]. В [3] в качестве суперлинзы использовался тонкий слой серебра толщиной 50 нм. Эксперимент состоял в получении изображения со сверхразрешением в слое резиста амплитудной решетки с периодом 145 нм вблизи тонкой пленки серебра. После освещения ультрафиолетовым светом с длиной волны 365 нм трехслойной структуры (амплитудная решетка + пленка серебра + резист) и после проявления, в слое резиста «записывалась» и была восстановлена с помощью преобразования Фурье решетка с периодом 170 нм. В [5] была реализована гиперлинза, представляющая собой полуцилиндр с кольцевым вырезом на подложке из кварца, состоящий из чередующихся слоев серебра (толщиной 35 нм) и диэлектрика А1203 (толщиной 35 нм). Такая многослойная структура имеет анизотропную диэлектрическую проницаемость. Объектом разрешения служили две амплитудные

0

п

0

п

линии, полученные литографическим способом, имеющие ширину 35 нм и разделенные расстоянием 150 нм. Этот объект освещался ультрафиолетовым светом с ТЕ-поляризацией и длиной волны 365 нм. С помощью обычной линзы было получено увеличенное изображение этих двух линий, разделенных уже расстоянием 350 нм. Таким образом, в описанных экспериментах [3,5] было достигнуто сверхразрешение равное 0,4Я.

Теоретически (без учета поглощения материала) с помощью суперлинзы может быть достигнуто любое разрешение. Так, в [9] с помощью моделирования показано, что гиперболическую линзу, формирующую в ближнем поле увеличенное субволновое изображение, можно реализовать и в виде плоскопараллельного слоя. Было показано, что слой анизотропного вещества с диэлектрическими проницае-мостями ех = 0,01-/0,01 и е7 = -100 толщиной 400 нм разрешает две щели шириной по 3 нм (вытянутые вдоль оси у) в металлическом экране с диэлектрической проницаемостью е = 1-/104, разделенных расстоянием 23 нм, если осветить их светом с ТМ-поляризацией и длиной волны 700 нм. Величину сверхразрешения можно оценить как 0,05Я.

В следующих подразделах будет численно показано (хорошо известным методом РБТО), что с помощью градиентной 2Б микрооптики также можно достичь сверхразрешения. Моделирование проводилось с помощью программы Ри1^АУБ фирмы Я^й (США).

4.1. Сверхразрешение с помощью микролинзы Микаэляна

Из-за дифракции света в однородном пространстве два близких точечных источника света перестают разрешаться на расстоянии много меньшем длины волны. Так, на рис. 2а показаны исходные (в плоскости 7 = 5 нм) пять гауссовых когерентных источников света шириной Я/200 каждый, которые разделены расстоянием Я/50. На рис. 2б показана интенсивность света от этих источников на расстоянии 7 = 30 нм для длины волны Я = 1550 нм. Из рис. 2 видно, что близкие точечные источники перестают разрешаться на расстоянии, примерно равном промежутку между ними (7 = Я/50).

На рис. 3 показан результат моделирования формирования изображения двух точечных источников света (ширина каждого из них по 35 нм), разделенных расстоянием 150 нм, с помощью линзы Микаэляна [24, 25].

Показатель преломления такой секонс-градиентной 2Б микролинзы имеет вид:

п(х) = п0оЬр) , (37)

где 2Ь - длина линзы. Ширина линзы 2Я находится из условия, что п(Я) = 1. Для линзы на рис. 3 показатель преломления на оси п0 = 2,1, ширина линзы 2Я = 1 мкм, длина линзы 2Ь = 1,144 мкм, длина волны Я = 365 нм. Выбранные параметры совпадают с параметрами эксперимента в работе [5].

гауссовых источников света шириной Я/200, разделенных интервалом Я/50 на разных расстояниях: 7 = 5 нм (а), 7 = 30 нм (б), Я = 1550 нм (по оси ординат отложены произвольные единицы)

МКМ 1,8

0,9

0,0

Рис. 3. Схема расположения 2П градиентной микролинзы Микаэляна в окне программы ЕыИЖЛУЕ для моделирования изображения двух близких точечных источников света.

Показана мгновенная амплитуда Еу(х,г) в момент времени, когда свет прошел расстояние в 18,2398 мкм

На рис. 3 показана мгновенная картина амплитуды электрической составляющей электромагнитной волны с ТЕ-поляризацией (свет распространяется на рис. 3 снизу вверх). Видно, что на противоположной по отношению к источникам света стороне линзы формируется их изображение. На рис. 4а показано распределение усредненной по времени интенсивно-

I I2

сти электрического поля I (х, 7 = 70) = \Еу (х, 7 = 70)

Изображения Микролинза

-1 0 +1 х, мкм

непосредственно на «задней» стороне линзы, то есть на расстоянии г0 = 2Ь от «передней» стороны линзы. Два источника света находились перед «передней» стороной линзы на расстоянии 20 нм.

0,035

а)

0,012

Рис. 4. Распределения усредненной интенсивности электрической составляющей ТЕ-поляризованной световой волны, рассчитанные на расстоянии г = С (а) и г = 1С нм (б) от «задней» поверхности линзы Микаэляна

(рис. 3), на входе у которой расположены два близких точечных источника света (горизонтальная линия делит интенсивность по полуспаду, по оси ординат отложены произвольные единицы)

Из рис. 4а видно, что два точечных источника, разделенных расстоянием 150 нм (между центрами источников 180 нм), уверенно разрешаются. Если плоскость регистрации интенсивности удалить от «задней» плоскости линзы на расстояние г = 10 нм (рис. 4б), то источники продолжают разрешаться, хотя величина «провала» в распределении интенсивности слегка уменьшается. Величина разрешения, достигнутая на рис. 4, равна 0,41Х. Из рис. 4а можно также определить, что ширина изображения линзой Микаэляна точечного источника по полуспаду интенсивности равна БЩИМ = 100 нм = 0,27Х. Полученное значение величины разрешения 0,41Х почти равно разрешению, полученному в [2,3,5] (0,4Х). Заметим, что полученная ширина изображения 0Д7Х точечного источника согласуется с минимальной шириной моды, способной распространяться в секонс-градиентном волноводе 0,4Х/пс = 0,27Х при пс = 1,5 (см. ур. (29)).

С целью увеличения разрешения линзы Мика-эляна, как следует из ур. (29) и (35), был увеличен показатель преломления на оси пс = 3,47 (кремний,

Другие параметры моделирования тоже изменились (рис. 5а): длина волны X = 1 мкм, ширина линзы 2К = 6 мкм, длина линзы 2Ь = 4,92 мкм.

г, мкм гг ^

Изо бражен ие

+4 х, мкм

0,3

0,0

аШД^

X, мкм

б) -3 0 +3

Рис. 5. Мгновенная картина амплитуды электрического поля ТЕ-поляризованной световой волны в линзе Микаэляна с точечным источником на передней плоскости линзы (а) и усредненное распределение интенсивности в задней плоскости линзы (б) (по оси ординат отложены произвольные единицы)

Шаг моделирования по пространственным осям во всех примерах равен Х/100. Ширина гауссового точечного источника света в начальной плоскости равна Х/20. На рис. 5а показана мгновенная картина амплитуды электрического вектора световой волны в линзе Микаэляна, рассчитанной в момент времени, когда световая волна от источника прошла 200 мкм. На рис. 5б показано распределение усредненной интенсивности света ТЕ-поляризованной волны на выходе из линзы (рис. 5а), на «задней» плоскости линзы Микаэляна. Расчет показал, что для интенсивности, показанной на рис. 5б, ширина центрального максимума по полуспаду в изображении точечного источника равна БЩИМ = 0,12Х. Однако интенсивность (или плотность мощности света) не является той величиной, которая показывает, какая часть мощности излучения источника распространяется в пространстве вдоль оси г. Поэтому на рис. 6 как раз показана такая величина: проекция на оптическую ось вектора Умова-Пойнтинга, рассчитанного на выходе из линзы Микаэляна (рис. 5а), на входе которой был расположен точечный источник света.

0,020

-0,005

X, мкм

-н

-3 0 +3

Рис. 6. Распределение вдоль оси х величины проекции на

ось 2 вектора Умова-Пойнтинга (в относительных единицах), рассчитанного на выходе линзы Микаэляна (рис. 5а) с входным точечным источником

Ширина центрального максимума потока мощности вдоль оси 2 на рис. 6 такая же, как и на рис. 56 и равна FWHM = 0,12Х.

Величина дифракционного предела, которого можно достичь при фокусировке света в однородной среде, известна и равна 0,5Х/п, где п - показатель преломления однородной среды. Для нашего случая п0 = 3,47, поэтому дифракционный предел равен FWHM = 0,144Х. Таким образом, из рис. 56 и рис. 6 следует, что с помощью микролинзы Микаэляна (рис. 5а) можно преодолеть дифракционный предел не только в свободном пространстве, но и в среде. Причина этого, на наш взгляд, заключается в том, что при фокусировке света на границе раздела двух сред в формирование фокуса дают вклад и поверхностные неоднородные волны, которые могут формировать интерференционные и дифракционные картины с периодом много меньшим длины волны [9, 10]. То, что поверхностные волны присутствуют при формировании изображения точечного источника, видно из рис. 6. Проекция вектора Умова-Пойнтинга на ось 2 в некоторой области на оси х отрицательна, то есть вблизи выходной поверхности линзы Микаэляна свет не только выходит из линзы, но и заходит в нее, то есть вдоль поверхности микролинзы распространяется поверхностная волна. При распространении света в свободном пространстве при удалении от поверхности линзы ширина изображения точечного источника быстро увеличивается, и на расстоянии от поверхности, примерно равным длине волны, ширина изображения достигает дифракционного предела 0,5 X.

Заметим, что смена знаков проекции вектора Умова-Пойнтинга, аналогичная рис. 6, ранее наблюдалась в работе [8], была названа оптическим вихрем и интерпретирована как интерференция распространяющейся волны и усиленной поверхностной волны. Однако в [8] рассматривалась Ш многослойная структура (Ш фотонный кристалл).

На рис. 7а показана мгновенная картина амплитуды электрического вектора световой волны в линзе Микаэляна (параметры такие же, как для рис. 5а), когда перед ее входной (нижней на рис. 7а) поверх-

ностью на расстоянии 10 нм находятся два точечных источника (шириной по 50 нм), разделенных расстоянием 300 нм. На рис. 76 показано усредненное распределение проекции на оптическую ось вектора Умова-Пойнтинга, рассчитанное за выходной поверхностью микролинзы на расстоянии 10 нм.

. г, мкм *

1 ' Изображения

0,30-

0,00

-0,20.

щ

г

6)

о

х, мкм -н

+5

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Рис. 7. Мгновенная картина амплитуды электрического вектора световой волны с ТЕ-поляризацией в линзе Микаэляна (рис. 5а), когда перед ее входной (нижняя горизонтальная линия на рисунке) поверхностью на расстоянии 10 нм находятся два точечных источника (шириной по 50 нм), разделенных расстоянием 300 нм (а); усредненное по времени распределение проекции вектора Умова-Поинтинга на оптическую ось, рассчитанное на расстоянии 10 нм от задней (верхняя горизонтальная линия на рисунке) плоскости линзы (по оси ординат отложены произвольные единицы) (6)

Из рис. 76 видно, что два источника разрешаются, при этом величина сверхразрешения равна 0,3Х, что меньше, чем в работах [2,4,5,11].

4.2. Сверхразрешение с помощью микролинзы «ры6ий глаз» Наряду с рассмотренной в предыдущем разделе микролинзой Микаэляна можно использовать и другие градиентные изображающие оптические элементы, зависимость показателя преломления от координат которых получена в явном аналитическом виде.

Одним из таких оптических элементов является «рыбий глаз» Максвелла [18,26,27], функция показателя преломления которого в двумерном случае в полярных координатах имеет вид:

п(Г) = п0

1+1 Я

(38)

где пс - показатель преломления в центре круга, Я -радиус элемента. Из (38) видно, что при г = Я показатель преломления уменьшается в 2 раза. Недостатком этого элемента является ограниченный перепад показателя преломления в центре и на краю оптического элемента: он не может изменяться больше, чем в 2 раза. В линзе Микаэляна перепад показателя преломления определяется только материалом. Например, для кремния показатель меняется от 3,47 до 1. Однако «рыбий глаз» обладает круговой симметрией, и поэтому любой точечный источник, расположенный на его поверхности, изображается «идеально» в диаметрально противоположенной точке на его поверхности.

На рис. 8а показано мгновенное распределение амплитуды электрического поля ТЕ-поляризованной световой волны внутри 2Э микролинзы «рыбий глаз», на поверхности которой находятся два точечных близких источника света.

г, мкм Тх ¿.

Изображения

Рис. 8. 2П микролинза «рыбий глаз» в окне программы

FullWAVE и мгновенная картина амплитуды электрического поля ТЕ-поляризованной волны внутри линзы, у которой на входе два точечных источника (свет

распространяется снизу вверх) (а), а также распределение усредненной интенсивности в плоскости изображения линзы (по оси ординат отложены произвольные единицы) (б)

Параметры моделирования: показатель преломления в центре линзы пс = 3,47, радиус линзы Я = 2,5 мкм, длина волны света Х = 1 мкм, ширина

точечных источников по 0,05Х, и они разделены расстоянием 440 нм или 0,44Х. На рис. 8б показано распределение усредненной по времени интенсивности в плоскости изображения (верхняя горизонтальная линия на рис. 8а). Видно, что два источника разрешаются (разрешение равно 0,44Х по критерию Рэлея 20%). Таким образом, данная микролинза незначительно (0,44Х) преодолевает дифракционный предел в разрешении (0,5Х).

При изображении отдельного точечного источника с помощью данной микролинзы (рис. 8а) получается распределение интенсивности, показанное на рис. 9. Ширина центрального максимума интенсивности по полуспаду равна БЩИМ = 0,3Х. Из сравнения рис. 9 и рис. 5б можно заключить, что микролинза «рыбий глаз» формирует более широкое изображение (примерно в 2 раза шире) точечного источника, чем микролинза Микаэляна при сравнимых параметрах. Хотя при этом величины разрешения 0,44Х (рис. 8б) и 0,3 Х (рис. 9) сравнимы с величинами сверхразрешения, полученными в [2, 4, 5, 8, 11, 12].

0,008

0,004

-1 0 +1 Рис. 9. Усредненное по времени распределение интенсивности электрического поля в плоскости изображения микролинзы «рыбий глаз» (рис. 8а),

у которой на поверхности расположен только один точечный источник света (по оси ординат отложены произвольные единицы)

Заключение

В работе получены следующие результаты:

- найдены модовые решения уравнения Гельм-гольца для произвольного градиентного планарного волновода; амплитуда моды представлена как экспонента с показателем в виде ряда Тейлора, коэффициенты которого находятся из рекуррентных соотношений (уравнения (9), (10) и (13));

- показано, что минимальная ширина моды в квадратичном и секонсном градиентных планарных волноводах равна 0,4 от длины волны в вакууме, деленной на значение показателя преломления вдоль оси волновода (уравнения (29) и (35));

- с помощью моделирования программой Ги11-ШАУБ показано, что градиентные микролинзы Мика-эляна и « рыбий глаз» Максвелла могут формировать изображения со сверхразрешением (рисунки 7 и 8);

- для 2Э микролинзы Микаэляна показано, что точечный источник света изображается вблизи задней поверхности линзы в виде светового пятна с шириной

-1

2

по полуспаду интенсивности 0,12 от длины волны в вакууме; это меньше, чем дифракционный предел для кремния с показателем преломления 3,47, который равен 0,144 от длины волны в вакууме, и меньше, чем получено в работах [8,11,12] (рисунки 5 и 6);

- также показано, что микролинза Микаэляна разрешает по полуспаду два близких точечных источника, разделенных расстоянием 0,3 А, это меньше дифракционного предела в свободном пространстве, равного

0.5. и меньше, чем в работах [2,4,5,11] (рисунок 7).

Благодарности

Работа выполнена при поддержке российско-американской программы «Фундаментальные исследования и высшее образование» (грант CRDF PG08-014-1), гранта Президента РФ поддержки ведущих научных школ (НШ-3086.2008.9) и гранта РФФИ 08-07-99007.

Литература

1. Pendry, J.B. Negative refraction makes a perfect lens / J.B. Pendry // Phys. Rev. Lett. 2000. - Vol. 85(18). - P. 3966-3969.

2. Blaikie, R.G. Imaging through planar silver lenses in the optical near field / R.G. Blaikie, D.O.S. Melville // J. Opt A: Pure Appl. Opt. 2005. - Vol. 7(2). - P. S176-S183.

3. Melvile, D.O.S. Super-resolution imaging through a planar silver layer / D.O.S. Melvile, R.J. Blaikie // Opt. Express 2005. - Vol. 13(6). - P. 2127-2134.

4. Fang, N. Sub-diffraction-limited optical imaging with a silver superlens / N. Fang, H. Lee, C. Sun, X. Zhang // Science 2005. - Vol. 308(5721). - P. 534-537.

5. Liu, Z. Far-field optical hyperlens magnifying sub-diffraction-limited object / Z. Liu, H. Lee, Y. Xiong, C. Sun, X. Zhang // Science 2007. - Vol. 315(5819). - P. 1686.

6. Thongrattanasiri, S. Hypergratings: nanophotonics in planar anisotropic metamaterials / S. Thongrattanasiri, V.A. Podolskiy // Opt Lett 2009. - Vol. 34(7). - P. 890-892.

7. Fu, Y. Experimental study of plasmonic structures with varint period for sub-wavelength focusing: analysis of characterization errors / Y. Fu, R.G. Mote, Q. Wang, W. Zhou // J. Mod. Opt. 2009. - Vol. 56(14). - P. 1550-1556.

8. Webb, K.J. Subwavelength imaging with a multilayer silver film structure / K.J. Webb, M. Yang // Opt. Lett. 2006.

- Vol. 31(14). - P. 2130-2132.

9. Liu, H. Submevelength imaging opportunities with planar uniaxial anisotropic lenses / H. Liu, Shivananad, К J. Webb // Opt Lett 2008. - Vol. 33(21). - P. 2568-2570.

10. Liu, H. Subwavelength imaging with nonmagnetic anisotropic bilayers / H. Liu, Shivanand, K.J. Webb // Opt Lett. 2009. - Vol. 34(14). - P. 2243-2245.

11. Wang, X. Unrestricted superlensing in a triangular twodi-mensional photonic crystal / X.Wang, Z.F. Ren, K Kempa // Opt Express 2004. - Vol. 12(13). - P. 2919-2924.

12. Casse, B.D.F. Imaging with subwavelength resolution by a generalized superlens at infrared wavelengths / B.D.F. Casse, W.T. Lu, R.K. Banyal, Y.J. Huang, S. Selvarasah, M.R. Dokmeci, C.H. Perry, S. Sridhar // Opt. Lett. 2009. -Vol. 34(13). - P. 1994-1996.

13. Tsukerman, I Superfocusing by nanoshells / I. Tsukerman // Opt Lett 2009. - Vol. 34(7). - P. 1057-1059.

14. Ingrey, P.C. Perfect lens with not so perfect boundaries / P.C. Ingrey, K.I. Hopcraft, O. French, E. Jakeman // Opt. Lett. 2009. - Vol. 34(7). - P. 1015-1017.

15. Ray, E.A. Simple demonstration of visible evanescent-wave enhancement with far-field detection / E. A. Ray,M.

J. Hampton, R. Lopez // Opt. Lett. 2009. - Vol. 34(13). -P. 2048-2050.

16. Cao, Z. Exact analytical method for planar optical waveguides with arbitrary index profile / Z. Cao, Y. Jiang, Q. Shen, X. Dou, Y. Chen // J. Opt. Soc. Am. A 1999. -Vol. 16(9). - P. 2209-2212.

17. Chung, M. General eigenvalue equations for optical planar waveguides with arbitrarily graded-index profiles / M. Chung, C. Kim // J. Lightwave Techn. 2000. - Vol. 18(6).

- P. 878-885.

18. Борн, М. Основы оптики / M. Борн, Э. Вольф - М.: Наука, 1973. - 719 с.

19. Миллер, У. Симметрия и разделение переменных / У. Миллер - М.: Мир, 1981. - 342 с.

20. Корн, Г. Справочник по математике / Г. Корн, Т. Корн

- М.: Наука, 1968. - 720 с.

21. Триандафилов, Я.Р. Фотонно-кристаллическая линза Микаэляна / Я.Р. Триандафилов, В.В. Котляр // Компьютерная оптика. - 2007. - Т. 31, № 3. - С. 27-31. -ISSN 0134-2452.

22. He, J. Wavelength switchable semiconductor laser using half-wave V-coupled cavities / J. He, D. Liu // Opt. Express 2008. - Vol. 16(6). - P. 3896-3911.

23. Lin, X. Design and analysis of 2x2 half-wave waveguide couplers / X. Lin, D. Liu, J. He // Appl. Opt. 2009. - Vol. 48(25). - P. F18-F23.

24. Микаэлян, АЛ. Применение слоистой среды для фокусирования волн / А. Л. Микаэлян // Доклады академии наук СССР. - 1951. - Т. LXXXI. - С. 569-571.

25. Котляр, В.В. Преобразование Абеля в задачах синтеза градиентных оптических элементов / В.В. Котляр,

A.С. Мелехин // Компьютерная оптика. - 2001. - №. 22. - С. 29-36. - ISSN 0134-2452.

26. Котляр, В.В. Преобразование Абеля для расчета градиентных оптических элементов со сферически-симметричным распределением показателя преломления /

B.В. Котляр, А.С. Мелехин // Компьютерная оптика. -2002. - № 24. - С. 48-52. - ISSN 0134-2452.

27. Котляр, В.В. Расчет обобщенных линз «рыбий глаз» Максвелла и Итона - Липмана / В.В. Котляр, А.С. Мелехин // Компьютерная оптика. - 2002. - № 24. - С. 5357. - ISSN 0134-2452.

References