УДК 535.2: 621.391
МОДИФИКАЦИЯ ДИСКРЕТНЫХ СПЕКТРАЛЬНЫХ ПАРАМЕТРОВ ОПТИЧЕСКИХ СОЛИТОНОВ В ВОЛОКНАХ С ПЕРЕМЕННОЙ ДИСПЕРСИЕЙ
К. С. Гочелашвили1, А. А. Сысолятин1, A. И. Конюхов1'2, Л.Л. Мельников1,3, П. А. Маврин2, М.Ю. Салганский4
Спектральные параметры солитонов являются результатом нелинейного преобразования Фурье для импульса, распространяющегося в волокне. В работе предложен новый подход для управления дискретной составляющей спектральных параметров. При наличии резонанса между осцилляциями поля и дисперсией волокна происходит, изменение действительной части спектральных параметров, что приводит к разделению многосолитоннного импульса на фундаментальные солитоны. Предложено использовать данный эффект для декодирования сигнала в линиях связи, работающих на основе нелинейного преобразования Фурье.
Ключевые слова: оптический солитон, нелинейное преобразование Фурье, обратная задача рассеяния.
Введение. Солитоны нелинейного уравнения Шредингера [1] характеризуются дискретными спектральными параметрами Лj обратной задачи рассеяния. Параметры Лj рассчитываются при помощи нелинейного преобразования Фурье [2, 3]. Кодирование сигналов на основе нелинейного преобразования Фурье стало возможным c развитием методов когерентного детектирования [4]. В настоящее время синтез и детектирование солитонов с заданными параметрами осуществляется методами цифровой обработки сигналов, что накладывает ограничения на скорость передачи. Для перехода к тераби-товым скоростям требуется использование оптических методов обработки сигнала.
1 ИОФРАН, 119991 Россия, Москва, ул. Вавилова, 38; e-mail: kai@optics.sgu.ru.
2 Саратовский государственный университет, 410012 Россия, Саратов, ул. Астраханская, 83.
3 Саратовский государственный технический университет, 410054 Россия, Саратов, ул. Политехническая, 77.
4 Институт химии высокочистых веществ РАН, 603950 Россия, Нижний Новгород, ул. Тропинина, 49.
Оптические солитоны, подчиняющиеся нелинейному уравнению Шредингера, взаимодействуют упруго, то есть при взаимодействии сохраняют свою групповую скорость и амплитуду. Эти особенности делают солитоны привлекательными для использования в высокоскоростных линиях связи, но в то же время затрудняют управление их параметрами. Под управлением мы подразумеваем возможность изменения числа солитонов, их скоростей, длительности и амплитуды. Отступление от классического уравнения Шредингера делает взаимодействие солитонов неупругим. Параметры солитонов могут изменяться при столкновениях, при взаимодействии с дисперсионными волнами. Достаточно обширный материал по управлению солитонами представлен в [5-7]. Для изменения солитонных состояний было предложено использовать нелинейность пятого порядка, нелинейность с насыщением, специального вида дисперсию, двухмодовые волокна. Солитоны могут образовывать структуры и взаимодействовать в лазерах, описываемых в рамках модели комплексного уравнения Гинзбурга-Ландау [8, 9]. Столкновение солитона с дисперсионной волной может приводить к его ускорению/замедлению, сдвигу его несущей частоты [10].
Использование специального типа нелинейности или дисперсии для управления параметрами солитонов трудно реализуемо на практике. Мы предлагаем способ, который может быть осуществлен при помощи одномодового оптического волокна. Модуляция диаметра волокна приводит к изменению дисперсии. В экспериментах [11, 12] было показано, что при совпадении периода осцилляций дисперсии и периода осцилляций двухсолитонного импульса, импульс разделяется на два фундаментальных солитона. Распространение солитонов в волокне с периодическим изменением дисперсии сопровождается изменением дискретных спектральных параметров А^- солитонов [13, 14]. Данная особенность может быль использована для кодирования/декодирования сигнала в системах связи на основе нелинейного преобразования Фурье.
Дискретный солитонный спектр. Распространение солитона в оптоволокне подчиняется нелинейному уравнению Шредингера (НУШ) [1].
дА в2 д2 А ., ..2.. .
Ж + 4 дТА = ^|А|2А(^Т), (1)
где А(г,т) - нормированная комплексная огибающая амплитуды электрического поля, г - дистанция распространения, т - время в бегущей системе координат — г/и), и - групповая скорость импульса, г - мнимая единица, в - коэффициент дисперсии второго порядка, 7 - коэффициент керровской нелинейности.
Для нахождения солитонных решений необходимо воспользоваться методом обратной задачи рассеяния солитона, т.е. перехода части его энергии в другие спектральные компоненты [1-3]. Для сигнала q(r) = A(z,T), который определен на интервале Ti < т < T2 при фиксированном значении z, решается система двух дифференциальных уравнений:
dv =( -гЛ q(Tм v(t, Л), (2)
дт ^ -q*(T) гЛ J К J W
где Л - солитонный параметр, * обозначает комплексное сопряжение. Начальное значение для системы (2) = (e-iXTl, 0)T. Два элемента матрицы рассеяния а(Л) и Ь(Л) определяются из решения (2) v(T2^) = (v1,v2)T как а(Л) = v1 ехр(гЛТ2) и Ь(Л) = v2 ехр(-гЛТ2). Солитонам соответствуют решения уравнения а(Л) = 0. Дискретные значения Лj являются комплексными. Действительная часть Re^j) определяет групповую скорость солитона, мнимая часть Im^j) - его амплитуду (см. главу 5 в [1]). Решение уравнения (1) находилось численно [15]. Для расчета матрицы (2) использовался "layer-peeling" метод [3]. Для нахождения комплексных корней уравнения а(Л) = 0 использовался метод Ньютона.
Разделение многосолитонного бризера. Рассмотрим распространение многосолитон-ного импульса в волокне с переменной дисперсией [11, 12]. В этом случае НУШ (1) имеет переменные коэффициенты
&(z) = <&> (1 + 0.2 sin(2n z/zm + Vm)), (3)
Y(z) = <Y>(1 - 0.028 sin(2n z/zm + <pm)), (4)
<в2> - среднее значение коэффициента дисперсии второго порядка, <y> - среднее значение коэффициента нелинейности, zm - период модуляции, - фаза модуляции. Выберем параметры волокна, соответствующие эксперименту [12] <в2> = -12.76 пс2/км, <Y> = 8.2 Вт ■ км-1.
Рассмотрим разделение многосолитонного импульса, заданного в виде сигнала
A(0,T) = sech(l) , (5)
где A(0, т) - амплитуда сигнала при z = 0, то - начальная длительность импульса. Для целого значения S импульс (5) имеет дискретные солитонные параметры Лj• = i(j —1/2), где j = 1, 2, 3 ... , j < S. Поскольку данный набор представлен чисто мнимыми числами, то все солитоны распространяются с одинаковыми групповыми скоростями и образуют бризер. При изменении действительной части Лj• несущая частота изменяется и за
счет дисперсии изменяется групповая скорость солитона. Такой солитон отделяется от основного пакета. Импульс (5) при распространении в волноводе с постоянной дисперсией e2(z) = const и нелинейностью 7(z) = const имеет периодичность по z, связанную с нелинейной интерференцией составляющих его солитонов. Период изменения формы импульса равен zo = пт0 |2в2|-1. Для импульса с длительностью т0 = 1.13 пс величина zo = 0.158 км.
Расстояние, км Re(/„) Re(A,)
Рис. 1: Разделение пятисолитоного импульса, Б = 5. (а) интенсивность I(г,т) = |А(г,т)|2; (б) линии уровня Яе(а(А)) = 0 - сплошная линия, 1ш(а(А)) = 0 - пунктир, z = 0.32 км; (в) то же что и на рис. (б) при z = 1.58 км. Цифрами 7 = 1, 2, 3, 4, 5 отмечены значения Аj• при z = 0.32 (б) и z =1.58 (в) км.
На рис. 1 показано разделение пятисолитонного импульса Б = 5 при zm = 0.16 км, <^>т = 0. Дистанция распространения равна 10zo. Период модуляции дисперсии близок к периоду солитона z0. Поэтому разделение импульса происходит на короткой дистанции распространения, после трех периодов осцилляции дисперсии (рис. 1(а)). Нули а(А^) можно найти на пересечении линий И,е(а(А)) = 0 и 1ш(а(А)) = 0. По характеру изменения этих линий можно отследить движение солитонных параметров Аj в комплексной плоскости И,е(А), 1ш(А). Начальный импульс (5) при Б = 5 имеет пять величин, находящихся на линии И,е(А) = 0: Аj = ¿(7 — 0.5), 7 = 1, 2, 3,4,5. При моделировании было обнаружено, что разделение многосолитонного импульса происходит через изменение действительной части пары "соседних" параметров. При z = 0.32 км появляется действительная часть А3 и А4 (рис. 1(б)), начальные значения которых А3 = ¿2.5, А4 = ¿3.5. Разделение солитонов происходит симметрично И,е(А3) = — И,е(А4), что диктуется законом сохранения импульса [1, 5]. Величина А5 = ¿4.5 уменьшается до А5 = ¿3.66 (рис. 1(б)), что связано с излучением дисперсионной волны. Изменение пары других
параметров А1 и Л2 происходит при г > 0.79 км по аналогичному сценарию (рис. 1(в)). После прохождения отрезка волокна длиной 1.5 км пятисолитонный импульс распадается на пять фундаментальных солитонов с различными групповыми скоростями, определяемыми Ке(А^). Мы предполагаем, что по соотношению пиковых интенсивно-стей выходных импульсов и по временному интервалу между ними можно будет детектировать порядок исходного солитона.
Влияние вынужденного комбинационного рассеяния. В экспериментах [11, 12] было показано, что вынужденное комбинационное рассеяние (ВКР) приводит к асимметричному разделению двухсолитонного импульса. Пиковые значения выходных импульсов при наличии ВКР становятся отличными друг от друга, нарушается закон сохранения импульса. При наличии дисперсии третьего порядка и вынужденного комбинационного рассеяния динамика импульса подчиняется уравнению
дА .&(*) д2 А вз (г) д3 А . ,.,|2, . , ,
(6)
где вз(г) = (вз> (1 + 0.01 81п(2п г/гт + рт)), (вз> = 0.076пс3км-1, 1к(г) = <1к> (1 — 0.028 вт(2п г/гт + <^т)), <7К) = 1.8Вт-1км-1. Нелинейный отклик вынужденного комбинационного рассеяния определен через затухающие осцилляции [16] (д2Я/дг2) + 2Т2-1 (дЯ/дг) + П2д(М) = П2|А(г, т)|2, Т2 = 32фе, П(2п)-1 = 13.1 ТГц.
Рис. 2: Разделение пятисолитонного импульса при наличии ВКР. (й) динамика со-литона при отсутствии модуляции дисперсии гт = то; (б) динамика солитона при наличии модуляции дисперсии гт = 0.016 км, = п; (в) дискретный солитонный спектр при гт = то (окружности) гт = 0.016 км (ромбы) после прохождения волокна длиной г = 0.79 км.
Известно, что ВКР приводит к разделению многосолитонного импульса [6]. На рис. 2(а) показано распространение пятисолитонного импульса Б = 5 при постоянной дисперсии и нелинейности = то). Солитоны разделяются по групповым скоростям. Для интенсивных импульсов ВКР оказывает сильное влияние на достаточно коротких дистанциях распространения импульса. Поэтому для эффективного воздействия модуляции дисперсии на динамику солитонов следует использовать малый период модуляции. Динамика солитонов при zm = zo/10 = 0.016 км показана на рис. 2(б). Начальные значения z = 0) солитонных параметров находятся на оси И,е(А) = 0 : А1 = ¿0.5, А2 = ¿1.5, А3 = ¿2.5, А4 = ¿3.5, А5 = ¿4.5. Параметры солитонов Аj после прохождения волокна с постоянным диаметром (положение Аj обозначено окружностями) и волокна с периодическим изменением дисперсии (положение Аj обозначено ромбом) показаны на рис. 2(в).
При постоянной дисперсии распад многосолитонного импульса происходит через изменение групповой скорости (И,е^-)) в зависимости от амплитуды солитона, определяемой ). В волокне с периодическим изменением дисперсии моделирование показывает возможность выравнивания амплитуд солитонов (рис. 2(в)). Один из солитонов исчезает за счет излучения дисперсионной волны.
Заключение. В данной работе мы рассмотрели возможность управления параметрами солитонов в рамках модели классического нелинейного уравнения Шредингера с переменными коэффициентами. Показано, что многосолитонный импульс можно разделить на отдельные составляющие, определяемые фундаментальными солитонами. Такое разделение может использоваться для детектирования солитонных состояний в линиях связи, использующих нелинейное преобразование Фурье для кодирования сигнала. После прохождения волокна с периодическим изменением дисперсии многосолитонный импульс, несущий полезный сигнал, даёт индивидуальный набор импульсов, используя который можно детектировать состояние исходного многосолитонного сигнала. Важным преимуществом предложенного подхода является использование полностью волоконной системы для детектирования солитонных состояний. Помимо практических применений резонансные явления в оптических волокнах с переменной дисперсией предлагают новую платформу для изучения физики солитонных взаимодействий.
Работа поддержана грантом РНФ № 17-12-01564.
ЛИТЕРАТУРА
[1] С. А. Ахманов, В. А. Выслоух, А. С. Чиркин, Оптика фемтосекундных лазерных импульсов (М., Наука, 1988).
[2] M. I. Yousefi, F. R. Kschischang, IEEE Transactions on Information Theory 60(7), 4312 (2014).
[3] M. I. Yousefi, F. R. Kschischang, IEEE Transactions on Information Theory 60(7), 4329 (2014).
[4] S. K. Turitsyn, J. E. Prilepsky, S. T. Le, et al., Optica 4(3), 307 (2017).
[5] Н. Н. Ахмедиев, А. Анкевич, Солитоны. Нелинейные импульсы и пучки (М., Физматлит, 2003). N. N. Akhmediev, A. Ankiewicz, Solitons - Nonlinear Pulses and Beams, 1st edn (Chapman & Hall, London, 1997).
[6] Ю. С. Кившарь, Г. П. Агравал, Оптические солитоны. От волоконных световодов до фотонных кристаллов (М., Физматлит, 2005). Yu. S. Kivshar, G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic Press, San. Diego, 2003).
[7] B. A. Malomed, Soliton Management in Periodic Systems (Springer, 2006).
[8] P. Grelu, N. Akhmediev, Nature Photonics 6(2), 84 (2012).
[9] Jose M. Soto-Crespo, N. Devine, and N. Akhmediev, J. Opt. Soc. Am. B 34(7), 1542 (2017).
[10] I. Oreshnikov, R. Driben, and A. V. Yulin, Opt. Lett. 40(23), 5554 (2015).
[11] A. A. Sysoliatin, E. M. Dianov, A. I. Konyukhov, et al., Laser Physics 17(11), 1306 (2007).
[12] A. A. Sysoliatin, A. K. Senatorov, A. I. Konyukhov, et al., Opt. Express 15(25), 16302 (2007).
[13] A. A. Sysoliatin, A. I. Konyukhov and L. A. Melnikov, in: Numerical Simulations of Physical and Engineering Processes, Ed. by Prof. Jan Awrejcewicz, (InTech, London, 2011), Chap. 13, р. 277.
[14] A. I. Konyukhov, A. S. Plastun, M. A. Dorokhova, L. A. Melnikov, Laser Phys. Lett. 12(5), 055103 (2015).
[15] O. V. Sinkin, R. Holzlohner, J. Zweck, C. R. Menyuk, J. Lightwave Technol. 21(1), 1 (2003).
[16] Э. М. Беленов, А. В. Назаркин, И. П. Прокопович, Письма в ЖЭТФ 55(4), 223 (1992). E. M. Belenov, A. V. Nazarkin, I. P. Prokopovich, JETP Letters 55(4), 218 (1992).
Поступила в редакцию 16 октября 2017 г.