CC BY
64
ЧИСЛЕННЫЕ МЕТОДЫ КОМПЬЮТЕРНОЙ ОПТИКИ
ПРИМЕНЕНИЕ ГАМИЛЬТОНОВА ФОРМАЛИЗМА ДЛЯ РЕШЕНИЯ НЕЛИНЕЙНОГО УРАВНЕНИЯ ВОЛОКОННОЙ ОПТИКИ
Алименков И.В.
Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет) (СГАУ)
Аннотация
Применён гамильтонов подход для решения расширенного уравнения распространения оптических импульсов в волоконных световодах для произвольных функций отклика нелинейной среды на внешнее гармоническое возмущение. Показано, что солитонное решение является функцией гиперболического косинуса при любой нелинейности.
Ключевые слова: волоконный световод, расширенное уравнение распространения, гамильтоновы системы, произвольная нелинейность, солитонное решение.
Введение
Поле оптического импульса, распространяющегося в одномодовом волоконном световоде, поддерживающем состояние линейной поляризации, имеет вид [1].
Е(г,0 = в^(х,у)А(г,Оехр{/-wot)} , (1)
где Е(х,у) - Гауссова функция вида ехр{-(х2+у2)/м>2} с радиусом моды V, А(г, ^ - комплексная огибающая импульса, юо - несущая частота, Ро = юоп(юо)/ с -центральное волновое число.
Для огибающей оптического импульса выведено [2] уравнение
' дА р аА| д^А удг 1 дt Г 2Ро дг2 др( Л\2) А = 0,
Р^ дА
2 а^
(2)
названное расширенным уравнением распространения, которое существенно отличается от традиционного, называемого основным уравнением распространения [1], наличием второй производной по координате. Здесь Р; = 1/ vg - величина, обратная групповой скорости, р2 - дисперсия групповой скорости, ДР(| А|2) - нелинейная поправка к постоянной распространения моды в линейном приближении, которая выражается через нелинейную часть Ап показателя преломления:
др =
ko ЦДп^ (х, у)|2ах ду Е (х, y)|2dx &у
(3)
где ko =юо / с.
В области прозрачности волновода Ар является вещественной функцией. Так как р2 в области прозрачности отрицательно, то тип уравнения (2) -эллиптический. Если в уравнении (2) отбросить одну из вторых производных, то полученное уравнение параболического типа сводится к нелинейному уравнению Шрёдингера, которое является на сегодняшний
день самым изученным из нелинейных уравнений, допускающих солитонные решения.
Полные уравнения второго порядка рассматривались и ранее. Так, в [3] выведено уравнение гиперболического типа для неограниченной керровской среды в упрощённой модели поля, не зависящего от поперечных координат, но затем это уравнение упрощается до нелинейного уравнения Шрёдингера. Следует упомянуть и статью [4], в которой методом секанса и tanh-методом решаются некоторые полные уравнения второго и даже третьего порядка, однако полученные решения сингулярны.
При небольших пиковых значениях интенсивности вводимого излучения нелинейную часть показателя преломления представляют степенным рядом
Дп = п2 Е + п. Е Н—,
(4)
что согласно (3) приводит к степенному разложению функции Др:
ДР = у| А|2 +ц| А|4
(5)
где параметры у и д зависят от характеристик световода. При достаточно малых интенсивностях оптического поля вторым слагаемым в (5) пренебрегают.
Полученная нелинейность называется керровской и солитонное решение уравнения (2) для этого случая найдено в [5].
С ростом интенсивности вводимого излучения наблюдается [6] отклонение от керровской зависимости показателя преломления и необходимо учитывать второе слагаемое в формуле (5). Такая нелинейность называется конкурирующей (иногда этот термин применяют только в случае, если слагаемые в (5) различаются знаком).
Экспериментальные исследования в нелинейной оптике [6] подтверждают такую зависимость нелинейного показателя преломления от интенсивности оптического поля в полупроводниковых волноводах, стёклах, допированных полупроводниками, и органических полимерах. Точное солитонное решение уравнения (2) для случая конкурирующей нелинейности найдено в [7].
Целью настоящей работы является сведение уравнения (2) к нормальной системе гамильтоновых уравнений для оптической интенсивности поля импульса в случае произвольной функции нелинейного отклика среды на внешнее гармоническое возмущение. Теория гамильтоновых систем со своим мощным аппаратом производящих функций на сегодняшний день, несомненно, является самым завершённым разделом математической физики.
Основной формализм От уравнения (2) для комплексной огибающей импульса перейдём к уравнению для оптической интенсивности I = |А|2 /2 с помощью подстановки
А(z, 0 = -¡21 ехр {ibz},
(6)
где Ь - произвольный параметр, являющийся поправкой к центральному волновому числу р0. Этот параметр всегда выражают через пиковое значение интенсивности или пиковое значение напряжённости.
Подстановка (6) в (2) приводит после отделения вещественной и мнимой частей полученного уравнения к следующим двум уравнениям для оптической интенсивности огибающей:
(
1+Ь Ро
51 1 Ы „
— +--= 0 ,
дz V„ дt
1
210 ^
Ь + ---ДВ(1)
-Р2
812.
21II - -
(7)
(8)
Таким образом, для одной неизвестной функции получено два уравнения. Уравнение (7) является линейным и однородным, что позволяет написать его общее решение, являющееся произвольной дифференцируемой функцией I = I (5(z, t)), где
5(2, ^ = 2 - 20 - Vt ,
V = V (1 + Ь / рв).
(9) (10)
Здесь учтено, что Р; = 1/ vg. Иными словами,
уравнение (7) определяет аргумент искомой функции, оставляя её вид произвольным. Таким образом, искомая функция представляет собой бегущую волну неизменного профиля, движущуюся с постоянной скоростью (10).
Профиль этой волны определяется уравнением (8), которое превращается в обыкновенное дифференциальное уравнение
1 -Р^2
Ро
21 £ - -
Ь + ---ДР(1)
2Р„ И
(11)
812.
Для дальнейшего удобства перейдём в (11) к без-
размерной независимой переменной т:
X = 5
Ь(Ь + 2Ро)
1 -Р„Р2V2
Тогда уравнение (11) принимает вид 211 -12 = 412 - 8у!2ДР(1),
где
У = -
Ро
(12)
(13)
(14)
Ь(Ь + 2Ро)
Точкой обозначена производная по т. Как легко проверить, уравнение (13) эквивалентно нормальной системе двух гамильтоновых уравнений первого порядка:
I = дН / дР
Р1 =-дН / II,
с функцией Гамильтона Н = (Р/ -1)I + 2yB(I);
(15)
(16)
(17)
где В(I) - первообразная для функции Др (I), т.е.
В (I) = | Др (I)<1/ . (18)
Из классической механики хорошо известно, что гамильтонова система (15), (16) описывает также одномерную точечную частицу с координатой I , следовательно, уравнение (13) является также уравнением Эйлера-Лагранжа для этой одномерной механической частицы. Роль времени в классической динамике частицы играет безразмерная переменная т, определённая формулой (12). Таким образом, гамильтонова система (15), (16), полученная для интенсивности огибающей оптического импульса, имеет ещё и механический смысл, если под I понимать не полевую функцию, а координату точки. Так как функция Гамильтона (17) не зависит явно от времени, то механическая энергия этой точки сохраняется. Хорошо известно, что решения !(т), асимптотически стремящиеся к нулю на бесконечности, существуют (если они, конечно, существуют) при нулевой механической энергии точки. Это означает, что огибающая оптического импульса является в этом случае локализованной функцией. Именно такие локализованные решения уравнения (13) и представляют интерес. Применение механической аналогии позволяет далеко продвинуться в теории солитонов [8], [9].
Решение уравнений Гамильтона (15), (16), а следовательно, и уравнения (13) легко проводится с помощью теории канонических преобразований. Как следует из теоремы Якоби-Пуанкаре [10], если существует дважды дифференцируемая функция £(I, р, т), такая что
|д2£ / д!др| ф 0, то преобразование (I,Р1) ^ р)
генерируемое этой функцией:
Р, =д£ / д!; q = д£ / др , (19)
2
2
2
является каноническим, а новая функция Гамильтона имеет вид
Н(д, р, т) = Н (I(д, р, т), Р1 (q, р, т), т) +
д^ (. ( ) ) (20)
+—и (q, p, т), p, т).
дт
В теории канонических преобразований для двумерного фазового пространства можно исходить из четырёх типов производящих функций, зависящих от одной новой и одной старой переменной [10]. Для дальнейших целей наиболее подходящей является функция, зависящая от старого импульса и новой координаты: F = R2(Р1,q,т) .Тогда явный вид преобразований найдётся из разрешения уравнений
, дF дF
I =—; р = —, (21)
дР1 дq
а новая функция Гамильтона
Н(д, р, т) = Н (I(д, р, т), Р7 (q, р, т), т)-
дF/ (22)
(р (q, p, т), q, т).
дт
Будем считать в гамильтониане (17) первое слагаемое невозмущённым гамильтонианом, а второе слагаемое, обусловленное нелинейным откликом среды, - возмущением. Перейдём от динамической системы {I, р, Н} к системе {д, р, Н} в представлении
взаимодействия с помощью производящей функции F (р, д, т), такой, чтобы в новых переменных новая функция Гамильтона определялась бы только возмущением. Это делается с помощью производящей функции вида
F = дАгЛР} - дт .
(23)
которая, как легко проверить, является полным интегралом невозмущённого уравнения Гамильтона-Яко-
би: дF/дт = (р2 -/др.
Из уравнений (21) следует явный вид преобразований:
I = д^2(р + т); р. = Щр + т).
(24)
а из (22) - явный вид новой функции Гамильтона: Н = 2yB(qch2 (р + т)). (25)
Новые уравнения Гамильтона
д = дН / др; р = -дН / дд (26)
определяют динамику новых канонических переменных д и р:
д = 4уд^( р + р +т)ДР(qch2( р + т)),
р = -2ych2( р +т)ДР(д^2( р + т)).
Так как в представлении взаимодействия гамильтониан (25) зависит от времени в комбинации (р + т), то первый интеграл системы легко угадывается:
(27)
2уВ(дсИ1 (р + т)) - д = Е = 0 .
(28)
Действительно, взяв полную производную по времени от этого выражения, получим тождество с учётом (25) и (26). Механическая энергия Е принята за ноль с учётом сделанного выше замечания о локализованных решениях (не путать Е с энергией оптического импульса - это просто произвольная постоянная, которая в механической интерпретации имеет смысл механической энергии одномерной частицы).
Для упрощения решения динамической задачи в представлении взаимодействия разделим первое уравнение (27) на второе. В результате получим
д / д = -2 р^( р + т).
(29)
Заметим, что вид последнего уравнения не зависит от вида функции нелинейного отклика, а определяется только каноническим преобразованием (24), которое, в свою очередь, определяется невозмущённым гамильтонианом Но = (р2 - Щ . Этот невозмущённый гамильтониан определяется в конечном итоге линейной частью уравнения (2).
Знание одного первого интеграла системы понижает её порядок на единицу, следовательно, до одного уравнения, в качестве которого удобнее всего взять уравнение (29). Вторым уравнением для решения задачи следует взять первый интеграл (28), из которого следует, что д является некоторой функцией гиперболического косинуса, зависящего от аргумента (р+т). Подстановка этой функции в (29) даёт временную зависимость р(т), что и довершает решение поставленной задачи. Так, к примеру, для керровской нелинейности В^) = ц/2 реализация этой схемы даёт: р= -2т, д = 1/ 2уцсЛ4 (т), что является точным решением.
Заключение
Таким образом, задача нахождения локализованных решений уравнения (2) сведена с помощью теории гамильтоновых систем к двум уравнениям (28) и (29), первое из которых является алгебраическим. Примечательно, что решение представляет собой некоторую функцию гиперболического косинуса независимо от конкретного вида функции нелинейного отклика.
Благодарности
Автор выражает благодарность В.В. Котляру за полезное обсуждение решаемой задачи.
Работа выполнена при поддержке Министерства образования и науки РФ в рамках реализации мероприятий Программы повышения конкурентоспособности СГАУ среди ведущих мировых научно-образовательных центров на 2013-2020 годы.
Литература
1. Агравал, Г.П. Нелинейная волоконная оптика. - М.: Мир, 1996. - 324 с.
2. Алименков, И.В. Решение расширенного уравнения распространения импульсов в оптических волокнах / И.В. Алименков, Ю.Ж. Пчёлкина // Компьютерная оптика. - 2014. - Т. 38, № 1. - С. 28-30.
3. Chen, Ch.-L. Foundations for Guided-Wave Optics / Chin-Lin Chen. - Wiley, 2007. - 462 p.
4. El-Wakil, S.A. New periodic and soliton solutions of nonlinear evolution equations / S.A. El-Wakil // Applied Mathematics And Computation. - 2008. - Vol. 197. -P. 497-506.
5. Алименков, И.В. Интегрирование в элементарных функциях двунаправленного уравнения распространения импульсов в оптических волокнах для степенной нелинейности / И.В. Алименков, Ю.Ж. Пчёлкина // Компьютерная оптика. - 2014. - Т. 38, № 3. -С. 377-379.
6. Кившарь, Ю.С. Оптические солитоны. От волоконных световодов к фотонным кристаллам / Ю.С. Кившарь, Г.П. Агравал. - М.: Физматлит, 2005. - 648 с.
7. Алименков, И.В. Решение расширенного уравнения распространения импульсов в оптических волокнах для конкурирующей нелинейности / И.В. Алименков, Ю.Ж. Пчёл-кина // Компьютерная оптика. - 2014. - Т. 38, № 4. -С. 686-689.
8. Раджараман, Р. Солитоны и инстантоны в квантовой теории поля / Р. Раджараман. - М.: Мир, 1985. - 416 с.
9. Тахтаджян, Л.А. Гамильтонов подход в теории солитонов / Л.А. Тахтаджян, Л.Д. Фаддеев. - М.: Наука, 1986. - 528 с.
10. Шмутцер, Э. Основные принципы классической механики и классической теории поля. - М: Мир, 1976. - 160 с.
References
1. Agrawal, G.P. Nonlinear Fiber Optics / G.P. Agrawal. -Moscow: "Mir" Publisher, 1996. - 324p. - (In Russian).
2. Alimenkov, I.V. Solution of expanded pulse-propagation equation for optical fiber / I.V. Alimenkov, Yu.Z. Pchelkina // Computer Optics. - 2014. - Vol. 38(1). - P. 28-30.
3. Chen, Ch.-L. Foundations for Guided-Wave Optics / Chin-Lin Chen. - Wiley, 2007. - 462 p.
4. El-Wakil, S.A. New periodic and soliton solutions of nonlinear evolution equations / S.A. El-Wakil // Applied Mathematics and Computation. - 2008. - Vol. 197. -P. 497-506.
5. Alimenkov, I.V. Integration in elementary functions of two-way pulse-propagation equation in optical fiber for power nonlinearity / I.V. Alimenkov, Yu.Z. Pchelkina // Computer Optics. - 2014. - Vol. 38(3). - P. 204-206.
6. Kivshar, Y.S. Optical solitons. From Fibers to Photonic Grystals / Y.S. Kivshar, G.P. Agrawal. - Moscow: "Fizmat-lit" Publisher, 2005. - 648 p. - (In Russian).
7. Alimenkov, I.V. The solution of expanded pulse-propagation equation in optical fiber for competing nonlinearity / I.V. Alimenkov, Yu.Z. Pchelkina // Computer Optics. -2014. - Vol. 38(4). - P. 686-689.
8. Rajaraman, R. Solitons and instantons in quantum field theory / R. Rajaraman. - Moscow: "Mir" Publisher, 1985. - 416 p. -(In Russian).
9. Takhtajan, L.A. Hamilton approach in theory of solitons / L.A. Takhtajan, L.D. Faddeev. - Moscow: "Nauka" Publisher, 1986. - 528 p.
10. Schmutzer, E. Grundprinzipien der klassischen Mechanik und der klassischen Fieldtheorie (kanonischer Apparat) / E. Schmutzer. - Berlin:VEB Deutscher Verlag der Wissenschaften, 1973. - 160 p.
HAMILTONIAN APPROACH IN NONLINEAR FIBER OPTICS
I. V. Alimenkov Samara State Aerospace University
Abstract
The Hamilton approach is applied to the solution of an expanded pulse-propagation equation in optical fibers for an arbitrary nonlinearity. It is shown that the solitonic solution is a function of the hyperbolic cosine at any nonlinearity.
Keywords: optical fiber, expanded pulse-propagation equation, Hamilton systems, arbitrary nonlinearity, solitonic solution.
Сведения об авторе
Алименков Иван Васильевич, 1949 года рождения. В 1977 году с отличием окончил Куйбышевский государственный университет по специальности «Физика». Кандидат физико-математических наук, работает в должности доцента кафедры прикладной математики СГАУ. Область научных интересов - нелинейная физика.
E-mail: [email protected] .
Ivan Vasilyievich Alimenkov, 1949 year of birth. In 1977 graduated with honours from Kuibyshev State University on a speciality "Physics". Candidate in Physics and Mathematics, works as associated professor of Applied Mathematics sub-department of SSAU. Research interests - nonlinear physics.
Поступила в редакцию 21 декабря 2014 г.
