Расчет спектральных дифракционных решеток Текст научной статьи по специальности «Физика»

ЧИСЛЕННЫЕ МЕТОДЫ КОМПЬЮТЕРНОЙ ОПТИКИ

РАСЧЕТ СПЕКТРАЛЬНЫХ ДИФРАКЦИОННЫХ РЕШЕТОК

Л.Л. Досколович Институт систем обработки изображений РАН, г. Самара

ДОЭ создаются, как правило, для работы с излучением одной заданной длины волны. Для большого количества прикладных задач, включающих задачи распознавания, обработки изображений, оптической связи, измерения температурных характеристик объектов, цветной печати и фотографии, требуются ДОЭ, позволяющие работать с излучением различных длин волн. Для работы с излучением различных длин волн известны фазовые дифракционные решетки (color separation gratings), позволяющие разделить 3 плоских пучка с длинами волн N N

А), А+1 = АоN ’ А = АN ’ ^

N + p N - p

где N, p - взаимно простые числа, по -1, 0 и +1 дифракционным порядкам [1-4]. В дальнейшем будем называть такие решетки спектральными. Спектральная решетка имеет на периоде N ступенек равной ширины. Высота ступенек определяется по формуле

di =-

Я0

-modN (а • i), i = 0, N -1,

(2)

(по - О

где а - целое число, определяемое по формуле

а = ( ■ 2 + 1)/(( + р), г є 2 . (3)

Поскольку числа Ми р взаимно простые, то уравнение (3) всегда имеет решение в целых числах по теореме о наибольшем общем делителе. Пренебрегая дисперсией материала решетки, определим связь фазового набега ф с высотой рельефа решетки для длины волны Я в виде:

р=Т( -1).dl

i = 0, N -1.

(4)

Согласно (4), решетка (2) для плоских пучков с длинами волн (1) формирует следующие фазовые набеги по модулю 2п:

2п ■ і ----------

Фо = 0, ф±1 = +—, і = о, N -1. (5)

Согласно (5), для длины волны Х0 фазовый набег равен нулю и, следовательно, данная спектральная компонента направляется в нулевой порядок. Для длин волн Х±1 фазовые набеги (5) соответствуют фазовым функциям призм, квантованных по N уровням. Это обеспечивает разделение длин волн Я±1 по порядкам ±1. Вследствие квантования доля энергии спектральных компонент Х±ь направляемая в порядки ±1, описывается функцией

Е(ы )= Шпс2 (ж/ N). (6)

Рассмотрим расчет спектральных дифракционных решеток, обобщающих известные решетки (2)

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

1 N

F(P1,...,PN )= II l=-1 j=1

\.2п .

exp| г 1 • j

N

N + p • l

-exp| i—n— P

(7)

где (pj = 2n ■ modN (a1 j), j = 0, N -1 (8)

соответствуют фазовым набегам решетки (4) для длины волны А0. Учитывая вышеприведенную связь функции невязки (7) с решеткой (2), предлагается искать решетку для разделения трех произвольных длин волн

А0, А+1 = А0c+1, А-1 = А0с-1 (9)

по порядкам ±1 из условия минимума функции невязки

1 N

F(P1,...,PN )= II

l=-1J=1

exp i ~Пі •j |-

- exp

—Pi

(10)

= 1

Естественным технологическим ограничением при синтезе решетки является ограничение на максимальную высоту профиля. Отметим, что высота профиля решетки (2) превышает в N раз высоту рельефа обычной решетки для работы с излучением одной длины волны Х0. Поэтому при больших значениях N в (1) или при общем соотношении длин волн целесообразно искать профиль спектральной решетки из условия минимума функции невязки (10) при ограничении ре [0,2пК], ] = 0, N , где значение К определяет максимальную высоту профиля. При этом значения р несложно найти простым перебором из условия минимума невязок

I sin2

p=-1

ж

N

j

p

2c,

(11)

где р}- е [0,2пК], ] = 1, N, с0 = 1. В качестве примера была рассчитана решетка для разделения длин волн Х0=1,55 мкм, Х+1=1,31 мкм, Х-1=1,625 мкм, использующихся в системах телекоммуникаций. Сле-

2

2

c

0

v cl

дует отметить, что приведенные длины волн не удовлетворяют соотношению (1);

А0

А

-+1

1,55 ^ 6 Ад_

1,31 ~ 5 , А-1

1,55 20

1,625 21

Расчет решетки проводился при К=16, N=8 и Р- е12п'1-,1- =1,2,... , К}. Рассчитанные высоты ступенек профиля равны

(,...,ё8)= /,55 ч(0,1,7,8,8,14,15,16)мкм . (12) («0 -1)

Выбранный вид значений р- = 2п • , 1- < К

дает для длины волны А0 фазовый набег по модулю 2п равный нулю. Поэтому пучок с длиной волны А0 направляется в нулевой порядок. Для оценки энергетической эффективности решетки были рассчитаны интенсивности ±1 порядков для длин волн А±1.

I (А±1 ) =

1 I i

J expl-----p(x)± 2nix

dx

(13)

где р(х)- фазовая функция решетки при А=А0. Согласно расчетам, решетка (12) направляет 100% энергии спектральной компоненты А0=1,55 мкм в 0-й порядок, 86,7% энергии спектральной компоненты Х+1=1,31 мкм в 1-й порядок и 80,1% энергии спектральной компоненты Х-1=1,625 мкм в -1-й порядок.

В качестве следующего примера была рассчитана решетка для разделения длин волн А0=1,55мкм, А+1=А0(21/23)=1,4152мкм, Х-1=Х0(21/19)=1,7132мкм,

удовлетворяющих соотношению (1) при N=21 и р=2. При расчете по аналитической формуле (2) решетка будет иметь 21 ступеньку на периоде при макси-

1,55

мальной высоте ступеньки а тах = 20^---------г мкм. С

(«0 - 1)

целью упрощения технологических параметров была рассчитана решетка из условия минимума невязок (11) с 8-ю ступеньками на периоде при К=8 и р- = 2п • I-, I- < К . Рассчитанные высоты ступенек

профиля равны

(,...,а8)= /,55 ч(0,1,3,4,5,7,8,8)мкм . (14)

(«0 -1)

Согласно результатам расчета интенсивностей ±1 порядков (13), решетка (14) направляет 100% энергии спектральной компоненты А0=1,55 мкм в 0-й порядок и по 87,5% энергии спектральных компонент А±1 в порядки 1-й и -1-й. Полученная решетка, по сравнению с аналитически рассчитанной решеткой (2), теряет 12,5% энергии для спектральных компонент А±ъ однако имеет в 2,5 раза меньшую высоту и более простую структуру периода. Приведенный пример показывает, что при N>8,9 для снижения технологических требований решетки для разделения длин волн (1) также следует рассчитывать из условия минимума невязок (11).

Приведенный подход легко обобщается на случай разделения более трех длин волн. При этом для

разделения длин волн А0

рядкам n0,...,nL значения pj

К = V,

n = 1, L по по-[0,2пК ] ищутся из

L 2

£ sin2

p=0

п

---npj

N p

Pj

2ct

= 1, j = 1, N .

(15)

В качестве примера была рассчитана решетка для разделения пяти длин волн А0=1,55 мкм, А-2= А0(23/19), Х-1=Я.0(21/19), Х+1=Хо(21/23), А+2=А0(21/25) по порядкам 0, -2, -1, +1, +2. Расчет решетки проводился при р- = 2п • I-, I- < К, К=10 и N=16. Рассчитанные высоты ступенек профиля равны 1,55

-10,1,1,2,3,3,4,

(16)

(,...,d16)= (1,55 )(0,1,1,2,3,3,4,

(n0 -1) 5,5,6,7,7,8,9,9,10)мкм.

Для решетки (16) доли энергии спектральных компонент А±1,А±2, направляемых в порядки -1, +1,-2,+2 составляют 96,3%, 96,3%, 80,4% и 85,9%, соответственно.

Использованный метод расчета спектральных решеток может быть также использован для расчета сложных дифракционных оптических элементов (ДОЭ), выполняющих заданные преобразования волновых фронтов для излучения различных длин волн. Будем считать, что микрорельеф ДОЭ определен на двумерной сетке отсчетов (, -) I, I = 1, N . Предположим, что для длин волн

А0, Ап =А0сп, п = 1, Ь ДОЭ должен сформировать пучки с фазовыми функциями р0(/,-),...,рЬ0,-). В

этом случае значения рг

= -тМп0 -^dj

А

опреде-

ляющие высоту микрорельефа ДОЭ, определяются из условия минимума невязок

(

£ sin2

p=0

p

Aj)

Рг

2с

^ min, i, j = 1, N . (17)

p

■. [0,2пК], j = 0, N, c0 = 1.

где pj 6

Литература

1. Досколович Л.Л. Расчет дифракционных оптических элементов для разделения и фокусировки излучения различных длин волн // Компьютерная оптика. 1998. В. 18. С. 16-24.

2. Досколович Л.Л., Петрова О.И. Расчет спектральных ДОЭ // Компьютерная оптика. 1999. В. 19. С. 29-32.

3. L.L. Doskololovich, M. Repetto Design of Does for Focusing Different wavelengths // Optical Memory and Neural Network. 2000. V. 9, № 1. P. 13-23.

4. Досколович Л.Л. Расчет дифракционных оптических элементов для фокусировки различных длин волн // Автометрия. 2000. В. 3. С. 99-108.

условия минимума невязок

2

е