CC BY
42
ФАЗОВЫЙ ФИЛЬТР ДЛЯ СЕЛЕКЦИИ УГЛОВЫХ ГАРМОНИК
В.В. Котляр, С.Н. Хонина Институт систем обработки изображений РАН
Аннотация
В данной работе рассматриваются ДОЭ, предназначенные для селекции угловых гармоник, присутствующих в лазерном пучке. Угловая гармоника - это комплексная экспонента, аргумент которой линейно зависит от полярного угла, причем коэффициент пропорциональности - целое число - ее номер. Рассчитан фазовый 8-ми канальный ДОЭ для селекции угловых гармоник с номерами ±1, ±2, ±3, ±4 и проведено численное моделирование для пучков Бесселя и Гаусса-Лагерра.
Введение
Известны световые лазерные моды, которые формируются в лазерных резонаторах и световых волокнах, и фаза которых линейно зависит от полярного угла в цилиндрической системе координат. Это моды Гаусса-Лагерра, Бесселя. В поперечном сечении распределения интенсивности таких мод имеются точки с нулевой интенсивностью. Световое поле в окрестности этих нулей интенсивности ведет себя как гидродинамическая воронка.
Комплексную амплитуду любого поля можно разложить в цилиндрических координатах по угловым гармоникам. Угловая гармоника (УГ) - это комплексная экспонента, аргумент которой линейно зависит от полярного угла, причем коэффициент пропорциональности - целое число - номер УГ.
УГ описывает сингулярность фазы волнового фронта, т.к. при обходе вокруг центра полярной системы координат фаза приобретает набег равный 2т, п - номер УГ.
В работах [1-7] УГ генерировались с помощью специальных оптических элементов цифровой голографии [1,2,6], амплитудных дифракционных решеток [5,7] и фазовых дифракционных оптических элементов (ДОЭ) [3,4].
УГ генерировались в составе мод Бесселя [8] и Гаусса-Лагерра на выходе резонатора лазера [9], в составе квазимодовых пучков [10,11]. Линейная комбинация из УГ формировалась с помощью амплитудных и фазовых оптических элементов для демонстрации эффектов вращения световых пучков [12-16]: бесселевых многомодовых пучков [12], многомодовых пучков Гаусса-Лагерра [13-15] и спиральных гауссовых пучков [16].
Оптическое разложение светового поля по ортогональному базису, согласованному с модами Га-усса-Лагерра и Бесселя рассматривалось в [17].
В данной работе рассматриваются ДОЭ, предназначенные для селекции УГ, присутствующих в лазерном пучке. В [18] такая задача была решена с помощью амплитудного 3-х канального фильтра, который представляет из себя дифракционную решетку со сбоем и выделяет только УГ ±1-го порядков. В данной работе рассчитан фазовый 8-ми канальный ДОЭ для селекции УГ с номерами ±1, ±2, ±3, ±4 и проведено численное моделирование для различных модовых пучков.
1. Теоретическое основание
Комплексную амплитуду светового поля /(Т,ф) в цилиндрических координатах {х=гсо$ф, у=г$>іпф } можно аппроксимировать усеченным функциональным рядом ортогональных базисных функций угловых гармоник [17,18] ехр(іпф), п=0,±1,±2,... :
N
где
/ (г,р) = £ /п (г )ехр(іпр),
п=-N
1 2п
/п (г) = — 11(г,ф)ехр(-іпф)&ф .
(1)
(2)
Пропускание многоканального фильтра, осуществляющего оптическое разложение светового поля (1) по базису угловых гармоник имеет вид:
т(г,р) = £ сп
т=-М
х ехр
- тф+ і -/ гРт СОЪ(р-вт )
(3)
где £=2 ПА - волновое число света с длиной волны А, / - фокусное расстояние линзы, выполняющей Фу-рье-преобразование и в плоскости пространственного спектра которой формируется комплексная амплитуда:
.£ К 2п
р(Р,в) =П 11 /(г, ^)г(г, Р) х 2/ 0 0
£
- і/грсо^ф-в)
х ехр
і£ М N ;т-п
Г d Г dф =
£ £ С т ~
2/ т=-Мп=-N
2п
(4)
х ехр[- і (т - п)$т р 0)]| ]п (г)--т-п Х 0
х 0ТГ-\/ Р2 +Р2т - 2РРт СО$(в-вт ) ^ г d Г,
где
X
х
3m (P,e) = arCtg
Pm sinem -Psine
Pm C0sem -PC0se Jm(x) - функция Бесселя первого рода m-го порядка.
(5)
В точке Фурье-плоскости с полярными координатами (Рт, вт) амплитуда светового поля (4) будет иметь вид:
F(Pm , вm ) =
ikC„
(2п)2 f
j fm (r) r d
r d r +
ik
iCm'Hy
(2п)2 f
m’^m
x exp
R f к I----------------------------------------'Л
p[i(m'-m)3m'(Pm',em')])j fm (r)Jm'-m Ij^tJpI + Pm - 2PmPm C0s(0 m - в m' ) J r d r
(6)
x
Второе слагаемое в ур. (6) много меньше первого, т.к. аргумент функции Бесселя в несколько раз больше 2п
к
Рт + Рт' - 2ртрт'С05(@т -@т') >> 2п .
Это неравенство следует из того, что точки (рт, вт) и (рт-, вт) принадлежат разным дифракционным порядкам (диаметр которых примерно равен Д/Ж), пространственно разделенным в Фурье-плоскости. С учетом этого замечания вместо ур. (6) запишем приближенное выражение для интенсивности в центре т-го порядка:
1 (Pm ,вm ) = |F(Pm ,вm ^
j fm (r) r d)
/ \ 2
f л
4п2 f.
(7)
Может оказаться, что интеграл в ур. (7) равен нулю при /т(г)^0, тогда датчик не определит наличие в пучке т-ой УГ, хотя она будет присутствовать в поле с амплитудой /(г, р). Чтобы избежать этого случая можно регистрировать последовательно несколько распределений интенсивности 1(рт, вт), #=1,2,... при разных радиусах диафрагмы Я#.
Если fm(r)=0, т.е. в световом поле f(r,p) не присутствует гармоника exp(inp), то в центре дифракционного порядка при (p, e=(Pm, em) значение интенсивности, как следует из ур. (7), будет равно нулю (с учетом замечаний, относящихся ко второму слагаемому в ур. (6)).
Для известных световых мод, описываемых функциями Бесселя или Гаусса-Лагерра, радиальными многочленами Цернике, функции из ур. (2) примут конкретный аналитический вид:
fn (r) = Jn (arX
fn (r) = exp
V 2w У
vw2 У
(S)
fn (r) = Rn (r)
где (г) - присоединенный многочлен Лагерра,
К (г) - радиальный многочлен Цернике.
Если имеется радиально-симметричное поле /(г, (р)=/о(г), то вместо ур. (6) получим:
F (P,e) =-
ik
(2п)2 f
I Cm (-i)m exp[im3m (P,e)]j ]0 (r) Jm f fr-JP2 + pI - 2PPm C0s(P - вm ) | r dr .
m 0 Vf )
(9)
В разностных координатах, связанных с центром т-го дифракционного порядка ДРт = |Р - Рт |, где Р = (р,9) и Рт = (Рт,вт), интенсивность будет пропорциональна преобразованию Ханкеля т-го порядка:
1 m
kCm
(2п)2 f
j f0(r ) Jm [V frP j r d r
(10)
То есть фильтр (3) для радиально-симметричного поля выполняет параллельно несколько преобразований Ханкеля заданных порядков.
3. Численное моделирование.
При расчете 8-ми порядкового фильтра (3) пространственные несущие (рт, вт), соответствующие каждому из порядков, были выбраны так, чтобы получить бинарную фазу. На рис. 1а показана бинарная (белый - 0, черный - п) фаза 8-ми порядкового фильтра (3). Причем коэффициенты Сп были подобраны таким образом, чтобы модуль функции Тг, р) был близок к постоянному значению. На рис. 1б дана диаграмма дифракционных порядков с номерами УГ ±1, ±2, ±3, ±4 (радиус Я=6 мм) порядков (б).
r
n
L
2
2
(б)
4 3 2
-1 1
- -3 4 -
Рис. 1. Бинарная фаза 8-ми порядкового фильтра УГ (а) и диаграмма расположения дифракционных (рис. 2в).
На рис. 2 показан результат работы фильтра с фазой, приведенной на рис. 1а, при освещении его модовым пучком Гаусса-Лагерра, состоящим из од-
ной моды с номером (1,п)=(4,0) ехр
.„2 Л
2^
С ,.2 Л
Т
V У
(рис. 2а) и композиции мод с номерами (1,1)+(4,-1)
ехр
..2 Л
2^
С 2 Л С 2 Л
Г 1 Г
ехр(/» + Т4
ехр(-/»[>.
Перетяжка гауссова пучка равна те=1 мм. На рис. 2б,г показаны соответствующие картины распределения интенсивности в фокальной плоскости.
Т.к. пучок Гаусса-Лагерра с модой (4,0) является радиально-симметричным полем, то в каждом из порядков на рис.2б можно видеть результат соответствующего преобразования Ханкеля. Пучок Гаусса-Лагерра с композицией мод (1,1)+(4,-1) имеет УГ с номерами т=±1, которые и были обнаружены фильтром: стрелками на рис. 2г отмечены дифракционные порядки, в центральных точках которых имеются не нулевые значения интенсивности (при т=±1), во всех других порядках на рис. 2б,г в центрах наблюдается нулевое значение интенсивности.
Рис. 2. Анализируемые пучки Гаусса-Лагерра: мода (4,0) (а) и композиция мод (1,1)+(4,-1) (в) и соответствующее им распределение интенсивности в фокальной плоскости после фильтра (б) и (г).
рис. 3б,г стрелками показаны дифракционные порядки, в центрах которых имеется ненулевая интенсивность. В центрах всех других порядков на рис. 3 б,г наблюдается нулевая интенсивность.
О О' О ▲
О О
Ш Ш (б>
▲
О о о
►
Рис. 3. Анализируемые пучки Бесселя: мода J3 (5r) ехр(/3ф) (а) и композиция мод J4(3r)еХр(/4ф) + J4(3r) еХр(—/4ф) (в) и соответствующее им распределение интенсивности в фокальной
плоскости после фильтра (б) и (г).
Видно, что фильтр работает независимо от подаваемого на фильтр размера пучка, при этом меняется только размер дифракционного порядка. Хотя для того, чтобы приближенное выражение (7) было более точным, нужно чтобы как можно дальше были разнесены.
Заключение.
В данной работе получены следующие результаты:
• Рассмотрен фазовый ДОЭ для обнаружения угловых гармоник, содержащихся в лазерном пучке.
• Проведено теоретическое рассмотрение и численное моделирование работы данного ДОЭ, и показано, что он одинаково эффективно обнаруживает угловые гармоники в многомодовых пучках Бесселя и Гаусса-Лагерра.
• В последующих работах планируется изготовление данного ДОЭ и его экспериментальная апробация.
Благодарность
Работа поддержана Российским фондом фундаментальных исследований
N 98-01-00894, 99-01-39012, 00-15-96114, 0001-00031.
Литература
1. Березный А.Е., Прохоров А.М., Сисакян И.И., Сойфер В.А., Бессель-оптика, ДАН СССР, 274(4), 802-805 (1984)
2. Vasara A., Turunen J., Friberg A.T., Realization of general nondiffracting beams with computergenerated holograms, J. Opt. Soc. Am. 6, 1748-1754 (1989)
3. Khonina S.N., Kotlyar V.V., Soifer V.A., Shinkaryev
M.V., Uspleniev G.V., Trochoson, Optics ^ммит-cations 91(3,4), 158-162 (1992)
4. Khonina S.N., Kotlyar V.V., Shinkaryev M.V., Soifer
V.A., Uspleniev G.V., The phase rotor filter, J. Mod. Opt. 39(5), 1147-1154 (1992)
На рис. 3 показан результат работы фильтра (рис.1а, радиус Я=5 мм) при освещении его модо-вым пучком Бесселя, состоящим из одной моды З3(5г )ехр(і 3ф) (рис. 3а) и композиции мод
З4(3г)ехр(і4ф) + З4(3г)ехр(-і 4ф) (рис. 3в). На
5. Bazhenov V.Yu., Soskin M.S., Vasnetsov M.V., Screw dislocations in light wavefronts, J. Mod. Opt. 39(5), 985-990 (1992)
6. Heckenberg N.R., McDuff R., Smith C.P, White A.G., Generation of optical phase singularities by computer-generated holograms, Opt. Lett., 17, 221223 (1992)
7. Basistiy I.V., Bazhenov V.Yu., Soskin M.S., Vasnet-
sov M.V., Optics of light beams with screw dislocations, Opt. Commun. 103, 422-428 (1993)
8. Paakkonen P., Turunen J., Resonators with Bessel-Gauss modes, Opt. Commun. 156, 359-366 (1998)
9. Abramochkin E., Losersky N., Volostnikov V., Generation of spiral-type laser beams, Opt. Commun. 141, 59-64 (1997).
10. Spektor B., Piestun R., Shamir J., Dark beams with a constant notch. Opt. Lett. 21(7), 456-458 (1996)
11. Piestun R., Schechner Y.V., Shamir J., Self-imaging with finite energy, Opt. Lett. 22(4), 200-202 (1997)
12. Paakkonen P., Lautanen J., Honkanen M., Kuittinen M., Turunen J., Khonina S.N., Kotlyar V.V., Soifer V.A., Friberg A.T., Rotating optical fields: experimental demonstration with diffractive optics, J. Mod. Opt. 45(11), 2355-2370 (1998)
13. Khonina S.N., Kotlyar V.V., Soifer V.A., Honkanen M., Lautanen J., Turunen J., Generation of rotating Gauss-Laguerre modes with binary-phase diffractive optics, J. Mod. Opt., 42(2), 227-238 (1999)
14. Schechner Y.V., Piestun R., Shamir J., Wave propagation with rotating intensity distributions, Physical Review E 54(1), 50-53 (1996)
15. Piestun R., Schechner Y.V., Shamir J., Self-imaging with finite energy, Opt. Lett. 22(4), 200-202 (1997)
16. Abramochkin E., Volostnikov V., Spiral-type beams: optical and quantum aspects, Opt. Commun. 125, 302-323 (1996)
17. Kotlyar V.V., Khonina S.N., and Soifer V.A., Light field decomposition in angular harmonics by means of diffractive optics, J. Mod. Opt., 45(7), 1495-1506 (1998)
18. Trapani P.D., Berzanskis A., Minardi S., Sapone S., Chinaglia W., Observation of optical vorticies and J0 Bessel-like beams in quantum-noise parametric amplification, Phys. Rev. Lett., 81(23), 5133-5151 (1998)
