ДИФРАКЦИОННАЯ ОПТИКА, ОПТИЧЕСКИЕ ТЕХНОЛОГИИ
НЕПАРАКСИАЛЬНОЕ РАСПРОСТРАНЕНИЕ ВЕКТОРНОГО ГАУССОВА ОПТИЧЕСКОГО ВИХРЯ
С НАЧАЛЬНОЙ РАДИАЛЬНОЙ ПОЛЯРИЗАЦИЕЙ
Алексей Андреевич Ковалев12 (научный сотрудник, e-mail: alanko@smr.ru), Виктор Викторович Котляр1'2 (заведующий лабораторией, e-mail: kotlyar@smr.ru) 1 Самарский государственный аэрокосмический университет имени академика С.П. Королева, 2 Учреждение Российской академии наук Институт систем обработки изображений РАН
Аннотация
Исследована непараксиальная дифракция гауссовых оптических вихрей с начальной радиальной поляризацией. Получены явные аналитические выражения для радиальной, азимутальной и продольной компонент вектора напряженности электрического поля. Моделирование показало, что непараксиальные формулы намного точнее при удалении от оптической оси, чем параксиальные. Дифракционное фокусное пятно по полуспаду интенсивности при фокусировке гауссового вихревого поля с топологическим единичным зарядом и радиальной начальной поляризацией получается меньше, чем значение, предсказываемое скалярной теорией.
Ключевые слова: непараксиальная дифракция, радиальная поляризация, гауссов оптический вихрь, фокусировка лазерного света.
Введение
Радиально поляризованные пучки могут быть сгенерированы с помощью пространственно-переменных диэлектрических решеток [1], вращаю-щих изображение резонаторов [2], интерференционными методами [3] или с помощью оптических волокон [4]. Интерес к лазерным пучкам с радиальной или азимутальной поляризацией появился из-за таких потенциальных применений, как лазерная резка [5] и оптическая микроманипуляция [6], поскольку фокусировка таких пучков приводит к формированию круглого субволнового фокусного пятна.
Распространение радиально-поляризованных пучков рассматривается в [7-9]. В [7] не была изучена продольная компонента поля. В [8] исследована непараксиальная дифракция, но рассматривалось поле без оптических вихрей. В [9] рассмотрена дифракция гауссова пучка с радиальной поляризацией, но в рамках параксиальной теории и также без оптических вихрей. Распространение параксиальных и непараксиальных вихревых пучков рассматривалось в [10-15]. В [10-13] рассматривалась скалярная параксиальная дифракция вихревых пучков, в [14] -
векторная непараксиальная дифракция вихревого пучка с линейной поляризацией и в [15] - скалярная непараксиальная дифракция гипергеометрических лазерных пучков. Острая (непараксиальная) фокусировка невихревых лазерных пучков с радиальной поляризацией и вихревых пучков с эллиптической поляризацией рассматривалась соответственно в [16] и [17]. В основе рассмотрения в [16, 17] были приближенные теории Дебая и Ричардса-Вольфа.
В данной работе получены явные аналитические выражения для радиальной , азимутальной и продольной проекций вектора напряженности электрического поля, описывающие непараксиальную дифракцию гауссовых оптических вихрей (вихревых лазерных пучков) с начальной радиальной поляризацией.
1. Интегральные преобразования для описания распространения радиально и азимутально поляризованных лазерных пучков в свободном пространстве
Распространение света в свободном пространстве вдоль оси 02 описывается хорошо известными интегралами Рэлея-Зоммерфельда:
Ex (u, v, z )=- Я Ex (x, y,0)fz
Ey(u,v,z ) = - 2P ff Ey(x, y,° £
Ez(u,v,z )=2P ff b(x, *0)|
exp (ikR) R
exp (ikR) R
exp (ikR) R
dxdy, dxdy,
+ E (u,v,0)— У dy
(1)
exp (ikR)
R
>dxdy,
где R2 = 21 +(х - и) +(у - V) , (х, у) - декартовы
координаты во входной плоскости, а (и, V) - декартовы координаты в выходной, Ех,у(х, у, 0) - проекции вектора напряженности электрического поля электромагнитной волны на оси х и у в плоскости 2 = 0,
Ех,у(и, V, 2) - проекции вектора напряженности электрического поля на оси х, у, 2 но в произвольной плоскости 2.
От декартовых составляющих электрического вектора вдоль осей х и у перейдем к радиальной и азимутальной составляющим:
или
где
" Er'
_ Ej _
" Ex'
E _ y _
cos j sin j - sin j cos j
cos j - sin j sin j cos j
E
E
cos j =
Í
2 , 2 x + y
, sin j =
y
V
22 x + y
(2)
(3)
(4)
Тогда для радиальной и азимутальной составляющих верны следующие выражения (в цилиндрических координатах):
E (р,0,z) = -JJ[Er (r, j,0)cos(j-0)-
-Ej (r, j,0)sin(j-0)]|-
exp (ikR)
R
(5)
rdrdj,
Ej (р, 0, z )=-JJ[ Ej (r, j,0) cos (j-0) +
+ Er (r, j,0)sin(j-0)]— J dz
exp (ikR)
R
(6)
rdrdj.
Видно, что если в начальной плоскости поле имело радиально-симметричное распределение и было радиально или азимутально поляризовано, то оно останется таким же радиально или азимутально поляризованным при распространении. Это следует из того, что в полярных координатах Я2 = z2 +р2 + r2 -2prcos(ф-0) и интегрирование
по переменой ^дает
J sin(j-0)|-
dz
exp (ikR)
R
dj =
2p -л
= -J "
J„ dz
exp (ikR)
R
dcos (j-0) =
= -C [cos (j-0)]| j"o = a где c(X) - первообразная от функции
dz
exp (ik\jz2 +р2 + r2 -
2prX)
4z2 +р2 + r2 - 2prX
Общий множитель в подынтегральных выражениях (6) можно приближенно упростить:
dz
exp (ikR)
R
x exp
ikr2
= ^i+zp2exp Ш-
z +р \
Л
z2 +р2 )x
2Л/
2 , „2 z +р
exp
ikpr cos (j-0)
J
22 z +р
(7)
Это приближение непараксиальное, так как в знаменателе дроби, стоящей в показателе экспоненты, присутствует не только продольная координата г, но и поперечная координата р в плоскости наблюдения.
Тогда выражения для радиальной и азимутальной составляющих упрощаются, и вместо (6) получим:
Er fo 0 z) = - /'k,Z 2\exp ('Wz2 +р2 )x
2p( z +р2) v >
xJJ[ Er (r, j+ 0,0) cos j- Ej (r, j+ 0,0) sin j] x (8)
x exp
ikr
\
24 z2 +р2
exp
ikz
Г Л
ikpr cos j
rdrdj,
/
-exp (ik^z2 + p2 )x
Ej(P, 0, z )=- - ( 2 +p2 V
2p(z +p )
xJJ [Ej (r, j + 0,0) cos j + Er (r, j + 0,0) sin j] x (9)
x exp
ikr
24 z2 +p2
exp
Г Л
ik pr cos j
VzW
rdrdj.
(10)
2. Распространение оптических вихрей с начальной радиальной поляризацией Если в начальной плоскости радиально поляризованное поле с вихревой составляющей, т.е. IX (г, Ф,0) = Лг (г)ехр(гпф), ф, ф,0) = 0,
то из (5), (6) следует, что в общем случае (п ф 0) компоненты Ег и Еф будут отличны от нуля. В этом разделе приводится подробный вывод выражения для Ег , и только в конце дается без вывода выражение для Еф:
кг
E (р, 0, z) = (-i)n
2 (z2 +p2)'
x exp (ik4z2 +p2 + in0) x
xJ Ar (r) exp
Г ikr2 Л
2
(11)
J
kpr
22 Wz +P J
Л
JzT+P2
( 1 ^ kpr
rdr.
В работе [9] рассмотрено распространение радиально поляризованного гауссова пучка с таким распределением в начальной плоскости:
' ^
Ex (x,y,0) = E0 - xexp[-a(x2 + y2)]
V w0 у Г /ТЛ
л/2
Ey (x,y,0) = E0 — y exp [-a( x2 + y2)]
(12)
V w0 у
где а = 1/ а>2 - гк/ (2/), (О0 - радиус перетяжки гауссова пучка, к - волновое число, а/- фокусное расстояние линзы.
Поле с декартовыми составляющими (12) соответствует полю с радиальной составляющей
R
(
x
R
(
\
R
0
x
0
x
(
f „д л
Er (r, j,°) = E° — r exp(-ar2).
V w° У
(13)
Подставляя это выражение вместо Ar(r) в (11), получим:
, , (—i)n kzE° Er (р, 9, z )= ^ ! . x
л/2ю° (z2 + ~2
P2)
x exp (ik^z2 +р2 + in 9) x
xf r2 exp(-pr2) [ Jn_x (gr)- Jn+i (Г)]dr,
(14)
где
g =
kp
i
2 , „2 z +p
b = a-
ik
aq
24 z2 + p2 V
22 z +p
ik
Г2 2 lk q =y z2 +p2--
2a
Интеграл в (18) можно вычислить с помощью следующего справочного интеграла [9, выражение 6.618.1]:
f exp (-ax2) Jv (bx) dx = —^ exp (-y) (y), (15) J 2v a -
где Rev >-1, Re a> 0, b> 0, y = b2/(8a), In (x) -функция Бесселя второго рода. Дифференцируя (15) по a, можно получить:
f x2 exp (-ax2) Jv (bx) dx =
0
= —exp (- y )x
4«a
(v-1 + 2y)Iv(y)- 2yIn_2 (y)
2
(16)
Используя этот интеграл, а также рекуррентное соотношение для функций 1п (2) [19, выражение 9.6.26], получим:
Er (р, 9, z) = -
(-i)n 4%kzE0
W^bV^ z2 +р2)
(n - 4x) In_-1 (x) + (n + 4x) In±1 (x) 2 2 .
xexp (ik^/z2 +p2 - x + in9),
(17)
где х = у7 (8р).
Аналогично можно получить выражение для азимутальной составляющей:
(—)"+' ^л/якгЕ0
Ej(ip 9, z )="W2w,b/b( z2 +p2 )'
In+1 (x) + 4-1 (x)
2
xexp (ik^z2 +p2 - x + in9),
(18)
Выражения (17) и (18) являются общими для описания гауссова оптического вихря произвольного порядка п с радиальной начальной поляризацией (10).
В [9] рассматривался случай, когда п = 0. Его также можно получить из уравнения (17) (как следует из (18), в этом случае Еф = 0):
Er (p, 9, z) =
-k2 zE0p
2л/2|
w0a2q2 л/ z2 + p2
x exp
ik^z1 + p2
22 k p
Л
(19)
4aq^z2 +p2
В параксиальном приближении q » 2 — 1к! (2а),
(22 +р2 )12 » 2 +р7(22) » 2 . Следовательно,
Ег (р, 9,2) = —кЕр ехр(1к2)ехрМ. (20) 2^J2q а2Ю0 ^ 2q )
Подобное выражение было получено и в работе [9] (формулы (27)), оно отличается от (20) только комплексной сопряженностью.
3. Выражения для продольной составляющей
Продольная составляющая вычисляется с помощью третьего интеграла в (1).
При переходе к радиальной и азимутальной компонентам получится следующее выражение (в цилиндрических координатах):
Ez (p, 9, z ) = -Р ff {Er (r, j,0) [r-p cos (<p-9)] + E<(r, j,0 )p sin (j-9)}R ^R
2р R^
exp (ikR)
R
rdrdj.
(21)
В частности, если в начальной плоскости поле было радиально поляризовано, то с использованием (7) можно получить:
Е2 (р,9,2)» ^ / ^ехр(¡к^22 + р2 )х
2р( z2 +p2)
xff Er (r, j,0)[r -pcos(j-9)]exp
ikr
bj z2 +p2
exp
ikpr cos (j-9)
(22)
rdrd< .
x
0
2
x
с
0
с
\
R
Если в начальной плоскости радиально поляризованное вихревое поле (10), то по аналогии с выражениями для поперечных составляющих можно показать, что для продольной составляющей получится выражение:
Ez (р,9,z) = (-/)" —2--exp(ikjz2 +р2 + /«е)х12Jг2ЛГ (г)exp
2(z +р ) I о ^2\]
' ikr2 ^ '
2J z2 +р2
крг
Л
i
2 , „2 z +р
dr +
+ i^JгЛг (г)
(
г) exp
ikr2
2 , „2 z +р
( 1 ^ крг
22 z +р
I
- Jn-1
( 1 ^ крг
22 z +р
I
dr!
(23)
Рассмотрим распространение радиально поляризованного гауссова пучка (13). Тогда
Ez (р,9,z) = (-i)n
ikE„
V2rn0 (z2 + ~2
р2 )
exp
(ik^/z2 +р2 + ine) х
хЫг3 exp(-рг2) Jn (уг)dг + г2 exp(-рг2)[Jn+X (уг)-Jn_x (уг)]dr
где
ß = a-
ik
\I z2 +р2
У =
kр
4.
z2 +р2
Второй интеграл имеет ту же форму, что и (16). Первый интеграл берется при дифференцировании по а следующего справочного интеграла [18, выражение 6.631.7]:
J x exp (-ax2) Jv(ßx) dx = ^apr х
х exp I
.ßl 8a
/ i-ß1! - / 418a ) ^ 18a
(25)
где Яе а > 0, Яе V > -2.
С помощью (25) можно получить:
¥ х3ехр (-рх2) Jv(gx) ёх = ехр (- )х
0 8р V р
I v-1 ( t )- ( t )
t
+— 2
Iv-3 (t)-3/v-! (t) + 3/v+1 (t)-Iv + 3 (t)
2
(26)
где t = у 7 (8ß).
Подставляя выражения (25), (26) в (24), получим:
Ez (р, 9, z) = (-г ^ikE;- х
, ' { ' W^ß2^ (z2+р2)
хexp (ik^Jz2 +р2 -t + in9) х
In-1 (t)- In + 1 (t)
2 2 _
In-3 (t)- 3C, (t) + 3I£±i (t)- In+3 (t) 2 2 2 2
-,ßP (n - 4t)In-1 (t) + (n + 4t)In+1 (t)
+ Г 2
n + 1
2
(27)
(24)
Заметим, что выражение (27) похоже по своей структуре на выражение, полученное авторами в [14] для гауссового оптического вихря с эллиптической поляризацией.
Используя рекуррентные соотношения для модифицированных функций Бесселя и приводя подобные слагаемые, получим окончательно:
(-г )п \1%1кЕ0
E (р e, z )= W2ß2Vß»o (z2+Р2 )
хexp (ik^z2 + р2 -1+ш9)х у (n + 2)
(28)
2
У( n - 2) 2
- 2yt - /'Ьр(п - 4t)
In-1 (t)
+ 2yt - /'Ьр( n + 4t)
In+1 (t)
Если пучок во входной плоскости невихревой (т.е. п = 0), то
Ez (р,9,z) =
ikEn
V2
q 2a2rn0
1 + M ! х
2q
х exp
ik^z2 + р2
k 2р2
(29)
4qa\j z2 +р2
В параксиальном приближении верно следующее равенство:
к 2р2
ik^/z2 + р2
4qa/ z2 +р2
.ib+- Ap_ = +ik р
2 z 4qaz 2q
Используя его, получим:
Ez (р,9,z) =
ikE
s,
q2a 2w0
х I 1+
ikр
exp | ikz +
ikр
(30)
2д ) ^ 2д
Это выражение так же совпадает (с точностью до комплексного сопряжения) с выражением для продольной компоненты радиально-по ляризован-
Л
0
0
0
2
х
+
2
(
\
2
2
х
2
2
2
х
ного гауссова пучка, полученным в работе [9] (формула (32)).
4. Численное моделирование
При численном моделировании интегралы (1) вычислялись по методу прямоугольников и сравнивались со значениями, получаемыми при использовании формул (19) и (29) и формул, получаемых из (19) и (29) для параксиального приближения.
Для расчетов были использованы следующие значения параметров: длина волны 1 = 633 нм, мощность пучка Е0 = 1 Вт/см2, радиус перетяжки гауссова пучка а = 1 мкм, фокусное расстояние параболической линзы /= 5 мкм, расстояние, пройденное светом г = 5 мкм (т.е. поле рассматривалось в фокусе линзы), полярный угол в выходной плоскости в= 0, пределы интегрирования для вычисления поля по формулам (1) х, у е [-4а0, +4а0 ] х [-4а0, +4а0 ], шаг дискретизации
по обеим переменным 8а№, N = 300.
На рис. 1а показана зависимость радиальной составляющей Ег (вещественной части) от радиальной координаты р. Кривая 1 соответствует полю, вычисленному по формуле (1), кривая 2 - по непараксиальной формуле (19), кривая 3 - по параксиальной формуле (20). На рис. 1б показаны те же зависимости для компоненты Ег: кривая 1 соответствует полю, вычисленному по формуле (1) (третий интеграл), кривая 2 - по непараксиальной формуле (29), кривая 3 - по параксиальной формуле (30). На рис. 1 в показана интенсивность I = Е2 + Е2, рассчитанная по формулам (1).
Из рис. 1 видно , что ближе к оптической оси параксиальные формулы (20) и (30) дают распределение амплитуды, близкое к значениям, даваемым интегралами (1), однако при удалении от оси непараксиальные формулы намного точнее.
Аналогичные сравнительные исследования проводились и для пучка (13) с добавлением оптического вихря с п = 1. На рис. 2 также показаны распределения амплитуд компонент Ег и Ег, а также интенсивности I. Кривая 1 соответствует полю, вычисленному по формуле (1), кривая 2 - по непараксиальным формулам (17) и (28).
Расчет интенсивности (рис. 2в) проведен по формуле (1), в ней автоматически учтена составляющая Еф, не показанная на рис. 2.
Интересно отметить, что невихревой гауссов пучок (п = 0) формирует в фокусе кольцевое распределение интенсивности (хотя I ф 0 при р = 0), а вихревой гауссов пучок (при п = 1) собирается в круглое фокальное пятно.
Заключение
В работе исследована непараксиальная дифракция гауссовых оптических вихрей с начальной радиальной поляризацией. Получены выражения для радиальной, азимутальной и продольной компонент вектора напряженности электрического поля. При п = 0 в пара-
ксиальном случае полученные выражения (17), (18) сводятся к выражению (20), которое совпадает с уравнением (27) в [9], а (28) сводится к (30), которое совпадает с уравнением (32) в [9]. Численное моделирование показало, что ближе к оптической оси параксиальные формулы (20) и (30) дают распределение амплитуды, близкое к значениям, даваемым интегралами (1), однако при удалении от оси полученные непараксиальные формулы намного точнее. Дифракционное фокусное пятно по полуспаду интенсивности для рассмотренного случая (рис. 2в) FWHM » 2,21, хотя числовая апертура линзы NA » + 52 » 0,196. Согласно скалярной теории диаметр фокусного пятна (по полуспаду интенсивности) должен быть равен FWHM » 0,51/^ . В нашем случае он должен быть равен FWHM » 2,51.
7 р,1хт
Рис. 1. Амплитуда радиально поляризованного гауссового пучка в плоскости геометрического фокуса линзы: радиальная составляющая Ег (а), продольная составляющая Ег (б) и интенсивность I (в)
Рис. 2. Амплитуда радиально поляризованного вихревого гауссового пучка (n = 1) в плоскости геометрического
фокуса линзы: радиальная составляющая Er (а), продольная составляющая Ez (б) и интенсивность I (в)
Благодарности Работа выполнена при поддержке российско-американской программы «Фундаментальные исследования и высшее образование» (грант CRDF PG08-014-1), гранта Президента РФ поддержки ведущих научных школ (НШ-3086.2008.9).
Литература
1. Bomzon, Z. Formation of radially and azimuthally polarized light using spacevariant subwavelength metal strip grating / Z. Bomzon, V. Kleiner, E. Hasman // Appl. Phys. Lett. 2001. - Vol. 79. - p. 1587-1589.
2. Armstrong, D.J. Generation of radially polarized beams with an image rotating resonator / D.J. Armstrong, M.C. Philips, A.V. Smith // Appl. Opt. 2003. - Vol. 42. - p. 3550-3554.
3. Passilly, N. Simple interferometric technique for generation of a radially polarized light beam / N. Passilly [and other] // J. Opt. Soc. Am. A 2005 22. - p. 984-991.
4. Volpe, G. Generation of cylindrical vector beams with few-mode fibers by Laguerre-Gaussian beams / G. Volpe, D. Petrov // Opt. Commun. 2004. - Vol. 237. - p. 89-95.
5. Niziev, V.G. Influence of beam polarization on laser cutting efficiency / V.G. Niziev, A.V. Nesterov // J. Phys. D 1999. - Vol. 32. - p. 1455-1461.
6. Gahagan, K.T. Simultaneous trapping of low-index and high-index microparticles observed with an optical-vortex trap / K.T. Gahagan, G.A. Swartzlander, Jr. // J. Opt Soc. Am. B 1999. - Vol. 16. - p. 533-537.
7. Wu, G. Focal shift in focused radially polarized ultrashort pulsed laser beams / G. Wu [and other] // Appl. Opt. 2007. - Vol. 46. - p. 6251-6255.
8. Deng, D. Nonparaxial propagation of radially polarized light beams / D. Deng // J. Opt. Soc. Am. B 2006. - Vol. 23. - p. 1228-1234.
9. Banerjee, P. A q-parameter approach to analysis of propagation, focusing, and waveguiding of radially polarized Gaussian beams / P. Banerjee, G. Cook, and D. Evans // J. Opt Soc. Am. A 2009. - Vol. 26(6) . - p. 1366-1374.
10. Kotlyar, V.V. Generation of phase singularity through diffracting a plane or Gaussian beam by a spiral phase plate / V. Kotlyar [and other] // J. Opt. Soc. Am. A 2005. -Vol. 22. - p. 849-861.
11. Котляр, В.В. Дифракция конической волны и гауссового пучка на спиральной фазовой пластинке / В.В. Котляр [и др.] // Компьютерная оптика, 2005. - № 28. -
C. 29-36. - ISSN 0134-2452.
12. Котляр, В.В. Дифракция Гауссового пучка на спиральном аксиконе / В.В. Котляр [и др.] // Компьютерная оптика, 2006. - № 30. - С. 30-35. - ISSN 0134-2452.
13. Ковалев, А.А. Параксиальные гипергеометрические лазерные пучки с особенностью в центре перетяжки / А.А. Ковалев [и др.] // Компьютерная оптика, 2007. -№ 31(1). - С. 9-13. - ISSN 0134-2452.
14. Ковалев, А.А. Непараксиальная векторная дифракция гауссового пучка на спиральной фазовой пластинке / А.А. Ковалев, В.В. Котляр // Компьютерная оптика, 2007. - № 31(4). - С. 19-22. - ISSN 0134-2452.
15. Котляр, В.В. Непараксиальные гипергеометрические моды / В.В. Котляр, А.А. Ковалев // Компьютерная оптика, 2008. - № 32(3). - С. 222-225. - ISSN 0134-2452.
16. Rashid, M. Focusing of high order cylindrical vector beams / M Rashid, O M Marago, P H Jones // J. Opt. A: Pure Appl. Opt. 2009. - Vol. 11. - p. 065204 (7pp).
17. Chen, B. Tight focusing of elliptically polarized vortex beams / B. Chen, J. Pu // Appl. Opt 2009. - Vol. 48. - p. 1288-1294.
18. Table of Integrals, Series, and Products / I.S. Gradshteyn and I.M. Ryzhik.- Elsevier, 2007.
19. Handbook of Mathematical Functions / ed. M. Abramowitz, I.A. Stegun.- National Bureau of Standards, Washington, DC, 1964.- 1044p.
References
1. Bomzon, Z. Formation of radially and azimuthally polarized light using spacevariant subwavelength metal strip grating / Z. Bomzon, V. Kleiner, E. Hasman // Appl. Phys. Lett. 2001. - Vol. 79. - p. 1587-1589.
2. Armstrong, D.J. Generation of radially polarized beams with an image rotating resonator / D.J. Armstrong, M.C. Philips, A.V. Smith // Appl. Opt. 2003. - Vol. 42. - p. 3550-3554.
3. Passilly, N. Simple interferometric technique for generation of a radially polarized light beam / N. Passilly [and other] // J. Opt. Soc. Am. A 2005 22. - p. 984-991.
4. Volpe, G. Generation of cylindrical vector beams with few-mode fibers by Laguerre-Gaussian beams / G. Volpe,
D. Petrov // Opt. Commun. 2004. - Vol. 237. - p. 89-95.
5. Niziev, V.G. Influence of beam polarization on laser cutting efficiency / V.G. Niziev, A.V. Nesterov // J. Phys. D 1999. - Vol. 32. - p. 1455-1461.
6. Gahagan, K.T. Simultaneous trapping of low-index and high-index microparticles observed with an optical-vortex trap / K.T. Gahagan, G.A. Swartzlander, Jr. // J. Opt. Soc. Am. B 1999. - Vol. 16. - p. 533-537.
7. Wu, G. Focal shift in focused radially polarized ultrashort pulsed laser beams / G. Wu [and other] // Appl. Opt. 2007. - Vol. 46. - p. 6251-6255.
8. Deng, D. Nonparaxial propagation of radially polarized light beams / D. Deng // J. Opt. Soc. Am. B 2006. - Vol. 23. - p. 1228-1234.
9. Banerjee, P. A q-parameter approach to analysis of propagation, focusing, and waveguiding of radially polarized Gaussian beams / P. Banerjee, G. Cook, and D. Evans // J. Opt. Soc. Am. A 2009. - Vol. 26(6) . - p. 1366-1374.
10. Kotlyar, V.V. Generation of phase singularity through diffracting a plane or Gaussian beam by a spiral phase plate / V. Kotlyar [and other] // J. Opt. Soc. Am. A 2005. -Vol. 22. - p. 849-861.
11. Kotlyar, V.V. Diffraction of conic wave and Gaussian beam by a spiral phase plate / V.V. Kotlyar [and other] // Computer Optics. - 2005. - Vol. 28. - P. 29-36. - ISSN 0134-2452. - (in Russian).
12. Kotlyar, V.V. Diffraction of Gaussian beam by a spiral axicon / V.V. Kotlyar [and other] // Computer Optics. -
2006. - Vol. 30. - P. 30-35. - ISSN 0134-2452. - (in Russian).
13. Kovalev, A.A. Paraxial hypergeometric laser beams with singularity at the center of the waist / A.A. Kovalev [and other] // Computer Optics. - 2007. - Vol. 31(1). - P. 9-13. - ISSN 0134-2452. - (in Russian).
14. Kovalev, A.A. Nonparaxial vertorial diffraction of Gaussian beam by a spiral phase plate / A.A. Kovalev, V.V. Kotlyar // Computer Optics. - 2007. - Vol. 31(4). - P. 1922. - ISSN 0134-2452. - (in Russian).
15. Kotlyar, V.V. Nonparaxial hypergeometric modes / V.V. Kotlyar, A.A. Kovalev // Computer Optics. - 2008. - Vol. 32(3). - P. 222-225. - ISSN 0134-2452. - (in Russian).
16. Rashid, M. Focusing of high order cylindrical vector beams / M Rashid, O M Marago, P H Jones // J. Opt. A: Pure Appl. Opt. 2009. - Vol. 11. - P. 065204 (7pp).
17. Chen, B. Tight focusing of elliptically polarized vortex beams / B. Chen, J. Pu // Appl Opt 2009. - Vol. 48. - p. 1288-1294.
18. Table of Integrals, Series, and Products / I.S. Gradshteyn and I.M. Ryzhik. - Elsevier, 2007.
19. Handbook of Mathematical Functions / ed. M. Abramowitz, I.A. Stegun. - National Bureau of Standards, Washington, DC, 1964.- 1044p.
NONPARAXIAL PROPAGATION OF VECTORIAL GAUSSIAN OPTCAL VORTEX
WITH RADIAL POLARIZATION
Alexey Andreevich Kovalev1'2 (researcher, e-mail: alanko@smr.ru), Victor Victorovich Kotlyar1'2 (head of laboratory, e-mail: kotlyar@smr.ru) 1 S.P. Korolyov Samara State Aerospace University, 2 Image Processing Systems Institute of the RAS
Abstract
Nonparaxial diffraction of Gaussian optical vortexes with initial radial polarization has been investigated. We have obtained explicit analytical expressions for radial, azimuthal, and longitudinal components of the electric vector. The simulation showed that nonparaxial formulae are more accurate than paraxial ones for areas located far from optical axis. The focal spot (area with intensity exceeding half of its maximal value) of focused Gaussian vortex field with unit topological charge and radial polarization appears to be less than that predicted by scalar theory.
Key words: nonparaxial diffraction, radial polarization, Gaussian optical vortex, focusing of laser light.
В редакцию поступила 11.09.2009г.