CC BY
36
ПРИБЛИЖЕННЫЕ ПАРАМЕТРИЧЕСКИЕ МОДЕЛИ ДИНАМИКИ САМО ВОЗДЕЙСТВИЯ ИМПУЛЬСОВ В НЕЛИНЕЙНЫХ ОПТИЧЕСКИХ СРЕДАХ С МОДОВОЙ ДИСПЕРСИЕЙ
Введение
Возрастающий интерес к анализу процессов самовоздействия импульсов когерентного излучения в нелинейных оптических средах связан с возможностями их использования ДЛЯ получения импульсов особо малой (фемтосекундной) длительности, повышения скоростей передачи информации по волоконно-оптическим каналам, создания различных новых оптических приборов, элементов оптической обработки информации со сверхвысоким быстродействием, а в перспективе - оптических компьютеров [1-9]. Доказана возможность самостабилизации импульсов и образования солитонов не только в одномодовых, но и в многомодовых оптических волокнах [2, 10, 11]. Хотя для создания высокоскоростных систем дальней оптической связи многомодовые волокна мало пригодны (главным образом потому, что из-за большого диаметра сердцевины для получения в них солитонных режимов требуются мощности на несколько порядков выше, чем в одномодовых волокнах), во многих других нелинейных оптических системах необходимо принимать во внимание модовую дисперсию. В первую очередь это относится к различным устройствам оптической обработки информации [12, 13]. Кроме того, и в одномодовых световодах при технологических нарушениях симметрии сечения волокна и неоднородностях появляются различия постоянных распространения двух ортогонально поляризованных мод (двойное лучепреломление), приводящие к модовой дисперсии [14]. Сходные задачи (но уже для разных частот, а не постоянных распространения) возникают при анализе вынужденного комбинационного рассеяния, которое предполагается использовать для компенсации потерь в волокне и увеличения дальности солитонных систем оптической связи [3, 4, 9].
Основу для решения задач анализа и расчета указанных оптических систем и устройств дает общая теория нелинейной динамики интенсивных сверхкоротких импульсов в многомодовых световодах, развитая в последние годы [15-17]. В общем случае уравнения эволюции модовых компонент импульса оказываются сложными и не допускают аналитического решения. Однако в большинстве практических задач расчета оптических систем (в частности, волоконно-оптических систем связи и устройств обработки информации в цифровой форме) нет необходимости в детальном анализе эволюции формы импульсов, достаточно знать характер изменения их отдельных параметров. В таких случаях анализ динамики самовоздействия может быть выполнен на основе гораздо более простых параметрических моделей, некоторые из которых рассматриваются ниже.
1. Нелинейные уравнения эволюции импульсов в многомодовом световоде
В общем случае такие уравнения должны быть записаны в векторной форме и требуют учета изменений поперечной структуры поля [15-17]. Если, однако, ограничиваться учетом только тех нелинейных явлений, которые связаны с зависимостью показателя преломления от интенсивности поля, и при этом не учитывать инерционность его нелинейной части, а также ряд других нелинейных эффектов, в том числе дисперсию выше второго порядка, то можно использовать скалярные модели, в которых взаимосвязь мод определяется только зависимостью показателя преломления от их суммарной интенсивности [10]. Такие упрощенные модели достаточно хорошо описывают динамику самовоздействия небольшого числа мод (наиболее точно - двух мод [17]).
При указанных условиях и обычных допущениях метода медленно меняющихся амплитуд процесс распространения импульса когерентного излучения в регулярном световоде с кубической нелинейностью поляризации описывается системой уравнений Шредингеровского типа для N взаимосвязанных ортогональных мод [10]
(/я = 1,2,...,ЛО
где ц=х/Ь , *=Ц-Ф0)/Т0 - продольная пространственная и смещенная временная
координаты^ нормированные с некоторыми масштабами по длине времени ?(Г0) и
скорости (у0); - комплексная огибающая импульса т-ой моды, нормированная на ее максимальное по модулю значение на входе световода (Дэ)> Ут> Кт1
коэффициенты, характеризующие дисперсию, затухание, относительную скорость взаимного смещения и нелинейность по каждой моде и определяемые выражениями
Здесь уш, Ит, рт - групповая скорость, коэффициент затухания и постоянная распространения т-ой моды (в размерных единицах), Вт, - коэффициенты,
определяемые в результате усреднения т-ой и /-ой мод по формулам, приведенным в [2], п0, /*2 - невозмущенное значение показателя преломления и коэффициент в его нелинейном члене.
При описании одномодовот режима (N=1) обычно выбирают в качестве Ь0 "дисперсионную длину" соответствующую удвоению ширины импульса в линейном режиме, в качествее Т0 - его начальную полуширину, у0 - 1рупповую скорость волн. Тогда в эволюционном уравнении у=0, а=1 [2]. Для системы нескольких уравнений вида (1) выполнить такое масштабирование в общем случае невозможно (хотя для двух мод можно учесть различие в скоростях с помощью дополнительных фазовых множителей, а в уравнениях положить У1=У2=® [11])-
Система уравнений (1) описывает процесс распространения импульса в световоде при начальных условиях
где ч/то(т) - комплексные огибающие импульсов, возбуждаемых на входе световода по отдельным модам.
К аналогичному виду приводятся уравнения динамики двухчастотного взаимодействия импульсов [7], а также импульсов в волокнах с двойным лучепреломлением (с тем отличием, что перекрестные члены уравнений содержат и линейные слагаемые) [14].
Точное аналитическое решение уравнений вида (1) затруднено, а численное решение требует больших затрат машинного времени. Однако в большинстве технических задач, связанных с расчетом нелинейных оптических элементов систем обработки и передачи информации, знание полной картины эволюции импульсов, обеспечиваемое такими решениями, как уже отмечено выше, не требуется, достаточно знать некоторые общие параметры импульсных сигналов. Уравнения их эволюции могут быть выведены из исходных уравнений (1) на основе автомодельных приближений, рассматриваемых ниже.
Импульсы ц/т(г|,т) , являющиеся автомодельными решениями уравнений (1), обладают свойством сохранять свою форму неизменной в процессе распространения по волокну, т.е. при всех значениях т| . Получение такого решения позволяет свести задачу анализа преобразования импульса к определению эволюции его параметров.
Существует два частных вида уравнений (1), для которых известны точные решения такого рода. В линейном режиме, когда все *^=0, система (1) распадается на N не связанных между собой линейных уравнений отдельных мод, для каждого из которых при начальном воздействии (3) гауссовской формы решение также имеет вид гауссовского импульса с измененными параметрами. Второй частный случай - одномерное (Л^!)
(2)
Ч'т(°’т) = Уто^)’ т= 1,2,—,ЛГ;
(3)
Автомодельные решения эволюционных уравнений
нелинейное уравнение Шредингера (НУШ), в котором (3=0, а а и к имеют критические значения, соответствующие распространению солитона (например, а=1 и к=2). Последний, как известно, имеет форму гиперболического секанса и без изменения этой формы распространяется по волокну.
В остальных режимах импульсы указанного вида можно рассматривать лишь как некоторые естественные приближения к точному решению уравнений (1). В [18] получено и исследовано автомодельное решение гауссовской формы для одномерного НУШ. Показано, что такое приближение дает в целом правильную качественную картину эволюции импульса, позволяет выявить ее основные режимы и рассчитать изменения параметров. В [19] эти результаты обобщены на многомодовые световоды, описываемые системами НУШ вида (1).
Однако, естественно ожидать, что для нелинейных режимов, близких к солитонному, более точным окажется автомодельное приближение в форме гиперболических секансов. Рассмотрим такое решение для системы НУШ (1) и сравним его с гауссовским приближением.
Пусть на входе световода, описываемого уравнениями (1), действует импульс когерентного излучения, комплексные огибающие которого по отдельным модам имеют вид
Ч'то(х)=4поДт; аШ0)ехр[Ю(х;
^то> сшо)1
(4)
где Лто - начальная комплексная амплитуда импульса т-ой моды,
/7(х;йя0) = 8єсЛ
/ т
(5)
- модуль его нормированной огибающей,
®(т’ ^то> сто)-^тох2+стох ^
- начальная фаза.
Здесь ато, £т0, сто - начальные значения полуширины, коэффициенты линейной частотной модуляции (ЛЧМ) и сдвига средней частоты импульса т-ой моды (при одночастотном воздействии следует положить все сто=0).
В полученном ранее автомодельном решении гауссовской формы [19] функция Дх;ато) была задана в виде
2
(7)
ехр
т0 '
Приближенное решение уравнений (1) при начальных условиях (3), в которых ч/то(т) определены (4), будем искать в форме функции такого же вида, но с параметрами, зависящими от г\ :
ЯтМ^хрМт-^тО; Мп); СщМЬРл} • (8)
Их начальные значения определяются параметрами входного воздействия:
Ап^)=Апо’ дто> ^то’ ст^)~Сто> * (^)
Объединим входящие в (8) неизвестные функции в блочный вектор
Дп)=№(л)Д2(т1)»-А(л)}» 0°)
где
(11)
Для их определения воспользуемся вариационным методом, который базируется на гамильтоновом представлении уравнений динамики [21] и широко испол^уется в теории нелинейных волн [22]. В рассматриваемой здесь задаче аналогично [18, 19] наилучшее (в
смысле принципа Гамильтона) приближение автомодельного представления (8) к точному
решению уравнений (1) достигается при выборе вектора Дц), удовлетворяющего
вариационным уравнениям
— = 0; к = 1,2,...,6# (12)
іде - компоненты векторов (11), а
Цф,ф;,ф;,)л из)
- лагранжиан с плотностью Л(Ф,Ф х,Ф л), зависящей от вектора функций
ф(лд)={фі(л/0^і(л/0,...,фм(л>*)>ф^1і>'0}
И его Производных ПО Т и Т], Фт и Ф^.
Здесь
Фт(ти'0=Ч'т(тьФхР(Ртті)
(14)
Для уравнений (1) с учетом подстановки (14) плотность лагранжиана, как показано [19], представима в виде
N ф,ф;.«> =£• 1 / ( • > Л ар» • ^ф- Ф-Т'-Ф- — / +~Ут ' * асР» Ф, „ ~Ф ” +ат
\ . » м/ **=1 .4 2 2 ^ дх дх , дх
-^ІХкҐкҐ ехр(-2м)
1 м
(15)
Подставляя функцию (14) с учетом (5), (6) и (8) в (15) и выполнив интегрирование согласно (13), получаем
/ . \ ЛАт . йАт
л.—-'*.—*
. |2| х 2 с1Ь аа 2ат 2 , ,
^ ІТа-^^-V гл+1?+7“-‘,а+
<*Ч
2а_ 2
2 а с
т т
(16)
1«есЛ:
'-я.
У а-
эесЛ3
Х~Я,
СІХ
я| /
После определения вариационных производных выражения (16) по входящим в него неизвестным параметрам, зависящим от т| , подставновки их в (12) и несложных преобразований получается следующая система обыкновенных дифференциальных уравнений для указанных параметров модели (8):
с!а_
<іт\
Г \
к
<1г\
•=4<МА ;
■ = а.
<-4Ь2я -Єхр(-2Рвті)2кя/Ж//(Ая,а/^ж>^) ;
Ы\
СІС
АЯ
; т=ъг,...,ы ■
где
Аат >а! >ЧШ,Ч)= 2 7 * і
п а_а
Щ ~ И;оІ2аіо і. * = 4/л2
\ Г \I / >
х'Чт
весЛ2
*-9.
вес Л'
( \ Х~Я,
< аі )
Ах
(17)
(18)
(19)
(20)
(21)
(22)
При этом
2ата,
т
г \ х~я,
5есЛ
т“9.
Бес/г'
4.М-
ИГ
а«М
(23)
(24)
Полученная система из 4N нелинейных обыкновенных дифференциальных уравнений (17)-(20) описывает эволюцию параметров импульса секанс-гиперболической формы (5), аппроксимирующего точное решение системы НУШ (1). Сравнивая ее с аналогичной системой, полученной в [19] для гауссовского импульса (7), нетрудно заметить, что они отличаются только коэффициентом к и видом функций в правой части. Для гауссовской модели к=1,
2 2 ( 2 ^
А*т >а1^т,я,) = °т1 Чт1 5 еХР _ Чт1 2 ъ 3 ►а II (-¥
а,п1 - ат1) “1 V
(25)
ще а2т]=а-т+а-, ;
При использовании секанс-гиперболической модели в правую часть уравнений (18) и (19) входят функции, аналитические выражения которых, если выполнить интегрирование в (22) и (23), оказываются очень громоздкими (и поэтому здесь опущены). При моделировании на ЭВМ их целесообразно определять путем численного интегрирования. В отличие от этого для гауссовской модели аналогичные функции, как видно из (25), имеют сравнительно простой вид, что ускоряет моделирование. Более подробное сравнение указанных двух моделей рассматривается ниже.
Результаты и выводы
*
Рассмотренный подход к анализу процессов самовоздействия импульсов в нелинейных оптических средах с модовой дисперсией позволяет перейти от их исходной модели в форме системы нелинейных уравнений в частных производных вида (1), т.е. в бесконечномерном пространстве состояний, к приближенному представлению в конечномерном пространстве параметров в виде системы обыкновенных дифференциальных уравнений (17)-(20). Свойства таких систем детально исследованы и известны многочисленные методы их приближенного и численного решения. При моделировании их на ЭВМ затраты времени и требования к объему памяти по крайней мере на порядок меньше по сравнению с исходной системой уравнений (1). Это особенно важно при статистическом моделировании с большим числом испытаний, которое необходимо при оценках вероятностей ошибочного приема цифровых сигналов и других статистических характеристик солитонных систем [20].
Полученные уравнения отражают и физический механизм нелинейного взаимодействия мод: за счет нелинейности развивается самомодуляция импульсов, проявляющаяся в изменениях как коэффициентов ЛЧМ Ьт, так и средних частот мод (смещении их на ст) в соответствии с (18) и (19), а эти параметры в свою очередь влияют на ширину ат и взаимное смещение <7т компонент импульса по отдельным модам согласно уравнениям (17) и (20).
Зависимости характера эволюции указанных величин от параметров нелинейности при использовании секанс-гиперболической аппроксимации аналогичны тем, которые были исследованы ранее для гауссовской модели [19].
При расчетах преобразований импульсных сигналов в нелинейных оптических средах с модовой дисперсией, как правило, необходимо знать параметры суммарного импульса, а не его компонент по отдельным модам (исключение составляют лишь о ласти применения, где используется селекция мод). Наиболее важным из таких параметров
является длительность (ширина) импульса, поскольку именно от нее зависят предельно достижимые скорость передачи информации, быстродействие и т.п. характеристики систем и устройств, в которых используются солитонные импульсы.
Поскольку рассматриваемые здесь импульсы не являются финитными, их необходимо характеризовать некоторой эквивалентной шириной (или полушириной). Из различных известных определений последней в данном случае наиболее удобной для расчетов является среднеквадратическая (по интенсивности) полуширина, определяемая как
где
~1(т-то)2Мт1^)|2^
(26)
= 1Мт1д)Г
с1х
- энергия импульса,
г/т
- среднее значение г .
Для одиночного импульса секанс-гиперболической формы (5) среднеквадратическая полуширина по интенсивности т8=а0л/2^3 для гауссовского импульса (7) т8=а0/>Й.
Применяя это определение к импульсу |у(т],т)|2 , образованному путем сложения интенсивностей всех мод, каждая из которых описывается выражением вида (8), получаем
х2 = т=1
5>.
2>.<
2>„
(27)
где для гауссовской модели у=0,5 , для секанс-гиперболической \>=л2/12.
В целях сравнения результатов, полученных на основе рассмотренных параметрических моделей и прямою решения исходных НУШ (1), было проведено моделирование на ЭВМ. Уравнения (17)-(20) решались численно методом Рунге-Кутга четвертого порядка, а система уравнений (1) - методом расщепления по физическим факторам с применением быстрого преобразования Фурье. Результаты моделирования показывают, что форма модовых компонент импульса в режимах, близких к солитонному, при использовании секанс-гиперболического приближения несколько лучше соответствует точному решению, чем при гауссовской модели. Для суммарного импульса указанное различие меньше.
На рис.1 показаны графики эволюции среднеквадратической полуширины т5(г|), определяемой по формуле (27), для импульса, распространяющегося в двухмодовом световоде, полученные на основе секанс-гиперболической (сплошные линии) и гауссовской (штриховые линии) моделей при различных значениях параметра нелинейности к . В этом примере все = к , а остальные параметры имели следующие значения:
= 1/^; а„ = 1‘> г, = -г2 =2; ьт<>=сто = =0; л*0=2л/з/л;
для гауссовского импульса ато=л/2.
Как видно из этих зависимостей, ширина импульса вследствие дисперсионных эффектов (материальной, волноводной и модовой дисперсии) в линейном режиме (к=0) быстро возрастает вдоль волокна. С увеличением мощности импульса и, соответственно, параметра к , это дисперсионное расширение замедляется благодаря нелинейным эффектам: частота и групповая задержка каждой модовой компоненты импульса в соответствии с уравнениями (17)-(20) изменяются за счет интенсивности второй моды (в общем случае - нескольких других мод) так, что волновой пакет начинает сжиматься, а дисперсия его - компенсироваться. Значение параметра к соответствующее самостабилизации импульса за счет полной компенсации дисперсии (т.е. образованию
солитона) в этом примере равно приблизительно 8. При к > 8 действие нелинейных эффектов еще сильнее и вызывает сжатие импульса, так что его общая ширина осциллирует в результате конкуренции нелинейных и дисперсионных эффектов.
Сравнение сплошных и штриховых кривых показывает, что различия в эволюции секанс-гиперболических и гауссовских импульсов с одинаковой начальной среднеквадратической полушириной незначительны. Таким образом, при прочих равных условиях гауссовское приближение может оказаться предпочтительным, как более простое. Разумеется, главным фактором, который следует принимать во внимание при выборе модели, является реальная форма входного импульса.
Аналогичные графики эволюции полуширины суммарного импульса двух мод, полученные при тех же значениях всех параметров путем прямого численного решения уравнений (1) указанным выше методом, представлены на рис.2. Сравнение их с рис.1 обнаруживает их достаточно хорошее соответствие.
В целом полученные результаты показывают: в тех случаях, когда детали эволюции формы импульсов (такие, как появление боковых лепестков, "пьедесталов", деформация . вершины и т.п.) не существенны для работы исследуемых оптических систем и устройств, при их расчетах с успехом могут быть использованы рассмотренные параметрические модели.
Рис.2
ЛИТЕРАТУРА
1. Хасэгава А., Кодама Ю./fТИИЭР, 1981, т.69, N 9, стр.57-63.
2. Сисакян И.Н., Шварцбург А.В.//Квантовая электроника, 1984, т.11, N 9, стр.1703-1721.
3. Беланов А.С., Головченко Е.А., Дианов Е.М., Никонова З.С., Прохоров А.М., Серкин Е.Н. //Труды ИОФАН, 1986, т.5, стр.35-59.
4. Ахманов С А., Выслоух В А., Чиркин А. С. Оптика фемтосекунды лазерных импульсов.-М.: Наука, 1988.
5. Blow K.J., Doran N.J., Nayar В.К.//Opt.Lett., 1989, v.14, pp.754-756.
6. Islam M.N.//Opt.Lett., 1989, v.14, N 22, pp.1257-1259.
7. Торчигин В.П.//В сб."Вычислительные машины с нетрадиционной архитекгурой’УПод ред.В.С.Бурцева. Вып.1.- М.: Наука, 1990.
8. Новые принципы оптической обработки информации/Под ред. СЛ.Ахманова и МА.Воронцова. -М.: Наука, 1990.
9. Гордон Г.И., Заркевич ЕА., Мишнаевский ПА., Оввян П.П., Смирнов ^.//.//Электросвязь, 1992, N 12, стр.14-17; 1993, N 2, стр.11-13.
10. Hasegava А.ЦOpt.Lett., 1980, v.5, N 10, рр.416-417.
11. Выслоух ВА.ЦУФНб 1982, т.136, N 3, стр.519-531.
12. Зуев МА., Молостов А А., Соколов В.Н., Торчигин В.П., Шварцбург ^.^.//Компьютерная оптика: Сб.МЦНТИ. -М., 1992, вып. 12, стр.51-56.
13. Зуев МА., Колчанов И.Г., Шварцбург А5.//Компьютерная оптика: Сб.МЦНТИ. -М., 1992, вып.12, стр.56-60.
14. Blow К.J., Doran N.J., Wood ^.//Opt.Lett., 1987, v.12, N 3, pp.202-204.
15. Колчанов И.Г., Шварцбург у4.£.//Компьютерная оптика: Сб.МЦНТИ. -М., 1989, вып.5, стр.95-100.
16. Водяницкий С.Я., Зуев М.А., Шапинский В.В., Шварцбург А.Б.//Компьютерная оптика: Сб.МЦНТИ. -М., 1989, вып.6, стр.37-42.
17. Зуев МА., Шварцбург ^.^.//Компьютерная оптика: Сб.МЦНТИ. -М., 1990, вып.7, стр.72-81.
18. Anderson D.//Phys.Rev.A: Gen.Phys., 1983, v.27, N 6, pp.3135-3145.
19. Широков C.M.//B сб.'Тадиотехнические и оптические системы связи", ЛЭИС, 1988, стр.27-34.
20. Широков С.М.//Электросвязь, 1992, N 3, стр.6-8.
21. Ландау ЛД., Лифшиц Е.М. Теория поля. -М.: Наука, 1973.
22. Заславский Г.М., Сагдеев Р.З. Введение в нелинейную физику. -М.: Наука, 1988.
