Численное моделирование и компьютерный дизайн оптических свойств наноструктурированных материалов тема автореферата и диссертации по физике, 01.04.05 ВАК РФ

084681731

^(а правах рукописи

Дейнега Алексей Вадимович

Численное моделирование и компьютерный дизайн оптических свойств наноструктурированных материалов

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

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

Москва - 2010

004601731

Работа выполнена в Институте водородной энергетики и плазменных технологий Российского Научного Центра „Курчатовский Институт".

Научный руководитель:

кандидат физико-математических наук, доцент Потапкин Борис Васильевич

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

доктор физико-математических наук, доцент Ёлкин Николай Николаевич (Троицкий институт инновационных и термоядерных исследований, г. Троицк)

доктор физико-математических наук, профессор Козарь Анатолий Викторович (Физический факультет МГУ имени М. В. Ломоносова, г. Москва)

Ведущая организация:

Институт спектроскопии РАН (г. Троицк)

Защита состоится „19" мая 2010 г. в 15 часов на заседании совета по защите докторских и кандидатских диссертаций Д 501.001.45 при Московском Государственном Университете имени М. В. Ломоносова по адресу: Россия, 119991, Москва, Ленинские горы, дом 1, строение 5 (19 корпус НИИ ядерной физики МГУ), ауд. 2-15.

С диссертацией можно ознакомиться в научной библиотеке НИИ ядерной физики имени Д. В. Скобельцына МГУ имени М. В. Ломоносова.

' )Г

Автореферат разослан апреля 2010 г.

Ученый секретарь совета по защите докторских и кандидатских диссертаций Д 501.001.45 при МГУ имени М. В. Ломоносова кандидат физико-математических наук

Вохник О. М.

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

Диссертация посвящена численному моделированию оптических свойств наноструктурированных материалов, а именно, в ней исследуются следующие вопросы:

• Возможность применения металлических фотонных кристаллов в качестве новых высокоэффективных источников света.

• Оптимизация размеров и формы антиотражающих нанотекстурирован-ных покрытий.

Основным используемым численным методом является метод решения уравнений Максвелла в конечных разностях (Finite-Difference Time-Domain, FDTD). Для расчетов применяется специально написанная параллельная программа, включающая в себя ряд новых оригинальных численных методов, которые также описываются в данной диссертации.

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

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

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

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

В этой работе основным используемым численным методом является метод решения уравнений Максвелла в конечных разностях (англ. - Finite-Difference Time-Domain method, FDTD) [6]. Выбор этого метода вызван присущей ему высокой параллельной эффективностью, что позволяет применять его для расчета больших задач на кластерных вычислительных системах.

Как и всякий другой разностный метод, FDTD сталкивается с проблемой неточного отображения на прямоугольную вычислительную сетку тел, обладающих произвольной формой поверхности. Решение этой проблемы особенно значимо для моделирования исследуемых в данной работе фотонных кристаллов и текстурированных поверхностей. Одним из наиболее удачных способов ее решения является метод подсеточного сглаживания [7]. Однако этот метод в его первоначальной постановке может применяться только для диэлектриков, в связи с чем актуальна его модификация, которая позволила бы расширить его на дисперсные среды.

Исследование оптических свойств наноматериалов включает в себя получение спектров прохождения и отражения от периодических структур при падении на них плоской волны. Их расчет FDTD для случая нормального падения предполагает применение периодических граничных условий. В

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

Цель работы состоит в 1) исследовании излучательных характеристик металлических фотонных кристаллов, а также возможности их практического применения в качестве источников света с высокой селективностью излучения в видимом диапазоне; 2) исследовании антиотражающих свойств текстурированных покрытий во всем диапазоне размеров составляющих текстуру рассеивателей, а также установлении оптимального размера для заданного диапазона длин волн; 3) дальнейшем развитии метода РБТБ, включающем, в частности, разработку метода подсеточного сглаживания для дисперсных сред и метода расчета наклонного падения плоской волны на периодическую структуру.

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

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

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

Получено количественное согласие результатов расчетов с результатами специально поставленного эксперимента.

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

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

5. Предложен новый подход программной реализации метода решения уравнений Максвелла в конечных разностях (РОТБ), основанный на интегральном представлении уравнений Максвелла. В рамках этого подхода написана параллельная программа на языке программирования С++. В ходе ее написания реализован ряд оригинальных решений: оптимизация использования памяти с помощью упаковки используемых в разностных уравнениях коэффициентов; увеличение производительности посредством разбиения алгоритма обновления сеток на основной цикл, поддающийся векторизации, и специализированные поправки к нему; увеличение параллельной эффективности путем балансировки доменов. Специально проведенные тесты продемонстрировали линейный характер масштабируемости программы вплоть до тысячи вычислительных ядер. Программа находится в открытом доступе в интернете (http://fdtd.kintechlab.com).

6. В рамках метода РБТО впервые предложен и реализован алгоритм подсеточного сглаживания для дисперсных сред, позволяющий существенно уменьшить величину численной ошибки, связанной с неточным отображением тел, обладающих произвольной формой поверхности, на прямоугольную вычислительную сетку. Его применение приводит к увеличению порядка точности метода ГВТБ по сравнению с обычным „лестничным" приближением.

7. Впервые предложен и реализован итерационный метод для моделирования наклонного падения плоской волны на периодическую структуру в чис-

ленном расчете РБТБ. С его помощью можно получать результаты во всем диапазоне углов от 0° до 90°, а также наблюдать реальное поведение поля во времени. Помимо этого, он лишен таких недостатков предложенных ранее методов, как неустойчивость при отвесных углах падения и увеличение требуемого размера сетки. Особенностью метода является то, что получение решения для одного угла требует нескольких отдельных итераций-экспериментов.

Положения, выносимые на защиту.

1. Геометрические параметры вольфрамовых фотонных кристаллов типа прямого опала, при которых КПД источника света на основе фотонных кристаллов выше в несколько раз по сравнению с обычными лампами накаливания.

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

3. Асимптотические зависимости величины отражения от геометрических параметров текстурированной поверхности для длинно- и коротковолнового пределов.

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

5. Параллельная программа ЕБТО, которая может использоваться для предсказательного моделирования в разнообразных оптических приложениях.

6. Новые численные методы в рамках ЕБТЮ: метод подсеточного сглаживания для дисперсных сред и итерационный метод для моделирования наклонного падения плоской волны на периодическую структуру.

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

1. Найдены геометрические параметры вольфрамовых фотонных кристаллов типа прямого опала, при которых КПД источника света на основе фотонного кристалла оказывается выше в несколько раз по сравнению с обычными лампами накаливания.

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

3. Найдена оптимальная геометрия текстурированной поверхности, при

которой достигается наименьшая величина отражения для заданного диапазона длин волн.

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

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

6. Предложены и реализованы численные методы, позволяющие существенно увеличить быстродействие и точность расчетов FDTD.

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

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

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

Апробация работы. Основные результаты диссертационной работы были доложены на научно-технических конференциях „The International Conference on Computational Science and Its Applications" (Малайзия, Kuala Lumpur, 2007), „Conference on Computational Physics" (Бразилия, Ouro Preto, 2008), „Conference on Computational Physics" (Тайвань, Kaohsiung, 2009), „Многомасштабное моделирование процессов и структур в нанотехнологиях" (Москва, 2009), „VII Курчатовская молодежная научная школа" (Москва, 2009). 52-ая научная конференция МФТИ (Москва, 2009).

Реализованная программа FDTD стала призером на конкурсе „Максимальная масштабируемость", проведенном в рамках РОСНАНОФОРУМА-2009 компаниями Intel и "Роснано".

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

и тезисы российских и международных конференций (список в конце автореферата).

Структура и объем диссертации. Диссертация состоит из введения, четырех глав и заключения, изложена на 157 страницах, включая 74 рисунков и 125 наименований цитируемой литературы.

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

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

В первой главе обсуждается метод решения уравнений Максвелла в конечных разностях (Finite-Difference Time-Domain, FDTD). Эта глава содержит в себе описание метода FDTD и изложение новых предложенных численных методов: метод подсеточного сглаживания для дисперсных сред и итерационный метод для моделирования наклонного падения плоской волны на периодическую структуру.

Метод подсеточного сглаживания предназначен для уменьшения численной ошибки, связанной с неточным отображением тел с произвольной формой поверхности на вычислительную сетку. Он основан на введении тензора обратной эффективной диэлектрической проницаемости ё-1 вблизи границы между телами [7]:

Г1 = Р<£~1 >+(1-Р) <£>-\ (1)

где Р - это матрица Рц = щщ соответствующая проектору на вектор нормали п к границе между двумя средами, а о есть усреднение по контрольному объему, окружающему выбранный узел сетки.

До сих пор метод подсеточного сглаживания применялся только для моделирования диэлектриков. Мы развили его также и для случая дисперсив-ных сред. Наша идея заключается в разбиении электрического поля на вспомогательные компоненты таким образом, что получающиеся уравнения для каждой компоненты могут быть решены во временном представлении традиционным для FDTD путем. Например, если зависимость диэлектрической проницаемости от частоты представляется в виде суммы произвольного числа членов Друде или Лоренца, эти уравнения могут быть решены методом вспомогательных дифференциальных уравнений [6].

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

Рис. 1: Слева: Зависимость фактора эффективности рассеяния С}вса от радиуса для свинцовой сферы, А = 1.5 (яа.

Справа: Угловая зависимость элемента матрицы рассеяния 5ц для свинцовой сферы радиуса г = 250 пт, Л = 1.5 /лп.

методов: 1) лестничное приближение, в котором значение диэлектрической проницаемости берется в узле сетки; 2) усреднение диэлектрической проницаемости по контуру или по объему; 3) усреднение обратной диэлектрической проницаемости; 4) УР-ЕР метод [7], который был нами реализован для дис-персивных сред; 5) собственно применение полного выражения (1).

Рассмотрим результаты расчетов для рассеяния плоской волны Л = 1.5дгп на свинцовой сфере (радиус г варьируется), рассчитанное на сетке с шагом Дг = 50 пт (Рис. 1). Отметим, что на рассматриваемой длине волны действительная часть диэлектрической проницаемости свинца отрицательна: £ « —81 + 18г. Поскольку решение уравнений Максвелла существенно различается для областей с положительной и отрицательной диэлектрической проницаемостью, методы, не учитывающие направление границы между телами и основанные на частичной аппроксимации выражения (1), показывают существенное снижение точности по сравнению новым методом, в котором примененяется тензор обратной эффективной диэлектрической проницаемости.

Далее в первой главе мы предлагаем новый итерационный метод для моделирования наклонного падения плоской волны на периодическую планарную структуру. Наш метод предполагает проведение нескольких (5-10) отдельных численных экспериментов РБТБ, которые мы в дальнейшем будем называть итерациями. Идея метода заключается в итерационном приближении к искомому решению, для которого поле Р(х, в граничных точках периодической ячейки жив любой момент времени Ь удовлетворяет условию:

Р{х,Ь) = Р{х±а^ + 5П,2), 5и = ¿а^тЯ^/с, (2)

PML.

a Sx

Рис. 2: Геометрия численного эксперимента. 1,2- генерирующие (TF/SF) границы, для которых в качестве источника берется сигнал, переносимый с точек образов 1', 2'; 3 - TF/SF граница для генерации падающей плоской волны; 4 - детекторы для записи прошедшего сигнала.

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

Геометрия численного эксперимента, соответствующего предлагаемому нами методу, представлена на Рис. 2 (в целях простоты изложения мы рассматриваем двумерный случай с а2 = а, аг = 0 и временным сдвигом St = asin0/c). Мы используем метод общего и рассеянного поля (Total Field / Scattered Field Method, TF/SF) [6] для генерации плоской волны в область общего поля (эта область заштрихована на Рис. 2). Зависящие от времени граничные условия TF/SF на границе 3 соответствуют падающей наклонно плоской волне и известны аналитически. Основная идея нового метода заключается в применении дополнительных граничных условий на границах 1 и 2, которые работают подобно границе TF/SF, где в качестве генерируемой волны используются значения полей в точках-образах на границах 1' и 2' (xi< = х\ + а и х? = Х2 — а) с соответствующим временным сдвигом. Для отрицательного временного сдвига в прошлое могут быть использованы значения полей с текущей г-ой итерации, а для положительного сдвига в будущее мы используем значения полей, записанные на предыдущей i — 1-ой итерации:

$(f2,t) = -Л), Fi(xut) = Fi.l(xv,t + 6t). (3)

Расстояние между границами 1 и 2 задается больше периода а на неко-

Рис. 3: Слева: Зависимость коэффициента прохождения однослойного фотонного кристалла, состоящего из упакованных в квадратную решетку металлических шариков (е = 4ео, а = 2, период решетки а = 1, радиус сфер 0.375) от частоты падающей волны, в = 45°. Справа: Зависимость коэффициента прохождения от угла падения, частота / = 0.6/а.

торое расстояние 5х, равное нескольким сеточным шагам, для того чтобы разделить точки-образы и дополнительную границу ТР/БР. Для записи значений Е{х\,4) и Н(х 1,4) на границе 1 используется дополнительный буфер памяти размером 6NT/At, где Т - продолжительность одной итерации, Д£ - шаг по времени, а N = Л^А^ - количество сеточных узлов, приходящихся на границу 1 для каждой из трех компонент полей Е и Н). Для

границы 2 требуется буфер размера 6Л№/Д4. Для 4 < 0 и Ь > Т поля в точках-образах полагаются равными нулю ^{хуу^) = 0. На первой итерации в качестве граничного условия на границах 1 и 2 берется аналитическое значение для падающей плоской волны или нулевой сигнал.

Далее в первой главе приводится доказательство сходимости предложенного метода. Идея доказательства заключается в представлении генерируемого на г-ой итерации на границах 1 и 2 сигнала в виде суммы Е = Ец + 5Е, где удовлетворяет условию (2), а есть сигнал-ошибка. По при-

чине линейности уравнений Максвелла и операции временного сдвига по отношению к полям /¡(ж, 4), эволюция ошибки 6Ёг при увеличении г может рассматриваться отдельно от решения Ео(хь,{). Нами доказано и численно проверено, что сигнал-ошибка затухает при увеличении числа итераций.

Метод был протестирован при моделировании однослойного фотонного кристалла, состоящего из упакованных в квадратную решетку металлических шариков (е = 4ео, с = 2). Период квадратной решетки полагался равным единице а = 1, радиус сферы 0.375, а шаг по сетке 0.05. На Рис. 3 представлено сравнение результатов для спектра прохождения, полученных с помощью описанного метода и с помощью метода ЬККИ [8]. Результаты РБТБ находятся в хорошем согласии с результатами ЬККИ, уже после пяти итераций.

Новый метод лишен таких недостатков предложенных раннее в литерату-

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

Во второй главе описывается использовавшийся нами подход для написания параллельного кода FDTD. В самом широком смысле метод FDTD включает в себя множество самых разнообразных возможностей: применение различных типов сеток, использование методов постпроцессорной обработки результатов и т. д. Для реализации всех этих возможностей требуется написание больших программ, которые обладали бы высокой производительностью, были бы легко расширяемы и обладали доступным, читаемым кодом. К сожалению, распространенные некоммерческие программы FDTD лишены этих свойств и не предназначены для широкого круга задач. Это явилось причиной разработки собственной модели FDTD, которую мы будем в дальнейшем называть EMTL (Electromagnetic Template Library). EMTL реализована в виде библиотеки шаблонов С++, которые могут быть использованы в качестве строительных компонент при создании конкретных численных приложений. Для реализации возможности проведения параллельных расчетов используется интерфейс MPI (Message Passing Interface).

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

Ключевым строительным блоком EMTL является вычислительная сетка. Сетка организует размещение данных в памяти и определяет способ интерполяции значений поля по сеточным узлам. В рамках этих требований могут быть реализованы сетки самого разного типа, включая не ортогональные и неоднородные. Проведение численного расчета с использованием сетки и внутреннее устройство самой сетки рассматриваются в EMTL как две независимые задачи. Это позволяет реализовать многие процедуры EMTL в более общем виде, так что они могут быть выполнены с помощью произвольной сетки. Такой подход обуславливает естественную легкость добавления новых типов сеток.

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

С целью измерения параллельной эффективности нашего кода при поддержке компании Intel были проведены тесты на системах Harpertown и

— Xeon, mvapich, ррп=1. диэл., сетка 1003

— Xeon, mvapich, ррп=8.ди9Л.. сетка 1003

— Xeon, OMPI.ppo^t, диэл., сетка 1003 y -+ HT. даалеюрик. сетка 2003

и—I—I—i—i—i—i—i—i—I

2 4' 8 16 32 64 128 256 512 1024 Количество ядер

EMTL (free network! -EMTL (busy network) -MEEP (free network) ■■ MEEP (busy network) -EMTL data MEEP data

10 15 20

number of processors

Рис. 4: Слева: Ускорение на единицу моделируемого объема для различных архитектур и реализаций MPI, нормировано на время исполнения на 64 ядрах. HT -Harpertown+Infiniband, HT - Nehalem+Infiniband (измерения в ИНТЕЛ), Xeon - кластер Т-тматформы, узлы 2хХеоп 5345 2.33 GHz 8 GB RAM, Infiniband.

Справа: Зависимость времени исполнения задачи от числа процессоров в серии масштабируемых тестов на кластере ФМБФ МФТИ (Athlon MP 2400, Gigabit). Начальный размер задачи (сетка 1403) масштабируется в одном измерении пропорционально числу процессоров. Количество временных шагов равно 100.

Nehalem. В этих тестах размер сетки полагался постоянным или зависел от числа процессоров таким образом, чтобы общий объем задачи сохранялся при его увеличении. Задача представляла собой численное решение уравнений Максвелла в диэлектрике и металле. Результаты тестов приводятся на Рис. 4, где показано ускорение на единицу моделируемого объема, характеризующее масштабируемость программы. Ускорение вычисляется относительно задачи, исполняемой на 64 процессорах при размере сетки 2003. Оба типа тестов показали линейный или суперлинейный характер масштабируемости приложения для большого числа вычислительных ядер. Также были проведены тесты на малом кластере, состоящем из 24 2-процессорных узлов с четырехядерными процессорами Xeon (2хХеоп 5345 2.33 GHz 8 GB RAM). Результаты тестов показывают, что характер производительности приложения зависит от архитектуры системы и реализации MPI.

Нами было проведено сравнение производительности с программой МЕЕР, разработанной в Массачусетском технологическом институте (MIT). Для этого была проведена серия масштабируемых тестов на кластере ФМБФ МФТИ (Athlon MP 2400, Gigabit). Результаты тестов показали, что производи-

тельность разработанной нами программы примерно в два раза выше, чем у МЕЕР (см. Рис. 4).

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

Вначале литературного обзора излагаются существующие методы расчета излучательной способности ФК. Затем рассматриваются имеющиеся в литературе результаты относительно излучательной способности металлических ФК. Отметим, что в настоящее время имеется лишь ограниченное число работ (см., напр., [9]), в которых исследуется излучательная способность ФК в видимом диапазоне. Вместе с тем, значительная часть результатов, касающихся влияния геометрических параметров ФК на спектр его излучения, получена для инфракрасного (ИК) диапазона в рамках исследований возможности применения ФК для управления спектром теплового излучения (см., напр. [3] и ссылки в этой работе). В этих работах рассматривались ФК разных геометрий: одномерные ФК, двумерные массивы рассеивателей и микрополостей, а также трехмерные ФК типа „поленница". Было показано, как влияют на интенсивность излучения ФК величина периода, количество слоев ФК, оптические свойства используемого материала.

Еще раз отметим, что число работ для излучения ФК в видимом диапазоне довольно ограничено. Также не было проведено какого-либо исследования возможности использования ФК в качестве источников видимого излучения. Этому исследованию посвящена основная часть третьей главы.

В качестве материала ФК мы используем вольфрам. Высокая селективность излучения вольфрама в видимом диапазоне выгодно отличает его от других металлов и вместе с высокой температурой плавления (3695fC) делает его лучшим материалом для применения в качестве источника света.

Для расчетов нами использовался метод FDTD. Для расчета излучательной способности ФК мы пользуемся способом, предполагающим расчет по-глощательной способности ФК и применение закона Кирхгофа [3], связывающего ее с излучательной способностью.

Для верификации написанной программы, реализующей метод FDTD, было проведено сравнение с экспериментом, специально проведенным в компании General Electric. Экспериментальный образец представлял собой вольфрамовый ФК монослой с периодом треугольной решетки а = 0.55 мкм. Поперечный размер рассеивающих элементов равнялся 0.2 мкм. Эксперимент проводился при комнатной температуре.

Сравнение спектра поглощения экспериментального образца с расчетами

Рис. 5: Слева: Поглощение вольфрамового ФК монослоя при нормальном падении при комнатной температуре (Т = 298К). Параметры ФК: а = 0.55мкм, поперечный размер рассеивающих элементов 0.1 мкм. 1 - экспериментальные данные; 2 - расчет идеального ФК с учетом реальной формы элементов; 3 - расчет с учетом дефектов решетки ФК и дисперсии элементов по размерам.

Справа: Поглощение ФК при нормальном падении для различных значений периода решетки а при Л = 0.07/0.45а и N = 2, Т = 2400ЛГ.

РБТБ показывает, что амплитуда пика поглощения в видимом диапазоне определяется степенью неупорядоченности решетки, а ИК-часть спектра поглощения определяется только размером элементов (Рис. 5). Последний факт связан с тем, что в ИК-диапазоне длина волны больше константы решетки, и поэтому электромагнитное излучение взаимодействует с ФК, как с эффективной средой.

Перейдем к результатам численных расчетов для температуры Т = 2400К, являющейся характерной рабочей температурой вольфрамовой нити в лампах накаливания. Мы рассматриваем ФК типа гранецентрированного кубического (ГЦК) прямого опала. Такой выбор обусловлен двумя причинами. Во-первых, прямой опал в отличие от альтернативных геометрий ФК (обратные опалы или ФК типа „поленница") обладает наилучшей термической устойчивостью при высоких температурах. Во-вторых, прямые опалы позволяют достигать малых значений фактора заполнения металлом. При этом, как нами было показано, именно при малых значениях фактора заполнения металлом достигается больший КПД светового источника.

Рассмотрим идеальный случай прямого опала, когда вольфрамовые шарики помещены в вакуум. На Рис. 5 представлена зависимость поглощатель-ной способности двухслойного ФК от длины волны падающего света Л для различных значений периода решетки а при сохранении величины отношения Я/а постоянной. Мы видим, что в спектре поглощения ФК присутствует ярко выраженный пик. Его положение связано с наименьшим вектором обратной решетки, величина которого равна 47г/\/За- Это соответствует длине волны Л = а\/3/2, которая определяет левый край пика.

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

Как показали дополнительные расчеты, при увеличении угла падения рассматриваемый пик перемещается в сторону ИК области. Увеличение радиуса шариков Д и числа слоев N в ФК приводят к увеличению амплитуды пика, что связано с увеличением доли вольфрама в ФК. С одной стороны это приводит росту излучаемой энергии в видимом диапазоне, однако при этом падает

кпд.

Нами были проведены расчеты для различных значений периода О.Змкм < а < 0.7мкм, радиуса шариков 0.02мкм < И < 0.2мкм и числа слоев 1 < N < 4. Нашей целью было нахождение таких параметров, при которых достигается максимальный КПДрс при сохранении величины светового потока, излучаемого через единичную площадку, Ррс не меньшим аналогичной величины для полубесконечного вольфрамового слоя Ршк- Наши расчеты показали, что искомыми параметрами ФК являются а = 0.45 мкм, Л = 0.07 мкм и N = 2 (фактор заполнения вольфрамом в этом случае примерно равен 2.2%). При этих значениях достигается КПДРС « 4.7%, а Ррс ~ 3.8 ■ 106лм/м2. Отметим, что, КПД полубесконечного вольфрамового слоя равен КПДЫА. = 1.6%, а световой поток, излучаемый через единичную площадку, Рык = 3.7- 10блм/м2. Таким образом, применение ФК может позволить увеличить КПД примерно в 3 раза.

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

При изготовлении ФК образца должна использоваться внешняя матрица, фиксирующая ФК элементы. Известно лишь небольшое количество материалов, которые были бы устойчивы к высоким температурам и которые можно использовать в качестве такой матрицы (в их числе: НЮг, Zr02 стабилизированный УгОз). Однако, при нагревании эти материалы начинают поглощать, а, следовательно, и излучать в ИК диапазоне, что может привести к существенному снижению КПД. Ниже мы исследуем этот вопрос на примере НЮ2.

Будем вначале считать, что поглощение отсутствует, и рассматривать ди-

Рис. 6: Слева: Поглощение ФК при нормальном падении (Г = 2400&Г) при N - 2. 1 - ФК. погружен в НЮг при комнатной температуре, а = 0.45 мкм и Л = 0.07 мкм; 2 - то же самое, а = 0.214 мкм и Л = 0.033 мкм; 3 - учитывается конечная толщина матрицы НГОг, (2 = 2Л + (1 + 1/\/3)а.

Справа: Поглощение ФК при нормальном падении (Г = 2400К) при а = 0.214 мкм, Л = 0.033 мкм, N = 2, й = 2Л + (1 + 1/\/3)а. Рассматриваются различные значения коэффициента затухания в формуле Лоренца для поглощающего НГОг (см- текст).

электрическую проницаемость НЮг при комнатной температуре (е = 4.41). Изменение диэлектрической проницаемости внешней среды от 1 до е приводит к тому, что спектральные кривые для поглощения, полученные нами для ФК в вакууме, масштабируются на величину ^/е. Это связано с уменьшением длины волны во внешней среде в \fi раз. Кроме того, конечность матрицы, в которую помещен ФК, приводит к появлению дополнительных осцилляций, вызванных интерференцией волн, многократно отраженных от границ между матрицей и вакуумом (Рис. 6).

При высоких температурах диэлектрическую проницаемость НЮг с хорошей точностью можно аппроксимировать с помощью формулы Лоренца

Дешо

Wq — 2¿W7 - ьР-

где Де = 3, ш0 = 1.2 • 1016 рад/сек, -у = 1.0 • 1016 рад/сек. Конкретные значения используемых параметров взяты нами из результатов экспериментов, проведенных в компании General Electric.

Для того, чтобы последовательно учесть влияние поглощения, мы полагаем 7 = s ■ 1016 рад/сек, и постепенно меняем параметр s от нуля до единицы (0 соответствует отсутствию поглощения, а 1 соответствует действительной диэлектрической проницаемости НЮг). Как видно из Рис. 6, при малых значениях s спектр поглощения в основном определяется вольфрамовыми шариками. При увеличении s существование матрицы НЮг вносит значительный вклад в поглощение в ИК диапазоне. При погружении ФК с найденными нами ранее оптимальными параметрами в матрицу толщины

Ф) = 1 + ,,2 п,,..° ,,2, (4)

й = 2Я + (1 + 1/\/3)а КПД снизился до 1.7%.

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

Четвертая глава посвящена изучению антиотражающих свойств тесту-рированных покрытий. Она состоит из литературного обзора и изложения новых результатов.

В имеющихся работах показано, что текстурированные покрытия обладают антиотражающими свойствами как в длинно-, так и в коротковолновом пределах. В первом случае это связано с тем, что оптические свойства таких покрытий подобны свойствам переходного слоя с плавно меняющимся эффективным коэффициентом преломления [10]. В другом предельном случае это связано с тем, что лучу требуются много переотражений внутри такой поверхности прежде чем окончательно отразиться обратно [11]. На примере кремниевых текстурированных покрытий было показано, что достижение низких значений коэффициента отражения возможно и в случае, когда характерный размер текстуры сравним с длиной волны [5]. Был эмпирически установлен ряд общих закономерностей касательно влияния размера и геометрии текстуры на ее антиотражающие свойства. Однако не было сделано какой-либо аналитической оценки характера убывания отражения при изменении геометрии рассеивателей. Отметим также, что в имеющихся работах, как правило, либо отдельно рассматривались длинноволновый или коротковолновый пределы, либо оптический диапазон для размеров рассеивателей порядка длины волны. Таким образом представляет интерес единое исследование текстурированных поверхностей во всем диапазоне длин волн.

В основной части главы мы приводим результаты исследования антиотражающих свойств текстурированных поверхностей во всем диапазоне размеров текстуры. Вначале мы оцениваем характер изменения отражения при изменении размеров текстуры в длинно- и коротковолновом пределах с помощью приближений эффективной среды и геометрической оптики соответственно. Затем мы исследуем антиотражающие свойства текстурированной поверхности во всем диапазоне размеров текстуры с помощью метода РБТО.

В конкретных численных расчетах мы рассматриваем подложку с нанесенными на нее стеклянными пирамидками с высотой ё,. Основаниями пирамидок являются треугольники, квадраты и круги (в этом случае пирамидка вырождается в конус). Пирамидки плотно упакованы в квадратную

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

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

Нами аналитически показано, что в случае полного замощения, для которого фактор заполнения (доля площади сечения пирамидок в элементарной ячейке на данной высоте) у вершин пирамидок /(0) = 0, а у оснований f(d) = 1, степень убывания R с увеличением d/X зависит от степени гладкости функции фактора заполнения f(z) в точках z = 0 и z = d (ось z направлена вдоль оси пирамидок). Так, если f(z) задавать многочленом степени (2к' - 1) таким, что /«(0) = f®(d) = 0,0 < г < /с', то R ~ (d/X)~2k'.

Нами предложен специальный профиль f(z) = С t S '(¿С, где значено

ние константы С побирается таким, чтобы f(d) = 1. В случае такого профиля коэффициент отражения при увеличении d/X уменьшается экспоненциально (см. Рис. 7). Выводы находятся в согласии с численными расчетами, проведенными с помощью приближения эффективной среды и метода FDTD, для пирамидок с различными профилями (Рис. 7).

В случае неполного замощения подложки основаниями пирамидок между этими основаниями имеется зазор (0 < /(d) < 1). Вследствие этого в точке z = d у функции эффективной диэлектрической проницаемости имеется разрыв. Отсюда следует, что при увеличении d/X отражение выходит на постоянное значение, равное коэффициенту отражения между средами со значениями диэлектрической проницаемости внешней среды и текстуры (см. Рис. 7).

Перейдем теперь к рассмотрению оптических свойств текстурированной поверхности в коротковолновом пределе Л Л. На Рис. 7 представлены результаты расчетов для плотно упакованных треугольных и квадратных пирамидок (полное замощение) и конусов (неполное замощение). Результаты, полученные с помощью FDTD находятся в хорошем согласии с результатами, полученными методом трассировки лучей.

В случае полного замощения коэффициент отражения экспоненциально

Рис. 7: Слева: Коэффициент отражения переходного слоя с оптическими свойствами, соответствующими (а) дифракционной решетке с „интегральным" профилем /(г) ~ ^еС-'О-СГ1^ (ТЕ поляризация); (Ь) плотно упакованным в квадратную решетку квадратным пирамидкам с (с) линейным к = 1 профилем и (с!) профилем к' = 3 (см. текст); (е) плотно упакованным в треугольную решетку конусам. Выражения для эффективной диэлектрической проницаемости квадратных пирамидок и конусов взяты из работ [12,13]. Результаты ГБТО (точки) соответствуют значениям параметров Л = 1, в. = 16,4 < А < оо, сеточный шаг Ах = 0.01.

Справа: Коэффициент отражения плотно упакованных пирамидок в приближении геометрической оптики. Результаты РБТП (точки) получены для Л/А = 15.

убывает при увеличении й/А. Мы объясняем это тем, что число переотражений внутри текстуры, необходимое для того, чтобы падающий луч поменял свое направление и отразился, линейно растет с увеличением высоты. В случае неполного замощения (например, плотно упакованные в треугольную решетку конусы на Рис. 7) при увеличении с(/Л отражение вначале падает, а после достижения некоторого минимума поднимается до постоянного значения. Мы связываем это с тем, что при ¿/А -» оо часть лучей после первоначального отражения от боковой поверхности движется почти параллельно ей и неминуемо попадает в зазор между основаниями, от которого эти лучи окончательно отражаются обратно.

Перейдем теперь к результатам РОТБ для всего диапазона размеров текстуры. Рассмотрим вначале случай полного замощения подложки основаниями пирамидок (Рис. 8). Зафиксируем какое-либо значение ¿/А и исследуем зависимость отражения Я от Л/Л (см. соответствующие черные кривые на Рис. 8). При увеличении Л/Л от 0 до 1 отражение падает и достигает локального минимума при Л/Л = 1, что качественно объясняется тем, эффективная диэлектрическая проницаемость в нулевом приближении не зависит от Л/Л, а отражение соответствующего промежуточного слоя падает при увеличении высоты й/Х = (¿/Л)(Л/Л). При Л/Л > 1 отражение продолжает падать, проходя ряд локальных минимумов, соответствующих тем значениям Л, при которых появляются очередные дифракционные порядки. Отметим, что при

— -10

-15

-20

-25

-30

-35 СО

":"'Г 0.1 "О

'Г 02 а -40 ф

/ л о

/ 0.5 Й .5 от

/ о

1 -50 ш 5=

2ь/д 1 ,5 0) СИ

-Ш

Рис. 8: Слева: Коэффициент отражения от плотно упакованных квадратных пирамидок при различных значениях Л/Ли ¿/Л.

Справа: Коэффициент отражения от плотно упакованных в треугольную решетку конусов при различных значениях Л/Л и «¿/Л. Результаты получены с помощью метода РОТО.

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

В случае неполного замощения, как, например, для случая плотно упакованных в треугольную решетку конусов (Рис. 8), наименьшее значение отражения достигается при А/Л ~ 1.

Итак, нами установлено, что оптимальные размеры текстуры зависят от характера замощения подложки основаниями пирамидок. В случае полного замощения, когда между основаниями нет зазоров, при фиксированной величине отношения ¿/Л наименьшее отражение достигается при Л —5- оо. Это согласуется с тем, что, как мы показали, скорость убывания отражения при росте (1 и фиксированных Л и Л в приближении геометрической оптики экспоненциальна, в то время как в приближении эффективной среды она в общем случае имеет степенной характер. При этом малая величина отражения достижима также и при Л порядка длины волны. В случае неполного замощения оптимальное значение Л имеет порядок длины волны. Наличие минимума в этом случае также находится в согласии с полученными нами оценками в приближениях эффективной среды и геометрической оптики.

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

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

3 Основные результаты и выводы работы

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

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

• С использованием трех различных методов найдены асимптотические зависимости величины отражения от геометрических параметров составляющих текстуру рассеивателей для длинно- и коротковолнового пределов. Установлено, что ключевым фактором, влияющим на оптимальный размер рассеивателей, является характер замощения подложки их основаниями. В случае полного замощения наименьшее значение отражения

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

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

• Предложен новый подход программной реализации метода решения уравнений Максвелла в конечных разностях (РОТБ), основанный на интегральном представлении уравнений Максвелла. В рамках этого подхода написана параллельная программа на языке программирования С++. В ходе ее написания реализован ряд оригинальных решений: оптимизация использования памяти с помощью упаковки используемых в разностных уравнениях коэффициентов; увеличение производительности посредством разбиения алгоритма обновления сеток на основной цикл, поддающийся векторизации, и специализированные поправки к нему, увеличение параллельной эффективности путем балансировки доменов. Специально проведенные тесты показали, что быстродействие программы примерно в два раза выше, чем у программы МЕЕР, разработанной в Массачусетском технологическом институте. Также был продемонстрирован линейный характер масштабируемости программы вплоть до тысячи вычислительных ядер. Развитый при написании программы подход обеспечивает возможность расчета с одновременным использованием нескольких сеток, каждая из которых максимально приспособлена к соответствующей ей области моделируемой структуры. Это является особо важным при моделировании оптических свойств структур сложной формы с неоднородностями, много меньшими длины волны. Программа находится в открытом доступе в интернете.

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

ную вычислительную сетку. Его применение приводит к увеличению порядка точности метода FDTD по сравнению с обычным „лестничным" приближением. Это является значимым для самых разнообразных оптических приложений: при моделировании оптических свойств металлических ФК, при расчете антиотражающих свойств текстурированных покрытий и т. д.

• Предложен и реализован итерационный метод для моделирования наклонного падения плоской волны на периодическую структуру в численном расчете FDTD. С его помощью можно получать результаты во всем диапазоне углов от 0° до 90°, а также наблюдать реальное поведение поля во времени. Это позволяет эффективно исследовать излуча-тельные характеристики металлических фотонных кристаллов, а также антиотражающие свойства текстурированных покрытий при любом угле падения.

Основные результаты диссертации опубликованы в следующих работах:

1. I. Valuev, A. Deinega, A. Knizhnik, В. Potapkin, "Creating Numerically Efficient FDTD Simulations Using Generic С++ Programming Lecture Notes in Computer Science, 4707, 213-226 (2007).

2. A. Deinega and I. Valuev, "Subpixel smoothing for conductive and dispersive media in the finite-difference time-domain method Optics Letters 32, 3429-3431 (2007).

3. I. Valuev, A. Deinega, and S. Belousov, "Iterative technique for analysis of periodic structures at oblique incidence in the finite-difference time-domain method Optics Letters 33, 1491-1493 (2008).

4. A. Deinega, I. Valuev, B. Potapkin, and Y. Lozovik,"Antireflective properties of pyramidally textured surfaces Optics Letters 35, 106-108 (2010).

5. А. В. Дейнега, И. В. Конистяпина, М. В. Богданова, И. А. Валуев, Ю. Е. Лозовик, Б. В. Потапкин, "Оптимизация антиотражающего слоя в солнечных батареях на основе первопринципных расчетов Известия высших учебных заведений. Физика., 11, 13-19 (2009).

6. С. А. Белоусов, М. В. Богданова, И. А. Валуев, А. В. Дейнега, С. Л. Эй-дерман, А. А. Книжник, И. Я. Полищук, Ю. Е. Лозовик, Б. В. Потапкин, Ю. А. Успенский, Э. Т. Кулатов, А. А. Титов, S. Zalyubovsky, В. Ramamurthi, "Предсказательное моделирование оптических свойств металло-диэлектрических

метаматериаловИзвестия высших учебных заведений. Физика., 11, 20-27 (2009).

Цитируемая литература

[1] Андриевский Р. А., Рагуля А. В. Наноструктурные материалы, — Москва: Издательский Центр Академия, 2005.

]2] Joannopoulos J. D., G.Johnson S., Winn J. N., Meade R. D. Photonic Crystals: Molding the Flow of Light. — Princeton Univ. Press, 2008.

[3] Fleming J. G., Lin S. Y., El-Kady I., Biswas R., Но К. M. All-metallic three-dimensional photonic crystals with a large infrared bandgap // Nature. — 2002. — Vol. 417. — Pp. 52—55.

[4] Improved broadband and quasi-omnidirectional anti-reflection properties with biomimetic silicon nanostrnctures / Y.-F. Huang, S. Chattopadhyay, Y.-J. Jen, C.-Y. Peng, T.-A. Liu, Y.-K. Hsu, C.-L. Pan et al. // Nature Nanotechnology. — 2007. — Vol. 2. — Pp. 770-774.

[5] Llopis F., Tobias I. Texture profile and aspect ratio influence on the front reflectance of solar cells U J. Appl. Phys. - 2006. - Vol. 100. - P. 124504.

[6] Taflove A., Hagness S. H. Computational Electrodynamics: The Finite Difference Time-Domain Method. — Boston: Artech House, 2005.

[7] Farjadpour A., Roundy D., Rodriguez A., Ibanescv M., Bermel P. Improving accuracy by subpixel smoothing in fdtd // Optics Letters. - 2006. — Vol. 31. — Pp. 2972-2974.

[8] Stefanou N., Yannopapas N., Modinos A. Heterostructures of photonic crystals: frequency bands and transmission coefficients // Comput. Phys. Commun. — 1998. —Vol. 113. —Pp. 49-77.

[9] Narayanaswamy A., Chen G. Thermal emission control with one-dimensional metallodielectric photonic crystal // Phys. Rev. B. - 2004, - Vol. 70. - P. 125101.

[10] Raguin D. Д., Morris G. M. Analysis of antireflection-structured surfaces with continuous one-dimensional surface profiles 11 Appl. Opt.- 1993.- Vol. 32,- Pp. 2582-2598.

[11] Sopori B. L., Pry or R. A. Design of antireflection coatings for textured silicon solar cells // Sol. Cells. - 1983. - Vol. 8. - Pp. 249-261.

[12] Brauer R., Bryngdahl O. Design of antireflection gratings with approximate and rigorous methods // Appl. Opt. - 1994. - Vol. 33. - Pp. 7875-7882.

[13] Wu F., Whites K. W. Computation of static effective permittivity for a multiphase lattice of cylinders 11 Electromagnetics. — 2001. — Vol. 21. — Pp. 97-114.

Подписано к печати ■ИМ.-/О Тираж 400 Заказ ¿9.

Отпечатано в отделе оперативной печати физического факультета МГУ

Введение

1 Развитие метода решения уравнений Максвелла в конечных разностях

1.1 Обзор литературы.

1.2 Метод нодсеточного сглаживания

1.3 Итерационный метод для моделирования наклонного падения плоской волны на периодическую структуру.

1.4 Выводы.

2 Применение шаблонного метапрограммирования для эффективной реализация FDTD

2.1 Контурный подход к дискретизации уравнений Максвелла

2.2 Код как совокупность взаимодействующих компонент

2.3 Стадии расчета.

2.4 Оптимизация использования памяти.

2.5 Проведение параллельных расчетов.

3 Металлические фотонные кристаллы как источники светового излучения

3.1 Обзор литературы.

2.6 Выводы.

3.2 Оптимизация излучательных характеристик металлических фотонных кристаллов в видимом диапазоне.

3.3 Учет внешней матрицы, необходимой для фиксации фотонного кристалла.

3.4 Выводы.

4 Антиотражающие текстурированные покрытия

4.1 Обзор литературы.

4.2 Численное моделирование текстурированных покрытий во всем диапазоне размеров текстуры.

4.3 Кремниевые текстурированные покрытия.

4.4 Выводы.

Диссертация посвящена численному моделированию оптических свойств наноструктурированных материалов, а именно, в ней исследуются следующие вопросы:

• Возможность применения металлических фотонных кристаллов в качестве новых высокоэффективных источников света.

• Оптимизация размеров и формы антиотражающих нанотекстуриро-ванных покрытий.

Основным используемым численным методом является метод решения уравнений Максвелла в конечных разностях (Finite-Difference Time-Domain, FDTD). Для расчетов применяется специально написанная параллельная программа, включающая в себя ряд новых оригинальных численных методов, которые также описываются в данной диссертации.

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

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

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

В этой работе основным используемым численным методом является метод решения уравнений Максвелла в конечных разностях (англ. - Finite-Difference Time-Domain method, FDTD) [6]. Выбор этого метода вызван присущей ему высокой параллельной эффективностью, что позволяет применять его для расчета больших задач на кластерных вычислительных системах.

Как и всякий другой разностный метод, FDTD сталкивается с проблемой неточного отображения на прямоугольную вычислительную сетку тел, обладающих произвольной формой поверхности. Решение этой проблемы особенно значимо для моделирования исследуемых в данной работе фотонных кристаллов и текстурированных поверхностей. Одним из наиболее удачных способов ее решения является метод подсеточного сглаживания [7]. Однако этот метод в его первоначальной постановке может применяться только для диэлектриков, в связи с чем актуальна его модификация, которая позволила бы расширить его на дисперсные среды.

Исследование оптических свойств наноматериалов включает в себя получение спектров прохождения и отражения от периодических структур при падении па них плоской волны. Их расчет FDTD для случая нормального падения предполагает применение периодических граничных условий. В случае наклонного падения периодические граничные условия содержат временной сдвиг, что затрудняет их применение. В литературе было предложено несколько методов для решения этой проблемы, однако у каждого из них есть свои недостатки: нестабильность для углов близких к 90°, увеличение объема используемой памяти и т. д. В связи с этим актуальна разработка альтернативного метода, который был бы лишен недостатков предыдущих. Такой метод позволил бы эффективно исследовать излуча-тельные характеристики фотонных кристаллов, а также антиотражающие свойства текстурированных покрытий при любом угле падения.

Цель работы состоит в 1) исследовании излучательных характеристик металлических фотонных кристаллов, а также возможности их практического применения в качестве источников света с высокой селективностью излучения в видимом диапазоне; 2) исследовании антиотражающих свойств текстурированных покрытий во всем диапазоне размеров составляющих текстуру рассеивателей, а также установлении оптимального размера для заданного диапазона длин волн; 3) дальнейшем развитии метода FDTD, включающем, в частности, разработку метода подсеточного сглаживания для дисперсных сред и метода расчета наклонного падения плоской волны на периодическую структуру.

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

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

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

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

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

5. Предложен новый подход программной реализации метода решения уравнений Максвелла в конечных разностях (FDTD), основанный на интегральном представлении уравнений Максвелла. В рамках этого подхода написана параллельная программа на языке программирования С++. В ходе ее написания реализован ряд оригинальных решений: оптимизация использования памяти с помощью упаковки используемых в разностных уравнениях коэффициентов; увеличение производительности посредством разбиения алгоритма обновления сеток на основной цикл, поддающийся векторизации, и специализированные поправки к нему; увеличение параллельной эффективности путем балансировки доменов. Специально проведенные тесты продемонстрировали линейный характер масштабируемости программы вплоть до тысячи вычислительных ядер. Программа находится в открытом доступе в интернете (http://fdtd.kintechlab.com).

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

7. Впервые предложен и реализован итерационный метод для моделирования наклонного падения плоской волны на периодическую структуру в численном расчете FDTD. С его помощью можно получать результаты во всем диапазоне углов от 0° до 90°, а также наблюдать реальное поведение поля во времени. Помимо этого, он лишен таких недостатков предложенных ранее методов, как неустойчивость при отвесных углах падения и увеличение требуемого размера сетки. Особенностью метода является то, что получение решения для одного угла требует нескольких отдельных итераций-экспериментов.

Положения, выносимые на защиту.

1. Геометрические параметры вольфрамовых фотонных кристаллов типа прямого опала, при которых КПД источника света на основе фотонных кристаллов выше в несколько раз по сравнению с обычными лампами накаливания.

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

3. Асимптотические зависимости величины отражения от геометрических параметров текстурироваиной поверхности для длинно- и коротковолнового пределов.

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

5. Параллельная программа FDTD, которая может использоваться для предсказательного моделирования в разнообразных оптических приложениях.

6. Новые численные методы в рамках FDTD: метод подсеточного сглаживания для дисперсных сред и итерационный метод для моделирования наклонного падения плоской волны на периодическую структуру.

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

1. Найдены геометрические параметры вольфрамовых фотонных кристаллов типа прямого опала, при которых КПД источника света на основе фотонного кристалла оказывается выше в несколько раз по сравнению с обычными лампами накаливания.

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

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

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

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

6. Предложены и реализованы численные методы, позволяющие существенно увеличить быстродействие и точность расчетов FDTD.

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

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

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

Апробация работы. Основные результаты диссертационной работы были доложены на научно-технических конференциях „The International Conference on Computational Science and Its Applications" (Малайзия, Kuala Lumpur, 2007), „Conference on Computational Physics" (Бразилия, Ouro Preto, 2008), „Conference on Computational Physics" (Тайвань, Kaohsiung, 2009), „Многомасштабное моделирование процессов и структур в нанотех-иологиях" (Москва, 2009), „VII Курчатовская молодежная научная школа" (Москва, 2009). 52-ая научная конференция МФТИ (Москва, 2009).

Реализованная программа FDTD стала призером на конкурсе „Максимальная масштабируемость", проведенном в рамках РОСНАНОФОРУМА-2009 компаниями Intel и "Роснано".

Публикации. По материалам диссертации опубликовано 5 работ в реферируемых научных изданиях [8-12], 1 в нерефериуемом научном издании [13] и тезисы российских и международных конференций.

Основные результаты и выводы работы

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

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

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

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

• Предложен новый подход программной реализации метода решения уравнений Максвелла в конечных разностях (FDTD), основанный на интегральном представлении уравнений Максвелла. В рамках этого подхода написана параллельная программа на языке программирования С++. В ходе ее написания реализован ряд оригинальных решений: оптимизация использования памяти с помощью упаковки используемых в разностных уравнениях коэффициентов; увеличение производительности посредством разбиения алгоритма обновления сеток на основной цикл, поддающийся векторизации, и специализированные поправки к нему; увеличение параллельной эффективности путем балансировки доменов. Специально проведенные тесты показали, что быстродействие программы примерно в два раза выше, чем у программы МЕЕР, разработанной в Массачусетском технологическом институте. Также был продемонстрирован линейный характер масштабируемости программы вплоть до тысячи вычислительных ядер. Развитый при написании программы подход обеспечивает возможность расчета с одновременным использованием нескольких сеток, каждая из которых максимально приспособлена к соответствующей ей области моделируемой структуры. Это является особо важным при моделировании оптических свойств структур сложной формы с неоднородностями, много меньшими длины волны. Программа находится в открытом доступе в интернете.

• В рамках метода FDTD предложен и реализован алгоритм подсеточ-ного сглаживания для дисперсных сред, позволяющий существенно уменьшить величину численной ошибки, связанной с неточным отображением тел, обладающих произвольной формой поверхности, на прямоугольную вычислительную сетку. Его применение приводит к увеличению порядка точности метода FDTD по сравнению с обычным „лестничным" приближением. Это является значимым для са мых разнообразных оптических приложений: при моделировании оптических свойств металлических ФК, при расчете антиотражающих свойств текстурированных покрытий и т. д.

• Предложен и реализован итерационный метод для моделирования наклонного падения плоской волны на периодическую структуру в численном расчете FDTD. С его помощью можно получать результаты во всем диапазоне углов от 0° до 90°, а также наблюдать реальное поведение поля во времени. Это позволяет эффективно исследовать из-лучательные характеристики металлических фотонных кристаллов, а также антиотражающие свойства текстурированных покрытий при любом угле падения.

5.2. Благодарности

Хочу выразить особую благодарность моему руководителю - Потап-кину Борису Васильевичу за постановку темы и плодотворное обсуждение результатов. Отдельную благодарность выражаю моим соавторам Юрию Ефремовичу Лозовику, Илье Валуеву, Марии Богдановой, Сергеям Бело-усову и Эйдерману за стимулирующие дискуссии.

5. Заключение

1. Андриевский Р. А., Рагуля А. В. Наноструктурные материалы.— Москва: Издательский Центр Академия, 2005.

2. Joannopoulos J. D., G.Johnson S., Winn J. N., Meade R. D. Photonic Crystals: Molding the Flow of Light. — Princeton Univ. Press, 2008.

3. Fleming J. G., Lin S. Y., El-Kady I., Biswas R., Но К. M. All-metallic three-dimensional photonic crystals with a large infrared bandgap // Nature. 2002. - Vol. 417. - Pp. 52—55.

4. Llopis F., Tobias I. Texture profile and aspect ratio influence on the front reflectance of solar cells // J. Appl. Phys. 2006. - Vol. 100. — P. 124504.

5. Taflove A., Hagness S. H. Computational Electrodynamics: The Finite Difference Time-Domain Method. — Boston: Artech House, 2005.

6. Farjadpour A., Roundy D., Rodriguez A., Ibanescu M., Bermel P. Improving accuracy by subpixel smoothing in fdtd // Optics Letters. — 2006.— Vol. 31.-Pp. 2972-2974.

7. Deinega A., Valuev I. Subpixel smoothing for conductive and dispersive media in the fdtd method // Optics Letters. — 2007. — Vol. 32. — Pp. 3429-3431.

8. Valuev I., Deinega A., Belousov S. Iterative technique for analysis of periodic structures at oblique incidence in the finite-difference time-domain method // Optics Letters. 2008. - Vol. 33. - Pp. 1491-1493.

9. Deinega A., Valuev I., Potapkin В., Lozovik Y. Antireflective properties of pyramidally textured surfaces // Optics Letters. — 2010. — Vol. 35. — Pp. 106-108.

10. Дейнега А. В., Конистяпина И. В., Богданова М. В., Валуев И. А., Лозовик Ю. Е.} Потапкин Б. В. Оптимизация антиотражающего слоя в солнечных батареях на основе первопринципных расчетов // Известия Вузов. Физика. — 2009. — Т. 11. — С. 13-19.

11. Valuev I., Deinega A., Knizhnik A., Potapkin В. Creating numerically ■ efficient fdtd simulations using generic C++ programming // Lecture Notesin Computer Science. 2007. - Vol. 4707. — Pp. 213-226.

12. Yee K. S. Numerical solution of inital boundary value problems involving maxwell's equations in isotropic media // IEEE Trans. Antennas and Propagation. — 1966. — Vol. 14. — Pp. 302-307.

13. Okoniewski M., Mrozowski M., Stuchly M. A. Simple tratment of multi-term dispersion in fdtd // IEEE Microware and Guided Wave Lett. — 1997. Vol. 7. - Pp. 121-123.

14. Joseph R. M., Hagness S. C., Taflove A. Direct time integration of maxwell's equations in linear dispersive media with absorption for scattering and propagation of femtosecond electromagnetic pulses // Opt. Lett. — 1991, —Vol. 16.-Pp. 1412-1414.

15. Umashankar К. R., Taflove A. A novel method to analyze electromagnetic scattering of complex objects // IEEE Transactions on Electromagnetic Compatibility. 1982. - Vol. 24. — Pp. 397—405.

16. Berenger J. P. A perfectly matched layer for the absorption of electromagnetic waves // J. Comput. Phys.— 1994, — Vol. 114. — Pp. 185-200.

17. Sacks Z. S., Kingsland D. M., Lee R., Lee J.-F. A perfectly matched anisotropic absorber for use as an absorbingboundary condition // IEEE Trans. Antennas and Propagation. — 1995. — Vol. 43. — Pp. 1460-1463.

18. Gedney S. D. An anisotropic perfectly matched layer absorbing media for the truncation of fdtd latices // IEEE Trans. Antennas and Propagation. 1996. - Vol. 44. — Pp. 1630-1639.

19. Chew W. C., Weedon W. H. A 3d perfectly matched medium from modified maxwell's equations with stretched coordinates // Microwave Optical Tech. Letters. 1994. - Vol. 7. - Pp. 590-604.

20. Teixeira F. L., Chew W. C. General closed-form pml constitutive tensors to match arbitrary bianisotropic and dispersive linear media // IEEE Microwave and Guided Wave Letters. — 1998. Vol. 8. - Pp. 223-225.

21. Bohren C. F., Huffman D. R. Absorption and Scattering of Light by Small Particles. — New York: Wiley-Interscience, 1983.

22. Zivanovic S. S., Yee K. S., Mei К. K. A subgridding method for the time domain finite-difference method to solve maxwell's equations // IEEE Trans. Microware Theory Tech. — 1991. — Vol. 38. — Pp. 471-479.

23. Jurgens T. G., Taflove A., Umashankar K., Moore T. G. Finite-difference time-domain modeling of curved surfaces // IEEE Trans. Antennas Propag. 1992. - Vol. 40. - Pp. 357-366.

24. Shankar VMohammadian A., Hall W. F. A time-domain finite-volume treatment for the maxwell equations // Electromagnetics. — 1990. — Vol. 10. Pp. 127-145.

25. Mohammadi A., Nadgaran H., Agio M. Contour-path effective permet-tivities for the two-dimensional finite-difference time-domain method // Optics Express. 2005. - Vol. 13. — Pp. 10367-10381.

26. Lee J.-Y., Myung N.-H. Locally tensor conformal fdtd method for modeling arbitrary dielectric surface // Microw. Opt. Technol. Lett. — 1999. — Vol. 23. Pp. 245-249.

27. Nadobny J., Sullivan D., Wlodarczyk W., Deuflhard P., Wust P. A 3d tensor fdtd-formulation for treatment of slopes interfaces in electrically inhomogeneous media // IEEE Trans. Antennas Propag.— 2003.— Vol. 51.-Pp. 1760-1770.

28. Dey S., Mittra R. A coonformal finite-difference time-domain technique for modeling cylindrical dielectric resonators // IEEE Trans. Microware Theory Tech. 1999. - Vol. 47. - Pp. 1737-1739.

29. Kaneda N., Houshmand В., Itoh T. Fdtd analysis of dielectric resonators with curved surfaces // IEEE Trans. Microware Theory Tech. — 1997. — Vol. 45. Pp. 1645-1649.

30. Okoniewski M., Mrozowski M., Stuchly M. A. Simple tratment of multi-term dispersion in fdtd // IEEE Microware and Guided Wave Lett. — 1997. Vol. 7. - Pp. 121-123.33. van de Hulst H. C. Light Scattering by Small Particles. — New York: Wiley, 1957.

31. Taflove A., Umashankar K. R. Radar cross-section of general three-dimensional structures // IEEE Trans. Electromagn. Compat.— 1983,— Vol. 25. Pp. 433-440.

32. Ordal M. A., Long L. L., Bell R. J., Bell S. E. Optical properties of the metals al, со, au, fe, pb, ni, pd, pt, ag, ti, and w in the infrared and far infrared // Appl. Opto. — 1983. —Vol. 22.-Pp. 1099-1119.

33. Kao Y. С. A., Atkins R. G. A finite-difference time-domain approach for frequency selective surfaces at oblique incidence // Proc. 1996 IEEE Antennas and Propagation Society International Symposium, MD, July 1996. 1996. - Vol. 2. - Pp. 1432-1435.

34. Roden J. A., Gedney S. D., Kesler M. PMaloney J. G.} Harms P. II. Time-domain analysis of periodic structures at oblique incidence: orthogonal and nonorthogonal fdtd implementations // Microwave Theory and Techniques. 1998. - Vol. 46. - Pp. 420-427.

35. Harms P., Mittra R., Ко W. Implementation of the periodic boundary condition in the finite-difference time-domain algorithm for fss structures // IEEE Trans. Antennas and Propagation. — 1994. — Vol. 42. — Pp. 1317-1324.

36. Aminian A., Rahmat-Samii Y. Spectral fdtd: a novel technique for the analysis of oblique incident plane wave on periodic structures // IEEE Trans. Antennas and Propagation. — 2006, — Vol. 54.— Pp. 1818-1825.

37. Ren J., Gandhi 0. P., Walker L. R., Fraschilla J., Boerman C. R. Floquet-based fdtd analysis of two-dimensional phased array antennas // IEEE Microwave and Guided Wave Letters.— 1994.— Vol. 109,— Pp. 109-111.

38. Stefanou N., Yannopapas N., Modinos A. Heterostructures of photonic crystals: frequency bands and transmission coefficients // Comput. Phys. Commun. 1998. - Vol. 113. - Pp. 49-77.

39. Stroustrup B. What is object-oriented programming? // IEEE Software. — 1988.-Vol. 5.-Pp. 10-20.

40. Корнеев В. Д. Параллельное программирование в MPI. — Новосибирск: Издательство ИВМиМГ СО РАН, 2002.

41. Страуструп Б. Язык программирования С++. — Москва: Издательство Бином, 2005.

42. Yablonovitch E. Inhibited spontaneous emission in solid-state physics and electronics // Phys. Rev. Lett. 1987. — Vol. 58. — Pp. 2059—2062.

43. Sakoda K. Optical Properties of Photonic Crystals.— Berlin: Springer, 2001.

44. Белотелое В. И., Звездин А. К. Фотонные кристаллы и другие мета-материалы.— Москва: Наука, 2006.

45. Narayanaswamy A., Chen С. Thermal emission control with one-dimensional metallodielectric photonic crystal // Phys. Rev. B. — 2004. — Vol. 70.-P. 125101.

46. Li Z.-Y. Modified thermal radiation in three-dimensional photonic crystals // Phys. Rev. B. 2002. - Vol. 66. - P. 241103.

47. Celanovic I., Perreault D., Kassakian J. Resonant-cavity enhanced thermal emission // Phys. Rev. B. 2008. — Vol. 72. — Pp. 075127—075132.

48. Florescu M., Lee H., Stimpson A. J., Bowling J. Thermal emission and absorption of radiation in finite inverted-opal photonic crystals // Phys. Rev. A. 2005. - Vol. 72. - P. 033821.

49. Yannopapas V. Thermal emission from three-dimensional arrays of gold nanoparticles // Phys. Rev. В. — 2009.-Vol. 73.—P. 113108.

50. Chan D. L. C., Soljacic M., Joannopoulos J. D. Thermal emission and design in 2d-periodic metallic photonic crystal slab // Optics Express. — 2006. Vol. 14. - Pp. 8785-8796.

51. Enoch S., Simon J.-J., Escoubas L., Elalmy Z., Lemarquis F., Torchio P., Albrand G. Simple layer-by-layer photonic crystal for the control of thermal emission // Appi Phys. Lett.— 2005.— Vol. 86.— Pp. 261101— 261103.

52. Lee В. J., Fu С. J., Zhang Z. M. Coherent thermal emission from one-dimensional photonic crystals // Appl. Phys. Lett. — 2005.— Vol. 87.— P. 071904.

53. Gall J. L., Olivier M., Greffet J.-J. Experimental and theoretical study of reflection and coherent thermal emissionby a sic grating supporting a surface-phonon polariton // Phys. Rev. B. — 1997. — Vol. 55. — Pp. 10105-10114.

54. Laroche M., Carminati R., Greffet J.-J. Coherent thermal antenna using a photonic crystal slab // Phys. Rev. Lett. — 2006. — Vol. 96. P. 123903.

55. Tsai M.-W., Chuang T.-H., Meng C.-Y., Chang Y.-T., Lee S.-C. High performance midinfrared narrow-band plasmonic thermal emitter // Appl. Phys. Lett. 2006. - Vol. 89. - P. 173116.

56. Cornelius С. M., Dowling J. P. Modification of planck blackbody radiation by photonic band-gap structures // Phys. Rev. A.— 1999. — Vol. 59. Pp. 4736-4746.

57. Pigeat P., Rouxel D., Weber B. Calculation of thermal emissivity for thin films by a direct method // Phys. Rev. B. — 1998. — Vol. 57. — Pp. 92939300.

58. Luo C., Narayanaswamy A., Chen G., Joannopoulos J. D. Thermal radiation from photonic crystals: A direct calculation // Phys. Rev. Lett. — 2004. Vol. 93. - Pp. 213905—213908.

59. Ландау Л. Д., Лифшиц Е. М. Теоретическая физика. — Москва: Наука, 1976.-Т. 5.

60. Рытое С. М. Введение в статистическую радиофизику. — Москва: Наука, 1978.

61. Шёпф Х.-Г. От Кирхгофа до Планка. — Москва: Мир, 1981.

62. Greffet J.-J., Nieto-Vesperinas M. Field theory for generalized bidirectional reflectivity: derivation of helmholtz's reciprocity principle and kirch-hoff's law // J. Opt. Soc. Am. A. — 1998. Vol. 15. - Pp. 2735-2744.

63. Сивухин Д. В. Общий курс физики.— Москва: Физматлит, 2006.— Т. 4.

64. Lin S. Y., Moreno J., Fleming J. G. Three-dimensional photonic-crystal emitter for thermal photovoltaic power generation // Appl. Phys. Lett.— 2003. Vol. 83. - Pp. 380-382.

65. Coutts T. J., Fitzgerald M. C. Thermophotovoltaics // Sci. Am.— 1998.-Vol. 279.-P. 90.

66. Sai H., Yugami H., Akiyama Y., Kanamori Y., Hane K. Spectral control of thermal emission by periodic microstructured surfaces in the near-infrared region // J. Opt. Soc. Am. A. — 2001.— Vol. 18.— Pp. 1471— 1476.

67. Maruyama S., Kashiwa Т., Yugami H., Esashi M. Thermal radiation from two-dimensionally confined modes in microcavities // Appl. Phys. Lett. — 2001. Vol. 79. - Pp. 1393-1395.

68. Sai H., Kanamori Y., Yugami H. High-temperature resistive surface grating for spectral control of thermal radiation // Applied Physics Letters. — 2003. Vol. 82. - Pp. 1685-1687.

69. Photonic crystal enhanced narrow-band infrared emitters / M. U. Pralle, N. Moelders, M. P. McNeal, I. Puscasu, A. C. Greenwald, J. T. Daly, E. A. Johnson et al. // Appl. Phys. Lett. — 2002. — Vol. 81. — Pp. 4685—4687.

70. Lin S. Y.; Fleming J. G., Li Z. Y., El-Kady I., Biswas R., Но К. M. Origin of absorption enhancement in a tungsten, three-dimensional photonic crystal // J. Opt. Soc. Am. В.- 2003. Vol. 40.-Pp. 1538-1541.

71. Lin S. Y., Fleming J. G., El-Kady I. Experimental observation of photonic-crystal emission near a photonic band edge // Appl. Phys. Lett. 2003. - Vol. 83. - Pp. 593-595.

72. Roberts S. Optical properties of nickel and tungsten and their interpretation according to drude's formula // Phys. Rev. — 1959.— Vol. 114,— Pp. 104-115.

73. Ландсберг Г. С. Оптика. — Москва: Физматлит, 2003.

74. Лифшиц И. М., Гредескул С. А., Пастур Л. А. Введение в теорию неупорядоченных систем // ЖЭТФ.— 1982,— Т. 83, № 6,— С. 23622376.

75. Ворн М., Вольф Э. Основы Оптики. — Москва: Наука, 1973.

76. Knittl Z. Optics of thin films. — London: John Wiley, 1976.

77. Macleod H. A. Thin film optical filter. — New York: McGraw-Hill, 1989.

78. Minot M. J. Single-layer, gradient refractive index antireflection films effective from 0.35 to 2.5 mkrn // J. Opt. Soc. Am.— 1976.- Vol. 66,— Pp. 515-519.

79. Jacobson R. Inhomogeneous and coevaporated homogeneous films for optical applications // Physics of Thin Films / Ed. by G. Hass, M. H. Francombe, R. W. Hoffman. — New York: Academic, 1975. — Vol. 8. — Pp. 51-98.

80. Bernhard C. G. Structural and functional adaptation in a visual system // Endeavour. 1967. - Vol. 26. - Pp. 79-84.

81. Bernhard C. G.; Miller W. H., Moller A. R. The insect corneal nipple array // Acta Physiol. Scand. — 1965, — Vol. 63. — Pp. 1-25.

82. Wilson S. J., Hutley M. C. The optical properties of 'moth eye' antireflection surfaces // Opt. Acta. — 1982. — Vol. 29. Pp. 993-1009.

83. Grann E. В., Moharam M. G., Pommet D. A. Optimal design for antire-flective tapered two-dimensional subwavelength grating structures // J. Opt. Soc. Am. A. — 1995. — Vol. 12.- Pp. 333-339.

84. Raguin D. H., Morris G. M. Analysis of antireflection-structured surfaces with continuous one-dimensional surface profiles // Appl. Opt. — 1993. — Vol. 32. Pp. 2582-2598.

85. Sopori B. L., Pryor R. A. Design of antireflection coatings for textured silicon solar cells // Sol. Cells. 1983. - Vol. 8. - Pp. 249-261.

86. Campbell P., Green M. A. Light trapping properties of pyramidally textured surfaces // J. Appl. Physics.— 1987. — Vol. 62. — Pp. 243-249.

87. Chen H. L., Chuang S. Y., Lin С. H., Lin Y. H. Using colloidal lithography to fabricate and optimize sub-wavelength pyramidal and honeycomb structures in solar cells // Opt. Express. — 2007. — Vol. 15. — Pp. 1479314803.

88. Kanarnori Y., Sasaki M., Plane K. Broadband antireflection gratings fabricated upon silicon substrates // Opt. Lett. — 1999. — Vol. 24. — Pp. 1422-142.

89. Kanamori Y., Hane K., Sai H., Yugami H. 100 nm period silicon antireflection structures fabricated using a porous alumina membrane mask // Appl. Phys. Lett. 2001. - Vol. 78. - Pp. 142-143.

90. Yu Z., Gao H.; Wu W., Ge H., Chou S. Y. Fabrication of large area subwavelength anti-reflection structures on si using trilayer resist nanoim-print lithography and lift-off // Journal of Vacuum. Science and Tech. B. 2003. - Vol. 21. - Pp. 2874-2877.

91. Song Y. M., Bae S. Y., Yu J. S., Lee Y. T. Closely packed and aspect-ratio-controlled antireflection subwavelength gratings on gaas using a lenslike shape transfer // Opt. Lett. — 2009. Vol. 34. - Pp. 1702-1704.

92. Hadobas К., Kirsch S., Carl A., Acet M., Wassermann E. F. Reflection properties of nanostructure-arrayed silicon surfaces // Nanotechnology. — 2000.-Vol. 11.-Pp. 161-164.

93. Lee C., Bae S. Y., Mobasser S., Manohara H. A novel silicon nanotips antireflection surface for the micro sun sensor // Nanoletters. — 2009. — Vol. 5. Pp. 2438-2442.

94. Boden S. А., В agnail D. M. Tunable reflection minima of nanostructured antireflective surfaces // Appl. Phys. Lett. — 2008. — Vol. 93. — P. 133108.

95. Lalanne P., Lemercier-Lalanne D. On the effective medium theory of subwavelength periodic structure //J. Mod. Opt. — 1994.— Vol. 43,— Pp. 7875-7882.

96. Tao R., Chen Z., Sheng P. First-principles fourier approach for the calculation of the effective dielectric constant of periodic composites // Phys. Rev. B. 1990. - Vol. 41. - Pp. 2417-2420.

97. Haggans C. W.? Li L.; Kostuk R. K. Effective-medium theory of zeroth-order lamellar gratings in conical mountings //J. Opt. Soc. Am. A.— 1993. Vol. 10. - Pp. 2217-2225.

98. Grann E. ВMoharam M. G., Pommet D. A. Artificial uniaxial and biaxial dielectrics with use of two-dimensional subwavelength binary gratings // J. Opt. Soc. Am. A. 1994. - Vol. 11. - Pp. 2695-2703.

99. Brauer R., Bryngdahl O. Design of antireflection gratings with approximate and rigorous methods // Appl. Opt. — 1994. — Vol. 33. — Pp. 78757882.

100. Wu F., Whites K. W. Computation of static effective permittivity for a multiphase lattice of cylinders // Electromagnetics. — 2001,— Vol. 21.— Pp. 97-114.

101. Lalanne P., Lemercier-Lalanne D. Depth dependence of the effective properties of subwavelength gratings // J. Opt. Soc. Am. A. — 1997. — Vol. 14. Pp. 450-458.

102. Raguin D. H., Morris G. M. Antireflection structured surfaces for the infrared spectral region // Appl. Opt.— 1993.— Vol. 32.— Pp. 1154— 1167.

103. Abouelsaood A. A., El-Naggar S. A., Ghannam M. Y. Shape and size dependence of the anti-reflective and light-trapping action of periodic grooves // Prog. Photovolt: Res. Appl. — 2002. — Vol. 10. — Pp. 513-526.

104. Franceschetti G. Scattering from plane layered media // IEEE Trans. Antennas Propag. — 1964. — Vol. 12. — Pp. 754-763.

105. Kraus J. D. Electromagnetics. — New York: McGraw-Hill, 1984.

106. Moharam M. G., Gay lord Т. K. Diffraction analysis of dielectric surface-relief gratings // J. Opt Soc. Am.— 1982,- Vol. 72. Pp. 1385-1392.

107. Ono Y., Kirnura Y., Ohta Y., Nishida N. Antireflection effect in ultrahigh spatial-frequency holographic relief gratings // Applied Optics. — 1987. — Vol. 26.-Pp. 1142-1146.

108. Southwell W. H. Gradient-index antireflection coatings // Opt. Lett.— 1983. Vol. 8. - Pp. 584-586.

109. Southwell W. H. Pyramid-array surface-relief structures producing antireflection index matching on optical surfaces //J. Opt. Soc. Am. A.— 1991. Vol. 8. - Pp. 549-553.

110. Smith A. W., Rohatgi A., Neel S. C. Texture: a ray tracing program for the photovoltaic community // Photovoltaic Specialists Conference, 1990., Conference Record of the Twenty First IEEE. — 1990. — Vol. 1. — Pp. 426-431.

111. Thorp D., Wenham S. R. Ray-tracing of arbitrary surface textures for light-trapping in thin silicon solar cells // Solar Energy Materials and Solar Cells. 1997. - Vol. 48. - Pp. 295-301.

112. Correig Z., Gomez J. L., Calderer J. Optical analysis of textured surfaces for photovoltaic solar cells // Eleetrotechnieal Conference, 1991. Proceedings., 6th Mediterranean. — 1991. — Vol. 1, — Pp. 202-205.

113. Bucci 0., Franceschetti G. Scattering from wedge-tapered absorbers // IEEE Trans. Antennas Propag.— 1971, —Vol. 19. — Pp. 96-104.

114. Aroutiounian V. M., Martirosyan K., Soukiassian P. Almost zero reflectance of a silicon oxynitride/porous silicon double layer antireflection coating for silicon photovoltaic cells //J. Phys. D. — 2006. — Vol. 39. — Pp. 1623-1625.

115. Striemer С. C., Fauchet P. M. Dynamic etching of silicon for broadband antireflection applications // Appl. Phys. Lett. — 2002. — Vol. 81. — Pp. 2980-2982.

116. Lalanne P., Morris G. M. Antireflection behavior of silicon subwavelength periodic structures for visible light // Nanotechnology. — 1997. — Vol. 8. — Pp. 53-56.

117. Green M. A., Keevers M. Optical properties of intrinsic silicon at 300 к // Progress m Photovoltaics.— 1995, — Vol. 3, no. 3. — Pp. 189-192.

118. Hovel H. J. Solar Cells, Semiconductors and Semimetals. — New York: Academic, 1975. Vol. 11. — 38 pp.

119. Даффет-Смит П. Практическая астрономия с калькулятором. — Москва: Мир, 1982.

120. Вопап G. В. Ecological Climatology: Concepts and Applications. — Cambridge: Cambridge University Press, 2008.