Integrable Systems and Inverse Scattering

⚡ Key Takeaways

Integrable systems are nonlinear evolution equations (ODEs/PDEs) exactly solvable via the inverse scattering transform (IST) — a nonlinear analogue of the Fourier transform. The key result: equations like KdV and NLS possess infinitely many conserved quantities and can be linearized by mapping the solution onto scattering data of an associated linear spectral problem, which evolves trivially. Applies to shallow-water waves, nonlinear optics, plasma physics, and Bose–Einstein condensates. Pioneered by Gardner, Greene, Kruskal & Miura (1967) and Lax (1968).

💡 The Core Insight

A nonlinear PDE $u_t = K(u)$ is integrable if it is equivalent to the statement that two linear operators $L$ and $M$ satisfy the Lax equation $L_t = [M, L]$ . This forces the eigenvalues of $L$ to be time-independent — they are the conserved quantities. The nonlinear flow on $u$ becomes a linear, isospectral flow on the spectrum of $L$ , and the nonlinear problem dissolves into the linear machinery of scattering theory.

DefinitionLax Pair

A Lax pair for an evolution equation $u_t = K(u)$ is a pair of linear operators $(L, M)$ , depending on the field $u(x,t)$ , such that the equation is equivalent to the Lax equation:

1

\quad \frac{\partial L}{\partial t} = [M, L] = ML - LM.

Derivation hint: Differentiate the eigenvalue equation $L\psi = \lambda\psi$ in $t$ and demand $\lambda_t = 0$ ; the compatibility of $\psi_t = M\psi$ with isospectrality yields exactly this commutator equation. Here $L$ is the spatial (scattering) operator whose spectrum encodes conserved data, and $M$ generates the time evolution of the eigenfunctions $\psi$ .

DefinitionKorteweg–de Vries Equation and its Lax Operators

The canonical example is the KdV equation:

2

\quad u_t - 6\,u\,u_x + u_{xxx} = 0,

Derivation hint: Balance nonlinear steepening $-6uu_x$ against dispersion $u_{xxx}$ ; rescaling absorbs constants. The " $-6$ " is a normalization that makes the scattering coefficient of $u$ equal to one in $L$ . with Lax pair given by the Schrödinger operator and a third-order operator:

3

\quad L = -\partial_x^2 + u(x,t), \qquad M = -4\,\partial_x^3 + 6\,u\,\partial_x + 3\,u_x.

Derivation hint: Choose $M$ as the most general antisymmetric third-order differential operator; computing $L_t - [M,L]$ produces a multiplication operator equal to $u_t - 6uu_x + u_{xxx}$ , which must vanish.

TheoremIsospectrality

If $u(x,t)$ solves the KdV equation and $L = -\partial_x^2 + u$ , then the discrete eigenvalues $\lambda_n = -\kappa_n^2$ of $L$ (the bound-state energies) are constants of the motion:

4

\quad \frac{d\lambda_n}{dt} = 0, \qquad n = 1, \dots, N.

Derivation hint: From $L\psi = \lambda\psi$ , apply $\partial_t$ , use the Lax equation $L_t = [M,L]$ , and project onto the normalized eigenfunction $\psi$ ; the commutator term integrates to zero, leaving $\lambda_t \langle\psi,\psi\rangle = 0$ .

TheoremInverse Scattering Transform — Solution Scheme

For KdV with rapidly decaying initial data $u(x,0) \to 0$ as $|x|\to\infty$ , the scattering data

S(t) = \big\{ R(k,t),\; \{\kappa_n\}_{n=1}^N,\; \{c_n(t)\}_{n=1}^N \big\}

(reflection coefficient $R$ , bound-state momenta $\kappa_n$ , normalization constants $c_n$ ) evolve linearly:

5

\quad R(k,t) = R(k,0)\,e^{8ik^3 t}, \qquad c_n(t) = c_n(0)\,e^{4\kappa_n^3 t}, \qquad \kappa_n = \text{const}.

Derivation hint: Insert the asymptotic form of $\psi$ as $|x|\to\infty$ (where $u\to 0$ , so $M \to -4\partial_x^3$ ) into $\psi_t = M\psi$ ; matching the plane-wave coefficients $e^{\pm ikx}$ gives ODEs $\dot R = 8ik^3 R$ etc.

ProofThe Three-Stage IST Pipeline

We solve $u_t - 6uu_x + u_{xxx} = 0$ with given $u(x,0)$ .

Step 1 — Direct scattering ( $t=0$ ). Solve the time-independent Schrödinger equation

6

\quad -\psi_{xx} + u(x,0)\,\psi = k^2\,\psi.

Derivation hint: This is $L\psi = \lambda\psi$ with $\lambda = k^2$ ; the continuous spectrum $\lambda>0$ gives scattering states, $\lambda<0$ gives bound states. Define the Jost solutions $\psi \sim e^{-ikx}$ as $x\to -\infty$ , with asymptotics $a(k)e^{-ikx} + b(k)e^{ikx}$ as $x\to+\infty$ . The scattering data are the transmission/reflection coefficients $T = 1/a$ , $R = b/a$ , and the bound states are the zeros of $a(k)$ on the positive imaginary axis $k = i\kappa_n$ .

Step 2 — Time evolution. By the Isospectrality Theorem and the linear ODEs of equation [5], propagate the scattering data forward to time $t$ . The bound-state momenta $\kappa_n$ are frozen; only phases and normalizations move.

Step 3 — Inverse scattering (reconstruction). Form the kernel

7

\quad B(\xi; t) = \sum_{n=1}^N c_n^2(t)\,e^{-\kappa_n \xi} + \frac{1}{2\pi}\int_{-\infty}^{\infty} R(k,t)\,e^{ik\xi}\,dk,

Derivation hint: $B$ is the Fourier-plus-bound-state transform of the scattering data; it serves as the inhomogeneity of the Marchenko equation below. and solve the Gelfand–Levitan–Marchenko (GLM) integral equation for $K(x,y;t)$ , $y>x$ :

8

\quad K(x,y;t) + B(x+y;t) + \int_x^\infty K(x,z;t)\,B(z+y;t)\,dz = 0.

Derivation hint: $K$ is the kernel of the transformation operator $\psi = (1+K)e^{-ikx}$ mapping free solutions to Jost solutions; substituting into the Schrödinger equation and using completeness yields this linear equation. Finally, recover the solution by differentiating along the diagonal:

9

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

Derivation hint: Expand the Schrödinger equation for the transformation operator to leading order; the potential is recovered as twice the $x$ -derivative of the diagonal of $K$ . Because every step after the nonlinear Step 1 is linear, the nonlinear PDE has been completely solved. $\blacksquare$

∎

PL

Peter Lax

1926–

Hungarian-American mathematician at the Courant Institute. In his 1968 paper "Integrals of nonlinear equations of evolution and solitary waves," Lax abstracted the GGKM discovery into the operator formalism $L_t = [M,L]$ that now bears his name — the Lax pair. This single algebraic reformulation revealed integrability as a structural property: any equation expressible as an isospectral deformation is solvable by inverse scattering. Lax received the Abel Prize (2005) in part for this work, which unified soliton theory, spectral theory, and Hamiltonian mechanics.

🧠 Intuition

Think of the inverse scattering transform as a nonlinear Fourier transform. The ordinary Fourier transform diagonalizes linear constant-coefficient PDEs: each mode $e^{ikx}$ evolves independently. The IST does the same for special nonlinear PDEs, but the "modes" are now of two species: a continuous spectrum (the reflection coefficient $R(k)$ , behaving like dispersive radiation) and a discrete spectrum (the bound states $\kappa_n$ , which are the solitons). In this nonlinear Fourier space, the dynamics trivializes — modes just rotate by phases. The deep magic is that the discrete eigenvalues are conserved: each soliton carries a fixed eigenvalue $-\kappa_n^2$ like a quantum number, so solitons survive collisions and merely exchange phase shifts. The nonlinear PDE is thus secretly a collection of non-interacting oscillators in the right coordinates.

Figure

The IST commuting diagram: nonlinear evolution is conjugated to linear flow in scattering space

ExampleSingle-Soliton from a Reflectionless Potential

Given: Initial data with one bound state, $R(k,0)=0$ (reflectionless), $\kappa_1 = \kappa$ , normalization $c_1^2(0) = 2\kappa\,e^{2\kappa x_0}$ .

Step 1. Time-evolve: $c_1^2(t) = 2\kappa\,e^{2\kappa x_0}\,e^{8\kappa^3 t}$ .

Step 2. Kernel (no continuous part): $B(\xi;t) = c_1^2(t)\,e^{-\kappa\xi}$ .

Step 3. With a single exponential kernel, the GLM equation [8] has a degenerate (separable) solution $K(x,y) = -c_1^2 e^{-\kappa y}/(1 + \tfrac{c_1^2}{2\kappa}e^{-2\kappa x})$ . Setting $y=x$ :

K(x,x;t) = -\frac{c_1^2(t)\,e^{-2\kappa x}}{1 + \frac{c_1^2(t)}{2\kappa}e^{-2\kappa x}}.

Step 4. Apply equation [9], $u = -2\,\partial_x K(x,x;t)$ . After simplification:

10

\quad u(x,t) = -2\kappa^2\,\operatorname{sech}^2\!\big[\kappa(x - 4\kappa^2 t - x_0)\big].

Derivation hint: Substitute $c_1^2(t)$ , complete the algebra, and recognize $\operatorname{sech}^2$ via $\frac{d}{dx}\tanh = \operatorname{sech}^2$ ; the phase shift $x_0$ absorbs $e^{2\kappa x_0}$ and the velocity $4\kappa^2$ emerges from the $e^{8\kappa^3 t}$ factor.

\boxed{\;u(x,t) = -2\kappa^2\,\operatorname{sech}^2[\kappa(x - 4\kappa^2 t - x_0)]\;}

Physical reading: amplitude $= 2\kappa^2$ , speed $= 4\kappa^2$ — taller solitons travel faster. (Units: with $x$ in length $L$ and $t$ in time, $u \sim L^{-2}$ , $\kappa \sim L^{-1}$ , speed $\sim L^2/t$ after dimensional rescaling.)

ExampleConserved Quantities of KdV

Given: KdV $u_t = 6uu_x - u_{xxx}$ (sign-flipped convention) with $u\to 0$ at $\pm\infty$ . Find the first three conserved densities.

Step 1. Mass: integrate the PDE directly. Since $u_t = (3u^2 - u_{xx})_x$ ,

\frac{d}{dt}\int u\,dx = \big[3u^2 - u_{xx}\big]_{-\infty}^{\infty} = 0 \;\Rightarrow\; I_1 = \int u\,dx.

Step 2. Momentum: multiply by $u$ and integrate; the flux is again a perfect $x$ -derivative, giving

I_2 = \int u^2\,dx, \qquad \frac{dI_2}{dt}=0.

Step 3. Energy (Hamiltonian): the third invariant is

11

\quad I_3 = \int \Big( u^3 + \tfrac{1}{2}u_x^2 \Big)\,dx.

Derivation hint: Construct the density via the Miura-linked Gardner expansion or by guessing $\alpha u^3 + \beta u_x^2$ and fixing $\alpha,\beta$ so $\partial_t(\text{density})$ is a total $x$ -derivative; this Hamiltonian generates KdV under the Gardner bracket.

\boxed{\;I_1=\!\int\! u\,dx,\quad I_2=\!\int\! u^2\,dx,\quad I_3=\!\int\!\big(u^3+\tfrac12 u_x^2\big)dx\;}

These are the first three of an infinite hierarchy $I_1, I_2, I_3, \dots$ , equivalently encoded as the coefficients in the large- $k$ asymptotic expansion of $\ln a(k)$ (the trace formula).

⚠️ Common mistake

Integrability is rare and non-generic. Adding an arbitrary nonlinear term to KdV almost always destroys the Lax structure and the infinite conservation laws. Integrable PDEs are measure-zero islands in the space of nonlinear equations — do not expect a "nearby" perturbed equation to be solvable by IST. Likewise, "integrable" in the soliton sense (Liouville/Lax integrability, exact IST solvability) is distinct from numerical solvability or from "integrable" meaning merely "has a solution."

⚠️ Common mistake

The IST is not a Fourier transform of $u$ . Step 1 (direct scattering) is genuinely nonlinear in $u$ — the map $u \mapsto S$ depends on $u$ in a complicated way. Only the evolution and the inverse steps are linear. People often mistakenly think solitons can be superposed; they cannot be added linearly, though their scattering data combine additively in the discrete spectrum.

Counterexample

Burgers' equation is NOT integrable in the IST sense — but is still exactly solvable. Consider

u_t + u u_x = \nu\,u_{xx}.

The Cole–Hopf transformation $u = -2\nu\,(\ln\phi)_x$ linearizes it to the heat equation $\phi_t = \nu\phi_{xx}$ . This is exact solvability, yet Burgers has no solitons, no infinite conservation-law hierarchy, and no Lax pair with isospectral discrete spectrum — its dissipative term smears structure rather than balancing it. The lesson: "exactly linearizable" $\neq$ "integrable in the soliton/IST sense." KdV's conservative dispersion is what permits solitons; Burgers' viscosity forbids them.

connection

Integrable systems sit at the crossroads of analysis, geometry, and algebra. The Lax equation $L_t=[M,L]$ is a special case of the <span class="wiki-link-missing">Zero-Curvature Representation</span> $\partial_t U - \partial_x V + [U,V] = 0$ , which generalizes IST to the AKNS scheme covering the Nonlinear Schrödinger Equation and Sine-Gordon Equation. The conserved quantities and Poisson structure connect to Liouville Integrability and bi-Hamiltonian geometry (Magri's theorem). The Gelfand–Levitan–Marchenko equation is a Riemann–Hilbert Problem in disguise, which is the modern analytic engine of the field. On the algebraic side, the hierarchy structure links to <span class="wiki-link-missing">Kac–Moody Algebras</span> and tau functions of the Sato Grassmannian, while finite-gap (algebro-geometric) solutions connect to Riemann Theta Functions on hyperelliptic curves.

🔬 Where this appears

⏳ History

- 1834 — John Scott Russell observes the "Wave of Translation," a solitary hump persisting unchanged along the Union Canal, Scotland — the first soliton sighting.

- 1895 — Diederik Korteweg and Gustav de Vries derive the KdV equation and find its solitary-wave

\operatorname{sech}^2

solution, vindicating Russell.

- 1955 — Fermi, Pasta, Ulam & Tsingou numerically study an anharmonic chain expecting thermalization; instead observe recurrence — energy returns to initial modes (the FPUT paradox).

- 1965 — Norman Zabusky and Martin Kruskal simulate KdV as the continuum limit of FPUT, discover that solitary waves pass through each other intact, and coin the term "soliton."

- 1967 — Gardner, Greene, Kruskal, and Miura (GGKM) solve KdV exactly via the inverse scattering method, linking it to the Schrödinger operator's spectrum.

- 1968 — Peter Lax introduces the Lax pair formalism

L_t=[M,L]

, abstracting integrability into operator language.

- 1971–72 — Vladimir Zakharov and Aleksei Shabat solve the Nonlinear Schrödinger equation by IST; Ablowitz, Kaup, Newell & Segur (AKNS, 1973–74) unify these via the

2\times2

zero-curvature scheme.

- 1974 — Faddeev, Zakharov, and Gardner reveal KdV as a completely integrable Hamiltonian system with action-angle variables (the scattering data).

- 1970s–80s — Sato, Date, Jimbo, Kashiwara, Miwa develop the tau-function / Grassmannian theory; algebro-geometric (finite-gap) solutions via theta functions (Novikov, Dubrovin, Matveev, Krichever) emerge.

open-problem

Where the theory strains and what remains open. (1) Beyond rapidly decaying data: IST for KdV requires $u\to 0$ fast enough that the scattering problem is well-posed; for periodic or algebro-geometric initial data the appropriate objects are Bloch spectra and Riemann surfaces, and the analysis is far subtler. The Deift–Zhou nonlinear steepest descent method (1990s) rigorously extracts long-time asymptotics from the Riemann–Hilbert formulation, but uniform asymptotics in transition regions remain delicate. (2) Higher dimensions: genuine integrable PDEs in $2+1$ dimensions exist (KP, Davey–Stewartson) requiring the $\bar\partial$ (d-bar) method, but a clean IST in $3+1$ dimensions is essentially absent — integrability appears to be a low-dimensional miracle. (3) Perturbed and "near-integrable" systems: how solitons respond to small non-integrable perturbations (radiation, soliton decay) is handled by perturbation theory built on IST, but its rigorous justification and the fate of the FPUT recurrence under the thermodynamic limit connect to deep open questions in statistical mechanics and nonlinear wave turbulence. (4) Classification: there is still no complete, constructive classification of all integrable PDEs; the symmetry approach (Mikhailov, Sokolov) provides powerful necessary conditions but no exhaustive list.

Integrable Systems and Inverse Scattering

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also