Метод обратной задачи рассеяния: как решают КдФ и НУШ

Метод обратной задачи рассеяния (МОЗР, по-английски inverse scattering transform) - это способ точно решать целый класс нелинейных уравнений в частных производных, которые «вручную» не интегрируются. Идея почти парадоксальная: вместо того чтобы атаковать нелинейное уравнение в лоб, его связывают со вспомогательной линейной задачей на собственные значения и переводят всю эволюцию в пространство данных рассеяния, где она становится тривиальной. Метод родился в 1967 году из работы Гарднера, Грина, Крускала и Миуры по уравнению Кортевега - де Фриза и стал «нелинейным аналогом преобразования Фурье». Разберём логику метода по шагам: что такое данные рассеяния, зачем нужна пара Лакса и почему на выходе так естественно появляются солитоны.

Если нужно не просто прочитать, а разобрать конкретное уравнение или вывести солитонное решение - соберите задачу в форме ниже, и решение придёт пошагово.

Зачем вообще нужен метод

Линейные уравнения в физике решает преобразование Фурье: раскладываем начальное условие по гармоникам, каждая гармоника эволюционирует независимо и просто, а потом собираем ответ обратно. Для нелинейных уравнений этот приём ломается - гармоники начинают взаимодействовать, принцип суперпозиции не работает.

Метод обратной задачи рассеяния повторяет ту же трёхшаговую схему «прямое преобразование - простая эволюция - обратное преобразование», но роль гармоник Фурье играют данные рассеяния вспомогательного линейного оператора. Получается удивительное: целый класс нелинейных уравнений оказывается точно интегрируемым - у них есть бесконечно много законов сохранения, а начальная задача решается аналитически.

Рассчитать для своей задачи

Ключевые «клиенты» метода - уравнение Кортевега - де Фриза (КдФ) для волн на мелкой воде, нелинейное уравнение Шрёдингера (НУШ) в нелинейной оптике, уравнение синус-Гордона, уравнение Тоды. Все они описывают среды, где дисперсия уравновешивается нелинейностью, и именно в этом балансе рождаются солитоны.

Пара Лакса: мост к линейной задаче

Центральная конструкция метода - пара Лакса. Питер Лакс в 1968 году заметил, что интегрируемое нелинейное уравнение можно записать как условие совместности двух линейных операторов $L$ и $A$, зависящих от искомой функции $u(x,t)$:

$$\frac{\partial L}{\partial t} = [A, L] = AL - LA.$$

Здесь $L$ - оператор задачи на собственные значения (для КдФ это оператор Шрёдингера), а $A$ задаёт эволюцию его собственных функций по времени. Главное следствие: если $u$ удовлетворяет нелинейному уравнению, то собственные значения оператора $L$ не зависят от времени. Они - те самые сохраняющиеся величины, аналог амплитуд Фурье-гармоник.

Для уравнения КдФ

$$u_t - 6 u u_x + u_{xxx} = 0$$

оператором $L$ служит одномерный оператор Шрёдингера

$$L = -\frac{\partial^2}{\partial x^2} + u(x,t),$$

в котором решение $u(x,t)$ играет роль потенциала. Так нелинейная динамика волны превращается в квантово-механическую задачу рассеяния на изменяющемся со временем потенциале. О том, что описывает уравнение Шрёдингера, полезно вспомнить отдельно - здесь его оператор работает как технический инструмент.

Прямая задача рассеяния

Первый шаг схемы - прямое преобразование рассеяния. Берём начальное условие $u(x,0)$ как потенциал в операторе $L$ и решаем для него стационарную задачу рассеяния. У потенциала, спадающего на бесконечности, спектр распадается на две части.

Дискретный спектр - связанные состояния с собственными значениями $\lambda_n = -\kappa_n^2$. Каждое такое состояние отвечает отдельному солитону: чем глубже уровень, тем выше и быстрее солитон. Им сопоставляются нормировочные коэффициенты $c_n$.

Непрерывный спектр описывается коэффициентами отражения $r(k)$ и прохождения. Физически это «излучение» - диспергирующая часть волны, которая со временем растекается и убывает.

Совокупность $\{\kappa_n, c_n, r(k)\}$ называется данными рассеяния. Это и есть нелинейный аналог Фурье-образа: полное описание состояния системы в «спектральном» пространстве, где эволюция проста.

Рассчитать для своей задачи

Эволюция данных рассеяния

Второй шаг - самый изящный. Подставив пару Лакса, можно явно выписать, как данные рассеяния зависят от времени. Для КдФ собственные значения $\kappa_n$ остаются постоянными, а коэффициенты эволюционируют по простым линейным законам:

$$c_n(t) = c_n(0)\, e^{4\kappa_n^3 t}, \qquad r(k,t) = r(k,0)\, e^{8 i k^3 t}.$$

Вся нелинейность ушла: в спектральном пространстве каждый «режим» развивается независимо по элементарному экспоненциальному закону, ровно как гармоники в обычном преобразовании Фурье. Дискретные данные растут или осциллируют, непрерывные просто набирают фазу. Именно ради этой простоты и затевался весь обходной манёвр.

Обратная задача: уравнение Гельфанда - Левитана - Марченко

Третий шаг - обратное преобразование рассеяния: по данным рассеяния в момент $t$ восстановить потенциал $u(x,t)$, то есть само решение. Это классическая обратная задача, решаемая интегральным уравнением Гельфанда - Левитана - Марченко (ГЛМ).

Из данных рассеяния собирают ядро

$$F(x) = \sum_n c_n^2\, e^{-\kappa_n x} + \frac{1}{2\pi}\int_{-\infty}^{\infty} r(k)\, e^{i k x}\, dk,$$

затем решают линейное интегральное уравнение относительно функции $K(x,y)$:

$$K(x,y) + F(x+y) + \int_x^{\infty} K(x,z)\, F(z+y)\, dz = 0,$$

и наконец восстанавливают решение по диагонали ядра:

$$u(x,t) = -2\,\frac{d}{dx} K(x,x).$$

Уравнение ГЛМ - линейное, и в этом вся соль метода: нелинейная задача целиком сведена к цепочке линейных операций. Тот же безотражательный потенциал, кстати, лежит в основе оптического солитона в волокне - там НУШ играет роль КдФ.

Откуда берутся солитоны

Самый наглядный результат метода - солитоны. Если в начальном условии непрерывный спектр отсутствует ($r(k)\equiv 0$), остаются только дискретные уровни - это безотражательный потенциал, и он порождает чистое многосолитонное решение. Для одного уровня $\kappa$ уравнение ГЛМ решается элементарно и даёт классический солитон КдФ:

$$u(x,t) = -2\kappa^2\, \operatorname{sech}^2\!\big(\kappa\,(x - 4\kappa^2 t) - \delta\big).$$

Из формулы сразу видны фирменные свойства: амплитуда $2\kappa^2$ и скорость $4\kappa^2$ жёстко связаны (высокий солитон бежит быстрее), а форма $\operatorname{sech}^2$ устойчива. При нескольких уровнях получаются многосолитонные решения: солитоны сталкиваются, проходят друг сквозь друга и выходят неизменными - лишь со сдвигом фазы. Эта «частицеподобность» и дала им название.

Рассчитать для своей задачи

Связь с другими уравнениями

Сила метода в его универсальности. Каждое интегрируемое уравнение получает свою пару Лакса и свою линейную задачу рассеяния:

• НУШ $i\psi_t + \psi_{xx} + 2|\psi|^2\psi = 0$ - система Захарова - Шабата (матричный оператор $2\times2$); описывает огибающие волновых пакетов и оптические солитоны в волокне.
• Синус-Гордона $u_{tt} - u_{xx} + \sin u = 0$ - даёт топологические солитоны (кинки) и бризеры.
• Уравнение Тоды - дискретная цепочка, обратная задача формулируется для матриц Якоби.

Во всех случаях работает одна и та же схема трёх шагов, меняется лишь конкретный спектральный оператор. Поэтому МОЗР часто называют единой теорией интегрируемых систем.

Частые ошибки

• Путают прямую и обратную задачу. Прямая - по потенциалу найти данные рассеяния; обратная - наоборот. Эволюция во времени делается именно с данными рассеяния, между двумя этими шагами.
• Считают, что метод универсален. Он работает только для интегрируемых уравнений, у которых есть пара Лакса. Большинство нелинейных уравнений ею не обладают.
• Думают, что собственные значения меняются со временем. Для КдФ дискретные $\kappa_n$ строго постоянны - это и есть законы сохранения; меняются только нормировочные коэффициенты.
• Отождествляют солитон с любой уединённой волной. Солитон обязан сохранять форму при столкновениях; просто горб, расплывающийся со временем, солитоном не является.
• Забывают про непрерывный спектр. Часть энергии уходит в диспергирующее излучение; чисто солитонным решение бывает только у безотражательных потенциалов.

FAQ

Чем метод обратной задачи рассеяния отличается от преобразования Фурье? Структура та же - прямое преобразование, простая эволюция, обратное преобразование, - но Фурье линеен и раскладывает по гармоникам, а МОЗР нелинеен и раскладывает по данным рассеяния спектрального оператора. Фурье работает для линейных уравнений, МОЗР - для интегрируемых нелинейных.

Для каких уравнений он применим? Только для интегрируемых: КдФ, НУШ, синус-Гордона, модифицированного КдФ, уравнения Тоды и ряда других. Признак применимости - существование пары Лакса. Для произвольного нелинейного уравнения метод не работает.

Обязательно ли решать уравнение Гельфанда - Левитана - Марченко численно? Для безотражательных потенциалов (чистые солитоны) оно решается аналитически в замкнутом виде. В общем случае с непрерывным спектром обратную задачу решают численно или через задачу Римана - Гильберта.

Коротко

Метод обратной задачи рассеяния решает интегрируемые нелинейные уравнения по схеме «прямое преобразование - линейная эволюция - обратное преобразование», нелинейным аналогом Фурье в которой выступают данные рассеяния вспомогательного оператора. Связь нелинейного уравнения с линейной задачей на собственные значения обеспечивает пара Лакса; для КдФ роль спектрального оператора играет оператор Шрёдингера. Прямая задача даёт дискретный (солитоны) и непрерывный (излучение) спектр, эволюция данных во времени линейна, а обратная задача через уравнение Гельфанда - Левитана - Марченко восстанавливает решение. Дискретные уровни безотражательного потенциала дают точные многосолитонные решения, устойчивые при столкновениях.