Введение в теорию солитонов
.pdf
|
|
Уравнение ГЛМ |
|
51 |
Функции Φ(λ, x) и a(λ) аналитичны по λ в верхней полуплоскости |
||||
Im λ > 0, причем имеет место оценка |
|
|
||
|
Φ(λ, x) |
eiλy −e−iλ(x−y) = O λ−1e−iλ(x−y) |
, |
λ → ∞. |
|
|
|||
|
a(λ) |
|||
Деформируем контур интегрирования в левой части (9.9) с вещественной оси на параллельную ей прямую −∞ < Re λ < +∞, Im λ > 0. Поскольку |e−iλ(x−y)| = eImλ(x−y) → 0 при (x−y) < 0, Im λ → +∞, то в силу последней оценки экспонента пропадет и интеграл в левой части занулится.
Уравнение (9.9) принимает вид
|
+∞ |
|
|
|
|
0 =−∞ eiλy # |
Θ(λ, x)−e−iλx+r(λ)Θ(λ, x)$ dλ. |
(9.10) |
|||
Будем искать его решение в виде |
|
||||
|
|
|
|
∞ |
(9.11) |
|
Θ(λ, x) = e−iλx+ K(x, y)e−iλy dy, |
||||
|
|
|
|
x |
|
где K – новая неизвестная вещественная функция. Это представление справедливо в силу интегрального уравнения (7.12), выведенного выше при рассмотрении прямой задачи рассеяния:
|
∞ |
|
Θ(λ, x) = eiλx+ |
sin λ(x−y) |
u(y)Θ(λ, y)dy. |
|
λ |
|
|
x |
|
Подставляя эту функцию в уравнение (9.10) и заменяя переменную интегрирования y на ξ, получим
0 = +∞eiλy |
∞K(x, ξ)eiλξ dξ + r(λ) |
eiλx + |
∞K(x, ξ)eiλξ dξ dλ. |
|
−∞ |
x |
|
x |
|
(9.12) Первый двойной интеграл здесь – суть последовательно примененные прямое и обратное преобразования Фурье. Точнее, вводя функцию
Kˆ (x, y) = |
0, |
x > y, |
|
K(x, y), x < y, |
|
52 |
|
Лекция 9 |
|
|
имеем |
+∞ |
|
|
|
+∞ |
|
|
|
|
eiλy |
Kˆ (x, ξ)e−iλξ dξ dλ = 2πKˆ (x, y). |
|||
−∞ |
−∞ |
|
|
|
Уравнение (9.12) теперь перепишется при x < y так: |
|
|||
+∞ |
|
+∞ |
∞ |
|
0 = 2πK(x, y) + |
r(λ)eiλ(x+y)dλ + r(λ)eiλy |
K(x, ξ)eiλξ dξdλ. |
||
−∞ |
|
|
−∞ |
x |
Введем обозначение |
|
|
|
|
|
def |
1 |
+∞ |
|
|
|
|
||
|
F (z) = |
r(λ)eiλz dλ |
|
|
|
2π |
|
||
|
|
|
|
|
−∞
и в уравнении переставим порядок интегрирования в последнем двойном интеграле:
|
+∞ |
∞ |
|
eiλy r(λ) |
K(x, ξ)eiλξ dξdλ = |
|
−∞ |
x |
∞ |
+∞ |
∞ |
= |
K(x, ξ) r(λ)eiλ(ξ+y)dλdξ = K(x, ξ)2πF (ξ+y)dξ. |
|
x |
−∞ |
x |
Окончательно получается уравнение Гельфанда–Левитана–Марченко
∞
K(x, y)+F (x+y)+ K(x, ξ)F (ξ+y)dξ = 0, x < y. |
(9.13) |
x
2. Формула обращения
Нам осталось сделать последний шаг в схеме метода обратной задачи рассеяния – перейти от решения K(x, y) уравнения ГЛМ к искомому решению уравнения КдФ u(x, t). Этот переход дается модификацией формулы обращения (9.6), выведенной в начале этой лекции.
Перепишем выражение (9.6) для потенциала u в терминах функции
Θ = χe−iλx: |
|
|
|
|
|
|
|
|
u(x) = |
∂ |
lim 2iλ(χ(λ, x) |
− |
1) = |
∂ |
lim 2iλ(Θ(λ, x)e−iλx |
− |
1). |
|
|
|||||||
|
∂x λ |
|
∂x λ |
|
||||
|
|
→∞ |
|
|
|
→∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
Формула обращения |
|
|
|
|
|
|
|
|
|
53 |
||||||||||
Подставим сюда выражение (9.11) для Θ(λ, x) |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Θ(λ, x) = eiλx+ |
K(x, ξ)eiλξ dξ, |
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
что дает |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
u(x) = ∂x λlim |
2iλ |
∞ |
K(x, ξ)e |
iλ(ξ |
− |
x) |
dξ |
= |
|
|
|
|||||||||||||||||
|
|
|
|
|
∂ |
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
→∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
∂x |
λ→∞ |
∞ |
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
2 x |
K(x, ξ)d |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
= |
∂ |
|
|
lim |
|
|
|
|
|
eiλ(ξ |
|
|
x) |
|
= |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
= 2 ∂x λlim |
|
|
K(x, x) |
∞ |
|
|
|
|
− Kξ (x, ξ)dξ |
= |
|||||||||||||||||
|
|
|
|
|
|
e |
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
∂ |
|
→∞ − |
|
|
− x |
|
iλ(ξ |
x) |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
= |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
iλ |
Kξ (x, x) + |
∞ |
|
|
|
− |
|
|
|
||||||||
− |
∂xK(x, x) + 2 ∂x λlim |
|
|
e |
|
|
|
Kξξ (x, ξ)dξ . |
|||||||||||||||||||||||
|
|
|
∂ |
|
|
∂ |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
iλ(ξ |
|
x) |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь последнее слагаемое обращается в нуль, поскольку интеграл в силу гладкости K(x, ξ) имеет порядок O(λ−N ) при λ → ∞, N – любое положительное число (N раз проинтегрировать по частям).
Таким образом, окончательно получим следующую формулу, которая называется формулой обращения:
u(x) = −2 |
∂ |
|
∂x K(x, x). |
(9.14) |
Подводя итог всему сказанному, сформулируем стандартный алгоритм решения обратной задачи рассеяния для уравнения КдФ.
•По данным рассеяния a(λ) = a(λ, 0) и b(λ, t) = b(λ, 0)e−8iλ3t находим коэффициент отражения
r(λ, t) = b(λ, t) = b(λ, 0) e−8iλ3t. a(λ, t) a(λ, 0)
54 |
Лекция 9 |
|||
• Вычисляем ядро уравнения ГЛМ |
|
|||
F (z, t) = |
|
1 |
|
+∞ |
|
|
r(λ, t)eiλz dλ. |
||
|
2π |
|||
|
|
|
||
|
|
|
−∞ |
|
• Решаем уравнение ГЛМ (9.13)
∞
K(x, y, t)+F (x+y, t)+ K(x, ξ, t)F (ξ+y, t)dξ = 0, x < y.
x
• Находим решение уравнения КдФ по формуле обращения (9.14)
u(x, t) = −2 ∂x∂ K(x, x, t).
Лекция 10
1. Дискретный спектр уравнения Шредингера
При выводе уравнения Гельфанда–Левитана–Марченко (9.13)
∞
K(x, y) + F (x + y) + K(x, ξ)F (ξ + y)dξ = 0
x
было существенно использовано условие a(λ) = 0, Imλ > 0. Это условие необходимо, чтобы сдвинуть контур интегрирования по λ для интеграла от левой части равенства
|
Φ(x, λ) |
|
|
|
|
|
|
||
|
|
= Ψ(x, λ) + r(λ)Ψ(x, λ). |
(10.1) |
||||||
|
|
||||||||
|
a(λ) |
|
|
||||||
В силу аналитичности функции |
Φ(x, λ) |
|
можно сдвинуть контур интегри- |
||||||
a(λ) |
|||||||||
|
|
|
|
|
|
|
|||
рования на Im(λ) вверх, а затем Im(λ) устремить в бесконечность. Если a(λ) имеет конечное число нулей в точках iλn, λn > 0, n = 1, 2, . . . , N , то справедливо разложение Тейлора
a(λ) = (λ − iλn)a (iλn) + O(λ − iλn)2, λ → iλn.
Предполагая a (iλn) = 0, имеем по теореме о вычетах |
|
|
|||||||||
|
a(λ) |
n=1 |
λ=iλn |
a(λ) |
|
=2πi n=1 |
a (iλn) |
||||
∞Φ(x, λ)eiλy |
N |
|
Φ(x, λ)e−iλy |
N |
Φ(x, iλn)e−λny |
||||||
−∞ |
|
|
|
|
|
|
|
|
|
||
|
|
dλ = 2πi |
|
Res |
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
(10.2) |
|
В дальнейшем мы воспользуемся этой формулой для обобщения уравнения Гельфанда–Левитана–Марченко.
Введем важное определение дискретного спектра уравнения Шредин-
гера |
|
Φxx + (λ2 + u(x))Φ = 0. |
(10.3) |
56 |
Лекция 10 |
Оказывается, что в классе быстроубывающих потенциалов
∞
(1 + |x|)|u(x)|dx < ∞
−∞
уравнение (10.3), кроме знакомых нам осциллирующих решений
|
|
|
eiλx, |
x → −∞, |
|||
Φ |
→ |
||||||
|
|
(λ)eiλx + |
b |
(λ)e−iλx, |
x → +∞, |
||
|
a |
||||||
имеет экспоненциально убывающие по x решения Φ(x, λ) при некоторых (комплексных) значениях λ. Такие решения соответствуют разрешенным состояниям квантовомеханической системы, описываемой уравнением (10.3), которое, как уже отмечалось выше, является одномерным стационарным уравнением Шредингера.
Точное определение дискретного спектра состоит в следующем. ОПРЕДЕЛЕНИЕ. Значение λ C принадлежит дискретному спектру
уравнения (10.3) с потенциалом u(x), если Φ(x, λ) L2(R), то есть
∞
|Φ(x, λ)|2 dx < const.
−∞
Эти значения λ называются собственными значениями, а соответствующие решения – собственными функциями дифференциального оператора Шредингера d2/dx2 + u(x).
1.1. Свойства дискретного спектра
1.Точки дискретного спектра совпадают с нулями коэффициента про-
хождения, то есть λn - собственное значение тогда и только тогда, когда a(λn) = 0.
Для доказательства воспользуемся выражением a(λ) через вронскиан (см. Свойство 3 данных рассеяния)
W (Φ, Ψ) a(λ) = ,
2iλ
откуда при a(λ) = 0, Im λ > 0 получим W (Φ(x, λn ), Ψ(x, λn)) = 0. Это означает, что Ψ и Φ линейно зависимы, то есть существует константа bn такая, что
Φ(x, λn) = bnΨ(x, λn).
|
|
|
Дискретный спектр |
|
|
|
57 |
||
Оценивая эти функции при больших значениях x, имеем |
|
|
|
|
|||||
|
|Φ(x, λn)| e−iλnx |
= eImλnx,Im |
Imλn > 0, x → −∞, |
+ |
|
, |
|||
| |
Φ(x, λn) = bn |
|| |
Ψ(x, λn) |
bn e− λnx, Imλn > 0, x |
→ |
∞ |
|||
| | |
| | | |
|
|
|
|||||
откуда следует, по определению, что Ψ и Φ есть собственные функции,
а |
λn – собственное значение. Обратно, из условия убывания |
Ψ(x, λn) |
на |
бесконечности по x и постоянства вронскиана, следует |
равенство |
W (Φ(x, λn ), Ψ(x, λn)) = 0, что эквивалентно a(λ) = 0.
2. Теорема Марченко. Если u(x) – гладкий потенциал и
+∞
(1+ | x |) | u(x) | dx < const,
−∞
тогда множество собственных значений {λ√n} не более чем конечно.
В частности, если интеграл меньше ln(2 + 3), то дискретного спектра нет, то есть a(λ) нигде не обращается в нуль.
Доказательство этой теоремы достаточно сложно и выходит за рамки данного курса.
3.При вещественном u(x) все собственные значения лежат на мнимой оси:
λn = iκ n, κ n > 0, n = 1, 2, . . . , N.
Для доказательства воспользуемся тем, что если Φ(x, λn) – собственная
функция, то Φ(x, −λn) – тоже собственная функция,
Φxx + ((−λn)2 + u(x))Φ = 0,
поскольку u(x) = u(x). Граничные значения Φ(x, λn ) и Φ(x, −λn) также
совпадают:
Φ(x, λn) → eiλnx, x → −∞,
|
|
|
|
n) |
|
|
|
nx = eiλnx, x |
|
Φ(x, |
− |
λ |
→ |
e−iλ |
. |
||||
|
|
|
|
|
|
|
→ −∞ |
||
Таким образом, эти функции совпадают и, следовательно, совпадают их
собственные значения λn = −λn, а это может быть, если они лежат на мнимой оси.
58 |
Лекция 10 |
4.Для различных собственных значений соответствующие собственные функции ортогональны:
+∞
Φn, Φm = Φ(x, λn)Φ(x, λm)dx = 0, λn = λm.
−∞
Для доказательства определим Φn = Φ(x, λn), Φm = Φ(x, λm), λn = = iκ n, λm = iκ m. Функции Φn и Φm удовлетворяют уравнению (10.3)
(Φn)xx + (−κ n2 |
+ u(x))Φn = 0, |
(10.4) |
||
(Φm)xx + (−κ m2 |
+ u(x))Φm = 0. |
|
||
Умножим (10.4) на Φm и |
проинтегрируем |
обе части на |
интервале |
|
(−∞, +∞). |
|
|
|
|
Из свойства 3 следует, что собственные функции на бесконечности |
||||
экспоненциально убывают |
|
|
|
|
Φn → e−κ n|x|, |
Φm → e−κ m|x|, |
x → ±∞, |
|
|
поэтому все интегралы ниже корректно определены. Интегрируя по частям, имеем
+∞
0 = Φm(x) #(Φn)xx + (−κ 2n + u(x))Φn(x)$ dx =
−∞
+∞
=#Φm(x)d (Φn)x + (−κ 2n + u(x))ΦmΦndx$ =
−∞ |
|
|
|
+∞ |
|
|
|
+∞ |
|
2 |
|
= Φm(Φn)x −∞ + |
# |
−(Φn)x(Φm)x + (−κ n + u(x))ΦmΦn |
dx = |
−∞ |
|
$ |
|
|
|
|
|
|
|
|
|
+∞
=#(−Φm)xd(Φn) + (−κ 2n + u(x))ΦmΦndx$ =
−∞ |
|
|
+∞ |
+∞ |
$dx = |
2 |
||
= (Φm)xΦn −∞ + |
Φn #(Φm)xx + (−κ n + u(x))Φm |
|
|
|
|
|
|
|
−∞
Дискретный спектр |
59 |
+∞
=Φn(x)[(Φm)xx + (−κ 2n + u(x))Φm + (−κ 2n + κ 2m)Φm]dx =
−∞
= (−κ n2 + κ m2 ) |
+∞ |
ΦnΦmdx. |
|
|
−∞ |
Поскольку по условию κ n = κ m, последний интеграл равен нулю, что и требовалось.
5. Скалярный квадрат собственной функции равен
Φn(x, iκ n) 2 = Φn, Φn = ibna (iκ n),
где |
|
Φn → eκ nx, |
|
x → −∞, |
|
|
|
||
|
Φn → bne−κ nx, x → +∞. |
|||
Для доказательства продифференцируем уравнение (10.3) по λ: |
||||
|
∂ |
2 |
|
|
0 = |
|
(Φn)xx + (λ |
|
+ u(x))Φn = |
∂λ |
|
|||
|
= (Φnλ)xx + (λ2 + u(x))Φnλ + 2λΦn. |
||
Умножим последнее уравнение на Φn и проинтегрируем по x на ин- |
|||
тервале (−∞, +∞). Интегрируя по частям, имеем |
|
||
+∞ |
+∞ |
|
|
−2λ−∞ Φn2 dx =−∞ Φn #(Φnλ)xx + (λ2 + u(x))ΦnΦnλ)$dx = |
|||
|
+∞ |
|
|
|
+∞ |
|
|
= Φn(Φnλ)x −∞ + |
−(Φn)x(Φnλ)x + (λ2 + u(x))ΦnΦnλ) dx = |
||
|
−∞ |
# |
$ |
|
|
|
|
|
|
|
|
|
+∞ |
|
|
= bnκ na (iκ n) + ∞ |
#(Φn)xd(Φnλ) + (λ2 + u(x))ΦnΦnλ)$dx = |
||
|
= bnκ na (iκ n) − (Φn)x(Φnλ) +∞ |
+ |
|
|
|
|
|
−∞
60 |
Лекция 10 |
+∞
+#(Φn)xxΦnλ + (λ2 + u(x))ΦnΦnλ)$dx =
−∞
+∞
=Φnλ #(Φn)xx + (λ2 + u(x))Φn$dx + 2κ nbna (iκ n) =
−∞
= 2κ nbna (iκ n).
