Автомодуляция одномерных волн на основе нелинейного уравнения Шредингера с некерровской нелинейностью

Тип работы:
Реферат
Предмет:
Физико-математические науки


Узнать стоимость

Детальная информация о работе

Выдержка из работы

АВТОМОДУЛЯЦИЯ ОДНОМЕРНЫХ ВОЛН НА ОСНОВЕ НЕЛИНЕЙНОГО УРАВНЕНИЯ ШРЕДИНГЕРА С НЕКЕРРОВСКОЙ НЕЛИНЕЙНОСТЬЮ
Иван Васильевич Алименков (доцент кафедры прикладной математики, e-mail: i-alimenkov@mail. ru) Самарский государственный аэрокосмический университет имени академика С.П. Королева
Аннотация
Показано, что нелинейное уравнение Шредингера с некерровской нелинейностью имеет решение в виде локализованного импульса, движущегося с постоянной скоростью без дисперсионного уширения. Данное решение найдено прямым методом, основанном на теории гамильтоновых систем, и содержит в себе, как частный случай, известное решение кубичного нелинейного уравнения Шредингера.
Ключевые слова: нелинейное уравнение Шредингера, нелинейность пятой степени, теория гамильтоновых систем, канонические преобразования, уравнение Гамильтона-Якоби, солитонные решения для степенной нелинейности.
Введение
Перечень физических приложений стандартного нелинейного уравнения Шредингера (НУШ)
idyl dt + Э 2yl dx2 + 2hy = 0, описывающего в
общем случае эволюцию огибающей несущей квазимонохроматической волны в слабо нелинейной системе, и связанных с ним нелинейных уравнений в настоящее время чрезвычайно обширен, и вряд ли можно сомневаться в их физической значимости. Нелинейный кубический член \ в
различных физических моделях, описываемых этим уравнением, возникает обычно из степенно —
го разложения некоторой физической величины
2
nNL, зависящей от интенсивности I = y. Это разложение имеет вид:
n
(I) = n21 + n412 +… = n 2y& quot- + n 4y4 +…
2
NL *
К примеру, в нелинейной оптике пщ^ - нелинейная часть показателя преломления п = п0 + пЩц. Стандартный безразмерный вид НУШ, приведенный выше, записан при учете наинизшего члена разло-
При этом коэффициент
жения n
¦¦n 21 = П2 У
нелинейности г пропорционален п2.
Функция нелинейного отклика системы яж (7) на внешнее воздействие несущей квазимонохро -матической волны в общем случае имеет сложный вид, определяемый конкретным физическим механизмом взаимодействия системы с полем несущей волны. Для нахождения явного вида яж (/) часто требуется квантовомеханический расчет, не позволяющий определить аналитическую зависимость функции отклика от интенсивности в широ -ком диапазоне. Функция отклика должна обладать двумя очевидными свойствами, а именно: обращаться в ноль при I = 0 и выходить на насыщение при I & gt->-1. Степенное разложение, указанное выше, применимо при малых значениях интенсивно -сти, но даже в этом случае, если ось I является правой касательной к графику функции пж (I) в
нуле, то разложение начинается только с члена n412 = n4y 4, что приводит к НУШ
idyl dt +d 2yl dx 2 + 2r\ = 0
с нелинейностью пятой степени.
Чтобы включить в решение и керровскую нелинейность третьей степени, рассмотрим более общее уравнение
idyldt + d2yldx2 + 2Гу = 0.
(1)
При V = 1 уравнение (1) очень хорошо изучено [1] и его решение имеет вид
ехр {[их/2 + (и2-ь2)/4 + ]}
y = A-
ch — (x — x0 — ut)
где и = 2А^. Здесь А, ь, х0,^0- свободные параметры.
Основной формализм Целью данной работы является решение уравнения (1) при произвольном V & gt- 0 (не обязательно целом) прямым методом, основанном на теории дифференциальных уравнений в частных производных первого порядка и теории гамильтоновых систем. Подстановка полевой функции у вида
у = Лехр{(хь/2 -1д + ф0)}, (2)
где ь, 3, ф0 — свободные параметры, в (1) приводит после отделения мнимой и вещественной частей к двум уравнениям
Ы Ы п
+ ь = 0, (3)
от ах
2IЦ — {^ I = a 212 — 8hV+2. dx2 I dx J '-
где
= u2 — 4d
(4)
(5)
Уравнение (3) является линейным однородным уравнением первого порядка с детально разработан-
2
a
2009
Компьютерная оптика, том 33, № 3
ной теорией [2]. Как известно из теории таких уравнений, общим решением (3) является любая дифференцируемая функция I = и («(х, /)}, где s (x,/) — левая часть первого интеграла уравнения характеристик, имеющая, как легко проверить, вид
5 = X — Х0.
Подстановка I = и («) в (4) приводит к обыкновенному дифференциальному уравнению второго порядка
…~ (6)
преобразование
2u (s)u '-(s) — u'- 2(s) = a 2 u 2(s) — 8ruv+2(s)
Упрощающее масштабное
(2 A a
1/v
f (t), t = as приводит (6) к виду
2f f — f 2 = f 2 — 8f v+2,
(7)
где / = 4/(г)/4г.
Как легко проверить, (7) следует из нормальной системы гамильтоновых уравнений
/ = дН / др/, р/ = -дН / д/ с функцией Гамильтона
H = - f|1 — 2pf I + fv+1/(v +1).
(8)
Решение уравнения (7) можно провести с помощью канонических преобразований и теории Га-мильтона-Якоби. Как следует из теоремы Якоби-Пуанкаре, если существует принадлежащая классу
С2 функция ?(/, р, г), такая, что д 2? / д/др Ф 0, то преобразование (/, р/) «(д, р), генерируемое этой функцией: р/ = д? / д/- д = д? / др, является
каноническим, а новая функция Гамильтона имеет вид
Н (д, р, г) = Н (/(д, р, г), р/ (д, р, г), г}+
±д? (/ (д, р, г), р, г}.
дт
В теории канонических преобразований для двумерного фазового пространства можно исходить из четырех типов производящих функций [3]. Для дальнейших целей наиболее подходящей является функция Е = Я2(р/, д, г). Тогда явный вид преобразований найдется из разрешения уравнений
/¦ дЕ дЕ
/=^-- р, (9)
др/ дд
а новая функция Гамильтона
Н (д, р, г) = Н (/(д, р, г), р/ (д, р, г), г}-
дЕ (() }. (10)
-^(р/д p, т), д, г}
дг
Перейдем от динамической системы {/, р/, Н } к представлению взаимодействия {д, р, Н }с помощью
производящей функции Е = 2дЛггИ4р/ - дг. Из уравнений (9) следует явный вид преобразований:
2 р + г.
f = 8qch 2
1 p + t
pf = - th--
f 4 2
(11)
(12)
а из (10) явный вид функции Гамильтона:
H = М^ ch 2(v+1) E±1. v +1 2
На заключительном этапе совершим второе преобразование от {q, p, H} к {& lt-2,P, H ° о}, откуда следует, что Q = const- P = const, а производящая функция преобразования F (p, Q, t) удовлетворяет уравнению H (dF / 5p, p, t) = dF / Эг. Теперь соотношения
q = dF / Bp, P = ЭF / Э2 (13)
в неявной форме задают первые интегралы гамиль-тоновой системы
q = ЭH /dp- p =-ЭH /Эq
(14)
Уравнение для производящей функции с учетом (12) имеет вид
(8ЭF / Эp) v+1 ch 2(v+1) Z±T = & lt-F
(15)
п +1 2 дг
Покажем, что частное решение системы (14) для случая нулевых значений произвольных постоянных Q = 0, Р = 0 приводит к солитонному решению исходного уравнения (1). В этом случае производящую функцию можно представить в виде
Е (р, Q/г) = я (р, г) + Qh (р, г), (16)
ограничившись первой степенью по Q, имея в виду, что после составления уравнений (13) следует положить Q = 0, Р = 0.
Подставляя (16) в (15) и удерживая только члены линейные по Q, находим
8v
v +1
f) + (v +1& gt- (f
Q Э±
Эp
2(v+1)
ch
p + г
=g + Q ^
Эг Эг
отсюда следуют два уравнения для g и h:
8V+1 (jg)V+1 ch 2(v+1) p +t Jg v + 1 (Эp) 2 Эг'-
8v+1|^gj dhch2(v+1) p + T =M
Э^ | Эp 2 Эг
Полагая g = g (p + г), из (17) находим
Эg (v +1)
1/v
Эp 81+1/v
ch
p + г
-2(1+1/v)
(17)
(18)
(19)
u =
r
2
2
Подстановка этого выражения в (18) приводит к линейному однородному уравнению первого порядка
дк дк
---(п +1) — = 0,
ат ар
простейшее нетривиальное решение которого имеет вид
к = p + (n + 1) t.
Соотношения (13) с учетом (16) дают:
q = *L + Q. dp dp
(20)
P = к.
Полагая здесь Q = P = 0, с учетом (19) и (20), находим:
p = -(n + 1) t,
)1/П
-, 1+1/n
q =
П1- (ck (nt/2))-2(1+1/n).
8
Подстановка последних двух формул в (11) дает решение исходной гамильтоновой системы с гамильтонианом (8):
f =
n +1
1/n
[ck (nt /2)]
-2/n.
(21)
pf = - - tk (nt/2)
Подстановка первой из этих формул в уравнение (7) обращает его в тождество.
Из приведенного выше масштабного преобразования и первой формулы (21) получается следующее выражение для интенсивности I:
I =
/ 2 1/n
a (n +1)
8h
[ck (asn /2)]
l-2/n
Учитывая (5) и вводя обозначение a -= 8^A2v / (n +1), окончательно находим решение (2):
exp & lt-i
y = A-
xu 2
V
A ln2^ n +1
t + & lt-
f
ck
1/n
nAl
2h n +1
(x — x0 — ut)
которое является гладкой функцией, локализованной вдоль направления x (t) = x0 +ut, и представляет собой волновой пакет, движущийся без дисперсионного уширения с постоянной скоростью u. Очевидно, что при v = 1 последняя формула превращается в известное решение НУШ с кубичной нелинейностью [1].
Заключение Таким образом, модуляция нелинейной системой (со степенным по интенсивности откликом на квазигармоническое возмущение) несущей волны в последовательность локализованных импульсов, имеющих форму гиперболического секанса, возведенного в некоторую степень, является общей чертой НУШ с любой положительной степенью нелинейности.
Литература
1. Тахтаджян, Л. А., Гамильтонов подход в теории соли-тонов / Л. А. Тахтаджян, Л. Д. Фаддеев // М.: «Наука», 1986. — 528с.
2. Степанов, В.В., Курс дифференциальных уравнений / В. В. Степанов // М.: Государственное издательство технико-теоретической литературы, 1953. — 468с.
3. Шмутцер, Э., Основные принципы классической механики и классической теории поля / Э. Шмутцер, // М.: Мир, 1976. — 160с.
References
1. Takhtajan L.A., Faddeev L.D. Hamilton approach in theory of solitons. — Moscow: Nauka, 1986, — 528p.
2. Stepanov V.V. Course of differential equations. — Moscow: GITTL, 1953. — 468p.
3. Schmutzer E. Grundprinzipien der klassischen Mechanik und der klassischen Feldtheorie (kanonischer Apparat). — Berlin: VEB Deutscher Verlag der Wissenschaften, 1973. — 160p.
4
8
AUTOMODULATION OF ONE-DIMENSIONAL WAVES BASED ON NONLINEAR SCHREDINGGER EQUATION WITH NON-KERR NONLIANERITY
Ivan V. Alimenkov (associatedprofessor, e-mail: i-alimenkov@mail. ru) S.P. Korolyov Samara State Aerospace University
Abstract
It is shown that nonlinear Schredinger equation with non-Kerr nonlinearity has a localized solution moving with constant velocity without dispersion. This solution is found by straight method based on Hamilton systems theory and it contains the well-known solution of cubic nonlinear Schredinger equation.
Key words: nonlinear Schredinger equation, nonlinearity 5Aorder, theory of Hamilton systems, canonical transformations, Hamilton-Jacoby equation, soliton solutions for degree nonlinearity.
В редакцию поступила 28. 05. 2009 г.

ПоказатьСвернуть
Заполнить форму текущей работой