Распространение волн в трехмерном периодическом диэлектрическом волноводе тема автореферата и диссертации по математике, 01.01.02 ВАК РФ

САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

на правах рукописи

<

Денисова Ирина Валерьевна

РАСПРОСТРАНЕНИЕ ВОЛН В ТРЁХМЕРНОМ ПЕРИОДИЧЕСКОМ ДИЭЛЕКТРИЧЕСКОМ ВОЛНОВОДЕ

специальность 01.01.02 - дифференциальные уравнения

АВТОРЕФЕРАТ

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

Санкт-Петербург 2003

Работа выполнена на кафедре математической физики математико-механического факультета Санкт-Петербургского государственного университета.

НАУЧНЫЙ РУКОВОДИТЕЛЬ доктор физико-математических наук, профессор Дергузов Виктор Иванович

ОФИЦИАЛЬНЫЕ ОППОНЕНТЫ доктор физико-математических наук, профессор Пламеневский Борис Алексеевич

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

ВЕДУЩАЯ ОРГАНИЗАЦИЯ Санкт-Петербургское отделение Математического института им. В.А. Стеклова РАН

Защита состоится «30 » 2003г. в^$ часов на заседании

диссертационного совета Д 212.232.49 по защитам диссертаций на соискание учёной степени доктора наук при Санкт-Петербургском государственном университете по адресу: 198504, Санкт-Петербург, Петродворец, Университетский пр., дом 28, математико-механичесий факультет СПбГУ. Аудитория 4526.

С диссертацией можно ознакомиться в Научной библиотеке им. М. Горького Санкт-Петербургского государственного университета по адресу: 199034, Санкт-Петербург, Университетская набережная, дом 7/9.

Автореферат разослан «_»_2003г.

Учёный секретарь диссертационного совета Д 212.323.49, профессор

А. А. Архипова

^ооЗИА

2.75"

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

Актуальность работы. Диэлектрические волноводы с периодической структурой находят широкое применение при модуляции, демодуляции. и фильтрации световых сигналов в различных устройствах интегральной оптики, включая фильтрующие устройства и решёточные элементы связи [1]. Волноводы этого типа продольно неоднородны. Поэтому теория направляемых волн для них, строго говоря, неприменима. Однако, когда пространственный период достаточно мал, волновой процесс в хорошем приближении предстаёт как обычная направляемая волна, длина которой намного больше пространственного периода [2]. Таким образом, свободный волновой процесс в периодической системе можно рассматривать как наложение бесконечной совокупности плоских неоднородных волн с разными постоянными распространения. Совокупность этих постоянных образует спектр волновода.

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

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

ШУ(дУЙ) + ЬМ = 0 (1)

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

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

Научная новизна и значимость работы. Представленные в диссертации результаты получены в период с 2000 по 2003 год. Все они являются новыми. Работа носит теоретический характер. Результаты работы позволяют более полно оценить спектральные свойства диэлектрического периодического волновода. Она дополняет теорию спектральных свойств диэлектрических периодических волноводов, развитую в работах [1] - [9]. Достоинством работы является более точное качественное описание непрерывного спектра трёхмерного диэлектрического периодического волновода.

Практическая значимость работы определяется возможностью применения полученных результатов для посташзмеи-^сзТйш^влЯ^ по

з 1 ^^/Г]

определению дискретного и непрерывного спектров трёхмерных диэлектрических периодических волноводов.

Апробация работы. Результаты диссертационной работы докладывались на общегородском семинаре по математическим вопросам теории дифракции в Санкт-Петербургском отделении Математического института РАН им. В.А. Стеклова.

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

Структура и объём работы. Диссертационная работа объёмом 133 страницы состоит из четырёх глав и списка литературы, содержащего 32 наименования.

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

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

ди"

[и] = 0 и

а дп

= 0> (2)

где [у]— скачок функции V на поверхности Г, —— производная функ-

дп

ции и по внешней нормали к поверхности Г.

Коэффициенты а и Ъ уравнения (1) удовлетворяют условиям а>0 и Ъ>.0. Вне области, ограниченной поверхностью Г, они равны постоянным Оо>0 и ¿„>0.

Разрешимость изучается в классе функций вида Флоке, то есть пред-ставимых в виде

и = е(3) где £— комплексное число, а функция V 2к - периодична по переменной

г.

Множество значений при которых существует решение задачи (1) - (3), образует спектр волновода.

Спектр волновода делится на непрерывную и дискретную составляющие. К дискретному спектру относятся те точки £, для которых решение поставленной задачи квадратично суммируемо на полосе П =

{(*,>>, г)г: 0 < г < 2л-}', а к непрерывному — те точки £, для которых оно ограничено на полосе П, но не квадратично суммируемо.

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

Задача (1) - (3) сводится к системе интегральных уравнений

и(х,у,г) = -(¡Б - |-УаУы/с (К1 +

ni «о aJ

ик(К1 + и0{х,у, z),

(4)

a(x,y,z)-a0

fi VeVn—<« + (

i a dn i

\s

ff-ds-I dn

4 ^

dk 8u0 8n dn

в пространстве W¡(Q.)xL2(S), где П = Пл/',5 = ГпП, P офаниченная поверхностью Г,

j. = 2a!l-adu

область,

а0 + а дп

Щ=е í;v0, v0-

решение уравнения

A)(ÍK = Av0-2<?^ + <f v0 Лу, = 0 • dz а0

Ядро К = в системе интегральных уравнений (4) равно:

*€ Z

где #<'> — функция Ханкеля первого рода,

(5)

и ветвь корня в выражении для Лк(^) выбрана так, что при 1тЯк(<^) Ф О выполнено неравенство 1т \ (£) > 0.

Множество тех £, для которых 1т Я^) = 0 равно множеству

Б = {//* + ¡к: к е Ъ,ц е < /а0 }и {у + ¡к: к е я}. (6)

1 Здесь 2я- — период. В случае, когда период волновода равен произвольному числу р>0 в качестве П следует взять множество ^х,у,г)': О <2<р].

и

f)

+ К

'и^

и.

(7)

При £ е С ряд (5) и х2 + у2 * 0 ряд (5) сходится и определяет непрерывно дифференцируемую по переменным x,y,z и 2тг - периодическую по переменной z функцию K{x,y,z\%), заданную при всех (х,у, z)T е R3, исключая счётное множество особых точек (0,0,2лк)т, к sZ. Функция K(x,y,z;£) и её производные по переменным х, у, z при любых фиксированных (x,y,z)' *(0,Q,2xk)T, fceZ— однозначные аналитические функции параметра £ на множестве комплексных чисел.

Система интегральных уравнений (4) представима в виде векторного уравнения

2(a(x,y,z)-a0)du0 ч a(x,y,z)-a0 дп у где К — компактный оператор в пространстве Ж2' (Q) х L2 (S). Кроме того, однородное уравнение (7) имеет только нулевое решение. Следовательно, по теореме Фредгольма [10] оно однозначно разрешимо для всех исключая дискретное множество точек.

Теорема 1 Точечный спектр <ур волновода состоит из изолированных точек, которые не сгущаются на конечном расстоянии. Непрерывный спектр сгс волновода равен множеству Е, описанному формулой (6).

Во второй главе вводится волноводный оператор А в пространстве Z^(n) как замыкание дифференциального оператора

О

. д А= — + dz

0

где у, =

дх' ду

WaVs) + 6 , заданного на множестве

0

D = \ 0,,v2)r: v, е IV22(Q)х fV22(П\Q), v2 е W\(£2)x W\(П\Q)

[v,] = 0,

a| — cos(M,;C) + — COS(M,,V) дх dy

+ v2 cos(«, z)

= 0,

öv,

vwUsvuL- &

, dz

где соб(и, х), соз(п,у), соб(п,г) — компоненты внешней нормали п к поверхности Г в точке (х,у,г)Т.

Оператор А обладает следующим свойством:

если векторная функция V е является решением уравнения

А К = 0, то её первая компонента является решением уравнения (1). В диссертации доказан следующий результат:

Теорема 2 Точечный ар(А) и непрерывный ас(А) спектры оператора А совпадают соответственно с точечным и непрерывным спектрами волновода. Остаточный спектр сгг(А) оператора А пуст. С оператором А связано индефинитное скалярное произведение

о Г ^ 1-1 о]

1 / 1 пч

и,У

сЮ.

в пространстве ¿^(П) , относительно которого оператор А антиэрмитов.

Кроме того, во второй главе были введены собственные функции непрерывного спектра (СФНС) как ограниченные 2л- периодические по

переменной г решения уравнения

/ /■

0

дг 1-1

1

-4

У= о

(8)

\У*№я)+ь о

с и найдено индефинитное скалярное произведение в пространстве обобщённых функций двух СФНС, отвечающих плоским волнам:

<1/(£Ш£)> =

-1

/л* (к, £ ф) В± (к, £ ср) а<р + — /''Щ^А^Ье, & +| +1У (к, Г], Э) +

1л

-е 4

8 т

Л(<?) Г 2 + 2)

<У(Яе(£ + т) + 1т+ пШу, - 5),

где Л1, В1, Ф, Т— асимптотические коэффициенты в первом члене асимптотического разложения СФНС при г - -^х2 + у2 -> +оо. Для нахождения СФНС используется сведение уравнения (8) к системе интегральных уравнений (4).

2 (,) — алгебраическое скалярное произведение в пространстве С2

7

В третьей главе вводятся проекторы Р^ на инвариантные подпространства оператора А, соответствующие отрезкам Д непрерывного спектра. Промежуток А здесь не содержит особых точек, то есть при всех £ е А система интегральных уравнений (4) однозначно разрешима в пространстве ^2'(0)х!2(5)-

Эти операторы определяются в виде ряда

кеЪ д

/еЛГ(Д) *

где Л^Д)— множество номеров к, при которых 1тЛД£) = 0 для всех ;

£ е А, {Уи(%)}кег ыг/(4)— набор СФНС, соответствующий решениям сис- ,

темы (4) с функцией и0, равной

где лг^еК, |дгЛ| = 1 -

Полнота этого набора СФНС следует из теоремы: Теорема 3 Предположим, что открытый промезкуток Д непрерывного спектра лежит на мнимой оси, и его замыкание А не содержит особых точек непрерывного спектра. Тогда спектр оператора А в инвариантных подпространствах РАЬ\{П) и (/ - /*Д)1^(П) равен соответственно Д и ст(Я)\Д.

Пусть открытые промежутки Д, и Д2 непрерывного спектра не лежат на мнимой оси, и Д2 ={-£ :£е А,}. И пусть замыкания промежутков Д, и Д2 не содержат особых точек непрерывного спектра. Тогда спектр оператора А на инвариантных подпространствах р& .

РА Ь\(П) и (1 - Р& - Р& )Ь22(П) равен соответственно А,, Д2 и

В четвёртой главе рассмотрен случай волновода с поглощением, то есть к коэффициенту Ь уравнения (1) добавлено чисто мнимое слагаемое ¡е, £>0.

В этом случае дискретный спектр оператора А является точечным множеством, не сгущающимся на конечном расстоянии, а непрерывный спектр равен множеству

ыг

В уравнении (А - V = 0 удобно сделать подстановку V = ,

где

1^0 г'/ ^ о ту

Т = и-1(\м\~1 +\}и, Г4 =;у-'(|//|ч и — преобразование Фурье

функций пространства ¿2(П) по переменным хну, //2)г— двой-

ственная по отношению к переменная при преобразовании Фурье

и. Операторы 5 и 5"' заданы на множествах Л(Т)х Ь2(П) и

^(П)хДП. П(Т) = |/ е 122(П):+ \}и/е 4(П)| и являются взаимно обратными, то есть их композиции 55"' и 5",5' являются ограниченными операторами, продолжающимися по непрерывности до тождественного оператора в пространстве ¿2(П).

Эта подстановка приводит к уравнению:

о,

где В = — замкнутый оператор в пространстве 12(П), заданный на плотном множестве. Спектр оператора В совпадает со спектром оператора А.

Для оператора В в пространстве 12(Г1) определены ограниченные

проекторы Р±, соответствующие частям спектра оператора В, лежащим слева и справа от мнимой оси. Для их определения вводятся вспомогательные операторы А, и В, по тем же формулам, что и операторы А и В,

но с заменой коэффициентов а и Ь операторами

a,^а0 + Е,(а-а0)Е1 = (1 - Е,)а0(Г - Е,) +Е,аЕ„

b, =Ь0 + Е,(Ь-Ь0)Е, =(/-£,)£„(/-£,) + где Е,— проектор в пространстве вида

ХЛМ) (0 при \/г\ > / Операторы Аг и В, являются замкнутыми операторами в пространстве 12(П), заданными на плотном множестве. Непрерывные спектры операторов А, и В, равны множествам

и{£еС:0<^)<*2}=

к Л

М|£еС:1т£ =----

ЫГ 2Яе£ 4(Ле#)2 а0 ]

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

точечными множествам, не имеющими на конечном расстоянии точек сгущения.

Пусть Г*,— расширяющиеся ориентированные ограниченные контуры, лежащие в резольвентном множестве оператора в, и охватывающие множества:

1. контур Г*, охватывает множество

к^-п

2. контур Г;, охватывает множество

€ С: < 0, Лк(0 е Я,0 < А(<?) ^ /}•

к=-п

Тогда формулы

2т л^+=сгх

* и

определяют ограниченные проекторы в пространстве ¿^(П).

Проекторы Р* имеют сильный предел в пространстве ¿^(П), пред-ставимый в виде:

+00

Р* = К + ¡{(В -¡Т]!)-1 -(Вй -/77/Г'К,

-ос

где

1 ^ 0 2/ ,

Ь,С1—ограниченные операторы в пространстве ¿2(П), I — тождественный оператор в пространстве ¿2(П), В0— частный случай оператора В, с коэффициентами а и Ь, тождественно равными постоянным а0 и

С помощью проекторов Р± определяются проекторы />±(г), ге [0,2л-] в пространстве ¿^(Я2):

Возможность такого определения следует из свойств локальности проекторов Р±:

р±{ф)0(х,у)) = ф){р±с)(х,у,г), <ре 1,(0,2^),Се В последнем параграфе четвёртой главы определяются разрешающие операторы задачи (1) - (3) по формулам:

™ . _

-оо

где операторы X„(z,s) равны e"(z~s)B" Р£ •

Для этих операторов справедливы теоремы:

Теорема 4 Операторы X±(z,s) являются ограниченными операторами в пространстве L2(R2). Их образы непрерывно зависят от переменных z и s в пространстве R2) и удовлетворяют равенству < X±(z,s) = Pt(z)X±(z,s) = Xt(z,s)Pt(s) = Pi(z)X±(z,s)Pt(s).

' Кроме того, для всех F е L\(R2)

i BX±(z,s)F = о3.

Теорема 5 Пусть U(z)— решение дифференциального уравнения \ BU(z) = 0 на промежутке (0,2;г), непрерывное в пространстве на

замкнутом промежутке [0,2яг]. Тогда функции U±(z) = P±(z)U(z) также являются решениями этого уравнения на промежутке (0,2;г). ; Решение U(z) однозначно определяется парой элементов

I P*(z)U(z) и P'(z)U(z) по формуле

U(z) = X* (z,0) Р* (0) U(0) + X-(z,2x)p-(2ж) 0{2л).

Список цитированной литературы

1. Интегральная оптика. Сборник под редакцией Т. Тамира. М. 1978. 344с.

2. Никольский В. В. Электродинамика и распространение радиоволн. М. 1973. 608с.

3. Сайханов И.Б. Спектр диэлектрического волновода с периодической границей. Вестник Ленинградского университета. Сер. 1. 1982. Вып. 19. С. 103-104.

4. Детгузов В.И., Сайханов И.Б. Точечный спектр волоконного волновода в окрестности особых точек спектрального параметра. Вестник Ленинградского университета. Сер. 1. 1985. Вып. 8. С. 9 -12.

5. Дегрузов В. И. Инвариантные подпространства периодического волновода. Вестник ЛГУ. 1985. №8.

3 Здесь В — дифференциальное выражение, определяющее оператор В.

I 11

S

6. Дергузов В.И., Лобанова С.Р. Резольвента периодического волновода на непрерывном спектре. Проблемы мат. анализа. 1992. Вып. 13. С. 79-89.

7. Дергузов В.И. Нормировка собственных функций непрерывного спектра трёхмерного периодического световода. Проблемы мат анализа. 1997. Вып. 16. С. 68-87.

8. Дергузов В.И. Ортогонализация собственных функций трёхмерного световода. Вестник С.-Петербургского университета. Сер. 1. 1998. Вып. 1. С. 6-9.

9. Дергузов В.И. Интегральные уравнения трёхмерного световода. Вестник Санкт-Петербургского университета. Сер. 1. 1999. Вып. 4. С. 22-25.

10. Рид М., Саймон Б. Методы современной математической физики. М. Мир. 1977.

Список работ по теме диссертации

1. Дергузов В.И., Денисова И.В. Характеристический спектр приведённого волнового уравнения с периодическими коэффициентами в трёхмерном пространстве. Проблемы мат. анализа., 2000. Вып. 21. С. 110-137.

2. Дергузов В.И., Денисова И.В. Инвариантные подпространства трёхмерного периодического световода. Проблемы мат. анализа., 2001. Вып 23. С. 3-13.

3. Дергузов В.И., Денисова И.В. Спектр периодического световода в трёхмерном пространстве. Проблемы мат. анализа., 2002. Вып. 24. С. 103-131.

4. Дергузов В.И., Денисова И.В. Проекционные и разрешающие операторы трёхмерного периодического световода. Проблемы мат. анализа., 2003. Вып . 26. С. 87 - 139.

5. Денисова И.В. Разрешающие операторы трёхмерного периодического световода. Проблемы мат. анализа., 2003. Вып. 26. С. 67 -87.

!

1

ч

I

I ЛР№ 040815 от 22.05.97

Подписано к печати ттг. Формат бумаги 60X84 1/16. Бумага офсетная. Печать неизофафическая. Объем 1усл. п.л. Тираж 100 экз. Заказ<ЗР<(# Отпечатано в отделе оперативной полиграфии НИИХ СПбГУ I с оригинал-макета заказчика.

к 198504, Санкт-Петербург, Старый Петергоф, Университетский пр., 26.

1 I

1

•1

I (

P 15 2 7 5

/У

4 //

/

Введение.

Глава 1 Точечный и непрерывный спектры волново^да.

§ 1 Вывод интегральных уравнений.

§2 Условия разрешимости интегрального уравнения (1.1.18).

Глава 2 Волноводный оператор для приведённого волнового уравнения.

§ 1 Определение волноводного оператора.

§2 Собственные функции непрерывного спектра и их асимптотические представления.

§3 Индефинитное скалярное произведение собственных функций непрерывного спектра.

Глава 3 Исследование непрерывного спектра приведённого волнового уравнения.

§1 Ортогонализация собственных функций непрерывного спектра.

§2 Проекторы на инвариантные подпространства оператора А, соответствующие конечным отрезкам непрерывного спектра.

§3 Спектры оператора А на инвариантных подпространствах соответствующих конечным отрезкам непрерывного спектра (Формулировка теоремы).

§4 Изучение системы интегральных уравнений (3.3.2).

§5 Выбор решения системы (3.3.2).

§6 Зависимость функции у(£) от правой части уравнения (2.1.16).

§7 Доказательство принадлежности функции пространству Ь2(П). .;.:.

Глава 4 Проекционные и разрешающие операторы для волновода с поглощением./.

§ 1 Преобразование оператора А.80

§2 Проекторы для случая невозмущённого волновода.

§3 Аппроксимация операторов^ и В более простыми.

§4 Определение проекторов Р1 и ^(г).

§5 Построение разрешающих операторов.

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

Для некоторого класса задач этот спектр изучался в работах [2.1] - [2.9]. Данное исследование является их продолжением.

Диссертация посвящена изучению граничной задачи для приведённого волнового уравнения

Цу(дУи) + йи = 0 (1) в трёхмерном пространстве с 2л - периодическими по переменной г коэффициентами, имеющими разрыв на некоторой 2л - периодической по переменной г и ограниченной по переменным х и у на поверхности Г. Вне области, ограниченной поверхностью Г, коэффициенты а и Ъ предполагаются постоянными.

В качестве граничных условий для этого уравнения берутся условия ди и] = 0 и а— дп 0, (2) где [у] — скачок функции V на поверхности Г, —— производная функции и по внешней дп нормали к поверхности Г.

Задачи такого типа имеют много физических применений. Ниже приведены некоторые из них.

Физическое применение полученных результатов

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

Распространение волн в этом волноводе описывается системой уравнений Максвелла: frot E = -icojuH,

1 (3) rot H = icosE, где EYi.fi— соответственно диэлектрическая и магнитная проницаемости вещества волновода, со — частота распространяемой волны. Функции е и ¡л являются 2ж - периодическими по переменной z и равны постоянным е0 и ju0 вне замыкания области Р. Векторные функции Ей Н— электрическая и магнитная составляющие электромагнитного поля.

Система (3) рассматривается в областях РиШ3\Р .

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

Векторная функция Н удовлетворяет уравнению div (////) = 0, вытекающему из первого уравнения системы (1). Следовательно, существует такая достаточно гладкая функция А, что Н = /и~] rot А.

Подставляя это равенство в первое уравнение системы (3), получаем, что rot(£' + icoA) = 0. Следовательно, Е = -icoA -Vq>, где (р— дифференцируемая скалярная функция.

Назовём А векторным потенциалом электромагнитного поля, а (р— скалярным потенциалом электромагнитного поля.

Подставим полученные выражения для Е и Н во второе уравнение системы (3).

Из этого уравнения и равенства div (sE) = 0 получим следующие уравнения для скалярного и векторного потенциалов Аир: di\(EV(p) + ia>div(£A) = 0, (4) rot(/T' rotA) = (o2sA-io)£V<р. (5)

Уравнения (4) и (5) содержат в себе обе неизвестные величины А и (р. Чтобы этого избежать наложим на скалярный и векторный потенциалы условие div (еА) = С(р, (6) где С — некоторая функция, которая будет определена ниже. Такое условие не ограничивает общность решения системы уравнений Максвелла, так как векторный потенциал А определён с точностью до произвольного слагаемого вида Уу/, где у/— произвольная достаточно гладкая функция.

При выполнении условия (6) уравнения (4) и (5) приводятся к виду: div{EW(p) + icL>Сер = 0 , (7) rot (/Г1 rot А) = со2еА - io)£4(C~x div (еА)). (8)

Уравнение (8) является уравнением в частных производных второго порядка со старшими членами

Г1 [V(div А)-АА] + ico£2C~'V(div А). При С = -ia>e2 /л они сводятся к выражению АА . Уравнение (7) при этом становится таким: div(£ Vq>) + со2е2/л<р = 0. Полученные уравнения для скалярного и векторного потенциалов в случае постоянных е и ju совпадают с классическими. Условия для скалярного и векторного потенциалов описаны, например, в работе [1.13].

Условие (4) в постоянной среде совпадает с традиционным.

При учёте инерционности процессов намагничивания магнитная проницаемость ¡л среды становится комплекснозначной функцией, и рассматриваемый волновод становится волноводом с поглощением. В диссертации рассмотрен случай, когда чисто мнимая часть коэффициента а>2е2/л постоянна. Общий случай требует дальнейшего исследования.

2. Пусть, также как и раньше, область Р является волокном диэлектрического In - периодического волновода и параметрами, описываемыми 2л- периодическими по переменной z функциями.

Свободные колебания электромагнитных волн в этом волноводе описывается системой

3).

Будем искать решение уравнений (3) в виде: е = е0 + уф(е),

9) н = н0 + уф(н), где Е0 и Н0 — векторные функции, равные нулю на границе области Р и при z = 2лк, к еЖ и удовлетворяющие в областях Р и Iß3 \ Р условиям: div(f£o) = 0, div (juH0) = 0. (10)

Такой способ решения предложен в работе [1.14].

Несложно убедиться, что пары функций Е0, УФ(Е) и Н0, УФ(Н) ортогональны в области Q = {(Ху, z)T е Р: 0 < z < 2яг| с весами е и /л:

J*£0VO(E)c/Q= JVE)div(££0)6fQ =0, n n j// я0 уф(н)<ю = |ф(н) div(juH0)dQ. = 0, n n а также в областях Qk, получающихся сдвигом области Q вдоль оси z на 2лк, к е Ж.

Более того, любая функция F, принадлежащая пространству L\ E(Q.) векторных функций, квадратично суммируемых с весом е, единственным образом представима в виде суммы F = Fx + F2, где функция F, принадлежит замыканию в пространстве Ú2 е (Q) множества к е Cq(Q) : div(sV) = о}, a F2 — множеству v<p:<pel2/oc(Q),V<peL32( Q)}.

Это доказывается также, как аналогичный результат для случая £• = 0 в работе [1.15]. При этом попутно доказывается, что если F е W\ (Q), то F2 е W\ (Q).

То же самое верно и для функций, заданных в областях П \ Q, где П = :0<z <27г).

Таким образом, любое решение Е, Я системы уравнений (3) представимо в виде (9). Для решений Е и Я системы уравнений (1) выполняются равенства div(¿r£) = 0, div(>tf) = 0.

Следовательно, согласно формуле (8), функции Ф(Е) и Ф(Н) должны являться решениями уравнений div(¿-V<D(E)) = 0, div(/¿V<D(H)) = 0. (11)

Пусть теперь функции Ф(Е) и Ф(Н) являются классическими решениями уравнений (11). Построим по ним решение системы уравнений (3) вида (9).

Для того чтобы векторные функции £иЯв формуле (9) были решением системы уравнений (3), необходимо и достаточно, чтобы векторные функции Е0 и Н0 удовлетворяли условию: rot Е0 = -iú)/uH0 - ку//УФ(Н), ^ rot Н0 = icoeE0 + icosV Ф(Е).

Поделим первое уравнение системы (12) на функцию ¡л и применим к обеим частям rot.1 Тогда, подставляя в правую часть выражения для rot Н0, полученное из второго уравнения системы (12), получим следующее выражение для векторной функции Е0: -т/л [iú)eE0 + icos УФ(Е) ]- ico/л УФ( rot rot Е0 М н) У

1 1 --АЕ0+ — VdivE0+V И 1

KMJ хrotЕ0 -со1 /леЕ^ =<v2peVФ(Е) -/¿у//УФ(Н). (13)

Это можно сделать, так как Е,Н Е W\ (Q) П W2 (П \ Q) , и, следовательно, ф 6 W2 (Q) П W2 (П \ Q)

Верно и обратное: если векторная функция Е0 является решением уравнения (11), то векторные функции Е0 и Н0 = -(/¿у//)-1 гог^ - УФ(Н) удовлетворяют системе уравнений (12),. оо .

Разрешимость уравнения (11) в пространстве \У) (О) гл (П \ О) была изучена в работе [1.14]. Таким образом, по всякому решению уравнений (11) строится решение системы (3), причём составляющая (УФ(Е), УФ(Н))г в разложении (9) этого решения есть составляющая, переносящая энергию вдоль волновода. •

3. Уравнения рассматриваемого вида можно использовать при моделировании упорядоченных полупроводниковых структур на поверхности твёрдого тела. Такие задачи возникают в современных технологиях получения наноструктур2. Более подробно эти задачи описаны в работах [1.16]-[1.17].

Характеристика работы по главам

В первой главе исследуется разрешимость поставленной граничной задачи в классе функций вида Флоке, то есть представимых в виде и = е'^у, где комплексное число, а функция V 2л - периодична по переменной г.

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

Дискретный спектр состоит из тех точек , для которых функция V квадратично суммируема на полосе П = г)т: 0 < г < 2л], а непрерывный спектр — это замыкание множества точек £, для которых функция V ограничена, но не квадратично суммируема на полосе П.

Доказан следующий результат.

Теорема Спектр волновода является дискретным множеством точек, не сгущающихся на конечном расстоянии, а непрерывный спектр совпадает с множеством: + Ис:кеЖ, ц е | /л |< / а0 }и {у + ¿к : к е Ъ, V е }.

Эти результаты похожи на результаты, полученные в работах [2.1], [2.2] и [2.9] для другого уравнения. Изложенные здесь результаты опубликованы в работе [2.10].

Во второй главе вводится волноводный оператор А, действующий в пространстве Ь] (П), непрерывный и дискретный спектры которого совпадает с непрерывным и дискретным спектрами волновода.

2 Наноструктуры это сверхпроводящие полупроводниковые структуры атомных размеров.

Оператор А определяется дифференциальным выражением А, обладающим следующим свойством: если векторная функция V е Ь22 (П) является решением уравнения А V = 0, то её первая компонента является решением уравнения (1).

Кроме того, во второй главе определяются собственные функции непрерывного спектра оператора А как 1п - периодические по переменной 2 ограниченные решения уравнения (А-£ 1)У = 0 , а также вводится естественное для оператора А индефинитное скалярное произведение собственных функций непрерывного спектра в пространстве обобщённых функций.

Результаты этой главы похожи на результаты работы [2.7] и опубликованы в работе [2.10].

Третья глава диссертации посвящена проекторам на инвариантные для оператора А подпространства пространства Ь\ (П), отвечающие конечным отрезкам непрерывного спектра волновода. Основной результат третьей главы сформулирован в теореме 3.3.1:

Теорема 3.3.1 Предположим, что открытый промежуток А непрерывного спектра лежит на мнимой оси, и его замыкание А не содержит особых точек непрерывного спектра, описанных в замечании 2.3.1. Тогда спектр оператора А в инвариантных подпространствах Р&Ь22(П) и (I-Ра)Ь22(П) равен соответственно А и сг(Л)\ А.

Пусть открытые промежутки А, и Д2 непрерывного спектра не лежат на мнимой оси, и Д2 = {—¿г : £ е Д,}. И пусть замыкания промежутков А, и Д2 не содержат особых точек непрерывного спектра, описанных в замечании 2.3.1. Тогда спектр оператора А на инвариантных подпространствах Р^ 1}2 (П), Р^ 1}2 (П) и (/ - Р^ — Р&)Ь\( П) равен соответственно А,, Д2 и сг(Л)\(Д, иА2).

Этот результат опубликован в работе [2.12].

В четвёртой главе рассматривается случай волновода с поглощением, то есть исследуется уравнение (1) с коэффициентом Ь, имеющим вид суммы вещественнозначной неотрицательной функции и чисто мнимого слагаемого ¿е, где £— ненулевое вещественное число. Результаты предыдущих глав легко переносятся на этот случай.

Непрерывный спектр волновода здесь имеет вид и|^е€ : 2Яе£(1т £-к) = -е, (Яе£)2-(1т£-к)2 > 0 кеЪ I а0

3 Тогда векторная функция II = е — решение уравнения А V = 0 вида Флоке.

В четвёртой главе дополнительно предполагается, что коэффициент а в уравнении (1) непрерывная функция во всём пространстве 0&3. Это предположение носит технический характер и упрощает некоторые доказательства. В общем случае аналогичных результатов пока не получено.

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

Точечный и непрерывный спектры операторов А и В совпадают.

Для оператора 5 определяются ограниченные проекторы Р± на инвариантные подпространства оператора В, соответствующие частям спектра оператора В, лежащим соответственно слева и справа от мнимой оси. Сама мнимая ось не содержит точек спектра оператора В.

Кроме того, определяются ограниченные проекторы Р±(г) в пространстве 1?2 (Н&2) как композиции операторов Р± , заданных в пространстве (К2), и оператора сужения функции, заданной в пространстве К , на плоскость [К х {г}. Проекторы Р (г) следующим образом связаны с введённым индефинитным скалярным произведением:

Следствие 4.5.1 Сужение индефинитного скалярного произведения (,) на подпространство

Р+(г)Ь\ (0&2) для всех г е [0,2я] является скалярным произведением в гильбертовом пространстве.

Точно так же и -(,) является скалярным произведением на подпространстве Р~ (г) (К2) для всех г б [0,2я].

В последнем параграфе четвёртой главы определяются разрешающие операторы в пространстве /^(П) для уравнения (1). Это означает, что первые компоненты их образов являются решениями уравнения (1) вида Флоке. Полнота этой системы решений следует из теоремы:

Теорема 4.5.2 Пусть и(г) — решение дифференциального уравнения В и {г) = 0, где В дифференциальное выражение, определяющее оператор В, на промежутке (0,2л-), непрерывное в пространстве Ь22 (¡К2) на замкнутом промежутке [0,2л-]. Тогда функции 11± {г~) = Р± {г) и (г) также являются решениями этого уравнения на промежутке (0,2л-).

Решение и (г) однозначно определяется парой элементов Р+( 0) С/ (0) и Р (2тг)£/(2тг) по формуле и {г) = Х+ (г,0) Р+ (0) 17(0) + Х~ (г,2я) Р' (2 л) Щ2я).

Результаты четвёртой главы опубликованы в работах [2.13] - [2.14].

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

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

1. Михлин С.Г. Линейные уравнения в частных производных. М.: Высшая школа. 1977.

2. Коренев Б.Г. Введение в теорию бесселевых функций. М. 1971.

3. Гельфанд И.М., Шилов Г.Е. Обобщённые функции и действия над ними. М.: Физмат-гиз, 1959.

4. Ладыженская O.A., Уральцева H.H. Линейные и квазилинейные уравнения эллиптического типа. М. Наука. 1964

5. Стейн И. Сингулярные интегралы и дифференциальные свойства функций. М. Мир, 1973, 342 стр.

6. Рид М., Саймон Б. Методы современной математической физики, М. Мир. 1977.

7. Федорюк М.В. Метод перевала, М., 1977, 368 стр.

8. Владимиров B.C. Уравнения математической физики. М. 1971.

9. Бирман М.Ш., Соломяк М.З. Спектральная теория самосопряжённых операторов в гильбертовом пространстве. Л. 1980. С. 264.

10. Лионе Ж.Л., Мадженес Э. Неоднородные граничные задачи и их приложения. М. 1977. С. 371.

11. Хилле Э., Р. Филлипс Р. Функциональный анализ и полугруппы. М. 1962. С. 431.

12. Интегральная оптика. Сборник под редакцией Т.Тамира. М. 1978.

13. Тамм И.Е. Основы теории электричества. М. 1956.

14. Вайнштецн Л.А. Электромагнитные волны. М. 1957.

15. Ладыженская O.A. Математические вопросы динамики вязкой несжимаемой жидкости. М. 1961.

16. Дубровский Г.В., Бауман Д.А., Дубровский В.Г., Козачёк В.В., Мареев В.В., Марков Ю.Г. Детальная кинетика многослойной адсорбции. СПб. 1998.

17. Марков Ю.Г. Эвристические подходы в моделировании сложных систем. СПб. 2002.

18. Никольский В.В. Электродинамика и распространение радиоволн. М. 1973.2. Статьи.

19. Сайханов И.Б. Спектр диэлектрического волновода с периодической границей // Вестник Ленинградского университета, Сер. 1.1982. Вып. 19. С. 103-104.

20. Дергузов В.И., Сайханов И.Б. Точечный спектр волоконного волновода в окрестности особых точек спектрального параметра // Вестник Ленинградского университета, Сер. 1.1985. Вып. 8. С. 9-12.

21. Дергузов В.И. Инвариантные подпространства периодического световода // Вестник ЛГУ. 1985. №8.

22. Дергузов В.И., Лобанова С.Р. Резольвента периодического световода на непрерывном спектре. // Проблемы мат. анализа. 1992. Вып. 13. С. 79-89.

23. Дергузов В.И. Решение некорректной задачи для одного класса линейных уравнений с периодическими коэффициентами // Проблемы мат. анализа. 1972. Вып. 3. С. 3-28.

24. Дергузов В.И. Операторные гамильтоновы и антиканонические уравнения с периодическими коэффициентами // Проблемы мат. анализа. 1973. Вып. 4. С. 9-36.

25. Дергузов В.И. Нормировка собственных функций непрерывного спектра трёхмерного периодического световода // Проблемы мат анализа. 1997. Вып. 16. С. 68-87.

26. Дергузов В.И. Ортогонализация собственных функций трёхмерного световода // Вестник С.-Петербургского университета. Сер. 1. 1998. Вып. 1. С. 6-9.

27. Дергузов В.И. Интегральные уравнения трёхмерного световода. // Вестник С.-Петербургского университета. Сер. 1. 1999. Вып. 4. С. 22-25.

28. Дергузов В.И., Денисова И.В. Характеристический спектр приведённого волнового уравнения с периодическими коэффициентами в трёхмерном пространстве. // Проблемы мат. анализа., 2000. Вып. 21. С. 110-137.

29. Дергузов В.И., Денисова И.В. Инвариантные подпространства трёхмерного периодического световода. // Проблемы мат. анализа., 2001. Вып 23. С. 3—13.

30. Дергузов В.И., Денисова И.В. Спектр периодического световода в трёхмерном пространстве. // Проблемы мат. анализа., 2002. Вып. 24. С. 103—131.

31. Дергузов В.И., Денисова И.В. Проекционные и разрешающие операторы трёхмерного периодического световода. // Проблемы мат. анализа., 2003. Вып . 25.

32. Денисова И.В. Разрешающие операторы трёхмерного периодического световода. // Проблемы мат. анализа., 2003. Вып. 25.