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

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

Рубрика Физика и энергетика
Вид автореферат
Язык русский
Дата добавления 02.03.2018
Размер файла 536,8 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

54

Размещено на http://www.allbest.ru/

1

Автореферат

диссертации на соискание ученой степени

доктора физико-математических наук

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

Специальность 01.04.05 - Оптика

Головашкин Димитрий Львович

Самара 2007

Работа выполнена в Самарском государственном аэрокосмическом университете имени академика С.П. Королева и

Институте систем обработки изображений Российской академии наук

Научный консультант: член-корреспондент РАН В.А. Сойфер

Официальные оппоненты:

доктор физико-математических наук, член-корреспондент РАН,

Крыжановский Б.В.

доктор физико-математических наук, профессор Ивахник В.В.

доктор физико-математических наук, профессор Степанов С.А.

Ведущая организация: Самарский филиал Физического Института имени П.Н. Лебедева РАН.

Защита состоится 2 ноября 2007 г. в 10 часов на заседании диссертационного совета Д 212.215.01 в Самарском государственном аэрокосмическом университете имени академика С.П. Королева по адресу: 443086, Самара, Московское шоссе, 34.

С диссертацией можно ознакомиться в библиотеке университета.

Автореферат разослан " " г.

Ученый секретарь

диссертационного совета к. т. н., профессор

В.Г. Шахов

Общая характеристика работы

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

Актуальность темы.

Дифракционная компьютерная оптика развивается более 25 лет, начиная с основополагающих работ А.М. Прохорова, И.Н. Сисакяна и В.А. Сойфера. За прошедшие годы решены фундаментальные задачи фокусировки и селекции мод лазерного излучения, формирования бездифракционных пучков и т.п. Созданные ДОЭ нашли широкое применение в лазерных технологических установках, оптических приборах и устройствах хранения и поиска информации. Количество публикаций по данной тематике отечественных и зарубежных авторов к настоящему времени весьма велико и продолжает бурно расти. В силу этого тематика работы является актуальной в широком смысле.

Необходимо отметить тенденции к миниатюризации ДОЭ и их интеграции с другими оптическими или электронными компонентами различных устройств. Кроме улучшения массогабаритных характеристик, миниатюризация позволяет использовать для изготовления ДОЭ материалы, применение которых при производстве элементов рефракционной оптики либо невозможно, либо связано со значительными техническими трудностями и финансовыми затратами. Миниатюризация приводит к уменьшению линейных размеров зон Френеля на ДОЭ и выдвигает соответствующее требование субволнового разрешения технологических установок, применяемых для их создания. При этом методы расчета, основанные на геометрическом и скалярном приближениях, становятся неадекватными, что приводит к постановке задачи решения уравнений Максвелла в векторной форме и обусловливает актуальность основного направления данной работы - расчета дифракции лазерного излучения в рамках строгой электромагнитной теории. В частности, дифракции на алмазных поликристаллических пленках (впервые примененных в 1999 году Коновым В.И., Кононенко В.В., Павельевым В.С. и Сойфером В.А. для формирования ДОЭ) с нанесенным субволновым микрорельефом.

Интеграция элементов микрооптики открывает широкие возможности для создания гибридных оптических структур, сочетающих достоинства ДОЭ и элементов традиционной оптики. Примером тому может служить формирование фокусатора в прямоугольник на торце кварцевого волновода (Prasciolu M., Cojoc D., Cabrini S и др., 2003 г.). С развитием исследований в данном направлении связано нанесение светоделительной дифракционной решетки на торец галогенидного ИК-волновода, произведенное в 2005 году Волковым А.В., Павельевым В.С. и Моисеевым О.Ю. Совмещение ДОЭ и волновода позволяет снизить потери на френелевское отражение (показатель преломления материала волновода n=2,15) и избежать юстировки объединенной оптической системы. Следовательно, актуально изучение распространения света через дифракционную решетку, сформированную на торце галогенидного ИК-волновода - также требующее применения строгой теории дифракции в силу субволнового профиля решетки.

Развитие микротехнологии позволяет создавать ДОЭ технологическими автоматами с субволновым разрешением. Однако, несовершенства технологий приводят к отклонению получаемых характеристик дифракционного микрорельефа от расчетных. Учетом влияния технологических погрешностей изготовления на работу дифракционных микролинз и бинарных решеток из оптического стекла (n=1,5) занимались Волков А.В., Скиданов Р.В., 2000; Скиданов Р.В., Хонина С.Н., 2004; Досколович Л.Л., Тявин Е.В. 2005 г, оставаясь, однако, в рамках геометрической и скалярной оптики.

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

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

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

Появившись в середине прошлого века (G. Cron, 1944 год), разностный метод решения уравнений Максвелла прошел несколько стадий развития. Ранее всего (S. K. Yee 1966 год) были записаны разностные уравнения, обладающие высокими порядками аппроксимации исходной дифференциальной задачи по времени и пространству. Однако, предложенные S. K. Yee явные схемы характеризуются условной устойчивостью, в силу чего актуальна задача разработки безусловно устойчивых разностных схем для уравнений Максвелла.

Десять лет назад J. - P. Berenger удовлетворительно решил проблему численного описания поглощения излучения, покидающего границы вычислительной области. Однако расположение поглощающих слоев в вычислительной области и традиционные граничные условия препятствует успешной векторизации вычислений по методу. Актуальна разработка новой компоновки поглощающих слоев в вычислительной области.

Задача моделирования работы источника падающей волны, поставленная еще Yee в 1966 г., решается с разной точностью вплоть до настоящего времени. Первый способ задания падающей волны, позволяющий ограничить вычислительную область исследуемым объектом и его ближайшей окрестностью, был сформулирован в 1975 году (Taflove A., Brodwin M.). Более точный метод опубликован в 1980 году (Taflove A.) в рамках TF/SF методики (Total-Field/Scattering-Field technique). С повышением точности данного подхода в области, заключенной в оболочку из однородной среды, связана работа 1999 года (D. W. Prather, S. Shi), авторы которой предпочли задавать падающую волну численно, вместо аналитической формы, как это ранее предлагалось. Однако, при моделировании работы дифракционных оптических элементов конструирование такой оболочки приводит к многократному росту вычислительной сложности алгоритма, так как в расчетную область приходится помещать весь оптический элемент. В силу чего актуальна адаптация методики задания падающей волны для решения задач дифракционной оптики.

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

Прием декомпозиции хорошо известен применительно к уравнению теплопроводности (Самарский А.А., Вабищевич П.Н., 2003), в силу чего актуальна его адаптация для уравнений Максвелла.

Параллельные алгоритмы решения явных разностных схем Yee (Perlik A. T., 1989) широко используются, однако алгоритмы вычислений по неявным схемам разработаны лишь в самом общем виде (монографии Голуба Дж., Ван Лоуна Ч., 1999 и Ортеги Д. М., 1991), без учета специфики задачи дифракции на ДОЭ. Решение этой задачи представляется актуальным.

Целью работы является решение задачи дифракции лазерного излучения на элементах микрооптики с применением разностных схем для уравнений Максвелла, декомпозиции сеточной области, реализацией параллельных вычислений и исследование на этой основе алмазных поликристаллических ДОЭ и дифракционных решеток на торце волновода.

В соответствии с поставленной целью определены основные задачи диссертации:

1. Создание математических методов и алгоритмов для проведения вычислительных экспериментов, обеспечивающих исследование дифракции лазерного излучения на элементах микрооптики; в том числе: внутри оптического элемента; на ДОЭ с высокой числовой апертурой; на ДОЭ, сформированных на диэлектрических и проводящих материалах.

2. Анализ погрешностей формирования дифракционных решеток на торце ИК-волновода и их влияния на оптические характеристики решетки.

3. Исследование прохождения света через микрорельеф ДОЭ на алмазных пленках с технологическими погрешностями изготовления и разработка методов уменьшения влияния таких погрешностей.

4. Сокращение вычислительных затрат при расчете дифракции электромагнитных полей по разностным схемам посредством декомпозиции сеточной области.

5. Синтез параллельных алгоритмов для реализации вычислений по неявным разностным схемам при решении задачи дифракции на оптическом микрорельефе.

Научная новизна работы.

дифракция лазерное излучение оптический

1. Для расчета дифракции на субволновом оптическом микрорельефе предложены неявные разностные схемы для уравнений Максвелла, отличающиеся от явных схем Yee способом выражения сеточных функций.

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

3. Разработаны методы задания падающего поля в рамках технологий TF/SF и TFF, отличающиеся от ранее известных расположением границы разделения полей и более точным учетом отраженной от ДОЭ волны, что позволяет от 3 до 5 раз снизить погрешность вычислительного эксперимента.

4. Численно исследовано влияние технологических погрешностей формирования дифракционной решетки на торце галогенидного ИК-волновода (перетрав/недотрав канавки профиля, образование клина травления, прогиб профиля при штамповке) на распределение энергии прошедшей волны по дифракционным порядкам. Показано, что наибольшее искажение в дифракционную картину (в несколько раз превышающее влияние остальных технологических погрешностей) вносит клин травления, уводя из рабочих порядков до четверти энергии падающей волны. Даны рекомендации о применении иной технологии для формирования матрицы штамповки, что позволяет снизить на 25% потери энергии в рабочих порядках.

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

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

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

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

На защиту выносятся:

неявные разностные схемы для уравнений Максвелла, основанные на методах расщепления и переменных направлений;

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

два метода задания падающей на ДОЭ волны: основанные на использовании вспомогательных одномерных задач и разделении полей;

результаты исследований технологических погрешностей нанесения дифракционного микрорельефа на торец галогенидного ИК-волновода и анализ их влияния на распределение интенсивностей дифракционных порядков за торцом;

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

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

метод декомпозиции сеточной области при разностном решении уравнений Максвелла, основанный на учете локально устоявшегося поля и принципе суперпозиции;

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

Практическая ценность работы.

Практическую значимость имеет разработанный разностный метод решения уравнений Максвелла, позволяющий моделировать распространение электромагнитного излучения через ДОЭ; в том числе через элементы с высокой числовой апертурой, внутри ДОЭ (изготовленных из различных материалов).

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

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

1) штамповки ДОЭ на торце галогенидного ИК-волокна,

2) абляции алмазной поликристаллической пластины. Предложенные рекомендации позволят повысить дифракционную эффективность элементов на торце галогенидного ИК-волокна и алмазных оптических элементов.

Апробация работы. Основные положения и результаты диссертационной работы докладывались на следующих конференциях и семинарах: VIII Всероссийское совещание по проблемам построения сеток для решения задач математической физики (Абрау-Дюрсо, 1998г.); IX Межвузовская конференция "Математическое моделирование и краевые задачи" (г. Самара, 1999 г.); Международная конференция "Diffractive Optics" (г. Йена, Германия, 1999 г.); VI Международная конференция "Электродинамика и техника СВЧ и КВЧ" (г. Самара, 1999г.); Международная школа для молодых ученых и студентов по оптике, лазерной физике и биофизике (г. Саратов, 1999 г.); Байкальская молодежная научная школа по фундаментальной физике (г. Иркутск, 1999 г.); Всероссийская научная конференция "Высокопроизводительные вычисления и их приложения" (г. Черноголовка, 2000 г.); Международная конференция "Математическое моделирование - 2001" (г. Самара, 2001 г.); Международный симпозиум "Optical Science and Technology" (США, г. Сан-Диего, 2003); Международная конференция “Science and computing" (г. Москва, 2003г.); Конференция "Photon Management" международного симпозиума “Photonics Europe” (Франция, г. Страсбург, 2004 г.); Четвертый международный научно-практического семинар и Всероссийская молодежная школа "Высокопроизводительные параллельные вычисления на кластерных системах" (г. Самара, 2004 г.); II Международная научная конференция “Сетевые компьютерные технологии” (Беларусь, г. Минск, 2005 года); Конференция "Photon Management II" международного симпозиума “Photonics Europe” (Франция, г. Страсбург, 2006 г.); Международная конференция “Laser Beam Shaping VII" (США, г. Сан-Диего, 2006); Международная конференция ICO “Topical Meeting on Optoinformatics/Information Photonics" (г. Санкт-Петербург, 2006 г.); Международный Российско-Китайский семинар по дифракционной оптике (КНР, г. Сиань, 2007 г.); Всероссийская научная конференция с международным участием “Математическое моделирование и краевые задачи” (г. Самара, 2007 г.).

Связь с государственными и международными программами

Работы по теме диссертаций проводились при поддержке грантов Президента Российской Федерации (МК-2568.2005.9) и РФФИ, Российской академии наук (программы Президиума РАН “Математическое моделирование и интеллектуальные системы”, “Поддержка молодых ученых”), программы фундаментальных научных исследований Отделения информационных технологий и вычислительных систем РАН “Новые физические и структурные решения в инфотелекоммуникациях”, федеральной целевой научно-технической программы Федерального агентства по науке и инновациям “Исследования по приоритетным направлениям науки и техники на 2002-2006 годы”, ведомственной научной программы Федерального агентства по образованию “Развитие научного потенциала высшей школы" и российско-американской программы “Фундаментальные исследования и высшее образование" (“BRHE”).

Публикации

Материалы диссертации опубликованы в 3-х монографиях и 47 основных научных работах, список которых приведен в конце автореферата.

Структура и объем работы. Диссертация состоит из Введения, четырех Глав, Заключения, Приложения, списка использованных источников из 214 наименований. Содержит 276 страниц, 78 рисунков, 41 таблицу.

Содержание работы

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

В первой Главе приводится математический аппарат, используемый для решения задачи дифракции на оптическом микрорельефе.

Для исследования ДОЭ с субволновыми размерами характеристических неоднородностей, при которых приближения геометрической и скалярной оптики теряют адекватность, предлагается применение строгой теории дифракции, основанной непосредственно на решении уравнений Максвелла. Эти уравнения и составляют основу представленной в первой Главе математической модели дифракции на ДОЭ. Для записи уравнений Максвелла выбирается декартова система координат, в которой наиболее естественно описывается квантованный ступенчатый микрорельеф ДОЭ.

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

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

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

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

В случае исследования ограниченного оптического элемента указанный прием позволяет произвести объединение восьми слоев (рис.1а), поглощающих по направлениям Y (слои 1,2,3), - Y (слои 5,6,7) и Z (слои 3,4,5), - Z (слои 1,7,8) в три (рис.1б): слои B и С поглощающие излучение, распространяющееся по направлению Y и - Y, слои A, B - по Z и - Z. Такие слои не позволяют рассеянному полю вернуться в область вычислительного эксперимента и исказить дифракционную картину.

а б

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

Переход от традиционной вычислительной области, изображенной на рис.1а, к области на рис.1б производится смещением начала координат в центр слоя B. Такое смещение возможно в силу цикличности граничных условий и сопровождается переносом электрических стенок в середину слоев A,B,C. Располагаясь в толще поглощающих слоев стенки не будут оказывать влияние на формирование дифракционной картины и могут быть удалены (рис.1б).

Достоинством предложенного перехода признается универсальность новой вычислительной области. С ее помощью (в отличие от области, представленной на рис.1. а) возможно моделирование не только ограниченных, но и бесконечных (пусть по направлению Y) периодических ДОЭ. Для этого обнуляются соответствующие электрические проводимости в слоях C и B, после чего поглощение по Y и - Y становится невозможным, а циклические условия ограничивают один период изучаемого элемента. При моделировании симметричного ограниченного оптического элемента ось симметрии располагается в середине слоев C и B (вдоль направления Z); проводимости в верхней части этих слоев обнуляются, а на ось накладывается магнитная стенка. Аналогично для симметричного неограниченного оптического элемента. Ось симметрии располагается в середине слоев C и B (вдоль направления Z); обнуляются проводимости в верхней и нижней части указанных слоев, а на ось накладывается магнитная стенка.

Подчеркивается возможность более удачной векторизации вычислений на новой области. Так, если минимальный размер вектора, состоящего из значений сеточных функций, ранее задавался толщиной необъединенных поглощающих слоев, то после объединения эта величина возросла вдвое. Как известно, эффективность векторизации возрастает с увеличением длины векторов. Для области, размерами 44 ( - длина волны) при толщине необъединенных слоев в , переход к области с объединенными слоями позволил на четверть снизить длительность вычислений, производившихся в среде Matlab 7.0 на процессоре AMD Opteron 244.

За начальное условие предлагается принимать отсутствие поля в вычислительной области в начальный момент времени. Не рекомендуется задавать поле аналитически в тех случаях, когда его формирование в части вычислительной области до начала вычислений возможно (например, дифракцию падающей волны в свободном пространстве до оптического элемента можно рассчитать заранее). Это приведет с снижению точности вычислительного эксперимента в силу неполного согласования аналитического и численного решений.

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

Построение неявных схем производится методами покоординатного расщепления и переменных направлений. Первый подход связан с представлением о раздельном распространении электромагнитной волны на каждом временном слое сначала по направлению Y, затем по Z. В результате три дифференциальных уравнения (случай распространения H-волны) аппроксимируются четырьмя сеточными: двумя неявными, выражающими электрическое поле, и двумя явными - для выражения магнитного через найденное электрическое.

Метод переменных направлений также подразумевает двухэтапный переход к следующему временному слою, при этом распространение на каждом этапе моделируется по обоим направлениям - Y и Z. Первое разностное уравнение для выражения электрического поля является неявным по Y и явным по Z. Второе наоборот, явным по Y и неявным по Z. Уравнения для нахождения магнитного поля остаются явными. Добавление явной части в уравнения, из которых выражается электрическое поле, приводит к увеличению порядка аппроксимации по времени.

Особое внимание уделяется методам задания падающей волны. Различаются два подхода к решению этой задачи. Первый (TF/SF методика) связывается с разделением электромагнитного поля (рис.2а) на результирующее (в подобласти оптического элемента) и рассеянное (в расположении поглощающих слоев). Такой подход реализуется записью на границе раздела полей разностных уравнений, содержащих напряженности падающего электрического и магнитного полей. Аналитическое задание падающей волны приводит к увеличению погрешности вычислительного эксперимента. Для случая расположения оптического элемента в вычислительной области целиком, известна более совершенная методика численного задания (D. W. Prather, S. Shi, 1999 г.) падающей волны, основанная на решении вспомогательной одномерной задачи распространения T-волны в вакууме. Напряженности полей, рассчитанные во вспомогательной задаче, используются при моделировании дифракции на ДОЭ в качестве напряженностей падающей волны.

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

Таблица 1. Погрешности моделирования распространения электромагнитного поля при реализации различных методик формирования падающей волны в рамках подхода, предусматривающего разделение полей. Параметры дискретизации сеточной области (Q,Qt) соответствуют количеству отсчетов, приходящихся на одну длину волны () по пространству и числу узлов во временной области, разбивающих интервал, за который плоская электромагнитная волна проходит в вакууме расстояние равное .

(Q,Qt)

способы задания падающей волны

TF/SF, аналитический

TF/SF, численный

TF/RF, численный

10/20

13,9267

12,3299

4,387

20/40

4,3193

3,8751

0,9745

50/100

1,4352

1,3202

0,1515

100/200

1,2614

1, 1947

0,0383

Применение двух вспомогательных задач позволяет снизить погрешность вычислительного эксперимента (Табл.1, третий столбец) при моделировании распространения волны через границу раздела диэлектрик (n=1,5) /вакуум.

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

Наложение сеточной области с циклическими краевыми условиями (рис.1б) позволяет при задании падающей волны не заключать оптический элемент в оболочку, образованную границей разделения результирующего и рассеянного полей, а разделять поля на результирующее и отраженное (методика TF/RF).

а б

Рис.2. Разделение полей в методике TF/SF (случай а) и TF/RF (случай б). Граница раздела нанесена штрих-пунктирной линией.

При этом граница такого раздела представляет собой отрезок, параллельный оси Y и расположенный до оптического элемента (рис.2б).

Реализация методики TF/RF связана с решением одной вспомогательной задачи (моделирование распространения T-волны в однородном диэлектрике) и характеризуется повышением точности задания падающей волны (Табл.1, четвертый столбец). Кроме того, снижается вычислительная сложность задания в силу сокращения протяженности границы раздела полей. При исследовании распространения излучения в волноводах, вытянутых вдоль оси Z, упомянутое снижение будет особенно заметным.

Продемонстрирована невозможность строгого следования концепции TF/SF. При аналитическом задании падающей волны неизбежно появление всплеска модуля комплексной амплитуды напряженности поля за правой границей разделения (рис.3).

Наблюдаемый всплеск объясняется временной задержкой, необходимой излученной на левой границе разделения волне для достижения правой границы. В течение этого времени волна, излученная на правой границе, распространяется в направлении Z.

Численная постановка падающей волны в методике TF/SF позволяет исправить это несоответствие. Однако, переход от исследования распространения излучения в свободном пространстве (случая, полезного для верификации разностного решения) к дифракции на ДОЭ, опять связан с нарушением принципа разделения полей. В этом случае причиной нарушения служит временная задержка, вызванная уменьшением скорости распространения волны в материале ДОЭ по сравнению со свободным пространством.

Рис.3. Распределение модуля комплексной амплитуды электрического поля () в области результирующего (0z2) и рассеянного (z<0 и z>2) полей при распространении T-волны в вакууме.

Второй подход к заданию падающей волны не характеризуется разделением полей: во всей вычислительной области присутствует только результирующее поле (TFF методика). При переходе к новому слою по времени предлагалось (Taflove A., Brodwin M., 1975 г.) определять результирующее поле в области задания падающей волны как сумму напряженностей электрических компонент такой волны (задаваемую аналитически) и результирующего поля (определенного на предыдущем временном слое). Результаты моделирования распространения T-волны в свободном пространстве по предложенному правилу представлены во втором столбце Табл.2.

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

Таблица 2. Погрешности моделирования распространения электромагнитного поля при реализации двух вариантов методики TFF.

(Q,Qt)

способы задания падающей волны

TFF, первый вариант

TFF, “прозрачный" источник

10/20

8,4609

1,3948

20/40

2,2177

0,1023

50/100

0,3741

0,0369

100/200

0,1019

0,0104

Тогда отраженную волну можно найти как разность результирующих полей основной и вспомогательной задач. Результирующее поле на новом временном слое в основной задаче задается сложением известной (аналитически выраженной) падающей волны и найденной отраженной. Предложенный метод получил название “прозрачного” источника, при котором, в отличие от “жесткого" источника (hard source, Taflove A., Hagness S., 1995 г.), отраженная от ДОЭ волна, возвращаясь к области задания падающего поля, не испытывает дополнительного переотражения.

Показано, что применение “прозрачного” источника при задании падающей волны не сопряжено со снижением точности моделирования: погрешности вычислительных экспериментов совпали с погрешностями, вызванными заменой производных разностными отношениями и наложением поглощающих слоев (Табл.2, третий столбец). Таким образом, в первоначальной формулировке методики TFF (Taflove A., Brodwin M., 1975 г.) за отраженную волну фактически принималась результирующая, что приводило к появлению дополнительной погрешности (Табл.2, второй и третий столбцы).

Представлено сравнение изложенных подходов TF/RF и TFF (на примере лучших модификаций) к заданию падающей волны. Показано, что оба не вносят дополнительной погрешности в результаты моделирования (за счет численного задания падающей либо отраженной волны), однако второй характеризуется меньшей вычислительной сложностью в силу того, что оперирует только электрическим полем, не нуждаясь в определении магнитного.

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

Во второй Главе методом вычислительного эксперимента проведено исследование влияния технологических погрешностей формирования дифракционной решетки на торце галогенидного ИК-волновода на распределение энергии прошедшей волны по дифракционным порядкам.

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

Исследование дифракции на бинарной идеально проводящей решетке с периодом 2,5, шириной канавки 1,25 и глубиной профиля при =1 мкм проводилось модовым и разностным методами. Результаты сравнения решений задачи дифракции представлены в Табл.3 и свидетельствуют о допустимости применения изложенного в первой Главе математического аппарата при исследовании идеально проводящих решеток.

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

интенсивности

дифракционных порядков

отраженной волны

метод моделирования

модовый

разностный

J0

0,6753

0,6824

J1=J-1

0,1093

0,1064

J2=J-2

0,053

0,05205

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

Для верификации на примере диэлектрической решетки привлекался метод связанных волн. Тестовая решетка характеризовалась теми же геометрическими размерами, что и предыдущая, показатель преломления материала задавался как n=1,5. Приведенные в Табл.4 результаты моделирования подтверждают состоятельность разностного решения и в этом случае.

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

интенсивности

дифракционных порядков

прошедшей волны

метод моделирования

связанных волн

разностный

I0

0,0552

0,05577

I1= I-1

0,3498

0,36845

I2= I-2

0,0756

0,07484

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

Рис.4 Зависимость отношения интенсивности отраженной волны к падающей (p) от коэффициента заполнения q для теории эффективных сред нулевого (o), второго () порядка и разностного решения (*).

Третьим тестовым примером служила антиотражающая бинарная алмазная решетка с периодом л/4 и глубиной профиля , где n=2,4, л=10,6 мкм. Отношение ширины ступеньки бинарной решетки к периоду (fill factor, коэффициент заполнения) варьировался от 0,1 до 0,9. Традиционно, для исследования антиотражающих свойств подобных структур применяется теория эффективных сред, известная в двух вариантах, как теория нулевого и первого порядков. На рис.4 представлены результаты моделирования дифракции света на предложенной решетке. Хорошее соответствие данных позволяет говорить о разностном методе решения уравнений Максвелла, как об инструменте моделирования работы антиотражающих структур.

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

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

Рассмотрены технологические погрешности формирования дифракционной светоделительной решетки на торце галогенидного ИК-волновода. Выделены два типа погрешностей: локальные, присущие каждому периоду решетки, и глобальная погрешность, характеризующая профиль в целом.

К “локальным” отнесены ошибки формирования, заметные при исследовании одного периода решетки (рис.5): в получении заданной высоты микрорельефа (недотрав или перетрав) и наличие клинов травления по обе стороны ступеньки формируемой бинарной решетки.

Рис.5 Профилограмма участка решетки, полученная с помощью интерферометра “Zygo Newview 5000. Ось абсцисс направлена поперек штрихов решетки, по оси ординат откладывается высота профиля. Пунктирная линия соответствует нулевой отметке высоты микрорельефа.

Отклонения высоты профиля решетки от расчетной связываются с ошибками в определении длительности травления матрицы, с которой штамповался микрорельеф. При недостаточном времени травления высота оказывалась ниже расчетной (недотрав), при избыточном - выше (перетрав). Если расчетная высота h профиля составляет 4,6 мкм, то высота формируемого рельефа варьируется от 4,0 до 5,2 мкм.

Образование клина травления при производстве матрицы объясняется изотропией химических процессов и аморфностью подложки. Острый угол между плоскостью торца и подъемом клина (характеризующий клин травления) для представленной технологии варьировался в интервале 36-440.

Глобальной технологической погрешностью, выявленной при изучении профилограммы всей решетки (рис.6) является прогиб профиля.

Рис.6 Профилограмма решетки, полученная с помощью интерферометра “Zygo Newview 5000. Ось абсцисс направлена поперек штрихов решетки, по оси ординат откладывается высота профиля.

Прогиб профиля решетки объясняется различием в структуре материалов сердечника и оболочки волновода и деформацией матрицы при штамповке. Величина прогиба в экспериментах сильно зависела от температурного режима штамповки и оказываемого на торец давления, варьируясь от 0,3 до 2 мкм в центре волновода.

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

При изучении влияния недотрава и перетрава ступеньки бинарной решетки, рассматривались результаты моделирования, представленные в Табл 5.

Таблица 5. Распределение энергии прошедшей волны по дифракционным порядкам в зависимости от высоты профиля прямоугольной решетки.

Высота профиля

h, мкм

интенсивности дифракционных порядков

I0

I1=I-1

I2=I-2

4,0

0,0716

0,3244

0,0111

4,2

0,046

0,3349

0,0126

4,4

0,0275

0,3417

0,0143

4,6

0,0163

0,345

0,0163

4,8

0,0125

0,3448

0,0184

5,0

0,016

0,341

0,0208

5,2

0,0267

0,3339

0,0234

Для решетки, лишенной указанных недостатков (h0=4,6), в первый и минус первый прошедшие порядки направляется суммарно 69% энергии падающего пучка, в нулевой, второй и минус второй порядки - по 1,63%. Отмечаются хорошие светоделительные показатели идеально изготовленной решетки. Исследуя случаи недотрава (h<h0) и перетрава (h>h0), делается вывод о меньших (по сравнению с недотравом) потерях, вызванных перетравом, характеризующимся не только меньшей суммарной энергией первого и минус первого порядков (65% для h=4,0 по сравнению с 67% для h=5,2), но и существенно большим количеством энергии, направляющейся в нулевой порядок (7,16% по сравнению с 2,67%). Таким образом, хотя суммарная энергия рабочих порядков (первого и минус первого) снижается несущественно (на 4% от общей энергии падающей волны), при максимальном недотраве становится заметным влияние нулевого порядка (возрастает с 1,63% до 7,16%).

К рассмотренной технологической погрешности добавляется клин травления с характеризующим его углом, равным 400. Результаты моделирования, представленные в Табл.6, свидетельствуют о существенном влиянии клина травления на распределение интенсивностей.

Так, в случае точно выдержанной высоты профиля (h0=4,6) добавление клина уводит из рабочих порядков четверть энергии падающей волны, а в нулевом порядке появляются лишние десять процентов такой энергии. Максимальный недотрав не ухудшает распределение энергии по первому и минус первому порядкам, но доводит долю энергии в нулевом порядке до 17,28%, что уже не позволяет говорить о такой решетке как о светоделительной. При предельном перетраве доля энергии, ушедшей в нулевой порядок, не увеличивается, однако, сокращается до 37% суммарная энергия рабочих порядков (при h0 она составляла 45%), что опять не позволяет считать светоделительное свойство решетки удовлетворительным.

Таблица 6. Распределение энергии прошедшей волны по дифракционным порядкам в зависимости от высоты профиля прямоугольной решетки с клином травления.

высота

профиля

h, мкм

интенсивности дифракционных порядков

I0

I1=I-1

I2=I-2

4,0

0,1728

0,2324

0,02796

4,2

0,1453

0,236

0,02827

4,4

0,13

0,2317

0,02896

4,6

0,1196

0,2241

0,02874

4,8

0,1141

0,2134

0,02778

5,0

0,1142

0,2

0,02601

5,2

0,1177

0,1858

0,02433

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

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

Таблица 7. Распределение энергии прошедшей волны по дифракционным порядкам в зависимости от глубины прогиба профиля решетки, изготовленной с идеальной матрицы.

глубина

прогиба

G

интенсивности дифракционных порядков

I0

I1=I-1

I2=I-2

0,2

0,0142

0,3487

0,0164

0,4

0,0142

0,347

0,0163

0,6

0,0141

0,3435

0,0162

0,8

0,014

0,34

0,016

1,0

0,0138

0,3328

0,0156

1,2

0,0136

0,326

0,0152

1,4

0,0134

0,3183

0,0147

1,6

0,0131

0,3082

0,0142

1,8

0,0127

0,2993

0,0136

2,0

0,0123

0,2859

0,0128

Даже в случае предельной величины прогиба (g=2,0) суммарные энергетические потери в первом и минус первом порядках составили всего 12%, что вдвое меньше потерь от появления клина травления при точно выдержанной высоте профиля и отсутствии прогиба. Решетка продолжает демонстрировать хорошие светоделительные характеристики: интенсивности нулевого, второго и минус второго порядков не меняются (Табл.7) при выбранных значениях g.

Добавление клина травления и варьирование высоты профиля при максимальном прогибе (Табл.8) существенно сказываются на распределении энергии по порядкам по сравнению с предыдущей серией экспериментов (Табл.7).

Таблица 8. Распределение энергии прошедшей волны по дифракционным порядкам в зависимости от высоты профиля прямоугольной решетки с клином травления и предельной глубиной прогиба.

высота

профиля

H

интенсивности дифракционных порядков

I0

I1=I-1

I2=I-2

4,0

0,1295

0, 1921

0,0211

4,2

0,1116

0, 1912

0,0206

4,4

0,1025

0,1844

0,0204

4,6

0,0974

0,175

0,0197

4,8

0,0969

0,1631

0,0186

5,0

0,0996

0,1506

0,0176

5,2

0,1048

0,1384

0,0169

Однако, оценивая влияние именно прогиба и сравнивая с результатами аналогичной серии вычислительных экспериментов с решеткой без прогиба (Табл.6), видим, что из рабочих порядков потеряно дополнительно не более десяти процентов энергии (независимо от высоты профиля). При этом поведение значений интенсивностей нулевого порядка (Табл.8) по сравнению с результатами из Табл.6 не ухудшилось.

Решетка, сочетающая глобальную и локальные погрешности изготовления, обладает низкими светоделительными характеристиками. Особенно это относится к случаю предельного перетрава, когда прошедшая энергия почти равномерно распределяется между нулевым, первым и минус первым порядками (10,48%, 13,84% и 13,84%). Более того, от такого профиля внутрь волновода отражается 58,46% падающей энергии.

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

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

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

Задача дифракции H-волны на бесконечном диэлектрическом цилиндре круглого сечения имеет аналитическое решение в виде разложения в ряд по цилиндрическим функциям: Бесселя внутри и Ханкеля вне цилиндра. Радиус тестового цилиндра задавался равным r=/2 (=1 мкм), показатель преломления материала выбирался как n=1,5. Отличие разностного решения от аналитического искалось на отрезке длиной 2 (с серединой в центре кругового сечения цилиндра и расположенном параллельно направлению падения плоской однородной электромагнитной волны) и характеризовалось (Табл.9) равномерной погрешностью разностного решения

и величиной (1)

, (2)

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

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

Таблица 9. Зависимости погрешностей вычислительных экспериментов при исследовании дифракции на диэлектрическом цилиндре круглого сечения от дискретизации сеточной области (Q,Qt) для “прозрачного” источника (колонки А) и аналитического задания падающей волны в методе TF/SF (колонки В).

(Q,Qt)


Подобные документы

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.