Формирование и анализ спектров оптических вихрей сингулярных пучков с аномалиями орбитального углового момента

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

Формирование и анализ спектров оптических вихрей сингулярных пучков с аномалиями орбитального углового момента

А.В. Воляр1, М.В. Брецько1, Я.Е. Акимова1, Ю.А. Егоров1 1 КФУ им. В.И. Вернадского, Физико-технический институт,

295007, Республика Крым, проспект Академика Вернадского, 4, Симферополь, Россия

Аннотация

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

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

Цитирование: Воляр, А.В. Формирование и анализ спектров оптических вихрей сингулярных пучков с аномалиями орбитального углового момента / А.В. Воляр, М.В. Брецько, Я.Е. Акимова, Ю.А. Егоров // Компьютерная оптика. - 2019. - Т. 43, № 4. - С. 517-527. -10.18287/2412-6179-2019-43-4-517-527.

Введение

Реальный сингулярный пучок, как правило, содержит множество оптических вихрей, даже если он сформирован специальным голографическим фильтром [1]. Вопрос только в том, как световой поток распределяется среди множества мод, переносящих единичные оптические вихри [2]. Существует множество способов измерения вихревого состава параксиальных пучков. Основное внимание уделяется измерению спектра орбитального углового момента (ОУМ), который непосредственно связан со спектром оптических вихрей. Действительно, если комплексную амплитуду Т(х, у, ¿) сложного пучка представить в виде суперпозиции комплексных амплитуд Тт (х, у, ¿) вихревых мод в виде

Т

ш=-Ы

где 2N - полное число мод в пучке, Ст - амплитуды

мод и X 1,-_ыСт = 1, то ОУМ на фотон 4 вдоль

направления распространения 2 задается простой формулой [3]

N 2

т=-N

В выражении (2) в неявной форме содержится условие численного равенства топологического заряда m оптического вихря в моде и ОУМ: tZm) = m .

Методы измерения спектра ОУМ можно условно разделить на три большие группы по способу влияния на структуру пучка в процессе эксперимента. К первой группе следует отнести методы голографиче-ской фильтрации вихревых мод, впервые детально рассмотренные в работах [4-6]. В основу метода положена голографическая фильтрация сложного пучка, представленного в базисе пучков Лагерра-Гаусса LGm таким образом, чтобы моды с орбитальным индексом (топологическим зарядом) m рассеивались по направлению азимутального угла фт. Радиальный индекс n отвечает за рассеяние вдоль полярного угла 0И. Второй способ предполагает использование специального рефракционного элемента и диафрагмы [7-9], которые осуществляют так называемое log-polar преобразование, при этом каждому значению топологического заряда m соответствует максимум интенсивности в горизонтальном (или вертикальном) направлении вдоль угла am. К третьей группе можно отнести метод, основанный на анализе преобразования Фурье с помощью сферической [10] или цилиндрической линзы [11-13]. В работах [11, 12] рассмотрен метод

измерения поля пучка с дробным ОУМ в фокальной плоскости цилиндрической линзы. Для измерения дробного ОУМ использовался подход, основанный на анализе перекрёстных элементов Jxy или Jyx матрицы интенсивности второго порядка [14]. Однако такой подход можно использовать, только если моды сложного пучка обладают осевой симметрией. Авторы работы [13] усовершенствовали метод за счет использования двух цилиндрических линз, оси которых повёрнуты на угол 90 ° относительно друг друга. Это позволило измерять дробный ОУМ пучков с любым типом симметрии.

Вместе с тем измерение дробного ОУМ не дает полной информации о состоянии сингулярного пучка. Как видно из выражения (2), одному и тому же значению ОУМ может соответствовать различный состав оптических вихрей с квадратами амплитуд Ст. Например, в слаботурбулентной среде [15] или оптических волокнах [16] имеет место интенсивный обмен энергией между оптическими вихрями сингулярного пучка и среды. При этом изменяется как величина так и спектр вихрей СЦ, (т). Для измерения спектра оптических вихрей был разработан специальный метод, основанный на анализе моментов интенсивности высших порядков [10] в фокальной плоскости сферической линзы. Характерной особенностью этого метода является возможность измерять как амплитуды, так и начальные фазы мод. Однако указанный подход имеет существенное ограничение. Метод позволяет измерять спектры вихрей только с одинаковыми знаками топологических зарядов. Дело в том, что моменты интенсивности в фокальной плоскости сферической линзы вырождены относительно знака топологического заряда. Поэтому авторы работы [10] отмечают, что метод применим только для невырожденных спектров сложных сингулярных пучков. Таким образом, целью данной статьи явилась демонстрация метода моментов интенсивности высших порядков для измерения вырожденных спектров вихрей и ОУМ сложных сингулярных пучков с действительными амплитудами.

1. Теоретическое обоснование

Представим комплексную амплитуду сложного сингулярного пучка в виде разложения по неполному базису мод Лагерра-Гаусса (ЬОот) с радиальным индексом п = 0 [10]

Т(Я, ф, z)= ]Г Ст1в0т)(Я, ф, z) =

т=-N

" СтЯте'т фО

где О = ехр (-Я2), Я = г/wо, г, ф - радиальные и азимутальные координаты соответственно, ст^) = 1 - iz/z0, z0 = к ч>0 /2, - радиус перетяжки каждой моды в

сечении z = 0, Мт = ^/2НтИ |т|!л - нормировочный

множитель, т - топологический заряд оптического вихря моды, 2Ы+1 - число мод в сингулярном пучке.

Основой наших исследований явилось измерение квадратов амплитуд мод Ст без разрушения структуры сингулярного пучка с помощью анализа момента интенсивности

ЛМрл

где 3(х, у) - распределение интенсивности поля (3) в плоскости z = 0, Мм (х, у) - функция момента интенсивности (см., например, [14] и ссылки на литературу), вид которой выбирается в зависимости от конкретных условий эксперимента (см. следующий раздел), р, q = 0, 1, 2, ... Как отмечалось в работе [10], моменты интенсивности Jp,q вырождены относительно знака топологического заряда оптических вихрей т в случае осевой симметрии мод в выражении (3). Чтобы снять вырождение, требуется изменить симметрию пучка без разрушения его сингулярной структуры таким образом, чтобы моды с противоположными знаками топологических зарядов имели отличительные признаки. Этому требованию отвечает астигма-тичное преобразование пучка Лагерра-Гаусса цилиндрической линзой, подробно рассмотренное в работе [17] (общую теорию параксиальных астигматичных преобразований можно найти в работе [18]).

Считаем, что цилиндрическая линза с фокусным расстоянием / расположена в плоскости z = 0. Рабочая ось линзы ориентирована вдоль оси у лабораторной системы координат. На вход линзы падает параксиальный монохроматический пучок с комплексной амплитудой (3). Плоскость перетяжки пучка совмещена с плоскостью линзы. Кроме того, ограничим наше исследование отсутствием начальных фаз модовых пучков, т.е. амплитуды Ст являются действительными величинами. Тогда, следуя работе [17], запишем комплексную амплитуду Т(х, у, z,/) сложного пучка в волновой зоне дифракции в виде

Т

т=-Ы

Здесь Нт(Ец) - полиномы Эрмита; ц = ± 1 - знак топологического заряда вихрей.

Введём новые нормированные координаты, оси которых ориентированы под углом л /4 относительно осей линзы: и = х / wo + у / wo, V = х / wo - у / wo. Будем рассматривать поле пучка в плоскости двойного фокуса г = 2/, а также потребуем, чтобы

*>/2/ = 1.

Тогда распределение интенсивности запишется в виде

1+

1 N 3= | С-2С-т 1 т,п=0 А N Нп и )Н- и 32

N + У СС т ,2=1 А N х ' т,2 Н2 ^т

N + У С С 2=0, А ^т,2 N 1 * т,2 Нп и )Нт [ 72

N

+Х с с-

п= 1, т=0

Мт

АтП-НА -41 Нп1 11

42) т ^л/2

х ехр

= 3 + з + з + з

л / 1\п + т 1 л

где Ат ,п = (-1)1 11 е 4

Ы2пт = Л( у) 2п+т+1 п!т!.

Как видно из соотношения (7), члены, соответствующие оптическим вихрям с положительными и отрицательными топологическими зарядами, только частично разделились. Последние два члена характеризуют перекрёстные амплитуды С-пСт и СпС-т. Кроме того, первые два члена также содержат перекрёстные амплитуды С-пС-т и СпСт(|п| Ф |т|). С другой стороны, наша задача - измерить только квадраты амплитуд модовых пучков С2т. Заметим, что при разности индексов |т - п| = 2 множитель Ат,2 = -Ап,т = / мнимый в (7) и, следовательно, перекрёстные члены с такой разностью индексов исчезают.

На первый взгляд, кажется, что для решения этой задачи целесообразно выбрать в качестве функции моментов интенсивности в выражении (4) многочлены Эрмита Мр (и) = Нр (и / >/2).

Далее можно воспользоваться условием ортогональности в (4) и записать систему линейных уравнений, число которых равно числу переменных, как это делалось в работе [10] для невырожденного массива вихрей. Однако, как показали оценочные расчёты, по мере увеличения индексов р и q, наряду с 2М неизвестными, для переменных С1т добавляется ещё 2Ы- 6 перекрёстных членов С±пС±т и С±2С+т. Систему уравнений невозможно замкнуть. Однако можно поступить иначе. Сначала сделаем два Фурье-преобразования распределения интенсивности 3(и, V) в виде

= | | 008I р(х,у)6хАу

I л/2)

да да / Л

^(л) = | {оо8У)АхАу.

Вид подынтегральных множителей в (8) показывает, что это преобразование уничтожает все перекрёстные члены в распределении интенсивности (7) с нечётной разностью индексов |т| - |п|. Значения интегралов в (8) найдем в работе [19]:

{-=2) нп ^ г п+2 „ ^

= (-1)" 22+1/2л/Л ю2" ехр [- — | ¿2

е - (и2/2)аи =

Г 008 I ю IН1I |е 2 Аи = -1 1 >/2)

= 22+1/2л/Л ехр| -^1 ¿2 2

Г 008 [«-^1Н 22 |-7=| е~ 2Аи =

л/2) 2 |л/2"

= (-1)2 22+1/2 ю22 ехр (-

ю 2

Г Н„

-да -да

В частности, члены в распределении интенсивности имеют вид

(£) = Г Г 008у)Ахёу =

N 22+1/2 (е2 I

= у С-2. -М- (V) е т +

2 =0 N

X С-2С

(-1)2+ 2т 22+1/2

-(2+2т) "

т,п =0, |т - п| ф0,2

х^2 тЩт (у ) е "4,

N 1)2 + 2т 22+1/2

2 = 0, т=1

N2.

Юе% }Н„ е ^ = 0.

I

М)

Таким образом, распределение интенсивности л) после Фурье-преобразования содержит все квадраты амплитуд, и перекрёстные члены с разностью индексов |т| - |2| = 4. Перекрёстные члены 3-т,т и 3т, -„, ответственные за вихри с положительными и отрицательными зарядами, исчезают. Причем члены 3_т,чи+2т)© содержат амплитуды С-п с одинаковыми коэффициентами, а члены 3т,(и+2т)© содержат амплитуды С,2 с различными коэффициентами. Аналогичные выражения получаем для членов 3-т-(и+2ш)(л) и 3

т,(и+2т)(л), если поменять знаки индексов.

т=1

Чтобы найти уравнения для квадратов амплитуд, воспользуемся разложением моментов интенсивности по полному базису многочленов Лагерра:

I

(+) _

р,2р

ОО^ /1о<

Ц / 2 \

I= |п Т Ъ [^]з(ц)1ц /I

где 100 - полная интенсивность пучка. Дополним эти уравнения интегралами [19]

| х2т+112„т (с х2 ))т (с х2) е—^ йх=

—да

(2т + п )! ^

п!с2т+' п,к'

да

| х2т+хЬ2пт (с х2) (сх2 ) е—сх'йх =

—да

—1) [ 2т + п + Х

| х2д+1е—ьр (сх)с1х =

(2т + Х)!.

2т

(2р + р )!(1/2 )р

р !4—

С учётом интегралов (16)-(18) уравнения для амплитуд мод (15) Сп запишутся в виде

X С—2

(2р + р )!(1/2)Р

р!4

-

1 р,р =

X (2т + п)()п

^ п!4—(2т+П+1) N

т,п =1 11 • ~ пт

X С

п=0

N

I + =

100

( + р )!(1/2)р

р !4—

(2т + п )!(1/2 )п

п!4—(2т+п+1)пт

10(

Система уравнений (19) и (20) позволяет найти значения квадратов амплитуд С1„ по измеряемым значениям моментов интенсивностей . Каждое уравнение содержит N членов для амплитуд С1п и N - 3 перекрёстных членов. Изменение индексов р и р не увеличивает число переменных. Полное число уравнений равно 4N - 6. Систему уравнений легко решить и найти амплитуды модовых пучков в реальном масштабе времени с помощью современного компьютерного обеспечения.

2. Эксперимент и обсуждение результатов

2.1. Экспериментальная установка и калибровочные измерения

Для решения системы линейных уравнений (19) и (20) относительно амплитуд С2т требуется измерить моменты интенсивности 1рр высших порядков [10, 14].

Для этого использовалась экспериментальная установка, схема которой приведена на рис. 1.

¿л

О £ , 5ьр о

яьм

ССВ1

Рис. 1. Схема экспериментальной установки. Ьз - Ие^е лазер, Р - поляризатор, БЬЕ - пространственный фильтр пучка, БЬМ - пространственный модулятор пучка, Ь1 - сферическая линза, ВБ - делительный кубик, СЬ - цилиндрическая линза, БТ - диоптрийная трубка, ССБ1, ССБ2 - ССБ-камеры

Свет от Не-№ лазера (X = 0,63 мкм) на основной продольной моде ТЕМ00 проходит поляризатор Р и пространственный фильтр пучка SLF. Затем свет направлялся на пространственный модулятор пучка SLM (ТЪог1^ ЕХиЬШ-4К1), на ПЗС-матрице которого формировалась цифровая голографическая решетка, сформированная компьютером в соответствии с комплексной амплитудой (3). Кроме того, SLF система совместно со сферической линзой L1 позволяла получить требуемый радиус перетяжки в плоскости цилиндрической линзы СЬ. Восстановленный пучок разделялся на два плеча. В первом плече линза L1 преобразовывала пучок таким образом, чтобы его плоскость перетяжки совпадала с плоскостью цилиндрической линзы СЬ. Фокусное расстояние цилиндрической линзы /су1 и длина Релея пучка г0 = к^ /2 соответствовали условию (6) 20/2/= 1. Измерение моментов интенсивности 1рр осуществлялось ССD1 камерой (Thorlabs DCC1645), расположенной в плоскости г = 2fyl. Контроль формы перетяжки пучка осуществлялся во втором плече, оптическая схема которого построена по принципу микроскопа бокового вида. Радиус перетяжки пучка в плоскости цилиндрической линзы вычислялся посредством измерения моментов интенсивности второго порядка 102 = ^у; и 120 = и ^02 = ^ + ^2)/2 [10, 20] на ССD2.

Измерение квадратов амплитуд Ст проводилось следующим образом. Сначала на CCD1 камере регистрировалось распределение интенсивности 3(х, у), а массив данных заносился в память компьютера. Затем компьютер производил преобразования (8), и записывался массив {3±т,±п(|), 3±т,±п(ц)}. Далее над этим массивом производились преобразования (15). Результаты компьютерной обработки информации в виде моментов интенсивности 1т,р подставлялись в

X!

левую часть уравнений (19) и (20). Решение уравнений выводились на монитор в виде спектра вихрей С2 (т) и спектра ОУМ 4.

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

Калибровочные измерения 1. Сначала отметим, что в общем случае аргумент многочлена Эрмита в (7) является комплексным за исключением плоскости г = 2 /су1 и выполнения условия (6). Это значит, что такой пучок переносит оптический вихрь и обладает ОУМ. Аналогичный вопрос о суперпозиции пучков Эрмита-Гаусса уже рассматривался в работе [23]. Погрешность измерений спектра вихря при небольшом отклонении от плоскости измерений г = 2/су1 и при нарушении условия г0 = 2 /су1 можно оценить, измеряя отклонение ОУМ единичных модовых пучков Лагерра-Гаусса от целочисленных значений. Для этого предположим, что величина отклонения Дг очень мала, т.е. 2 = 10/2/су1 = 1+Дг, Дг << 1. Тогда, используя соотношение

Hm

I tx + y I

л/1

+12

ml

+ t2 j=0 V J

I m I

XI . I tkHk

и уравнение (5), ограничиваясь только членами первого порядка малости Д г, получаем для ОАМ простое соотношение 4 ~ Дхт. Следовательно, можно ожидать возрастание погрешности измерений спектра вихрей по мере увеличения топологического заряда вихря т в модовом пучке (3). Таким образом, на SLM модуляторе формировалась дифракционная решётка мод Ла-герра-Гаусса высших порядков, входящих в сумму в исходном выражении (3), с распределением интенсивности, представленным на рис. 2. Поскольку отклонение Д г вызывает появление т мод для каждого ЬО0т пучка, измерялся весь спектр вихрей для используемых в эксперименте пучков.

Измерительный процесс начинался с формирования Гауссова пучка в плоскости цилиндрической линзы СЬ с заданным радиусом перетяжки w0, который удовлетворяет условию (6): w0 = 2^[/Тк . В нашем

случае фокусное расстояние цилиндрической линзы есть /су1 = 3,5 см, что соответствует радиусу перетяжки w0 = 37,4 мкм. На эксперименте радиус перетяжки регулировался с помощью пространственного фильтра SLF и микроскопа бокового вида на рис. 1 и измерялся в соответствии с методом, описанным в работе [10]. Нам удалось установить диаметр перетяжки равным 2w0 « 78 ± 1,5 мкм с оптимальной точностью для данного эксперимента. Следующим шагом явилось измерение спектра вихрей и ОУМ для каждого модового пучка последовательности на рис. 2. Ре-

зультаты измерений представлены на рис. 3. Хотя при измерениях вклад сателлитов в спектрах составляет доли процента, погрешность измерений ОУМ вихрей высших порядков может быстро возрастать с увеличением индекса т. В нашем эксперименте максимальная ошибка не превышала Д4/4 ~ 2 %.

т=±5 т=±10 т=±15 т=±20 т=±25 т=±30

LGom О о о о о с

HGom "t ч *ч \

HGom i* ✓ и / /

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

С2

^т 0,8 0,6 0,4 0,2

0

-30 -20 -10 0 10 20 т

Рис. 3. Спектры вихрей \Cm(m)\2 единичныхLGo,m модовых пучков, измеренных в ходе калибровочного эксперимента

Калибровочные измерения 2. Вторая серия контрольных экспериментов предполагала измерение сложных спектров вихрей со случайным выбором амплитуд мод в соответствии с выражением

I cos2 p (am + c)m :Л

С2т (m) = -

exp

X cos2p (an + c)"

exp

где О, а, Ь, с - параметры пучка. Небольшие изменения индекса р = 0, 1, 2, 3... и управляющих параметров {а, Ь, с} вызывают резкие изменения вида спектра вихрей (22). Параметры спектра выбирались так, чтобы охватить область как положительных, так и отрицательных топологических зарядов. Для определения квадратов амплитуд в каждом спектре потребовалось провести 60 измерений, составить и решить 60 линейных уравнений (19) и (20). Оптимальные результаты измерений представлены спектром С|т (т) на рис. 4.

Заметим, что для выбора оптимальных результатов потребовалось проводить контроль степени корреля-

1

ции исходного распределения интенсивности 3(г, ф) и интенсивности Зс (г, ф), восстановленной в соответствии с измеренным спектром С1т (т) на рис. 4:

K = jj 3„ (г, ф) Зс, (г, ф) dS / J* JC

где J0"o и J00 - полные интенсивности исходного и восстановленного пучков. Если степень корреляции была ниже К = 0,9, то узлы экспериментальной установки подвергались дополнительной юстировке. В данном случае степень корреляции составила К = 0,93, а погрешность измерения квадратов амплитуд Ст (т) была не хуже 4,5 %.

Ст •

0,2

0,1

- — - 1 1 1 > |

V*

y я 6 [

1

Л JL -1. JL _ Jj L _i

0

-30 -20 -10 0 10 20 т

Рис. 4. Рассчитанные (серый) и измеренные (чёрный) спектры вихрей \Ст\2для выражения (23): а = 14, Ь = 0,004, с = -5, р = 20, О = 1,05, N = 30. Выноска: распределения интенсивности Зс(г,ф):

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

2.2. Срыв лавины вихрей при слабом возмущении голографической решётки

В этом подпараграфе мы применим метод измерения спектра вихрей и ОУМ в новом типе комбинированных пучков с выраженными провалами в спектре ОУМ. Свойства такого пучка невозможно было исследовать методом, предложенным в работе [10]. Действительно, в работах [10, 21, 22] было показано, что при соответствующем выборе амплитуд Ст вихревых мод ЬО0т в спектре ОУМ возникают глубокие резонансные всплески и провалы при условии совпадения целочисленного значения параметра возмущения р с топологическим зарядом т одной из вихревых мод: р = т. Чтобы понять физический механизм возникновения резонансов, в работе [21] было показано, что при малых локальных возмущениях 5р рельефа голографической решётки, настроенной на формирование т порядка вихря, возникает лавина оптических вихрей с широким спектром как положительных, так и отрицательных топологических зарядов. В области провалов ОУМ вклад вихрей с т < 0 и т > 0 при р > 0 был сравним между ними, в то время как в области всплесков ОУМ вклад вихрей с т > 0 был существенно большим.

Поэтому можно было использовать метод измерения в работе [10] только для исследования всплесков ОУМ. В этом подпараграфе мы сфокусируем внимание на теоретическом и экспериментальном исследовании физического механизма провалов в спектре ОУМ. Для анализа структуры возмущенного пучка удобно ввести параметр слабого возмущения 5р рель-

ефа голографической решетки как отклонение целочисленного индекса числа ветвей вилочной дислокации М: р = М + 5р, где р - параметр возмущения голо-графической решетки. Тогда комплексная амплитуда пучка, восстановленного голографической решеткой, запишется в виде

Y (г,ф,z) = sin(p%)e'p% X

cos

/

d ml

Rme

|m|!(( - m + 5p) CT(z)H+1

exp

где d - масштабный параметр пучка; множители cos(m% /2) или sin (m% /2) прореживают спектр вихрей, оставляя моды либо с чётными, либо с нечётными индексами m.

При этом провалы в спектре ОУМ становятся более отчётливыми. Комплексной амплитуде (24) соответствует ОУМ в виде

cos2 (m % /2)) d2lh

sin2 (m%/2) |m|!2 (( -m +5p)2

£ z =-

N

X m

h=-N

X

2

cos2

Л

d2l h

sin2 (m % /2)) |m| !2 (( - m + 5p)2

На рис. 5 представлены спектральные кривые 4(т) для случая чётных топологических зарядов. Спектр ОУМ на рис. 5а описывает общий вид областей резонансных всплесков и провалов, соответствующих параметру М. Масштабный параметр О позволяет сдвигать резонансные области вдоль оси р = М + 5р. Их характерный вид иллюстрируют кривые 1 и 2 на рис. 5б. Все резонансные провалы локализуются вблизи чётных топологических зарядов вихревых мод. Замена косинуса на синус в выражении (25) сдвигает резонансы на единицу вдоль ось р. Уменьшение индекса М сопровождается резким сжатием резонансов, при этом их глубина составляет примерно 2М при О = 14 и р < 6. При больших значениях параметров О >> 20 наблюдается удвоение ОУМ 4 = 2М в широком интервале исходных значений топологического заряда М.

Как было показано в работе [21], возмущение 5р вносит локальное искажение структуры вилочного дефекта при больших значениях параметра р. В рассматриваемом случае, как показывает форма голо-графических решёток на рис. 6, состояния с малыми значениями М подвержены существенным искажениям всего рельефа голограммы даже при возмущении 5р = -0,0001 (см. рис. 6б). Такое различие в рельефе решёток неизбежно сказывается на форме спектра вихрей С1т (т).

Действительно, спектр вихрей Ст (т) комбинированных пучков, рассмотренных в работе [21], содержит сателлиты в области как положительных, так и отрицательных топологических зарядов. Однако в данном случае перекачка энергии из основного со-

S

стояния осуществляется в сателлиты с меньшими топологическими зарядами, что и вызывает провалы в спектре ОУМ. Перестройка спектра С2т (т) показана на рис. 7.

20

15

10

5

0

а)

к. 8765 4 3 2

1 О

б)

Рис. 5. Теоретические спектральные кривые ОУМ ^

Рис. 6. Возмущение вилочного дефекта голографической

решётки сМ = 4 и П= 14: §р = -0,0001 (а), §р = -0,5 (б), П= 14. Выноска: распределения интенсивности в пучках, восстановленных соответствующими голограммами

Уже малые возмущения голографической решётки (см. рис. 7а) вызывает срыв лавины вихрей как с положительными, так и отрицательными зарядами таким образом, что перекачка энергии осуществляется в сателлиты с более высокими топологическими зарядами, чем заряд исходного пучка. По мере роста возмущения 8р вклад сателлитов в ОУМ справа от исходной моды (рис. 7б) становится значительно больше, чем от сателлитов с отрицательными топологическими зарядами. Такая перекачка энергии формирует глубокие провалы в спектре ОУМ I г (¿), показанные на рис. 8.

0,7 0,6 0,5 0,4 0,3 0,2 0,1 О

а) -20 -10 0 10 т

0,20

0,15

0,10

0,05

О

б) -20 -10 0 10 т

Рис. 7. Экспериментальные спектры вихрей С2т возмущённого пучка с топологическим зарядом М = 4: (а) 8р = - 0,03, (б) 8р = - 0,5, П= 14;

исходное состояние

4 7

б

5 4 3 2

1 О

Рис. 8. Спектр ОУМ ¿г

В области топологического заряда моды т = 2 глубина провала происходит от 4=4 до 4=2, а в области т = 4 провал составляет 4 = 4(6,8) - 4(4) « - 2,8. Дальнейшее возрастание параметра р сопровождается сглаживанием провалов 0, а затем провалы сменяются резкими всплесками ОУМ.

Заключение

Мы предложили и применили новый метод измерения спектров оптических вихрей и ОУМ комбини-

4 (п= 14 у

V /

Ж 1(,

А

г* /

/ Т —\

£ У г

V

г

ь У

у р

0 5 10 15 20

а=14

£

/

"К

П=0,4

О 1 2 3 4 5 б р

в

8рЬ су1 -

- - - -- -

_____ -1-1- 1____ _____ 1-1-1 -1- -—

'с* -

гт _

$

Ш Ч- Ш

% *

V % _

зрк су1 -

.1 СГ 1

1.

рованных сингулярных пучков, содержащих вихри как с положительными, так и отрицательными топологическими зарядами. Метод основан на измерении моментов интенсивности высших порядков комбинированного пучка в плоскости двойного фокуса цилиндрической линзы при условии равенства длины Релея исходного Гауссова пучка в плоскости цилиндрической линзы удвоенному фокусному расстоянию этой линзы. На основании измерений составлялось 4N - 6 линейных уравнений (где 2N - число вихревых мод в пучке) и строились спектры квадратов амплитуд мод и ОУМ. Предварительные калибровочные измерения показали, что точность измерений не превышает 4,5 %, а степень корреляции исходного и восстановленного пучков не ниже К = 0,93.

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

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

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 19-2901233; Грант №ВГ24/2018 «Крымский федеральный университет имени В.И. Вернадского».

Литература

1. Soifer, V.A. Laser beam mode selection by computergenerated holograms / V.A. Soifer, M.A. Golub. - Boca Raton: CRC Press, 1994. - 215 p. - ISBN: 978-0-8493-2476-5.

2. Kotlyar, V.V. Vortex laser beams / V.V. Kotlyar, A.A. Kovalev, A.P. Porfirev. - Boca Raton, London, New York: Taylor & Francis Group, 2018. - 404 p. - ISBN: 9781-138-54211-2.

3. Berry, M.V. Paraxial beams of spinning light. Proceedings of SPIE. - 1998. - Vol. 3487. - P. 6-11. - DOI: 10.1117/12.317704.

4. Khonina, S.N. Measuring the light field orbital angular momentum using DOE / S.N. Khonina, V.V. Kotlyar,

V.A. Soifer, P. Paakkonen, J. Turunen // Optical Memory and Neural Networks (Information Optics). - 2001. -Vol. 10, Issue 4. - P. 241-255.

5. Khonina, S.N. Generation and selection of laser beams represented by a superposition of two angular harmonics / S.N. Khonina, V.V. Kotlyar, V.A. Soifer, K. Jefimovs, J. Turunen // Journal of Modern Optics. - 2004. - Vol. 51, Issue 5. - P. 761-773. - DOI: 10.1080/09500340408235551.

6. Khonina, S.N. Optical vortices in a fiber: mode division multiplexing and multimode selfimaging / S.N. Khonina, N.L. Kazanskiy, V.A. Soifer. - In: Recent progress in optical fiber research / ed. by M. Yasin, S.W. Harun, H. Arof. -Chap. 15. - Rijeka, Croatia: InTech Open, 2012. - P. 327352. - DOI: 10.5772/28067.

7. Gbur, G.J. Singular optics / G.J. Gbur. - New York: CRC Press, 2016. - 545 p. - ISBN: 978-1-4665-8077-0.

8. Lavery, M.P.J. Measurement of the light orbital angular momentum spectrum using an optical geometric transformation / M.P.J. Lavery, G.C.G. Berkhout, J. Courtial, M.J. Padgett // Journal of Optics. - 2011. - Vol. 13, Issue 6. - 064006.

9. Berkhout, G.C.G. Efficient sorting of orbital angular momentum states of light / G.C.G. Berkhout, M.P.J. Lavery, J. Courtial, M.W. Beijersbergen, M.J. Padgett // Physical Review Letters. - 2010. - Vol. 105. - 153601.

10. Volyar, A.V. Measurement of the vortex spectrum in a vortex-beam array without cuts and gluing of the wavefront / A. Volyar, M. Bretsko, Ya. Akimova, Yu. Egorov // Optics Letters. - 2018. - Vol. 43, Issue 22. - P. 5635-5638. - DOI: 10.1364/OL.43.005635.

11. Alperin, S.N. Quantitative measurement of the orbital angular momentum of light with a single, stationary lens / S.N. Alperin, R.D. Niederiter, J.T. Gopinath, K.E. Siements // Optics Letters. - 2016. - Vol. 41, Issue 21. - P. 50195022. - DOI: 10.1364/OL.41.005019.

12. Alperin, S.N. Angular momentum of topologically structured darkness / S.N. Alperin, K.E. Siements // Physical Review Letters. - 2017. - Vol. 119. - 203902.

13. Котляр, В.В. Методы определения орбитального углового момента лазерного пучка / В.В. Котляр, А.А. Ковалёв, А.П. Порфирьев // Компьютерная оптика. - 2019. -Т. 43, № 1. - С. 42-53. - DOI: 10.18287/2412-6179-201943-1-42-53.

14. Flusser, J. Moments and moment invariants in pattern recognition / J. Flusser, S. Tomás, B. Zitová. - New York: John Wiley & Sons, 2009.

15. Fu, S. Influences of atmospheric turbulence effects on the orbital angular momentum spectra of vortex beams / S. Fu, Ch. Gao // Photonics Research. - 2016. - Vol. 4, Issue 5. -P. B1-B4.

16. Kaiser, T. Complete modal decomposition for optical fibers using CGH-based correlation filters / T. Kaiser, D. Flamm, S. Schröter, M. Duparré // Opt Express. - 2009. - Vol. 17. -9347-9356.

17. Kotlyar, V.V. Astigmatic transforms of an optical vortex for measurement of its topological charge / V.V. Kotlyar, A.A. Kovalev, A.P. Porfirev // Applied Optics. - 2017. - Vol. 56, Issue 14. - P. 4095-4110. - DOI: 10.1364/AO.56.004095.

18. Abramochkin, E. General astigmatic transform of Her-mite-Laguerre-Gaussian beams / E. Abramochkin, E. Razueva, V. Volostnikov // Journal of the Optical Society of America A. - 2010. - Vol. 27, Issue 11. - P. 2506-2513.

19. Prudnikov, A.P. Integrals and series, elementary functions / A.P. Prudnikov, Y.A. Brychkov, O.I. Marichev. - New York: Gordon and Breach, 1986. - 708 p.

20. Воляр, А.В. По ту сторону интенсивности или моменты интенсивности и измерение спектра оптических вихрей слож-

ных пучков / А.В. Воляр, М.В. Брецько, Я.Е. Акимова, Ю.А. Егоров // Компьютерная оптика. - 2018. - Т. 42, № 5. -С. 736-743. - Ш1: 10.18287/2412-6179-2017-42-5-736-743.

21. Воляр, А.В. Лавинная неустойчивость орбитального углового момента оптических вихрей высших порядков / А.В. Воляр, М.В. Брецько, Я.Е. Акимова, Ю.А. Егоров // Компьютерная оптика. - 2019. - Т. 43, № 1. - С. 14-24. - DOI: 10.18287/2412-6179-2019-43-1-14-24.

22. Volyar, A.V. Super pulses of orbital angular momentum in fractional-order spiroid vortex beams / A.V. Volyar, Yu.A. Egorov // Optics Letters. - 2018. - Vol. 43, Issue 1. -P. 74-77. - DOI: 10.1364/OL.43.000074.

23. Котляр, В.В. Лазерные пучки Эрмита-Гаусса с орбитальным угловым моментом / В.В. Котляр, А.А. Ковалёв, А.П. Порфирьев // Компьютерная оптика. - 2014. - Т. 38, № 4. - С. 651-656.

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

Воляр Александр Владимирович - 1949 года рождения. Выпускник Симферопольского государственного университета (1972), а после окончания - его сотрудник. Доктор физико-математических наук (1992). С 1993 г. - профессор, с 1999 г. - заведующий кафедрой общей физики. Сфера научных интересов: оптика, лазерная физика. E-mail: volyar.sinsular.optics@smail.com .

Брецько Михаил Владимирович - 1994 года рождения. С 2012 года студент Крымского федерального университета им. В.И. Вернадского. На данный момент аспирант первого курса кафедры общей физики физико-технического института КФУ им. В.И. Вернадского. Сфера научных интересов: сингулярная оптика, обработка графических изображений, программирование. E-mail: mikhailbretcko4@smail.com .

Акимова Яна Евгеньевна - 1994 года рождения. С 2012 года студент Крымского федерального университета им. В.И. Вернадского. На данный момент аспирант первого курса кафедры общей физики, физико-технического института КФУ им. В.И. Вернадского. Сфера научных интересов: сингулярная оптика, дифракционная оптика. E-mail: ak1mova.yana@yandex.ru .

Егоров Юрий Александрович - 1975 года рождения. Выпускник Таврического национального университета (в настоящее время КФУ). Получил степень кандидата физико-математических наук в области оптики и лазерной физики в Таврическом национальном университете в 2005 году. Начиная с 2006 года доцент кафедры общей физики. Сфера научных интересов: сингулярная оптика, дифракционная оптика, оптика анизотропных сред. E-mail: yuriyegorov@cfuv.ru .

ГРНТИ: 29.31.15

Поступила в редакцию 28 марта 2019 г. Окончательный вариант - 11 июня 2019 г.

Shaping and processing the vortex spectra of singular beams with anomalous orbital angular momentum

A.V. Volyar1, M.V. Bretsko1, Ya.E. Akimova1, Yu.A. Egorov1 1 Physics and Technology Institute of V.I. Vernadsky Crimean Federal University, Academician Vernadsky 4, 295007, Republic of Crimea, Simferopol, Russia

Abstract

The article examines physical mechanisms responsible for shaping the vortex avalanche induced by a weak perturbation of the holographic lattice of a combined vortex beam. For this, we have developed a new technique for measuring the degenerate spectra of optical vortices and orbital angular momentum of combined singular beams. The technique is based on measuring the intensity moments of higher orders of a beam containing vortices with both positive and negative topological charges. The appropriate choice of the mode amplitudes in the combined beam enables us to form orbital angular momentum anomalous spectral regions in the form of resonance dips and bursts. Since the intensity moments of a vortex mode with positive and negative topological charges are the same (the moments are degenerate) for an axially symmetric beam, the measurements are carried out in the plane of the double focus of a cylindrical lens. The calibration measurements show that the experimental error is not higher than 4.5 %. We also reveal that the dips and bursts in the orbital angular momentum spectrum are caused by the vortex avalanche induced by weak perturbations of the holographic grating relief responsible for the beam shaping. The appearance of the orbital angular momentum dips or bursts is controlled by the relation between the energy fluxes in the vortex avalanche with positive or negative topological charges.

Keywords: diffractive optics, image processing, optical vortices, orbital angular momentum, moments of intensity.

Citation: Volyar AV, Bretsko MV, Akimova YaE, Egorov YuA. Shaping and processing the vortex spectra of singular beams with anomalous orbital angular momentum. Computer Optics 2019; 43(4): 517-527. DOI: 10.18287/2412-6179-2019-43-4-517-527.

Acknowledsements: The work was funded by the Russian Foundation for Basic Research under RFBR grant No. 19-29-01233 and under grant Br24/2018, "V.I. Vernadsky Crimean Federal University Development Program in 2015 - 2024".

References

[1] Soifer VA, Golub MA. Laser beam mode selection by computer-generated holograms. Boca Raton: CRC Press; 1994. ISBN: 978-0-8493-2476-5.

[2] Kotlyar VV, Kovalev AA, Porfirev AP. Vortex laser beams. Boca Raton, London, New York: Taylor & Francis Groop LLC; 2018. ISBN: 978-1-138-54211-2.

[3] Berry MV. Paraxial beams of spinning light. Proc SPIE 1998; 3487: 6-11. DOI: 10.1117/12.317704.

[4] Khonina SN, Kotlyar VV, Soifer VA, Paakkonen P, Turunen J. Measuring the light field orbital angular momentum using DOE. Optical Memory and Neural Networks 2001; 10(4): 241-255.

[5] Khonina SN, Kotlyar VV, Soifer VA, Jefimovs K, Turunen J. Generation and selection of laser beams represented by a superposition of two angular harmonics. J Mod Opt 2004; 51(5): 761773. DOI: 10.1080/09500340408235551.

[6] Khonina SN, Kazanskiy NL, Soifer VA. Optical vortices in a fiber: mode division multiplexing and multimode self-imaging. In Book: Yasin M, Harun SW, Arof H, eds. Recent progress in optical fiber research. Rijeka, Croatia: InTech Open; 2012; Ch 15: 327-352. DOI: 10.5772/28067.

[7] Gbur GJ. Singular optics. New York: CRC Press; 2016. ISBN: 978-1-4665-8077-0.

[8] Lavery MPJ, Berkhout GCG, Courtial J, Padgett MJ. Measurement of the light orbital angular momentum spectrum using an optical geometric transformation. J Opt 2011; 13(6): 064006.

[9] Berkhout GCG, Lavery MPJ, Courtial J, Beijersbergen MW, Padgett MJ. Efficient sorting of orbital angular momentum states of light. Phys Rev Lett 2010; 105: 153601.

[10] Volyar A, Bretsko M, Akimova Ya, Egorov Yu. Measurement of the vortex spectrum in a vortex-beam array without cuts and gluing of the wavefront. Opt Lett 2018; 43(22): 5635-5638. DOI: 10.1364/OL.43.005635.

[11] Alperin SN, Niederiter RD, Gopinath JT, Siements KE. Quantitative measurement of the orbital angular momentum of light with a single, stationary lens. Opt Lett 2016; 41(21): 5019-5022. DOI: 10.1364/OL.41.005019.

[12] Alpeiin SN, Siements KE. Angular momentum of topological structured darkness. Phys Rev Lett 2017; 119: 203902.

[13] Kotlyar VV, Kovalev AA, Porfirev AP. Methods for determining the orbital angular momentum of a laser beam. Computer Optics 2019; 43(1): 42-53. DOI: 10.18287/2412-6179-2019-43-1-42-53.

[14] Flusser J, Tomás S, Zitová B. Moments and moment invariants in pattern recognition. New York: John Wiley & Sons; 2009.

[15] Fu S, Gao Ch. Influences of atmospheric turbulence effects on the orbital angular momentum spectra of vortex beams. Photon Res 2016; 4(5): B1-B4.

[16] Kaiser T, Flamm D, Schröter S, Duparré M. Complete modal decomposition for optical fibers using CGH-based correlation filters. Opt Express 2009; 17: 9347-9356.

[17] Kotlyar VV, Kovalev AA, Porfirev AP. Astigmatic transforms of an optical vortex for measurement of its topologi-cal charge. Appl Opt 2017; 56(14): 4095-4104. DOI: 10.1364/A0.56.004095.

[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.

[19] Prudnikov AP, Brychkov YA, Marichev OI. Integrals and series, elementary functions. New York: Gordon and Breach; 1986.

[20] Volyar AV, Bretsko MV, Akimova YaE, Egorov YuA. Beyond the light intensity or intensity moments and measurements of the vortex spectrum in complex light beams. Computer Optics 2018; 42(5): 736-743. DOI: 10.18287/2412-6179-2017-42-5-736-743.

[21] Volyar AV, Bretsko MV, Akimova YaE, Egorov Yu.A. Avalanche instability of the orbital angular momentum higher order optical vortices. Computer Optics 2019; 43(1): 14-24. DOI: 10.18287/2412-6179-2019-43-1-14-24.

[22] Volyar AV, Egorov YuA. Super pulses of orbital angular momentum in fractional-order spiroid vortex beams. Opt Lett 2018; 43(1): 74-77. DOI: 10.1364/OL.43.000074.

[23] Kotlyar VV, Kovalev AA, Porfirev AP. Hermite-Gaussian laser beams with orbital angular momentum. Computer Optics 2014; 38(4): 651-656.

Authors' information

Alexander Vladimirovich Volyar - born in 1949. A graduate of the Simferopol State University (1972), and after graduation - his employee. Doctor of Physical and Mathematical Sciences (1992). Since 1993, he has been a professor, since 1999, he has become the head of the General Physics department. The sphere of scientific interests is optics, laser physics. E-mail: volyar.singular.optics@smail.com .

Mikhail Vladimirovich Bretsko - born in 1994. Since 2012, a student of the V.I. Vernadsky Crimean Federal University. At the moment, the first-year postgraduate student of the General Physics department, the V.I. Vernadsky Phys-ico-Technical Institute of KFU. Scope of scientific interests - singular optics, processing of graphic images, programming. E-mail: mikhailbretcko4@smail.com .

Yana Evgenievna Akimova - born in 1994. Since 2012, a student of the V.I. Vernadsky Crimean Federal University. At the moment, the first-year postgraduate student of the General Physics department, the V.I. Vernadsky Physico-Technical Institute of KFU. Scope of scientific interests - singular optics, difraction optics. E-mail: ak1mova.yana@yandex.ru .

Yuriy Alexandrovich Egorov - born in 1975. A graduate of Taurida National University (2001). Received the degree of Candidate of Physical and Mathematical Sciences in the field of Optics and Laser Physics at the Taurida National V.I. Vernadsky University (Crimean Federal University is now) in 2005. Since 2006, he has been a lecturer at the Physics and Technology Institute of V.I. Vernadsky Crimean Federal University. The sphere of scientific interests is singular optics, diffraction optics. E-mail: vuriyesorov@cfuv.ru .

Received 28 March, 2019. The final version - 11 June, 2019.