CC BY
31
ОПТИМИЗАЦИЯ БИНАРНЫХ ДОЭ ДЛЯ ФОРМИРОВАНИЯ ПРОДОЛЬНЫХ РАСПРЕДЕЛЕНИЙ ИНТЕНСИВНОСТИ
Денис Георгиевич Качалов2 (инженер, e-mail: kachalov@ssau.ru), Владимир Сергеевич Павелъев1'2(главный научный сотрудник, профессор, e-mail: pavelyev@smr. ru), Светлана Николаевна Хонина12 (ведущий научный сотрудник, профессор, e-mail: khonina@smr. ru)
1 Учреждение Российской академии наук Институт систем обработки изображений РАН,
2 Самарский государственный аэрокосмический университет им академика С.П. Королева
Аннотация
Рассмотрен стохастический подход к оптимизации квантованных дифракционных оптических элементов (ДОЭ), предназначенных для формирования заданных продольных распределений интенсивности.
Ключевые слова: дифракционные оптические элементы, углубленный фокус, генетические алгоритмы.
Введение
Элементы дифракционной микрооптики, формирующие заданные продольные распределения интенсивности (соосные отрезки, наборы последовательных осевых фокусов и т.д.), находят применение во многих областях техники. Например, при выравнивании и совмещении конструкций на больших расстояниях [1]; при формировании изображений протяженных или движущихся объектов, что особенно актуально в медицине [2-4], а также для неразрушающего исследования материалов [5] и устройств [6], в метрологии для сканирующих [7] и интерференционных [8, 9] оптических систем.
Классическим рефракционным оптическим элементом для генерации осевой световой нити является аксикон [10]. Однако распределение интенсивности вдоль оси, формируемое аксиконом, очень неравномерное - кроме существенных флуктуаций наблюдается линейный рост интенсивности при увеличении расстояния от оптического элемента.
Оптимизация фазовой функции дифракционного оптического элемента, фокусирующего лазерное излучение в продольный отрезок, в рамках геомет-рооптического подхода рассматривалась в работах [11-13]. Отметим, что геометрооптический расчет принципиально не позволяет учесть дифракционных эффектов.
Также известны итерационные методы расчета ДОЭ [14], основанные на алгоритме Герчберга-Сакстона [15], которые позволяют учитывать дифракционные эффекты, однако не обладают сходимостью (может быть показана лишь нерасходимость). Кроме того, необходимость применения обратного оператора на каждой итерации метода затрудняет получение квантованной фазы оптического элемента, в то время как значительная часть технологий (главным образом литографических) ориентирована как раз на получение квантованного микро- нанорельефа [16]. Поэтому в работе [17] для синтеза бинарного рельефа был применен метод последовательных приближений с бинаризацией фазовой функции на каждом шаге. Полученный ДОЭ имел вид кольцевой решетки из 186 колец,
причем на периферии период был близок к 10 мкм при радиусе элемента 4 мм.
В случае, когда шаг поперечной дискретизации должен быть не меньше определенной величины, и соответственно, количество колец в структуре ДОЭ желательно свести к минимуму, предпочтительнее применение стохастических алгоритмов прямого поиска, так как поиск функции высоты микрорельефа можно организовать на множестве квантованных решений, к тому же алгоритмы прямого поиска обладают сходимостью.
В работах [18, 19] обсуждался расчет бинарных ДОЭ, формирующих множество фокусов, а также фокусирующих в заданные кривые в электромагнитном приближении. Однако при формировании осевых световых распределений, удовлетворяющих параксиальному приближению, применение скалярной теории дифракции является вполне корректным и ресурсосберегающим подходом.
Данная работа посвящена стохастической оптимизации квантованных дифракционных оптических элементов, предназначенных для формирования заданных продольных распределений интенсивности.
1. Основной формализм
Рассмотрим постановку задачи расчета ДОЭ, формирующего продольное осевое распределение интенсивности из лазерного пучка с гауссовым распределением интенсивности в сечении (рис. 1).
т\
ДОЭ
X, а
R
V
Рис. 1. Оптическая схема формирования продольного распределения интенсивности
Для расчета комплексной амплитуды Е (р, г) на
осевом расстоянии г и радиальном расстоянии р используем интеграл Френеля [20]
^ (р, 2 )=-^ек ехр (
2 I 2 2
€ ехр I -02)ехр
'. , л кг2 ф( г)+1—
Л
кгр
(1)
г бт,
где Я - радиус ДОЭ, Л0 (х) - функция Бесселя первого рода нулевого порядка, ф (г) - фаза ДОЭ, о - радиус перетяжки освещающего гауссова пучка, к = 2р1 - волновое число.
Запишем формулу для интенсивности на оптической оси
I (0, 2) = (0, 2)Г
1к
^ (0,2) =--ехр(/кг) х
1оеХР I )ехр
кг2" /ф( г) + / —
(2)
(3)
г &г.
Для формирования соосного отрезка представляется целесообразным выбрать в качестве минимизируемого функционала невязку:
Ф(ф( г )) = £|1 ( )-I ( 2,. )• С^ , (4)
,=1
где I (21) - интенсивность, сформированная ДОЭ с фазовой функцией ф(г) на оптической оси на расстоянии , I () - требуемый уровень интенсивности на оптической оси на расстоянии 2, , рассчитанный из соображений закона сохранения энергии, N - число контролируемых точек на оптической оси, Ст - энергетический коэффициент, определяемый количеством уровней квантования т.
Для оценки оптимизируемого ДОЭ применялись два критерия: энергетическая эффективность и среднеквадратичное отклонение (СКО). Энергетическая эффективность рассчитывалась по следующей формуле:
е = (£I()I()) '. (5)
СКО рассчитывалось по формуле:
5 =
_N. __о
Е (I ( )-1 ( ))2
X12 ( )
(6)
2. Выбор алгоритма оптимизации
Существует несколько подходов к организации оптимизации микрорельефа оптического элемента методами прямого поиска [21]. Прежде всего, это алгоритм полного перебора, когда производится поиск по всей области допустимых решений без каких-либо ограничений. Этот метод хорошо применим для плохо формализованных задач. Минусы
его очевидны - это большие вычислительные затраты. Популярен также алгоритм имитации отжига [22], однако его особенности таковы, что по мере выполнения итераций все менее вероятным становится выход из локального минимума. Также существуют эволюционные алгоритмы, основанные на моделировании процесса естественной эволюции. В данной работе был использован генетический алгоритм, предложенный Джоном Холландом [23]. Этот алгоритм работает с группами решений, что является его отличительной особенностью и дает возможность выхода из локальных минимумов при разумных вычислительных затратах.
Принцип работы алгоритма заключается в создании начальной группы решений - популяции (родительского пула), состоящей из кодированных наборов параметров оптимизации задачи, называемых хромосомами, и применении к популяции хромосом генетических операторов. Наиболее популярные генетические операторы - это рекомбинация и мутация, которые увеличивают количество хромосом в популяции, и селекция, выбирающая по определенному принципу лучшие хромосомы. Основанием для оценки качества хромосомы служит значение функции приспособленности или функционала, подлежащего оптимизации. Ранее [24] было показано, что генетический алгоритм эффективен для решения задач подобного класса.
В данном исследовании была применена одноточечная и двухточечная рекомбинация (взаимоисключающие равновероятные) с вероятностью 85%. При этом точка (точки) среза выбиралась произвольно на протяжении всей длины хромосомы. Выбор хромосом для применения оператора рекомбинации был принят случайным, что обеспечило большую рассредоточенность направления оптимизации, тем самым снизив вероятность попадания в локальные минимумы.
Применялся стандартный оператор мутации с вероятностью 3% для каждого гена хромосомы, выбранной для применения данного оператора.
Оператор селекции требует более подробного рассмотрения в силу того, что существует множество стратегий селекции. Наиболее распространенная и простая в реализации - это элитарная селекция [21], приводящая к обеднению генетического материала популяции и, как следствие, к стагнированию алгоритма, попаданию в локальный минимум. Более предпочтительной для решения вышеобозначенной задачи выглядит турнирная селекция, которая заключается в случайном выборе нескольких хромосом в отдельную группу и анализе этих хромосом с точки зрения функции приспособленности. Такая стратегия, несмотря на несколько большую вычислительную нагрузку, обеспечивает высокую эффективность работы алгоритма, предотвращает вырождение генетического материала и, следовательно, снижает вероятность стагнации. В качестве дополнительного средства предотвращения стагнации среди хромосом группы с небольшой веро-
х
г
х
,=1
,=1
ятностью (2%) в следующее поколение отбиралась хромосома, следующая за наилучшей.
Для оценки приспособленности хромосомы можно применить функционал (4), но, учитывая специфику данной задачи (в частности, снижение интенсивности излучения вдоль оптической оси по мере удаления от источника света), функционал был изменен:
Ф(ф(г )) = I
(
ехр
I(г)-1 (г)• С I (г.)• С
\ I / т
\
-1
(7)
3. Результаты численных экспериментов
В данной работе решалась задача оптимизации микрорельефа бинарного радиально-симметричного ДОЭ с радиусом Я = 3,6 мм и количеством радиальных отсчетов М = 128, фокусирующего гауссовый пучок с параметрами X = 1,06 мкм, о = 2 мм в десять равномерно распределенных фокусов на расстоянии от Е = 200 мм до Е2 = 220 мм, составляющих высокоэнергетический световой отрезок.
Изменение значения функционала (7) с ходом итераций показано на рис. 2.
|Ф(г), ед
45 40 35 30 25 20 15 10 5 О
1000 3000 5000 7000 1,итер.
Рис. 2. Зависимость значения функционала от номера итерации
В процессе оптимизационной процедуры был получен микрорельеф ДОЭ, обеспечивающий следующие показатели качества распределения интенсивности, приведенной на рис. 3: е = 4,58 %, 5=2,07 %.
Вид радиального сечения рассчитанного бинарного фазового распределения приведен на рис. 4.
202 206 210 216 220 г, мм
Рис. 3. Распределение интенсивности вдоль оптической оси
4. Сравнение полученных результатов с традиционными оптическими элементами Для решения задачи фокусировки в продольный осевой отрезок также применяются ДОЭ, имеющие
аналитические представления фазовой функции. К таким можно отнести конический аксикон [10]. К преимуществам можно отнести простоту расчета фазовой функции ДОЭ, а к недостаткам - неизбежные погрешности квантования при использовании литографических технологий [25] и неравномерное распределение интенсивности вдоль оптической оси. Фазовая функция конического аксикона имеет вид:
ф(г) = шоа2р (-аг), (8)
где а - параметр аксикона.
111111111
•о •о "О "О 4-1 г
ГЧ ГЧ гч Гч <4 <ч
"-ч 'о "■ч "ч "-ч ММ
ОО «Л «г> 20 "О °0 го
^ > "о 0\ о
С о IX «п
С? о" сГ »-Г "-Г •-Г к
Рис. 4. Фаза бинарного ДОЭ, формирующего осевое распределение интенсивности
Известен также дифракционный аксикон [25], фазовая функция которого соответствует фазе бесселевой моды:
ф(г) = а^[ Jo (аг)],
(9)
где J0 (г) - Бессель функция первого рода нулевого порядка, а - масштабирующий множитель. Основное преимущество ДОЭ (9) - это бинарный профиль фазы, а недостаток такой же, как у конического аксикона -неравномерное распределение на осевом отрезке.
На рис. 5 приведено распределение интенсивности, сформированное коническим аксиконом (8) с параметром а = 30 и квантованной по двум уровням фазой; распределение интенсивности, сформированное дифракционным аксиконом (9) с масштабирующим множителем а = 30 (следует отметить, что фазовый профиль данного оптического элемента по определению является двухуровневым, принимая значения 0 или р, так что квантование не требуется); и распределение интенсивности, сформированное рассчитанным бинарным ДОЭ.
Из рис. 5 видно, что распределение интенсивности вдоль оптической оси, сформированное рассчитанным ДОЭ, является более энергетически выгодным, а также обеспечивает более равномерное распределение, чем распределения, сформированные другими оптическими элементами.
Заключение
Полученные результаты численного эксперимента позволяют сделать вывод об эффективности применения алгоритмов прямого поиска к задаче оптимизации квантованных оптических элементов, формирующих продольное распределение интенсивности. Сравнение с известными оптическими элементами, такими как бинарный вариант конического аксикона и дифракционный аксикон, показывает выгодное отличие рассчитанного в работе бинарного ДОЭ как по энергетической эффективности, так и по равномерности распре-
¿=1
деления. Также выгодное отличие стохастической дискретной оптимизации от прочих методов расчета в том, что решение ищется на пространстве технологически реализуемых решений, что делает поиск действительно прямым, исключая промежуточное кванто-
120 100 80 60 40 20 О
Литература
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, Jian-yu 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. Stein-mann, and T. Lasser // Opt. Lett. 2006, - Vol. 31(16). - P. 2450-2452.
4. Lee, Kye-Sung 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(15). - P. 1696-1698.
5. Lu, Jian-yu 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, José J. Photographing by means of a diffractive axicon / José J. Lunazzi, Daniel S. F. Magalhâes // 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(31). - P. 6653-6657.
8. Fortin, M. Optical tests with Bessel beam interferometry / Mathieu Fortin, Michel Piché and Ermanno F. Borra // Optics Express - 2004. - Vol. 12(24). - P. 5887-5895.
9. Reichelt, S. Self-calibration of wavefront testing interferometers by use of diffractive elements / Stephan Reichelt, Hans Tiziani, and 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. Kazanskii N.L. Correction of focuser phase function by computer-experimental methods // Computer Optics. -Pergamon Press, 1989, Vol. 1, No. 1. - P.69-73.
13. Analysis of quasiperiodic and geometric optical solutions of the problem of focusing into an axial segment / Dosko-lovich L.L., Kazanskiy N.L., Soifer V.A., Tzaregorodtzev A.Ye. // Optik. - 1995. - Vol.101, № 2. - P.37-41.
14. Котляр В.В. Дифракционный расчет фокусаторов в продольный отрезок / Котляр В.В., Сойфер В.А., Хонина С.Н. // Письма в ЖТФ. - 1991. - Т. 17. - № 24. -С. 63-66
15. 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.
16. Методы компьютерной оптики (Издание второе, исправленное) / под ред. В. А. Сойфера - М.: Физматлит, 2003. - 688 с.
17. Волков А.В. Бинарный дифракционный оптический элемент для фокусировки гауссового пучка в продольный отрезок / Волков А. В., Котляр В. В., Моисеев О.Ю., Рыбаков О.Е., Скиданов Р.В., Сойфер В.А., Хо-нина С.Н. // Оптика и спектроскопия. - 2000. - Т. 89. -№ 2. - С. 347-352.
18. Soifer V.A., Doskolovich L.L., Kazanskiy N.L. Multifocal diffractive elements // Optical Engineering. - 1994. - Vol.33, № 11. - P.3610-3615.
19. Сойфер В.А., Казанский Н.Л., Харитонов С.И. Синтез бинарного фокусатора в произвольную кривую в электромагнитном приближении // Компьютерная оптика, 1996, № 16, с. 22-27.
20. Дифракционная компьютерная оптика / под ред. В.А. Сойфера - М.: Физмалит, 2007. - Гл. 4.
21. Meister, M. Novel approaches to direct search algorithms for the design of diffractive optical elements / M. Meister, R.J. Winfield // Optics Communications 2002. - Vol. 203. - P. 39-49.
22. Kirkpatrick S. Optimization by Simulated Annealing / S. Kirkpatrick, C. D. Gelatt and M. P. Vecchi // Science -1983. - Vol. 220. - No. 4598. - P. 671-680.
23. Holland, J. H. Adaptation in Natural and Artificial Systems / University of Michigan Press, M.I., 1975.
100 120 140 160 180 200 220 240 260 280 300 320 z, мм
Рис. 5. Результат сравнения рассчитанного ДОЭ с известными аналогами
24. Pavelyev, V.S. Technology-oriented stochastic optimization of radially symmetric DOEs//Proceedings SPIE 2005. - Vol. 5962 - P. 891-897.
25. Котляр В.В. Фазовые оптические элементы для формирования квазимод свободного пространства / Котляр В.В., Сойфер В.А., Хонина С.Н. // Квантовая электроника. - 1991. - Т. 18. № 11. - С. 1391-1394.
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, Jian-yu 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. Stein-mann, and T. Lasser // Opt. Lett. 2006, - Vol. 31(16). - P. 2450-2452.
4. Lee, Kye-Sung 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(15). - P. 1696-1698.
5. Lu, Jian-yu 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, José J. Photographing by means of a diffractive axicon / José J. Lunazzi, Daniel S. F. Magalhâes // 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(31). - P. 6653-6657.
8. Fortin, M. Optical tests with Bessel beam interferometry / Mathieu Fortin, Michel Piché and Ermanno F. Borra // Optics Express - 2004. - Vol. 12(24). - P. 5887-5895.
9. Reichelt, S. Self-calibration of wavefront testing interferometers by use of diffractive elements / Stephan Reichelt, Hans Tiziani, and 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. Golub M.A., Karpeev S.V., Prokhorov A.M., Sisakyan I.N., Soifer V.A. // JETP Letters. - 1981. -Vol.7. - No. 10. - P.618-623. - (in Russian).
12. Kazanskii N.L. Correction of focuser phase function by computer-experimental methods // Computer Optics. -Pergamon Press, 1989, Vol. 1, No. 1. - P.69-73.
13. Analysis of quasiperiodic and geometric optical solutions of the problem of focusing into an axial segment / Dosko-lovich L.L., Kazanskiy N.L., Soifer V.A., Tzaregorodtzev A.Ye. // Optik. - 1995. - Vol.101, № 2. - P.37-41.
14. Kotlyar V.V. Calculation of the focusators into a longitudinal linesegment / Kotlyar V.V., Soifer V.A., Khonina S.N., JETP Letters. - 1991. - Vol.17. - No. 24. - P.63-66. - (in Russian).
15. 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.
16. Methods of Computer Optics (Secondary Edition) / edited by V.A. Soifer - Moscow: Fizmatlit, 2003. - 688 p. - (in Russian).
17. 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. Skida-nov,V. A. Soifer, and S. N. Khonina // Optics and Spectros-copy. - 2000. - Vol. 89. - No.2. - P.347-352. - (In Russian)
18. Soifer V.A., Doskolovich L.L., Kazanskiy N.L. Multifocal diffractive elements // Optical Engineering. - 1994. - Vol.33, № 11. - P.3610-3615
19. Soifer V.A., Kazanskiy N.L., Kharitonov S.I. Synthesis of a Binary DOE Focusing into an Arbitrary Curve Using the Electromagnetic Approximation // Computer Optics, 1996, № 16, pp. 22-27 (In Russian)
20. Diffractive Computer Optics / edited by V.A. Soifer -Moscow: Fizmatlit, 2007. - 736 p. - (in Russian).
21. Meister, M. Novel approaches to direct search algorithms for the design of diffractive optical elements / M. Meister, R.J. Winfield // Optics Communications 2002. - Vol. 203.
- P. 39-49.
22. Kirkpatrick S. Optimization by Simulated Annealing / S. Kirkpatrick, C. D. Gelatt and M. P. Vecchi // Science -1983. - Vol. 220. - No. 4598. - P. 671-680.
23. Holland, J. H. Adaptation in Natural and Artificial Systems / University of Michigan Press, M.I., 1975.
24. Pavelyev, V.S. Technology-oriented stochastic optimization of radially symmetric DOEs//Proceedings SPIE 2005.
- Vol. 5962 - P. 891-897.
25. V.V. Kotlyar Phase optical elements for generation of free-space quasimodes / V.V. Kotlyar, V.A. Soifer, S.N. Khonina // Quantum Electronics. - 1991. - Vol.11. -No.11. - P.1391-1394. - (In Russian).
OPTIMIZATION OF BINARY DOE FORMING INTENSITY DISTRIBUTION ALONG AN AXIAL FOCAL ZONE
Denis Georgievich Kachalov2 (engineer, e-mail: kachalov@ssau.ru), Vladimir Sergeevich Pavelyev1'2 (chief scientist, professor, e-mail: pavelyev@smr.ru), Svetlana Nikolaevna Khonina 12 (lead. scientist, professor, e-mail: khonina@smr.ru) 1 Image Processing Systems Institute of the RAS, 2 S.P. Korolyov Samara State Aerospace University
Abstract
This paper deals with application of the known genetic stochastic procedure to determine the optimum of the function of many variables to designing quantized DOEs forming pre-given intensity distribution along an axial focal zone. Computer simulation results are presented. Key words: diffractive optical elements, light focusing, stochastic optimization procedure.
В редакцию поступила 7.12.2009г
