УДК 535.42
ДИФРАКЦИЯ ПРОСТРАНСТВЕННО-ОГРАНИЧЕННОГО ПУЧКА НА РАДИАЛЬНО-СИММЕТРИЧНЫХ ДИФРАКЦИОННЫХ ОПТИЧЕСКИХ ЭЛЕМЕНТАХ
© 2008 С.И. Харитонов1, Н.Л. Казанский1,2, А.Ю. Дмитриев1
Институт систем обработки изображения РАН Самарский государственный аэрокосмический университет
В работе предложен метод решения задачи дифракции на оптических элементах с радиально-симметричным, не зависящим от продольной координаты распределением диэлектрической проницаемости. Предлагаемый метод является обобщением метода связанных волн (RCWA) для радиально-симметричных структур. В качестве базиса для разложения решения были выбраны конические волны, представляющие собой решения системы уравнений Максвелла в среде с постоянной диэлектрической проницаемостью.
Фокусатор, дифракционный оптический элемент, уравнения Максвелла, дифракция
Введение
Наряду со сферической линзой, фокусирующей свет в точку, широко применяются на практике оптические элементы, фокусирующие свет в кольцо.
Имеется ряд работ, посвященных исследованию фокусировки в кольцо в рамках геометрической оптики [1-7]. Как известно, приближение геометрической оптики не позволяет оценить ширину кольца, энергетическую эффективность фокусировки. В статье [8] исследуется фокусировка в кольцо с помощью оптического элемента с бинарной
фазовой структурой вида sign ^ J0 (krr0 / f0 )j,
который является дополнением к линзе. В работах [9,10] предложены несколько видов функций комплексного пропускания, описывающих оптический элемент с кольцевым импульсным откликом, и получены интегральные представления для интенсивности светового поля в фокальной плоскости вблизи кольца. В статье [11] было получено распределение интенсивности в фокальной плоскости (вблизи кольца) пары аксикон-линза. Дифракционный оптический элемент (ДОЭ) с повышенной глубиной фокуса или фокусатор в отрезок на оптической оси был впервые предложен в работе [12]. В дальнейшем было предложено несколько вариантов расчета, совершенствования и исследо-
вания фазовых функций таких ДОЭ [4, 1322]. Использование оптических элементов с повышенной глубиной фокуса актуально для лазерных проигрывателей компакт-дисков [23], для получения оптического разряда в газе [24], для создания опорной световой линии в метрологии [25], лазерных технологических установках [26].
Во всех изложенных работах прохождение поля через оптический элемент рассчитывался в приближении геометрической оптики. Для расчета поля после оптического элемента использовался интеграл Кирхгофа. Цель данной работы - разработать метод расчета поля внутри ДОЭ с радиальносимметричным распределением диэлектрической проницаемости. Предложенный метод является аналогом метода RCWA (Rigorous Coupled Wave Analysis), используемого для расчета дифракции полей на периодических структурах [27-29]. В работе [30] рассмотрено использование компактной записи для решения системы уравнений Максвелла с помощью многомерных матриц. В работе [30] рассмотрен наиболее общий способ решения системы уравнений Максвелла с помощью решения системы интег-ро-дифференциальных уравнений. В настоящей работе рассматривается применение данного метода к решению уравнений Максвелла в цилиндрической системе координат. Задача сводится к решению системы
обыкновенных дифференциальных уравнений с постоянными коэффициентами. В качестве базиса для описания результирующего поля были выбраны конические волны с определенным орбитальным угловым моментом, представляющие собой решения системы уравнений Максвелла в среде с постоянной диэлектрической проницаемостью. В скалярном приближении конические волны с определенным угловым моментом подробно рассмотрены в работах [32-34]. Эти пучки использовались авторами для вращения микрочастиц [35-36].
1. Решение эволюционных уравнений в операторной форме
Напомним основные элементы теории представлений, изложенной в работе [30]. Пусть имеется операторное уравнение, записанное в эволюционной форме:
і- =Щш),
k Эх3 1 1
Н =
0
В
А
0
где Н! - матричные элементы оператора Н в Б-представлении.
В результате мы получили систему ин-тегро-дифференциальных уравнений. Следует отметить, что вследствие специфической структуры оператора Н можно выбрать базис 1%^), в котором будут отличны
от нуля только следующие матричные элементы:
Н „1=Г “'|А|К„3), Н Ю4 = % Ю'|АЮ% 4), н Ю2 =% “А їй, н“ =% "2іа; %,), н Ю'3 =% Ю31 в|%„^, н 2 =% "3ів% л,
н “ =% ю4ів
(1) ню4 =%-|В%
(5)
(6)
(7)
(8) (9)
(10)
(11)
(12)
где Н - матричный дифференциальный оператор размерностью 4х4. В данной работе использованы обозначения, введенные П. Дираком в работе [31].
Пусть матричный дифференциальный оператор имеет вид
В данном случае используются матрицы 2 х 2, а в качестве базисных векторов используются столбцы из двух элементов, содержащие ненулевые компоненты векторов с размерностью, равной 4.
Система интегро-дифференциальных уравнений распадается на две части:
где А и В - матричные дифференциальные операторы размерности 2x2.
Представим решение уравнения в виде разложения по базису |РШ):
k ЭХ k=3,4
для 5 = 1, 2 ,
і Э!
k Эх3
= I Г Н! (Х3) V“ (х3) ёа
(13)
= ЦГ Н! (х3) V * (Х3) ёа
(14)
k=1,2
(3)
Vю5 | х3 | решением
Будем называть уравнения (1) в Б-представлении.
Уравнение для функции Vш ( х3 j имеет
вид
^ = 21. ^ IV (х3 )аа, (4)
для 5 = 3,4.
Можно также получить систему интег-ро-дифференциальных уравнений второго порядка:
1 д^
k2 Э(х3 )2
ЦТ рт (хз )1 vpm (хз) ар,
т=1,2
-р!
рт
ЦІН!1 (х5) Нрт (х3) ёа
По аналогии
1 Э
^ Э( х3 у
I Гот (хз )1 vрm (хз) ар,
а
юk
ЦТ Н^ (х3) Нрт (х3) ёа
(16)
(17)
(18)
2. Уравнения Максвелла в цилиндрических координатах
В данном разделе используем физические компоненты векторов. Для обозначения цилиндрических координат используем обозначения (х1, х2, х3 ^, где х1 - радиальная ко-
2 „3
ордината, х - полярный угол, х - расстояние вдоль оптической оси.
В цилиндрической системе координат систему уравнений Максвелла для комплексных амплитуд можно представить в виде
і ЭЮ = Н\ш), k Эх
(19)
Н
0
— А k Эх1
і 1 Э k х1 Эх
2 2
-і Э
А + ikє
k Эх —і 1 Э k х1 Эх
ц
2 1
і Э
єk Эх1 2 і 1 Э
-Ц —/^є
1 1
єk х1 Эх2 2 0
0
±_ _Э єk Эх
D,+ik
11
Ц2 — ik —1Д-Ц
єk х1 Эх2 1 0
0
(20)
|Ш>:
Н1
Н
(21)
где Ц1, Ц2 - дифференциальные операторы следующего вида:
1 Э| х1
Ц
ц2
1
Эх1
Э
х1 Эх2
(22)
(23)
Оператор Ц1 действует на функцию следующим образом:
А/ =
1 Э(х'/
Эх1
(24)
Здесь х1, х2, х3 - цилиндрические координаты, связанные с обычными декартовыми координатами (х, у, г) следующим преобразованием:
: х1 СОБI х2
у = х1 БІП I" х2
(25)
(26)
г = х3. (27)
Е1, Е2, Е3, Н1, Н2, Н3 - компоненты векторов электрического и магнитного полей в цилиндрической системе координат; є -диэлектрическая проницаемость среды;
k = : - волновое число в вакууме; : - круговая частота; с - скорость света.
2.1. Распространение света в среде с постоянной диэлектрической проницаемостью
В качестве простого примера рассмотрим распространение волны в среде с постоянной диэлектрической проницаемостью. Представим решение системы уравнений Максвелла в виде
E fx1, x2, x3 j = F (x1, x2)exp fikgx3 j, (2В)
H I" x1, x2, x3 1 = G f x1, x2 0 exp(ikgx3 1. (29)
ЭF
ik tF - ikG2 =^r,
* 1 2 ax1
-ik eF + ik tG2 =
x Эx
*=д-
G2
ikTF2 + ikGl = xl Эx2 ,
ik eF2 + ik tG = ^% Эx
Решения С1 и Р2 выражаются через определители Крамера следующим образом:
F2 = АЦ 2А
G = А4
G-=T'
(36)
(37)
Подставляя полученные выражения для £1, Е2, Н1, Н2 в третью пару уравнений Максвелла, получаем уравнения для компо-
нент F3, G3:
x Эx
+
Эx1
+
Такое представление справедливо в случае, если свойства среды не зависят от координаты х3; у- константа, возникающая при разделении переменных.
Для определения р и 02 используется следующая пара уравнений Максвелла:
Эx2
V
1 Э '
1 Э2 F3
-11 Э| x2
(3В)
+ k 2p2 F3 = G,
x1 Эx1
+-
ЭG3
Эx1
+
Эx2
1 Э 2G3
x112 ЭСx212
+ k 2p2G3 = G,
(3G) p2 =(e-T
(31)
(39)
(4G)
Решения выражаются через определители Крамера А1, А, А2 следующим образом:
где Ь - конический параметр.
Для решения полученных уравнений для функций G3, Р3 разделяем переменные следующим образом:
(32)
(33)
F3 С x1, x21 = Z ( x1 I exp f imx
G3 (x1, x2
Y (x1 ] exp (imx
(41)
(42)
Выражения для определителей приведены в приложении.
По аналогии для определения G1 и Р2 используется следующая система уравнений:
1 дР3
где т - целое число, определяющее орбитальный угловой момент количества движения поля, I ^ х1) и У ^ х10 - цилиндрические
функции. Подставляем это представление в уравнения для функций G3 и Р3 и получаем
уравнение для функций IС х10 и У(х1 ]. При-
(34)
(35)
веденные решения, а также линейные комбинации с одинаковыми индексами т будем называть решениями с определенным угловым моментом:
1 Э
x Эx
д ЭZ С Эx1
m
Z + k 2p2 Z = G.
(43)
1 Э
x1 Эx1
і Э7
' Эx1
m
Y + k p2Y = G.
(44)
Функции Z С x1 ] и YI x1 ], как и следо
+
Эx1 k2 p2 (x
Эx1
v /
д ЭZ + x —— +
Э^
і / \2 I
2q2 L„1 Г M„2
m
Z = G.
Эz
ґд2_л
Решением этого уравнения будут цилиндрические функции.
Запишем теперь решение для Е3 и Н3:
Е3 = Zm [k^x11 exp [ikgx3 lx
ik
F = —
ik
g2 = — 2А
ik
F2 = — 2А
T
Эр 1 ЭG3
-+
Э^ x1 Эx2
tlG ,рэр
x1 Эx2 Эx1 T Эр ЭG3
G1
ik
i
x1 Эx2 Эx1
ЭG3 e ЭF3
T
Эx1 x1 Эx2
1 Э[xG2 ) 1 Э^
1
G3
-i
x‘ Эx1 x1 Эx2
V J
/ -4 і 1F 1 \
2
1 Э( x'F2 I _Д Э<Д_
x1 Эx2
1
Эx1
получаем выражение для конических волн с определенным орбитальным угловым моментом:
F =-
вало ожидать, удовлетворяют одному и тому же уравнению.
Перепишем уравнение для функции
2 ^ х10 в виде
2 Э (Э2Л
. Эр mG3
/g
(45)
G2
F2
G,
ik
migG3 ЭF3
ч x1 Эx1 у
ik і migF3 ЭG3 Л
А ( x1 Эx1
k і
j't-
ЭG3 + meF3
Эx
(55)
(56)
(57) (5В)
Сделаем замену переменных kpx1 = z и получим уравнение Бесселя
1. Случай поперечной магнитной волны ( G3 = G ):
+ z у- + [ z2 - m21Z = G. (46)
F1
G2
ik і ЭрЛ
т 3
V
ik
V
дк1
Эр Л Эx1
, 21 , 2 3 (47) xexp(imx = Q(x ,x ,x ),
H3 = Zm fkpx11 exp (ikgx3 1exp fimx21. (4В)
ik
g-=d
ik і T Эр Л
x1 Эx2 /
eG Э^3 x1 Эx2
Для дальнейших рассуждений будем использовать выражения для компонент р,
р, ^, G2 через Р3, G3:
(59)
(6G)
(61)
(62)
(63)
А = k 2e- k2 т2.
2. Случай поперечной электрической волны ( F3 = G ):
(49)
(5G)
(51)
(52)
(53)
(54)
F1
G2
ik
3G3 x1 Эx2
x1 Эx2 ЭG3
ik
F = —
/
\
ЭG3
Эx1
(64)
(65)
(66) (67)
Подставляем в данные выражения решения с разделяющимися переменными и
2.2. Матричное представление решений с определенным значением орбитального углового момента
Выражение для компонент конических волн с определенным орбитальным угловым моментом можно переписать в матричном виде:
А
О,
02
ехр (ік ух31 = —
т
У
_дУ_ Эх1 ± ЭУ
±уЭ7
±у т у
х
ехр І ік ух3 I ехр
ішх2
(68)
р2
01
02
ехр І ік ух3
ік
1
±У
Э1_
Эх1
, іт
±у— 2
х
ітєп _
Э7
Эх1
ехр І ік ух3 І ехр І" ітх
(69)
Следует отметить, что конические волны ортогональны между собой при различных конических параметрах, ортогональны между собой при различных индексах, определяющих угловой момент, и также ортогональны для волн различных типов (имеются в виду волны поперечно-электрические и поперечно-магнитные). Конические волны
можно выбрать в качестве базиса. Однако данный базис не очень удобный, так как не является ортогональным. Это приводит к тому, что выражения для матричных элементов становятся достаточно сложными. В качестве базиса можно также выбрать систему ортогональных бивекторов
ік
1
'у|т Л (кЬх'
у т(крх-
х 1
т
х
Э
і Jm [крх1
------J (кВх1
Эх1 т (
О
О
О
О
_е0 х Jm (кРх' ) у^т [ кРх'
Є0 Эх1 ^ ^ ) уіт'1т(^
(70)
где 3т ^к$х10 - функция Бесселя.
Приведенные базисные векторы являются линейной комбинацией конических волн, распространяющихся в различных направлениях. ||ж|| - диагональная матрица, обеспечивающая нормировку базисных
функций. Под нормой вектора будем понимать скалярное произведение П=П П.
Под скалярным произведением будем понимать выражение
<П||П) = 1 П+( х10П (х11 х1^1, (71)
где VІ х10 - вектор-столбец, V х11 - комплексно сопряженная вектор-строка.
3. Распространение волн с определенным угловым орбитальным моментом в радиально-симметричной среде
Для дальнейших выкладок представим оператор Г амильтона-Максвелла в виде
і О А1
В О
н
(72)
где
■-А1 А
к Эх е
і 1 Э
ек х1 Эх
А _ ік
.2 2
1-Э-1 Б+ік к Эх1 е
—-—— Б ек х1 Эх2 1
(73)
В =
- -ЭТ б2 к Эх
і 1 Э „ ..
кх Эх-+іке
_і Э
Д _іке
к Эх1
—-—— Б к х1 Эх2 1
(74)
Следует отметить, что в пространстве состояний с определенным орбитальным
Э
моментом действие оператора
Эх2
іт , где
т - число, определяющее орбитальный момент.
В качестве базиса будем использовать (7О). Для удобства приведем матрицу из ненулевых компонент базисных векторов
(V V V V )
\ а1 а 2 а3 а 4/
ік
у~г Jm {ках11 — J [ках11 _еО — J (ках11 у—J (ках1
Эх т 1 0 х1 т ( 0 О х1 т І і Эх1 т І
у—Jm |І к ах11 _—Jm ( к ах1 ] еО —, Jm {к ах11 imуJm (к ах1
“О Эх1 т
(75)
На практике для упрощения вычислительной процедуры вместо непрерывного базиса необходимо использовать базис дискретный. Это будет соответствовать замене интегральных выражений интегральными суммами. В свою очередь это означает, что базисные вектора имеют вид
V™ = £ 0(а-а„„). (76)
п
Здесь 8(а-апт) - функция Дирака. Величины апт можно выбрать следующим образом:
Л (к а п,Д )= 0, (77)
Л, (к ап 2 Д ) = 0, (78)
5ХТ Л, (к а пз Д ) = 0, (79)
Л Л (к а „4 Д ) = 0. (80)
ах
С физической точки зрения такой выбор базиса означает, что оптическая система помещена в цилиндр радиуса Д, обладающий абсолютной проводимостью. Наличие абсолютно проводящего цилиндра вносит искажения в первоначальную постановку задачи дифракции. Однако в каждом конкретном случае можно выбрать значение Д таким образом, чтобы это влияние было незначительным. Например, в задаче фокусировки гауссова пучка в кольцо необходимо выбрать Д таким образом, чтобы он был значительно больше поперечного размера гауссова пучка. В этом случае влияние цилиндра на поле в фокальной плоскости будет минимально. При дальнейшем распространении пучка будут наблюдаться многократные отражения
от цилиндрическом поверхности, что приведет к значительному отличию поля. Однако, если увеличить Я, можно добиться, чтобы при данном расстоянии от ДОЭ искажение поля, вызванное наличием проводящего цилиндра, были минимальны.
При таком выборе счетного базиса все интегралы заменяются интегральными суммами, а интегро-дифференциальные уравнения превращаются в систему обыкновенных дифференциальных уравнений.
Система интегро-дифференциальных
уравнений распадается на две части:
і
к дх3
•=11 «п (х5) V* (х5),
п к=1,2
(82)
і
к дх3
II 47 (х5) V* (х3),
(81)
к=3,4
для 5 = 1, 2 ,
для ^ = 3,4.
Полученную систему необходимо свести к системе дифференциальных уравнений второго порядка.
4. Дифракция на ДОЭ с радиальносимметричным распределением диэлектрической проницаемости
Рассмотрим решение задачи дифракции на дифракционном оптическом элементе с радиально-симметричным распределением диэлектрической проницаемости. Пусть входной пучок падает из области 1 (рис. 1) на дифракционный оптический элемент.
Х2
XI
— /
— 1 ( ^ хз
—— '
——
— 1
— '
р
Рис. 1. Оптическая схема фокусировки с помощью радиально-симметричного ДОЭ
Распределение комплексной амплитуды входного пучка имеет вид
■■II г ю ехр і
п 5=1,2
(83)
где е1 - диэлектрическая проницаемость среды 1, 1п определяют вклад различных конических волн во входном поле, |£ж) - вектор-столбец из четырех элементов. Столбцы с различным вторым индексом запишем в виде матрицы
^ ^ (крх‘) 7 ^ (^'1 *п (крх‘) 7 л (крх')
^7 ^ (кЬ х'] -^ 'п (крх‘) -Т7 ^ -дх ^( ^
-£1 х( кРх' ] ^ ( кРх‘ ) ( кРх' ] -У^ ( ( (
^ ^ ^ ] т7 (кРх') е ^( ч*) -У1т зп (крх1 х V
Поле в области 1 представляется в виде суммы падающей и отраженной волн и имеет вид
Щ =11 /"Ю ехр (ік^/е^х51+
п 5=1,2
+1 I Кт\Рш) еХР \-1Ч Є -ашх3 1,
п 5=3,4
Щз)
X
п 5=1,2
хехр(іку]е1 -аю (х3 -ё)),
1 д2уп (х3)
- к2 д( х3 )2
где і = 1,2,
М7 = I А^Н’:.
д=3,4
--мпу
п1^,т
уУ (х3) = I Е7п (а+тпех’ (ікт ппх3
пп
+ а-ппех’(-іктпп (х3 -ё))),
МЩЕ’1 = Ц2 Ещ
’I тп г^тп тп '
+
.’і
Остальные V* можно найти по формуле Э (91)
і I о
7, 1
’1 _________
т5 дх3
(85)
где Яп5 определяют вклад различных конических волн в сумму, описывающую отраженное поле.
Поле в области 3 за оптическим элементом представляется в виде
(86)
к „
I = 3,4.
Многомерная матрица О^ удовлетворяет соотношению
I О’1Аш: = 8* (92)
т,5=1,2
Систему дифференциальных уравнений первого порядка можно также свести к системе уравнений второго порядка следующего вида:
1 дVі (х3)
т1^,т
2 К*
где ё - толщина дифракционного оптиче-
ГП5
определяют вклад различных конических волн в сумму, описывающую прошедшее поле.
Рассмотрим распространение света внутри ДОЭ с радиально-симметричным, не зависящим от х3 распределением диэлектрической проницаемости. Решение описывается с помощью системы дифференциальных уравнений первого порядка. Полученную систему уравнений можно свести к системе уравнений второго порядка. При этом сокращается размерность системы уравнений
к д( х3)
где і = 3,4,
Ыт1 = I Нтілп
т рд т '
(93)
(94)
д=3,4
Решение системы дифференциальных уравнений второго порядка имеет вид
V’1 ( х3 )^ Р7і (4+тпехр ( ікЦті^х‘
тп
+ Ъ~тпех’ ( ікЦтп (х3 - ё))
ту5Чр’1 = Ц 2 рщ ’I тп г^тп тп '
+
(95)
(96)
(87)
Остальные Vр1 в этом случае можно найти по формуле
- I К
д*
рї
т дх3 ’
(97)
(88)
Решение системы дифференциальных уравнений второго порядка имеет вид
Л т,5=3,4
где I = 1,2.
Многомерная матрица От5 удовлетворяет соотношению
I кт5А7 = 8„1 . (98)
(89)
(90)
т, 5=3,4
Поле внутри оптического элемента представляется в виде
(99)
п 5=1,2,3,4
где | ¥т) - вектор-столбец из четырех элементов (70). Набор коэффициентов а+тп, Ъ+тп, Ятп, Ттп находится из условия непрерывности полей на границах оптического элемента. Далее, используя полученные коэффициенты, получим выражение для компонент электромагнитного поля на выходе оптического элемента. Для того чтобы получить поле за оптическим элементом, необходимо использовать теорию Кирхгофа-Котлера.
5. Заключение
В работе разработан метод решения задач дифракции полей с определенным значением углового момента. Предложенный метод является аналогом метода RCWA, используемого для расчета дифракции полей на периодических структурах. Получен вид системы обыкновенных дифференциальных уравнений для нахождения решения системы уравнений Максвелла в V -представлении в пространстве функций с определенным угловым моментом
Приложение А
Используя правила Крамера, решаем систему линейных уравнений относительно Е и Н2. Определители Крамера имеют вид
2 У к - е 2 к N А (100)
дК ік дG3 А, = ік у—3 + 3, 1 1 дх1 х1 дх2 (101)
ік удG3 дД А 2 = х дх2 +ік ^ (102)
Аз = ^ - ік *?, 3 х1 дх2 дх1 (103)
А = ік уд°з ік едрз 4 у дх1 х1 дх2 ' (104)
Приложение В
Для того чтобы вычислить матричные символы необходимо найти действие следующих операторов на вектора
л\уа з) = -к-1 “3/ А(а)
і д 1 _ ( іт Т (, 1 0 (і д 1 _
е^-Т Зт (каХ + - ТУ" Д1 + ік
х ( 0 к дх е
к дх1 е 2 ґ і 1 д
ек х дх
-А
Д - ік
л л
е0 —- Зт ( к ах1
0 дх1 т'
ек х дх
е0 —- З т(к ах1
0 дх1 т (
(105)
Л\¥.
ік
А(а)
- - -дТ 1Д2 к дх е
^ і 1 д
У—1 Зт I к ах1 дх т (
+
і д 1 „
------—Д + ік
к дх е
Л( іту З„ (к ах1 '
(
ек х дх
,1 П„2 2
Д2 - ік
У—1 Зт \ к ах1 дх т (
ек х дх
Д
(іту
Зт | к ах1
(106)
( • "\ Л( "\
і д д д
ВК
а1
ік
А(а)
к дх1 2
ґ і 1 д
У—1 Зт I к ах1 дх т (
+
-і д
к дх
Д - ік е
1-^1
іт
У~Зт !ках
(
к х дх
Д, + ік е
(
У—1 Зт I ках1 дх т (
А
+
-і 1 д
,1 П„2 1
к х дх
Д
(107)
і д
ВК 2) = ~іГТ 1 1 А (а)
к дх1 2
/ і 1 д
“ т Зт і к ах'
х
к х дх
Д Д2 + ік е
+
-і д
к дх
Д - ік е
л л
у іт ( 1 л! і-і 1 д
хГ Зт іках 1 +
(
к х дх
1 3,.2 Д
—1 Зт і к ах1
дх1 т'
аіх
1 т
Напомним выражение для производных цилиндрических функций:
Jm ( k aX'
m
k ax1
Jm (k ax11 - Jm (kax1
(109)
Действие операторов на состояния с определенным значением углового момента определяется следующим образом:
_Э_
Эх1
(1 д (11
- DH 2
1 2 - =3x1 V-0
DiH2 ) +
_д_
дх1
( 1
- d2
-
д (1Л
дх1
e
ГЛ 1 д ТЛ
D +---- D2 e дх
= im
= im
1A -
x1 дх1 -1 д
1 д 1
+ -
x1 дx1
(11 1
— + -
- -
V У V
- дx1 x1 1 д
x1 дx1
//
-i-D=±
ax' ' 3x‘
1 д^ x1 E2
x1 дx1
1 д
x1 дx1 1 д
x дx
D
4
3x'
im
1
x д( x1
112 дx1
im
V
3x1
д ^ д 1 -im im д
—tD = im—;—т =--------- + -
3x ‘
3x' x1
2 (x1
3x1
1 д d = -m2
1 -Ч..2 2
x 3x
(110)
(111)
(112)
(113)
(114)
(115)
Благодарности
Работа выполнена при поддержке фонда «Фундаментальные исследования и высшее
образование» (RUXO-014-SA-06) и грантов РФФИ № 07-07-97601-р_офи, 07-01-96602-р_Поволжье_а, 07-07-91580-АСП_а, 07-0700210, 08-07-99005-р_офи и Фонда содействия отечественной науке.
Библиографический список
1. Belanger, P.A. Ring pattern of a lens-axicon doublet illuminated by a Gaussian beam [текст] / P.A. Belanger, M. Rioux // Applied Optics. - 1978. - Vol. 17, № 7. - P.1080-1086.
2. Belanger, P.A. Diffraction ring pattern at the focal plane of a spherical lens-axicon doublet [текст] / P.A. Belanger, M. Rioux // Journ. Canadien de Physique. - 1976. - Vol. 54.
- P.1774-1780.
3. Farn, M.W. Effect of VLSI fabrication errors on kinoform efficiency [текст] / M.W. Farn, J.W. Goodman // Proceedings SPIE.
- 1990. - Vol. 1211. - P.1256-136.
4. Голуб, М.А. Вычислительный эксперимент с элементами плоской оптики [текст] / М. А. Голуб [и др.] // Автометрия. -1988. - № 1. - С. 70-82.
5. Голуб, М.А. Дифракционный расчет оптического элемента, фокусирующего в кольцо [текст] / М. А. Голуб [и др.] // Автометрия - 1987. - № 6. - С. 8-15.
6. Казанский, Н.Л. Исследование дифракционных характеристик фокусатора в кольцо методом вычислительного эксперимента [текст] / Казанский Н.Л. // Компьютерная оптика. - М.: МЦНТИ, 1992. - Вып.10-11. - С.128-144.
7. Doskolovich, L.L. Focusators into a ring [текст] / L.L. Doskolovich [and others] // Optical and Quantum Electronics. - 1993. -Vol. 25. - P.801-814.
8. Fedotowsky, A. Optimal filter design for annular imaging [текст] / A. Fedotowsky, K. Lehovec // Applied Optics. - 1974. - Vol. 13, № 12. - P.2919-2923.
9. Коронкевич, В.П. Киноформные оптические элементы: методы расчета, технология изготовления, практическое применение [текст] / В.П. Коронкевич [и др.] // Автометрия. - 1985. - № 1. - С.4-25.
1G. Коронкевич, В.П. Киноформные оптические элементы с кольцевым импульсным откликом [текст] / В.П. Коронкевич [и др.] // Препринт № 265 ИАиЭ СО АН СССР.
- Новосибирск, 19S5. - 23 с.
11. Belanger, P.A. Diffraction ring pattern at the focal plane of a spherical lens - axi-con doublet [текст] / P.A.Belanger, M. Rioux // Journ. Canadien de Physique. - 1976. - Vol. 54.
- P.1774-17SG.
12. Голуб, М.А. Фокусировка когерентного излучения в заданную область пространства с помощью синтезированных на ЭВМ голограмм [текст] / М. А. Голуб [и др.] // Письма в ЖТФ. - i9Si. - Т. 7, вып. 1G. -C61S-623.
13. Казанский, Н.Л. Анализ характеристик фокусаторов лазерного излучения методом вычислительного эксперимента [текст] / Н.Л. Казанский // Диссертация на соискание ученой степени кандидата технических наук. - Куйбышев: КуАИ, 19SS. - 1S3 с.
14. Казанский, Н.Л. Процедура корректировки фазовой функции фокусатора по результатам вычислительного эксперимента [текст] / Н.Л. Казанский // Компьютерная оптика. - М.: МЦНТИ, 19S7. - Вып.1. - C9G-96.
15. Пальчикова, И.Г. Киноформные оптические элементы с увеличенной глубиной фокуса [текст] / И.Г. Пальчикова // Компьютерная оптика. - М.: МЦНТИ, 19S9. -Вып. 6. - С. 9-І9.
16. Васин, А.Г. Расчет и исследование когерентного волнового поля в фокальной области радиально-симметричных оптических элементов [текст] / А.Г. Васин [и др.] // Препринт № 3G4 ФИАН СССР. - М.: ФИАН, 19S3. - 3S с.
17. Сойфер, В.А. К расчету фокусато-ра в соосный отрезок [текст] / В. А. Сойфер // Оптическая запись и обработка информации.
- Куйбышев: КуАИ, 19SS. - С. 45-52.
1S. Doskolovich, L.L. Analysis of qua-siperiodic and geometric optical solutions of the problem of focusing into an axial segment [текст] / L.L. Doskolovich [and other] // Optik.
- 1995. - Vol. 1G1, № 2. - P.37-41.
19. Kazanskiy, N.L. Correction of fo-cuser phase function by computer-experimental method [текст] / N.L. Kazanskiy // Computer Optics. - 19S9. - Vol. 1, № 1. - P.69-73.
20. Khonina, S.N. Calculation of the fo-cusators into a longitudinal line-segment and study of a focal area [текст] / S.N. Khonina, V.V. Kotlyar, V.A. Soifer // Journal of Modern Optics. - 1993. - Vol. 40. - P. 761-769.
21. Michaltsova, I.A. Kinoform axicon [текст] / I.A. Michaltsova, V.I. Nalivaiko, I.S.Sol-datenkov // Optik. - 1984. - Vol. 67, № 3. -P. 267-270.
22. Kolodziejczyk, A. The light sword optical element - a new diffraction structure with extended depth of focus [текст] / A. Kolodziejczyk [and other] // Journal of Modern Optics. - 1990. - Vol.37, № 8. - P.1283-1286.
23. Brenden, B.B. Optical playback apparatus focusing system for producing a prescribed energy distribution along an axial focal zone [текст] / B.B. Brenden, J.T. Russel // Applied Optics. - 1984. - Vol. 23, № 19. -P. 3250-3253.
24. Tremblay, R. Laser plasmasoptically pumped by focusing with axicon a CO2-TEA laser beam in a high-pressure gas [текст] / R. Tremblay [and other] // Optics Communications. - 1979. - Vol. 28, № 2. - P. 193-196.
25. Воронцов, М.А. К расчету фокуса-
торов лазерного излучения в дифракционном приближении [текст] / М.А. Воронцов,
А.Н. Матвеев, В.П. Сивоконь // Компьютерная оптика. - М.: МЦНТИ, 1987. - Вып.1. -С.74-79.
26. Rioux, M. Linear, annular and radial focusing with axicons and applications to laser machining [текст] / M. Rioux, R. Tremblay, P.A. Belanger // Applied Optics. - 1978. -Vol.17, № 10. - P.1532-1536.
27. Moharam, M.G. Stable implementation of the rigorous coupled-wave analysis for surface-relief gratings: enhanced transmittance matrix approach [текст] / M.G. Moharam [and other] // J. Opt. Soc. Am. A. - 1995. - Vol. 12(5). - P. 1077-1086.
28. Moharam, M.G. Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings [текст] / M.G. Moharam [and other] // J. Opt. Soc. Am. A. - 1995. - Vol. 12(5). - P. 10681076.
29. Li, L. Use of Fourier series in the analysis of discontinuous periodic structures [текст] / L. Li // J. Opt. Soc. Am. A. - 1996. -Vol. 13(9) - P. 1870-1876.
3G. Казанский, Н.Л. Компактная запись решений системы уравнений Максвелла в пространственно-частотном представлении [текст] / Н.Л. Казанский, М. Л. Каляев, С.И. Харитонов // Антенны. - 2GG7. - № 1G. -C. 13-21.
31. Дирак, П.А.М. Принципы квантовой механики [текст] / П.А.М. Дирак - М.: Наука, 1979.
32. Khonina, S.N. An analysis of the angular momentum of a light field in terms of angular harmonics [текст] / S.N. Khonina [and other] // Journal of Modern optics. - 2GG1. -4S(1G). - P. 1543-1557.
33. Котляр, B.B. Измерение орбитального углового момента светового поля с помощью дифракционного оптического эле-
References
1. Belanger, P.A. Ring pattern of a lens-axicon doublet illuminated by a Gaussian beam / P.A. Belanger, M. Rioux // Applied Optics. - 197S. - Vol. 17, № 7. - P.1GSG-1GS6.
2. Belanger, P.A. Diffraction ring pattern at the focal plane of a spherical lens-axicon doublet / P.A. Belanger, M. Rioux // Journ. Ca-nadien de Physique. - 1976. - Vol. 54. - P.1774-17SG.
3. Farn, M.W. Effect of VLSI fabrication errors on kinoform efficiency / M.W. Farn, J.W. Goodman // Proceedings SPIE. - 199G. -Vol. 1211. - P.1256-136.
4. Golub, M.A. Computing experiment with elements of flat optics / M.A. Golub [and other] // Avtometria. - 19SS. - N 1. - P. 7G-S2.
- [in Russian].
5. Golub, M.A. Diffractial calculation of the optical element focusing in a ring / M.A. Golub [and other] // Avtometria. - 19S7. -N 6. - P. S-15. - [in Russian].
6. Kazanskiy, N.L. Research diffractial characteristics of focusator in a ring using a method of computing experiment / N.L. Kazanskiy // Computer Optics. - 1992. - N. 1G-11. -P.12S-144. - [in Russian].
7. Doskolovich, L.L. Focusators into a ring / L.L. Doskolovich [and others] // Optical and Quantum Electronics. - 1993. - Vol. 25. -P.SG1-S14.
мента [текст] / В.В. Котляр [и др.] // Автометрия. - 2GG2. - S(3). - C. 33-44.
34. Khonina, S.N. Astigmatic Bessel laser beams [текст] / S.N.Khonina [and other] // Journal of Modern optics. - 2GG4. - 51(5). -P. 677-6S6.
35. Khonina, S.N. Rotation of microparticles with Bessel beams generated by diffractive elements [текст] / S.N.Khonina [and other] // Journal of Modern optics. - 2GG4. - 51(14). -P. 2167-21S4.
36. Skidanov, R.V. Micromanipulation in Higher-Order Bessel Beams [текст] / R.V. Skidanov [and other] // Optical Memory & Neural Networks (Information Optics), Allerton Press. - 2GG7. - 16(2). - P. 91-9S.
8. Fedotowsky, A. Optimal filter design for annular imaging / A. Fedotowsky, K. Lehovec // Applied Optics. - 1974. - Vol. 13, № 12. - P.2919-2923.
9. Koronkevich, V.P. Kinoform optical elements: calculation methods, manufacturing techniques, practical application / V.P. Koronkevich [and other] // Avtometria. - 1985.
- N 1. - P.4-25. - [in Russian].
10. Koronkevich, V.P. Kinoform optical elements with the ring pulse response / V.P. Koronkevich [and other] // Preprint 265 HA^ CO Academy of Science of USSR. - Novosibirsk, 1985. - 23 p. - [in Russian].
11. Belanger, P.A. Diffraction ring pattern at the focal plane of a spherical lens - axi-con doublet / P.A.Belanger, M. Rioux // Journ. Canadien de Physique. - 1976. - Vol. 54. -P.1774-1780.
12. Golub, M.A. Focusing of coherent radiation in the set area of space by means of the hologrammes synthesised on the computer / M.A. Golub [and other] // Letters for GTF. -1981. - V. 7, N 10. - P.618-623. - [in Russian].
13. Kazanskiy, N.L. The analysis of characteristics of laser radiation focusators by a method of computing experiment / N.L. Kazanskiy // The dissertation on competition of a scientific degree of a Cand. Tech. Sci. - Kujbyshev: KuAI, 1988. - 183 p. - [in Russian].
14. Kazanskiy, N.L. Procedure of updating of phase function of a focusator by results of
computing experiment / N.L. Kazanskiy //
Computer Optics. - 19S7. - N 1. - P. 9G-96. -[in Russian].
15. Paljchikova, I.G. Kinoform optical elements with the increased depth of focus / I.G. Paljchikova // Computer Optics. - 19S9. - N 6. -P. 9-І9. - [in Russian].
16. Vasin, A.G. Calculation and research of a coherent wave field in focal area of ra-dially-symmetric optical elements / A.G. Vasin [and other] // Preprint 3G4 FIAN USSR. -Moscow, FIAN, 19S3. - 3S p. - [in Russian].
17. Soifer, V.A. To focusator calculation in a coaxial piece / V.A. Soifer // Optical recording and information processing. - Kuibyshev, KuAI, 19SS. - P. 45-52. - [in Russian].
1S. Doskolovich, L.L. Analysis of qua-siperiodic and geometric optical solutions of the problem of focusing into an axial segment / L.L. Doskolovich [and other] // Optik. - 1995. -Vol. 1G1, № 2. - P.37-41.
19. Kazanskiy, N.L. Correction of fo-cuser phase function by computer-experimental method / N.L. Kazanskiy // Computer Optics. -19S9. - Vol. 1, № 1. - P.69-73.
2G. Khonina, S.N. Calculation of the fo-cusators into a longitudinal line-segment and study of a focal area / S.N. Khonina, V.V. Kotlyar, V.A. Soifer // Journal of Modern Optics. - 1993. - Vol. 4G. - P. 761-769.
21. Michaltsova, I.A. Kinoform axicon / I.A. Michaltsova, V.I. Nalivaiko, I.S.Soldatenkov // Optik. - 19S4. - Vol. 67, № 3. - P. 267-27G.
22. Kolodziejczyk, A. The light sword optical element - a new diffraction structure with extended depth of focus / A. Kolodziejczyk [and other] // Journal of Modern Optics. - 199G.
- Vol.37, № S. - P.12S3-12S6.
23. Brenden, B.B. Optical playback apparatus focusing system for producing a prescribed energy distribution along an axial focal zone / B.B. Brenden, J.T. Russel // Applied Optics. - 19S4. - Vol. 23, № 19. - P. 325G-3253.
24. Tremblay, R. Laser plasmasoptically pumped by focusing with axicon a CO2-TEA laser beam in a high-pressure gas / R. Tremblay [and other] // Optics Communications. - 1979. -Vol. 2S, № 2. - P. 193-196.
25. Vorontsov, M.A. To calculation of fo-cusators of laser radiation in diffracting approach / M.A. Vorontsov, A.N. Matveev, V.P. Sivok-onj // Computer Optics. - 19S7. - N 1. - P.74-79. - [in Russian].
26. Rioux, M. Linear, annular and radial focusing with axicons and applications to laser machining / M. Rioux, R. Tremblay, P.A. Belanger // Applied Optics. - 197S. - Vol.17, № 1G. - P.1532-1536.
27. Moharam, M.G. Stable implementation of the rigorous coupled-wave analysis for surface-relief gratings: enhanced transmittance matrix approach / M.G. Moharam [and other] // J. Opt. Soc. Am. A. - 1995. - Vol. 12(5). - P. 1G77-1GS6.
2S. Moharam, M.G. Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings / M.G. Moharam [and other] // J. Opt. Soc. Am. A.
- 1995. - Vol. 12(5). - P. 1G6S-1G76.
29. Li, L. Use of Fourier series in the analysis of discontinuous periodic structures / L. Li // J. Opt. Soc. Am. A. - 1996. - Vol. 13(9) -P. 1S7G-1S76.
3G. Kazanskiy, N.L. Compact record of decisions of system of Maxwell equations in spatially-frequency representation / N.L. Kazanskiy, M.L. Kalyaev, S.I. Kharitonov // Aerials. -2GG7. - N 1G. - P. 13-21. - [in Russian].
31. Dirak, P.A.M. Principles of quantum mechanics / P.A.M. Dirak - Moscow, Science, 1979. - [in Russian].
32. Khonina, S.N. An analysis of the angular momentum of a light field in terms of angular harmonics / S.N. Khonina [and other] // Journal of Modern optics. - 2GG1. - 4S(1G). -P. 1543-1557.
33. Kotlyar, V.V. Measurement of the orbital angular moment of a light field with the help an diffractive optical element / V.V. Kotlyar [and other] // Avtometria. - 2GG2. - S(3). -P. 33-44. - [in Russian].
34. Khonina, S.N. Astigmatic Bessel laser beams / S.N.Khonina [and other] // Journal of Modern optics. - 2GG4. - 51(5). - P. 677-6S6.
35. Khonina, S.N. Rotation of microparticles with Bessel beams generated by diffractive elements / S.N.Khonina [and other] // Journal of Modern optics. - 2GG4. - 51(14). -P. 2167-21S4.
36. Skidanov, R.V. Micromanipulation in Higher-Order Bessel Beams / R.V. Skidanov [and other] // Optical Memory & Neural Networks (Information Optics), Allerton Press. -2GG7. - 16(2). - P. 91-9S.
DIFFRACTION OF THE SPACE-LIMITED BEAM ON THE DIFFRACTIVE OPTICAL ELEMENTS WITH RADIAL SYMMETRY
© 2008 S.I. Kharitonov, N.L. Kazanskiy, A.Yu. Dmitriev
Image Processing Systems Institute of the RAS,
S. P. Korolyov Samara State Aerospace University
We discuss a method for solving the diffractive problem on the optical elements with the symmetrical distribution of the permittivity. The method under discussion is the generalization of the rigorous coupled wave analysis (RCWA) for structures with the radial symmetry. As the basis for the decomposition of the solution the conic waves are choosen. This conic waves are the solution of the Maxwell’s equations in the medium with the constant permittivity.
Focusator, diffractial optical element, Maxwell equations, diffraction
Сведения об авторах
Харитонов Сергей Иванович, Учреждение Российской академии наук Институт систем обработки изображений РАН, Самара, Россия, старший научный сотрудник (prognoz@smr.ru). Область научных интересов - физическая оптика, компьютерная оптика, теория оптимизации, оптика, фазовая проблема
Казанский Николай Львович, Учреждение Российской академии наук Институт систем обработки изображений РАН, Самара, Россия, заместитель директора; Самарский государственный аэрокосмический университет имени академика С.П. Королева, профессор, email: ka-zansky@smr.ru. Область научных интересов - лазерные информационные технологии, компьютерная оптика, методы оптимизации, распределенные вычислительные системы, математическое моделирование, обработка изображений.
Дмитриев Антон Юрьевич, Учреждение Российской академии наук Институт систем обработки изображений РАН, Самара, Россия, стажер-исследователь (tonydm@yandex.ru). Область научных интересов - оптика, дифракционная оптика, геометрическая оптика.
Kharitonov Sergey Ivanovich, Establishment of the Russian Academy of Sciences Image Processing Systems Institute of the Russian Academy of Sciences, Samara, Russia, the senior scientific employee (prognoz@smr.ru). Area of scientific interests: physical optics, computer optics, the optimisation theory, optics, a phase problem.
Kazanskiy Nikolay Ljvovich, Establishment of the Russian Academy of Sciences Image Processing Systems Institute of the Russian Academy of Sciences, Samara, Russia, the Deputy Director; S. P. Korolyov Samara State Aerospace University, the professor, email: kazansky@smr.ru. Area of scientific interests: laser information technologies, computer optics, the optimisation methods, the distributed computing systems, mathematical modelling, processing of images.
Dmitriev Anton Yurjevich, Establishment of the Russian Academy of Sciences Image Processing Systems Institute of the Russian Academy of Sciences, Samara, Russia, the trainee-researcher (tonydm@yandex.ru). Area of scientific interests - optics, diffraction optics, geometrical optics.