Введение в теорию солитонов
.pdf
|
|
2. Свойства данных рассеяния |
41 |
|
Ψ11 |
→ |
e4itλ3 2λbeiλx |
, |
x → +∞. |
Ψ21 |
e4itλ3 (ae−iλx + beiλx |
|||
Уравнение Ψt = AΨ для первого столбца при x → +∞ переписывается в виде
|
e4itλ3 (2λbt + 4iλ3 2λb)eiλx |
= |
||
e4itλ3 (at + 4iλ3a)e−iλx +·(bt + 4iλ3b)eiλx |
||||
= −4λ2 |
iλ |
0 |
e4itλ3 2λbeiλx |
|
i |
−iλ e4itλ3 (ae−·iλx + beiλx) , |
|||
откуда получаем |
at + 4iλ a |
·= 4iλ a. |
|
|
|
|
iλ3 |
3 |
|
|
2bt + 4 3 |
2b = 3−4iλ · 2b, |
|
|
Таким образом, эволюция от t данных рассеяния очень проста и составляет содержание первого из свойств данных рассеяния, формулируемых ниже.
2. Свойства данных рассеяния
a(λ, t) = a(λ, 0),
1.b(λ, t) = b(λ, 0)e−8itλ3 .
Следующие свойства описывают зависимость данных рассеяния от λ. Будем опускать аргумент t, если это не вызывает недоразумений.
2. |a(λ)|2 − |b(λ)|2 = 1, λ R.
Доказательство следует из закона сохранения вронскиана W = const, где W (Φ1, Φ2) = Φ1Φ2 − Φ1Φ2, а Φ1 и Φ2 - два произвольных решения уравнения Шредингера.
Φxx+(λ2+u(x))Φ = 0. |
(7.6) |
В самом деле, из вещественности λ2 + u следует, что Φ также является решением уравнения (7.6), а в силу нормировки (6.5)
|
|
|
|
|
|
eiλx, |
|
|
|
|
|
x → −∞, |
|
|
|
|
Φ |
→ |
|
|
|
|
|
||||
|
|
|
|
|
(λ)eiλx + |
b |
(λ)e−iλx, |
x → +∞. |
|||||
|
|
|
|
a |
|||||||||
Тогда |
|
|
|
|
|
|
|
|
|
|
|||
W (Φ, |
|
|
= 2iλ = W (Φ, |
|
|
|
+ = 2iλ(|a|2 − |b|2), |
||||||
Φ) x→−∞ |
Φ) x |
|
|||||||||||
|
|
|
|
|
|
|
|
||||||
что и доказывает свойство 2. |
|
→ ∞ |
|||||||||||
42 |
Лекция 7 |
3. a(λ) = −21iλ W (Φ, Θ), где Φ и Θ – два решения уравнения (7.6), заданные условиями
Φ → e−iλx, |
x → −∞, |
Θ → eiλx, |
x → +∞. |
Доказательство следует из вычислений вронскиана решений Φ и Θ
W (Φ, Θ)|x→+∞ = (−iλae−iλx+iλbeiλx)eiλx−
−(ae−iλx+beiλx)iλeiλx = W (Φ, Θ)|x→−∞ = −2iλa.
4.Функция a(λ) аналитична при Im λ > 0, причем a(λ) → 1 при λ → ∞. Аналогично, a(λ) – аналитическая функция при Im λ < 0, и a(λ) → 1,
λ → ∞.
Доказательство: Покажем, что Φ и Θ – аналитические в полуплоскости Im (λ) > 0. Функция Φ удовлетворяет следующему интегральному уравнению:
|
|
x |
|
|
Φ(x, λ) = e−iλx |
− |
sin λ(x−ξ) |
u(ξ)Φ(ξ, λ) dξ. |
(7.7) |
|
λ |
|
||
−∞
Оно получается из нормировки (7.4) и функции Грина G для уравнения Шредингера (7.6), где G(λ, x, ξ) = sin λ(x −ξ)/λ при ξ < x и G(λ, x, ξ) = 0 при ξ > x. Впрочем, уравнение (7.7) можно проверить и непосредственно, что удобно сделать для функции
χ(x, λ) = Φ(x, λ)eiλx. |
(7.8) |
Подставляя (7.8) в уравнение (7.7), получим эквивалентное интегральное уравнение
|
|
x |
|
|
χ(x, λ) = 1 |
− |
e2iλ(x−ξ)−1 |
u(ξ)χ(ξ, λ) dξ. |
(7.9) |
|
2iλ |
|
||
−∞
С другой стороны, подставляя (7.8) в уравнение Шредингера (7.6), получим
Φxx = (χxx−2iλχx−λ2χ)e−iλx = u(x)χe−iλx, |
|
то есть |
|
χxx−2iλχx+uχ = 0. |
(7.10) |
Свойства данных рассеяния |
43 |
Продифференцируем теперь обе части интегрального уравнения (7.9)
x
χx = − e2iλ(x−ξ)u(ξ)χ dξ, χxx = −u(x)χ(λ, x) −
−∞ x
−2iλ e2iλ(x−ξ)u(ξ)χ(λ, ξ)dξ = −u(x)χ+2iλχx,
−∞
так что выполнено дифференциальное уравнение (7.10) для функции χ(λ, x) и формула (7.7) доказана.
Выберем Im λ > 0. Тогда в интегральном уравнении (7.11) ядро аналитично и ограничено по λ, поскольку
| e2iλ(x−ξ) |= e−2Im λ(x−ξ) < 1, (x−ξ) > 0, Imλ > 0. |
(7.11) |
Тогда его решение χ(λ, x) – суть аналитическая и ограниченная во всей верхней полуплоскости функция по λ. Отметим, что существование и единственность этого решения следует из того факта, что уравнение (7.11) является уравнением Вольтерра, причем абсолютная сходимость интеграла на бесконечности обеспечена оценками (7.3), (7.11) и условием χ(x, λ) = = 1 + o(1) при x → −∞.
Аналогично, докажем аналитичность в верхней полуплоскости второй функции Θ(λ, x), входящей во вронскиан в свойстве 3:
Θ(λ, x) → eiλx, λ → +∞.
Повторяя вышеприведенные рассуждения, легко получить интегральное уравнение
|
+∞ |
|
|
Θ(λ, x) = eiλx+ |
sin λ(x−ξ) |
u(ξ)Θ(λ, ξ)dξ, |
(7.12) |
|
λ |
|
|
x
которое анализируется аналогично уравнению (7.11). Детали этого доказательства оставляем читателю в качестве упражнения.
Поскольку функции Φ и Θ аналитичны по λ в верхней полуплоскости, то таковым же будет и их вронскиан W , а значит, и функция a(λ) = = −W (Φ, Θ)/2iλ (свойство 3).
Лекция 8
1. Свойства данных рассеяния (продолжение)
Проверим ограниченность a(λ) при λ → ∞. Из интегральных урав-
нений для Φ и Θ имеем: |
|
|
|
|
Θ(λ, x) = e−iλx + o(1), |
Θx = iλeiλx+O(1), |
|||
Φ(λ, x) = e−iλx + o(1), |
Φx = −iλe−iλx+O(1), |
|||
так что |
1 |
|
|
|
a(λ) = − |
W (Φ, Θ) → 1, λ → ∞. |
|||
2iλ |
||||
Свойства функции a(λ) проверяются аналогичным образом, выполняя комплексное сопряжение в интегральных уравнениях для Θ и Φ.
Свойство 4 доказано.
ПРИМЕР 1: ДЕЛЬТА-ОБРАЗНЫЙ ПОТЕНЦИАЛ.
Рассмотрим пример вычисления данных рассеяния для одного специального потенциала – дельта-функции Дирака u(x) = δ(x). Это обобщенная функция, которая может быть определена так:
|
|
0, x < −2ε , |
|
|
δ(x) = lim |
1ε , ε < x < ε , |
(8.1) |
||
→ |
|
0, x > ε . |
|
|
ε 0 |
|
−2 |
2 |
|
|
|
2 |
|
|
|
|
|
|
|
Можно представлять себе дельта-функцию как «бесконечно высокую стенку бесконечно малой толщины». В физике с ее помощью удобно записывать пространственную плотность различных величин (массы, заряда, интенсивность источника тепла, силу и т.п.), сосредоточенных в точке x0. Так, плотность точечной массы m, находящейся в точке a, записывается в виде mδ(x − a). В теории обобщенных функций доказывается, что для любой непрерывной функции f (x)
∞
f (x)δ(x − a)dx = f (a). |
(8.2) |
−∞
Свойства данных рассеяния |
45 |
Это равенство можно принять и за определение дельта-функции, оно эквивалентно формуле (8.1).
Воспользуемся свойством (8.2) для вычисления данных рассеяния. Подставляя u(x) = δ(x) в интегральное уравнение (7.7), получим
0, |
x < 0, |
Φ(λ, x) = e−iλx− sin λx |
Φ(λ, 0), x > 0. |
λ |
|
Переходя здесь к пределу при x → 0, немедленно получаем Φ(λ, 0) = 1, откуда имеем окончательно
|
|
|
|
|
|
|
0, |
x < 0, |
|
|
|
|
|
|
|
|
|||
Φ(λ, x) = e−iλx − sin λx , x > 0. |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
λ |
|
|
|
|
|
|
|
|
|
|
|
|
Отсюда при любом x > 0 получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Φ(λ, x) = e−iλx |
− |
eiλx−e−iλx |
= (1+ |
1 |
|
)e−iλx |
− |
1 |
eiλx, |
||||||||||
|
|
|
|||||||||||||||||
|
|
|
2iλ |
|
2iλ |
|
2iλ |
|
|
||||||||||
так что по определению (Φ |
→ |
a(λ)e−iλx + b(λ)eiλx, |
x |
→ |
+ |
∞ |
) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
a(λ) = 1+ |
|
, b(λ) = − |
|
. |
|
|
|
|
|
|
|
|
|||||||
2iλ |
2iλ |
|
|
|
|
|
|
|
|
||||||||||
В заключение данного примера отметим, что данные рассеяния можно вычислить в явном виде для любого «ступенчатого» потенциала вида
|
|
|
u(x) = |
|
0, |
x < x0, |
|
|
(8.3) |
||||||
|
|
|
u |
, x |
|
< x < x , |
|
||||||||
|
|
|
|
|
|
.u.1., x0 |
< x < x1 |
, |
|
|
|
||||
|
|
|
|
|
. .k. k−1 |
|
|
k |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0, |
xn < x. |
|
|
|
|
||||
В самом деле, на каждом |
интервале |
x |
|
|
< x < x |
мы имеем для функ- |
|||||||||
|
|
|
|
k−1 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
k |
2 |
|
|
||
ции Φ уравнение (7.6) с постоянными коэффициентами λ + uk, так что ре- |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
√ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
шение находится в виде линейной комбинации экспонент exp(i |
λ + ukx) |
||||||||||||||
√ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и exp(−i λ |
+ uk x). Коэффициенты при этих экспонентах подбираются из |
||||||||||||||
условий согласования Φ и ее производной на границе соседних интервалов:
Φ(λ, xk − 0) = Φ(λ, xk + 0), Φx(λ, xk − 0) = Φx(λ, xk + 0).
В итоге полученные формулы для a и b будут представлять собой композицию экспонент с рациональными по λ коэффициентами. К сожалению, они достаточно громоздки и мы их здесь не приводим.
46 |
Лекция 8 |
ПРИМЕР 2: ДАННЫЕ РАССЕЯНИЯ СОЛИТОННОГО ПОТЕНЦИАЛА. Вычислим величины a(λ) и b(λ) для односолитонного потенциала
2a2 u(x) = .
ch2(ax)
Решение уравнения Штурма-Лиувилля
2 |
2a2 |
||
Φxx + λ + |
|
|
Φ = 0 |
ch2 |
|
||
|
(ax) |
||
будем искать в виде |
|
|
|
Φ(x, λ) = |
C − th(ax) |
e−iλx, C = const. |
(8.4) |
|
|||
|
C + 1 |
|
|
Выбор такой функции диктуется стремлением получить члены вида композиции экспонент eax и e−ax в уравнении, а также необходимостью добиться нужной нормировки на бесконечности
|
e−iλx, |
x |
→ −∞ |
, |
|
||||||
Φ → |
C − 1 |
|
|
|
|
|
|||||
e−iλx, |
|
x |
→ |
+ |
∞ |
. |
|||||
|
|
C + 1 |
|
|
|
|
|
|
|
||
Дифференцируя функцию (8.4), имеем |
|
|
!e− |
|
|
|
|||||
Φxx = |
ch3(ax) |
− ch2(ax) |
|
|
− λ Φ. |
||||||
|
2a2sh(ax) |
|
2iλa |
|
|
iλx |
2 |
||||
iλ
Подстановка в уравнение немедленно дает C = a , откуда по определению данных рассеяния имеем
a(λ) = |
C − 1 |
= |
iλ − a |
, b(λ) |
≡ |
0. |
|
C + 1 |
|
iλ + a |
|
|
2.Коэффициент отражения и интегральная формула Коши
Обсудим вопрос о величине произвола в данных рассеяния. Они состоят из двух комплекснозначных функций a и b, то есть параметризуются четырьмя вещественными функциями на прямой λ R. С другой стороны,
Коэффициент отражения |
47 |
данные рассеяния полностью определяются одной вещественной функцией – потенциалом u(x). Этот факт означает, что должны существовать соотношения на функции a(λ) и b(λ), оставляющие только один вещественный функциональный произвол.
Одно из этих соотношений, |a|2 − |b|2 = 1, мы уже вывели в Лекции 7 (см. свойство 2). Сейчас будут указаны остальные функциональные зависимости между данными рассеяния.
Покажем сначала, что данные рассеяния a и b однозначно выражаются через одну комплекснозначную функцию
b(λ)
r(λ) = |
|
, |
(8.5) |
|
a(λ)
называемую коэффициентом отражения.
Для модулей a, b и r справедливы соотношения:
2 |
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
1− |
| b | |
= |
|
|
, |
1− | r |2= |
|
|
, |
|
||||
2 |
|
2 |
| a | |
2 |
|
|||||||||
|
| a | |
|
|
| a | |
|
|
|
|
|
|||||
| a |2= |
|
1 |
|
, |
| b |2= |
| r |2 |
. |
(8.6) |
||||||
|
1− | r | |
2 |
2 |
|||||||||||
|
|
|
|
|
|
|
1− | r | |
|
|
|
||||
Для восстановления самой функции a(λ) воспользуемся ее аналитичностью в верхней полуплоскости. Как известно из теории функций комплексного переменного, аналитичность влечет связь между вещественной и мнимой частями функции (условия Коши–Римана).
Напомним, что для любой аналитической в области D функции g(λ) имеет место формула Коши
|
|
|
g(λ) = 2πi "Γ |
|
ξ−λ dξ, |
|
|
(8.7) |
|||||||||
|
|
|
|
|
1 |
|
|
|
g(ξ) |
|
|
|
|
|
|
|
|
где контур интегрирования Γ – граница области D – проходится против ча- |
|||||||||||||||||
совой стрелки. Далее, для граничного значения g(λ+0) = |
lim |
g(ξ) |
|||||||||||||||
справедлива формула Сохоцкого–Племеля |
|
|
|
|
|
ξ→λ, ξ D, λ Γ |
|
||||||||||
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
g(λ) |
1 |
|
|
g(ξ) |
|
|
|
|||||
g(λ + 0) = |
|
+ |
|
|
|
|
|
dξ, |
|
(8.8) |
|||||||
2 |
2πi |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
Γ |
|
ξ−λ |
|
|
|
||
где интеграл понимается в смысле главного значения |
|
|
|||||||||||||||
|
g(ξ) |
dξ = lim |
|
|
|
|
|
g(ξ) |
dξ. |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
ξ |
− |
λ |
|
|
|
|
|
− |
|
|
||||||
|
|
ε→0 |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
ξ λ |
|
|
|
||||||
Γ |
Γ/{|λ−ξ|>ε} |
|
48 |
Лекция 8 |
Выберем теперь в качестве D верхнюю полуплоскость, а в качестве g(λ) – логарифм функции a(λ) (мы предполагаем, что a(λ) = 0 всюду при λ D)
g(λ) = ln a(λ), D = {λ |Imλ > 0}, Γ = R.
Поскольку в силу свойства 4 a(λ) → 1, имеем ln a → 0 при λ → ∞ и интеграл сходится
+∞
ln a(λ+0) = |
ln a(λ) |
− |
1 |
2 |
2πi |
−∞
ln a(ξ)
ξ−λ dξ.
В силу свойства логарифма комплекснозначной функции
ln a(λ+0) = ln a(λ) = ln |a(λ)| + i arg a(λ),
так что
|
|
+∞ |
|
|
|
|
ln a(λ) + i arg a(λ) = |
1 |
|
ln |a(ξ)|+i arg a(ξ) |
dξ, |
||
|
|
|||||
| | |
πi |
|
ξ |
− |
λ |
|
|
|
|||||
−∞
или, окончательно, приравнивая чисто мнимые части,
arg a(λ) = 1 |
∞ |
ln |a(ξ)| |
dξ = |
1 |
∞ |
ln(1 − |r(ξ)|2 ) |
dξ. |
||||
|
|
|
|||||||||
−π |
|
ξ |
− |
λ |
2π |
|
ξ |
− |
λ |
||
|
|
|
|
||||||||
|
−∞ |
|
|
|
|
−∞ |
|
|
|
||
Таким образом, функция a(λ) однозначно выражается через |r(λ)|. Вспоминая формулы (8.6), (8.5), убеждаемся, что b(λ) также полностью выражается через r(λ):
| b |2= |
| r |2 |
|
, arg b(λ) = arg r(λ) + arg a(λ). |
1− | r | |
2 |
||
|
|
|
Используем далее вещественную редукцию, налагаемую на решение уравнения Шредингера условием вещественности потенциала u(x):
Φ(x, −λ) = Φ(x, λ).
В самом деле, взятие комплексного сопряжения и замена λ → −λ не изменяют уравнение (7.6) и граничное условие на Φ. Отсюда непосредственно вытекает такая же редукция на данные рассеяния:
a(−λ) = a(λ), b(−λ) = b(λ), r(−λ) = r(λ).
Последнее соотношение означает, что r(λ) параметризуется, например, своей вещественной частью при λ < 0 и своей мнимой частью при λ > 0.
Лекция 9
1.Уравнение Гельфанда–Левитана–Марченко
Вэтой лекции будет выведено уравнение обратной задачи рассеяния в случае отсутствия дискретного спектра, то есть всюду в этой лекции будем предполагать, что
a(λ) = 0 для всех Imλ 0. |
(9.1) |
Наша задача – восстановить потенциал u(x) в уравнении Шредингера (7.6)
|
|
Φxx+(λ2 |
+u(x))Φ = 0 |
(9.2) |
по его данным рассеяния a(λ) и b(λ): |
|
|||
|
→ |
a(λ)e−iλx |
+b(λ)eiλx, x → +∞. |
|
Φ |
|
e−iλx, |
x → −∞, |
(9.3) |
Нам понадобится также каноническое решение Θ, определенное выше уравнением Шредингера (7.6)
Θxx+(λ2+u(x))Θ = 0, |
(9.4) |
Θ → eiλx, x → +∞. |
(9.5) |
Выведем сначала формулу обращения для потенциала u(x) по известной функции Φ.
Поскольку интегральное уравнение (7.7) для Φ эквивалентно уравнению (7.9) для χ(λ, x) = Φ(λ, x)eiλx,
|
+∞ |
χ(λ, x) = 1 − |
e2iλ(x−ξ) |
2iλ u(ξ)χ(λ, ξ) dξ. |
|
|
−∞ |
Используя аналитичность χ по λ при Im λ > 0, перейдем к пределу при λ → ∞. Воспользовавшись оценкой
| e2iλ(x−ξ) |= e−2Im λ(x−ξ) → 0, λ → ∞, x > ξ,
50 |
|
|
Лекция 9 |
|
|
|
||
получаем |
|
|
|
|
|
|
|
|
|
∂ |
2iλ(χ(λ, x)−1) = |
∂ |
2iλ(Φ(λ, x)e |
iλx |
−1). (9.6) |
||
u(x) = |
|
λlim |
|
λlim |
|
|||
∂x |
∂x |
|
||||||
|
|
→∞ |
|
|
→∞ |
|
|
|
Приступая непосредственно к выводу уравнения обратной задачи, напишем соотношение рассеяния (9.3) несколько по-другому. Считая λ веще-
ственным, введем третье решение Θ(x, λ) уравнения Шредингера (7.6):
|
|
|
= 0, |
|
Θ |
xx+(λ2+u(x))Θ |
(9.7) |
||
|
→ e−iλx, x → +∞. |
(9.8) |
||
Θ |
||||
Тогда между тройкой функций Φ, Θ и Θ существует линейное соотношение, поскольку они являются решениями одного и того же уравнения (7.6) второго порядка. Коэффициенты этого соотношения достаточно вычислить в какой-нибудь одной точке x, так как они от x не зависят. В силу асимптотик (9.3), (9.4) и (9.7) имеем
Φ(λ, x) = a(λ)Θ(λ, x)+b(λ)Θ(λ, x) для любых x, λ R.
Разделим на a(λ) обе части этого соотношения и умножим на экспоненту exp(iλy), где y – вещественный параметр.
a(λ) |
eiλy = |
Θ(λ, x)+ a(λ) Θ(λ, x)! eiλy = |
|
Θ(λ, x)+r(λ)Θ(λ, x) eiλy . |
|||
Φ(λ, x) |
|
|
|
b(λ) |
# |
|
$ |
Проинтегрируем это соотношение по λ, вычтя из обеих частей экспоненту exp(−iλ(x − y)):
+∞ |
eiλy −e−iλ(x−y)!dλ = |
|
|||
|
a(λ) |
|
|||
|
Φ(λ, x) |
|
|
||
−∞ |
|
|
|||
+∞ |
|
|
|
|
|
=−∞ eiλy # |
Θ(λ, x)−e−iλx+r(λ)Θ(λ, x)$ dλ. |
(9.9) |
|||
Интегралы в последней формуле определены корректно, поскольку их сходимость на бесконечности следует из оценки подынтегральных функций O λ−1 exp(iλ(x − y)) , которая, в свою очередь, обеспечена асимптотиками (9.3), (9.4), (9.7) и свойством 4.
