ИССЛЕДОВАНИЕ ФОРМИРОВАНИЯ ОСЕВОГО ОТРЕЗКА С ПОМОЩЬЮ ОПТИМИЗИРОВАННОГО БИНАРНОГО ДОЭ
Хонина С.Н. 12, Скиданов Р.В. 12, Качалов Д.Г. 2, Павельев В.С. 12, Моисеев О.Ю. 1
1 Учреждение Российской академии наук Институт систем обработки изображений РАН,
2 Самарский государственный аэрокосмический университет имени академика С.П. Королёва
Аннотация
Проведено численное и экспериментальное исследование действия оптимизированного бинарного ДОЭ, предназначенного для формирования продольного отрезка с равномерной интенсивностью, в условиях изменения параметров освещающего пучка.
Ключевые слова: бинарный аксикон, дифракционный оптический элемент (ДОЭ), продольный световой отрезок, лазерное излучение, фазовый и амплитудный ДОЭ.
Введение
Дифракционные оптические элементы, предназначенные для формирования заданных распределений интенсивности вдоль оптической оси (дифракционные линзы с продольно-протяженным фокусом, многофокусные элементы и др.), находят применение во многих областях техники. Например, при выравнивании и совмещении конструкций на больших расстояниях [1]; при формировании изображений протяжённых или движущихся объектов, что особенно актуально в медицине [2-4], а также для неразрушающего исследования материалов [5] и устройств [6], в метрологии для сканирующих [7] и интерференционных [8, 9] оптических систем.
Классическим оптическим элементом (как рефракционным, так и дифракционным) для генерации осевой световой нити является аксикон [10]. Однако распределение интенсивности вдоль оси, формируемое аксиконом, очень неравномерно - кроме существенных флуктуаций наблюдается линейный рост интенсивности при увеличении расстояния от оптического элемента.
Расчёт фазовой функции ДОЭ, фокусирующего лазерное излучение в продольный отрезок, в рамках геометрооптического подхода рассматривался в работах [11-14].
Также известны итерационные методы расчёта ДОЭ [15], основанные на алгоритме Герчберга-Сакстона [16], которые позволяют учитывать дифракционные эффекты, однако не обладают сходимостью (может быть показана лишь нерасходимость). Кроме того, необходимость применения обратного оператора на каждой итерации метода затрудняет получение квантованной фазы оптического элемента, в то время как значительная часть технологий (главным образом литографических) ориентирована как раз на получение квантованного микрорельефа [17]. Поэтому в работе [18] для синтеза бинарного рельефа был применён метод последовательных приближений с бинаризацией фазовой функции на каждом шаге. Полученный ДОЭ имел вид кольцевой решётки из почти двух сотен колец, причём на периферии период был близок к 10 мкм при радиусе элемента 4 мм.
В работах [19, 20] обсуждался расчёт бинарных ДОЭ, формирующих множество фокусов, а также
фокусирующих в заданные кривые в электромагнитном приближении. Однако при формировании осевых световых распределений, удовлетворяющих параксиальному приближению, применение скалярной теории дифракции является вполне корректным и ресурсосберегающим подходом.
В случае, когда шаг поперечной дискретизации должен быть не меньше определённой величины, и соответственно, количество колец в структуре ДОЭ желательно свести к минимуму, предпочтительнее применение стохастических алгоритмов прямого поиска, так как поиск функции высоты микрорельефа можно организовать на множестве квантованных решений. В работе [21] с помощью стохастического алгоритма прямого поиска был рассчитан оптимизированный бинарный ДОЭ, предназначенный для формирования продольного отрезка с равномерной интенсивностью из освещающего га-уссового пучка.
В данной работе проведено численное исследование действия данного оптического элемента в условиях изменения параметров освещающего пучка. Рассмотрено три типа освещающего пучка: ТЕМ00, ТЕМ10 и плоский пучок с равномерным распределением интенсивности, ограниченный круглой апертурой радиуса ДОЭ. Также проведено моделирование изменения характеристик ДОЭ в случае исполнения его в виде фазового или амплитудного оптического элемента в зависимости от длины волны освещающего пучка.
Для предварительной экспериментальной верификации численных исследований по маске, рассчитанной в [21], был изготовлен амплитудный ДОЭ. Анализ полученных экспериментальных результатов при освещении данного оптического элемента излучением от различных лазеров показал хорошее согласование с численными результатами.
1. Моделирование действия бинарного формирователя осевого отрезка в различных условиях В данном разделе с помощью численного моделирования исследуется действие бинарного формирователя осевого отрезка, рассчитанного в [21] в различных условиях: при освещении гауссовым и равномерным пучками, излучением с различными
длинами волн, при исполнении ДОЭ в виде фазового или амплитудного оптического элемента.
Изначально рассматриваемый ДОЭ радиусом Я = 3,6 мм и числом отсчётов по радиусу 128 был рассчитан для формирования осевого отрезка протяжённостью Л = 20 мм на расстоянии от х1 = 200 мм до х2 = 220 мм от ДОЭ при освещении гауссовым пучком с радиусом в перетяжке о = 2 мм и длиной волны 1 = 1,06 мкм.
На рис. 1 показано радиальное сечение фазы рассматриваемого ДОЭ и сравнение результатов численного моделирования с помощью преобразования Френеля при освещении его расчётным гауссовым и равномерным пучками. ф, рад 3
2
1
а) 0
3 г, мм
в)
0
0,04 0,08 0,12 0,16 г, мм
Рис. 1. Радиальное сечение фазы (а), осевое распределение интенсивности (б) и распределение интенсивности в плоскости z = 210 мм (в) при освещении рассматриваемого ДОЭ расчётным гауссовым (сплошная линия) и равномерным (точечная линия) пучками
Как видно на рис. 1б, на отрезке от х\ = 200 мм до х2 = 220 мм при использовании для освещения расчётного гауссового пучка формируется равномерное распределение осевой интенсивности (среднеквадратичное отклонение от прямой линии составляет 5%). При освещении ДОЭ плоским пучком с равномерным распределением интенсивности, ограниченным круглой диафрагмой радиуса Я = 3,6 мм, осевое распределение несколько изменяется, осевое распределение интенсивности на заданном отрезке становится менее равномерным (среднеквадратичное отклонение от прямой линии составляет 19%), но значение интенсивности увеличивается примерно в полтора раза.
В последнем случае для наглядности можно сравнить действие оптимизированного ДОЭ с бинарным аксиконом вида:
t(r) = exp [i arg (J0(ar)] . (1)
При использовании чисто фазовых аксиконов (1) формирование осевой световой нити происхо -дит не сразу, а на некотором, иногда довольно значительном расстоянии от плоскости ДОЭ [22]. Максимальное расстояние, до которого воспроизводится световая нить, пропорционально радиусу ДОЭ и обратно пропорционально параметру бес-селевого пучка a:
kR
z =— , (2)
max ’ v ■'
a
где k = 2p /1, 1 - длина волны излучения.
Используя выражение (2), можно предположить, что при a = 85 мм-1 в диапазоне z е [200 мм, 220 мм]
осевая нить фазовым ДОЭ уже будет сформирована, но ещё не начнёт разрушаться. Как видно по рис. 2, это действительно так, однако, чтобы аккуратно прописать фазовую функцию такого аксикона, нужно хотя бы 1000 отсчётов на радиус. При значительном уменьшении числа отсчётов будут возникать сильные осцилляции осевой интенсивности (рис. 2).
На рис. 3 показано сравнение равномерности осевого распределения интенсивности на более протяжённом отрезке z е [190 мм, 230 мм] для рассчитанного ДОЭ и бинарного аксикона с a = 85 мм-1. Основное преимущество оптимизированного стохастическим методом ДОЭ перед бинарными акси-конами, частота периода которых пропорциональна параметру a, состоит в возможности формирования осевого отрезка близко от плоскости ДОЭ при достаточно крупном шаге дискретизации микрорельефа оптического элемента, что важно, в частности, при реализации фазовых распределений на динамических транспарантах. ф, рад
3
2
1
а)
0,1 0,2 0,3 0,4 0,5 0,6 г, мм
Рис. 2. Результаты моделирования действия бинарного аксикона с а = 85 мм-1: дискретизация фазы (показан фрагмент) с шагом 3,6 мкм (сплошная линия) и 28 мкм (точечная линия) (а) и соответствующие осевые распределения интенсивности (б)
Аксикон формирует очень неравномерное распределение интенсивности - среднеквадратичное отклонение от прямой линии на отрезке от 200 мм
до 220 мм составляет 44%. По рис. 3 видно, что в отличие от аксикона, рассматриваемый ДОЭ производит довольно «толстый» световой отрезок. В этом состоит особенность метода расчёта, т.к. отрезок формируется из набора фокусов, расположенных близко друг к другу. В этом случае, варьируя размеры и положение фокальных областей, можно создавать практически произвольные трёхмерные фигуры.
Рис. 3. Сравнение по равномерности осевого распределения интенсивности рассчитанного ДОЭ (сплошная линия на (а)), (6) и бинарного аксикона с а = 85 мм-1 (точечная линия на (а)), (в) при дискретизации фазы с шагом 28 мкм.
Поперечный размер на рис. (б) и (в) составляет 0,1 мм (негативное изображение)
При изменении длины волны освещающего излучения расстояние (2) изменится. В частности, если вместо 1 = 1,06 мкм использовать 1 = 0,633 мкм, то расстояние максимального самовоспроизведения световой нити увеличится в 12 » 1,67 раз. Однако при освещении ДОЭ излучением с длиной волны отличной от той, для которой был этот элемент изготовлен, высота микрорельефа становится неоптимальной. Оценить изменение высоты рельефа можно по формуле:
Ф2М =ф1(г) Л/12. (3)
Результаты моделирования для бинарного аксикона (1) с а = 85 мм-1 приведены на рис. 4а, по которому видно удлинение расстояния (2) и существенные изменения (в частности, уменьшение эффективности) в осевом распределении в связи с измене-
нием высоты микрорельефа по отношению к новой длине волны. Аналогичный результат получается для рассматриваемого ДОЭ (рис. 4б, для удобства восприятия показаны только распределения, соответствующие формируемому отрезку).
При использовании амплитудного оптического элемента негативный эффект, связанный с несогласованностью высоты рельефа и длины волны освещающего излучения, будет отсутствовать, однако общая доля энергии, поступающая в полезную область, уменьшится.
Результаты численного моделирования для оптимизированного ДОЭ, выполненного в виде амплитудного элемента и освещённого гауссовым пучком с радиусом перетяжки а = 2 мм, показаны на рис. 5 . Видно, что осевая интенсивность уменьшилась почти в 8 раз и несколько изменилась форма отрезка по сравнению тем, что формирует фазовый элемент. Также видно, что в этом случае изменение длины волны освещающего излучения приводит лишь к смещению отрезка и его масштабированию.
Рис. 4. Влияние уменьшения длины волны по сравнению с расчётной при освещении фазового ДОЭ: осевое распределение для расчётной длины волны 11 = 1,06 мкм (сплошная линия) и 12 = 0,633 мкм (точечная линия) для бинарного аксикона (а) и оптимизированного ДОЭ (б) (показаны только распределения, соответствующие формируемому отрезку)
Рис. 5. Влияние уменьшения длины волны по сравнению с расчётной при освещении амплитудного оптимизированного ДОЭ гауссовым пучком: осевое распределение для расчётной длины волны 11 = 1,06 мкм (сплошная линия), для 12 = 0,633 мкм (точечная линия) и для 13 = 0,532 мкм (пунктирная линия)
2. Экспериментальное формирование осевого отрезка с помощью амплитудного ДОЭ По радиальным отсчётам, рассчитанным в [21], было синтезировано изображение с шагом дискретизации 20 мкм (рис. 6а) и был изготовлен амплитудный фотошаблон (рис. 6б) на стеклянной подложке с нанесённым слоем хрома толщиной 100 нм.
Запись фотошаблона производилась на станции лазерной записи CLWS-200 путём создания защитной окисной плёнки окиси хрома на засвечиваемых участках с последующим удалением слоя хрома в жидкостном травителе.
Принцип действия станции основан на формировании в светочувствительном материале скрытого изображения будущих элементов. Система CLWS-200 работает в полярной системе координат. Запись топологии элементов ведётся в режиме непрерывного кругового сканирования под управлением компьютера. Во время записи рабочая заготовка, покрытая пленкой светочувствительного материала (в данном случае хрома), вращается относительно сфокусированного пучка света, который перемещается в радиальном направлении от периферии к центру. Рисунок на поверхности плёнки формируется с помощью компьютера путём программного управления мощностью излучения на каждом элементе поверхности.
Микрофотография центрального участка фотошаблона представлена на рис. 6б. Оптическая схема эксперимента представлена на рис. 7.
Рис. 6. Центральная часть бинарного элемента: амплитудная маска (а), вид изготовленного амплитудного ДОЭ (б), полученный с помощью высокоразрешающего сканера
Луч лазера L естественным образом расширяется (без коллиматора) до размеров дифракционного оптического элемента DOE, затем полученное распределение фиксируется CCD-камерой. В экспериментах было использовано два лазера с разными длинами волн: твердотельный лазер на неодимовом стекле с длиной волны 1,0б4 мкм и гелий-неоновый лазер длиной волны 0,б33 мкм.
При этом твердотельный лазер на неодимовом стекле генерировал лазерную моду TEM01, а гелий-неоновый лазер - моду TEM00.
2.1. Эксперимент с инфракрасным излучением
Распределения интенсивности, полученные при дифракции лазерного пучка с длиной волны 10б4 нм на разных расстояниях, представлены на рис. В, а изображения (диаметром 2 мм), полученные при моделировании, приведены на рис. 9. Из сравнения экспериментальных и модельных результатов видно общее согласование картины поперечной интенсивности.
Отметим, что элемент обладает оптической силой - освещение элемента пучком TEM10 привело к воспроизводству пучка TEM10 вдоль фокального отрезка несмотря на то, что элемент рассчитывался для фокусировки пучка TEM00.
Падающий на оптический элемент пучок имел расходимость в несколько градусов. В связи с этим в эксперименте протяжённость отрезка получилась несколько больше, чем предполагалось. Также данный факт привёл к некоторому смещению начала формирования отрезка дальше от плоскости оптического элемента.
РС
Рис. 7. Оптическая схема Расходящийся сферический фронт можно описать выражением:
V(r) = exP (ibr 2) . (4)
Для 1 = 1,0б мкм расчётный центр формируемого отрезка должен находиться на расстоянии
zc = 210 мм от оптического элемента, а в эксперименте центр сместился в плоскость ze = 235 мм . Определить параметр расходящегося сферического фронта b в (4) можно но разности Az = ze - zc:
к Az
b =
2z (z +Az)
(5)
Из выражения (5) получается Ь = 1,6 мм 2.
На рис. 10 показано сравнение продольных распределений интенсивности при использовании пучка ТЕМ01 излучения твердотельного лазера на неодимовом стекле в расчётном случае и при наличии слаборасходящегося сферического фронта.
• " - " -
200 мм 205 мм 210 мм 215 мм 220 мм 225 мм
* - * -
230 мм 235 мм 240 мм 245 мм 250 мм 255 мм Рис. 8. Экспериментально зафиксированные распределения интенсивности на разных расстояниях от ДОЭ (негатив)
" ((« »)
180 мм 190 мм 200 мм 210 мм 220 мм 230 мм 240 мм
Рис. 9. Поперечные картины интенсивности при моделировании для 1 = 1,064 мкм (негативное изображение)
а)
б)
Рис. 10. Продольные распределения интенсивности (негатив) при освещении амплитудного ДОЭ лазерным пучком ТЕМ01 с 1 = 1,064 мкм в расчётном случае (а) и при наличии слаборасходящегося сферического фронта Р = 1,6мм-2 (б) на отрезке оптической оси г є [180 мм, 280 мм], х є [-0,1 мм, 0,1 мм]
Для определения интенсивности на оптической оси проводилось измерение средней интенсивности в круге с радиусом, который бы полностью захватывал центральную часть пучка. Такое усреднение необходимо выполнять в связи со сложностью точного определения положения оптической оси в натурных экспериментах. На рис. 11 представлен график зависимости усреднённой интенсивности от расстояния от ДОЭ.
2.2. Эксперимент с красным излучением видимого диапазона Для гелий-неонового лазера распределения интенсивности, полученные при дифракции лазерного пучка, представлены на рис. 12, а соответствующие результаты моделирования - на рис. 13.
Для определения интенсивности на оптической оси в эксперименте проводилось измерение интенсивности в центре пучка. На рис. 13 представлен график зависимости этой интенсивности от расстояния от ДОЭ.
Рис. 11. Зависимость усреднённой интенсивности в центре пучка от расстояния при 1= 1,064 мкм: экспериментальные результаты (толстая линия) и результаты моделировании при наличии слаборасходящегося сферического фронта Р = 1,6ммГ2 (тонкая линия)
Как видно из рис. 13, продольный отрезок для длины волны 0,633 мкм удалился от ДОЭ и стал существенно длиннее, чем при использовании излучения с большей длиной волны.
Для 11 = 1,06 мкм (расчётная длина волны) середина отрезка приходилась на = 210 мм и протяжённость отрезка составляла Лг = 30 мм, для 12 = 0,633 мкм (длина волны освещающего излучения) увеличится в 1^ 12 » 1,67 раз, т.е. будет
г0 = 351 мм и Лг = 50 мм (рис. 5).
Отличие экспериментальных результатов от расчётных (рис. 13) также связано с наличием слабо-расходящегося сферического фронта. По формуле (5) параметр расходящегося сферического фронта примерно равен Р = 2,5 мм-2.
На рис. 14 показано сравнение продольных распределений интенсивности при формировании осевого отрезка при использовании излучения гелий-неонового лазера в расчётном случае и при наличии слаборасходящегося сферического фронта.
0 © \3sJ- Ь />-ч\ -• '
250 мм 260 мм 270 мм 280 мм 290 мм 3 00 мм
V Г • \ Ф) .
310 мм 320 мм 330 мм 340 мм 350 мм 360 мм
Ф' '© г ♦ // , • Ф ^ | ‘ А. *
370 мм 380 мм 390 мм 400 мм 410 мм 420 мм
кл 1 <3>
430 мм 440 мм 450 мм 460 мм 470 мм 480 мм
Рис. 12. Экспериментально зафиксированные распределения интенсивности на разных расстояниях от ДОЭ для длины волны 0,633 мкм (негативное изображение)
300 мм 320 мм 340 мм 360 мм 380 мм 400 мм
Рис. 13. Поперечные картины интенсивности при моделировании для 12 = 0,633 мкм (негативное изображение)
6)________________________________________________________________________________________________________________________
Рис. 14. Продольные распределения интенсивности (негатив) при освещении амплитудного ДОЭ лазерным пучком ТЕМ00 с 1 = 0,633 мкм в расчётном случае (а) и при наличии слаборасходящегося сферического фронта Р = 2,5 мм- (б) на отрезке оптической оси z е [300мм, 500мм], х е [-0,1 мм, 0,1 мм]
Для определения в эксперименте интенсивности на оптической оси проводилось измерение средней интенсивности в круге с радиусом, который бы полностью захватывал центральную часть пучка. На рис. 15 представлен график зависимости усреднённой интенсивности от расстояния от ДОЭ.
Рис. 15. Зависимость усреднённой интенсивности в центре пучка от расстояния при 1 = 0,633 мкм: экспериментальные результаты (толстая линия) и результаты моделировании при наличии слаборасходящегося сферического фронта Р = 2,5 мм2 (тонкая линия)
Заключение В работе проведено численное и экспериментальное исследование характеристик бинарного ДОЭ, предназначенного для формирования осевого отрезка, в условиях изменения параметров освещающего пучка.
Основное преимущество оптимизированного стохастическим методом ДОЭ для получения равномерного распределения интенсивности на заданном отрезке перед бинарными аксиконами состоит в возможности формирования осевого отрезка близко от плоскости ДОЭ при достаточно крупном шаге дискретизации микрорельефа оптического элемента, что важно, в частности, при реализации фазовых распределений на динамических транспарантах.
При исследовании оптимизированного бинарного ДОЭ рассмотрено три типа освещающего пучка: ТЕМ00, ТЕМ01 и пучок с равномерным распределением интенсивности, ограниченный круглой апертурой радиуса ДОЭ. При замене гауссового (расчётного) пучка равномерным интенсивность на заданном отрезке становится менее равномерной (среднеквадратичное отклонение от прямой линии увеличивается с 5% до 19%), но осевая интенсивность увеличивается (в рассмотренном случае в полтора раза). Отметим, что элемент обладает оптической силой - освещение элемента пучком ТЕМ10 привело к воспроизводству пучка ТЕМ10 вдоль фокального отрезка несмотря на то, что элемент рассчитывался для фокусировки пучка ТЕМ00.
Также проведено моделирование изменения длины волны в случае исполнения ДОЭ в виде фазового или амплитудного оптического элемента. При отличии длины волны освещающего излучения (12) от расчётной (11) протяжённость формируемого светового отрезка изменится 12 раз, т.е. увеличится при уменьшении длины волны (осевая интенсивность будет пропорционально уменьшаться). Однако при использовании фазового ДОЭ высота микрорельефа становится неоптимальной. При ис-
пользовании амплитудного оптического элемента негативный эффект, связанный с несогласованностью высоты рельефа и длины волны освещающего излучения, будет отсутствовать.
Для предварительной экспериментальной верификации проведённых численных исследований был изготовлен амплитудный ДОЭ. Анализ полученных экспериментальных результатов (при использовании лазерного излучения с различными длинами волн) с учётом небольшой расходимости волнового фронта освещающего пучка показал хорошее согласование с результатами моделирования.
Благодарности
Работа выполнена при поддержке российско-американской программы «Фундаментальные исследования и высшее образование» (грант CRDF PG08-014-1), грантов РФФИ 10-07-00109-а, 10-07-00438-а, гранта Президента РФ поддержки ведущих научных школ НШ-7414.2010, гранта президента РФ МД-8026.2010.2 и Государственного контракта № 02.740.11.0805
Литература
1. Wang, K Influence of the incident wave-front on intensity distribution of the nondiffracting beam used in large-scale measurement / K. Wang, L. Zeng, and Ch. Yin // Opt. Commun. - 2003. - Vol. 216. - P. 99-103.
2. Lu, J.-Y. Diffraction-limited beams and their applications for ultrasonic imaging and tissue characterization / Jian-yu Lu, J. F. Greenleaf // Proc. SPIE. - 1992. - Vol. 1733. -P. 92-119.
3. Leitgeb, R.A. Extended focus depth for Fourier domain optical coherence microscopy / R.A. Leitgeb, M. Villiger,
A.H. Bachmann, L. Steinmann, and T. Lasser // Opt Lett.
- 2006. - Vol. 31, N 16. - P. 2450-2452.
4. Lee, K.-S. Bessel beam spectral-domain high-resolution optical coherence tomography with micro-optic axicon providing extended focusing range / Kye-Sung Lee, Jannick P. Rolland // Opt Lett - 2008. - Vol. 33, N 15. - P. 1696-1698.
5. Lu, J.-Y. Producing deep depth of field and depth-independent resolution in NDE with limited diffraction beams / Jian-yu Lu, J. F. Greenleaf // Ultrason. Imag. -1993. - Vol. 15, N 2. - P. 134-149.
6. Lunazzi, J.J. Photographing by means of a diffractive axicon / Jose J. Lunazzi, Daniel S. F. Magalhaes // XXIX ENFMC - Annals of Optics. - 2006. - P. 1-4.
7. Arimoto, R. Imaging properties of axicon in a scanning optical system / R. Arimoto, C. Saloma, T. Tanaka, and S. Kawata // Appl. Opt. - 1992. - Vol. 31, N 31. -P. 6653-6657.
8. Fortin, M. Optical tests with Bessel beam interferometry / Mathieu Fortin, Michel Piche and Ermanno F. Borra // Optics Express. - 2004. - Vol. 12, N 24. - P. 5887-5895.
9. Reichelt, S. Self-calibration of wavefront testing interferometers by use of diffractive elements / Stephan Reichelt, Hans Tiziani, Hans Zappe // Proc. of SPIE. - 2006. -Vol. 6292. - P. 629205.1-629205.10.
10. McLeod, J.H. The axicon: a new type of optical element // J. Opt. Soc. Am. - 1954. - Vol. 44. - P. 592-597.
11. Голуб, М.А. Фокусировка излучения в заданную область пространства / М. А. Голуб, С.В. Карпеев, А.М. Прохоров, И.Н. Сисакян, В.А. Сойфер // Письма в ЖТФ. - 1981. - Т. 7, № 10. - С. 618-623.
12. Голуб, М. А. Вычислительный эксперимент с элементами плоской оптики / М.А. Голуб, Н.Л. Казанский, И.Н. Сисакян, В. А. Сойфер // Автометрия. - 1988. -№ 1. - С. 70-82.
13. Kazanskii, N.L. Correction of focuser phase function by computer-experimental methods // Computer Optics. Per-gamon Press. - 1989. - Vol. 1, N 1. - P. 69-73.
14. Doskolovich, L.L. Analysis of quasiperiodic and geometric optical solutions of the problem of focusing into an axial segment / L.L. Doskolovich, N.L. Kazanskiy, V.A. Soi-fer, A.Ye. Tzaregorodtzev // Optik. - 1995. - Vol. 101, N 2. - P. 37-41.
15. Котляр, В.В. Дифракционный расчёт фокусаторов в продольный отрезок / В.В. Котляр, В.А. Сойфер, С.Н. Хонина // Письма в ЖТФ. - 1991. - Т. 17, № 24. -С. 63-66.
16. Gerchberg, R.W. A practical algorithm for the determination of the phase from image and diffraction plane pictures / R.W. Gerchberg, W.O. Saxton // Optik. - 1972. -Vol. 35. - P. 237.
17. Методы компьютерной оптики / под ред. В.А. Сойфе-ра. - 2-е изд., испр. - М.: Физматлит, 2003. - 688 с.
18. Волков, А.В. Бинарный дифракционный оптический элемент для фокусировки гауссового пучка в продольный отрезок / А.В. Волков, В.В. Котляр, О.Ю. Моисеев, О.Е. Рыбаков, Р.В. Скиданов, В.А. Сойфер, С.Н Хонина // Оптика и спектроскопия. - 2000. - Т. 89, № 2. -С. 347-352.
19. Soifer, V.A. Multifocal diffractive elements / V.A. Soifer, L.L. Doskolovich, N.L. Kazanskiy // Optical Engineering.
- 1994. - Vol. 33, № 11. - P. 3610-3615.
20. Сойфер, В.А. Синтез бинарного фокусатора в произвольную кривую в электромагнитном приближении /
В.А. Сойфер, Н.Л. Казанский, С.И. Харитонов // Компьютерная оптика. - 1996. - № 16. - C. 22-27.
21. Качалов, Д.Г. Оптимизация бинарных ДОЭ для формирования продольных распределений интенсивности / Д.Г. Качалов, В.С. Павельев, С.Н. Хонина // Компьютерная оптика. - 2009. - Т. 33, № 4. - С. 441-445.
22. Хонина, С.Н. Сравнительный анализ распределений интенсивности, формируемых дифракционным аксиконом и дифракционным логарифмическим ак-сиконом / С.Н. Хонина, С.А. Балалаев // Компьютерная оптика. - 2009. - T. 33, № 2. - С. 162-174.
References
1. Wang, K. Influence of the incident wave-front on intensity distribution of the nondiffracting beam used in large-scale measurement / K. Wang, L. Zeng, and Ch. Yin // Opt. Commun. - 2003. - Vol. 216. - P. 99-103.
2. Lu, J.-Y. Diffraction-limited beams and their applications for ultrasonic imaging and tissue characterization / Jian-yu Lu, J. F. Greenleaf // Proc. SPIE. - 1992. - Vol. 1733. - P. 92-119.
3. Leitgeb, R.A. Extended focus depth for Fourier domain optical coherence microscopy / R.A. Leitgeb, M. Villiger, A.H. Bachmann, L. Steinmann, and T. Lasser // Opt Lett.
- 2006. - Vol. 31, N 16. - P. 2450-2452.
4. Lee, K.-S. Bessel beam spectral-domain high-resolution optical coherence tomography with micro-optic axicon providing extended focusing range / Kye-Sung Lee, Jan-nick P. Rolland // Opt. Lett. - 2008. - Vol. 33, N 15. -P. 1696-1698.
5. Lu, J.-Y. Producing deep depth of field and depth-independent resolution in NDE with limited diffraction beams / Jian-yu Lu, J. F. Greenleaf // Ultrason. Imag. 1993. Vol. 15, № 2. P. 134-149.
6. Lunazzi, J.J. Photographing by means of a diffractive axicon / Jose J. Lunazzi, Daniel S. F. Magalhaes // XXIX ENFMC - Annals of Optics. - 2006. - P. 1-4.
7. Arimoto, R. Imaging properties of axicon in a scanning optical system / R. Arimoto, C. Saloma, T. Tanaka, and
S. Kawata // Appl. Opt. - 1992. - Vol. 31, N 31. -P. 6653-6657.
8. Fortin, M. Optical tests with Bessel beam interferometry / Mathieu Fortin, Michel Piche and Ermanno F. Borra // Optics Express. - 2004. - Vol. 12, N 24. - P. 5887-5895.
9. Reichelt, S. Self-calibration of wavefront testing interferometers by use of diffractive elements / Stephan Reichelt, Hans Tiziani, Hans Zappe // Proc. of SPIE. - 2006. -Vol. 6292. - P. 629205.1-629205.10.
10. McLeod, J.H. The axicon: a new type of optical element // J. Opt. Soc. Am. - 1954. - Vol. 44. - P. 592-597.
11. Golub, M.A. Radiation focusing in the set area of space / M.A. Golub, S.V. Karpeev, A.M. Prokhorov, I.N. Sisakyan, V.A. Soifer // JETP Letters. - 1981. - Vol. 7, No. 10. -P.618-623. - (in Russian).
12. Golub, M.A. Computing experiment with elements of plain optics / M.A. Golub, N.L. Kazanskii, I.N. Sisakyan, V.A. Soifer // Avtometria. - 1988. - N 1. - P. 70-82. - (in Russian).
13. Kazanskii, N.L. Correction of focuser phase function by computer-experimental methods // Computer Optics. Per-gamon Press. - 1989. - Vol. 1, N 1. - P. 69-73.
14. Doskolovich, L.L. Analysis of quasiperiodic and geometric optical solutions of the problem of focusing into an axial segment / L.L. Doskolovich, N.L. Kazanskiy, V.A. Soifer, A.Ye. Tzaregorodtzev // Optik. - 1995. - Vol. 101, N 2. - P. 37-41.
15. Kotlyar V.V. Calculation of the focusators into a longitudinal linesegment / V.V. Kotlyar, V.A. Soifer, S.N. Kho-nina , JETP Letters. - 1991. - Vol. 17, No. 24. - P. 63-66.
- (in Russian).
16. Gerchberg, R.W. A practical algorithm for the determination of the phase from image and diffraction plane pictures / R.W. Gerchberg, W.O. Saxton // Optik. - 1972. -Vol. 35. - P. 237.
17. Methods of Computer Optics (Secondary Edition) / edited by V.A. Soifer - Moscow: “Fizmatlit” Publisher, 2003. -688 p. - (in Russian).
18. Volkov, A.V. Binary Diffraction Optical Element Focusing a Gaussian Beamto a Longitudinal Segment / A.V. Volkov, V.V. Kotlyar, O.V. Moiseev, O.E. Rybakov, R.V. Skidanov,V.A. Soifer, and S.N. Khonina // Optics and Spectroscopy. - 2000. - Vol. 89. - No.2. -P.347-352. - (In Russian).
19. Soifer, V.A. Multifocal diffractive elements / V.A. Soifer, L.L. Doskolovich, N.L. Kazanskiy // Optical Engineering.
- 1994. - Vol. 33, № 11. - P. 3610-3615.
20. Soifer, V.A. Synthesis of a binary focusator in any curve in electromagnetic approach / V.A. Soifer, N.L. Kazanskiy, S.I. Kharitonov // Computer optics. - 1996. - N 16. -P. 22-27. - (In Russian).
21. Kachalov, D.G. Optimization of binary DOE forming intensity distribution along an axial focal zone / D.G. Kachalov, V.S. Pavelyev, S.N. Khonina // Computer optics. -2009. - V. 33, N 4. - P. 441-445. - (In Russian).
22. Khonina, S.N. The comparative analysis of the intensity distributions formed by diffractive axicon and diffractive logarithmic axicon / S.N. Khonina, S.A. Balalaev // Computer optics. - 2009. - V. 33, N 2. - P. 162-174. -(In Russian).
INVESTIGATION OF THE AXIAL LINE FORMATION BY OPTIMIZED BINARY DOE
S.N. Khonina 12, R.V. Skidanov 12, D.G. Kachalov 2, V.S. Pavelyev 12, О.У. Moiseev 1 1 Institution of Russian Academy of Sciences, Image Processing Systems Institute RAS,
2 Samara State Aerospace University named after S.P. Korolev
Abstract
Numerical and experimental investigation of the optimized binary DOE intended to form an axial line with uniform intensity is realized for various parameters of an illuminating beam.
Key words: binary axicon, diffractive optical element (DOE), axial light line, laser beam, phase and amplitude DOE.
Сведения об авторах
Хонина Светлана Николаевна, доктор физико-математических наук, профессор Самарского государственного аэрокосмического университета имени академика С.П. Королёва; ведущий научный сотрудник Учреждения Российской академии наук Институт систем обработки изображений РАН. E-mail: khonina@smr.ru .
Область научных интересов: дифракционная оптика, сингулярная оптика, модовые и поляризационные преобразования, оптическое манипулирование, оптическая и цифровая обработка изображений.
Svetlana Nikolaevna Khonina, Doctor of Physical and Mathematical Sciences; Professor of the Samara State Aerospace University named after S.P. Korolev. Leading researcher of the Image Processing Systems Institute of the RAS. E-mail: khonina@smr.ru.
Research interests: diffractive optics, singular optics, mode and polarization transformations, optical manipulating, optical and digital image processing.
Скиданов Роман Васильевич, 1973 года рождения. В 1990 году с отличием окончил Самарский государственный университет (СамГУ) по специальности «Физика». Доктор физико-математических наук (2007 год), работает старшим научным сотрудником лаборатории лазерных измерений Института систем обработки изображений РАН (ИСОИ РАН), профессором кафедры технической кибернетики СГАУ. Скиданов Р.В. - специалист в области дифракционной оптики, математического моделирования, обработки изображений и нанофотоники. В списке научных работ Р.В. Скиданова 60 статей, 3 монографии. E-mail: romans@smr.ru .
Roman Vasilevich Skidanov (b. 1973) graduated with honours (1990) from the Samara State University (SSU)), majoring in Physics. He received his Doctor in Physics & Maths (2007) degrees from Samara State Aerospace University. He is the senior researcher laboratory of laser measurement at the Samara Image Processing Systems Institute of the Russian Academy of Sciences (IPSI RAS), holding a part-time position of professor at SSAU’s Technical Cybernetics sub-department. He is co-author of 60 scientific papers, 3 monographs. His current research interests include diffractive optics, mathematical modeling, image processing, and nanophotonics. E-mail: ro-mans@smr.ru.
Качалов Денис Георгиевич, инженер кафедры наноинженерии, м.н.с. научнообразовательного центра «Нанотехнологии» Самарского государственного аэрокосмического университета имени академика С.П. Королёва. E-mail: kachalov@ssau.ru. Область научных интересов: дифракционная оптика, методы стохастической оптимизации.
Denis Georgievich Kachalov, engineer of Nanoengineering Department and researcher of Research and Educational Center “Nanotechnology” of Samara State Aerospace University named after S.P. Korolev. E-mail: kachalov@ssau.ru. Research interests: diffractive optics, methods of stochastic optimization.
Павельев Владимир Сергеевич, доктор физико-математических наук, заведующий кафедрой наноинженерии Самарского государственного аэрокосмического университета имени академика С.П. Королва; главный научный сотрудник Учреждения Российской академии наук Институт систем обработки изображений РАН.
E-mail: pavelyev@smr.ru.
Область научных интересов: дифракционная микрооптика и оптические волноводы, нанофотоника и синтез оптических метаматериалов, технологии формирования микро- и наноструктур.
Vladimir Sergeevich Pavelyev, Doctor of Physical and Mathematical Sciences; Head of Nanoengineering Department of the Samara State Aerospace University named after S.P. Korolev, chief researcher of the Image Processing Systems Institute of the RAS. E-mail: pavelyev@smr.ru. Research interests: diffractive microoptics, optical waveguides, nanophotonics, optical metamaterials synthesis, technologies for forming of micro- and nanostructures.
Моисеев Олег Юрьевич, кандидат технических наук, старший научный сотрудник Учреждения Российской академии наук Институт систем обработки изображений РАН. E-mail: moiseev@smr.ru.
Область научных интересов: технологические процессы создания микроструктур дифракционной оптики.
Oleg Yurevich Moiseev, Candidate of Technical Sciences; senior staff scientist of the Image Processing Systems Institute of the RAS. E-mail: moiseev@smr.ru.
Research interests: technological process of diffractive optics.
Поступила в редакцию 31 августа 2010 г.