Анализ дифракции лазерного излучения на аксиконе с числовой апертурой выше предельной

АНАЛИЗ ДИФРАКЦИИ ЛАЗЕРНОГО ИЗЛУЧЕНИЯ НА АКСИКОНЕ С ЧИСЛОВОЙ АПЕРТУРОЙ ВЫШЕ ПРЕДЕЛЬНОЙ

Устинов А.В., Хонина С.Н. Институт систем обработки изображений РАН, Самарский государственный аэрокосмический университет имени академика С.П. Королёва

Аннотация

В работе выполнен теоретический анализ дифракции плоской волны на аксиконе в зоне затухающих волн. Полученные аналитические оценки качественно согласуются с результатами численного моделирования в рамках скалярной непараксиальной волновой модели. Ключевые слова: аксикон, предельная числовая апертура, затухающие волны.

Введение

Уже несколько десятков лет аксикон [1] считается классическим оптическим элементом наравне с линзой. Одним из основных его применений является формирование Бесселевых пучков, обладающих бездифракционными свойствами [2]. В этом случае создаётся световой пучок, сохраняющий постоянный поперечный размер на очень большом расстоянии (значительно превышающем расстояние нерасходимости Гауссова пучка). Такие пучки нашли применение во многих областях [3, 4], включая оптическое манипулирование [5 - 7], оптическую когерентную томографию [8, 9], метрологию [10, 11].

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

При увеличении размеров оптического элемента до нескольких десятков длин волн использование разностных методов расчёта, таких как конечно-разностный метод решения уравнений Максвелла, становится затруднительным в силу их огромной ре-сурсозатратности. Для быстрого приближённого расчёта можно использовать более простые модели, в частности интегральные операторы распространения [21, 22], а также аналитические оценки [23], которые затем уточняются более строгими методами [24, 25].

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

1. Аксикон в рамках скалярной непараксиальной волновой модели

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

Отсутствию внутреннего отражения соответствует условие, что синус угла преломления не превосходит единицы, то есть угол при вершине удовлетворяет неравенству

a0 > atír = arccos

где n - показатель преломления материала аксикона.

Значение числовой апертуры (NA) аксикона определяется углом, под которым лучи, выходящие из аксикона, пересекают оптическую ось:

NA = p = sin b.

Очевидно, что выходящий из рассматриваемого аксикона луч не может идти под более крутым углом к оптической оси, чем atir, соответствующим наступлению полного внутреннего отражения (1). В рамках волновой теории действие аксикона описывается через комплексную функцию пропускания:

t

где k = 2p /1 - волновое число, 1 - длина волны излучения, падающего на элемент.

Если условие (1) не выполняется, то происходит полное внутреннее отражение и в рамках геометро-оптической теории лучи из аксикона выходят сбоку от элемента [26] и формируют расходящиеся, не пересекающие оптическую ось лучи. Однако из волновой теории известно, что при полном внутреннем отражении часть излучения проникает за границу оптического элемента и образует так называемые затухающие волны. Поэтому формально запишем выражение для параметра аксикона в виде комплексного числа (при отсутствии полного внутреннего отражения оно становится вещественным и меньшим единицы):

p = cos a0 (nsin a0 - i<Jn2 cos

Отметим, что при освещении аксикона с острой стороны его параметр всегда будет вещественным, но может превзойти единицу [26]:

p = ncos

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

Заметим, что геометро-оптическая модель является очень полезной для понимания хода лучей и направления распространения энергии. Однако в этом случае проблематично получить общую картину распределения интенсивности и оценки размера формируемых пучков. В этом случае необходимо применять волновую теорию.

'a0 -1

Непараксиальная скалярная модель, основанная на теории Рэлея-Зоммерфельда [27], позволяет получать корректные результаты на очень близких расстояниях (всего несколько длин волн) от апертуры [21, 22]. Интегральное преобразование Рэлея-Зоммерфельда первого типа в декартовых координатах имеет следующий вид: E (и, z) =

е^А-^у,

;Я Е„

2я _

где Е(хУ - входное поле, I= ^¡(U—X)^+(v—У)^+ZF , Е„ - область, в которой задано входное поле.

2. Анализ продольного распределения

При рассмотрении поля на оптической оси (и = V = „) выражение (6) значительно упрощается, и в условиях радиальной симметрии можно записать:

я ехр ('к^г^+г2 Е(„,„, z) = -zj Е„(г)-

'к -

1

7777

гдтг = Е„(„)е'к - Е„( Я)е'к

х

z

Е„ (г)-

^ёг.

%Я77 Т„ 7777

Второе выражение в (7) получено интегрированием по частям, его удобство показано в [28].

Рассмотрим дифракцию плоской волны, ограниченной круговой апертурой радиуса Я, на аксиконе (3). Результаты дифракции на аксиконе существенно зависят от того, превышает параметр аксикона р предельную числовую апертуру для данной среды или нет. Для свободного пространства ЛЛ, соответствующая распространяющимся волнам, должна быть меньше единицы. Для р>1 в ближней зоне важную роль играют затухающие волны. Их влияние мы учтём в дальнейшем анализе.

2.1. Числовая апертура меньше предельной

В работе [28] было получено следующее выражение для комплексного распределения на оптической оси при использовании пропускающего линейного аксикона (р < 1): Е(„,„, z) = ехр ('кг)-

77

ехр ('к7я2 + г2 - 'крЯ) -

-'Лкгр

я ехр ('кл/г2 + г2 - 'крг) -'крг Г -, -- д.г,

Л 7777

где г„ = рг- р2 , 1 = -7772 р(1 - р2) ,

2 г2

2 = А7к/(2г) (1 - р2) , А2 =-

(1 - р2 )3

Последнее слагаемое в (8) соизмеримо с остальными (не считая первого) только при довольно малых г, поэтому в работе [28] его вклад по сравнению с основным членом не рассматривался. В данном исследовании нам необходимо это сделать.

Если г мало, то в интеграле последнего слагаемого выражения (8) показатель экспоненты приближённо представим в виде линейной функции:

я ехр('к7777-'крг) } --

ехр('кр„) } ехр1-к

г„ +А Л/г +

-&г.

7777

ехр

-2 + г2

Если выполнено одно из неравенств Я> 1,5г или Я > 2 я /(кР1), то верхний предел в (9) можно заменить на бесконечный и сделать преобразование:

Я ¥ г„ +А

Г

г„ +А „ „

Интеграл в бесконечных пределах является табличным [29]:

Г ехр(-'кР,г) =

Г„ л/г2 + г2

= К„ (кр^)- '„, 5р[ 1„ (кр^)- Ь„ (кр1 г)], где К„ (х) - функция Макдональда (модифицированная функция Бесселя 3-го рода) нулевого порядка; 1„ (х) -модифицированная функция Бесселя 1 -го рода нулевого порядка; Ь„( х) - модифицированная функция Струве „-го порядка.

Второй интеграл можно, если ограничиться первым лепестком и заменить косинус и синус линейными приближениями, приближённо оценить в виде:

^ехр (-'кУ) ^

ехр (-'кр1г) ^

Лп

3

7г2 + г2 [

я

+ .

2кф1 ^

7г2 + г2 —Л

2кф1

+1

2кр1

+ г2-.

Таким образом, если г мало, то с учётом функций К„; /„;Ь„ (8) можно приближённо записать в следующем виде:

Е(„, „, г) »1 - ^ехр ['кЯ(1 - р)] -[ -

-'крг

2

1

2 + р

р

+ 'крг ехр

х{„, 571' + [„,647 + „,331п (кр1г)]}.

Теоретически при очень малом г слагаемое, содержащее логарифм, должно давать заметный вклад.

„

х

е

„

„

2

я

2

я

Попытаемся найти эту границу. Для этого проведём некоторые упрощения. Если Я~ 1 / (3яр), то в (13) можно опустить второе слагаемое; также учтём, что при р третье слагаемое неограниченно растёт, поэтому нас интересует область небольших р, когда в третьем слагаемом выражение в скобках можно заменить на 2р. В итоге получим:

Е(0,0, ¿) »1 - 2йр2г + ikpzехр (кЬ0) х х{0,57Ь' + [0,647 + 0,331п (к01 г)]}.

Анализируя выражение (14), можно сказать с учётом его приближённости, что слагаемое с логарифмом начинает сказываться (но ещё не преобладать) только на довольно малых расстояниях г < 0,11. В то же время поправочное слагаемое (9) при малом г имеет величину, близкую к основному слагаемому. Поясним это на примере. Значение р1 близко к единице (величина в фигурных скобках также близка к единице по модулю): при г = е1/(2я) ~ 0,431 оно равно 1,131, а при г = 1/(2я) ~ 0,1591 - равно 0,863. Получаем, что при р < 1/2 поправочное слагаемое больше основного (если основное взять без единицы).

2.2. Числовая апертура больше или равна предельной (р > 1) Очевидно, что вычисления при р > 1 существенно отличаются от описанных ранее, т.к. использование метода стационарной фазы невозможно в силу отсутствия стационарной точки. При использовании классического метода стационарной фазы пришлось бы сказать, что здесь на любом расстоянии г будет тень, хотя и не совсем ноль, т.к. в (7) имеются и внеинтегральные слагаемые. В геометро-оптическом подходе такой ак-сикон считается не пропускающим излучение.

В силу заметного отличия вида фазовой функции при р = 1 и р > 1 рассмотрим эти случаи отдельно. Хотя первый случай имеет меру нуль, он является хорошим приближением для ситуаций, когда р близко к единице.

2.2.1. Числовая апертура равна предельной (р = 1) Для р = 1 формула (7) принимает следующий вид:

гй , „2

Е(0,0, г) = ехр^) -

г • ехр

я ехр (ik7777 - ikг) -ikz|-'—, -'- &г.

7777

Чтобы вычислить интеграл в (15), удобно разбить отрезок интегрирования на две части и на каждой сделать своё приближение фазовой функции. Деление будем производить в точке максимальной кривизны фазовой функции.

Можно доказать, что поиск точки максимальной кривизны эквивалентен поиску точки графика функции г - 7777, ближайшей к началу координат.

Квадрат расстояния равен 3г2 + г2 - 2г7777 . Минимум этого выражения достигается при значении г,

удовлетворяющем уравнению 5г4 + 5 Л? - г4 = 0, решение которого имеет следующий вид:

2 г2 I 3 , . Г2 =-1 -;= - 1 I ^ Г

2

0,413г.

Если гс > Я, что выполняется при не очень малых г, то второго отрезка не будет. Для общности обозначим правый край первого отрезка гр = тш(гс, Я).

На первом отрезке функцию г - 7777 приблизим параболой Ь212 + Ь1г + Ь> Для нахождения коэффициентов используем коллокацию: приравняем значения функций на обоих концах отрезка и значения производных на левом конце. В результате получим:

Ь0 = -^ А = 1; Ь2 =

Если второй отрезок есть (гр = гс), то Ь2 » -0,48 / г. Подставив параболическое приближение в интеграл из (15), получим

гр ехр

1  -1

7г2 + г2 гр ехр I -ik (Ь2г2 + г - г) I

') -- аг =

7г2 + г2

= ехр

/к

г + -

2Ь

ехр

-кЬ2 (г- г0 )2

7г2 + г2

аг,

где г0 = -1/ (2Ь2) = 1,042 г.

Далее возможно применение метода стационарной фазы, и тогда вместо (18):

гр ехр (-ikг + ik^/777) | 1 , -^

7777

ехр

1к

г +-

2Ь

7-Ь^7 г2 + г2

'2

| ехр (^) ^,

где 1 = ^2 = (грг - г0 )4~Ь2 .

Искусственно полученная стационарная точка г0 лежит за пределами отрезка интегрирования, поэтому интеграл (19) будет иметь малую величину.

На втором отрезке, если он есть, функцию

г- 777 г2 приблизим дробью -£3 / (г + Ь4). Для нахождения коэффициентов можно использовать различные типы коллокации. Если приравнять значения функций на левом конце и асимптотическое поведение при большом г, то получится простое выражение для Ь3 = ^/2.

Метод стационарной фазы здесь неприменим, но можно использовать некий его аналог, пользуясь тем,

что функция 1/7г2 + г2 на этом участке меняется сравнительно медленно и её можно заменить значением в срединной точке. Тогда:

0

0

р /

1

0

0

0

,exp

ikb3 r + b

7777

dr =

VTf+z'

J exp

ik

b,

r + b

dr =

4 J

i

r2 + z2

exp

ikb3 R + b

(R + b.) -

4J

exp

ikb, V f + b4

'2

(r, + b.) + ikb, J

äft,x

dx ^

где х = У( ^ + Ь), х = У(гр + Ь).

Оставшийся в (20) интеграл тривиально выражается через интегральные синус и косинус.

Исходя из того, что функция г-7г2 + г2 при увеличении г стремится к нулю, что соответствует исчезновению осцилляций, можно сделать вывод, что, увеличивая R, мы будем иметь всё возрастающую амплитуду. К сожалению, приближение в (20) не показывает этого - амплитуда стремится к конечному постоянному значению.

Более грубую оценку можно получить значительно проще, заменяя фазу нулём:

exp (ikb,/ (г + b.)) sjr2 + z2

dr

dr

7777

= ln (R + 7 R2 + z2)- ln (rp +J p + z2).

Оценка (21) будет более точной, если нижний предел заменить на такое значение, которое соответствует началу последнего лепестка, но на зависимость от R это влияния не оказывает. Как видно из (21), рост амплитуды намного слабее, чем в случае параболической линзы в приближении Френеля [30], где он был пропорционален R2, но там сама величина R была ограничена в силу требования параксиальности.

2.2.2. Числовая апертура больше предельной р>1

Выражение (7) в этом случае принимает следующий вид:

E (0,0, z) = exp(ikz) -

>( ik4

z • exp l

7R2 + z2

-ikpz • J

R exp ik

r2 + z2 -

ikpr)

7r2 + z2

dr.

Заметим, что фазовая функция рг - 777 г2 в интеграле (22) при р > 1 имеет совсем иную зависимость, чем было при р = 1. После деления отрезка интегрирования на две части можно выполнить аппроксимацию фазовой функции линейными функциями вида Ь1г + Ь0. Деление производится так же, как при

р = 1: в точке графика функции рг - г2 , ближайшей к началу координат, хотя визуально она не очень выражена. Как и при р = 1, она совпадает с точкой максимальной кривизны.

Найдём эту точку: квадрат расстояния от начала координат равен (2 + р2)г2 + г2 - 2 рг7г2 + г2 . Минимум достигается при значении г , удовлетворяю-

4 4 4 2 2 2 4

щем уравнению (4+ р ) г + (4+ р ) г г - рг =0, решение которого имеет следующий вид:

z2 2

2 + p2

f4

г- 1

+ p4

При увеличении p выражение в квадратных скобках сначала возрастает до 72 -1 при p = 72, а затем асимптотически убывает до нуля.

Так же, как в 1.2.1, обозначим правый край первого отрезка rp = min (rc, R). Подставив в интеграл из (22) приближение фазовой функции, получим:

p J

exp l

7777

dr =

t -1U exp(-ikb1r)л - exp(-ikb0) J —. dr.

о 7r

: + z2

Если гр хотя бы в 2-3 раза превышает я/(Щ), то

верхний предел можно заменить бесконечным, что позволит использовать справочные выражения. Учитывая, что величина Ь1 близка к р (обычно немного меньше), и используя (23), при среднем коэффициенте 2,5 получим, что замена возможна, если

2,572

z >

2 nl (2 + pVTiTp4 -1

•l.

В частном случае при р = 1,2 имеем г > 2,341, а при увеличении р в пределе получим г> 1,251. Разделив вещественную и мнимую части, после использования табличных формул [29], получим ответ:

¥ ехр

7777

-dr =

= K0 (kb z) - i - [ I0 (kb z) - L0 (kb z)].

Вычисление интеграла по второму отрезку можно выполнить аналогично (заменить предел на бесконечный и вычесть интеграл по первому отрезку от функции для второго отрезка) только при достаточно большом радиусе R. Это значение заметно больше, чем требуется для замены верхнего предела на бесконечный. Дело в том, что первый отрезок не вмещает приблизительно целое число периодов с параметрами второго отрезка. Поэтому в левой части формулы, аналогичной (25), под корнем будет не г2, а (г - г50г)2 , а такой формулы нет в справочнике. При большом же радиусе этот эффект можно игнорировать.

Полученные формулы показывают, что при больших г амплитуда убывает экспоненциально, и это не противоречит априорно известному убыванию амплитуды как г4, так как такое убывание обеспечивают первые два слагаемых в (8).

1

J

r

r

p

p

1

r2 =

Xi

0

R

R

J

J

<

r

p

p

0

0

Если применить формулу (25) невозможно, то всё же есть способ вычислить интеграл, фигурирующий в (22). Он требует больше вычислений, чем обычно требуют приближённо-аналитические методы; но при небольших г и Я проигрыш невелик. Его отличие в том, что фазовая функция в (22) на начальном шаге сохраняется в исходном виде без аппроксимации. Опишем этот способ.

Шаг 1.

Преобразуем экспоненту в (22) следующим образом:

•>(iks[r

exp I ik^jr2 + z2 -

= eikz exp I -ik

ikpr) = pr + z -V r2 + z2

= eikz (cos(kf2 (r)) - i sin(kf2 (r))).

Такое представление удобно тем, что даёт одинаковую фазу /2 (г = 0) = 0 независимо от г, и тем, что появляется такой же фазовый множитель в'кг, как у первого слагаемого в (22); второе удобство проявится ниже. После этого находим границы лепестков (нули) у косинуса и синуса.

Нули косинуса:

rc =

m

V( Р2 - !) z2 + L - pzm ,

p2 -1

zm = z-p(m + 0,5)/k; m = 0,1,2,....

Нули синуса т вычисляются по той же формуле,

но меняется гт = г - рт / к; т = 1,2,.....

Если подставить значения гт и обозначить Р1, то получим «безразмерные» границы лепестков:

гС

V

p2p2 - b( m + 0,5) + (m + 0,5) - pb + 0,5 p(m + 0,5) 4

2 , -'

Р - 1

т 7р2Р2-Рт+ т /4 - рР + 0,5рт

Т= рг-1 .

Шаг 2.

После того как границы лепестков найдены, первый лепесток аппроксимируется параболой. А именно: первый лепесток косинуса (от 0 до г^, строго говоря, это полулепесток) приближаем кривой 1 - (г / г^ )2. При малом г первый лепесток синуса (от 0 до г]5) заметно несимметричен относительно максимума, поэтому левая и правая половины лепестка аппроксимируются отдельно: левый кривой 1 - ((г- г^)/ гЦ)2, а правый кривой 1 - ((г - г0с )/(г^ - г0с ))2. При не очень

малом г первый лепесток синуса становится примерно симметричен, поэтому его можно приблизить единой кривой 1 - ((г - 0,5,1 )/0,5г;5 )2. Остальные лепестки приближаются лепестками обычной симмет-

ричной синусоиды. Поэтому площадь всех лепестков, кроме первого, приближённо вычисляется следующим образом:

^ = (-1Г1 J

! rn cos

Vr2 + z2

dr-

(-1)n

J cos

>/((r-1 + П )/2)2 + z2

= 2(-1)n-1 _ r^-r^

(rD_1 + rn )/2)2 + z2 p .

Если последний лепесток не является приближённо полным, то в данном выражении множитель 2 заменится на 1 - cos 8, где угол 8 получаем из пропорции: rn - rn-1 соответствуют 180 градусам, а R- rn-1 соответствуют 8 . Если z мало, то можно положить z = 0 .

Шаг 3.

При вычислении площади первого лепестка знаменатель сохраняется в исходном виде и используются формулы для интегралов

dr

Vr2 + z2' J Vr2 + z2

J-

r dr

r dr

1 Wr2 + z2 - — in (r+V r2 + z2).

2 2 v I

__- гу г2 + г2

у1г2 + г2 2 2

Если г мало, то можно использовать следующие приближения. Первый лепесток косинуса имеет площадь 1п(2г,с / г) - 0,5, левая половина лепестка синуса имеет площадь 1,5 - (2 г / г^).

3. Анализ поперечного распределения

В данном разделе рассмотрим поперечное распределение комплексной амплитуды. Для этого перепишем (6) в полярных координатах:

Е (р, 8, г) =

= JJ E0(r cos j, r sin j)eXp(ikl) | ik -1J rdrdj

exp

I2 I I.

где I2 = р2 + г2 + г2 - 2ргсоБ(ф- 8). Пусть входное поле радиально-симметрично, тогда полярный угол в подынтегральном выражении будет только под знаком косинуса; поэтому интеграл по полному периоду не зависит от сдвига угла 8 . То есть выходное поле будет тоже радиально-симметрично и при вычислениях можно положить 8 = 0 :

E (р, z) = - 2Р J E0(r)rdr J

I2

exp

I

' 0 0

Ограничимся распределением вблизи оптической оси, что позволит сделать упрощения: в выражении для I положим р = 0, кроме показателя экспоненты, где разложим I по степеням р , ограничившись дву-

мя слагаемыми

(I «>/r2 + z2-prcos j/Vr2 + z2 ). В

)

l

S

этом случае интегрирование по полярному углу производится аналитически, а конечный результат будет следующим:

¥ ( \

1к -

4

Е

0 V х exp (1кл1г2 + г2) крГ

1

г2 + г2 ;

2 2 ' г + г

7777

с1г.

Для получения хотя бы начального приближения используем классический метод стационарной фазы. Тогда зависимость от поперечной координаты будет только в функции Бесселя, что позволит легко найти радиусы нескольких первых колец.

В случае аксикона Е0 (г) = exp(-ikpг), ограничиваясь случаем апертуры меньше предельной, находим стационарную точку г5 = рг/-^1 - р2 .

В результате получаем, что Е(р,г) ~ ]0(крр), т.е. получаем пучок Бесселя, масштаб которого зависит от числовой апертуры.

4. Численный расчёт на основе интеграла Рэлея-Зоммерфельда

В данном разделе приведены результаты моделирования дифракции ограниченной плоской волны на аксиконе (3) с помощью численного интегрирования по формуле (6). Параметры расчёта: длина волны 1 = 0,532 мкм, радиус ограничивающей апертуры Я = 2,5 мкм.

На рис. 1 показано распределение интенсивности вдоль оптической оси для аксиконов с различными числовыми апертурами.

90

50

10

а) 0

Л\ : / / ^ / \

/ г / / 1 \ \ \ \ \ ч\ ч,

.......... .....

4 2, МКМ

0 1 2 3 4 2, мкм

Рис. 1. Осевое распределение интенсивности для аксиконов с числовыми апертурами ЫА=0,7

Как видно из рис. 1а, с увеличением числовой апертуры до предельной ЫА = 1 происходит постепенное уменьшение концентрации энергии на оптической оси. Это значит, что энергия уходит во внеосевую область. В случае, когда числовая апертура становится больше предельной (ЫА = 1,2) (рис. 1б), энергия практически перестаёт проходить через оптический элемент, однако явно становится заметным вклад затухающих волн на расстоянии менее длины волны.

Вид графика на рис. 1б можно, по крайней мере качественно, объяснить при помощи рассуждений из конца предшествующего пункта. Если в формуле (22) опустить второе слагаемое, подставить выражение интеграла через площади лепестков, положить р = 1,2 и г = Р1, то без учёта фазового множителя е1кг получим выражение:

Е

Производя вычисление площадей лепестков, как указано выше, можно найти местоположение всех экстремумов. Даже не делая слишком частый перебор значений Р, можно обнаружить экстремумы, хотя и не так точно, как при прямом вычислении по квадратурным формулам. Но в отличие от квадратур, можно произвести доказательство наличия экстремумов.

При очень малом Р первый лепесток превалирует над всеми остальными, причём его площадь с ростом Р убывает (это очевидно геометрически и следует из приближений для площади первых лепестков). Мы получаем произведение Р на убывающую функцию, что означает наличие максимума. Вычисляя площади всех лепестков, получим значение максимума 1,758, которое соответствует z = 0,121. Рассчитанное на основе численного интегрирования значение получилось 0,241, что в 2 раза больше приведённой оценки.

С увеличением Р площади всех лепестков, очевидно, убывают, но убывание тем быстрее, чем меньше номер лепестка. Поэтому настанут моменты, когда площади обратятся в нуль. Нуль площади косинусов даёт нуль мнимой части амплитуды, нуль вещественной части будет немного раньше, чем обратится в нуль площадь синусов (см. формулу (33)). В этой области получится минимум. Он получается в области 0,351 < г < 0,391. Далее площади становятся отрицательными, но их модуль на некотором участке будет расти, поэтому начнёт расти и амплитуда.

В то же время из (25) и наглядных соображений следует, что при больших г амплитуда монотонно убывает к нулю, поэтому неизбежно должен быть и второй максимум. Отметим, что в области применимости (25) опускать второе слагаемое в (22) нужно с осторожностью.

Таким образом, методика, предложенная в разделе 2 для аксикона с числовой апертурой выше предельной, позволяет качественно предсказывать поведение интенсивности вдоль оптической оси.

Структура внеосевого распределения показана в табл. 1.

Таблица 1. Результаты моделирования дифракции ограниченной плоской волны на аксиконе (3) _с помощью интеграла Рэлея-Зоммерфельда (6)_

р=0,7 р=0,9 р = 1 р = 1,2

Вид фазы аксикона (3) # ф ф #

Распределение амплитуды в продольной плоскости z е [0,1 мкм, 5 мкм], х е [-3 мкм, 3 мкм] 1 1

Распределение интенсивности в поперечной плоскости с максимальной интенсивностью • • • •

х, уе [-1 мкм, 1 мкм] z = 1,3 мкм FWHM = 0,26 мкм р0 = 0,27 мкм z = 0,6 мкм FWHM = 0,21 мкм р0 = 0,22 мкм z = 0,37 мкм FWHM = 0,2 мкм р0 = 0,21 мкм z = 1,1 мкм FWHM = 0,21 мкм р0 = 0,22 мкм

Также в табл. 1 показаны распределения интенсивности в поперечной плоскости с максимальной интенсивностью и приведены полученные численно размеры центрального пятна. Видно, что с увеличением числовой апертуры размер центрального светового пятна уменьшается, вплоть до достижения предельного значения FWHM = 0,2 мкм = 0,371.

Оценки, полученные в разделе 3, предсказывают

2,4 0,203

радиус центрального пятна р0 =-=-мкм , что

кр р

согласуется с результатами численных расчётов вплоть до р=1. При р > 1 полученные теоретические оценки следует применять с осторожностью, т.к. при их выводе предполагалось р<1 (см. конец пункта 3).

Дальнейшее увеличение числовой апертуры не приводит к уменьшению размера светового пятна. Аналогичный результат был получен в работе [26].

Заметим, что использование скалярной волновой модели в ближней зоне дифракции правомерно только для одной из поперечных компонент электрического поля. Причём с увеличением числовой апертуры роль продольной компоненты электрического поля становится очень важной, её вклад может превышать вклад поперечных компонент [31], хотя для однородно-поляризованного излучения, не имеющего фазовой сингулярности, энергия продольной компоненты электрического поля концентрируется вне оптической оси.

Известны ситуации [32-34], когда вещество или устройство является чувствительным только к поперечным или продольной компонентам электрического

поля. Таким образом, полученные аналитические оценки становятся актуальны не только для расчёта отдельных компонент, но и картины в целом.

Заключение

В работе выполнен теоретический анализ дифракции плоской волны на аксиконе в зоне затухающих волн.

Для распределения комплексной амплитуды вдоль оптической оси получены аналитические оценки для высоких значений числовой апертуры, в том числе превышающих предельное.

При числовой апертуре меньше предельной использован модифицированный метод стационарной фазы с поправками, которые существенны только в ближней зоне дифракции. Поправки выражены аналитически через модифицированные функции Бесселя.

При числовой апертуре, равной предельной или больше неё, использование метода стационарной фазы невозможно. Для анализа в этом случае использована аппроксимация подынтегральной функции, обеспечивающая сохранение точных нулей вещественной и мнимой частей. Такая модель позволяет качественно предсказывать поведение интенсивности вдоль оптической оси.

Оценки, полученные для поперечного распределения, хорошо согласуются с результатами численных расчётов вплоть до предельной числовой апертуры.

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

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

Работа выполнена при финансовой поддержке

гранта РФФИ 13-07-97004-р_поволжье_а.

Литература

1. McLeod, J.H. The axicon: a new type of optical element / J.H. McLeod // Journal of the Optical Society of America. -1954. - V. 44. - P. 592-597.

2. Durnin, J. Diffraction-free beams / J. Durnin, J.J. Miceli, Jr. and J.H. Eberly // Physical Review Letters. - 1987. -V. 58. - P. 1499-1501.

3. McGloin, D. Bessel beams: diffraction in a new light / D. McGloin and K. Dholakia // Contemporary Physics. -2005. - V. 46(1). - P. 15-28.

4. Jaroszewicz, Z. Axicon - the most important optical element / Z. Jaroszewicz, A. Burvall, A.T. Friberg // Optics & Photonics News, April 2005.

5. Arlt, J. Optical micromanipulation using a Bessel light beam / J. Arlt, V. Garces-Chavez, W. Sibbett and K. Dholakia // Optics Communications. - 2001. - V. 197. - P. 239-245.

6. Garces-Chavez, V. Simultaneous micromanipulation in multiple planes using a self-reconstructing light beam / V. Garces-Chavez, D. McGloin, H. Melville, W. Sibbett and K. Dholakia // Nature. - 2002. - V. 419. - P. 145-147.

7. Сойфер, В.А. Оптическое манипулирование микрообъектами: достижения и новые возможности, порождённые дифракционной оптикой / В.А. Сойфер, В.В. Кот-ляр, С.Н. Хонина // Физика элементарных частиц и атомного ядра. - 2004. - Т. 35(6). - P. 1368-1432.

8. Ding, Z. High-resolution optical coherence tomography over a large depth range with an axicon lens / Z. Ding, H. Ren, Y. Zhao, J.S. Nelson and Z. Chen // Optics Letters. - 2002. - V. 27. - P. 243-245.

9. Lee, K. Bessel beam spectral-domain high-resolution optical coherence tomography with micro-optic axicon providing extended focusing range / K. Lee and J. Rolland // Optics Letters. - 2008. - V. 33. - P. 1696-1698.

10. Котляр, В.В. Бесконтактное прецизионное измерение линейных смещений с использованием ДОЭ, формирующих моды Бесселя / В.В. Котляр, Р.В. Скиданов, С.Н. Хонина // Компьютерная оптика. - 2001. - № 21. - C. 102-104.

11. Fortin, M. Optical tests with Bessel beam interferometry / M. Fortin, M. Piché and E.F. Borra // Optics Express. -2004. - V. 2, issue 24. - P. 5887-5895.

12. Turunen, J. Holographic generation of diffraction-free beams / J. Turunen, A. Vasara and A.T. Friberg // Applied Optics. - 1988. - V. 27. - P. 3959-3962.

13. Arlt, J. Optical dipole traps and atomic waveguides based on Bessel light beams / J. Arlt [et al.] // Physical Review. -2001. - V. 63. - P. 063602.

14. Cizmar, T. An optical nanotrap array movable over a mili-metre range / T. Cizmar, M. Siler, P. Zemanek // Applied Physics B. - 2006. - V. 84. - P. 197-203.

15. Khonina, S.N. Narrowing of a light spot at diffraction of linearly-polarized beam on binary asymmetric axicons / S.N. Khonina, D.V. Nesterenko, A.A.Morozov, R.V. Skidanov, V.A. Soifer // Optical Memory and Neural Networks (Information Optics). - 2012. - V. 21(1). - P. 17-26.

16. Khonina, S.N. Experimental demonstration of the generation of the longitudinal E-field component on the optical axis with high-numerical-aperture binary axicons illuminated by linearly and circularly polarized beams / S.N. Khonina, S.V. Karpeev, S.V. Alferov, D.A. Savelyev, J. Laukkanen, J. Turunen // Journal of Optics. - 2013. -V. 15. - P. 085704 (9 p).

17. Vahimaa, P. Electromagnetic analysis of nonparaxial Bes-sel beams generated by diffractive axicons / Pasi Vahimaa,

Ville Kettunen, Markku Kuittinen, Jari Turunen and Ari T. Friberg // Journal of the Optical Society of America A. -1997. - V. 14(8). - P. 1817-1824.

18. Zhang, Y. Vector propagation of radially polarized Gaussian beams diffracted by an axicon / Y. Zhang, L. Wang, C. Zheng // Journal of the Optical Society of America A. -2005. - V. 22(11). - P. 2542-2546.

19. Котляр, В.В. Моделирование острой фокусировки ра-диально-поляризованной лазерной моды с помощью конического и бинарного микроаксиконов /

B.В. Котляр, С.С. Стафеев // Компьютерная оптика. -2009. - Т. 33, № 1. - С. 52-60.

20. Хонина, С.Н. Дифракция на бинарных микроаксиконах в ближней зоне / С.Н. Хонина, Д.А. Савельев, П.Г. Серафимович, И.А. Пустовой // Оптический журнал. -2012. - Т. 79, № 10. - С. 22-29.

21. Totzeck, M. Validity of the scalar Kirchhoff and Rayleigh-Sommerfeld diffraction theories in the near field of small phase objects // Journal of the Optical Society of America A. - 1991. - V. 8(1). - P. 27-32.

22. Tsoy, V.I. The use of Kirchhoff approach for the calculation of the near field amplitudes of electromagnetic field / V.I. Tsoy, L.A. Melnikov // Optics Communications. -

2005. - V. 256. - P. 1-9.

23. Устинов, А.В. Обобщённая линза: анализ осевого и поперечного распределения / А.В. Устинов, С.Н. Хонина // Компьютерная оптика. - 2013. - Т. 37, № 3. - С. 307-315.

24. Хонина, С.Н. Распространение радиально-ограниченных вихревых пучков в ближней зоне. Часть I. Алгоритмы расчёта / С.Н. Хонина, А.В. Устинов, А.А. Ковалёв,

C.Г. Волотовский // Компьютерная оптика. - 2010. -Т. 34, № 3. - С. 317-332.

25. Хонина, С.Н. Дифракция на бинарных микроаксиконах в ближней зоне / С.Н. Хонина, Д.А. Савельев, П.Г. Серафимович, И.А. Пустовой // Оптический журнал. -2012. - Т. 79, № 10. - С. 22-29.

26. Устинов, А.В. Расчёт комплексной функции пропускания рефракционных аксиконов / А.В. Устинов, С.Н. Хонина // Компьютерная оптика. - 2011. - Т. 35, № 4. - С. 480-490.

27. Born, M. Principles of Optics / M. Born, E. Wolf - 6th ed.

- Oxford: Pergamon, 1980. - Chap. 8.3.

28. Устинов, А.В. Анализ дифракции плоского пучка на рассеивающем фраксиконе в непараксиальном режиме / А.В. Устинов, С.Н. Хонина // Компьютерная оптика. -2014. - Т. 38, № 1. - С. 42-50.

29. Диткин, В.А. Интегральные преобразования и операционное исчисление / В.А. Диткин, А.П. Прудников. -М.: Физматлит, 1961. - Главы VII, VIII.

30. Устинов, А.В. Сравнительный анализ параболической линзы и аксикона в моделях геометрической и скалярной параксиальной оптики / А.В. Устинов, А.В. Карсаков, С.Н. Хонина // Вестник Самарского государственного аэрокосмического университета имени академика С.П. Королёва (национальный исследовательский университет). - 2012. - № 4(35). - C. 230-239.

31. Khonina, S.N. Controlling the contribution of the electric field components to the focus of a high-aperture lens using binary phase structures / S.N. Khonina, S.G. Volotovsky // Journal of the Optical Society of America A. - 2010. -V. 27(10). - P. 2188-2197.

32. Huse, N. Z-polarized confocal microscopy / N. Huse, A. Schonle, and S. W. Hell // Journal of Biomedical Optics.

- 2001. - V. 6. - P. 273-276.

33. Grosjean, T. Photopolymers as vectorial sensors of the electric field / T. Grosjean and D. Courjon // Opt. Express. -

2006. - V. 14. - P. 2203-2210.

34. Dedecker, P. Orientational effects in the excitation and de-excitation of single molecules interacting with donut-mode

laser beams / P. Dedecker, B. Muls, J. Hofkens, J. Enderlein, and J.-I. Hotta // Optics Express - 2007. - V. 15. -P. 3372-3383.

References

1. McLeod, J.H. The axicon: a new type of optical element / J.H. McLeod // Journal of the Optical Society of America. -1954. - V. 44. - P. 592-597.

2. Durnin, J. Diffraction-free beams / J. Durnin, J.J. Miceli, Jr. and J.H. Eberly // Physical Review Letters. - 1987. -V. 58. - P. 1499-1501.

3. McGloin, D. Bessel beams: diffraction in a new light / D. McGloin and K. Dholakia // Contemporary Physics. -2005. - V. 46(1). - P. 15-28.

4. Jaroszewicz, Z. Axicon - the most important optical element / Z. Jaroszewicz, A. Burvall, A.T. Friberg // Optics & Photonics News, April 2005.

5. Arlt, J. Optical micromanipulation using a Bessel light beam / J. Arlt, V. Garces-Chavez, W. Sibbett and K. Dholakia // Optics Communications. - 2001. - V. 197. - P. 239-245.

6. Garces-Chavez, V. Simultaneous micromanipulation in multiple planes using a self-reconstructing light beam / V. Garces-Chavez, D. McGloin, H. Melville, W. Sibbett and K. Dholakia // Nature. - 2002. - V. 419. - P. 145-147.

7. Soifer, V.A. Optical microparticle manipulation: advances and new possibilities created by diffractive optics / V.A. Soifer, V.V. Kotlyar, S.N. Khonina // Physics of Particles and Nuclei. - 2004. - Vol. 35 (6), P. 733-766.

8. Ding, Z. High-resolution optical coherence tomography over a large depth range with an axicon lens / Z. Ding, H. Ren, Y. Zhao, J.S. Nelson and Z. Chen // Optics Letters. - 2002. -V. 27. - P. 243-245.

9. Lee, K. Bessel beam spectral-domain high-resolution optical coherence tomography with micro-optic axicon providing extended focusing range / K. Lee and J. Rolland // Optics Letters. - 2008. - V. 33. - P. 1696-1698.

10. Kotlyar, V.V. Contactless precision measurement of a linear displacement using DOEs forming Bessel beams / V.V. Kotlyar, R.V. Skidanov and S.N. Khonina // Computer Optics. - 2001. - V. 21. - P. 102-104. - (In Russian).

11. Fortin, M. Optical tests with Bessel beam interferometry / M. Fortin, M. Piché and E.F. Borra // Optics Express. -2004. - V. 2, Issue 24. - P. 5887-5895.

12. Turunen, J. Holographic generation of diffraction-free beams / J. Turunen, A. Vasara and A.T. Friberg // Applied Optics. - 1988. - V. 27. - P. 3959-3962.

13. Arlt, J. Optical dipole traps and atomic waveguides based on Bessel light beams / J. Arlt [et al.] // Physical Review. -2001. - V. 63. - P. 063602.

14. Cizmar, T. An optical nanotrap array movable over a mili-metre range / T. Cizmar, M. Siler, P. Zemanek // Applied Physics B. - 2006. - V. 84. - P. 197-203.

15. Khonina, S.N. Narrowing of a light spot at diffraction of linearly-polarized beam on binary asymmetric axicons / S.N. Khonina, D.V. Nesterenko, A.A.Morozov, R.V. Skidanov, V.A. Soifer // Optical Memory and Neural Networks (Information Optics). -2012. - V. 21(1). - P. 17-26.

16. Khonina, S.N. Experimental demonstration of the generation of the longitudinal E-field component on the optical axis with high-numerical-aperture binary axicons illuminated by linearly and circularly polarized beams / S.N. Khonina, S.V. Karpeev, S.V. Alferov, D.A. Savelyev, J. Laukkanen, J. Turunen // Journal of Optics. - 2013. - V. 15. - P. 085704 (9 p).

17. Vahimaa, P. Electromagnetic analysis of nonparaxial Bes-sel beams generated by diffractive axicons / Pasi Vahimaa, Ville Kettunen, Markku Kuittinen, Jari Turunen and Ari

T. Friberg // Journal of the Optical Society of America A. -1997. - V. 14(8). - P. 1817-1824.

18. Zhang, Y. Vector propagation of radially polarized Gaussian beams diffracted by an axicon / Y. Zhang, L. Wang,

C. Zheng // Journal of the Optical Society of America A. -2005. - V. 22(11). - P. 2542-2546.

19. Kotlyar, V.V. Modeling sharp focus radially-polarized laser mode with conical and binary microaxicons / V.V. Kotlyar, S.S. Stafeev // Computer Optics. - 2009. - V. 33(1). -P. 52-60. - (In Russian).

20. Khonina, S.N. Diffraction at binary microaxicons in the near field / S.N. Khonina, P.G. Serafimovich,

D.A. Savelyev, I.A. Pustovoi // Journal of Optical Technology. - 2012. - V. 79(10). - P. 626-631.

21. Totzeck, M. Validity of the scalar Kirchhoff and Rayleigh-Sommerfeld diffraction theories in the near field of small phase objects // Journal of the Optical Society of America A. - 1991. - V. 8(1). - P. 27-32.

22. Tsoy, V.I. The use of Kirchhoff approach for the calculation of the near field amplitudes of electromagnetic field / V.I. Tsoy, L.A. Melnikov // Optics Communications. -2005. - V. 256. - P. 1-9.

23. Ustinov, A.V. Generalized lens: Calculation of distribution on the optical axis / A.V. Ustinov, S.N. Khonina // Computer Optics. - 2013. - V. 37(3). - P. 307-315.

24. Khonina, S.N. Propagation of the radially-limited vortical beam in a near zone. Part I. Calculation algorithms / S.N. Khonina, A.V. Ustinov, A.A. Kovalev, S.G. Volotovsky // Computer Optics. - 2010. - V. 34, No. 3. - P. 317-332. - (In Russian).

25. Khonina, S.N. Diffraction at binary microaxicons in the near field / S.N. Khonina, P.G. Serafimovich, D.A. Savelyev, I.A. Pustovoi // Journal of Optical Technology. - 2012. - T. 79, № 10. - P. 626-631.

26. Ustinov, A.V. Calculating the complex transmission function of refractive axicons / A.V. Ustinov and S.N. Khonina // Optical Memory and Neural Networks (Information Optics). - 2012. - V. 21(3). - P. 133-144.

27. Born, M. Principles of Optics / M. Born, E. Wolf - 6th ed. - Oxford: Pergamon, 1980. - Chap. 8.3.

28. Ustinov, A.V. Analysis of flat beam diffraction by divergent fracxicon in nonparaxial mode / A.V. Ustinov, S.N. Khonina, // Computer Optics. - 2014. - V. 38(1). - P. 42-50.

29. Ditkin, V.A. Integral transforms and operational calculus / V.A. Ditkin and A.P. Prudnikov. - Moscow: "Fizmatlit" Publisher, 1961. - Charters VII, VIII. - (In Russian).

30. Ustinov, A.V. Comparative analysis of parabolic lens and axicon in geometric and scalar paraxial optical models / A.V. Ustinov, A.V. Karsakov, S.N. Khonina // Bulletin SSAU. - 2012. - № 4(35). - C. 230-239.

31. Khonina, S.N. Controlling the contribution of the electric field components to the focus of a high-aperture lens using binary phase structures / S.N. Khonina, S.G. Volotovsky // Journal of the Optical Society of America A. - 2010. -V. 27(10). - P. 2188-2197.

32. Huse, N. Z-polarized confocal microscopy / N. Huse, A. Schonle, S.W. Hell // Journal of Biomedical Optics. -2001. - V. 6. - P. 273-276.

33. Grosjean, T. Photopolymers as vectorial sensors of the electric field / T. Grosjean and D. Courjon // Optics Express. - 2006. - V. 14. - P. 2203-2210.

34. Dedecker, P. Orientational effects in the excitation and de-excitation of single molecules interacting with donut-mode laser beams / P. Dedecker, B. Muls, J. Hofkens, J. Enderlein, and J.-I. Hotta // Optics Express - 2007. - V. 15. -P. 3372-3383.

ANALYSIS OF LASER BEAM DIFFRACTION BY AXICON WITH THE NUMERICAL APERTURE ABOVE LIMITING

A. V. Ustinov, S.N. Khonina Image Processing Systems Institute, Russian Academy of Sciences, Samara State Aerospace University

Abstract

The theoretical analysis of plane wave diffraction by axicon in the zone of evanescent waves is executed. The received analytical estimations are qualitatively in accordance with results of numerical calculations in nonparaxial scalar wave model.

Key words: axicon, limiting numerical aperture, evanescent waves.

Сведения об авторах Устинов Андрей Владимирович, 1968 года рождения, в 1991 году окончил Куйбышевский авиационный институт имени академика С.П. Королёва (КуАИ) по специальности «Прикладная математика», работает ведущим программистом в Федеральном государственном бюджетном учреждении науки Институте систем обработки изображений РАН; является аспирантом Самарского государственного аэрокосмического университета имени академика С.П. Королёва (национальный исследовательский университет). Область научных интересов: дифракционная оптика, разработка программ моделирования работы оптических элементов; обработка изображений, в частности, гидродинамических процессов и биомедицинских. E-mail: andr@smr.ru .

Andrey Vladimirovich Ustinov, (b. 1968) graduated from Kuibyshev Aviation Institute named after academician S.P. Korolyov (KuAI), majoring in Applied Mathematics, works as the leading programmer in the Image Processing Systems Institute of the RAS; postgraduate student of Samara State Aerospace University named after S.P. Korolyov. Research interests: diffractive optics, software design for modeling of optical elements operating; images processing, particularly images of hydrodynamic processes and biomedical images.

Сведения об авторе Хонина Светлана Николаевна - см. стр. 180этого номера.

Поступила в редакцию 31 марта 2014 г.