Параксиальные гипергеометрические лазерные пучки с особенностью в центре перетяжки

ПАРАКСИАЛЬНЫЕ ГИПЕРГЕОМЕТРИЧЕСКИЕ ЛАЗЕРНЫЕ ПУЧКИ С ОСОБЕННОСТЬЮ В ЦЕНТРЕ ПЕРЕТЯЖКИ

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

Аннотация

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

Введение

В современных научных исследованиях интерес к различным типам лазерных пучков заметно возрос. Продолжаются исследования хорошо известных мод Лагерра-Гаусса [1]. В [2] с помощью функции Вигнера рассмотрены непараксиальные моды Ла-герра-Гаусса. Продолжают исследоваться элегантные пучки Лагерра-Гаусса [3] и Эрмита-Гаусса [4]. Эти пучки описываются соответствующими полиномами с комплексными аргументами. В [5, 6] исследовались разными методами эллиптические пучки Лагерра-Гаусса. Лазерный пучок с эллиптической симметрией можно сформировать при наклонном падении плоской волны, например, на аксикон. В [7] исследуются лазерные пучки, сформированные при наклонном падении на аксикон. Интерес также вызывает линейная комбинация лазерных мод. Так, в [8] показано, что с помощью подбора мод в осевой линейной комбинации можно получить лазерные пучки с вращающейся поперечной интенсивностью, но не обладающие орбитальным угловым моментом. Можно также сформировать лазерный вихревой пучок, картина дифракции которого будет представлять собой только одно светлое кольцо, а все периферийные кольца будут подавлены [9].

В недавних работах авторов [10, 11] рассматривалось новое семейство параксиальных лазерных пучков - гипергеометрические моды. Эти моды близки к известным модам Бесселя [12, 13], они также обладают бесконечной энергией и могут быть сформированы с помощью ДОЭ только приближенно. В данной работе рассматривается некоторое обобщение гипергеометрических мод - гипергеометрические пучки, которые переходят в моды при замене исходного Гауссового пучка плоской волной или при стремлении радиуса перетяжки гауссового пучка к бесконечности. Показано, что комплексная амплитуда у некоторых из таких пучков в дальней зоне пропорциональна модифицированной функции Бесселя целого или полуцелого порядков.

1. Гипергеометрические лазерные пучки

Рассмотрим световое поле с начальной функцией

комплексного пропускания вида:

1^2 Л

/ Ч 1 ^ г г

Еуп

2п г 2а ^

где (г,ф) - полярные координаты в начальной плоскости (z = 0), w и у - действительные параметры логарифмического аксикона, а - радиус перетяжки гауссова пучка, n - целый порядок спиральной фазовой пластинки. Комплексная амплитуда (1) описывает световое поле с бесконечной энергией и с особенностью при r = 0 . Несмотря на это, в любой другой поперечной плоскости на расстоянии z от начальной плоскости комплексная амплитуда светового поля, порожденная функцией (1), уже не будет иметь особенности и будет конечной.

При параксиальном распространении светового поля (1) его комплексная амплитуда на расстоянии z будет определяться преобразованием Френеля, которое в полярных координатах имеет вид:

E (рв z ) = - П ИЕ (r -0 )х

Г ik R 1

хexpI — [p2 + r2 -2prcos(--в) j-rdrd-,

где (p, в) - полярные координаты в плоскости, перпендикулярной оптической оси, k - волновое число.

Подставляя (1) в (2) и используя справочный интеграл [14]

i

х" 1 exp

(-px2 )v (cx )dx = cvp

xrl^lr-1 (v + l)i Fi

v + a

„2 Л

v +1,--

4 p

получим следующее выражение для комплексной амплитуды в плоскости z :

Еуп

п+1 —

1-1/

п.п+2 I I _

2 п'.пу г

р *

<ехр

('кр

+ 'пв

'к

2а2 2г

п+1+'/ 2

*г| ^^ 11 Ъ

п +1 + 'у

п +1, —

кР

27

2 А

1

¡к

2а2 2г

где 1 (а, Ь, х) - вырожденная или конфлуэнтная гипергеометрическая функция (функция Куммера):

а I х

т Л

т

т=0 (Ь )т т!

(а)т =Г(а + т)/г(а) - символ Похгаммера. Из

выражения (4) видно, что комплексная амплитуда светового поля на расстоянии 7 > 0 имеет конечный модуль в каждой точке р .

При большом удалении от начальной плоскости (7 >> ка2) выражение (4) упростится:

(-1 )п+1 к А

Еуп (р,в,г >> ка2;а) = -

р ехр

('кр2 2г

+ 'пв

1

оп+2 „ I

2 п.п

п+1+'у

2

2а2

*Г| \ 1 ^

п +1 + 'у

п +1, -

^ кра^

V27.

2 А

Частным случаем пучков (6) являются пучки при у = 0 . Из (6) следует, что комплексная амплитуда таких пучков пропорциональна функции Бесселя второго рода (модифицированной) с целым и полу целым порядками:

Е0п (р,в,7 >> ка2;а) = КТ 7 'кр2 (кра

п+1 к а—

<ехр

2 г

2 г

+ ¡пв

г 242ж кра

2 г

2. Гипергеометрические моды

При устремлении радиуса перетяжки гауссово-го пучка к бесконечности (а ^ж, т.е. Гауссов пучок переходит в плоскую волну) выражение (4) существенно упрощается:

п +1 + 'ул

Еуп (р,в :

м>-прп 2жп\

'к—' 2 г

п+1-'у

.2 А

ехр

( Ч 2 'кр

2г

-¡пв

<1 ъ

2

п +1 + 'у , ¡кр

-п +1,--—

2 2г

Для функции Куммера известно соотношение:

1 ъ

Тогда (8) преобразуется:

ЕГп (р,в, г) =

— р [ п +1 + 'У | ,

-—г|-- Iехр('.

2пп\ 1 ^ 1 ^

^ 'к—2 А

п+1-'/ 2

2 г

1 ъ

^ п +1 - 'у , 'кр2 2

--, п +

2 2г

Следует отметить, что справочный интеграл (3) определен для случаев, когда Яе (р) > 0 , что не соблюдается при а ^ж. Однако при подстановке (10) в параксиальное уравнение распространения света

(

д д

2

2'к — +-+--+

2

1 д 1 д

дг др2 р др р2 дв2

и

получается верное тождество, а именно, дифференциальное уравнение Куммера:

Ж

д2

д

—у + (Ь-%)~--а

дХ дХ

Ъ

где

а =

п +1 - 'у 2 '

Ь = п +1, 'кр2

х = -

2 г

Гипергеометрическая мода (10) формируется из начальной комплексной амплитуды (г = 0) с особенностью в начале координат:

Е

2п г I — )

3. Гипергеометрические моды при у = 0

Рассмотрим однопараметрическое семейство гипергеометрических мод (10) при условии, что у = 0 .

В этом случае воспользуемся известным соотношением для конфлюэнтной функции [15]

п 1 п 'кр1

1 11 2 2 2 2г

ехр

21 4г

4г Д 8г

После подстановки (15) в (10) получается

Е°и (р,в, г)=— ( П Г ехр

<ехр| '-кр + 'пв 13 ' кр

-' 4

4 г

4 г

п

Из (16) видно, что комплексная амплитуда гипергеометрических мод при у = 0 пропорциональна

функции Бесселя первого рода (сравните с уравнением (7), где появляется функция Бесселя второго рода) целого и полу целого порядков. Заметим, что все рассмотренные выше поля (и гипергеометрические пучки (6) и моды (10) зависят от квадрата радиальной переменной. Это означает, что с ростом переменной р будет уменьшаться расстояние между соседними локальными максимумами или минимумами. Неограниченное увеличение пространственной частоты картины дифракции является следствием особенности (амплитудной и фазовой) комплексной амплитуды в центре начальной плоскости (в центре перетяжки). В реальности для формирования таких лазерных пучков следует использовать кольцевую диафрагму, «блокирующую» особенность в центре координат. При этом, правда, гипергеометрические пучки можно будет сформировать только приближенно, но зато пространственная частота картины дифракции будет конечной величиной.

Рассмотрим первую вихревую моду при п = 1. Тогда вместо (16) получим комплексную амплитуду:

Em (р,в, z) = — -JIz sin J Р | exp(i<9) • пр l 21z ,

Нули комплексной амплитуды (17) будут иметь координаты p^ =^2mXz . Первый максимум амплитуды находится из уравнения 2 x = tg (x) и имеет координату вблизи точки pmax <-Jxz . Заметим, что

расстояние между двумя соседними нулями интенсивности поля (17) с номерами m и m +1 при m >> 1 равно выражению Дрт (2m). То есть

пространственная частота картины дифракции стремиться к бесконечности при p^-». Устранить это можно, ограничив исходное поле(14) кольцевой диафрагмой, блокирующей особенность в центре.

4. Обобщенные гипергеометрические моды

Рассмотрим функцию следующего вида (более общего, чем (10):

Е

где F - некоторая функция.

После подстановки функции (18) в параксиальное уравнение распространения света (11) получается следующее соотношение:

s2 m 2p2m-2 z21F " +

+ [

+ [P (P -1) p- + pp"2 - n2p~2 + 2ikqz- ] F = 0.

При p = ±n уравнение упрощается. Для определенности рассмотрим случай, когда p = + n :

2 2 2т-2 21 т-'ГГ ,

Я т р 2 Г +

+ \_яш

+2^ = 0.

Чтобы получить из уравнения (20) уравнение Куммера, нужно чтобы коэффициент при производной второго порядка равнялся аргументу функции Г. Для этого разделим обе части уравнения на

8т2рт- 2 :

m I 2n + m 2iklp2z 1 sp z F +|-+-—2-

mm

F ' +

+ mqp2-mz l 1 f = 0. sm2

Далее нужно, чтобы коэффициент при производной первого порядка был равен разности константы и аргумента функции Г, это возможно только в случае, когда

m = 2, l = -1,

2ik1 ik = m2 ~ 2. Тогда вместо (21) получим уравнение вида:

у p2z-1F" +

Получено уравнение Куммера, в котором

a = -q, b = n +1.

Решением такого уравнения является конфлу-энтная функция:

j F I a, n +1

ikp2

т.е. решением параксиального уравнения распространения (11) является следующая функция:

Еа +п

При p = -n получается

Ea -„ (р,в, z) = p-nz-a exp(inff) 1F1I a,1 -

ikp2 ' 2 z

.

Объединяя (26) и (27) в одно выражение, получится семейство функций

Ean (р,в, z) = pnz-a exp(±тв) 1F11 a,n +1

ikp2 ' 2z

Комплексную амплитуду (28) можно рассматривать как обобщенный параксиальный гипергеометрический пучок, который переходит в модовые пучки при условии а = ( п +1 + ¡у)/2.

5. Численное моделирование

Формула (10) описывает комплексную амплитуду для идеальной гипергеометрической моды, сформированной с помощью начального поля вида (14). На практике это поле должно быть ограничено в размерах. Модуль функции комплексного пропускания также должен быть ограничен единицей. На рис. 1а и б показаны рассчитанные по формуле (14) радиальные распределения амплитуды в плоскости г = 0 : неограниченное (1а) и ограниченное кольцевой диафрагмой с радиусами Я1 = 0,2 мм и Я2 = 10 мм. На рис. 1в и г показаны рассчитанные соответствующие радиальные распределения амплитуды в плоскости г = 2000 мм: по формуле (10) (1в) и с помощью преобразования Френеля (1г). Параметры расчета: порядок СФП п = 10, параметры логарифмического акси-

кона — = 0,5 мм и у = 0 , длина волны X = 514,5 нм, дискретизация Др = 0,16 мм.

На рис. 2 показано рассчитанное радиальное распределение амплитуды на расстоянии г = 2000 мм от исходного поля (14) для двух различных значений параметра логарифмического ак-сикона у. Видно, что с ростом у радиус первого (главного) кольца увеличивается.

На рис. 3 показано радиальное распределение амплитуды на расстоянии г = 2000 мм от исходного поля (1) для двух различных значений радиуса перетяжки гауссового пучка а при у = 0 . Видно, что с ростом а уменьшается контраст дифракционной картины.

На рис. 2 и 3 остальные параметры расчета такие же, что и на рис. 1.

1,0

а)

в)

0,10*

0,08

0,06

0,04

0,02 )

0

\Еоло(г,<р,0)\ 1,0-

10 р, мм

Рис. 1. Радиальные распределения модуля комплексной амплитуды в плоскостях г = 0 мм (а и б) и г = 2000 мм (в и г) для светового поля (14) с бесконечной апертурой (а и в) и ограниченного кольцевой диафрагмой с радиусами ^ = 0,2 мм

и ^2 = 10 мм

0,07 0,06 0,05 0,04 0,03 0,02 0,01

а)

о

|Я5оло(рДг)| 0,05 0,04 0,03 0,02 0,01

8 10 р, мм

б)

0

8 10 р, мм

Рис

2. Радиальные распределения модуля комплексной амплитуды для начального поля

а)

ojo1

0,08

0,06

0,04

0,02 J

0

ll^

4

8 10 р, мм

б)

0,08

0,07

0,06

0,05

0,04

0,03

0,02

0,01 J

0

10 p, мм

Рис. 3. Радиальные распределения модуля комплексной амплитуды от исходного поля

Заключение

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

Гипергеометрические моды обладают одной особенностью: пространственная частота картины дифракции асимптотически стремиться к бесконечности.

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

Работа выполнена при поддержке российско-американской программы «Фундаментальные исследования и высшее образование» (грант CRDF RUX0-014-SA-06), а также грантов РФФИ 05-0S-5029S и 07-07-97600.

Литература

1. Gao M., Gao C., Lin Z. Generation and application of the twisted beam with orbital angular momentum, Chin.Opt. Lett., 2007. V.5. N.2. P.S9-92.

2. Duan K., Lu B. Application of the Winer distribution function to complex-argument Hermite-and laguerre-

Gaussian beams beyond the paraxial approximation // Opt. Las. Tech., 2007. V.39. P.110-115.

3. Seshardi S.R. Self-interation and mutual interaction of complex-argument Laguerre-Gauss beams // Opt. Lett., 2006. V.31. N.5. P.619-621.

4. Kostenbauder A., Sun Y., Siegman A.E. Eigenmode expansions using biorthogonal functions: complex-valued Hermite-Gaussians: reply to comment // J. Opt. Soc. Am. A, 2006.V. 23. N.6. P.1528-1529.

5. Kotlyar V.V., Khonina S.N., Almazov A.A., Soifer V.A., Jefimovs K., Turunen J. Elliptic Laguerre-Gaussian beams // J. Opt. Soc. Am. A, 2006. V.23. N.1. P.43-56.

6. Mei Z., Gu J., Zhao D. The elliptical Laguerre-Gaussian beam and its propagation // Optik, 2007. V.118. P.9-12.

7. Burvall A., Kolacz K., Goncharov A.V., Jaroszewicz Z., Dainty C. Lens axicons in oblique illumination // Appl. Opt., 2007. V.46. N.3. P. 312-318.

8. Kotlyar V.V., Khonina S.N., Skidanov R.V., Soifer V.A. Rotation of laser beams with zero of the orbital angular momentum // Opt. Commun., 2007 (in press).

9. Kotlyar V., Kovalev A., Soifer V., Tuvey C., Devis J. Sidelobes contrast reduction for optical vortex beams using a helical axicon // Opt. Lett., 2007. V.32 (in press).

10. Kotlyar V.V., Skidanov R.V., Khonina S.N., Soifer V.A. Hypergeometric modes // Opt. Lett., 2007. V.32. N.7. (in press).

11. Котляр В.В., Скиданов Р.В., Хонина С.Н., Балалаев С.А. Гипергеометрические моды // Компьютерная оптика, Самара, ИСОИ РАН, 2006. № 30. С. 16-22.

12. Durnin J., Miceli J.J., Eberly J.H. Difraction-free beams, Phys. Rev. Lett., 1987. V. 58. P. 1-1501.

13. Khonina S.N., Kotlyar V.V., Skidanov R.V., Soifer V.A., Jefimovs K., Simonen J., Turunen J. Rotation of micropar-ticles with Bessel beams generated by diffractive elements // J. Mod. Opt., 2004. V.51. N. 14. P. 2167-2184.

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

15. Справочник по специальным функциям // Под ред. М. Абрамовица, И. Стигана. - М., Наука, 1979.