Обобщенные эрмитовые световые пучки в свободном пространстве

ОБОБЩЕННЫЕ ЭРМИТОВЫЕ СВЕТОВЫЕ ПУЧКИ В СВОБОДНОМ ПРОСТРАНСТВЕ

Котляр В.В., Сойфер В.А., Хонина С.Н.

Институт систем обработки изображений РАН

Аннотация

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

Введение

Интерес к эрмитовым пучкам или, точнее, модам Гаусса-Эрмита, обусловлен тем, что в свободном пространстве они распространяются без изменения своей структуры, изменяясь только масштабно, а в световых волокнах с параболическим профилем показателя преломления они распространяются и без изменения масштаба.

В работах [1,2] приведены разные аналитические выражения для мод Гаусса-Эрмита. Это послужило толчком к тому, чтобы выяснить степень различия и сходства двух разных форм эрмитовых пучков. Оказалось, что оба типа пучков являются частными случаями более общих модо-вых эрмитовых пучков.

Для того, чтобы сформировать в пространстве такие световые пучки, следует синтезировать на компьютере и изготовить по технологии микролитографии специальные дифракционные оптические элементы. Имеется ряд работ [3-5], по-

священных расчету таких фазовых оптических элементов. В [3] они названы моданами. Для их расчета используется метод Кирка-Джонса кодирования с пространственной несущей частотой [6].

В [4] для расчета оптического элемента, генерирующего одномодовый пучок, используется итеративная процедура и вводится в рассмотрение вспомогательная область.

В [5] расчет ведется на основе итеративной аппроксимации функции пропускания транспаранта конечной суммой ортогональных мод.

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

Приведены также аналитические выражения, описывающие дифракцию Фраунгофера на таких фазовых элементах.

1. Два типа эрмитовых пучков

Известно, что частным решением параболического уравнения распространения

(

э э2

2ik---------+

+

э

2 Л

V

Эz Эx Эx2

E (x, у, z) = О

У

где к - волновое число света, (х,у) - поперечные и г - продольная координаты в пространстве, Е(х,у,г) - медленно меняющаяся вдоль г комплексная амплитуда светового поля, являются модовые функции вида 0:

-| т+п+2

Юо

e1 (x, у, z ) =

x exp

2 2 x2 + у2

со2

w(z). я

x

я

w

у

где

о( z) = оО 1+ —

О

k оО

z0 =

2

Нп(х) - многочлен Эрмита, Юо — минимальный радиус пучка при г = 0, 70 —конфокальный параметр.

Функция Е при изменении расстояния г изменяется только масштабно, сохраняя свой вид. Поэтому решение (2) уравнения (1) можно назвать световой эрмитовой модой. Заметим, что моды (2) взаимно не ортогональны, что следует из справочного интеграла 0:

j e-2x яп (x)ят (x)d x =

m+n m+n-1 /

= (-i^T^-T- rf m + n +1

2

где Г(х) - гамма-функция.

С другой стороны, известны взаимно ортогональные моды Гаусса-Эрмита, которые также являются частными решениями уравнения (1) [2]:

E(mn (x, у, z) = exp [i (n + m + 1)h z )]x

x exp

a

x я.

2R

л

я

az

V ' ' У

exp

л

22 x2 + у 2

a 2 (z

где

az

V ' ' У

a( z):

w,

C 2\2

z 2 2

1+ 72

V z0 J

h( z) = arctg —

R(z ) =

О

R(z) - радиус кривизны модового пучка.

Функции Е и E^l описывают разные световые пучки, так как порождаются разными граничными условиями. При z=0 вместо (2) и (6) получим:

Eml (Х >”,0) = exp

22 x 2 + у 2

Оо

x я.

V Оо У

я

f j Л

V Оо У

j = 1,2

q < 1, j = 1 где qj = U j = 2-

Эти два типа модовых эрмитовых пучков отличаются также разным видом дифракции Фраунгофера. Световое поле Е(тП является собственной функцией оператора Фурье [8]:

¥ I \

j e-x2 яп Уїx)e- 2xXd x =

= 4P

Это свойство мод Гаусса-Эрмита сохранять свою структуру не только в ближней зоне дифракции Френеля, но и в дальней зоне Фраунгофера (или в фокальной плоскости линзы) используется для эффективного ввода лазерного излучения в световые волокна [3].

Световое поле Е , сохраняя свою структуру в зоне дифракции Френеля, в дальней зоне «вырождается» в поле с нулевой интенсивностью в центре. Это следует из выражения:

{е-х2 Нп

2. Обобщенные эрмитовые пучки

Можно показать, что функции Е и Е из уравнений (2) и (6) являются частным случаем более общего решения уравнения (1), которое имеет вид:

Е<тП (Х ^ Г) = Атп (г) ехр[- В(г)(х2 + У 2 )1Х

x я„\ * \ я 1 у

nl C

где

2

x

Атп ( - ) =

, т+ п т+ п 1

іу- 2 , і-

1+ — 1+ —

0=

22

Ю0£

1 +

V

1+у

-о)

-1

+ ку

2 -г

-1

С(-) =

- ) = Юо£

1 +

. - 0)

V

у=1-е" , е - модовый параметр. При е=1 решение (13) переходит в решение Е<тПп, а при е

==№)- ‘

решение

переходит в

Е і2).

Обобщенные эрмитовые световые пучки (13) порождаются с помощью комплексной амплитуды при г=0 вида:

Е<тП ^ У,о)= ЄХР

2 2 X2 + у2

х И,

И

Выражение для комплексной амплитуды светового пучка (13) в дальней зоне следует из соотношения:

п И е -і 2 хХ

V е)

Ґ

-Х2 Ип и 1

п

х

X

где е<1. При е= (л/2) уравнение (18) переходит в уравнение (11). При е=1 используя ассимпто-тику Нп(юХ) ——— > 2п(юХ)п , получим:

(-у)2 И

Ип ЮХ)

( \ X

л/ї-

е2 )

=

переходит в уравнение (12).

3. Фазовый элемент для формирования эрмитова пучка

Для энергетически эффективного формирования эрмитовых пучков обоих типов предлагается в плоскости г=0 располагать фазовый оптический элемент с функцией пропускания:

?тп (^ У) = ^Шп (х)Нт (У)] =

= 5ёп Нп

где sgn(x) - знаковая функция.

Задача поиска оптимального фазового фильтра, формирующего эрмитовые пучки сходна

с аналогичной задачей для формирования бесселевых пучков [9].

Покажем в одномерном случае, что транспарант (19) эффективно формирует эрмитовый пучок п-го порядка. Для этого разложим функцию тп(х) по ортогональным многочленам:

sgn Ип

р=о

где

СРп) = А | sgn[Hn (х)]е х2 Ир (х)а х,

А = (2пп!л/Р) . Для четных номеров п=21 нули

многочлена хк/. Н21(хк>1)=0, будут расположены симметрично относительно точки х=0: х_к1=-хк1. Тогда интеграл в уравнении (21) можно представить в виде суммы интегралов:

{ 5Ип[Н21 (х)]е-х2 Нр (х)а х =

“ 2 I

= | е-х Нр

х .

к=1

хк ,1

л 2

| е - х Ир (х)d х

- хк ,і

С учетом справочного интеграла [7]:

/ е-х2 Нп

о

вместо (22) получим следующее выражение для коэффициентов разложения (21):

)!л/РГ х

С221) =[22"-1(25)!

X (-1)1X (-1)к е -хк,; Н2,-1(хК1) к=1

При этом нечетные коэффициенты равны нулю.

Аналогично можно показать, что для нечетных многочленов Эрмита п=2/+1 имеет место следующее выражение для коэффициентов разложения (21):

С22++1 = [2

= |22" (25 + Ц^Н)1 х

И2,

к=1

где

Н 21+1 (хк ,1 ) = 0,

Н2*(0) = (-1)* • 2* • 3 • 5 • 7...

Четные коэффициенты равны нулю.

Таким образом, освещая оптический элемент с пропусканием (19) коллимированным гауссовым пучком с амплитудой ехр[- (х2 + у2 )у2], получим поле

за элементом в виде суперпозиции эрмитовых пучков, один из которых будет энергетически преобладать среди прочих с эффективностью около 70%:

-

-

+

-

V о у

2

х

х

Юое

Юое

2

- х

е

х

п

2

е

п

х

е- х 72

sgn[ Ип

р= о

4. Дифракция Фраунгофера для транспаранта со знаковой функцией пропускания

Для расчета дифракции Фраунгофера плоской волны на участке [-а,а] оптического элемента с пропусканием sgn[Ип (х)], где а должно

быть больше максимального корня многочлена Ип(х), аналогично работе [1о], удобно использовать представления вида (п=21):

^п[и2і (х)]= гей! х

+ 2(-1)1X (-1)кгесі

к=1

( \ х

V Хк ,1 )

где И21(хк,1 )= о, гесі

( \ х

о, х > хп

Тогда преобразование Фурье от выражения

(27) будет иметь вид:

Р21 (Х) = 3^[Н 21 (х)]}= 2а 8шс(аХ) +

1 ' ,

+ 2

к=1

где

3{/ (х)}= | / (х)ехр(-іхХИ х,

БІПС(х) = 8ІП(х) / х

Параметр диафрагмы а выбирается из энергетических соображений.

Аналогично получаются выражения для нечетных полиномов Эрмита:

^[И 21 +1( х)] =

+ 2(-1)1X (-1)к-1гес1

( \ х

к=1

Кк ,1

sgn х

Р21+1 (Х) = 3^п[и 21+1 (х)]} =

= -и2а со8с(аХ)+ 2(-1)1 X (-1)к х,

[ к=1

х 2хк,і соЄс(хк,іх)}

где СО БС( х ) =

1- СО8х

5. Пример

Пусть требуется сформировать первую моду Гаусса-Эрмита при п=1, т=о. Для этого в плоскости -=о, следует расположить транспарант с

пропусканием Н 1(х)ехр(- х2/2) . Тогда в фокальной плоскости линзы согласно уравнению (11) сформируется световое поле с комплексной амплитудой, пропорциональной выражению:

- гл/лН 1

где Н1( Х) = 2Х.

Если в плоскости 2=0 разместить энергетически более выгодный фазовый элемент с пропусканием sgn[Н 1(х)]= sgn(х), то согласно уравнению (25) в пространстве за фильтром сформируется суперпозиция световых мод:

Sgn х = X С21+1И 2, +1( х) =

,=о

=7Р (Их) -12И 3( х)+1бо И 5( х) -

Освещая такой фазовый элемент коллимированным гауссовым пучком, получим, что согласно уравнению (32) в частотной плоскости сформируется световое поле с амплитудой, пропорциональной выражению:

| sgn(х)е х ^2е lхXd х-

= 4РХс21+1

*=0

= -г еХ2 {Н1 (Х) - — Н3 (Х) + — Н5 (Х) -...

^ 1 ^ 12 160 5

Отношение энергий световых пучков, описываемых вторым и первым слагаемыми в уравнении (33) равно:

|И2 (X) е-Х dX (12)2 |И 12(Х)е-x2dX

Таким образом, в первом слагаемом суммы (33), описывающем первую моду Гаусса-Эрмита, содержится более 80% всей энергии светового пучка.

Еще более эффективным способом формирования первой эрмитовой моды оказывается способ освещения фазового фильтра с пропусканием:

^пН\

не гауссовым пучком, а плоской световой волной. В этом случае в фокальной плоскости линзы сформируется световое поле с амплитудой, описываемой выражением:

а

-а

где функция соБс(х) определена в уравнении (31).

+

а

х

\^оу

х

Взаимная корреляция светового поля (36) и первой эрмитовой моды оценивается выражением:

¥ 2 |[2асоъс(аХ)\Н1 (Х)е-Х ^2dХ

л=-

Для сравнения на рис. 2 показан результат дифракции Фраунгофера для второго типа эрмитовых пучков. На рис. 2а показан сам пучок

„2/

е32)(х) = И3(х)ехр(- х2/2), х є [- 8,8],

И1(Х)е

-Х2/2

^Х .

Їл/Р

4 1 - е-а- )2

Размер апертуры фазового элемента (35) [-а,а] выбирается из условия максимума выражения (37), который достигается при а=2.2. При этом корреляция (37) будет равна ^=0.85.

Таким образом первая мода Гаусса-Эрмита содержится в световом поле (36) с эффективностью около 85%.

6. Численные результаты

При численном моделировании использовались следующие параметры: размерность массивов N=256, область изменения аргумента хе [-8,8], длина волны 1=0.63 |лт, фокусное расстояние /=100 тт.

На рис. 1а показан эрмитовый пучок

Е

На рис. 1 б,в показаны соответственно интенсивность (квадрат модуля) и фаза (аргумент) Фурье-преобразования от функции (38).

а на рис. 2б,в показаны интенсивность и фаза его Фурье-образа, соответственно.

На рис. 3а показан результат бинаризации эрмитового пучка (39): кривая 1 - это функция

^3

а кривая 2 - это функция (39). На рис. 3б,в показаны соответственно интенсивность и фаза Фу-рье-преобразования от бинарной функции (4о) при хє[-3,3]. Диапазон изменения аргумента х функции (4о) выбран таким образом, чтобы корреляция ^, аналогичная определенной в уравнении (37), была максимальной. В этом случае корреляция между функциями интенсивности, показанных на рис. 2б и рис. 3б составляет о.84. Поэтому можно заключить, что 84% световой энергии плоского пучка, освещающего фазовый элемент (4о), будет идти на формирование эрмито-вого пучка (39). На рис. 4 показана зависимость взаимной корреляции ^ от полудлины отрезка [-а,а] для эрмитового пучка третьего порядка п=3. В Таблице приведены оптимальные размеры отрезков [-а, а], на которых заданы бинарные функции sgn[Ип (х)], п = 1,7 и значения взаимной корреляции между Фурье-образами от функций (39) и (4о) при п = 1,7 .

Рис. 1.

3.14

-о.25 -0.12

-1.76

о.13 о.25

Рис. 2.

о

Таблица

Номер моды, n 1 2 3 4 5 6 7

Оптимальный размер апертуры, a (mm) 2,1 2,6 3,0 3,4 3,7 4 4,2

Корреляция, h 0,86 0,84 0,84 0,82 0,82 0,81 0,81

а) б) в)

Рис. 3.

Рис.4

Из Таблицы видно, что с ростом номера эр-митового пучка энергетическая эффективность ^ медленно спадает с 86% до 81% (п<7).

Заключение

В работе получены следующие результаты: получено выражение для обобщенных эрмитовых пучков, сохраняющих свою структуру при распространении в зоне дифракции Френеля (ур. (13)-(16)); получены разложения в ряд по многочленам Эрмита бинарных знаковых функций, которые можно рассматривать как функции пропускания оптических элементов, для эффективного формирования эрмитовых пучков (ур. (20), (24), (25)); получены выражения, описывающие дифракцию Фраунгофера плоской волны на бинарной знаковой функции от многочлена Эрмита (ур. (28), (31)); показано, что знаковые бинарные функции можно рассматривать как пропускание оптических элементов, формирующих эрми-товые пучки с эффективностью не меньше 80% (ур. (34), (37) и Таблица).

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

Эта работа выполнена при поддержке Российского фонда фундаментальных исследований

(№ 96-01-00021).

Литература

1. М.Б. Виноградова, О.В. Руденко, А. Л. Сухо-руков Теория волн, М., Наука, 1979.

2. А. Ярив. Введение в оптическую электронику, М., Высшая школа, 1983.

3. V.A. Soifer, M.A. Golub. Laser beam mode selection by computer-generated holograms, CRC Press, Boca Raton, 1994.

4. М.А. Голуб, В.С. Павельев, В.А. Сойфер. Построение итеративного алгоритма расчета фазовых дифракционных элементов, формирующих заданное одномодовое распределение на основе применения метода обобщенных проекций, Компьютерная оптика, М., МЦИТИ, 1995, вып. 14-15, часть 2, с. 85-93.

5. V.A. Soifer, V.V. Kotlyar, L.L. Doskolovich. Iterative methods for diffractive optical elements computation, Taylor and Francis, London, 1997.

6. J.P. Kirk, A.L. Jones. Phase-only complex-valued spatial filter, J. Opt. Soc. Am., 1971, v. 61, no. 8, p. 1023.

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

8. В.И. Смирнов. Курс высшей математики, М., Наука, 1971, т. IV, ч. 1

9. A. Fedotovsky, K. Lehovec. Optical filter design for annular imaging, Appl. Opt., 1974, v. 13, № 12, p. 2919-2923.

10. J. Mait. Design of binary-phase and multiphase Fourier grating for array generation, J. Opt. Soc. Am. A, 1990, v. 7, no. 8, p. 1514-1528.