СВОЙСТВА ВРЕМЕННЫХ ОПТИЧЕСКИХ СОЛИТОНОВ В ОПТИЧЕСКИХ ВОЛОКНАХ И ВОЗМОЖНОСТЬ ИХ ИСПОЛЬЗОВАНИЯ В ТЕЛЕКОММУНИКАЦИЯХ. ЧАСТЬ 1
Е.И. Андреева 1 М.С. Былина С.Ф. Глаголев 1 П.А. Чаймарданов1
^анкт-Петербургский государственный университет телекоммуникаций им. проф. М.А. Бонч-Бруевича, Санкт-Петербург, 193232, Российская Федерация *Адрес для переписки: [email protected]
Информация о статье
УДК 621.39, 530.182 Язык статьи - русский
Ссылка для цитирования: Андреева Е.И., Былина М.С., Глаголев С.Ф., Чаймарданов П.А. Свойства временных оптических солитонов в оптических волокнах и возможность их использования в телекоммуникациях. Часть 1 // Труды учебных заведений связи. 2018. Т. 4. № 1. С. 5-12.
Аннотация: В первой части работы приведено теоретическое обоснование существования и устойчивости солитонов в одномодовом оптическом волокне без потерь. Показана возможность применения солитонов в оптических системах связи. Рассмотрены проблемы, связанные с взаимодействием солитонов. Приведены результаты моделирования процессов распространения солитонов по оптическому волокну. Во второй и третьей частях будут рассмотрены вопросы сохранения солитонного режима в одномодовом оптическом волокне с помощью управления потерями и дисперсией в линейном тракте волоконно-оптической системы связи.
Ключевые слова: оптические солитоны, временные оптические солитоны, одномодовое оптическое волокно, волоконно-оптическая система связи, дисперсия групповых скоростей, фазовая самомодуляция.
Введение
Условия существования временных светлых оптических солитонов в одномодовых оптических волокнах и возможности их использования для телекоммуникаций были высказаны Хасегавой и Таппертом в 1973 году. В 1980 году Молленауэром, Столеном и Гордоном существование светлых со-литонов в одномодовых оптических волокнах было экспериментально подтверждено [1, 2]. В работе мы не будем обсуждать пространственные оптические солитоны и временные темные солитоны. Все это систематически представлено в [2]. В дальнейшем изложении временные светлые оптические солитоны мы будем называть просто соли-тонами.
Работа включает три части. Первая посвящена вопросам теоретического обоснования существования и устойчивости солитонов в оптических волокнах без потерь. В следующих частях реальные потери в оптических волокнах учитываются: во второй части рассматриваются практические вопросы сохранения солитонного режима в оптических волокнах с использованием техники управления потерями в линейном тракте волоконно-
оптической системы связи, а в третьей - вопросы сохранения солитонного режима с помощью управления дисперсией в линейном тракте волоконно-оптических систем связи (ВОСС). В работе для проверки теории и получения наглядных результатов широко используется имитационное моделирование процессов, происходящих в линейном тракте ВОСС.
1. Теоретическое обоснование существования солитонов в оптическом волокне без потерь
Солитоны формируются внутри одномодового оптического волокна вследствие баланса между дисперсией групповых скоростей (ДГС) 02 и фазовой самомодуляцией (ФСМ).
Теоретически существование солитонов в оптических волокнах может быть обосновано с помощью решения нелинейного уравнения Шредингера (НУШ), которое для оптического волокна без потерь имеет вид [1-4]:
. дА 1 „ д2А 1.1?.
^ = ;?Р2д75-УИ12А (1)
где А - нормированная комплексная амплитуда напряженности (огибающая импульса) электрического поля световой волны, распространяющейся в направлении оси 2 по оптическим волокнам, Вт05; Т - время, отсчитываемое от середины импульса, пс; 02 - ДГС, пс2/км; у - коэффициент нелинейности, (Вт-км)-1; 1 - расстояние от начала оптического волокна, км.
Введем обозначения:
и = N -и, г = г/Ьа, т = Г/Го, и = А/^Т0, И2 = , Ьв = Го2/|р2|, = 1/(уРо),
(2)
где Ьо - дисперсионная длина оптического волокна; То - полуширина импульса на входе в оптическое волокно, пс; Ьмь - нелинейная длина волокна; Ро - пиковая мощность импульса.
Учитывая, что в области аномальной дисперсии 02 < 0, и используя обозначения (2), преобразуем (1) к стандартному виду НУШ, в котором все параметры безразмерные:
ди 1 д2и
+ _ ._ +Ы2.и = 0.
(3)
НУШ позволяет решить прямую задачу об изменении формы и спектра оптического сигнала произвольной формы после прохождения его по оптическому волокну с известными параметрами.
Решением обратной задачи рассеяния является определение параметров сигнала на входе в оптическое волокно для достижения определенной цели. Например, одной из целей решения обратной задачи является определение формы и мощности (энергии) оптического импульса на входе оптического волокна, который сохраняет свою форму (не расширяется) или периодически восстанавливает свою форму при прохождении по оптическому волокну больших расстояний.
Особое значение имеет решение, получившее название фундаментального солитона, который при распространении по оптическому волокну без потерь не меняет свою форму, т.е. не подвержен влиянию хроматической дисперсии (ХД). Фундаментальный солитон имеет огибающую напряженности электрического поля в форме гиперболического секанса, и его распространение по оптическому волокну без потерь описывается функцией [1-4]:
и(г,т) = БесИ(т) - ехр([ - г/2). (4)
Фундаментальный солитон также называют со-литоном 1-го порядка, так как для него N = 1. Из (2) можно определить соотношение между пиковой мощностью и длительностью солитона 1-го порядка [1]:
Ро = |р21/(у - Т02) * 3.11 - |р2|/(у - гд¥НМ), (5)
где Тгк/нм - длительность импульса, имеющего форму гиперболического секанса, на полувысоте по интенсивности:
= 2Т0 - 1п(1 + VI) * 1.763 - Т0.
(6)
На рисунке 1 представлены результаты расчетов по (5). Из графиков видно, что для стандартного одномодового оптического волокна (0.652) на длине волны 1550 нм при 02 = -20 пс2/км пиковая мощность составляет Р0 = 5 Вт при Т0 = 1 пс и 50 мВт при Т0 = 10 пс. Для оптического волокна со смещенной дисперсией (0.653) с 02 = -2 пс2/км при Т0 = 10 пс значение Р0 уменьшается до 5 мВт.
Ро, мВт
105
104 5 Вт
103 500 мВт
102 50 мВт
10 5 мВт
1
0 10 20 30 40 50 |р2|, пс2/км
Рис. 1. Зависимости пиковой мощности, необходимой
для возбуждения фундаментального солитона, от ДГС
На рисунке 2 показаны результаты расчета огибающей по интенсивности и2(0, т) и по амплитуде напряженности поля и(0, т) фундаментального солитона в начале оптического волокна. Для сравнения на рисунке представлены результаты аналогичных расчетов для гауссовского импульса, огибающая которого в начале оптического волокна имеет вид:
ид(0, т) = ехр(-т2/2), ид2(0, т) = ехр(-т2). (7)
Отметим, что солитонный импульс имеет большую ширину по сравнению с гауссовским. Для гауссовского импульса длительность на полувысоте по интенсивности равна TFwнм = 1.665Т0.
Среди солитонов высших порядков особую роль играют солитоны, имеющие начальную форму (при 2 = 0) в виде гиперболического секанса [1]:
и(0,т) = N - БесИ(т), (8)
где N - целочисленный порядок солитона. Из (2) следует выражение для мощности, необходимой для формирования солитона ^го порядка:
|Р2|/(У - Го2) *М2 - 3.11 - |Р2|/(у - Т^ни) (9)
Ром = N2 ■
Из (5) и (9) следует, что значение пиковой мощности, необходимой для создания солитона ^го порядка, в N раз больше, чем для возбуждения фундаментального солитона.
1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
u2(0, т), u2g(0, т)
Z0, км
-4 -3 -2 -1
а)
u(0, т), ug(0, т)
1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
Рис. 2. Формы огибающих фундаментального солитона и гауссовского импульса по интенсивности (а) и напряженности электрического поля (б)
Форма и спектр солитонов высших порядков периодически изменяется при распространении по оптическому волокну, возвращаясь в конце периода к исходной форме. Период изменения составляет z = zo = п/2 или:
п п Z° = 2 Ld =2
Л.
iß2!
= 0.322-
n-7V2
2 • iß2
(10)
На рисунке 3 представлены результаты расчетов по (10). Из рисунка 3 видно, что для стандартного оптического волокна (G.652) с ß2 = -20 пс2/км на длине волны Я = 1550 нм при То = 10 пс период солитона составляет Z0 = 8 км. Для оптического волокна со смещенной дисперсией (G.653) и ß2 = -2 пс2/км на длине волны Я = 1550 нм при Т0 = 10 пс период солитона составляет Z0 = 80 км.
На рисунке 4 показаны изменения формы и спектра солитонов 2-го и 3-го порядков по интенсивности на расстоянии, равном одному периоду.
Для понимания физического смысла распространения солитонов высших порядков рассмотрим изменения формы и спектра солитона третьего порядка (см. рисунок 4), которые обусловлены совместным действием ФСМ и ДГС.
100
80 км
10
8 км
0.1
0.01
0
10
40
50
20 30
|ß2|, пс2/км
Рис. 3. Зависимости периода изменения формы солитонов высших порядков от ДГС
0 -10"
Т 10 о 0.2 0.4 0.6 0.8 1.0 z/z0
а)
0 -10
Т 10оо!Го.4 0.6 08 z/z0
1.2 0.8 0.4
oj 11.1
" z/z0 в)
kS
1.2 0.8 0.4
^ ^о.^-о.Г"0^" 10
° 2120 г)
Рис. 4. Изменение формы (а, б) и спектра (в, г) солитона второго (а, в) и третьего (б, г) порядков на одном периоде
2
4
Известно, что ФСМ и ДГС приводят к возникновению частотной модуляции или чирпинга внутри оптического импульса. Чирпинг, вызванный ФСМ, всегда положителен, то есть частота внутри импульса нарастает. При этом частота на переднем фронте импульса смещается в стоксову область (более низкие частоты), а на заднем фронте в ан-тистоксову (более высокие частоты). На рисунке 4г при 2 / 20 = 0.2 хорошо видны типичные для ФСМ расширение спектра и появление модуляции (два пика).
Чирпинг, вызванный ДГС, зависит от ее знака. Он положителен при нормальной дисперсии (02 > 0) и отрицателен при аномальной (02 < 0). Поэтому при 02 < 0 ДГС может сжимать импульс, имеющий положительный чирпинг. Сокращает свою длительность только центральная часть импульса (хорошо видно на рисунке 4б при 2 / 20 = 0.2), т.к. только там частотная модуляция является линейной. Из-за существенного увеличения интенсивности в центральной части импульса его спектр также значительно изменяется (см. рисунок 4г для 2 / 20 = 0.3).
Важное свойство оптических солитонов - это их высокая устойчивость к возмущениям. Поэтому, хотя для формирования солитона и требуется специальная форма импульса и определенная пиковая мощность, он может быть создан, даже если форма импульса и пиковая мощность отклоняются от идеальных условий.
Влияние начальной формы импульса на формирование солитона может быть исследовано путем численного решения уравнения (3).
На рисунке 5а показана динамика импульса с пиковой мощностью, соответствующей мощности фундаментального солитона N = 1), но имеющего гауссово начальное распределение поля (7). Импульс распространяется в среде без потерь с аномальной дисперсией (02 < 0). Из рисунка 5а видно, что форма импульса изменяется при его распространении, поскольку вначале она отличается от гиперболического секанса, характерного для фундаментального солитона. Гауссовский импульс здесь асимптотически стремится к фундаментальному солитону. Эволюция фактически заканчивается при 1 / Ьо = 5, что соответствует примерно трем периодам солитона.
Похожая картина имеет место и для импульсов с другими начальными формами, например, с супер-гауссовской. Длительность солитона в конечном состоянии и расстояние, необходимое для эволюции импульса в солитон, зависит от начальной формы, но качественно поведение остается одним и тем же. Ясно, что солитон может быть сформирован при пиковой мощности начального импульса, превышающей пороговую величину.
1.0 0.8 0.6 0.4 0.2
2 n = 1, ß2 < 0
-5 0 5 0
T / To
1.6 1.2 0.8 0.4
-10
** 0
z|Lp
б)
Рис. 5. Изменение формы импульса, распространяющегося в среде с аномальной дисперсией: а) гауссовского при N = 1, б) секансного при N > 1
Рассмотрим случай, когда значение пиковой мощности не точно соответствует условию (9), полученному для оптического волокна без потерь (величина N, определенная по выражению (2) не является целым числом). Использование теории возмущений для решения НУШ показывает, что такой импульс, распространяясь по оптическому волокну, асимптотически преобразуется в солитон, порядок которого есть целое число N, ближайшее к начальному значению N. Часть его энергии при этом рассеивается.
Представим изначальное значение N в виде N = N + £, где |е| < 0.5, тогда (8) примет вид:
и(0,т) = (N + 2е) • sech[(1 + 2e/N) • т]. (11)
Длительность импульса в (11) будет равна: N
70 = ■=--70. (12)
0 N + 2e 0
Из (10) следует, что для фундаментального солитона (N = 1) длительность возрастает при е < 0 (при N < 0.5 солитон не образуется) и уменьшается при е > 0.
На рис. 5б показана динамика секансного импульса (11) с N = 1.2. Длительность импульса и его пиковая мощность сначала осциллируют, но постепенно стабилизируются до того момента, когда наконец параметры импульса будут удовлетворять условию N = 1.
2. Общие вопросы использования солитонов в волоконно-оптических системах связи
Работа высокоскоростных ВОСС с бинарной амплитудной модуляцией обычно ограничена ХД, из-
5
за которой импульс расширяется и возникает межсимвольная интерференция. Поскольку солитоны могут сохранять свою форму благодаря балансу между нелинейными и дисперсионными эффектами, их использование могло бы увеличить расстояние и скорость передачи ВОСС с амплитудной модуляцией. Однако, для практического использования солитонных ВОСС, необходимо рассмотреть эффекты, ограничивающие их возможности. Наиболее важными из них являются: 1) потери в оптическом волокне, 2) наличие частотной модуляции в начальном импульсе, 3) взаимодействие соседних импульсов.
На рисунке 6 изображен солитонный поток двоичных сигналов с RZ кодированием. Обычно расстояние между солитонами превышает их полную длительность по уровню 0.5 от максимума (FWHM) в несколько раз. Математически солитон имеет бесконечную длительность. На практике солитонное решение справедливо только при больших скважностях солитонных импульсов. Запишем выражение для скорости передачи [2]:
J-J 1 ^ o.oo
D = - = - ~ -,
ТВ 2Ч0Т0 4oTFWHM
(13)
где Тв - длительность битового интервала; 2цо = Тв / То - нормированное расстояние между соседними солитонами.
Из выражений (2) и (8) следует, что мощность фундаментального солитона равна [2]:
р(т) = Р0 • sech2(i),
(14)
где Ро определяется выражением (5). Энергия фундаментального солитона получается интегрированием мощности по внутриимпульсному времени Т = т-То [2]:
= J Р(Т)
dT = 2Po • To.
(15)
В предположении равной вероятности нулей и единиц в передаваемом сигнале средняя мощность Ш сигнала становится равной:
Ps = Es• (В/2) = Ро/2Чо.
(16)
Например, для солитонной ВОСС со скоростью передачи 10 Гбит/с получим To = 10 пс, если выбрать qo = 5. Длительность импульса по уровню 0.5 составит 17.6 пс. Пиковая мощность входного импульса составит 5 мВт при значениях ß2 = -1 пс2/км и у = 2 (Вт-км)-1, типичных для оптического волокна со смещенной дисперсией (DSF). Это значение пиковой мощности отвечает энергии импульса 0.1 пДж и уровню средней мощности всего лишь 0.5 мВт.
Tb = 1/B 4-
Солитон
110 1 0 1 Рис. 6. Последовательность солитонов при двоичном кодировании сигналов в RZ формате
3. Взаимодействие солитонов
Рассмотрим проблемы, возникающие из-за взаимодействия соседних солитонов в ВОСС. Наличие импульсов в соседних битах искажает временно'й солитон уже потому, что суммарное поле не является решением НУШ. Из-за нелинейного взаимодействия соседние солитоны либо сближаются (притягиваются), либо расходятся (отталкиваются). Для понимания процессов взаимодействия солитонов можно численно решить НУШ с начальными условиями в виде пары солитонов [2, 4]:
и(0,т) = sech(r — q0) + г • sech[r • (т + q0)] х х exp(i • 0),
(17)
где r - отношение амплитуд двух солитонов; 0 -относительная фаза; 2q0 - начальное нормированное расстояние.
С помощью численного решения НУШ можно определить расстояние между солитонами Цб в зависимости от расстояния 2 вдоль оптического волокна для начального расстояния Цо между солитонами. При Цо >> 1 можно использовать теорию возмущений. В частном случае при г = 1 и 0 = о текущее расстояние между солитонами Цб можно определить по выражению:
2 • ехр[2 • (д8 - оо)] = = 1 + соб[4 • г • ехр(-^0)].
(18)
Это соотношение показывает, что расстояние Цб^) между соседними солитонами периодически колеблется с периодом:
п
zp = — exp^o).
(19)
На рисунке 7 приведены результаты расчета зависимости периода взаимодействия двух соли-тонов 2р от начального расстояния Ц0 между ними.
6 7
Чо
Рис. 7. Зависимость периода изменения взаимодействия двух солитонов от первоначального расстояния между ними
Выражение (19) обеспечивает достаточную точность при Ц0 > 3 и может использоваться при проектировании солитонных ВОСС. Если величина
¿7шах, км В2 = -1 пс2/км
106 I
10
B, Гбит/с
а)
100
ZpLü много больше длины линии передачи Lt, взаимодействием солитонов можно пренебречь, так как расстояние между солитонами будет лишь незначительно отклоняться от своего первоначального значения. При qo = 6 получим Zp = 634. Положив дисперсионную длину равной Ld = 100 км, получим, что Lt << Zp^Ld даже при Lt = 10000 км.
C учетом (2), (13) и (19) условие Lt << Zp^Ld можно записать в виде:
п • exp(qo)
B2LT << ■
8 IP2Í
(20)
На рисунке 8 показаны рассчитанные по (20) зависимости максимальной длины солитонной линии, ограниченной взаимодействием соседних солитонов, от скорости передачи при различных Ц0 и ДГС. Видно, что ограничения на длину солитон-ной линии существенны только при высоких скоростях передачи и малых относительных расстояниях между солитонами.
LTmax, км
qo = 5
10
10
B, Гбит/с
б)
100
Рис. 8. Зависимости максимальной длины солитонной линии от скорости передачи при различных нормированных начальных расстояниях между солитонами (а) и при различных значениях ДГС (б)
1
1
4. Моделирование процессов распространения импульсов секансной формы по оптическому волокну
Для подтверждения изложенной теории было проведено моделирование ВОСС со скоростью передачи В = 10 Гбит/с на основе оптического волокна со смещенной дисперсией длиной 200 км со следующими параметрами: 02 = -1.9 пс2/км, а = 0, у = 2.6 (Вт-км)-1.
На рисунке 9 представлены последовательности входных и выходных импульсов, имеющих параметры до = 5, Т0 = 11.3 пс, Р0 = 6 мВт. Видно, что реализован режим распространения фундаменталь-
ных солитонов, т.к. длительности и амплитуды импульсов одинаковы. ХД полностью компенсируется фазовой самомодуляцией на всем протяжении оптического волокна.
Для полноты картины на рисунке 10 представлены последовательности входных и выходных импульсов при входной мощности Р0 = 3 мВт в два раза меньшей той, которая нужна для формирования фундаментального солитона. Видно, что преобладает ХД, выходной импульс расширился, а его амплитуда равная 1.7 мВт значительно меньше амплитуды входного импульса.
Р, мВт
У У У
0
1
2
3
Р, мВт
У У У
0
1
2
3
г / Тв г / Тв
а) б)
Рис. 9. Последовательности импульсов на входе (а) и выходе (б) оптического волокна в режиме распространения фундаментальных
солитонов
Р, мВт
3 2
1
У V У V У V
Р, мВт
0
1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0
3
1 2 3 0 1 2
г / Тв г / Тв
а) б)
Рис. 10. Последовательности импульсов на входе (а) и выходе (б) оптического волокна, мощность которых в два раза меньше необходимой для формирования фундаментальных солитонов
Р, мВт
12 10 8 6 4 2 0
Р, мВт
У V У V У
18 16 14 12 10 8 6 4 2 0
У У )
12
г / Тв
12
г / Тв
а) б)
Рис. 11. Последовательности импульсов на входе (а) и выходе (б) оптического волокна, мощность которых в два раза больше необходимой для формирования фундаментальных солитонов.
На рисунке 11 представлены последовательно- солитона. Видно, что преобладает фазовая самомости входных и выходных импульсов при входной дуляция, выходные импульсы сузились, а их ампли-мощности Р0 = 12 мВт в два раза большей той, кото- туда стала равной 18 мВт (значительно больше ам-рая нужна для формирования фундаментального плитуд входных импульсов, составляющих 12 мВт).
0
3
0
3
Заключение
Проведенный анализ показал, что:
- во-первых, в оптическом волокне с аномальной дисперсией (отрицательной ДГС) может быть сформирован солитонный импульс, сохраняющий или периодически восстанавливающий свою форму в процессе распространения по волокну; период изменения формы солитона увеличивается при возрастании его длительности и уменьшении ДГС волокна;
- во-вторых, для формирования солитона требуется специальная форма импульса и определенная пиковая мощность, однако высокая устойчивость солитона к возмущениям позволяет создать его, даже если форма импульса и пиковая мощность отклоняются от идеальных условий;
- и в-третьих, дальность связи солитонной линии Lтmax теоретически ограничивается фактом взаимодействия соседних солитонов при распространении по волокну, однако эти ограничения существенны только при высоких скоростях передачи и малых относительных расстояниях между солитонами.
Таким образом, можно утверждать, что использование солитонных линий позволяет обеспечить высокое качество связи на большие расстояния при сравнительно небольших мощностях источников излучения. В следующих частях будут рассмотрены проблемы, возникающие при учете потерь в оптических волокнах, и методы их решения. Авторы надеются, что данная публикация будет способствовать возрождению интереса к использованию оптических солитонов в телекоммуникациях.
Список используемых источников
1. Агравал Г. Нелинейная волоконная оптика. М.: Мир, 1996. 323 с.
2. Кившарь Ю.С., Агравал Г.П. Оптические солитоны. От волоконных световодов до фотонных кристаллов: Пер. с англ. под ред. Н.Н. Розанова. М.: ФИЗМАТЛИТ, 2005. 648 с.
3. Ахмедиев Н.Н., Анкевич А. Солитоны. М.: ФИЗМАТЛИТ, 2003. 304 с.
4. Linn F. Mollenauer, James P. Gordon. Solitons in Optical Fibers: Fundamentals and Applications. Academic Press, 2006. 278 p.
* * *
PROPERTIES OF TEMPORARY OPTICAL SOLITONS
IN OPTICAL FIBERS AND THE POSSIBILITY OF THEIR USE IN TELECOMMUNICATIONS. PART 1
E. Andreeva1, M. Bylina1, S. Glagolev1, P. Chaimardanov1
*The Bonch-Bruevich State University of Telecommunications, St. Petersburg, 193232, Russian Federation
Article info
Article in Russian
For citation: Andreeva E., Bylina M., Glagolev S., Chaimardanov P. Properties of Temporary Optical Solitons in Optical Fibers and the Possibility of their Use in Telecommunications. Part 1 // Proceedings of Telecommunication Universities. 2018. Vol. 4. Iss. 1. PP. 5-12.
Abstract: In the first part of the paper, a theoretical justification is given for the existence and stability of solitons in an optical fiber without losses. The possibility of using solitons in optical communication systems is shown. The problems associated with the interaction of solitons are considered. The results of modeling the processes of soliton propagation along an optical fiber are presented. In the second and third parts, we will consider the conservation of soliton regime in an optical fiber by controlling losses and dispersion in the linear path of a fiber-optic communication system.
Keywords: optical soliton, temporary optical soliton, single-mode optical fiber, fiber optic communication system, group velocity dispersion, self-phase modulation.