Формирование интерференционных картин поверхностных электромагнитных волн с изменяемым периодом с помощью дифракционных решеток

Автор: Безус Евгений Анатольевич, Досколович Леонид Леонидович, Кадомин Иван Иванович, Казанский Николай Львович, Civera Pierluigi, Pizzi Marco

Журнал: Компьютерная оптика @computer-optics

Рубрика: Дифракционная оптика, оптические технологии

Статья в выпуске: 3 т.32, 2008 года.

Бесплатный доступ

Рассмотрено формирование интерференционных картин поверхностных электромагнитных волн с помощью диэлектрической дифракционной решетки с металлическим слоем. Моделирование в рамках электромагнитной теории показывает возможность получения контрастных интерференционных картин с периодом в несколько раз меньшим периода дифракционной решетки. При этом интенсивность поля в интерференционных максимумах в десятки раз превышает интенсивность падающей волны. Рассмотрены способы управления периодом интерференционной картины за счет изменения длины волны и угла падения.

Дифракция, дифракционная решетка, фотолитография, интерференционная картина, поверхностная электромагнитная волна

Короткий адрес: https://sciup.org/14058823

IDR: 14058823

Текст научной статьи Формирование интерференционных картин поверхностных электромагнитных волн с изменяемым периодом с помощью дифракционных решеток

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

В [1] показана возможность получения интерференционной картины затухающих волн, соответствующих -1,+1 затухающим порядкам субволновой дифракционной решетки (ДР). Интерференционная картина в [1] обладает высоким контрастом и высокой интенсивностью в 3-4 раза большей, чем интенсивность падающей волны. Период интерференционной картины в [1] в два раза меньше периода используемой ДР. В [2, 3] рассмотрен аналогичный метод, основанный на использовании интерференции ПЭВ, возникающих на поверхности перфорированной металлической пленки. В [4] предложен метод, основанный на интерференции затухающих волн, получаемых при полном внутреннем отражении. В [5, 6] для формирования интерференционной картины ПЭВ используется ДР, расположенная над металлической пленкой. В настоящей работе, как и в [6], для получения интерференционной картины ПЭВ используется диэлектрическая ДР с металлическим слоем в области подложки. В отличие от [5, 6], для возбуждения ПЭВ предлагается использовать высшие затухающие порядки дифракции. Приведенные расчеты показывают высокое качество формируемых интерференционных картин. Использование высших затухающих порядков позволяет формировать высокочастотные интерференционные картины с помощью низкочастотной ДР, имеющей на порядок меньший период. При этом коэффициент усиления поля на порядок больше, чем в [5].

Для данной структуры рассмотрены способы управления частотой интерференционной картины за счет изменения длины волны и угла падения волны.

Формирование интерференционной картины ПЭВ

Исследуемая структура состоит из бинарной диэлектрической ДР и металлической пленки, расположенной под решеткой (рис. 1). ДР предназначена для возбуждения на нижней границе металлической пленки двух встречных ПЭВ, которые формируют интерференционную картину.

Рис. 1. Геометрия структуры (один период) и формируемая интерференционная картина

При нормальном падении TM-поляризованной волны на структуру константы распространения дифракционных порядков имеют вид kxn = 2nnd , ky = 0,                          (1)

где d – период решетки. Условие возбуждения ПЭВ дифракционными порядками с номерами ± n на нижней границе раздела металлический пленки имеет вид

V k x. + k y = Re ( k spp ( ^ 0 ) ) ,                          (2)

где kspp (Х 0 ) = k 0 V S m S II K£m + S II ) - константа распространения ПЭВ, к 0 = 2 п / Х 0, Х 0 - длина волны, s II - диэлектрическая проницаемость материала под пленкой, s m - диэлектрическая проницаемость металлической пленки.

Согласно (2), при периоде решетки d = 2п n/ Re (kspp) (3) происходит возбуждение ПЭВ. При этом период формируемой интерференционной картины ПЭВ dint = d /2n , (4)

в 2 n раз меньше периода ДР. На рис. 1 схематично показана интерференционная картина при n =5.

В [5] ПЭВ возбуждались высшими незатухающими порядками, сформированными ДР в области над металлической пленкой. В этом случае диэлектрическая проницаемость материала решетки s gr должна быть больше, чем диэлектрическая проницаемость материала s II под пленкой. В данной работе рассматривается случай, когда порядки дифракции, используемые для возбуждения ПЭВ, являются затухающими. При этом значение s gr может быть равным s II или меньше его.

Расчет интерференционной картины ПЭВ проводился при следующих параметрах: Х 0 = 550 нм , s I = 1, s II = 2,56, s m = -12,922 + 0,44727 i . Значение s II соответствует фоторезисту, а s m - диэлектрической проницаемости серебра для выбранной длины волны. Диэлектрическая проницаемость материала решетки s gr была также выбрана равной 2,56. Период решетки d = 1539,1 нм был рассчитан из (3) при n =5. В этом случае период интерференционной картины d int = 154 нм в 10 раз меньше периода ДР. Значения остальных геометрических параметров структуры w = 0,5 d , h gr = 435,4 нм, h = 0 , h m = 65 нм были определены с помощью оптимизационной процедуры. Целевой функцией являлась мера близости расчетной интерференционной картины к «идеальной» интерференционной картине, формируемой при интерференции двух ПЭВ. Для расчета интерференционной картины использовался модовый метод (rigorous coupled wave analysis) в формулировке работ [7-9]. График нормированной интенсивности поля, формируемого непосредственно под металлическим слоем при указанных параметрах, представлен на рис. 2. Рис. 2 показывает формирование интерференционной картины с расчетным периодом d int = 154 нм .

Отметим, что период dint не в только в 10 раз меньше периода ДР, но и в 3,57 раз меньше длины волны. Коэффициент усиления поля, показывающий значение интенсивности в максимумах интерферен- ции относительно интенсивность падающей волны, близок к 50. Контраст полученной интерференционной картины равен 0,701. Таким образом, предложенная структура позволяет формировать интерференционную картину высокого качества.

Рис. 2. Интерференционная картина

Управление частотой интерференционной картины

Рассмотрим возможность управления частотой интерференционной картины за счет изменения длины волны и угла падения.

Рассмотрим использование различных длин волн для формирования интерференционных картин различного периода. Пусть период d в (3) определен из условия возбуждения ПЭВ порядками с номерами ±n при длине волны Х0. Константа распространения ПЭВ к „ = к „ (X) зависит от длины волны. Поэто-spp spp му возможно возбуждение ПЭВ при другой длине волны X * Х0 с использованием дифракционных порядков ±m, m * n . Длины волн, которые будут возбуждать ПЭВ порядками ±m, могут быть найдены из уравнения

2 п m / d = k spp (X) .                                   (5)

В частности, для рассмотренного выше случая X 0 = 550 нм, n = 5 , ПЭВ будут возбуждаться ±4 порядками при X =659 нм и порядками с номерами ±3 при Х = 852 нм. Соответствующие интерференционные картины будут иметь периоды 192 нм и 257 нм. Геометрические параметры структуры w = 0,37 d , h gr = 1000 нм, h , = 0 , hm = 69,8 нм , предназначенной для формирования трех различных интерференционных картин при трех указанных длинах волн, были рассчитаны с использованием оптимизационной процедуры. Как и в предыдущем случае, целевой функцией являлась мера близости расчетных интерференционных картин к «идеальным» интерференционным картинам, формируемым при интерференции двух соответствующих ПЭВ.

Расчетные графики нормированной интенсивности интерференционных картин приведены на рис.

3а-в. Значения контраста и коэффициентов усиления поля составляют (0,87;76), (0,78;37), (0,71;25), соответственно.

|£"№Jj

Рис. 3. Интерференционные картины, соответствующие длинам волн 852, 659 и 550 нм

где 9 - угол падения. Согласно (6), условие возбуждения ПЭВ (2) может выполняться при различных сочетаниях угла падения и номера порядка дифракции i .

Рассмотрим пример. На рис. 5 приведены графики модуля проекции волнового вектора kn (9) = ^k2n + ky (9) от угла падения для n = 3 (штрих-пунктирная линия), n = 4 (точечная линия) и n = 5 (пунктирная линия). Графики получены при еI = 2,31, £gr = 9, значения параметров Х0,d,£m,£II совпадают с предыдущими случаями. Условия kn (9) = Re(kspp) выполняются в точках 70,53° (n=3), 45° (n=4) и 0° (n=5). Это означает, что при указанных углах ПЭВ будут возбуждаться ±3, ±4 и ±5 порядками дифракции. Периоды соответствующих интерференционных картин составляют 257нм, 192нм и 154нм, соответственно. Геометрические параметры структуры w = 0,534d , hgr = 273 нм, hl = 95,9 нм, hm = 60 нм были найдены в процессе оптимизации из условия максимизации качества формируемых интерференционных картин.

0 10 20 30 40 50 60 70 80 90

Рис. 5. Зависимость k n^ + k y от угла падения

Расчетные графики нормированной интенсивности на нижней границе металлического слоя показаны на рис. 6 а-в . Значения контраста и коэффициента усиления составляют (0, 67; 25), (0, 79; 35), (0, 88; 41), соответственно.

Заключение

Рассмотрено формирование интерференционных картин ПЭВ с помощью диэлектрической ДР с металлическим слоем. Показана возможность формирования интерференционных картин с высоким контрастом и высокой интенсивностью при использовании высших затухающих порядков дифракции. Рассмотрены способы управления периодом интерференционной картины. Результаты исследования показывают возможность формирования интерференционных картин с высоким контрастом и интенсивностью при различных длинах волн и углах падения.

Работа выполнена при поддержке грантов РФФИ № 07-07-97601-р_офи, 07-01-96602-р_поволжье_а, 07-07-91580-АСП_а, гранта Президента РФ № НШ-3086.2008.9, фонда «Фундаментальные исследования и высшее образование» (RUXO-014-SA-06) и Фонда содействия отечественной науке.

Статья научная