Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Введение в теорию солитонов

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
874 Кб
Скачать

 

 

2. Свойства данных рассеяния

41

Ψ11

e4itλ3 2λbeiλx

,

x → +∞.

Ψ21

e4itλ3 (ae−iλx + beiλx

Уравнение Ψt = AΨ для первого столбца при x → +переписывается в виде

 

e4itλ3 (2λbt + 43 2λb)eiλx

=

e4itλ3 (at + 43a)e−iλx +·(bt + 43b)eiλx

= 4λ2

0

e4itλ3 2λbeiλx

 

i

−iλ e4itλ3 (ae·iλx + beiλx) ,

откуда получаем

at + 4iλ a

·= 4iλ a.

 

 

 

3

3

 

 

2bt + 4 3

2b = 34iλ · 2b,

 

Таким образом, эволюция от t данных рассеяния очень проста и составляет содержание первого из свойств данных рассеяния, формулируемых ниже.

2. Свойства данных рассеяния

a(λ, t) = a(λ, 0),

1.b(λ, t) = b(λ, 0)e8itλ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 ,

 

 

= 2= W ,

 

 

 

+ = 2(|a|2 − |b|2),

Φ) x→−∞

Φ) x

 

 

 

 

 

 

 

 

 

что и доказывает свойство 2.

 

→ ∞

42

Лекция 7

3. a(λ) = 21W , Θ), где Φ и Θ – два решения уравнения (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)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

e2(x−ξ)1

u(ξ)χ(ξ, λ) dξ.

(7.9)

 

2

 

−∞

С другой стороны, подставляя (7.8) в уравнение Шредингера (7.6), получим

Φxx = (χxx2iλχx−λ2χ)e−iλx = u(x)χe−iλx,

 

то есть

 

χxx2iλχx+= 0.

(7.10)

Свойства данных рассеяния

43

Продифференцируем теперь обе части интегрального уравнения (7.9)

x

χx = e2(x−ξ)u(ξ)χ dξ, χxx = −u(x)χ(λ, x)

−∞ x

2e2(x−ξ)u(ξ)χ(λ, ξ)= −u(x)χ+2iλχx,

−∞

так что выполнено дифференциальное уравнение (7.10) для функции χ(λ, x) и формула (7.7) доказана.

Выберем Im λ > 0. Тогда в интегральном уравнении (7.11) ядро аналитично и ограничено по λ, поскольку

| e2(x−ξ) |= e2Im λ(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 , Θ)/2(свойство 3).

Лекция 8

1. Свойства данных рассеяния (продолжение)

Проверим ограниченность a(λ) при λ → ∞. Из интегральных урав-

нений для Φ и Θ имеем:

 

 

 

Θ(λ, x) = e−iλx + o(1),

Θx = eiλx+O(1),

Φ(λ, x) = e−iλx + o(1),

Φx = −iλe−iλx+O(1),

так что

1

 

 

a(λ) =

W , Θ) 1, λ → ∞.

2

Свойства функции a(λ) проверяются аналогичным образом, выполняя комплексное сопряжение в интегральных уравнениях для Θ и Φ.

Свойство 4 доказано.

ПРИМЕР 1: ДЕЛЬТА-ОБРАЗНЫЙ ПОТЕНЦИАЛ.

Рассмотрим пример вычисления данных рассеяния для одного специального потенциала – дельта-функции Дирака u(x) = δ(x). Это обобщенная функция, которая может быть определена так:

 

 

0, x < −2ε ,

 

 

δ(x) = lim

1ε , ε < x < ε ,

(8.1)

 

0, x > ε .

 

 

ε 0

 

2

2

 

 

 

2

 

 

 

 

 

 

 

Можно представлять себе дельта-функцию как «бесконечно высокую стенку бесконечно малой толщины». В физике с ее помощью удобно записывать пространственную плотность различных величин (массы, заряда, интенсивность источника тепла, силу и т.п.), сосредоточенных в точке x0. Так, плотность точечной массы m, находящейся в точке a, записывается в виде (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λxsin λ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λxe−iλx

= (1+

1

 

)e−iλx

1

eiλx,

 

 

 

 

 

 

2

 

2

 

2

 

 

так что по определению (Φ

a(λ)e−iλx + b(λ)eiλx,

x

+

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

a(λ) = 1+

 

, b(λ) =

 

.

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

В заключение данного примера отметим, что данные рассеяния можно вычислить в явном виде для любого «ступенчатого» потенциала вида

 

 

 

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

Подстановка в уравнение немедленно дает C = a , откуда по определению данных рассеяния имеем

a(λ) =

C − 1

=

iλ − a

, b(λ)

0.

 

C + 1

 

+ 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(ξ)

= 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(ξ)|

=

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

e2(x−ξ)

2u(ξ)χ(λ, ξ) dξ.

 

−∞

Используя аналитичность χ по λ при Im λ > 0, перейдем к пределу при λ → ∞. Воспользовавшись оценкой

| e2(x−ξ) |= e2Im λ(x−ξ) 0, λ → ∞, x > ξ,

50

 

 

Лекция 9

 

 

 

получаем

 

 

 

 

 

 

 

 

 

2(χ(λ, x)1) =

2(Φ(λ, 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)!=

 

 

a(λ)

 

 

Φ(λ, x)

 

 

−∞

 

 

+

 

 

 

 

=−∞ eiλy #

Θ(λ, x)e−iλx+r(λ)Θ(λ, x)$ dλ.

(9.9)

Интегралы в последней формуле определены корректно, поскольку их сходимость на бесконечности следует из оценки подынтегральных функций O λ1 exp((x − y)) , которая, в свою очередь, обеспечена асимптотиками (9.3), (9.4), (9.7) и свойством 4.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]