Вращающиеся элегантные пучки Бесселя-Гаусса

ДИФРАКЦИОННАЯ ОПТИКА, ОПТИЧЕСКИЕ ТЕХНОЛОГИИ

ВРАЩАЮЩИЕСЯ ЭЛЕГАНТНЫЕ ПУЧКИ БЕССЕЛЯ-ГАУССА

Котляр В.В., Ковалёв А.А., Скиданов Р.В., Сойфер ВА.

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

Аннотация

Рассмотрено новое трёхпараметрическое семейство вращающихся асимметричных пучков Бесселя-Гаусса (аБГ-пучки) с целым и дробным орбитальным угловым моментом (ОУМ). аБГ-пучки описываются произведением функции Гаусса и функции Бесселя первого рода n-го порядка с комплексным аргументом и имеют конечную энергию. Степень асимметрии аБГ-пучка зависит от вещественного параметра c > 0: при с = 0 аБГ-пучок совпадает с обычным радиально-симметричным пучком Бесселя-Гаусса; с ростом c аБГ-пучок приобретает форму полумесяца и при c >> 1 вытягивается по вертикальной оси и смещается по горизонтальной оси. Распределение интенсивности асимметричных пучков Бесселя-Гаусса в начальной плоскости имеет счётное число изолированных нулей, расположенных на горизонтальной оси. На месте этих нулей находятся оптические вихри с единичными топологическими зарядами и противоположными знаками с разных сторон от начала координат. При распространении пучка центры этих вихрей вместе со всем пучком вращаются вокруг оптической оси с неравномерной скоростью (при большом с >> 1): на расстоянии, равном длине Рэлея, они повернутся на 45 градусов, и на остальном расстоянии - ещё на 45 градусов. При разных значениях параметра c нули интенсивности в поперечном распределении интенсивности пучка меняют свои местоположения и изменяют ОУМ пучка. Изолированный ноль интенсивности на оптической оси порождает оптический вихрь с топологическим зарядом n. Лазерный пучок в виде вращающегося полумесяца был сформирован с помощью модулятора света.

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

Введение

В 1987 году Гори [1] рассмотрел пучки Бесселя-Гаусса (БГ-пучки). Комплексная амплитуда таких пучков описывается произведением Гауссовой функции на функцию Бесселя первого рода п-го порядка и на фазовую функцию, описывающую угловую гармонику. Комплексная амплитуда БГ-пучка удовлетворяет параксиальному уравнению распространения. Эти пучки имеют радиально-симметричное распределение интенсивности и обладают орбитальным угловым моментом (ОУМ). Радиальная симметрия сохраняется при распространении пучка. Но БГ-пучки не являются модами свободного пространства, так как при распространении не только масштабно изменяются, но также перераспределяется световая энергия между различными кольцами в распределении интенсивности в поперечном сечении пучка. БГ-пучки обобщались в ряде работ [2 - 4]. Например, в [4] рассмотрены пучки Гельмгольца-Гаусса, частным случаем которых являются пучки Бесселя-Гаусса [1]. БГ-пучки имеют конечную энергию, но известны также моды Бесселя, не обладающие конечной энергией [5]. Моды Бесселя являются решением уравнения Гельмгольца [6] и при распространении в однородном пространстве сохраняют свою интенсивность, поэтому также называются бездифракционными пучками Бесселя [5]. Линейная комбинация мод Бесселя с произвольными коэффициентами также является решением уравнения Гельмгольца. В [7] предложен алгоритм

расчёта фазового оптического элемента, который формирует бездифракционные пучки Бесселя с заданным модовым составом. В [8] предложено рассматривать пучок Матье как альтернативу пучку Бесселя. В [9] показано, что линейная комбинация чётного и нечётного пучков Матье с комплексными коэффициентами обладает ОУМ. Но сами пучки Матье не обладают ОУМ. Интересно [10], что линейная комбинация двух мод Эр-мита-Гаусса с комплексными коэффициентами, которые не обладают ОУМ, обладает ОУМ. Периодические функции Матье можно разложить в ряд Фурье [6]. Например, чётные функции Матье раскладываются по косинусам от полярного угла в цилиндрической системе координат, а нечётные - по синусам. Поэтому бездифракционный пучок Матье можно представить как линейную комбинацию мод Бесселя. Такие пучки рассматривались в [11].

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

зано, что такой асимметричный пучок Бесселя-Гаусса (аБГ-пучок) имеет в начальной плоскости счётное число изолированных нулей интенсивности, расположенных на горизонтальной оси. Все эти нули (кроме осевого) порождают оптические вихри с единичным топологическим зарядом и противоположными знаками с разных сторон от начала координат. Нуль интенсивности на оптической оси порождает оптический вихрь с топологическим зарядом п. При распространении в свободном пространстве аБГ-пучки вращаются вокруг оптической оси. Показано, что такие пучки имеют ОУМ, который растёт с номером п, а также увеличивается с ростом параметра с асимметрии пучка. Причём ОУМ аБГ-пучков может быть как целым, так и дробным.

1. Линейная комбинация БГ-пучков

Запишем комплексную амплитуду БГ-пучка [1] в начальной плоскости г = 0:

( 2 Л

Еп (г, j, z = 0) = exp

г

--2 + inj

ю„

Jn (ar) >

где а = k sin 90 = (2n/^)sin 90 - масштабирующий множитель, ш0 - радиус перетяжки Гауссова пучка, k = 2n/X - волновое число света с длиной волны X, 90 -угол конической волны, формирующей пучок Бесселя. В любой другой плоскости z комплексная амплитуда (1) будет иметь вид:

En (г, j, z) = q-1 (z) exp

f

ikz -

i'a2 z ^

2kq

x exp

o q

+ inj

J„

ar

qz)

где д(г) = 1 + и/го, г0 = М^/О - длина Рэлея, ]п(х) -функция Бесселя первого рода п-го порядка. Пучки (0) не являются параксиальными модами свободного пространства, так как аргумент функции Бесселя комплексный. Рассмотрим следующую суперпозицию БГ-пучков ц = д(г):

Ев (г, j, z; с) = q 1 exp

ikz -

v

. 2

ia z

2kq

2

qw

- cp exp (in j+ipj) г i ar

Ii Jn+ p I _

p=o p! v q

Поле (3) формирует параксиальный асимметричный пучок Бесселя-Гаусса (аБГ-пучок) при любых целых n и любой комплексной постоянной с. Но для простоты в дальнейшем будем считать постоянную c вещественной положительной величиной: c > 0. Хотя если рассмотреть c как комплексную величину c = |c| arg c или учесть, что она может быть отрицательной c<0, то распределение интенсивности поля (3) повернётся на угол arg c вокруг оптической оси. При c = 0 в (3) отличным от нуля остаётся только одно слагаемое при p = 0, и аБГ-пучок становится обычным БГ-пучком (1), (2).

Здесь рассматривается линейная комбинация БГ-пучков в виде (3) потому, что этот ряд равен функции Бесселя с комплексным аргументом. Действительно, в [10] есть справочный ряд (выражение 5.7.6.1):

Ё ki Jk +v ( *) =

k=0 k !

= ^2

(* - 2í У2 Jv [JXF-lX).

С помощью (4) преобразуем (3) и получим Í ,„.2 „2 Л

En (г, j, z; с) = —exp

ikz -

aг

v

B 2

---- + inj

2kq qWo

aг - 2cq exp (ij) xJn {q ^aг[aг- 2cqexp (ij)]}.

Выражение (5) является замкнутой формой для комплексной амплитуды трёхпараметрического семейства параксиальных скалярных аБГ-пучков. Два непрерывных вещественных параметра аБГ-пучков управляют масштабом (а) и степенью асимметрии (с). Радиус Гауссова пучка ш0 считается одинаковым для всего семейства аБГ-пучков. Целый параметр п вместе с непрерывным с определяют величину ОУМ. Когда знаменатель в (5) равен нулю, тогда и аргумент функции Бесселя равен нулю. Неопределённость ноль на ноль раскрывается.

Пучок (5) имеет счётное число изолированных нулей интенсивности, которые порождают оптические вихри с единичным топологическим зарядом, кроме осевого нуля, который порождает оптический вихрь с п-м топологическим зарядом. Для получения полярных координат нулей аБГ-пучка приравняем аргумент функции Бесселя в (5) корню функции Бесселя Ър (Л (ЪР) = 0):

а2г2 - Оасцгехр

Выделим вещественную и мнимую части этого уравнения:

a2г2 - 2aс|q|гcos(j+Y) = y2p|q| cos(2Y), l^ac^T sin (j+Y) = y2p|q|2 sin (2Y),

где Y = атс^г^) - фаза Гоу. Из (7) находятся координаты точек с нулевой интенсивностью:

„0

1

jnp = 2arccos

cos

г

np

= qVgBp cos(2Y) + 2c2 ± 24d,

Из обоих

где D =

уравнений (8) следует: чтобы координаты нулей были вещественными, необходимо выполнение условия:

gp sin Y < с.

x

x

2

2

X

Из (9) следует, что если с> упр, то все нули интенсивности с номерами от д = 0 до д=р при распространении не исчезают (так как для этих нулей упд < с и неравенство (9) выполняется для любых расстояний 2), однако остальные нули с номерами д=р + 1, р + 2, ... (для которых с< упд) пропадают на некотором расстоянии 2 = С20 / (упд2 - с2)1/2.

При 2 = 0 вместо (8) запишем:

Гр+=а~1 (с +у! С2 +у2р), Фпр = 2 рк, гР-=а~1 (Vс2 +У2пр - с), Фпр =(2р+1)к.

Из (10) следует, что интенсивность в сечении аБГ-пучка несимметрична относительно начала координат, так как гр+ > гр-. Причём, когда с > 0 растёт, асимметрия аБГ-пучка увеличивается, а при с = 0 изолированные нули аБГ-пучка пропадают, а появляются кольца нулевой интенсивности радиально-симмет-ричного пучка Бесселя-Гаусса. Все нули аБГ-пучка (кроме осевого при г = 0) порождают оптические вихри с единичным топологическим зарядом и противоположными знаками с разных сторон от начала координат. Для нулей интенсивности с радиальными координатами гр+ топологический заряд оптических вихрей +1, для гр_ - топологический заряд -1. Чтобы сменить знаки этих оптических вихрей на противоположные, надо вместо (5) взять комплексно сопряжённое выражение. Из (8) видно, что нули интенсивности, а значит, и само распределение интенсивности в поперечной плоскости аБГ-пучка вращаются при распространении. Но вращение это достаточно сложное. Только когда с>> 1, в первом уравнении в (8) можно пренебречь вторым слагаемым в квадратных скобках, и тогда картина интенсивности в поперечном сечении пучка вращается как единое целое. При этом полярный угол меняется с расстоянием 2

Ф = агс12

Из (11) следует, что аБГ-пучок от начальной плоскости 2 = 0 до расстояния Рэлея 2 = 20 повернётся против часовой стрелки на л/4, а на остальном участке пути от z = 20 до z = го - ещё на п/4, и, таким образом, на всём протяжении пути пучок повернётся на п/2. Вращение аБГ-пучка не зависит от номера п, то есть вращаться будет и нулевой пучок при п = 0.

На рис. 1 показаны распределение интенсивности и фазы аБГ-пучка (5) при разных значениях параметра асимметрии с. При моделировании были использованы следующие значения параметров: длина волны 1 = 532 нм, радиус перетяжки Гауссова пучка Щ =101, масштабирующий множитель а = 1/(101), границы расчётной области -401 < х, у< 401. На рис. 1 видно, что с ростом с > 0 интенсивность из кольцевой становится похожей на полумесяц, выпуклый в правую сторону. Можно записать выражение для аБГ-пучка, имеющего вид полумесяца, выпуклого в левую сторону. аБГ-пучок, имеющий распределение интенсивно-

сти, зеркально-симметричное относительно оси х = 0 по сравнению с пучком (5), имеет амплитуду вида:

0>2

Еп (г, Ф, 2 = 0; с) = ехр

х

(-с) р ехр[;(п + р)ф] л + р (аг) =

р=0 р'

, г2

= еХР| --2

Щ

'0 у

аг

аг + 2 с ехр(ф) х Лп |^аг [аг+2с ехрф)]} ехр(шф).

Рис. 1. Распределение интенсивности

Также на рис. 1 видно, что нули интенсивности, лежащие на оси х< 0, с ростом с приближаются к центру координат: на рис. 1б, г их ещё не видно, а на рис. 1е появился один ноль. Такое поведение нулей следует из второго уравнения (10) при больших с >> 1. В центре рис. 1б, г, е расположен нуль интенсивности третьего порядка (п = 3).

На рис. 2 показаны распределения интенсивности и фазы аБГ-пучка третьего порядка (п = 3) при относительно большом значении параметра с = 10 на разных расстояниях вдоль оптической оси в плоскостях 2 = 0 (а, б), 2 = 20 (в, г), 2 = 1020 (д, е). На рис. 2 видно, что пучок вращается вокруг оптической оси. На расстоянии Рэлея 20 картина интенсивности повернулась против часовой стрелки на 45 градусов (рис. 2в), а на

п /2

х

0

расстоянии 10^о картина интенсивности повернулась почти на 90 градусов (рис. 2д).

Рис. 2. Распределение интенсивности

Размер картинок на рис. 2 разный: -801 < x, y< 801 (а, б); -1001 < x, y< 1001 (в, г) и -5001 < x, y< 5001 (д, е). Остальные параметры те же, что и для рис. 1. На рис. 2 также видно, что нули интенсивности, лежащие на оси x при z = 0, тоже начинают вращаться против часовой стрелки c увеличением z. На рис. 2б видны три изолированных нуля интенсивности на оси x слева от центрального нуля 3-го порядка. На рис. 2г видно, что при z = 2q осталось только два нуля, которые повернулись на 45 градусов, а на рис. 2е (z = 10zô) эти два нуля объединились в один ноль 2-го порядка и повернулись почти на 90 градусов.

Из неравенства (9) следует, что, так как gn1 > 1 для всех целых n, при c < 1 все изолированные нули интенсивности, лежащие на оси x (кроме центрального нуля), будут «пропадать» по мере распространения пучка, начиная с самых дальних нулей (yp >> 1). Последним «пропадёт» первый ноль интенсивности y1 при z, равном:

z = zn tg

arcsin

Y,

Нули интенсивности «пропадают» в том смысле, что значения их координат, согласно (8) и при условии (9), становятся комплексными, а не действительными.

На рис. 3 показаны распределения интенсивности и фазы аБГ-пучка при небольшом параметре асимметрии с = 1. Остальные параметры расчёта: длина волны 1=532 нм, радиус перетяжки Гауссова пучка Щ = 1001,

масштабирующий множитель а = 1/(101), границы расчётной области -3001 < х, у< 3001 (рис. 3а, б) и - 10001 < х, у< 10001 (рис. 3в, г). При этих параметрах вращение аБГ-пучка уже не описывается простой формулой (11), так как параметр асимметрии с небольшой. На рис. 3 видно, что картина интенсивности поворачивается на расстоянии 2= 20 почти на 90 градусов против часовой стрелки. А изолированные нули интенсивности, которые лежат на оси х и видны на рис. 3б, при распространении пучка «пропадают»: на рис. 3г при 2= 20 не видно нулей интенсивности, кроме центрального.

Рис. 3. Распределение интенсивности

На рис. 4 показаны интенсивность и фаза аБГ-пучка нулевого порядка при следующих параметрах: длина волны 1=532 нм, радиус перетяжки Гауссова пучка Щ) = 101, масштабирующий множитель а = 1/(101), параметр асимметрии с = 10, границы расчётной области -401 < х, у < 401.

г)

Рис. 4. Интенсивность

В случае нулевого порядка аБГ-пучок (рис. 4) обладает интересным свойством: имеет максимум интенсивности вблизи оптической оси и ОУМ, отличный от нуля. Это свойство можно использовать при манипулировании диэлектрическими микрочастицами. Частица, которая по размерам в несколько раз больше основного максимума интенсивности пучка на рис. 4, может удерживаться этим максимумом интенсивности и одновременно вращаться вокруг своей оси.

2. Фурье-спектр аБГ-пучка

Угловой спектр обычного БГ-пучка известен [4]:

¿рю0

An (р, f) = (-iГ exp

2 f

xI„

k f j exp (inf),

где Iп (x) - модифицированная функция Бесселя п-го порядка, р - радиальная координата в Фурье-плоскости, /фокусное расстояние сферической линзы, формирующей пространственный спектр БГ-пучка. Рассмотрим линейную комбинацию, аналогичную (3), но состоящую из функций (14), получим угловой спектр аБГ-пучка:

ч2~

кр«0

An (р, f) = exp

2 f

х£

p=0

( ic)"exP[i(n + p)f] fkарю2 j

p!

2f

Используя справочное соотношение [12]:

£ к! Ik +v ( =

k=0 k !

= xV2 (x + 2"у2 L иxf + 2"x |

,(V x2 + 2x),

получим окончательное выражение для Фурье-спектра аБГ-пучка:

An ^ f) = exP

J x

kPW 2 f

х Т

+ inf

yx [X + 2c ei(f-P/2) ] },

^ х+ 2се1'(ф-р/2), где Х = акрю2/(2^). Угловой спектр (17) асимметричный: при ф = п/2 и на кольце фиксированного радиуса р = р0 модуль амплитуды (17) имеет максимум, а при ф = -п/2 - минимум. Вид углового спектра (17) похож на вид амплитуды аБГ-пучка (5) (с точностью до замены функции Бесселя модифицированной функцией Бесселя), но повёрнут на 90 градусов. Поэтому асимметрия спектра (17) будет похожа на асимметрию аБГ-пучка (5).

3. Орбитальный угловой момент аБГ-пучка

Орбитальный угловой момент ]г (проекция ОУМ на оптическую ось) и суммарная интенсивность I светового пучка в плоскости, поперечной оптической оси, определяются по формулам [6]:

Jz = Im<j ff E ^drdj № Эф

I = JJ KEr dr dj.

Подставив в (18), (19) комплексную амплитуду (3) при г = 0, получим орбитальный угловой момент и суммарную интенсивность аБГ-пучка:

^2 p

Jz=2*£ ( л

p=o

(n + p).

xJ exp

0

I = 2p£

2r

2

w

J+ p (ar)rdr,

c

p=o

( p !)2

J exp

2r2

w

J2+ p (ar)rdr

o;

Интегралы в этих выражениях описаны в [12]: J xexp (-px2) Jv (bx) Jv (cx) dx =

= ( 2 p )-1

exp

b2 + c2j

4 p

I' I bp

Используя интеграл (22) и разделив (20) на (21), получим выражение для орбитального углового момента, нормированного на интенсивность:

-1

J=n+£

1 p=0

>4+ p (7) ( p!)2

- c2pi

£ n+ p

p=0

где 7 =a2o>2/4. Выражение (23) дальше упростить не удаётся. Из (23) следует, что ОУМ аБГ-пучков больше n, так как все слагаемые рядов в (23) положительные. То есть с ростом параметра с растёт асимметрия аБГ-пучка и почти линейно растёт его ОУМ, как видно из рис. 5. Так как параметры а и c являются действительными положительными числами, то второе слагаемое в (23) может быть как целым, так и дробным положительным числом. Поэтому из (23) также следует, что аБГ-пучок нулевого порядка (n = 0) может иметь любой ОУМ. Заметим также, что в (23) масштабные параметры пучка Бесселя а и Гауссова пучка ш0 входят в виде произведения, поэтому у разных аБГ-пучков, у которых аш0 = const, ОУМ будет одинаковым (при равных n и c).

Jz, отн. ед.

6

4 8 12 16 20

Рис. 5. Зависимость ОУМ от параметра асимметрии c при n = 3, ю0 = 10Я, а =1/

2

0

0

p

2

х

4. Взаимная неортогональность функций, описывающих семейство аБГ-пучков

Подобно тому, как был вычислен орбитальный угловой момент, можно вычислить скалярное произведение комплексных амплитуд двух ЭБ-пучков -пучка п-го порядка с параметрами а и с и пучка т-го порядка с параметрами (3 и ±

(

х exp

^nßd

0 (a2 +

ß2)

(cd *)p+

\n - m

p=0 p!(p+\n-m)!

p + max

w2 aß 4

где In(x) - модифицированная функция Бесселя. Из выражения (24) видно, что в отличие от мод Бесселя [5] аБГ-пучки не ортогональны ни по масштабирующему множителю, ни по порядку функции Бесселя, ни по параметру асимметрии. В (24) параметры асимметрии c и d считались комплексными, а «*» -комплексное сопряжение.

Заметим, что аБГ-пучки при стремлении радиуса Гауссова пучка к бесконечности w0 ^ да переходят в бездифракционные асимметричные элегантные моды Бесселя [13, 14].

5. Эксперимент

Схема экспериментальной установки представлена на рис. 6. В эксперименте был использован фазовый пространственный модулятор PLUTO VIS. Свет твердотельного лазера с длиной волны 532 нм расширяется коллиматором, ограничивается диафрагмой диаметром 8 мм. В результате получается равномерное распределение интенсивности, которое можно считать плоской волной. Далее свет проходит через светоделительный кубик и попадает на пространственный модулятор света, отражается от него и отклоняется светоделитель-ным кубиком на ПЗС-камеру. На полутоновой модулятор с компьютера передаётся фаза для формирования аБГ-пучка в комбинации с фазой параболической линзы с фокусным расстоянием 960 мм.

Рис. 6. Схема экспериментальной установки.

1 - твердотельный лазер с длиной волны 532нм, 2, 3 - коллиматор, 4 - диафрагма, 5 - светоделительный кубик, 6 - модулятор PLUTO VIS, 7 - зеркало, 9 - ПЗС-камера

ПЗС-камера перемещается на небольшом отрезке в непосредственной близости от фокуса. На рис. 7 показаны фазовое распределение (а), сформированное на

модуляторе (без добавленной фазы линзы), и распределения интенсивности, зарегистрированные ПЗС-матрицей на расстоянии 850 мм (б), 900 мм (в) и 950 мм (г) от модулятора. Размерность модулятора составляет 1920*1080 пикселов, размер одного чувствительного элемента -8 мкм. Фаза, представленная на рис. 7, имела размерность 1024*1024 пикселов и формировалась в центре модулятора. Таким образом, точный размер сформированного фазового распределения составлял 8,2 мм.

Рис. 7. Распределение фазы, сформированной на модуляторе PLUTO VIS

На рис. 7 видно, что сформированное модулятором распределение интенсивности в форме полумесяца не только вращается в сходящемся лазерном пучке (за 100 мм пути пучок повернулся почти на 90 градусов против часовой стрелки), но и искажается. Это искажением связано с тем, что фазовое распределение (рис. 7 а) только частично учитывает амплитудное распределение аБГ-пучка (5).

Заключение

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

- получено новое решение параксиального уравнения Гельмгольца, описывающее трёхпараметриче-ское семейство асимметричных пучков Бесселя-Гаусса (аБГ-пучков), обладающих конечной энергией и описываемых функциями Бесселя первого рода целого порядка с комплексным аргументом; при с = 0 аБГ-пучки совпадают с обычными пучками Бесселя-Гаусса [1];

- распределение интенсивности аБГ-пучков имеет счётное число изолированных нулей, лежащих на горизонтальной прямой в начальной плоскости (z=y = 0); в местах этих нулей расположены оптические вихри, имеющие единичные топологические заряды и разные знаки с разных сторон от оптической оси (от центра координат); осевой ноль интенсивности имеет топологический заряд, равный порядку функции Бесселя;

n-m

/ л

2

х

х

2

- аБГ-пучки при больших параметрах асимметрии c >> 1 вращаются как единое целое при распространении в свободном пространстве: на расстоянии длины Рэлея они поворачиваются на 45° и ещё на 45° - за весь остальной путь; при небольших c пучки вращаются более сложным образом: на определённом расстоянии от начальной плоскости оптические вихри (изолированные нули интенсивности), которые присутствовали в сечении пучка, пропадают (кроме центрального нуля интенсивности) и на расстоянии длины Рэлея пучок может повернуться почти на 90°;

- аБГ-пучки имеют асимметричный угловой спектр, выражающийся через произведение Гауссовой функции на модифицированную функцию Бесселя с комплексным аргументом; вид этого спектра имеет асимметрию, аналогичную асимметрии аБГ-пучка, но повёрнутую на 90°;

- аБГ-пучки могут иметь целый и дробный ОУМ, который растёт с ростом номера п и почти линейно с ростом параметра c; пучок нулевого порядка имеет максимум интенсивности, смещённый с оптической оси на величину с/a, изолированные нули интенсивности и может обладать любым ОУМ в зависимости от выбора параметра асимметрии с;

- аБГ-пучки не ортогональны ни по одному из параметров п, а, С;

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

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

Работа выполнена при поддержке грантов Президента РФ поддержки ведущих научных школ (НШ-3970.2014.9) и молодого доктора наук (МД-1929.2013.2), а также грантов РФФИ 13-07-97008 и 14-07-31092.

Литература

1. Gori, F. Bessel-Gauss beams / F. Gori, G. Guattari, C. Pa-dovani // Optics Communications. - 1987. - V. 64. - P. 491495.

2. Li, Y. New generalized Bessel-Gauss beams / Y. Li, H. Lee, E. Wolf // Journal of the Optical Society of America A. - 2004.

- V. 21. - P. 640-646.

3. Kisilev, A.P. New structures in paraxial Gaussian beams // Optics and Spectroscopy. - 2004. - V. 96. - P. 479-481.

4. Gutierrez-Vega, J.C. Helmholz-Gauss waves / J.C. Gutierrez-Vega, M.A. Bandres // Journal of the Optical Society of America A. - 2005. - V. 22. - P. 289-298.

5. Durnin, J. Exact solutions for nondiffracting beams. I. The scalar theory // Journal of the Optical Society of America A.

- 1987. - V. 4, Issue 4. - P. 651-654.

6. Miller Jr., W. Symmetry and separation of variables. - Ad-dison-Wesley Pub. Comp., 1977.

7. Kotlyar, V.V. Algorithm for the generation of non-diffracting Bessel beams / V.V. Kotlyar, S.N. Khonina, V.A. Soifer // Journal of Modern Optics. - 1992. - V. 42, Issue 6. - P. 1231-1239.

8. Gutierrez-Vega, J.C. Alternative formulation for invariant optical fields: Mathieu beams / J.C. Gutierrez-Vega, M.D. Itur-be-Castillo, S. Chavez-Cedra // Optics Letters. - 2000. -V. 25, Issue 20. - P. 1493-1495.

9. Chavez-Cedra, S. Elliptic vortices of electromagnetic wave fields / S. Chavez-Cedra, J.C. Gutierrez-Vega, G.H.C. New // Optics Letters. - 2001. - V. 26, Issue 22. - P. 1803-1805.

10. Kotlyar, V.V. Hermite-Gaussian modal laser beams with orbital angular momentum / V.V. Kotlyar, A.A. Kovalev // Journal of the Optical Society of America A. - 2014. -V. 31, Issue 2. - P. 274-282.

11. Dennis, M.R. Propagation-invariant beams with quantum pendulum spectra: from Bessel beams to Gaussian beam-beams / M.R. Dennis, J.D. Ring // Optics Letters. - 2013. -V. 38, Issue 17. - P. 3325-3328.

12. Прудников, А.П. Интегралы и ряды. Специальные функции / А.П. Прудников, Ю.А. Брычков, О.И. Мари-чев. - М.: Наука, 1983.

13. Котляр, В. В. Бездифракционные асимметричные элегантные пучки Бесселя с дробным орбитальным угловым моментом / В.В. Котляр, А.А. Ковалёв, В.А. Сойфер // Компьютерная оптика. - 2014. - Т. 38, № 1. - С. 4-10.

14. Kotlyar, V.V. Asymmetric Bessel modes / V.V. Kotlyar, A.A. Kovalev, V.A. Soifer // Optics Letters. - 2014. -V. 39, Issue 8. - P. 2395-2398.

References

1. Gori, F. Bessel-Gauss beams / F. Gori, G. Guattari, C. Padovani // Optics Communications. - 1987. - V. 64. - P. 491-495.

2. Li, Y. New generalized Bessel-Gauss beams / Y. Li, H. Lee, E. Wolf // Journal of the Optical Society of America A. - 2004.

- V. 21. - P. 640-646.

3. Kisilev, A.P. New structures in paraxial Gaussian beams // Optics and Spectroscopy. - 2004. - V. 96. - P. 479-481.

4. Gutierrez-Vega, J.C. Helmholz-Gauss waves / J.C. Gutierrez-Vega, M.A. Bandres // Journal of the Optical Society of America A. - 2005. - V. 22. - P. 289-298.

5. Durnin, J. Exact solutions for nondiffracting beams. I. The scalar theory // Journal of the Optical Society of America A.

- 1987. - V. 4, Issue 4. - P. 651-654.

6. Miller Jr., W. Symmetry and separation of variables. - Ad-dison-Wesley Pub. Comp., 1977.

7. Kotlyar, V.V. Algorithm for the generation of non-diffracting Bessel beams / V.V. Kotlyar, S.N. Khonina, V.A. Soifer // Journal of Modern Optics. - 1992. - V. 42, Issue 6. - P. 1231-1239.

8. Gutierrez-Vega, J.C. Alternative formulation for invariant optical fields: Mathieu beams / J.C. Gutierrez-Vega, M.D. Itur-be-Castillo, S. Chavez-Cedra // Optics Letters. - 2000. -V. 25, Issue 20. - P. 1493-1495.

9. Chavez-Cedra, S. Elliptic vortices of electromagnetic wave fields / S. Chavez-Cedra, J.C. Gutierrez-Vega, G.H.C. New // Optics Letters. - 2001. - V. 26, Issue 22. - P. 1803-1805.

10. Kotlyar, V.V. Hermite-Gaussian modal laser beams with orbital angular momentum / V.V. Kotlyar, A.A. Kovalev // Journal of the Optical Society of America A. - 2014. -V. 31, Issue 2. - P. 274-282.

11. Dennis, M.R. Propagation-invariant beams with quantum pendulum spectra: from Bessel beams to Gaussian beam-beams / M.R. Dennis, J.D. Ring // Optics Letters. - 2013. -V. 38, Issue 17. - P. 3325-3328.

12. Prudnikov, A.P. Integrals and series. Special functions / A.P. Prudnikov, J.A. Brychkov, O.I. Marichev. - Moscow: "Nauka" Publisher, 1983. - (In Russian).

13. Kotlyar, V.V. Diffraction-free asymmetric elegant Bessel beams with fractional orbital angular momentum / V.V. Kotlyar, A.A. Kovalev, V.A. Soifer // Computer Optics. - 2014. - V. 38(1). - P. 4-10.

14. Kotlyar, V.V. Asymmetric Bessel modes / V.V. Kotlyar, A.A. Kovalev, V.A. Soifer // Optics Letters. - 2014. -V. 39, Issue 8. - P. 2395-2398.

ROTATING ELEGANT BESSEL-GAUSSIAN BEAMS

V. V. Kotlyar, A.A. Kovalev, R. V. Skidanov, V.A. Soifer Image Processing Systems Institute, Russian Academy of Sciences, Samara State Aerospace University

Abstract

We considered a new three-parameter family of rotating asymmetric Bessel-Gaussian beams (aBG-beams) with integer and fractional orbital angular momentum (OAM). Complex amplitude of aBG-beams is proportional to a Gaussian function and to n-th order Bessel function with complex argument. These beams have finite energy. The asymmetry ratio of the aBG-beam depends on the real parameter c > 0. For c = 0, the aBG-beam coincides with the conventional radially symmetric Bessel-Gaussian beam, with increasing c intensity distribution of the aBG-beam takes the form of Crescent. At c >> 1, the beam is widening along the vertical axis while shifting along the horizontal axis. In initial plane, the intensity distribution of asymmetric Bessel-Gaussian beams has a countable number of isolated zeros located on the horizontal axis. Locations of these zeros correspond to optical vortices with unit topological charges and opposite signs on different sides of the origin. During the beam propagation, centers of these vortices rotate around the optical axis along with the entire beam at a non-uniform rate (for large c >> 1): they will rotate by 45 degrees at a distance of the Rayleigh range, and then by another 45 degrees during the rest distance. For different values of the parameter c zeros in the transverse intensity distribution of the beam change their location and change the OAM of the beam. Isolated intensity zero on the optical axis corresponds to an optical vortex with topological charge of n. Laser beam with the shape of rotating Crescent has been generated by using a spatial light modulator.

Key words: Bessel-Gaussian laser beam, orbital angular momentum, rotation of light beam, complex argument Bessel function.

Сведения об авторах

Котляр Виктор Викторович, 1957 года рождения, доктор физико-математических наук, профессор, заведующий лабораторией лазерных измерений ИСОИ РАН и по совместительству профессор кафедры технической кибернетики Самарского государственного аэрокосмического университета (СГАУ). В 1979 году окончил физический факультет Куйбышевского государственного университета, в 1988 году защитил кандидатскую диссертацию в Саратовском государственном университете, а в 1992 году - докторскую диссертацию в Центральном конструкторском бюро Уникального приборостроения РАН (г. Москва). Область научных интересов: нанофотоника, дифракционная компьютерная оптика. Публикации: 300 научных трудов, 5 монографий, 7 авторских свидетельств.

Victor Victorovich Kotlyar is a head of Laboratory at the Image Processing Systems Institute (Samara) of the Russian Academy of Sciences and professor of Computer Science department at Samara State Aerospace University. He received his MS, PhD and DrSc degrees in Physics and Mathematics from Samara State University (1979), Saratov State University (1988) and Moscow Central Design Institute of Unique Instrumentation, the Russian Academy of Sciences (1992). He is co-author of 300 scientific papers, 5 books and 7 inventions. His current interests are diffractive optics, gradient optics, nanophotonics, and optical vortices.

Ковалёв Алексей Андреевич, 1979 года рождения, в 2002 году окончил Самарский государственный аэрокосмический университет имени академика С.П. Королёва - СГАУ по специальности «Прикладная математика». Доктор физико-математических наук (2012 год), работает старшим научным сотрудником лаборатории лазерных измерений Института систем обработки изображений РАН. В списке научных работ более 80 статей. Область научных интересов: математическая теория дифракции, сингулярная оптика, фо-тонно-кристаллические устройства. E-mail: alanko@smr.ru .

Alexey Andreevich Kovalev (b. 1979), graduated (2002) from S.P. Korolyov Samara State Aerospace University (SSAU)), majoring in Applied Mathematics. He received his Doctor in Physics & Maths degree (2012). He is a senior researcher of Laser Measurements laboratory at the Image Processing Systems Institute of the Russian Academy of Sciences (IPSI RAS), holding a part-time position of assistant at SSAU's Technical Cybernetics sub-department. He is co-author of more than 80 scientific papers. Research interests are mathematical diffraction theory, singular optics, and photonic crystal devices.

Скиданов Роман Васильевич, 1973 года рождения. В 1990 году с отличием окончил Самарский государственный университет (СамГУ) по специальности «Физика». Доктор физико-математических наук (2007 год), работает заведующим лабораторией микро- и нанотехнологий Института систем обработки изображений РАН (ИСОИ РАН), профессором кафедры технической кибернетики СГАУ. Скиданов Р.В. - специалист в области дифракционной оптики, математического моделирования, обработки изображений и нанофотоники. В списке научных работ Р.В. Скиданова 90 статей, 5 монографий. E-mail: romans@smr.ru .

Roman Vasilevich Skidanov (b. 1973). Graduated with honours (1990) from Samara State University (SSU)), majoring in Physics. He received his Doctor in Physics & Maths (2007) degrees from Samara State Aerospace University. He is the head of Micro- and Nanotechnologies laboratory at the Samara Image Processing Systems Institute of the Russian Academy of Sciences (IPSI RAS), holding a part-time position of professor at SSAU's Technical Cybernetics sub-department. He is co-author of 90 scientific papers, 5 monographs. His current research interests include diffractive optics, mathematical modeling, image processing, and nanophotonics.

Сойфер Виктор Александрович - член-корреспондент РАН, доктор технических наук, директор Института систем обработки изображений РАН. Лауреат государственной премии и премии правительства РФ в области науки и техники. Диплом инженера получил в Куйбышевском авиационном институте (ныне Самарский государственный аэрокосмический университет имени академика С.П. Королёва) в 1968 году. Докторскую диссертацию защитил в 1979 году в Санкт-Петербургском государственном электротехническом университете «ЛЭТИ». Автор и соавтор значительного числа научных публикаций с высоким индексом цитируемости, 10 книг и 50 авторских свидетельств и патентов. E-mail: soifer@ssau.ru .

Victor Alexandrovich Soifer is a corresponding member of the Russian Academy of Sciences (RAS), Doctor of Engineering, Professor, and director of the Image Processing Systems Institute (IPSI) of the RAS. He is the recipient of the State Prize of the Russian Federation and the RF government prize in Science and Technology. He was awarded an Engineer qualification (1968) by Kuibyshev Aviation Institute (KuAI, presently, Samara State Aerospace University named after academician S.P. Korolyov, SSAU). He received a Doctor in Engineering degree (1979) from Saint Petersburg State Electrotechnical University "LETI" (ETU). He is a SPIE- and IAPR-member. He is the author and coauthor of a substantial number of scientific publications with high citation index, 10 books, and 50 author's certificates and patents.

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