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

0 Введение

1 Глава 1. Нелинейные краевые задачи на собственные значения

1.1 Постановка краевой задачи на собственные значения для системы уравнений Максвелла (ТЕ-поляризация)

1.2 Сведение к нелинейной краевой задаче на собственные значения для дифференциальных уравнений (ТЕ-поляризация)

1.3 Постановка краевой задачи на собственные значения для системы уравнений Максвелла (ТН-поляризация)

1.4 Сведение к нелинейной краевой задаче на собственные значения для дифференциальных уравнений (ТН-поляризация)

2 Глава 2.Исследование разрешимости задач на собственные значения

2.1 Функция Грина и ее свойства.

2.2 Сведение краевой задачи к нелинейным интегральным уравнениям.

2.3 Теоремы о существовании и единственности решений интегральных уравнений.

2.4 Теорема о непрерывной зависимости решения от спектрального параметра.

2.5 Теоремы о существовании решений дисперсионного уравнения и задачи на собственные значения.

2.6 Формулировка итерационного метода решения интегрального уравнения.

2.7 Теорема о сходимости итерационного метода.

2.8 О некоторых оценках параметров.

3 Глава З.Комплекс программ и результатов расчетов

3.1 Алгоритм решения задачи на собственные значения в нулевом и первом приближениях.

3.2 Алгоритм полного решения задачи на собственные знаг чения.

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

Распространение электромагнитных волн в волноводах с заполнением линейной средой (то есть когда диэлектрическая и магнитная проницаемости не зависят от электромагнитного поля) -тема классической электродинамики [9], [28], [29], [30].

В случае волновода кругового сечения и постоянных электрической и магнитной проницаемостей уравнения Максвелла решаются в цилиндрических координатах, при использовании метода разделения переменных появляется линейное обыкновенное дифференциальное уравнение, называемое уравнением цилиндрических функций или уравнением Бесселя, решение которого является комбинацией цилиндрических функций. Собственные функции и собственные значения определяются как решения краевых задач с дополнительными условиями на контуре для решений и их первых производных [9]. С появлением нелинейной оптики предметом изучения в электродинамике стали сильные волновые поля, в которых начинает проявляться нелинейность сред. Качественно новыми эффектами нелинейной оптики стали порождение средой высших гармоник, а также самовоздействие волнового процесса, распространяющегося в нелинейной среде [9]. При распространении через жидкость или газ волны, создаваемой лазером, учитывают нелинейность поляризации среды, вызываемую целым рядом факторов [9]. Помимо поведения электронов в сильном электромагнитном поле существенно механическое воздействие поля на вещество: возникает давление, пропорциональное средней мощности волны, в результате чего в областях сгущения увеличивается диэлектрическая проницаемость. При распространении резко неоднородной волны — луча лазера -можно сказать, увеличивается оптическая плотность среды в области сильного поля. Иными словами, в определенных условиях волновому процессу сопутствует образование канала, направляющего его энергию - нечто вроде диэлектрического волновода. Это называется самоканализацией, если канал сужается, наступает так называемая самофокусировка [9].

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

Квантовая теория самофокусировки была разработана Талановым для плазмы и группой американских ученых R.Y.Chiao, Е. Garmire, С.Н. Townes [22] для твердых тел в 1964 году, где исследуются условия, при которых электромагнитный луч может распространяться без затухания в средах с диэлектрической проницаемостью, возрастающей пропорционально интенсивности электромагнитного поля, и постоянной в его отсутствие. В работе [22] представлено теоретическое объяснение самофокусировки. Если при воздействии электромагнитного поля диэлектрическая постоянная обуславливает критический угол внутреннего отражения, превышающий дифракционный угол отклонения луча, то рассеяние путем дифракции не будет иметь места. Таким образом, при мощности луча (уменьшающейся пропорционально квадрату длины волны), выше некоторого критического уровня этот луч может быть сфокусирован при любом произвольном радиусе. В обычных диэлектрических материалах для достижения самофокусировки достаточно мощности обычного лазерного светового луча, т.е. порядка 108 Вт.

Авторами [22] получено точное решение волнового уравнения в форме линейно поляризованного плоского светового луча и в форме цилиндрически поляризованного луча кругового сечения с поперечной составляющей электромагнитного поля .

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

Расслоение среды на области, в одних из которых электромагнитное поле поперечно, а в других появляется также и продольная составляющая, было обнаружено также В.И.Елеонским и В.П.Силиным [20] при исследовании сходной задачи отражения наклонно падающих р - поляризованных волн от среды с нелинейной диэлектрической проницаемостью. Рассматривается отражение плоско-поляризованных волн (вектор электрической составляющей которых лежит в плоскости падения) от полупространства, заполненного средой с нелинейной диэлектрической проницаемостью. Чередование областей с поперечной и продольной составляющими поля объясняется чередованием продольной и поперечной степеней свободы электромагнитного поля в нелинейной среде.

Наряду с рассматриваемыми проблемами плоской геометрии, В.И. Елеонский представил теорию цилиндрических самофокусирующихся волноводов, в основном базирующуюся на качественном анализе фазовых траекторий [25], а также численные расчеты основной и низшей неосновной ТЕ-мод.

Y. Chen, основываясь на вариационной технике, получил аналитические решения для фундаментальной моды в свободном пространстве в нелинейных средах, используя аппроксимацию Гауссовой функцией [17], [24].

R.A.Sammut, C.Pask. [15] использовали вариационную формулировку волнового уравнения для волноводов с произвольной нелинейностью оптического волокна и представили аналитическую аппроксимацию для поля.

Позже D.Sjoberg [23] анализировал распространение электромагнитных волн в нелинейных цилиндрических волноводах методом возмущений, используя коэффициент нелинейности как параметр возмущения.

Для моделирования распространения волн в средах с нелинейностью Керра также используются методы распространяющегося луча [37], метод линий [38], конечно — разностные методы [8].

Наряду с вышеописанными разработками существуют также направления, посвященные проблемам устойчивости. Изучение устойчивости распространения волн необходимо при исследовании возможности практического применения эффекта самофокусировки. Важными в этой области являются исследования Н. Н. Ахмедиева [5], [6], [18], которые посвящены нестационарной проблеме нестабильности модуляции основной моды цилиндрического волновода с нелинейным заполнением среды по закону Керра. Эффект нестабильности модуляции в одномодовом оптическом волокне к этому времени был уже предсказан, исследован теоретически и практически наблюдаем [31] - [36]. Рассматривается распространение высоко- и низкочастотного поля в волноводе круглого сечения, диэлектрическая проницаемость среды внутри и вне волновода описывается законом Керра. В одномодовом оптическом волокне при нелинейных эффектах поперечная структура пучка определяется самим волноводом в линейном режиме, а продольная функция для относительно низких частот определена сочетанием нелинейности и дисперсии волокна. Для высоких частот или больших радиусов волновода должна также быть учтена дифракция, и ее влияние при высоких частотах на нестабильность модуляции исследовано автором. Эффект нестабильности модуляции вызывает проявление отдельных световых групп, однако, дальнейшее поведение этих групп является открытым вопросом. Наиболее их вероятное поведение в среде с нелинейностью по закону Керра - это пространственно - временное затухание или коллапс. Таким образом, в результате нестабильности модуляции на выходе оптического волокна имеем затухание в пространстве и времени вибраций. В среде с насыщенной нелинейностью, где также существуют подобные эффекты, световые группы при определенных условиях могут распространяться на достаточно длинные дистанции в виде трехмерных стационарных солитонов.

Вопросу нелинейных волновых взаимодействий посвящена работа А.П. Сухорукова [60]. Распространение волн в нелинейной среде в слое изучала группа исследователей под руководством В.П. Шестопалова [54] — [59], сходной проблемой о плоском трехслойном волноводе занимались H.W. Shurmann, B.C. Серов и Ю.В. Шестопалов [1] — [4], в работах которых описываются отражение и распространение плоской ТЕ-поляризованной волны. Аналогично случаю, рассматриваемому в настоящей диссертационной работе, ищутся решения уравнений Максвелла, удовлетворяющие условиям непрерывности на границе раздела сред. При этом также используется метод разделения переменных.

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

Изучая распространение волн в слое, разные авторы использовали различные подходы, вводя специальные условия по отношению к нелинейной системе Фабри - Перота: J.H.Marburger и

F.S.Felber [39] упростили анализ путем введения граничных условий, предполагающих, что нелинейный слой отделен от линейных сред идеальными зеркалами; J.Danikaert [40] рассматривает реакцию твердого резонатора Фабри — Перота, включающую нелинейное поглощение и наклонное падение поперечно — электрического и поперечно — магнитного полей; M.Haeltermann [41] и G.Vitrant [42] представили объединенную нелинейную теорию для поперечных эффектов резонаторов Фабри -Перота, упростив численные расчеты и хорошо объяснив понятие оптической бистабильности.

Решения нелинейных уравнений Гельмгольца были представлены частично, для специальных случаев) в виде эллиптических функций Якоби [43] — [47]. Как показано в работе одного из авторов [48], общее решение более эффективно может быть представлено в форме эллиптических функций Вейерштрасса. Хорошо известно [14], [27], что эллиптические функции, в общем случае, обладают полюсами, однако, в предыдущих исследованиях е общее решение может быть использовано непосредственно для определения критических значений коэффициента отражения в различных случаях. Также отсутствовал анализ неограниченных решений в случае отрицательной константы Керра, в частности, никаких условий для существования таких решений.

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

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

Проводя сравнительный анализ этой работы с предыдущими исследованиями [50], [45], [46], можно отметить, что построение диэлектрической проницаемости 8, процедура решения нелинейного уравнения Гельмгольца и формулировка граничных условий являются традиционными, в то время как подход представления решений в виде эллиптических функций Вейерштрасса вместо использующихся ранее Якобиановых является более эффективным и имеет явное преимущество, описывая более общий физический случай.

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

Работа посвящена решению нелинейной краевой задачи на собственные значения для нелинейного уравнения Гельмгольца, к постановке которой приводит вопрос об изучении распространения собственных волн в цилиндрических диэлектрических волноводах кругового сечения, заполненных средой с нелинейностью по закону Керра. Такого рода задачи так или иначе сводятся к изучению сложной оператор - функции, нелинейно зависящей от спектрального параметра, которое весьма затруднительно традиционными методами теории дифракции. В связи с этим перспективным является метод оператор - функций, впервые примененный к задачам электродинамики в работах А.С. Ильинского и Ю.В. Шестопалова [51] - [53] . Для сведения задачи о собственных волнах к задаче на собственные значения с нелинейным вхождением спектрального параметра для интегральных уравнений используется метод функции Грина. Впервые для доказательства существования решений краевой задачи на собственные значения используются методы, базирующиеся на теореме Шаудера и принципе сжимающих отображений. Следует отметить, что такие методы использовались ранее для решения интегральных уравнений [12], но не были использованы в целом при решении задачи на собственные значения. При использовании предлагаемого метода спетральный параметр X изначально не фиксируется, а для каждого его значения решается уравнение краевой задачи, а затем полученное решение подставляется в дисперсионное соотношение. Принцип сжимающих отображений применяется для доказательства существования решений интегральных уравнений. В качестве численного метода предлагается итерационый метод, сходимость которого теоретически обоснована в работе.

Итак, основными целями настоящей работы являются:

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

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

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

Работа состоит из трех глав и двух приложений.

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

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

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

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

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

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

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

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

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

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

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

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

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

4 Заключение

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

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

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

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

4. На основании анализа дисперсионного соотношения возможна продуктивная разработка анализа стабильности различных мод. Известно, что интеграл потока энергии Р = / dpp3z (р, 7) (где обозначает усрененную по времени z-компоненту ветора Пойтин-га) может быть использован в некоторых случаях для анализа стабильности. z может быть вычислено посредством итерационных решений un(7,/?).

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

1. Serov, V.S., Shestopalov, Yu.V., and Schiirmann, H.W. Existence of eigenwaves and solitary waves in lossy linear and lossless nonlinear layered waveguides. // Dokl. Maths. -1996. -P.98-100

2. Schiirmann, H.W., Serov, V.S., and Shestopalov, Yu.V. TE-polarized waves guided by a lossless nonlinear three-layer structure.Phys. Rev. E. 1998. P.1040-1050.

3. Serov, V.S., Shestopalov, Yu.V., and Schiirmann, H.W. Propagation of ТЕ waves through a layer having permittivity depending on the transverse coordinate and lying between two half-infinite nonlinear media. // Dokl. Maths. 1999. P.742-744.

4. Schiirmann, H.W., Serov, V.S., and Shestopalov, Yu.V. Solutions to the Helmholtz equation for TE-guided waves in a three-layer structure with Kerr-type nonlinearity. //J. Phys. A: Math. Gen. 2002. P.10789-10801.

5. Akhmediev N.N., Ankiewicz A. Solitons,Nonlinear Pulses and Beams. Chapman and Hall. London. - 1997.

6. Akhmediev N.N. Spatial solitons in Kerr and Kerr-like media. // Invited Paper,Optical and Quantum Electronics. 1998. V.30. - No 7-10. - P.535-570.

7. Davies J.R.,Mendonca J.T. Laser Propagation in Cylindrical Waveguides. // Physical Review E. 2002 E.66.046604

8. Romanova E.A., Melnikov L.A., Bekker E.V. Light Guiding in Optical Fibers with Kerr-like Nonlinearity. // Microwave and optical technology letters. 2002. V.30, No 3.

9. Никольский, B.B. Электродинамика и распространение радиоволн. М., Наука, 1978.

10. Snyder A. and Love J. Optical Waveguide Theory. Chapman and Hall. 1983 - London.

11. Stakgold, I. Green's Functions and Boundary Value Problems. -Wiley. 1998. New York.

12. Треногин B.A. Функциональный Анализ. M.: Наука, 1993.

13. Zeidler Е. Applied Functional Analysis. Springer, New York, Berlin, Heidelberg, 1997.

14. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров.- М.: Наука, 1968.

15. R.A.Sammut, C.Park. Gaussian and equivalent-step-index approximations for nonlinear waveguides. // Journal of the Optical Society of America B. 1991. V.8. - No.2. - P.395-402.

16. L.Novotny., C. Hafner. Light rpopagation in a cylindrical waveguide with a complex, metallic, dielectric function. // Physical Review E. 1995. V.50. - No.5. - P.4094-4106.

17. Yijiang Chen. Nonlinear fibers with arbitrary nonlinearity. // Journal of the Optical Society of America B. 1991. V.8. - No.ll. -P.2338-2341.

18. Akhmediev N.N., Korneev V.I., Nabiev R.F. Modulation instability of the lowest nonlinear mode of a cylindrical waveguide. // Physical Review A. 1991. V.46. - No.l. - P.430-435.

19. Eleonskii V.N., Silin V.P. Self-focusing of a vector field. // ZhETF Fis. Red. 1971. V.13. - No.3. - P.117-121.

20. Eleonskii V.N., Silin V.P. Nonlinear theory of penetration of polarized waves into a conductor. // Zh. Eksp. Teor. Fiz. 1971. -V.33. No.5. - P.1039-1044.

21. Herausgegeben von S.Fliigge. Elektrische felder und wellen. -Springer-verlag. 1958. Berlin.

22. R.Y.Chiao., E. Garmire., C.H. Townes. Self-trapping of optical beams.// Physical Review Letters. 1964. V.13. - No.15. - P.479-482.

23. Daniel Sjoberg. Guided waves in nonlinear media //

24. Yijiang Chen. ТЕ family of self-guided beams in saturable nonlinear media. // Journal of lightwave technology. V.9. - No.9. — P.1208-1213.

25. Eleonskii V.N., Oganes'yants L.G., Silin V.P. Cylindrical nonlinear waveguides. // Zh. Eksp. Teor. Fis. Red. 1972. V.35. - No.l. -P.44-47.

26. Abakarov D.I., Akopyan A.A., Pekar S.I. Contribution to the theory of self focusing of light in nonlinearly polarizing media. // Zh. Eksp. Teor. Fis. Red. 1967.- V.25. No.2. - P.303-305.

27. Зайцев, В.Ф., Полянин, А.Д. Справочник по нелинейным дифференциальным уравнениям, М., Факториал, 1997г.

28. Bartling, J.Q. Propagation of an electromagnetic wave in an infinite rectangular dielectric waveguide. // Journal of the Franclin . May 1969.- V.287. No.5. - P.303-305.

29. Borgnis, F., Papas,Stratton, J.A. Electromagnetic Theory C. Electromagnetic Waveguides. // Handbuch der Physic 16, Editor Flugge.

30. Stratton, J.A. Electromagnetic Theory.// Mc Graw Hill, New York 1941.

31. Hasegawa, A., Brinkman, W.F. IEEE J. // Quantum Electron, QE-16, 694 (1980)

32. Akhmediev, N.N., Eleonskii, V.M., Kulagin, N.E. // Zh. Eksper. Teor. Fiz.,89, 1542 (1985) Sov. Phys. JETP 62,894 (1985)].

33. Tai, K., Hasegawa, A., Tomita, A. // Phys. Rev. Lett. 56,132 (1986).

34. Tai, K., Tomita, A., Jewell, J.L., Hasegawa, A. // Appl. Phys. Lett. 49, 236 (1986).

35. Hasegawa, A., Kodama, Yu. // Proc IEEE 69, 1145 (1981).

36. Mollenauer, L.F., Stolen, R.H. // Laser Focus 18, 193 (1982).

37. Gonthier,, F., Henault, A., Lacroix, S., Black, R.J., Bures, J. Mode copling in nonuniform fibers: comparison between coupled mode theory and finite-difference beam-propagation method simulation. // J. Opt. Soc. Amer. В 8 (1991), P.416-421.

38. Bertolotti, M., Masciulli. P., Sibilia, C. Method of lines numercal analysis of nonlinear plener waveguide. //J Lightwave Technol 12 (1994), P.784-789.

39. Marburger, J. H., Felber, F.S. Theory of lossless nonlinear Fabry -Perot interferometer. // Phys. Rev. A 17 (1978) P.335.

40. Danikaert, J., Thienpont, H., Veretennikoff, I. Self-consistent stationary description of a nonlineary Fabry Perot. // Opt. Commun. 71 (1989) P.5.

41. Haeltermann, M. Transverse effects in nonlinear planar resonators.I. // Modal Theory, J. Opt. Soc. Am. В 7 (1990) P.1309.

42. Vitrant, G., Haeltermann, M., Reinisch, R. Transverse effects in nonlinear planar resonators. II. Modal analysis for normal and oblique incidence. // J. Opt. Soc. Am. В 7 (1990) P.1319.

43. Chen, W., Mills, D.L. Optical response of a nonlinear dielectric film. // Phys. Rev., В 35 (1986) P.524 532.

44. Leung, K.M. Scattering of transverse electric electromagnetic waves with a finite nonlinear film. // J. Opt. Soc. Am. В 5(2) (1988) P.571 - 574.

45. Chen, W., Mills, D.L., Optical behaviour of a nonlinear thin film with oblique S polarized incident wave. // Phys. Rev. В 38 (1988) P.12814 - 12822.

46. Leung, K.M. Exact results for the scaterring of electromagnetic waves with a nonlinear film. // Phys. Rev. В 39 (1988) P.3590 3598.

47. Peschel, Th., Dannberg, P., Lederer, F. Investigation of optical tunneling through nonlinear films. //J. Opt. Soc. Am. В 5 (1) (1988) P.29 36.

48. Schurmann, H.W., Schmoldt, R. On the theory of reflectivity and trassmissivity of a lossless nonlinear dielectric slab. //Z. Phys. В 92 (1993) P.179 -186.

49. Born, M., Wolf, E. Principles of optics. 3rd Edition, Pergamon Press, Oxford, 1965, P. 61-66.

50. Chen, W., Maradudin, A.A. S polarized guided and surface electromagnetic waves supported by a nonlinear dielectric film. // J. Opt. Soc. Am. В 5 (1988) P. 529 - 538.

51. Ильинский, А.С., Шестопалов, Ю.В. О спектре нормальных волн щелевой линии передачи, Радиотехника и электроника, 1981, т.26, 10, с. 2064 2073.

52. Ильинский, А.С., Шестопалов Ю.В. Математическая модель для задачи распространения волн в микроволновых устройствах, Вып. 32. М.: Изд во МГУ, 1980, с. 85 - 103.

53. Шестопалов, Ю.В. Существование дискретного спектра нормальных волн микрополосковых линии передачи со слоистым диэлектрическим заполнением. ДАН СССР, 1983, т. 273, 3, с. 594 596.

54. Виноградова, М.Б., Руденко, О.В., Сухоруков, А.П. Теория волн. Москва, Наука, 1990.

55. Akhmediev, N.N., Ankiewicz, A. Solitons. Nonlinear Pulses and Beams. London Weinheim - New York - Tokyo - Melbourne -Madras, Chapman and Hall, 1997.

56. Yatsyk, V.V. Own Modes Nonlinear Dielectric of a Layer and Algorithm of the Decision of Nonlinear Spectral Problems.// Dopovidi Nacionalnoi Akademii Nauk Ukrainy, no.11,2000, pp. 114 117.

57. Yatsyk, V.V. Numericallly Analytical Methods of Solution of the Diffraction Problem on the Nonlinear Dielectric Layer.// Computational Physics, no. 12,2003, pp. 1-6.

58. Yatsuk, V.V. Calculation of Own Oscillation Wave Modes Nonlinear Spectral Problems.// Europian Symposium on Numerical Methods in Electromagnetics, JEE"02, Toulouse, France, 6-8 March 2002, part 2, session CW, pp. 1-6.

59. Shestopalov, V.P., Yatsyk, V.V., Equations of Dispersional and Space Time Evolutions of Field of Quasiuniform Electrodynamic Structure Near Morse Critical Point.// Radiotechnica i Electronica, vol,45, no. 2, 2000, pp.157 - 164.

60. Сухоруков, А.П., Нелинейные волновые взаимодействия в оптике и радиофизике, М., Наука, 1975.

61. Ивлиева, С.Н. О точках дискретности спектра в нелинейной задаче о собственных волнах цилиндрического волновода. // Материалы ХПМеждународной школы семинара " Синтез и сложность упраляющих систем", Пенза, 15-21 октября, 2001. - С.104 - 109.

62. Ivleeva S.N., Smirnov Yu.G. Electromagnetic waves guided by a lossless nonlinear open structure.// Proceedings of European Symposium on Numerical Methods in Electromagnetics, 6-8 March 2002, Toulouse, France. P. 94 -97.

63. Ивлиева, С.Н. О некоторых решениях солитонного типа в дискретных моделях для волноведущих структур.// Материалы XIII Международной школы семинара " Синтез и сложность уп-раляющих систем", Пенза, 14-20 октября, 2002. - С.100 - 104.

64. Смирнов, Ю.Г., Куприянова, С.Н. Рапространение электромагнитных волн в цилиндрических волноводах, заполненных нелинейных средой. // Журнал вычислительной математики и математической физики — выпуск 10, 2004.

65. Куприянова, С.Н. Собственные электромагнитные волны в диэлектрических волноводах с нелинейной средой. // Труды международного юбилейного симпозиума (АПНО 2003), Пенза, 19-22 ноября, 200, т.1 - с. 30 -32.