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

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

и спиральной фазовой пластинке: сравнение

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

Д.А. Девис, С. Тувей, Д. Коттрел Государственный университет, Сан-Диего, Калифорния

Аннотация

Получены аналитические выражения, описывающие дифракцию Фраунгофера плоской волны конечного радиуса на спиральном аксиконе (СА) и спиральной фазовой пластинке (СФП). Решения получены в виде ряда из функций Бесселя для СА и в виде конечной суммы из функций Бесселя для СФП. Причем решение для СА переходит в решение для СФП, если параметр аксикона приравнять к нулю. На численных примерах показано, что добавление «слабого» аксикона к СФП приводит к снижению контраста концентрических колец в картине дифракции.

Введение

Впервые спиральный аксикон (СА) [1] и спиральная фазовая пластинка (СФП) [2] были изготовлены по технологии фотолитографии и экспериментально исследованы в 1992 году. СА используется для формирования бездифракционных лазерных пучков Бесселя, а СФП используется для формирования оптических вихрей и для оптического выполнения радиального преобразования Гильберта [2-4]. В последнее время возрос интерес к СА и СФП [519]. Это связано с тем, что улучшилось качество изготовления пространственных модуляторов света (ПМС), с помощью которых можно теперь формировать дифракционные оптические элементы, в том числе СА и СФП. Так, в [5] с помощью ПМС сформированы СФП высоких порядков (п > 30) и исследовались оптические вихри высоких порядков. С помощью ПМС можно сформировать составную СФП, которая будет генерировать лазерный пучок, состоящий из нескольких соосных оптических вихрей [6]. Также с помощью ПМС были сформированы бездифракционные пучки Бесселя [7, 8], эллиптические пучки Бесселя [9], пучки Гаусса-Айнса [10] и полые (трубчатые) пучки [11].

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

Теоретическое исследование параксиальной дифракции Френеля и Фраунгофера на СФП было проведено для освещающих гауссового пучка [16], неограниченной плоской волны [2, 12], ограниченной плоской волны [5, 13], эллиптического пучка [17]. Теоретическое исследование дифракции на СА проводилось для неограниченной плоской волны

[14] и гауссового пучка [18]. Интерес к исследованию СА и СФП связан также с их применением для манипулирования микрочастицами [5, 14, 19].

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

2. Описание дифракции Фраунгофера на СА с помощью гипергеометрических функций

Рассмотрим скалярную параксиальную дифракцию плоской волны конечного радиуса Я , комплексная амплитуда которой имеет вид

/ ч i Г i И г < Я

Е0 (г ) = агс!

0, r > R,

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

Дг,ф) = exp(/ar + тф),

где (г,р) - полярные координаты в плоскости СА при г = 0 , г - оптическая ось, а - параметр акси-кона, п = 0,±1,±2,... - номер спиральной фазовой пластинки (СФП).

Тогда комплексная амплитуда света Лп (р,в) в фокальной плоскости идеальной сферической линзы с фокусным расстоянием / имеет вид (опустим несущественный множитель ехр(кг)):

2лг

К 2п

I I ехр

0 0

1аг + 1пр-^угрсо^

гёгёр =

= НГ 7 ехр (гпв)^Щ- } гт+1 3,

I

рг

гёг,

1 т=0 0

где (р,в) - полярные координаты в плоскости дифракции Фраунгофера, к = 2я/2 - волновое число, 3п (х) - функция Бесселя.

Используя справочный интеграл [20]:

К Лт ( ъ (аК),+2+1 I х 3,, (ах )ах = —-—-—/--—:

0 л * aя+12v ( + 2 +1,

1^2

V + 2 +1 V + 2 + 3 , ( аК

-,-, V +1,-1 —

2 2 I 2

вместо выражения (3) получим:

Р„ (р, в) =

'кК2 ^

п!

I

кК р

21

х

\ Z'J /

х

X

(аК )'

=0 (т + п + 2)т!

х Л

т + п + 2 т + п + 4

п + 1, -

'кК Рх2

21

где 1Л2 (а, Ь, с, х) - гипергеометрическая функция:

(а) хт !Л(а,Ь,с,х)= _ ,

т=0 (Ь)т (с)тт!

(а )т = Г(а + т)/Г(а), (а)0 = 1, Г(х) - гамма-функция.

Из выражения (5) следует, что картина дифракции представляет собой набор концентрических колец. При р = 0 в центре картины дифракции при любом п Ф 0 будет нулевая интенсивность. Так как комплексная амплитуда (5) зависит от комбинации переменных кКр/(21), то радиусы р1 локальных максимумов и минимумов картины дифракции должны удовлетворять выражению:

21п

Р1 =

пК

где /1 - постоянные, зависящие только от номера кольца I = 1,2,... картины дифракции.

Из уравнения (5) при а = 0 (аксикон отсутствует) следует выражение для комплексной амплитуды дифракции Фраунгофера плоской волны конечного радиуса на СФП [13]:

Еп (р,в)=(~ ')п+'ехр(п

( кК 2 ^

Х1

(п + 2)п! п + 2 п + 4

I

кКр

2

2

п +1,-

кКр

~2Г

Интересно сравнить выражение (5) с выражением для комплексной амплитуды, описывающей дифракцию Френеля плоской волны конечного радиуса на СФП [13]:

Еп (р, в, г) = 2ехр

/  —2

'Г0 Р

+ тв

Р_ п!

-1Г0

г

^=0 (2т + п + 2)т!

X, Л

2т + п + 2 2т + п + 4

, п +1,

г0 р г

где г0 = кК2 ¡2, р = р/К .

Структура выражений (5) и (9) аналогична и описывает картины дифракции, состоящие из набора концентрических световых колец. Причем, как в выражении (5) при стремлении а к нулю из всей суммы остается только первое слагаемое, также в выражении (9) при стремлении г к бесконечности остается только одно слагаемое, равное (8) при г=/.

О 0,2 0,4 0,6 0,8 Р'ММ

Рис. 1. Радиальное сечение модуля амплитуды Л

дифракции Фраунгофера плоской волны конечного радиуса на СА

На рис. 1 показан рассчитанный по формуле (5) модуль амплитуды |Лп (р,в) в относительных единицах в зависимости от радиальной переменной. Эта кривая показывает радиальное сечение картины дифракции Фраунгофера плоской волны радиуса К = 2 мм и длиной волны 2 = 633 нм на СА с параметром а = 30мм-1 и п = 4. Фокусное расстояние линзы равно I = 140 мм.

п

X

2

2

3. Описание дифракции Фраунгофера на СА

с помощью ряда из функций Бесселя Вместо ряда (5), состоящего из гипергеометрических функций, можно получить для комплексной амплитуды, описывающей дифракцию Фраунгофера плоской волны радиуса Я на СА, выражение в виде ряда из функций Бесселя:

Р (р,в) =

- ехр( аЯ) 2

2п /

(

д1'п да

дГ

да

•/т+п

где

А" =1

ех

р(тр)с1ф

0 а+рсоБф р = кр/, Jm+n (х) - функция Бесселя. Выражения для интеграла I" и производной д1" /да из уравнения (11) можно получить в явном виде:

2п sgn а

Iп

VI —2

а -р

т ((*п -в") _ , , —-Р> 4

I

р

а

где

в =

—2 2 р -а

Р

х = -

7 2 —2

а - р

sgnа =

р

1,а > 0, -1, а < 0,

д1" да

-2пХ

п аsgnа + |п|уа2-р2

"Ф

,0 < р < Ы,

/

!/2

* п| + п И

-2 2 р -а

р > а.

При Я выражение (10) приводит к уравнениям, описывающим дифракцию неограниченной плоской волны на СА, полученным в [14]. Из уравнений (12) и (14) также следует, что поле Рп (р,0) описывается разными функциональными зависимостями при 0 < р < |а| и при р > |а|. Разный характер картины дифракции при изменении р от 0 до |а| и после |а| виден также из рис. 1. Интересно

отметить, что пространственная частота картины дифракции на рис. 1 на интервале 0,1 < р < 0,3 мм примерно в два раза больше, чем на интервале

0,5 < р < 0,7 мм. Заметим также, что глобальный максимум интенсивности 1п (р,в) = |Рп (р,в) возникает при р > |а| (на рис. 1 точка р = \а\ отмечена

вертикальной линией).

Для слабого СА (аЯ << 1) имеет место эффект снижения контраста периферийных колец дифракционной картины. Действительно, при малом аЯ в ряду (5) можно оставить только первых два слагаемых:

Рп

I кЯ 21

п + 2

п!

п + 2 п + 4

/

2

2

п +1,-х 2 1 +

п + 3

г аЯ Г п + 3 п + 5 , 2 . +-Р I -,-, п +1,-х 1+ ...

2

2

где х = кЯр/(2/). Первое слагаемое в уравнении (15) совпадает с выражением (8) - оно действительное, второе слагаемое в уравнении (15) - только мнимое. Поэтому интенсивность 1п (р,в) будет пропорциональна сумме квадратов этих двух слагаемых. Это означает, что нули интенсивности будут иметь место в тех точках р , в которых оба слагаемых в выражении (15) равны нулю. Это приведет к снижению контраста периферийных колец картины дифракции.

а)

б)

0 0,1 0,2 0,3 0,4 р, мм 0

0,1 0,2 0,3 0,4 Р, мм

в) г)

Рис. 2. Рассчитанные картины дифракции Фраунгофера |Рп (р, 0) плоской волны

конечного радиуса на СФП п = 4

хп X

т

т=-ы

1

X

На рис. 2 показаны две рассчитанных картины дифракции Фраунгофера (2Б модуль амплитуды |Еп (р,в), I = 140 мм) плоской волны конечного радиуса К = 2 мм с длиной волны 2 = 633 нм на СФП (а = 0, п = 4) (а) и СА (а = 1мм-1, п = 4) (б). Размер обеих картин дифракции 0,45х0,45 мм. На рис. 2вг приведены радиальные сечения модулей амплитуды картин дифракции, показанных на рис. 2аб.

Из рис. 2 видно, что если при а = 0 число локальных максимумов (кроме главного) функции модуля амплитуды (р,в) в диапазоне 0 < р < 0,45 мм равно 8, и большая часть минимумов равна нулю (рис. 2в), то при а = 1 мм-1 увеличивается число локальных максимумов амплитуды примерно вдвое, почти все минимумы не равны нулю, и поэтому контраст картины дифракции в области боковых лепестков существенно снижается.

4. Описание дифракции Фраунгофера на СФП с помощью конечных сумм функций Бесселя

Чтобы получить выражение для дифракции Фраунгофера плоской волны радиуса К на СФП с целым номером п в уравнениях (10)-(14) следует положить а = 0. Тогда интегралы (12) и (14) существенно упростятся:

0, п = 2т,

II =

дГх да

Р

2n i? I

n = 2m +1;

n = 2m,

Р

0, n = 2m + 1.

Подставляя выражения (16), (17) в уравнение (10), получим выражение для комплексной амплитуды дифракции Фраунгофера плоской волны конечного радиуса на СФП:

, „ч (~1 )п+1 к ехр (тв) Еп (р, в) = У ' _ Г ' X

f Р2

(n-2)/2

1-Jo {y)~ 2 2 ^ 2 m Ы

m=1

' У (n-1)/2

f Jo {t) dt - 2 ^ J2m—1 {У )

где y = Rp = kRp/ f,

— yJn—1 (y ), n = 2m, (18) yJn—1 Ы.n = 2m + 1,

13 0 = 2 {{ ((У У)Н 0 (у) + 3 0 (у)[2-пЯ1(у)] (19) 0 2

Н 01 (у) - функции Струве нулевого и первого порядков.

Выражение (18) описывает дифракцию на СФП в виде конечных сумм функций Бесселя. Для четных и нечетных номеров СФП дифракционная картина описывается разными выражениями. Для пер-

вых трех номеров n = 0,1, 2 из (18) следуют про стые выражения:

¿0 (p,0)-=-f .iW,

fy

E p)=-fe

f p

\ J 0

¿2

f p

\E\ 40

30

20

10

0,1 0,2 0,3 0,4 p.MM

в) 0 0,5 1,0 1,5 p,MM

Рис. 3. Рассчитанные модули комплексной амплитуды |En

Выражение (20) описывает обычную картину дифракции Фраунгофера плоской волны на круглой апертуре радиуса R . На рис. 3 показаны рассчитанные модули амплитуды |En (p,e) радиальных сече-

х

ний картин дифракции Фраунгофера ( f = 100 мм) плоской волны радиуса R = 1 мм и с длиной волны Я = 633 нм на СФП с номерами n = 2,16,50 .

Из рис. 3 видно, что с ростом n уменьшается контраст боковых лепестков модуля амплитуды. Минимумы модуля амплитуды боковых лепестков не равны нулю. Из уравнения (18) можно асимптотически оценить радиус первого кольца, на котором модуль амплитуды поля равен нулю. При p^-œ (n=2m) из уравнения (18) следует приближенное равенство:

E (p в)k[n-k E^ f p2 fpp

2 y . ( n +1 Sinl y--п

п К 2

В уравнении (23) воспользовались асимптотическим преставлением функции Бесселя при большом значении аргумента. Из уравнения (23) следует, что примерная зависимость радиуса р0, при котором

имеет место первый нуль амплитуды |Еп (р, 0), имеет вид (п>>1):

А/п2

p0

4R

Из уравнения (24) видно, что зависимость радиуса для первого нуля амплитуды от номера СФП квадратичная. Это согласуется с рис. 3, так как видно, что с ростом п растет значение р0, при котором имеет место первый нуль амплитуды.

Из уравнения (18) можно получить рекуррентное соотношение, с помощью которого быстро и удобно рассчитывать комплексные амплитуды при разных п :

E'+2 (p,0)=-

n+2

-exi

p(i20)E„ (p,ß)-

+ 2(- i Г exp[i(n + 2 J ( -k-Rpl.

p n К f J

\En (pm A = 2 A ( ^ J J'

pm К ' + 2 J

р п

Уравнение (25) верно для всех целых п . Из уравнения (25) можно получить уравнение, которому должны удовлетворять координаты рт точек, в которых амплитуда поля |Еп+2(р,0) равна нулю(п>0):

Я I п + 11 I кЯрт , п + 2 М /

Из рис. 4 видно, что максимумы модуля амплитуды |Еп (р,00) почти совпадают с максимумами

функции 2 Я/ р-(п +1)1 (п + 2) • /п+1 ((Яр /), и наименьшая координата р0, при которой максимумы этих двух функций совпадут точно, и будет первым нулем амплитуды Еп+2 (р,0). Из сравнения рис. 4а и рис. 4б видно, что с ростом номера п увеличивается значение р0.

Интересно сравнить уравнение (26) с уравнением [13], которому должны удовлетворять координаты р!1 минимумов и максимумов амплитуды поля Еп (р, 0):

\En М =

2

б) 0 0,05 0,10 0,15 0,20

Рис. 4. Графики функций (n=2 (а) и n=4 (б)) из левой и правой частей в уравнении (26): толстая кривая -график функции |En (р, 0), тонкая кривая - график

функции 2R¡р• (n +\)/(n + 2)- Jn+l (kRp¡f). Параметры расчета: R=1mm, f=100 мм, 1=633 нм

-10

Рис. 5. Графики функций из левой и правой частей в уравнении (27): толстая кривая - график функции |Еп (р,0), тонкая кривая - график функции

(кЯ 2 )/(2/)• /п (кЯр/). Параметры расчета: п=14, Я=1мм,/=100 мм, 1=633 нм

n

n

Из уравнения (27) видно, что для интенсивности основного (первого) максимума картины дифракции Фраунгофера имеет место неравенство:

\En (Ps ,ef *

< kR 2 ^

2 f

J,

0n

где /0п = тах|/п (У)2.

На рис. 5 показаны два графика функций из левой и правой частей уравнения (27) при следующих параметрах: п=14, Я=1 мм, /=100 мм, 1=633 нм. Видно, что графики пересекаются в точках, в которых амплитуда |Еп(р,,0) имеет экстремумы.

Из рис. 5 видно также, что координата ру первого максимума амплитуды |Епрх,0)

находится между первым нулем производной функции Бесселя п-го порядка и между первым нулем самой функции Бесселя п-го порядка:

fY'n ,1 fYn,1 -- < pv <--

kR v kR

где у„1 и у'„1 - первые корни функции Бесселя п -го порядка и ее производной: /п (упЛ) = /'„ (у'п,\) = 0 . Из неравенств (29) можно получить приближенное ра-

венство для определения радиуса первого максимума амплитуды |Еп (р,, 0) :

Pv

Ynf - =Yn,1 +fn ,1

7n

кЯ ' "' 2

Интересно сравнить выражение (30) с аналогичным выражением для радиуса воронки, полученным в [13]:

Гп-1,1/

Pv =-

kR

где у„.11 - корень функции Бесселя (п - 1)-го порядка: /п-1 (хп-1,1) = 0. В таблице показаны значения радиусов максимальной амплитуды (радиусов воронок) для п от 1 до 8, вычисленных с помощью уравнений (18) (первая строка), (30) (вторая строка) и (31) (третья строка). Радиус максимальной амплитуды с помощью уравнения (18) находился численными методами и считается более точным, а выражения (30) и (31) дают приближенное значение этого радиуса. Из Таблицы видно, что оба выражения дают примерно одинаковую оценку радиуса воронки (при Я=1 мм, /=100 мм, 1=633 нм) с максимальной относительной ошибкой 4% (п>1), хотя выражение (31) немного точнее (максимальная ошибка 3%).

Таблица. Приближенные значения радиуса pvмаксимальной амплитуды |En

N 1 2 3 4 5 6 7 8

Pv (мкм),ур. (18) 24,6 39,6 53,0 65,8 78,2 90,4 102,4 114,0

р„(мкм),ур. (30) 28,6 41,2 53,3 64,9 76,5 87,8 99,1 110,2

р„(мкм),ур. (31) 24,2 38,6 51,7 64,3 76,4 88,4 100,0 111,7

5. Сравнение теории и эксперимента

Эксперименты по дифракции Фраунгофера плоской волны конечного радиуса на СФП проводились с помощью программируемого фазового жидко-кристаллического пространственного модулятора света (ПМС). Свет от аргонового лазера с длиной волны 514 нм коллимировался и падал на поверхность ПМС. Электрически управляемый нематический ПМС фирмы Seiko Epson с разрешением 640x480 пикселов имеет размер по диагонали 3,3 см. Пропускание ПМС было пропорционально функции (2) СФП с номером n=10 и а=0 мм-1, а также СА с номером n=10 и а=11 мм-1. На рис. 6 показаны картины дифракции Фраунгофера плоской волны на СФП (n=10) (а) и ее вертикальное сечение (б) и на СА (n=10, а=11мм-1) (в), ее вертикальное сечение (г), зарегистрированные CCD-камерой в фокальной плоскости сферической линзы с фокусным расстоянием /=100 мм.

Сравнивая рисунки 66 и 6г можно сделать вывод, что добавление слабого аксикона к СФП приводит к снижению примерно в 2 раза интенсивности пери-

ферийных колец картины дифракции. Зарегистрированные интенсивности в вертикальном сечении обеих картин дифракции хорошо согласуются с теоретическими кривыми интенсивности.

Заключение В работе получены аналитические выражения для параксиальной дифракции Фраунгофера плоской волны конечно радиуса на СА и СФП с произвольным целым п. Для СА комплексная амплитуда описывается либо рядом из гипергеометрических функций (уравнение (5)), либо рядом из функций Бесселя (уравнение (10)), а для СФП - конечной суммой из функций Бесселя (уравнение (18)). Причем комплексные амплитуды картины дифракции для СА описываются разными аналитическими выражениями внутри и вне некоторого радиуса (уравнения (12), (14)). Картина дифракции для СФП также описывается разными суммами для четных и нечетных номеров п. Численно показано, что при добавлении слабого аксикона к СФП контраст периферийных колец картины дифракции Фраунгофера существенно снижается.

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

Работа выполнена при поддержке российско-американской программы «Фундаментальные ис-

следования и высшее образование» (грант СКОБ ИБС БЛ-014-02), а также гранта РФФИ 05-0850298.

а)

в)

Я

V

L

100

200

300

400

0,4

0,3

0,2

0,1

0

I

A \

\\ ]

ТлЛ/

100

200

300

400

б) г)

Рис. 6. Картины дифракции Фраунгофера плоской волны на СФП (п=10) (а) и ее вертикальное сечение (б) и на СА (п=10, а=11мм1) (в), ее вертикальное сечение (г), зарегистрированные ССБ-камерой в фокальной плоскости сферической линзы с фокусным расстоянием 1=100 мм. А - эксперимент, В - теория. По горизонтальной оси - число пикселей ПМС, по вертикальной оси - относительные значения интенсивности

Литература

1. Khonina S.N., Kotlyar V.V., Soifer V.A., Shinkaryev M.V., Uspleniev G.V. Trochoson // Opt. Commun., V.91. No.3-4. Pp.158-162 (1992).

2. Khonina S.N., Kotlyar V.V., Shinkarev M.V., Soifer V.A., Uspleniev G.V. The phase rotor filter // J. Mod. Opt.,V.39. No.5. P. 1147-1154 (1992).

3. Davis J.A., McNamara E., Cottrell D.M., Campos J., Image processing with the radial Hilbert transform: theory and experiments // Opt. Lett., V. 25. No.2. P. 99-101 (2000).

4. Guo C., Han Y., Xu J., Ding J.

5. Radial Hilbert transform with Laguerre-Gaussian spatial filters // Opt. Lett., V.31. No.10. Pp.1394-1396 (2006).

6. Curtis J.E., Grier D.G. Structure of optical vortices // Phys. Rev. Lett., V.90, No.13, P. 1.33901-1-3 (2003).

7. Lin J., Yuan X., Tao S.H., Burge R.E. Synthesis of multiple collinear helical modes generated by a phase-only element. J. Opt. Soc. Am. A., V.23. No.5. P. 1214-1218 (2006).

8. Chattrapiban N., Rogers E.A., Cofield D., Hill W.T., R. Roy Generation of non-diffracting Bessel beams by use of a spatial light modulator, Opt. Lett., V.28. No.22. P. 2183-2185 (2003).

9. Hakola A., Shevchenko A., Buchter S.C., Kaivola M., Tabiryan N.V. Creation of a narrow Bessel-like laser beam using a nematic liquid crystal // J. Opt. Soc. Am. B, V.23, No.4, P. 637-641 (2006).

10. Chakraborty R., Ghosh A. Generation of an elliptic Bessel beam II Opt. Lett., V.31. No.1. P. ЗВ-40 (2006).

11. Bentley J.B. , Davis J.A. , Bandres M.A., Gutiérrez-Vega J.C. Generation of helical Ince-Gaussian beams with a liquid-crystal display II Opt. Lett., V.31, No.5, Pp.649651 (2006).

12. Fatemi F.K., Bashkansky M. Generation of hollow beams by using a binary spatial light modulator II Opt. Lett., V.31. No. 7. P.864-866 (2006).

13. Kotlyar V.V., Almazov A.A., Khonina S.N., Soifer V.A., Elfstrom H., Turunen J. Generation of phase singularity through diffracting a plane or Gaussian beam by a spiral phase plate II J. Opt. Soc. Am. A., V.22, No.5, Pp.849-861 (2005).

14. Kotlyar V.V., Khonina S.N., Kovalev A.A., Soifer V.A., Elfstrom H., Turunen J. Diffraction of a plane, finite-radius wave by a spiral phase plate II Opt. Lett., V.31, No.11, Pp.1597-1599 (2006).

15. Kotlyar V.V., Kovalev A.A., Khonina S.N., Skidanov R.V., Soifer V.A., Elfstrom H., Tossavainen N., Turunen J. Diffraction of conic and Gaussian beams by a spiral phase plate II Appl. Opt., V.45, No.12, Pp.2656-2665 (2006).

16. Ahluwalia B.P.S., Cheong W.C., Yuan X.-C., Zhang L.-S., Tao S.-H., Bu J., Wang H. Design and fabrication of a double-axicon for generation of tailorable self-imaged three-dimensional intensity voids // Opt. Lett., V.31. No.7. P. 987-989 (2006).

17. Rozas D., Law C.T., Swartzlander G.A. Propagation dynamics of optical vortices // J. Opt. Soc. Am. B, Vol.14, No.11, Pp.3054-3065 (1997).

18. Dennis M.R. Rows of optical vortices from elliptically perturbing a high-order beam, Opt. Lett., V.31, No.9, Pp.1325-1327 (2006).

19. Ling D., J. Li, J. Chen Analysis of eigenfields in the axi-con-based Bessel-Gauss resonator by the transfer-matrix method // J. Opt. Soc. Am. A., V.23, No.4, Pp.912-918 (2006).

20. Cheong W.C., Lee W.M., Yuan X.-C., Zhang L.-S., Dholakia K., Wang H. Direct electron-beam writing of continuous spiral phase plates in negative resist with high power efficiency for optical manipulation // Appl. Phys. Lett., V. 85. No.23. P. 5784-5786 (2004).

21. Prudnikov A.P., Brichkov Y.A., Marichev O.I. Integrals and Series // Special Functions. Moscow, Nauka.1983.