УДК 537.86
Е.М. Громов, Т.В. Наседкина, В.В. Тютин
СОЛИТОНЫ В РАСШИРЕННОМ НЕЛИНЕЙНОМ УРАВНЕНИИ ШРЕДИНГЕРА С ПРОСТРАНСТВЕННОЙ НЕОДНОРОДНОСТЬЮ ДИСПЕРСИИ
В ВИДЕ ПЕРЕПАДА
Национальный исследовательский университет Высшая Школа Экономики
Цель работы: В данной работе рассмотрена динамика солитонов в рамках расширенного неоднородного уравнения Шредингера с пространственной неоднородностью дисперсии в виде плавного перепада и при учете псевдо индуцированного рассеяния на затухающих НЧ волнах. Научный подход: Исследование проведено как численно, так и аналитически.
Результат: Показана возможность баланса псевдоиндуцированного рассеяния, приводящего к смещению пространственного спектра волновых чисел солитона в длинноволновую область, и убывающей дисперсии, смещающей спектр солитона в коротковолновую область. В явном виде получено солитонное решение, возникающее в результате этого баланса.
Новизна: Результаты исследования новые и могут иметь приложение для разработки новых поколений оптоволоконных линий связи на базе коротких оптических солитонов.
Ключевые слова: расширенное нелинейное уравнение Шредингера, солитон, вынужденное рассеяние, неоднородность, дисперсия второго порядка, аналитическое решение, численный эксперимент
Введение
Интенсивные волновые пакеты довольно большой протяженности способны распространяются в нелинейных средах без изменения формы (солитоны). Для решения многих прикладных задач существует необходимость уменьшения протяженности солитонов до нескольких длин волн. Подобное уменьшение сопровождается, как правило, индуцированным рассеянием на низкочастотных (НЧ) возмущениях среды, разрушающим стабильность распространения солитонов. До настоящего времени довольно подробно изучено индуцированное рассеяние на однородных по пространственной координате временных НЧ модах среды (индуцированное рассеяние Рамана) [1]. Данное рассеяние описывается в расширенном нелинейном уравнении Шредингера (НУШ) дополнительным слагаемым с временным запаздыванием нелинейного керровского отклика. Для солитонов учет данного рассеяния приводит к смещению их частот вниз по спектру [2] и, как следствие, к потере устойчивости и разрушению. Возможность компенсации рамановского сдвига частоты изучена к настоящему времени довольно подробно [3, 4].
Однако для ряда сред распространение коротких солитонов сопровождается возбуждением затухающих НЧ волн. В плазме это ионно-звуковые волны, возбуждаемые электромагнитными или ленгмюровскими солитонами; в стратифицированной жидкости - внутренние волны, возбуждаемые солитонами поверхностных волн. Учет волнового характера возбуждаемых НЧ возмущений существенным образом меняет динамику коротких солитонов высокочастотного (ВЧ) поля в подобных средах. Впервые модель для описания индуцированного рассеяния ВЧ волн на затухающих НЧ волнах, названное как псевдоиндуцированное рассеяние, предложена в работе [5], в которой получено модельное уравнение в виде расширенного НУШ с псевдоиндуцированным рассеянием и показано, что псевдоиндуцированное рассеяние приводит к сдвигу пространственного спектра волновых чисел солитона в длинноволновую область. До настоящего времени рассмотрены случаи дисперсии, убывающей в пространстве по линейному [5] и экспоненциальному [6,7] законам.
В данной работе динамика солитонов в рамках расширенного НУШ с псевдоиндуци-рованным рассеянием рассмотрена при пространственной зависимости дисперсии в виде перепада. Показана возможность компенсации псевдоиндуцированного рассеяния.
© Громов Е.М., Наседкина Т.В., Тютин В.В., 2014.
Исходное уравнение и интегральные соотношения
Рассмотрим динамику ВЧ волнового пакета и (%, t )ехр (/ш, — /к%) в расширенном НУШ с псевдоиндуцированным рассеянием и неоднородной дисперсией [5]:
ди д Г
2/-+
дt
ди
д(^)ди| + 2и|и|2 + д\
дМ
= 0.
(1)
где д(%) - линейная дисперсия второго порядка; д - коэффициент псевдоиндуцированного рассеяния. Последнее слагаемое в (1) описывает индуцированное рассеяние ВЧ волнового поля на затухающих НЧ волнах и является пространственным аналогом индуцированного рассеяния Рамана в оптике (псевдоиндуцированное рассеяние). Уравнение (1) при нулевых условиях на бесконечности и^^ ^ 0 имеет следующие дифференциальные соотношения
для интегральных моментов волнового пакета:
0 ,2
1 го»
21,1 ки\2 *=—д1
й 1 эк) ■
1 |и| 0% = 0,
ГО 7
1
— ГО ~
и д%
+ к2 \и\2
7« 7 ГО +ГО
31и2 2 ,
(2)
(3)
(4)
где и = р\ ехр (/ф), К = дф / д% - добавочное волновое число пакета. Масштабы неоднородности дисперсии и волнового числа к примем много большими масштаба неоднородности огибающей волнового пакета ~ОчК >> . В этом случае, полагая добавочное волновое число не зависящим от пространственной координаты [5], К (%,, ) = к (,), система (2)-(4) для со-литоно подобных волновых пакетов сводится к виду:
2
Л
дЬ0 / П — д'(ф0 / п2 — д'(%)ык2,
йп Ж
где д'(%)=(0д / ,
п = д
(%)/ ?0
?0
= д(0); ¿0 = |[д(и (%,, = 0)2 )
)/ д%
(5)
0%,
го
= 1 [д|и(%,, = 0) / д%]2 0% - начальные интегральные моменты волнового поля. Система
—го
(5) имеет состояние равновесия
к = Д = Д* = — д'(%0 ^0 / ¿0 , (6)
при котором псевдоиндуцированное рассеяние компенсируется убывающей дисперсией д'(%0)< 0 : при Д = Д* начальный волновой пакет при , > 0 остается неизменным.
Солитонное решение
Рассмотрим решение уравнения (3) в виде стационарной волны и (%,, ) = у(%)ехр (/О,) при пространственном профиле дисперсии в виде перепада д(%) = д0 [1 — 1апЬ(% / Ъ)]:
д0 [1 — 1апЬ(% / Я)]
%
0у
0 у
Ъ собЬ2 (% / Я) 0%
+ 2у3 — 2Оу + ду
0 (у2) 0%
= 0.
(7)
—ГО
2
—го
—ГО
—ГО
ГО
— ГО
Масштаб неоднородности дисперсии примем много большим масштаба неоднородности огибающей волнового пакета Б >> Ьу . С учетом малости параметра в ~ Ьу / Б ~ ц << д0 и
аппроксимации 1апЬ(£, / / Б, решение уравнения (7) будем искать в виде у = у0 +у!,
где у! ~ еу 0 << у0. Удерживая члены порядка в , имеем
%
d Vo
d 2ш, % d^ +
dt,
(ба^0 - 2о)
2 +2^0-2Qy0=0,
D dt? 3 dt D dt
(8)
(9)
Уравнение (8) имеет классическое солитонное решение у0 = А0/ ео8Ь(Е, / А), где А = л[% / А, О = А2 /2. Уравнение (9) после замены ^ = / А и ^ = /(А) А) примет вид
' 6 ^
d
dr\ ^ cosh ц
-1
¥ — 2-
ц
ц
5 ц sinh ц sinh ц
(10)
cosh3 ц cosh ц 4 ц* cosh4 ц cosh2 ц
где ц* = -5%o /(8 A D) - величина псевдоиндуцированного рассеяния, отвечающая состоянию равновесия системы (5). При ^(o) = 0 уравнение (10) имеет точное решение
I ц2 Ц I 1 ¥(ц) = ¥'(0)ц--tanhцн--(tanhц) ln(cosh ц) sech^--
4ц*
12
-Ц-1 ц*
(tanh2 ц) sinh ц. (11)
Рис. 1. Распределение ¥(ц) при ц — ц* и различных ¥'(0)
При д = д*, отвечающем балансу псевдоиндуцированного рассеяния и убывающей дисперсии, решение (11) локализовано. Асимптотика (11) в этом случае на больших значениях аргумента следующая: ¥(| ^ ±го)« ±|2 ехр (—||). На рис. 1 показано распределение волнового поля ¥(|) при д = д* и различных значений ¥'(0).
Решение ¥(|) антисимметрично. Солитоны с антисимметричными «хвостами» возникают и в хорошо известной системе линейно связанных НУШ.
При д ф д* решение (11) не локализовано |¥(| ^ ±го) ^ го . На рис. 2 приведено распределение функции ¥(|) при начальных условиях ¥'(0) = 0 и различных значениях д.
Рис. 2. Распределение ¥(|) при ¥'(0) = 0 и различных д
Численное моделирование
Рассмотрим численно начальную задачу динамики пакета и (%,, = 0) = БесЬ% в рамках (1) при д(%) = 1 — 1апЬ(%/10) и различных значениях д . Величина равновесного значения псевдоиндуцированного рассеяния для выбранного профиля неоднородности дисперсии и данного начального импульса из аналитического соотношения (6) составляет д* = 1/16. На
рис.3 приведены численные результаты пространственно-временного распределения Ц при
различных д. При д = 1/16 = д (рис. 3, а) исходный волновой пакет не меняется во времени, что соответствует состоянию равновесия псевдо индуцированного рассеяния и неоднородной дисперсии. Отклонение параметра д от равновесного значения д* приводит к изменению параметров солитона (волновое число, амплитуда) во времени. На рис. 3, б показано пространственно-временное распределение модуля огибающей и\ при д = 1 / 32 = д* /2.
На рис. 4 приведена динамика во времени локального волнового числа в точке максимума модуля огибающей волнового пакета при различных значениях д. Непрерывные кривые отвечают результатам численного счета, пунктирные - аналитическое решение системы (5).
Результаты численного счета близки к аналитическим. Небольшое различие связано с
малыми полями излучения из ядра солитона, приводящими к уменьшению эффективного параметра псевдо индуцированного рассеяния.
а) б)
Рис. 3. Численные результаты пространственно-временного распределения | и\
при различных ц: ц = 1/16 = (а); - ц = 1/32 = ц*/2 (б)
Рис. 4. Численные и аналитические результаты (сплошные и пунктирные кривые) для локального волнового числа в точке максимума огибающей волнового пакета в зависимости от времени t при различных значениях ц
Заключение
В данной работе рассмотрена динамика солитонов в рамках расширенного НУШ с псевдоиндуцированным рассеянием и пространственно неоднородной дисперсией в виде перепада. Показана возможность баланса псевдоиндуцированного рассеяния, приводящего к смещению спектра волновых чисел солитона в длинноволновую область, и убывающей в пространстве дисперсии, смещающей спектр солитона в коротковолновую область.
Данное научное исследование (№ проекта 14-01-0023) выполнено при поддержке Программы «Научный фонд НИУ ВШЭ» в 2014/2015 гг.
Библиографический список
1. Kivshar, Y.S. Optical Solitons: From Fibers to Photonic Crystals / Y.S. Kivshar, G.P. Agraval // Academic, San Diego, 2003.
2. Gordon, J.P. Theory of the soliton self-frequency shift // Optics Letters 11 (1986) 662-664.
3. Malomed, B.A. Soliton Management in Periodic Systems / B.A. Malomed. Springer, New York, 2006.
4. Yang, Y. Solitons in Field Theory and Nonlinear Analysis / Y. Yang. Springer, New York, 2001.
5. Gromov, E.M. Soliton dynamics in an extended nonlinear Schrödinger equation with a spatial counterpart of the stimulated Raman scattering / E.M. Gromov, B.A. Malomed // J. Plasma Physics 79 (2013) 1057-1062.
6. Gromov, E.M. Damped solitons in an extended nonlinear Schrödinger equation with a spatial stimulated Raman scattering and decreasing dispersion / E.M. Gromov, B.A. Malomed // Opt. Comm. 320 (2014) 88-93.
7. Aseeva, N.V. The short envelope soliton dynamics in inhomogeneous dispersive media with allowance for stimulated scattering by damped low-frequency waves / N.V. Aseeva, E.M. Gromov, V.V. Tyutin // Radiophysics and Quantum Electronics 56 (2013) 157-166.
Дата поступления в редакцию 22.04.2014
E.M. Gromov, T.V. Nasedkina, V.V. Tyutin
SOLITONS IN AN EXTENDED NONLINEAR SCHRÖDINGER EQUATION WITH A SPATIAL KINK-TYPE DISPERSION DISTRIBUTION
National Investigate University Higher School of Economics
Purpose: Dynamics of solitons is considered in the framework of the extended nonlinear Schrödinger equation
including a pseudo-stimulated-Raman-scattering term, i.e., a spatial-domain counterpart of the stimulated-Raman-
scattering term, which is a known ingredient of the temporal-domain nonlinear Schrödinger equation in optics. Also
included is inhomogeneity of the spatial second-order dispersion in kink-type form.
Approach: Soliton's dynamic investigated as analytically as numerically.
Findings: An analytical soliton solution with a permanent shape is found in an approximate form.
Originality: The obtained results is original and can be important for optical fibers application.
Key words: Extended Nonlinear Schrödinger Equation, Soliton Solution, Stimulated Scattering, Inhomogeneity, Second-Order Dispersion, Analytical Solutions, Numerical Simulation