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

УДК 535.42

РАСЧЁТ ДИФРАКЦИИ ПЛОСКОЙ ВОЛНЫ НА РАССЕИВАЮЩЕМ ДРОБНОМ АКСИКОНЕ С УЧЁТОМ ЗАТУХАЮЩИХ ВОЛН

© 2014 А.В. Устинов, С.Н. Хонина Институт систем обработки изображений Российской академии наук, г. Самара

Поступила в редакцию 20.06.2014

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

ВВЕДЕНИЕ

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

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

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

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

Устинов Андрей Владимирович, ведущий программист. E-mail: andr@smr.ru

Хонина Светлана Николаевна, доктор физико-математических наук, ведущий научный сотрудник. E-mail: khonina@smr.ru

гарифмическому аксикону [10], обладает фокусирующими свойствами.

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

1. ОСЕВОЕ РАСПРЕДЕЛЕНИЕ ПРИ ДИФРАКЦИИ ОГРАНИЧЕННОЙ

ПЛОСКОЙ ВОЛНЫ НА ДРОБНОМ АКСИКОНЕ В НЕПАРАКСИАЛЬНОМ СКАЛЯРНОМ СЛУЧАЕ

Непараксиальная скалярная модель, основанная на теории Рэлея-Зоммерфельда [11], позволяет получать корректные результаты на очень близких расстояниях (всего несколько длин волн) от апертуры [12, 13]. Также в работах [14, 15] было показано качественное согласование результатов, полученных с помощью интегральных методов и конечно-разностных методов, с точностью до учета конкретной конфигурации оптического элемента.

Интегральное преобразование Рэлея-Зом-мер-фельда первого типа в декартовых координатах имеет следующий вид:

E

£0

где E0(x,y) - входное поле,t =,J(u-x)2 +(v-y)2 +Z2, £0 - область, в которой задано входное поле, k = 2к /X - волновое число, l - длина волны.

При рассмотрении поля на оптической оси (u=v=0) выражение (1) значительно упрощает-

ся, и в условиях радиальнои симметрии можно записать:

Е

ехр (к/

г2 + 22) (

I г + 2

к —

I

Г + 2

г сГ

•

у

1-у

0 <у< 1.

Е(0,0,2) = ехр(г'к2) —

2 • ехр

(к2 + 22 -'(ка0К)у

-У (к ао)у 2 •

к ехр пк-^г2 + 22 - г

у-1 г '

г

Сг,

которое удовлетворяет граничному условию при z=0 независимо от погрешности последующего вычисления оставшегося интеграла.

В работе [16] интеграл в (5) был приближенно вычислен модифицированным методом стационарной фазы [2]:

2 • ехр(к/л2 + 2

Е(0,0, 2) ~ ехр(г'кг) - ехр(—(ка0К)У)-, —

У ' л/К2 + 22

-а (2) ехр ^^¡^Щ^ ^ )С/.

где = -год/Иг,)|/2 , /2 = (К - го —г,)|/2 ,

^(го ) = (каого )У- ^го2 +; к

У (го ) = -

г

(1 - у) +-

Рассмотрим дифракцию плоской волны на дифракционном фраксиконе [1]:

г

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

Для 0 < у < 1 оптический элемент (3) на периферии имеет минимальное значение числовой апертуры, а в центре - максимальное. Значение параметра а0 можно определить в соответствии с заданной минимальной числовой апертурой на краю оптического элемента Я:

г0 - стационарная точка, которая определяется из уравнения:

(ка

У-1

Уг0У =

кг

4

2 , 2 г0 + 2

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

При падении на фраксикон (3) плоской волны, ограниченной радиусом Я, с единичной амплитудой, распределение на оси можно описать в приближении тонкого оптического элемента по формуле (2) при Е0 (г) = ехр(-г'(ка0г)У). В работе [16] было получено представление (2):

Выражение (6) достаточно удобно для численного расчета, но мало пригодно (по крайней мере, без дополнительных упрощений) для анализа зависимости комплексной амплитуды (или интенсивности) от расстояния и параметров фраксикона, так как в общем случае уравнение (9) не решается аналитически.

Далее проведём дополнительный анализ полученного решения для 0 <у< 1 в ближней зоне дифракции.

2. ПОПРАВКИ В БЛИЖНЕЙ ЗОНЕ

Можно доказать [16], что при 0 < у < 1 всегда имеется одна стационарная точка (максимум щ(г) ) и нет точек перегиба: щ"(г) < 0 .

В отличие от г0 (2), из (9) легко получается выражение обратной функции 2 (г0), поэтому рассмотрим её поведение.

,2 - 2у

2 2у 0

у О/

Из (10) можно определить границу тени (т.е. правую границу существенных значений на оптической оси), соответствующую г0(2) = К:

=к.1

у о/

В окрестности точки максимума, соответствующей г0 = К, имеется приближённое равенство:

22 ((

>К)« К2

(К

а

2-2у ^

+2К

(К

2-2у

-(2-У)-1

Отсюда имеем приближение для стационарной точки:

2 Г) 2

2 - К

(кК)2

2 2 у

У а0У

--1

2 К | (кК?)2-2У (2-у) -1

2 2 у

У а0У

Зависимость в ближней зоне при 2 ^ 0 можно получить, рассмотрев окрестность точки

\1/

гс = к 1 (у\у

2 (г0 ^гс ) = 2г (1-У) (г0 -гс ) ^ г0с = гс +

2(1 -У) г

2 г =

2

г

Так как при 0 <у< 1 интеграл (5) является несобственным, вклад окрестности начала координат может быть сопоставим с вкладом стационарной точки. Поэтому полученное в первом разделе общее решение (6) требует соответствующих поправок.

Вклад окрестности нуля с учётом осцилляций не превосходит площади нескольких первых лепестков. Оценим эту площадь. В центральной части (при г»0) фазовая функция (7):

y/{r) ~ -kz + ß(kaor)r ,

Четверть периода (15) достигается при r = rq:

1

ka

nr

я

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

exp í-i

Ecor (z) = -Y (k ao)r z exp (ikz) • J-X-f

0

= -iГ(kao)r z exp (ikz)• EIcor (z) .

rr 1 dr =.

Верхний предел в (17) равен нескольким гд. Из вида подынтегральной функции ясно, что поправка не будет пренебрежимо малой только при малых 2. Способ приближённого вычисления интеграла следующий. Запишем его в тригонометрическом виде:

r (z) = J

hr" COS

r

ldr-i J

* sin(kao 1 rr-d

.

!г +г 0 Ыг +г

С учётом последующего приближения делаем лепестки симметричными, что получается при замене г7 = у :

EI

r

r o yjy + z r o y y + z

где г = (Ьг^)r.

Теперь интеграл стал собственным. Отрезок интегрирования разбивается по границам лепестков, внутри каждого лепестка арка синусоиды заменяется аркой параболы (для этого и была нужна симметричность) из условия коллокации: совпадают нули и вершина арки. Знаменатель преобразуется следующим образом: он заменяется на z или на y11 r в зависимости от того, какая величина больше. Естественно, что лепестку, на котором происходит переход с одной замены на другую, будет соответствовать два интеграла. Из-за особого вида интеграла J ypdy при p = -1 ответ для r = 1/ 2 и Г = 1 /3 придётся получить отдельно.

Переписав(19)в виде

EI (z) = CI (z) -iSI (z),

cor ^ s cor ^ s cor ^ ^ 7

выпишем окончательный ответ, выраженный в исходной переменной г . Он, естественно, выражается разными формулами в зависимости от соотношения между г и гч. Ограничимся первыми двумя (при синусе) и тремя (при косинусе) лепестками: С1сог (г) = ^ (г) + Sc2 (г) + 53с (г) (граничные точки - значения, когда под знаком косинуса стоит 0, 90, 270, 450 градусов); Б1сог (г) = (г) + 512 (г) + 52 (г) (граничные точки -значения, когда под знаком синуса стоит 0, 90, 180, 360 градусов).

Лепестки с косинусом, г < гц:

S (z) =

2r

r

rr _1,o<r< 1,

1 4-Jz 4 -r+----r, r=1/2,

Vz 3r, Jrg

4 5 2 5 ln

--—1 r = 1/3

_2/3 .2/3 0..2/3 ' ' '

S2 (z) =

Sí (z) =

rr-1 -j 8

2/3 2/3 z rq 3r,2'3 '

32-1/r -1 32-1' 'r-3 33

2 r-1 r -1 :

8(1n 3-1) o,789

V^

4 - 3 ln 3 r2" ч o,7o4 ,.2,3 ' Я

32-1 /r 3 52-1/r

'-1 r_ 1

8 (21n(5/3)-1) o,173

&

У 1 V ?

8/5-31n(5/3) o, o675

3r- 1

r = 1/2,

r = 1/3,

, o<r<1

3r-1 r = 1/2,

r=1/3.

,o<r<1,

Лепестки с косинусом, r < z < 3 rr :

sc

-rT, o < r < 1,

3 rz q

4» F

r = 1/2,

2r

3z

1/3

r = 1/3,

S2 (z) =

3zr

2z21

r

2 1 ^ z3r-1

2T-T- -1+~

r;

4 1

__^ ^

3r 3r-1 у

---I, o <r< 1

4rr

__-r r- • 32-1/r

3-z q ^ 2r-1 r-1 3r-1

811 (q77) = 1/2

•Tz 3z 3r fi ' '

13,5

18

4,5 4rq'3 31n

--— +---—-", r = 1/3,

z

5^(г) совпадает с(20в). Лепестки с синусом, г < гч:

5Тх(г) =

г27-1 г3'-1 4у—1 л

-—Г- + —-7-Л—- г , 0 <'<1

да—1) 3у(3у—1)^ (2у—1)(3у—1) 4/г 2Цг, /г

3г

у = 1/2,

—, 9 21п(/г) 7=1/3

г?'3 ( Г3 Г ' '

51^2

' 23—1/7 — 2 23—1/7 —

^ 2 7 — 1 3 7 — 1

41п 2 — 2 0, 773

^ Я '

3(1 —1п 2) 0,921

г 2/3 9 г 2/3 , 9

, 0 < ' < 1,

%(г) =

7 = 1/3,

[ 4Ь"' —2—1/' 2-42—1/' —2*—1/' 43—1/' —Г5—1/' 1

-г,—1 -¡24----0<'< 1,

" { 2у—1 '—1 37—1 1

121п2—8 0,316

3(3/4—1п2) 0,171

7=1/2,

7=1/3.

Лепестки с синусом, г < г < 2 7г :

5'( г) =

-гу, 0 < 7 < 1,

37г 9

4/ГГ

3г ,

2г1

7 = 1/2, 7 = 1/3,

ЗД =

2274 (1

(1- 1

7 2у—1) Г 2 \ 3у 3у—1

2г7

_1 +г7— 1 23—1/7

3уг 9 '

1

1

2у—1 3у—1

2 4$ 4Тг 41п

+-+

г 3г 3,

4 2г 9

ЗЦ(/2

0 <у<1,

7 = 1/2, (23б) 7=1/3,

52(г) совпадает с(22в).

3. РЕЗУЛЬТАТЫ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ И АНАЛИЗ НЕОБХОДИМОСТИ ПОПРАВОК В БЛИЖНЕЙ ЗОНЕ.

Рассмотрим прямой расчёт (по квадратурным формулам) дифракции плоской волны в ближней зоне на фраксиконах с высокой числовой апертурой с использованием формулы (2).

Примем в (3) к = 5Л и в (4) = 1. Тогда значение МЛ = 1 будет достигаться на краю оптического элемента, а в центре значение числовой апертуры будет выше и это означает необходимость учитывать при расчете затухающие волны.

По формуле (4) в принятых нами условиях параметр дробного аксикона оценивается следующим образом:

1—7 ( 1 > 7

а0 =(10л) 7 I - I , 0 < у < 1.

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

Как видно из табл. 1 два метода практически не различаются по результатам расчета на расстояниях более длины волны. Однако на очень близких расстояниях ^<0,1 Л ), где влияние затухающих волн очень велико, результаты несколько различаются. Данный факт отмечался также в работе [14].

Заметим также, что если рассматривать осевой участок z>0,5Л (после ямы), то наблюдается поведение, аналогичное отмеченному ранее для фраксиконов большего размера [16]. А именно, максимумы и условная граница тени отодвигаются от элемента при уменьшении 7. Однако на начальном участке имеется принципиальное отличие дробного аксикона ( 0 < 7 < 1) от линейного аксикона (7 = 1): интенсивность не сразу начинает возрастать при увеличении расстояния, а вначале происходит убывание. Причём при 7 , не близком к единице, этот минимум достаточно глубокий. Вычисление по модифицированному

г

г

Таблица 1. Результаты расчёта дифракции для дробного аксикона

Вид фазы элемента

Распределение интенсивности на оптической оси г е [0,1Л, 5Л]

Распределение интенсивности на оптической оси г е [0,001Л, 0,1Л]

Y=1, а о = 1

40

1 2 3 4 5

Максимум (z=0,7A) = 40

0.01 0.03 0.05 0.07 0.09

Минимум (z=0,001A)=1

0.9

0.7

0.5

1 2 3 4 5

Максимум (z=1,7A) = 10

0.01 0.03 0.05 0.07 0.09

Минимум (z=0,04A)=0,48

1.0

0.1

1 2 3 4 5

Максимум (z=2,15A) =5,8

0.01 0.03 0.05 0.07 0.09

Серый: минимум (z=0,034A)=0,01 Черный: минимум (z=0,065X)=0,044

у=0,25, а0 = 7937607

/ /

У

1.0

0.6

0.2

12 3 4

Максимум (z=2,4A) = 4,5

0.01 0.03 0.05 0.07 0.09

Минимум (z=0,055A) = 0,018

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

Для проверки этой возможности выпишем в явном виде выражение (6) при малом z и добавим поправку (17):

Е(0,0,г) ~ ехр(кг) ——ехр\1Ш—1 (ка0К)71 —

ТТехрГ—¡кг^Ц

—¡кг-Г-^—Г | ехр

Сравним слагаемые между собой. Если К >Л/(2л), то второе слагаемое мало по сравнению с первым. Коэффициент

л/2л

М =

д/Г—7 (7,аг )1/(2(1—7)) также невелик, так как в

Д^у

1

знаменателе присутствует величина а0 в положительной степени (например, при 7 = 1/2 ц = 4^л/а0 ), а в данном диапазоне а0 велико. Кроме того, величина коэффициента завышена, так как в ней подставлено значение имеющегося в (25) ин-

5

3

1

9

1

5

1

4

1

5

теграла Френеля при бесконечных пределах.

Поэтому в первом приближении оставим только первое и последнее слагаемое:

E

Так как общий фазовый множитель не влияет на интенсивность, его опустим:

E1

В силу громоздкости выражений для поправочного слагаемого, в общем виде анализировать (26б) не будем. Ограничимся частным значением r = 1/2 . Результаты численного моделирования (третья строка Табл. 1) показывают, что области z < 9rq , для которой имеются формулы, вполне достаточно для выяснения, необходимы ли поправки.

Если воспользоваться формулами из предыдущего пункта, подставить в них значение rq = л2/(4kao) и ввести для удобства переменную 5 = 2yjkaoz /л, то получим следующее кусочно-непрерывное выражение для амплитуды поля:

Ej(o,o, s) --

1 + л (s2 ln s- o, 114s2 -s3/3) + +m(-s+1,154s2 -s3 /3), o < s < 1

1+л(0,58s2 + s21n(s /2)- s3 / 3) + v '

+.л(1/3—3s —o, o43s2 -2s21n(s /3)+s3 / 3) ,1 < s < 2

Анализируя это выражение, можно показать, что интенсивность имеет минимум, равный

0,0478 при s = 2,117 , т.е. при

1,12л

kao

o,o15l .

Этот результат получается средним между двумя вариантами, полученными численными методами. В случае прямого расчета по формуле (2) получается очень близкое значение минимума, а при использовании метода разложения по плоским волнам оказывается более точное положение минимума.

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

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

Поперечное распределение интенсивности в ближней зоне дифракции, рассчитанное на основе метода разложения по плоским волнам для рассмотренных выше оптических элементов, приведено на рис. 1.

Как видно, в плоскости максимальной интенсивности размер центрального светового пятна примерно одинаковый во всех случаях и равен по полуспаду от максимального значения: FWHM=0,385Л для 7=0,8, FWHM=0,395Л для 7=0,5 и FWHM=0,4Л для 7=0,25. Однако с уменьшением параметра g уменьшается интенсивность светового пятна.

Выигрыш использования рассеивающего дробного аксикона с 7<1 проявляется только на очень близких к оптическому элементу расстояниях. В частности, на расстоянии z=0,01Л размер центрального светового пятна равен FWHM=0,266Л для 7=0,8, FWHM=0,123Л для 7=0,5 и FWHM=0,118l для 7=0,25. В этом случае световое пятно оказывается в 5 раз меньше дифракционный предел.

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

1.0

0.6

0.2

-1.0 -0.5

0.0

0.7

0.4

0.1

0.5 1.0

-1.0 -0.5

0.0

0.5 1.0 -1.0

-0.5

0.0

0.5 1.0

Рис. 1. Распределение интенсивности на расстоянии (а) z=0,01 Л , (б) z=0,1 Л и (в) в плоскости максимума интенсивности для 7 =0,8 (серая линия), 7 =0,5 (тонкая черная линия),

7 =0,25

Z

9

5

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

ЗАКЛЮЧЕНИЕ

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

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

Численный расчет на основе интеграла Рэлея-Зоммерфельда и с использованием метода разложения по плоским волнам показал, что дробные аксиконы с 7 <1 и имеющие высокую числовую апертуру в плоскости максимальной интенсивности формируют центральное светового пятна примерно одинакового размера (FWHM=0,39Л ). Однако с уменьшением параметра 7 уменьшается интенсивность светового пятна. Это связано с тем, что элементы в этом случае становятся более высокоапертурными, и основная энергия идет в затухающие волны.

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

Работа выполнена при финансовой поддержке гранта РФФИ 13-07-97004-р_поволжье_а.

СПИСОК ЛИТЕРАТУРЫ

1. Хонина С.Н., Волотовский С.Г. Фраксикон - дифракционный оптический элемент с конической фокальной областью // Компьютерная оптика. 2009.

Т. 33, №4. С. 401-411.

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

3. Annular-aperture logarithmic axicon / J. Sochacki, Z. Jaroszewicz, L.R. Staronski and A. Kolodziejczyk //J. Opt. Soc. Am. A. 1993. V. 10. P. 1765-1768.

4. Устинов А.В., Хонина С.Н. Геометрооптический анализ обобщённой рефракционной линзы // Известия Самарского научного центра РАН, 14(4), 28-37 (2012).

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

6. Extended depth of field and aberration control for inexpensive digital microscope systems / S.C. Tucker, W.T. Cathey andE.R. Dowski// Opt. Express. 1999. Vol. 4(11). P. 467-474.

7. Extending the imaging volume for biometric iris recognition / R. Narayanswamy, G.E. Johnson, P.E.X. Silveira and H.B. Wach // Appl. Opt. 2005. Vol. 44. P. 701-712.

8. Asymmetric phase masks for extended depth of field / A. CastroandJ. Ojeda-Castaneda// Applied Optics. 2004. Vol. 43, N 17. P. 3474-3479.

9. Хонина С.Н. Фазовая аподизация изображающей системы с целью увеличения глубины фокуса в когерентном и некогерентном случаях // Компьютерная оптика, 36

10. Хонина С.Н., Балалаев С.А. Гипергеометрические пуч-

ки в ближней зоне дифракции в рамках скалярной модели // Компьютерная оптика, 33(4), 427-435 (2009).

11. Born M., Wolf E. Principles of Optics. 6th ed. Pergamon, Oxford, 1980. Chap. 8.3.

12. Totzeck M. Validity of the scalar Kirchhoff and Rayleigh-Sommerfeld diffraction theories in the near field of small phase objects // J. Opt. Soc. Am. A. 1991. V. 8, No. 1. P. 27-32.

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

14. Хонина С.Н., Устинов А.В., Ковалев А.А., Волотовский С.Г. Распространение радиально-ограниченных вихревых пучков в ближней зоне: I. Алгоритмы расчёта, Компьютерная оптика, 34(3), 317-332 (2010).

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

16. Устинов А.В., Хонина С.Н. Анализ дифракции плоского пучка на рассеивающем фраксиконе в непараксиальном режиме // Компьютерная оптика, 38(1),

42-50 (2014).

17. Хонина С.Н., Устинов А.В., Волотовский С.Г. Ананьин М.А. Алгоритмы быстрого расчета дифракции

радиально-вихревых лазерных полей на микроапертуре // Известия Самарского научного центра РАН, 12(4), 15-25 (2010).

CALCULATION OF PLANE WAVE DIFFRACTION BY DIVERGENT FRACTIONAL AXICON TAKING INTO ACCOUNT EVANESCENT WAVES

© 2014 A.V. Ustinov, S.N. Khonina

Image Processing Systems Institute of the Russian Academy of Sciences, Samara

We performed theoretical and numerical analysis of plane wave diffraction by fractional axicon whose degree is less than unity, in the near-field taking into account evanescent waves. For calculating the Rayleigh -Sommerfeld integral of first type was used. It is shown that application of the method of stationary phase to obtain analytical estimates in near diffraction zone is insufficiently accurate. For adjustment the obtained solution additional corrections are introduced to obtain qualitative agreement of theoretical and numerical results. It is shown that divergent fractional axicon with degree significantly less than unity can create a light spot having longitudinal and transverse size of one tenth of the wavelength in the zone of evanescent waves.

Keywords: near diffraction zone, divergent fractional axicon, Rayleigh-Sommerfeld integral, modified stationary phase, evanescent waves

Andrey Ustinov, Leading Programmer. E-mail: andr@smr.ru Svetlana Khonina, Doctor of Physics and Mathematics, Leading Research Fellow. E-mail: khonina@smr.ru