Решение расширенного уравнения распространения импульсов в оптических волокнах для конкурирующей нелинейности

Автор: Алименков Иван Васильевич, Пчлкина Юлия Жиганшевна

Журнал: Компьютерная оптика @computer-optics

Рубрика: Дифракционная оптика, оптические технологии

Статья в выпуске: 4 т.38, 2014 года.

Бесплатный доступ

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

Волоконный световод, расширенное уравнение распространения, точное решение, конкурирующая нелинейность, солитонное решение

Короткий адрес: https://sciup.org/14059293

IDR: 14059293

Текст научной статьи Решение расширенного уравнения распространения импульсов в оптических волокнах для конкурирующей нелинейности

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

E ( r , t ) = e x F ( x , У ) A ( z , t ) exp { i ( в o2 - to ot ) } ,        (1)

\                4,                    J- ( x 2 + y 2) 1

где f ( x, y ) - гауссовская функция вида exp <----- 2--- ^

[ w с радиусом моды w, A(z, t) – комплексная огибающая импульса, too - несущая частота, в o = to on (too)/ c -центральное волновое число.

Для огибающей оптического импульса выведено

  • [2]    уравнение

. ( d A „ d A ^ 1 d 2 A в 2 d 2 A t d z 1 d t J 2 в o d z 2 2 d t1 +Ap (| A 2) A = 0,

названное расширенным уравнением распространения, которое существенно отличается от традиционного, называемого основным уравнением распространения [1], наличием второй производной по координате. Здесь в 1 = 1/ v g - величина, обратная групповой скорости, β 2 дисперсия групповой скорости, Ap(| A |2) - нелинейная поправка к постоянной распространения моды в линейном приближении, которая выражается через нелинейную часть Δ n показателя преломления:

kot fA n F ( x , У ) d x d y

Ав =     -------й-------,

J] I F ( x , y )| d x d y

где k o = to o I c .

В области прозрачности волновода Δβ является вещественной функцией. Так как β2 в области прозрачности отрицательно, то тип уравнения (2) – эллиптический. Если в уравнении (2) отбросить одну из вторых производных, то полученное уравнение параболического типа сводится к нелинейному уравнению Шрёдингера, которое является на сегодняшний день самым изученным из нелинейных уравнений, допускающих солитонные решения. Сводку различных его решений можно посмотреть, например, здесь: ,

Полные уравнения второго порядка рассматривались и ранее. Так, в [3] выведено уравнение гиперболического типа для неограниченной керровской среды в упрощённой модели поля, не зависящего от поперечных координат, но затем это уравнение упрощается до нелинейного уравнения Шрёдингера. Следует упомянуть и статью [4], в которой методом секанса и tanh-методом решаются некоторые полные уравнения второго и даже третьего порядка, однако полученные решения сингулярны. Решения же в статьях [2], [5] являются гладкими локализованными функциями, как и в данной работе.

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

A n = n 2 | E |2 + n 4 | E |4 +••• , (4) что согласно (3) приводит к степенному разложению функции Δβ:

Ав = у| a| 2 +ц| a| 4 +••• , (5) где параметры γ и µ зависят от характеристик световода. При достаточно малых интенсивностях оптического поля вторым слагаемым в (5) пренебрегают.

Полученная нелинейность называется керровской, и солитонное решение уравнения (2) для этого случая найдено в [5].

С ростом интенсивности вводимого излучения наблюдается [6] отклонение от керровской зависимости показателя преломления и необходимо учитывать второе слагаемое в формуле (5). Такая нелинейность называется конкурирующей (иногда этот термин применяют только в случае, если слагаемые в (5) различаются знаком).

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

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

Основной формализм

В [5], [7] показано, что если искать локализованное решение уравнения (2) в виде

A ( z , t ) = R ( z , t )exp { iqz } , (6) где R действительная функция, а q произвольный параметр, являющийся поправкой к центральному волновому числу в o , то функция R определяется

пиковое значение напряжённости выражается следующим образом:

E max

двумя квадратурами:

I /    dR2      2 =^ ,(7)

J R1 - - B ( R 2)/ pR2

R 2

B(R2) = J Ap(I)dI.(8)

Здесь

^ = I 2 p Po (z — z — vt),(9)

V1 -Po P2 vo’ p = q(1 + q /2po),(10)

v = vg (1 + q / в o),(11)

2 7 p / Y

1 + + 71 + 16 p Ц / 3 y 2

Отсюда находим

Y E^ Ц Е 4ax max max p 2      3

Теперь из формулы (10) выражаем

q = p o

- 1 ±

V

или, подставляя сюда (18), окончательно получим

q = p o

(

- 1 ±

V

h e l ) з в o J

A

где z 0 – произвольная постоянная, определяющая начальное положение импульса.

Полагая в (8)

Ap = y I + Ц 1 2

и подставляя полученное выражение в (7), имеем

J            d R            = ^ .

J R 71 - ( Y R 2/2 p \u. R 4/3 P )

Интеграл в левой части (13) с помощью подстановки u =1/ R сводится к известному [8] интегралу, и в результате получим

1 . 1   1/R2-y/4 p.

— Arch .= ^

  • 2      7y2 /16 p 2 /3 p

Обращая это выражение, находим

2 J p / Y

R = ! I.

71 + 71 + 16 p ц /3 Y 2 ch 2 ^

Формула (15) описывает вещественную огибающую волнового пакета, локализованного вдоль направления z = z o + vt и движущегося с постоянной скоростью (11).

C учётом (6), солитонное решение уравнения (2) имеет вид:

A ( z , t )

27 p / Y exp{ iqz }

1 + 71 + 16 p Ц /3 Y 2 ch 2 ^

В том, что (16) является решением уравнения (2), можно убедиться прямой подстановкой. Последняя формула представляет собой однопараметрическое семейство решений, в котором свободным параметром является q .

От свободного параметра q (а вместе с ним и от p ) удобно перейти к новому свободному параметру [9] амплитудного характера. Действительно, из соображений размерности и из формулы (15) следует, что

Формулы (16), (9), (11), (18), (20) решают поставленную задачу.

Заключение

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

Работа выполнена при поддержке Министерства образования и науки РФ в рамках реализации мероприятий Программы повышения конкурентоспособности СГАУ среди ведущих мировых научно-образовательных центров на 2013–2020 годы.

Статья научная