CC BY
0
Управление всплесками орбитального углового момента в астигматических структурированных пучках в рамках ABCD
матричных преобразований
А.В. Воляр1, М.В. Брецько1, С.И. Халилов1, Я.Е. Акимова1 1КФУ им. В.И. Вернадского, Физико-технический институт, 295007, Россия, Республика Крым, г. Симферополь, проспект Академика Вернадского, д. 4
Аннотация
Разработан и экспериментально воплощен ABCD матричный формализм для описания структурных преобразований и орбитального углового момента структурированных Ла-герр-Гауссовых пучков в астигматической оптической системе, содержащей цилиндрическую и корректирующую сферическую линзу. Показано, что матричный формализм не только хорошо согласуется с методом интегральных астигматических преобразований, но и значительно расширяет область его применения. Выявлено, что корректирующая сферическая линза способна не только разделять вихревую и астигматическую составляющие орбитального углового момента, но и превращать структурно неустойчивый астигматический Лагерр-Гауссов пучок после цилиндрической линзы в структурно устойчивый, при условии сохранения формы быстрых осцилляций и всплесков орбитального углового момента.
Ключевые слова: вихревые пучки, ABCD-матрица, структурированный свет, орбитальный угловой момент.
Цитирование: Воляр, А.В. Управление всплесками орбитального углового момента в астигматических структурированных пучках в рамках ABCD матричных преобразований / А.В. Воляр, М.В. Брецько, С.И. Халилов, Я.Е. Акимова // Компьютерная оптика. - 2024. -Т. 48, № 2. - С. 171-179. - DOI: 10.18287/2412-6179-CO-1411.
Citation: Volyar AV, Bretsko MV, Khalilov SI, Akimova YE. Orbital angular momentum burst control in astigmatic structured beams in ABCD-matrix transforms. Computer Optics 2024; 48(2): 171-179. DOI: 10.18287/2412-6179-CO-1411.
Введение
Наглядность и относительная математическая простота унитарных АВСБ матричных преобразований, открытых еще в начале 70-х годов прошлого века [1, 2], оказывается сейчас наиболее удобным инструментом для инженерного моделирования сложных оптических устройств и систем на основе структурированных когерентных и частично когерентных пучков [3 - 8]. Основным элементом в этом подходе является унитарная ЛВСБ-матрица вида
M =
'A B
. C D
где detM=AD - BC = 1 для отдельного элемента системы, так что вся оптическая система описывается последовательным произведением этих матриц. Как правило, такие матричные преобразования используются для вычисления комплексного параметра q(z) фундаментального Гауссова пучка, включающего в себя радиус кривизны волнового фронта R(z) и радиус поперечного сечения пучка w(z) как 1 / q(z) = 1 / R(z) - i2 / kw2(z), где k - волновое число. Однако, когда приходится работать с пучками Эрми-та-Гаусса (ЭГ-), Лагерра-Гаусса (ЛГ-) или другими типами пучков высших порядков, обращаются к интегралу Коллинза [9], записанному для когерентных параксиальных пучков как
Т(г,z) =-[ 2Т(р)ехр(-/'-[А(z) | р |2 -
2тВ(z) -1«2 2B(z) (1)
-2(рг) + D( ^|г|2]}^ 2р,
где ^(р) - комплексная амплитуда пучка на входе оптической системы, р = г = (х, у), A(z), B(z), С^), D(z) - элементы ЛБСБ-матрицы.
Используя этот интеграл, Сейгман [2], а затем Бе-ленджер [10] наметили общий подход для описания единичных ЭГ- и ЛГ-пучков в простейших оптических системах, содержащих сферические линзы. Аналогичный расчет для единичного ЛГ-пучка был выполнен Таше [11]. Алиева и Бастианс [12, 13], рассматривающие пучки ЭГ- и Эрмита-Лагерра-Гаусса (ЭЛГ-) в оптических системах первого рода с каскадом линз, усилителем и ортосимплектической системой (система, которая является одновременно упрощенной и ортогональной), используют подход на основе 4^4 симплектических матриц. Преобразование ЭГ-пучков в более сложных системах - конверторах оптических вихрей, содержащих астигматические элементы, авторы статей [14, 15] применили скорее искусственный, чем стандартный ЛБСБ-подход, в то время как Хонина и соавторы [16] рассмотрели различные типы астигматических модовых конвертеров, используя для моделирования стандартные дифракционные методы. Дело в том, что, если в оптической системе присутствуют астигматические элементы
(например, цилиндрические линзы), ее анализ значительно усложняется. Теперь любой цилиндрический пучок должен быть представлен в базисе собственных ЭГ-мод, горизонтальные и вертикальные оси которых направлены по осям цилиндрической линзы с соответствующим масштабированием координат [17 -21], а в сложной астигматической системе они могут перепутываться.
Ситуация значительно усложняется при использовании методов ЛБСБ-матриц для структурированных пучков, которые в настоящее время широко используются в различных областях науки и техники, требующих удобных математических подходов для проектирования оптических систем [22]. Прежде всего, это связано с тем, что в простейшем случае требуется представить ЛГ-моду в терминах ЭГ-мод в собственной системе координат цилиндрической линзы. Это выглядит так [17, 18]:
ЬОп±± ,(г) =
(_1)п 2п+1 (2)
- 1 ' -X (±20]Р'Т1_]'"_] \0)И02п+1_],] (г), ' '
2 n!
j=o
где Р(п+е_]'"_])(.) - многочлен Якоби, п и I - радиальное и азимутальное числа. Если оси ЭГ- мод повернуты на угол л / 4 относительно оси астигматизма линзы, то следует использовать разложение в виде
НГ (X1Z ) = nn+'( V2 ' V2 )
(_1)и 2n+t
(3)
22n+3H2 n!
S 2 JP}"+t_j'"_j л(0) НГ2п+1_ j, j (r),
j=o
а ЛБСБ-преобразование осуществляется для каждого направления х и у осей астигматического элемента.
В более сложном случае произвольной ориентации осей пучка относительно осей линз необходимо
использовать громоздкие базисные преобразования, полученные Алиевой и Бастиансом [23]. Разработка оптической системы с единичным астигматическим преобразованием структурированных пучков требует применения этих базисных преобразований к каждой моде пучка, и в то же время необходимо внимательно следить за амплитудами и фазами мод, возникающих в результате преобразований. Когда речь идет об оптической системе, содержащей последовательность астигматических элементов с различной ориентацией осей, для каждого элемента следует использовать базисные преобразования. В результате основной подход к астигматическим преобразованиям структурированных пучков был сосредоточен на анализе траекторий на 2Б-сфере с использованием унитарности астигматических преобразований [17, 24, 25]. Но здесь исследователь сталкивается с длинной последовательностью громоздких нестандартных вычислений, что не оптимально для инженерного моделирования. Таким образом, для инструментальной реализации устройств, преобразующих неустойчивый астигматический структурированный пучок в структурно устойчивый без потери его уникальных свойств, требуется разработка базисных методов ЛБСБ-матриц.
В данной статье мы теоретически и экспериментально обсудим ЛБСБ-преобразования в одной астигматической системе, которая наиболее важна для практической инженерии. Система содержит цилиндрическую линзу, формирующую всплески орбитального углового момента (ОУМ) и корректирующую сферическую линзу, позволяющую разделять астигматическую и вихревую компоненту ОУМ, а также преобразующую структурно неустойчивый астигматический структурированный пучок в устойчивый без потери его уникальных свойств.
Vw,(2.)
Рис. 1. (а) Эскиз астигматической системы с цилиндрической (СЬ) и сферической (БЬ) линзой; (б) взаимная ориентация поперечного сечения астигматического пучка в координатах цилиндрической линзы (х, у) и лабораторных координат (х', у')
1. Правило ЛВСВ для астигматических структурированных ЛГ-пучков (простой астигматизм)
Основной целью этого параграфа является поиск оптимальных условий для сохранения формы всплесков и быстрых осцилляций ОУМ, подчеркивающих вклад радиального числа п при распространении структурированного Лагерр-Гауссового (сЛГ-) пучка
через цилиндрическую линзу и разделении вихревого и астигматического ОУМ после корректирующей сферической линзы. Схема оптической системы показана на рис. 1а. Мы предполагаем, что структурированный sZГn>l-пучок [24] с радиальным числом n, азимутальным числом I и комплексным параметром q(z) = z - izo (где z0 = kwO, /2 - длина Рэлея) падает на цилиндрическую линзу с фокусным расстоянием fx, расположенную в плоскости z = 0, так что начальный
комплексный параметр равен qo = - izo и с радиусом перетяжки w0 Гауссова пучка. Сферическая линза выполняет преобразование Фурье и позволяет не только разделить вихревую и астигматическую составляющие ОУМ, но и преобразовать структурно неустойчивый пучок в структурно устойчивый без потери всплеска ОУМ при изменении параметров оптической системы. В целом, сложную оптическую систему удобно представлять в виде произведения матриц каждого оптического элемента [2]. Однако, поскольку преобразование структурированных сЛГ-пучков астигматическим элементом ранее не рассматривалось, мы сначала подробно рассмотрим преобразование сЛГ-луча на цилиндрической линзе, определив его комплексные параметры qx и q}, а затем используем их для анализа распространения пучка через остальные оптические элементы.
2. Структура сЛГ-пучка после единственной цилиндрической линзы
Подчеркнем, что цилиндрическая линза вводит различный масштаб по собственным координатам (х, у) (см. рис. 16). Предположим, что оси цилиндрической линзы и лабораторные координаты совпадают, что соответствует случаю так называемого простого астигматизма [17, 18, 19]. Поскольку ЭГ-пучки являются собственными модами астигматического элемента, представим структурированный сЛГ-пучок в терминах стандартных ЭГ-мод в выражении (2), но с различным масштабом по координатам (х, у) (см. также формулу (3.8) из монографии [19])
sLGnM (r, zi) = -
(-1)"
22n+3l/2 "I 2n+l
< X (±2i)j Qj (1 + eeij0 )lf+1- j,"- j) (0) >
(4)
j=0
V2x )H 42y ) г ik x2 y2
H+l- j (^—)Hj exp[- (■
)],
-l- j V J±J- j 4 XXL V
Wx (zi) Wx (Zi) 2 qx (Zi) qz (Zi)
где H„ - многочлен Эрмита, Qn - амплитуда
ЭГ-моды, которая в рамках ABCD-правила равна [2]
Qx, У = (
qo,x
Ax, y • qo, x, y + Bx
-)".
(5)
Очевидно, что в общем случае нам придется использовать две группы матриц ЛБСБ для направлений х и у, которые задают комплексные параметры
_ Axqo (Zi =0) + Bx _ Ayqo (Zi =0) + By
Cxqo (Zi =0) + Dx Cq (Zi =0) + Dy
. (6)
Однако цилиндрическая линза не изменяет масштаб в направлении у, то есть в области (0, zl), и мы можем написать
qy(Zi) = Zi -iZo.
(7)
В направлении х действует матрица цилиндрической линзы и смещение на длину z, которое записывается в виде
f
i
Л
--i
. f*
Откуда qx параметр получаем в виде [Zi( к2 +1)-к x ] - i
(8)
qx ( z ) = Zo
К 2 +1
где
ry Z1 Z0
Z1 = —, к x = —.
Zo fx
(9)
(10)
Полученные результаты позволяют найти радиусы пучков для х и у направлений следующим образом
ik
) = -.
1
2qx (Zi)' w2( Zi) w2( Zi) = Wo2[(1 - ZiK x )2 + Zi2],2 W?( Zi) = w2(1 + Zi2),2
в то время как фазы Гуи мод равны [2]
Г x (Zi) = arg
(11) (12)
= ^arg(1 -KxZ1+iZi), (13) A • qx + B 2
1
Г y (Zi) = -arg(1 + iZi).
(14)
Таким образом, после подстановки (9) - (14) в выражение (4) и последующих преобразований получаем комплексную амплитуду астигматического сЛГ-(асЛГ-) пучка после цилиндрической линзы в виде
asLGn l (x, y, Z1 | e, 0) = -
1
^Wx (Zi) Wy (Zi)
< exp{-i[(2n +1 +1/ 2)Г x (Z1) + Г y (Z1)]} >
(-1)n 2n+l
22n+3l/2 "|
X (2i)j Pjn+l - j-l> (0)(1 + eeij0);
(i5)
j=0
H
( V2x ) H t Лy
2"+l-j ( -Z—)H j ( -Z—) ' Wx (Zi) Wy (Zi)
cexp(i/r xy (Zi))exp(i[=
x2 y2 — + ]),
qx (Zi) qy (Zi)
где Гу (Zi) = rv (Zi) - Гх (Zi), wo = yj2Zo /k , x ^ x / Wo,
y ^ y / Wo,
qx (Zi) =
[Z1(k2 +1)-к x ] - i
к 2 +1 .
(16)
Заметим, что комплексная амплитуда (15) совпадает с амплитудой, полученной Бекшаевым и др. в работе [26] для простейшего случая с нулевым радиальным
числом п = 0. Компьютерное моделирование и экспериментальные результаты эволюции картин интенсивности (а - в) стандартной ЛГ-моды и (г - е) сЛГ-пучка вдоль направления распространения 2\ показаны на рис. 2. Как и ожидалось, ЛГ-пучок (рис. 2а) преобразуется в ЭГ-пучок на длине 21 = 2/ (21 = 1), а количество нулей интенсивности вдоль осей х и у позволяет определить топологический заряд ЛГ-пучка [27, 28]. Эволюция картины интенсивности структурированного пучка иллюстрирует по крайней мере один интересный эффект (рис. 2г). Вблизи 21 = 1 картина интенсивности асЛГ-пучка превращается в почти типичную картину
ЭГ-моды. Однако, как мы покажем ниже, ОУМ пучка испытывает резкий всплеск. Следует также отметить, что полученные результаты хорошо согласуются с картинами интенсивности, полученными методом интегральных преобразований в статье [24]. Однако метод интегральных преобразований позволяет получать достоверные результаты только в дальней зоне дифракции или в плоскости двойного фокуса цилиндрической линзы. В то же время представленные результаты, основанные на матричном ЛВСБ-подходе, значительно упрощают вычисления, приводящие к оптимальной форме комплексной амплитуды.
/,,-й 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
Рис. 2. Распределение интенсивности: (а,г) теория и (в,е) эксперимент, (б, д) фазы, эволюции ЛГ- (а-в) п = 2, I = 3 и астигматического сЛГ- (г-е) п = 10,1 = 1, е = 1, 0 = 0,99п пучка, с 2о = 1 м,/х = 0,5 м
3. Преобразования ОУМ сЛГ-пучка после цилиндрической линзы
Орбитальный угловой момент структурированного пучка в базисе ЭГ-мод удобно представить как
4 = * (г)(хд у - уд х )¥ (г)ёхёу. (17)
лк2
Очевидно, что основной вклад в ОУМ вносят амплитуды мод, которые мы находим из выражения (15)в виде
(,\)= (18)
= (2/) 1 Р(п+£-1 ■п -£ > (0)(1 + ее10) ехр(уГ ху (2)).
Чтобы найти удельный ОУМ, вычислим поток энергии в направлении распространения пучка по формуле
^ = ^(г) ¥*(г )йхйу, (19)
так что удельный ОУМ равен
£ 2 = 4 / ^. (20) Используя выражения (17 - 20), получаем
2п+1-1
¿2 = 22п+£ п X (1 + 1)!(2п + £ -1 )!3(с,+с), (21)
1=0
2п+£
= 22п+£-1 пХ 1!(2п + £ -1 )!сс. (22)
1=0
Мы обнаружили, что полученные выражения могут быть сведены к выражениям (8), полученным Котля-ром и др. в [27], несмотря на то, что в их расчетах масштабирование комплексной амплитуды одинаково по осям х и у. Полученные результаты для ОУМ представлены на рис. 3 в виде быстрых осцилляций ОУМ при изменении управляющего параметра 0 в различных поперечных сечениях 21 пучка. Мы видим возникновение и подавление быстрых осцилляций ОУМ по мере смещения 21 вдоль пучка, в то время как вблизи 0 = 0,98п зарождается и растет всплеск ОУМ. Второй всплеск ОУМ с противоположным знаком возникает вблизи 0 = 1,02п. Важно, что ОУМ при 21 = 1 хорошо согласуется с результатами нашей работы [24], использующей метод интегральных преобразований, но ЛВСБ матричный метод позволяет проследить эволюцию быстрых осцилляций и всплесков ОУМ по всей длине 21 пучка. Мы также видим, что всплески ОУМ
незначительно изменяются в дальнейшей зоне дифракции. Такое сохранение максимума всплеска ОУМ наглядно проиллюстрировано на рис. 4а. ОУМ достигает своего максимального значения при 71 = 1, не-
смотря на то, что картина интенсивности пучка значительно деформирована (см. рис. 46 - д), в то время как ОУМ стремится к половине радиального числа lz ^ п / 2 в дальней зоне дифракции.
4-
Г
-4-
-6 —
г, о,4 г, о,5 г, 1 г, 5
Рис. 3. Осцилляции ОУМ&(&) при вариации 0-^управляющего параметра асЛГ-пучка с п = 20,£ = 1, е пучка Ъ1. Выноски: картины интенсивности для первого всплеска ОУМ
-1 на различных длинах
г д
Рис. 4. (а) ОУМ£2(7.1) вдоль асЛГ-пучка с п = 20, £ = 1, 20 = 1 м, е = 1, 0 = 0,098к; (б)-(д) картины интенсивности в различных сечениях пучка
4. Цилиндрическая и корректирующая сферическая линза
Еще в начале 90-х годов Ананьев и Бекшаев предсказали в [29], что ОУМ сингулярного пучка имеет астигматическую и вихревую составляющие, которые можно разделеить с помощью обычной сферической линзы в плоскости, где радиусы Wx и Wy сфокусированного пучка вдоль направлений х и у становятся одинаковыми. Их анализ был основан на методе моментов интенсивности второго порядка. В дальнейшем этот подход был реализован для анализа как простых [21], так и структурированных сингулярных пучков [30]. В этом параграфе мы не будем углубляться в математику моментов интенсивности, а про-
сто проанализируем условия разделения вихревой и астигматической составляющих на основе ЛБСБ матричного подхода.
Стоит отметить, что согласование радиусов пучков wx = wy происходит и в плоскости двойного фокуса цилиндрической линзы, однако здесь не происходит разделения ОУМ на вихревую и астигматическую составляющие. Для этого необходимо также согласование Гх и Гу фаз. Эту операцию делает сферическая линза при условии, что она выполняет Фурье-преобразование. Это значит, что плоскость наблюдения преобразованного пучка после цилиндрической линзы должна располагаться в плоскости заднего фокуса сферической линзы, как показано на рис. 1. Для расчета мы используем матрицу (8), в которой выполняется замена / ^/н, 2 ^ 22. Тогда qx - комплексный параметр пучка определяется по рекуррентной формуле
qx 2
20 qx
-(20/ /л )qx +1
+ 20 72 20 qx2,
(23)
где
- 1 + П71 +кх (-1 + и1кх)] qx =-
qx2
1+к2
qx + 7 2[-к^ +1]
1 -к^х
и кн = 20 //и, 7.2 = 22 / 20. Теперь асЛГ-пучка следующим образом
(24)
находим х-радиус
1
Wx2
= -Ж-^ ) = -
2qx 2
1
Ж—).
qx 2
(25)
а х-фазу Гуи определяем как [2]
г Г^
Г х2 = arg 1=-.
Нх2
(26)
Из правила ABCD находим qv комплексный параметр
= (1 - Z2кsh )(Zi - i) + Z2 = -
qy 2 z0 . z0 qy 2,
1 + iKsh (iZi + 1)
радиус Wy и фазу Гуи
1 ik 1 i = -*(—-) = -— И(=-),
Vy 2
2qy 2
Г y 2 =arg(-^).
W02
qy2
(27)
(28) (29)
Теперь комплексная амплитуда пучка записывается в виде
asLG„t (х, y, z | e, 6) =
1
•¡Wx2(z)Wy2(z)
< exp{-i[(2n+1+1 / 2)Гх (z) / 2 + Г( z) / 2]}-
(-1)n
' 22n+3£/2 n! 2n+l
<X(2i) jPf+1-j ,n- j)(0)(1+eej6 )exp(i/Tk 2( z)) x
(30)
1=0
л/2х 4rr / 42y x2 y2
xH2n+t- j (=-)H (^-4exp(i[=-]),
wx2(z) wy2(z) qx2(z) qy 2 (z)
где Гк2 = Гх2 - Ty2.
Кривые на рис. 5 определяют условия, при которых происходит разделение вихревой и астигматической составляющих ОУМ, а также преобразование неустойчивого асЛГ-пучка в устойчивый пучок. Кривые (а, г) задают условия Wx(Z1, Z2) = Wy(Z1, Z2), при которых астигматическая составляющая ОУМ исчезает и только вихревая составляющая вносит основной вклад в ОУМ [21]. Можно показать, что эти условия выполняются при любом соотношении фокусных расстояний цилиндрических и сферических линз. Однако это не означает, что асЛГ-пучок становится структурно устойчивым после сферической линзы. Все же необходимо выполнить два дополнительных условия. Первое из них требует, чтобы разница между радиусами в направлениях х и y wx(Z1, Z2) - wy(Z1, Z2) = const оставалась неизменной, как показано на рис. 56, д. Это требование накладывает ограничение на неизменность астигматического компонента при прохождении через сферическую линзу (см. рис. 5е, е). Например, условия структурной устойчивости выполняются для всех пучков с параметрами, указанными на рис. 5г - е, однако структурно устойчивые пучки с параметрами, указанными на рис. 5е, е Z2 >> 1, не допускают супервсплесков ОУМ (см. также рис. 6 и рис. 7). Вторым дополнительным
условием является требование к положению сферической линзы
Z1=2 fx + fsh
(31)
что выполняется для кривых на рис. 5е. Теперь посмотрим, как преобразуются всплески ОУМ после корректирующей линзы.
где Рис. 5. Теоретические зависимости (а, г) ■«х(7.1,22) = му(21, 22); (б, д) Ам = мх- му; (в - е) мх(22) (темная кривая), м>у(22) (светлая кривая), для различных фокусных расстояний: (а - в) -/х = 0,5 и /¡и = 0,3; (г - е) -/х = 0,5 и /и = 0,3; (б, в) -21 = 0,3, (д, е) -21 = 0,5; (а - е) -20 = 1 М
I
2,521,510,5-
фф 7.,=0,85
f § Z, 1,5
фф /.=0.75
2 а
98765-4L 32-I-
%
т
Г
7., 0,85
о
2 б
Рис. 6. Зависимость ОУМ£¿(6) от управляющего параметра 6 асЛГ-пучка с (а) п = 4, I = 1, (б) п = 20, I = 1 на длине 21 = 2/х +/¡и м,/х = 0,5 м, 20 = 1 м
5. ОУМ сЛГ-пучка после корректирующей линзы
Для того чтобы рассчитать ОУМ, мы используем описанный выше подход. Для этого достаточно записать амплитуды ЭГ-мод из уравнения (30) в виде
cj (zb z2) = (2i) jPjn
'+£ - j,n- j )
(0)>
<(1 + eej6 )exp(jr k 2( zb z2)).
(32)
После подстановки уравнения (30) в уравнения (21) и (22) получаем ОУМ на фотон в виде многопараметрической функции I = €¿(21, 22,/х,/и, е, ), п, I). Мы сосредоточимся на оптической системе, в которой сферическая линза расположена в соответствии с условием (31). Рассмотрим особые точки в зависимости ОУМ 12(22) от смещения вдоль асЛГ-пучка в состоянии с малым п = 4 и большим п = 20 радиальными
числами и минимальным азимутальным числом I = 1, в то время как фазовый параметр 0 соответствует всплеску ОУМ в двойном фокусе Zl = 2fx асЛГ-пучка. Кривая ОУМ на рис. 6, окруженная картинами интенсивности, имеет резкий провал в точках Wx = Wy. (а) Ъ2 = 0,74, п = 4, I = 1 и (б) Z2 = 0,74, п = 20, I = 1, соответствующих условию коррекции астигматического пучка. Сферическая линза, выполняющая преобразование Фурье пучка со всплеском ОУМ превращает асЛГ-пучок в неастигматический пучок сЛГ-, соответствующий пучку на входе цилиндрической линзы. Как известно [24], максимум ОУМ в структурированном сЛГ-пучке не может превышать азимутальное число I = 1, в то время как максимальный ОУМ астигматического асЛГ-пучка превышает половину радиального числа I > п / 2. Ширина провала ОУМ зависит от фокусного расстояния fsh и быстро уменьшается с уменьшением фокусного расстояния. Затем
4
а
I
/
\А
[ У
у
ОУМ резко возрастает до своего первоначального значения, и его величина не меняется при смещении плоскости наблюдения вдоль пучка, его структура интенсивности также не меняется. Пучок становится структурно устойчивым вплоть до масштабирования и вращения. Структурная стабильность распространяется как на быстрые колебания, так и на всплески ОУМ. На рис. 7 показаны изменения формы осцилля-ций ОУМ вдоль пучка после сферической линзы. Если непосредственно за сферической линзой форма осцилляций ОУМ точно соответствует колебаниям в плоскости двойного фокуса цилиндрической линзы, то в плоскости согласования * и у радиусов пучков характер осцилляций резко меняется и принимает форму осцилляций ОУМ перед астигматическими преобразованиями [24]. Небольшое смещение от этой плоскости вдоль пучка возвращает форму осцилля-ций к исходной форме, содержащей всплески ОУМ.
\ Г У
ч /
]д [\1
1АА 1
V
6-
42 0-
Г
2
6- 4 "я " Зж Ш ' .7-
Л
г
Ъ=0
г2=о,б
г2=о,74
г2=о,88
Ъ=1
Рис. 7. Эволюция всплесков ОУМ асЛГ-пучка с (а) п = 4, I = 1, 0С = 0,98к и (б) п = 20, I = 1, 0С = 0,99к после корректирующей сферической линзы /¡и = 0,5 м, го = 1 м, 11 = 2/х +/¡и ). Выноски: картины интенсивности, соответствующие характерным
областям (указанных на рисунке) всплесков ОУМ
Выводы
Использование ЛБСБ матричного подхода, рассмотренного в статье, значительно упростило и расширило математическое описание астигматических преобразований структурированных пучков по сравнению с интегральными преобразованиями [24]. Мы продемонстрировали эффективность использования матриц ЛБСБ для преобразования ЭГ-мод высокого порядка в структурированных сЛГ-пучках. Было показано, что астигматическое преобразование приводит к возникновению всплесков ОУМ, амплитуда которых превышает половину радиального числа асЛГ-пучка с минимальным азимутальным числом. Более того, в отличие от метода интегральных преобразований, мы продемонстрировали эволюцию тонкой структуры интенсивности и ОУМ по всей длине пучка. Было показано, что ОУМ незначительно изменяется в дальней области дифракции цилиндрической линзы, несмотря на то, что сЛГ-пучок теряет свою структурную устойчивость.
Исследуя трансформацию асЛГ-пучка в системе цилиндрических и сферических линз, мы не только подтвердили предсказание Ананьева и Бекшаева о разделении вихревой и астигматической составляющих ОУМ после сферической линзы [21], но и показали, что сферическая линза способна превратить структурно неустойчивый асЛГ-пучок в структурно устойчивый в дальней области дифракции сферической линзы при условии, что сферическая линза выполняет преобразование Фурье асЛГ-пучка в двойном фокусе цилиндрической линзы. Мы обнаружили, что в асЛГ-пучке наблюдается резкий провал ОУМ, где радиусы пучков Wx(Z2) = Wy(Z2) совпадают, что соответствует подавлению астигматической составляющей ОУМ, так что основной вклад в ОУМ вносят оптические вихри пучка. Корректируемый пучок преобразуется в свою первоначальную форму сЛГ-пучка перед астигматическим преобразованием. Небольшое смещение плоскости наблюдения из этой области приводит к резкому увеличению ОУМ. Астигматический пучок становится структурно устойчивым
вплоть до масштабирования и вращения, сохраняя при этом форму осцилляций и всплесков ОУМ. Экспериментальные исследования хорошо согласуются с нашим компьютерным моделированием.
Благодарности
Авторы благодарят Е.Г. Абрамочкина за полезную дискуссию.
References
[1] Kogelnik H, Li T. Laser beams and resonators. Appl Opt 1966; 5(10): 1550-1567. DOI: I0.1364/A0.5.001550.
[2] Siegman AE. Lasers. Mill Valley, CA: University Science Books; 1986. ISBN: 0-935702-11-5.
[3] Soifer VA, Golub MA. Laser beam mode selection by computer-generated holograms. Boca Raton: CRC Press; 1994. ISBN: 0-8493-2476-9.
[4] Kotlyar VV, Kovalev AA, Porfirev AP. Vortex laser beams. Boca Raton, FL: Taylor & Francis Groop LLC; 2019. DOI: 10.1201/9781351009607.
[5] Forbes A, de Oliveira M, Dennis MR. Structured light. Nat Photonics 2021; 15: 253-262. DOI: 10.1038/s41566-021-00780-4.
[6] Shen Y, Yang X, Naidoo D, Fu X, Forbes A. Structured ray-wave vector vortex beams in multiple degrees of freedom from a laser. Optica 2020; 7(7): 820-831. DOI: 10.1364/OPTICA.382994.
[7] Fadeyeva T, Rubass A, Egorov Y, Swartzlander JrG. Quadrefringence of optical vortices in a uniaxial crystal. J Opt Soc Am A 2008; 25(7): 1634-1641. DOI: 10.1364/JOSAA.25.001634.
[8] Fadeyeva T, Alexeyev C, Rubass A, Volyar A. Vector erf-Gaussian beams: Fractional optical vortices and asymmetric TE and TM modes. Opt Lett 2012; 37(9): 1397-1399. DOI: 10.1364/OL.37.001397.
[9] Collins SA. Lens-system diffraction integral written in terms of matrix optics. J Opt Soc Am 1970; 60(9): 11681177. DOI: 10.1364/JOSA.60.001168.
[10] Bélanger PA. Beam propagation and the ABCD ray matrices. Opt Lett 1991; 16(4): 196-198. DOI: 10.1364/OL.16.000196.
[11] Taché JP. Derivation of ABCD law for Laguerre-Gaussian beams. Appl Opt 1987; 26(14): 2698-2700. DOI: 10.1364/AO.26.002698.
[12] Alieva T, Bastiaans MJ. Alternative representation of the linear canonical integral transform. Opt Lett 2005; 30(24): 3302-3304. DOI: 10.1364/OL.30.003302.
[13] Alieva T, Bastiaans MJ. Mode mapping in paraxial lossless optics. Opt Lett 2005; 30(12): 1461-1463. DOI: 10.1364/OL.30.001461.
[14] Beijersbergen MW, Allen L, van der Veen HELO, Wo-erdman JP. Astigmatic laser mode converters and transfer of orbital angular momentum. Opt Commun 1993; 96(1): 123-132. DOI: 10.1016/0030-4018(93)90535-D.
[15] Courtial J, Padgett MJ. Performance of a cylindrical lens mode converter for producing Laguerre-Gaussian laser modes. Opt Commun 1999; 159(1): 13-18. DOI: 10.1016/S0030-4018(98)00599-9.
[16] Khonina SN, Karpeev SV, Ustinov AV. Functional enhancement of mode astigmatic converters on the basis of application of diffractive optical elements [In Russian]. Izvestia of Samara Scientific Center of the Russian Academy of Sciences 2009; 11(5): 13-23.
[17] Abramochkin E, Volostnikov V. Beam transformations and nontransformed beams. Opt Commun 1991; 83(1): 123-135. DOI: 10.1016/0030-4018(91)90534-K.
[18] Abramochkin E, Razueva E, Volostnikov V. General astigmatic transform of Hermite-Laguerre-Gaussian beams. J Opt Soc Am A 2010; 27(11): 2506-2513. DOI: 10.1364/J0SAA.27.002506.
[19] Abramochkin E, Volostnikov V. Modern optics of Gaussian beams [In Russian]. Moscow: "Fizmatlit" Publisher; 2010. ISBN: 978-5-9221-1216-1.
[20] Wada A, Ohtani T, Miyamoto Y, Takeda M. Propagation analysis of the Laguerre-Gaussian beam with astigmatism. J Opt Soc Am A 2005; 22(12): 2746-2755. DOI: 10.1364/JOSAA.22.002746.
[21] Bekshaev AYa, Soskin MS, Vasnetsov MV. Optical vortex symmetry breakdown and decomposition of the orbital angular momentum of light beams. J Opt Soc Am A 2003; 20(8): 1635-1643. DOI: 10.1364/JOSAA.20.001635.
[22] Shen Y, Wang X, Xie Z, Min C, Fu X, Liu Q, Gong M, Yuan X. Optical vortices 30 years on: OAM manipulation from topological charge to multiple singularities. Light Sci Appl 2019; 8(90): 90. DOI: 10.1038/s41377-019-0194-2.
[23] Alieva T, Bastiaans MJ. Orthonormal mode sets for the two-dimensional fractional Fourier transformation. Opt Lett 2007; 32(10): 1226-1228. DOI: 10.1364/OL.32.001226.
[24] Volyar A, Abramochkin E, Akimova Y, Bretsko M. Astigmatic-invariant structured singular beams. Photonics 2022; 9(11): 842. DOI: 10.3390/photonics9110842.
[25] Volyar A, Abramochkin E, Bretsko M, Khalilov S, Akimova Y. General astigmatism of structured LG beams: Evolution and transformations of the OAM super-bursts. Photonics 2023; 10(7): 727. DOI: 10.3390/photonics10070727.
[26] Bekshaev AYa, Soskin MS, Vasnetsov MV. Transformation of higher-order optical vortices upon focusing by an astigmatic lens. Opt Commun 2004; 241(4): 237-247. DOI: 10.1016/j.optcom.2004.07.023.
[27] Kotlyar VV, Kovalev AA. Orbital angular momentum of paraxial propagation-invariant laser beams. J Opt Soc Am A 2022; 39(6): 1061-1065. DOI: 10.1364/JOSAA.457660.
[28] Volyar A, Abramochkin E, Egorov Yu, Bretsko M, Akimova Ya. Fine structure of perturbed Laguerre-Gaussian beams: Hermite-Gaussian mode spectra and topological charge. Appl Opt 2020; 59(25): 7680-7687. DOI: 10.1364/AO.396557.
[29] Anan'ev YuA, Bekshaev AYa. Theory of intensity moments for arbitrary light beams. Opt Spectrosc 1994; 76(4): 558-568.
[30] Fadeyeva TA, Rubass AF, Aleksandrov RV, Volyar AV. Does the optical angular momentum change smoothly in fractional-charged vortex beams? J Opt Soc Am B 2014; 31(4): 798-805. DOI: 10.1364/JOSAB.31.000798.
Сведения об авторах
Воляр Александр Владимирович - доктор физико-математических наук (1992). С 1993 г. - профессор, с 1999 г. - заведующий кафедрой общей физики. Сфера научных интересов: оптика, лазерная физика. E-mail: volyar.singular.optics@smail.com
Брецько Михаил Владимирович - кандидат физико-математических наук (2022), ассистент кафедры общей физики Физико-технического института Крымского федерального университета им. В.И. Вернадского. Сфера научных интересов: компьютерное моделирование, сингулярная оптика. E-mail: mihailbretcko4@smail.com
Халилов Сервер Искандарович - ассистент кафедры бизнес-информатики и математического моделирования Физико-технического института Крымского федерального университета им. В.И. Вернадского. Сфера научных интересов: сингулярная оптика, лазерная физика. E-mail: server.khalilov.94@mail.ru
Акимова Яна Евгеньевна - кандидат физико-математических наук (2022), ассистент кафедры общей физики Физико-технического института Крымского федерального университета им. В.И. Вернадского. Сфера научных интересов: сингулярная оптика, лазерная физика. E-mail: aklmova. yana@yandex.ru
ГРНТИ:29.31.15
Поступила в редакцию 16 августа 2023 г. Окончательный вариант - 13 сентября 2023 г.
Orbital angular momentum burst control in astigmatic structured beams
in ABCD-matrix transforms
A.V. Volyar1, M.V. Bretsko1, S.I. Khalilov1, Y.E. Akimova1 1 Physics and Technology Institute of V.I. Vernadsky Crimean Federal University, 295007, Simferopol, Republic of Crimea, Russia, Academician Vernadsky 4
Abstract
We developed and experimentally implemented an ABCD matrix formalism for describing structural transformations and orbital angular momentum of structured Laguerre-Gaussian beams in an astigmatic optical system containing cylindrical and corrective spherical lenses. It was shown that the matrix formalism not only agrees well with the method of integral astigmatic transformations, but also significantly expands the scope of its application. It was found that a corrective spherical lens is capable of not only separating the vortex and astigmatic components of the orbital angular momentum, but also transforming a structurally unstable astigmatic Laguerre-Gaussian beam after a cylindrical lens into a structurally stable one, provided that the form of fast oscillations and bursts of the orbital angular momentum is preserved.
Keywords: vortex beams; ABCD-matrix, structured light; orbital angular momentum.
Citation: Volyar AV, Bretsko MV, Khalilov SI, Akimova YE. Orbital angular momentum burst control in astigmatic structured beams in ABCD-matrix transforms. Computer Optics 2024; 48(2): 171-179. DOI: 10.18287/2412-6179-CO-1411.
Acknowledgements: the authors thank E. Abramochkin for the useful discussion.
Authors' information
Alexander Vladimirovich Volyar, (b. 1949) graduated from Simferopol State University (1972). After graduation he is its employee. Doctor of Physical and Mathematical Sciences (1992). Since 1993 he is a professor. Since 1999 he is a head of the General Physics department. His research spans various fields of optics and laser physics. E-mail: volyar.sinsular.optics@smail.com
Mikhail Vladimirovich Bretsko, (b. 1994) graduated from V.I. Vernadsky Crimean Federal University (2022), majoring in Physics and Astronomy. He received his Ph. D. in Physics in 2022. He is assistant of General Physics department in the Physical-Technical Institute of V.I. Vernadsky Crimean Federal University. Co-author of about 40 scientific papers. His research interests include computer modeling, singular optics. E-mail: mihailbretcko4@smail.com
Server IskandarovichKhalilov, (b. 1994) assistant of Business Informatics and Mathematical Modeling department of the Physical-Technical Institute of the Crimean Federal University named after V.I. Vernadsky. Research interests: singular optics, laser physics. E-mail: server.khalilov.94@mail.ru
Yana Evgenievna Akimova, (b. 1994) graduated from V.I. Vernadsky Crimean Federal University (2022), majoring in Physics and Astronomy. She received her Ph. D. in Physics in 2022. She is assistant of General Physics department in the Physical-Technical Institute of V.I. Vernadsky Crimean Federal University. Co-author of about 40 scientific papers. Her research interests include diffractive and singular optics, laser physics. E-mail: akimova. yana@yandex.ru
Received August 16, 2023. The final version - September 13, 2023.
