Статьи Соросовского Образовательного журнала в текстовом формате
Излагаются основы создания дифракционных оптических элементов (ДОЭ) с широкими функциональными возможностями. Обсуждаются пути и методы получения зонированных пластинок со сложным профилем зон. Анализируются уравнения границ и профиля зон. В качестве примеров рассматриваются плоские сферические и цилиндрические линзы, плоская призма, отражающие фокусаторы.
КОМПЬЮТЕРНАЯ ОПТИКАЧасть 1. Дифракционные оптические элементы
В. А. СОЙФЕР
Самарский государственный аэрокосмический университет им. С.П. Королева
ВВЕДЕНИЕ
Предметом рассмотрения этой статьи является новая область знаний, терминология в которой еще до конца не установилась. Разные авторы используют различные названия: "компьютерная оптика", "оптика, синтезированная компьютером", "бинарная оптика", "плоская оптика", "дифракционная оптика" и т.д. В любом случае речь идет об оптических элементах, различные отличительные особенности которых раскрываются перечисленными названиями.
Внешне такой оптический элемент представляет собой пропускающую или отражающую пластинку с тонким фазовым микрорельефом, рассчитанным в рамках теории дифракции. Первым представителем этого класса оптических элементов является дифракционная решетка, созданная более 200 лет тому назад, задолго до появления компьютеров. Следующим по хронологии представителем указанного класса оптических элементов является зонная пластинка. Эти дифракционные оптические элементы (ДОЭ) имели бинарное амплитудное или фазовое пропускание. Если дифракционные решетки нашли широкое применение в приборостроении, то зонные пластины в основном использовались в учебном лабораторном практикуме по оптике для иллюстрации положений теории дифракции.
Настоящую революцию в создании ДОЭ произвело применение компьютеров. Причем вначале, в 70-х годах, компьютеры стали применять для синтеза голограмм математически заданных объектов. При этом была решена ключевая проблема кодирования - запись на физическую среду комплекснозначных функций и получение соответствующих амплитудно-фазовых транспарантов. Это послужило толчком к компьютерному синтезу бинарных амплитудно-фазовых и фазовых пространственных фильтров, для записи которых на физическую среду применялись графопостроители, а затем и фотопостроители. В результате прогресса в области микроэлектроники и лазерной техники, который привел к появлению в конце 80-х годов прецизионных фотопостроителей и электронных литографов, практически реализовалась задача создания дифракционных оптических элементов со сложным профилем зон. Возникло широкое поле деятельности для физиков, оптиков, специалистов в областях прикладной математики и электроники, технологов, специалистов по автоматизации. Появились ДОЭ с уникальными характеристиками, недостижимыми в рамках традиционной оптики, например фокусаторы лазерного излучения.
ФУНКЦИОНАЛЬНЫЕ ВОЗМОЖНОСТИ ЗОНИРОВАННЫХ ДИФРАКЦИОННЫХ ОПТИЧЕСКИХ ЭЛЕМЕНТОВ
Амплитудная одномерная (1D ) дифракционная решетка представляет собой плоский транспарант, на котором чередуются светлые и темные параллельные полоски (штрихи и щели). Полоски имеют одинаковую ширину. Если осветить такую решетку монохроматическим пучком света с длиной волны l, падающим нормально к плоскости решетки, то в результате дифракции света на периодической структуре щелей будет образовано множество пучков света, выходящих под различными углами a p и соответствующих различным порядкам дифракции. Угол a p зависит от периода решетки и в предположении малости углов определяется по формуле
T - период решетки. Интенсивность света убывает с ростом абсолютной величины p. Как следует из теории дифракционных решеток, при равной ширине штриха и щели характер убывания определяется формулой
Из-за того что значительная часть света поглощается амплитудным транспарантом, дифракционная эффективность в первом порядке не превышает 10%. У одномерной фазовой дифракционной решетки этот показатель приблизительно в 4 раза выше.
Функция пропускания в данном случае является чисто фазовой, периодической с периодом T, аргумент которой скачкообразно меняется на p. Функционально фазовая и амплитудная дифракционные решетки одинаковы.
Естественным развитием одномерных бинарных амплитудных и фазовых дифракционных решеток являются соответствующие двумерные (2D ) радиально-симметричные решетки (зонные пластинки). На рис. 1 показан центральный фрагмент зонной пластинки Рэлея-Сорэ.
Поясним понятие зоны. Зона - это ограниченная область ДОЭ, на которой функция пропускания света претерпевает однократное изменение от минимального до максимального значения. Для амплитудной дифракционной решетки зона представляет собой сочетание темной и светлой полос в пределах одного периода решетки. Граница зон при этом прямая линия. Для зонной пластинки Рэлея-Сорэ зона представляет собой совокупность темного и светлого колец переменной толщины (рис. 1). Граница зон при этом окружность переменного диаметра.
Радиусы окружности меняются пропорционально квадратным корням из последовательных целых чисел p:
где f - фокусное расстояние.
Фрагмент 2D фазовой зонной пластинки изображен на рис. 2.
Зонная пластинка выполняет функцию фокусировки света и ведет себя как линза с множеством фокусов f- 2 , f-1 , f0 , f1 , f2 , соответствующих различным порядкам дифракции и расположенных на оси распространения света:
Интенсивность света убывает с ростом номера порядка по закону
Совершенствование базовых ДОЭ, какими являются 1D дифракционные решетки, возможно за счет изменения профиля зон. При этом получаются так называемые решетки с блеском. Например, можно рассмотреть амплитудную 1D дифракционную решетку, функция пропускания которой меняется по закону При дифракции на такой решетке теоретически помимо нулевого образуются только ?1 порядки дифракции. Практически, с учетом погрешностей изготовления ДОЭ, дело обстоит гораздо сложнее.
Для двумерных (2D ) дифракционных решеток естественными являются построение зон, отличных от периодических щеле- и кольцеобразных, и выявление новых функциональных возможностей ДОЭ. Несмотря на кажущуюся очевидность, с момента создания зонной пластинки этот вопрос не поднимался в течение многих лет из-за отсутствия практической возможности реализации зонных пластинок с варьируемым характером зон. Такая возможность появилась в связи с созданием компьютеров и технических средств машинной графики, в частности графо- и фотопостроителей.
А.В. Ломан и Д.П. Парис в 1967 году предложили 2D дифракционные решетки с варьируемыми зонами Френеля, получаемые на основе муаровых эффектов в результате суперпозиции пары синтезированных масок. Синтезированные маски содержат периодические по одной из координат ( например, по x) бинарные изображения, описываемые специально подобранной алгебраической функцией j(x, y). Варьируя вид и параметры функции j(x, y), а также параметры относительно сдвига Dx и поворота масок можно получить множество амплитудных 2D дифракционных решеток, каждой из которых будет соответствовать своя картина дифракции. К сожалению, описанная процедура содержит регулярный метод синтеза только для цилиндрических и сферических зон Френеля. Результат фокусировки представляет собой прямую фокальную линию. Имеется также возможность построения конических зонных пластинок путем оптической регистрации результата физической суперпозиции цилиндрической и плоской волн. Результат фокусировки представляет собой наклонную прямую линию в плоскости, параллельной зонной пластинке.
Рассмотренные примеры демонстрируют возможность получения зонных пластинок с различной формой зон. В то же время представляет интерес построение 2D бинарных образов зон Френеля для более сложных случаев, например для зонной пластины, фокусирующей в продольный или поперечный отрезок заданной длины кольцо или какую-либо другую геометрическую фигуру. ДОЭ такого рода получили название фокусаторов. Рассмотрим для примера фокусатор в кольцо. Соответствующий образ зон Френеля решетки может быть получен путем комбинации 1D дифракционной решетки и зонной пластинки. Возьмем достаточно узкий сегмент 1D дифракционной решетки, который ведет себя так же, как целая дифракционная решетка, то есть отклоняет входной монохроматический пучок на определенный угол в плоскости (рассматривается первый порядок дифракции).
Вращая сегмент вокруг центра получим ДОЭ, зоны Френеля которого отображаются системой равноотстоящих концентрических черно-белых колец одинаковой ширины. Назовем его дифракционным аксиконом. Нетрудно сообразить, что такой ДОЭ будет отклонять входной пучок на определенный телесный угол в пространстве. Для того чтобы собрать дифрагированные пучки в фокальной плоскости, применим зонную пластинку, причем дифракционный аксикон и зонная пластинка могут быть совмещены в одной плоскости и записаны на одной подложке. При этом в плоскости ДОЭ образуется система концентрических, неравномерно отстоящих колец переменной ширины. Таким образом, мы синтезировали амплитудный ДОЭ, фокусирующий в кольцо.
Следуя описанной методике можно построить различные фокусирующие ДОЭ. Например, записывая на одной подложке две скрещенные цилиндрические линзы с разными фокусными расстояниями можно получить ДОЭ, фокусирующий в отрезок, лежащий в плоскости, перпендикулярной оси распространения (рис. 3).
По-видимому, нет особой необходимости говорить, что и фокусатор в кольцо, и фокусатор в поперечный отрезок могут быть реализованы как фазовые ДОЭ. Достичь этого в простейшем случае можно путем фотографического отбеливания соответствующих амплитудных масок (фотошаблонов).
По существу в двух только что рассмотренных примерах речь идет о фазовой модуляции, при которой одна бинарная функция, описывающая дифракционную решетку, умножается на другую. Однако процесс модуляции параметров дифракционной решетки может осуществляться и при более широком классе модулирующих функций. В частности, осуществляя модуляцию радиусов колец зонной пластинки можно построить фокусатор в соосный отрезок определенной длины, который фактически будет представлять собой зонную пластинку с протяженной продольной аберрацией. Идея получения соответствующего радиально-симметричного ДОЭ состоит в том, чтобы его периферийная часть соответствовала зонной пластинке с фокусным расстоянием F1 , а центральная часть - зонной пластинке с расстоянием F2 > F1 . Между центральной и периферийной частями должны быть записаны кольца, соответствующие зонным пластинкам, фокусное расстояние которых уменьшается от F2 к F1 . В результате будет получен фокусатор в соосный отрезок длины l = F2 - F1 . На рис. 4 показан процесс фокусировки.
Общие методы получения модулированных дифракционных решеток, обладающих заданными функциональными свойствами, требуют решения обратных задач теории дифракции. Заметим, что решение обратной задачи дифракции осуществляется относительно фазовой функции ДОЭ, который выполняет требуемое преобразование светового пучка. Эта фазовая функция должна быть записана на оптическую среду в виде зонированной структуры с фазовым микрорельефом, в результате чего и получается ДОЭ.
ГРАНИЦЫ ЗОН И ФАЗОВЫЕ ФУНКЦИИ ДИФРАКЦИОННЫХ ОПТИЧЕСКИХ ЭЛЕМЕНТОВ
Приведение функции к интервалу
Под фазовой функцией j(u, u) оптического элемента будем понимать функцию, описывающую зависимость от пространственных координат величины набега фазы, создаваемой оптическим элементом в точке (u, u).
Иногда границы зон ДОЭ могут быть получены исходя из фазовых функций соответствующих традиционных (рефракционных) оптических элементов.
Традиционный оптический элемент характеризуется гладкой фазовой функцией где - поперечные координаты точки, лежащей в плоскости оптического элемента. При освещении оптического элемента пучком с комплексной амплитудой непосредственно за оптическим элементом образуется поле
Толщина традиционного оптического элемента H может составлять тысячи длин волн. При этом значения фазы j лежат в интервале от 0 до тысяч единиц 2p. В то же время изменения фазы на рассмотренных ранее ДОЭ лежат в диапазоне от 0 до 2p.
C учетом того, что комплексная экспонента в выражении (6) является периодической функцией с периодом 2p, фаза может быть приведена к интервалу [0, 2p). Приведение к интервалу можно выполнить по формуле
F = mod2p j,
где
mod2p j = j - 2pj
при
j " 2p # j # ( j + 1)2p, j = 0, ?1, ? 2, _
Из практических соображений может оказаться целесообразным приводить фазу к интервалу [0, 2pm) по формулам
F = mod2pm j,
где обычно m ~ 1-102. Эквивалентность исходной гладкой фазы j и приведенной к интервалу кусочной фазы F выражается соотношением
G ╞ exp (iF) = exp (ij).
Плоская сферическая линза
Уравнение фазовой функции сферической линзы в параксиальном приближении имеет вид
где k = 2p / l - волновое число, f - фокусное расстояние, D - диаметр линзы,
Если материал линзы имеет коэффициент преломления n, то максимальная высота рельефа составляет
и имеет порядок длины волны. Высота микрорельефа определяется по формуле
Радиусы зон Френеля можно найти из соотношения
j(rj) = - 2pj,
откуда следует
Число полных зон j0 на линзе определяется из условия и удовлетворяет соотношению
где ] " [означают целую часть числа с округлением в меньшую сторону.
Ширина зон плоской сферической линзы
Dj = rj - rj - 1 , j = 1, 2, _, j,
является переменной и уменьшается к периферии линзы. Характерным параметром является ширина D самой узкой (в данном случае последней) периферийной зоны, определяющей требования к технологическому оборудованию. В более общем случае для ДОЭ с фазовой функцией j(u, u) можно также определить ширину наиболее узкой зоны микрорельефа с максимальной фазовой высотой 2pm. Используя в пределах этой зоны линейную аппроксимацию фазовой функции получим выражение
где max берется вдоль поверхности ДОЭ,
Сравним между собой уравнения радиусов колец зонной пластинки Рэлея-Сорэ (3) и сферической плоской линзы (12). Видим, что радиусы зон на линзе в раз больше радиусов колец на зонной пластинке. Соответственно при заданном диаметре число колец на зонной пластинке в 2 раза больше, то есть каждому кольцу сферической плоской линзы соответствуют два кольца зонной пластины: темное и светлое. Плоская линза имеет непрерывный (в пределах зоны) фазовый рельеф, а фазовая зонная пластинка имеет аппроксимирующий его бинарный фазовый рельеф.
Если подходить к созданию ДОЭ с общих позиций дискретной аппроксимации непрерывной фазовой функции, можно продолжить процедуру ступенчатой аппроксимации в пределах зоны Френеля исходной непрерывной фазовой функции. На рис. 5 показан пример равномерного квантования фазы плоской сферической линзы Френеля на четыре уровня.
Для практического получения амплитудной маски с N градациями яркости можно использовать N - 1 бинарных амплитудных масок, осуществляя последовательный фотопроцесс. Существует также возможность уменьшить число бинарных масок до n = log2 N.
Нетрудно рассчитать радиусы колец, соответствующих градациям фазы в различных зонах. Знание их необходимо для расчета бинарных амплитудных масок, однако они не несут такого ясного физического смысла, как в случае зонной пластинки. Отметим также, что уравнение колец зонной пластинки Рэлея-Сорэ невозможно получить с помощью равномерной дискретизации фазовой функции плоской сферической линзы.
Плоская цилиндрическая линза
Рассмотрим цилиндрическую линзу, описываемую фазовой функцией
Приведение фазовой функции к интервалу [0, 2p) аналогично проведенному выше для сферической линзы. Границы зон в данном случае - прямые линии, а расстояния между ними определяются формулой (12) при замене rj на uj .
Отражательная зонная пластина
Рассмотрим отражательную зонную пластину, предназначенную для работы под углом a с освещающим пучком и осуществляющую фокусировку в точку F.
При нормальном падении максимальная высота рельефа отражательной зонной пластины должна составить l /2 для того, чтобы набег фазы при прямом и обратном ходе составил 2p, при наклонном падении
Уравнения границ зон отражательной зонной пластины при нормальном падении в точности соответствуют уравнениям границ зон пропускающей пластинки (высота зон различается в 2/(n - 1) раза).
Уравнение границ зон при наклонном падении под углом можно получить путем простых рассуждений. Наклонное падение пучка эквивалентно введению поворотного зеркала. Спроецировав кольца зонной пластинки Рэлея-Сорэ на это зеркало мы получим уравнение краев зон отражательной зонной пластинки при наклонном падении под углом a. Понятно, что это будут эллипсы с коэффициентом растяжения по оси u, равным cos a.
Высота рельефа отражательного ДОЭ, фокусирующего в точку под углом a,
где j(u, u) определяется формулой (10).
При выводе (16) использовано приближение плоской волны для падающего и отраженного пучков.
Плоская призма
Рассмотрим призму с углом a. Призма обеспечивает фазовый сдвиг, линейно зависящий от координаты, и характеризуется фазовой функцией
где q = (n - 1) tg a.
Приводя фазовую функцию j(u) к интервалу [0, 2p), получим 1D дифракционную решетку с блеском (рис. 6).
Максимальная высота рельефа
и уравнение высоты микрорельефа имеет вид
Различным углам отклонения b соответствуют различные периоды решетки
где sin b = (n - 1) tg a.
Бинарная 1D амплитудная дифракционная решетка получается при замене линейно меняющейся фазовой функции в пределах одного периода на двоичную функцию, принимающую значения ? p /2.
Комбинированные и сегментированные ДОЭ
Выше уже говорилось о возможности получения комбинированных зонных пластинок, представляющих собой суперпозицию двух исходных зонных пластинок.
Если два расположенных соосно друг за другом оптических элемента имеют фазовые функции j1 и j2 соответственно, то они могут быть заменены комбинированным оптическим элементом с фазовой функцией j = j1 + j2 .
В качестве примера рассмотрим две скрещенные цилиндрические линзы с различными фокусными расстояниями
Фазовая функция комбинированного ДОЭ имеет вид
Границы краев зон являются эллиптическими. Такой ДОЭ осуществляет фокусировку в отрезок, перпендикулярный оптической оси.
Плоскость сегментированного ДОЭ в пределах апертуры разбивается на два или большее число сегментов, а в каждом из сегментов записывается зонная пластинка, соответствующая определенной фазовой функции. Например, таким образом можно получить фокусатор в букву R, фотошаблон которого приведен на рис. 7.
ЛИТЕРАТУРА
1. Борн М., Вольф Э. Основы оптики. М.: Наука, 1973. 720 с.
* * *
Виктор Александрович Сойфер, профессор, ректор Самарского государственного аэрокосмического университета, специалист в области обработки изображений и компьютерной оптики. Лауреат Государственной премии России в области науки и техники 1992 года за разработку лазерных технологий и их внедрение при создании новой авиационно-космической техники. Член международного оптического общества SPIE. Автор 310 научных работ, в том числе трех монографий, 50 работ в зарубежных изданиях.