🌀 nonlinear sciences

Perturbation theory for kinks of the defocusing modified Korteweg-de Vries equation

This paper develops an integrable perturbation theory for defocusing modified Korteweg-de Vries kinks using squared eigenfunction expansions to derive evolution equations for perturbed parameters and demonstrate that generic perturbations generate a radiative shelf, a finding validated by numerical simulations.

Original authors: Nicholas J. Ossi, Barbara Prinari, Theodoros P. Horikis, Dimitrios J. Frantzeskakis

Published 2026-06-29

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Nicholas J. Ossi, Barbara Prinari, Theodoros P. Horikis, Dimitrios J. Frantzeskakis

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a perfectly smooth, endless ocean where a single, massive wave is traveling. In the world of physics, this isn't just any wave; it's a "kink." Think of it as a permanent, stable ramp in the water that connects a low level to a high level. It doesn't crash or fade away; it just glides along, maintaining its shape forever. This is the "kink solution" of a famous mathematical equation called the defocusing modified Korteweg-de Vries (mKdV) equation.

In the real world, however, nothing is perfectly isolated. There is always a little bit of wind, a bit of friction, or some other tiny disturbance. The big question this paper asks is: What happens to our perfect, stable wave ramp when we poke it with a tiny stick?

Here is the breakdown of what the researchers discovered, using simple analogies:

1. The Problem: The Perfect Wave vs. Real Life

The math behind these waves is incredibly complex. Usually, when you disturb a perfect system, you have to throw out the old math and start over with messy, approximate calculations. The authors wanted a better way. They wanted to predict exactly how the wave would change its speed, where it would move, and what new ripples it would create, without losing the elegance of the original math.

2. The Tool: The "Squared Eigenfunction" Toolkit

To solve this, the team used a sophisticated mathematical toolkit called "squared eigenfunction expansion."

The Analogy: Imagine you have a complex sound (like a chord on a piano). To understand it, you break it down into individual notes (frequencies). In this paper, the "notes" are special mathematical shapes called squared eigenfunctions.
The Innovation: While scientists have used this toolkit for other types of waves for decades, this paper is the first to successfully adapt it specifically for the "kink" wave (the ramp) on a non-zero background. It's like finally figuring out how to tune a specific, rare instrument that everyone else had been ignoring.

3. The Discovery: The "Radiative Shelf"

When they applied their new toolkit to a perturbed kink, they found something surprising and consistent.

The Result: The wave doesn't just wiggle; it leaves a trail. Specifically, a flat, elevated "shelf" of water forms in front of the kink (on the left side, since the wave moves left).
The Metaphor: Imagine a snowplow moving down a street. As it pushes the snow, it doesn't just clear a path; it leaves a pile of snow in front of it. In this case, the "snow" is a flat, raised shelf of water that travels ahead of the main wave.
The Speed: This shelf moves faster than the kink itself. Because it's faster, the kink can never catch up to it or overtake it. The shelf is always running away ahead of the wave.

4. The Predictions: How the Wave Changes

The paper doesn't just say "a shelf appears." It gives a precise recipe for how the wave behaves under different types of "pokes" (perturbations):

If you add "diffusion" (like friction or spreading): The wave keeps its height, but it slows down and shifts its position. The shelf that forms in front has a specific, predictable height that doesn't depend on how big the wave was to begin with.
If you add "damping" (like energy loss): The wave actually shrinks in height over time, and the shelf that forms in front is different in size.
The "No-Shelf" Rule: The researchers found that if the disturbance is perfectly symmetrical (like certain types of higher-order ripples), the shelf might not form at all. It's like pushing a swing perfectly in rhythm; sometimes, nothing extra happens.

5. Why It Matters

The authors tested their math against computer simulations (essentially, a digital wind tunnel). The results matched perfectly.

The Takeaway: They have built a reliable "crystal ball" for these waves. If you know the type of tiny disturbance hitting the wave, you can now calculate exactly how the wave will speed up, slow down, shift position, and how big the "shelf" in front of it will be.

Summary

This paper is about taking a perfect, theoretical wave that exists only in math books and figuring out how it survives when the real world gets messy. They developed a new, precise method to predict that when you disturb this wave, it leaves a distinct, flat "shelf" of energy running ahead of it, and they can tell you exactly how big that shelf will be and how the wave will move. It turns a chaotic problem into a predictable, solvable one.

Technical Summary: Perturbation Theory for Kinks of the Defocusing Modified Korteweg-de Vries Equation

Problem Statement
The paper addresses the behavior of the kink solution to the defocusing modified Korteweg-de Vries (mKdV) equation under small, non-integrable perturbations. The defocusing mKdV equation, $q_t + 6\sigma q^2 q_x + q_{xxx} = 0$ with $\sigma = -1$ , admits a topological kink solution $q(x,t) = q_0 \tanh[q_0(x + 2q_0^2 t - x_0)]$ that interpolates between two distinct constant asymptotic states ( $q \to \pm q_0$ ). While the inverse scattering transform (IST) provides a complete solution for the unperturbed integrable case, systematic perturbation theories for the mKdV kink (as opposed to bright solitons or the KdV/NLS equations) have been notably sparse in the literature. The authors aim to develop an integrable perturbation theory to predict the slow-time evolution of the kink's parameters (velocity and center) and the emergence of radiative structures when the system is subjected to small perturbations such as diffusion or linear damping.

Methodology
The authors employ a squared eigenfunction expansion method based on the underlying Zakharov–Shabat scattering problem associated with the mKdV equation. The methodology proceeds through the following steps:

Scattering Framework: The study utilizes the IST for the defocusing mKdV with non-vanishing boundary conditions. The scattering problem is formulated using a uniformization variable $z$ to handle the branch points at $k = \pm q_0$ . For the pure kink solution, the scattering data contains a single discrete eigenvalue at $z = iq_0$ (and its conjugate), with no continuous spectrum reflection.
Completeness Relation: A central technical contribution is the derivation of the completeness relation for the squared eigenfunctions appropriate to the kink background. Adapting results from the defocusing Nonlinear Schrödinger (NLS) equation, the authors construct scalar squared eigenfunctions $\eta_s$ and adjoint squared eigenfunctions $\chi_s$ . They explicitly handle the singularities at the branch points $\pm q_0$ and the discrete eigenvalue $iq_0$ , establishing a closure relation that includes both a continuous integral over the real axis and a discrete residue term.
Perturbation Expansion: The perturbed equation $q_t - 6q^2 q_x + q_{xxx} = \epsilon F[q]$ is analyzed using a multiple-scales ansatz. The kink parameters ( $q_0$ and $x_0$ ) are allowed to evolve on a slow time scale $T = \epsilon t$ . The first-order correction $\tilde{q}$ is expanded in the basis of the squared eigenfunctions.
Modulation Equations: By projecting the linearized perturbation equation onto the adjoint squared eigenfunctions, the authors derive orthogonality conditions. These conditions eliminate secular growth in the first-order correction, yielding explicit evolution equations for the background amplitude $q_0(T)$ and the kink center $x_0(T)$ .
Radiative Shelf Analysis: The first-order correction integral is analyzed in the asymptotic limit. The authors demonstrate that the integral predicts the formation of a "radiative shelf" (a step-like structure) on the left side (front) of the kink, propagating with a velocity distinct from the kink itself.

Key Contributions and Results

Derivation of Completeness: The paper establishes the specific completeness and orthogonality relations for squared eigenfunctions on a kink background, accounting for the unique symmetries of the real mKdV equation (specifically the symmetry $z \to -z$ ) which distinguishes it from the complex NLS case.
Evolution Equations: Explicit leading-order evolution equations are derived for the kink's amplitude and position. The evolution of the amplitude is determined by the divergent part of the inner product with the perturbation, while the position evolution is determined by the convergent part.
Radiative Shelf Prediction: The study shows that generic perturbations produce a prominent shelf in front of the kink (the region $x < x_0$ ). The shelf propagates with a velocity of $-6q_0^2$ , which is strictly faster in magnitude than the kink's velocity ( $-2q_0^2$ ), preventing the kink from overtaking it. The height of the shelf is explicitly calculated as $\tilde{q} \sim -\alpha / (4q_0^3)$ , where $\alpha$ is a functional of the perturbation $F$ .
Specific Applications: The theory is applied to two physically relevant cases:
- Diffusive Perturbation ( $F = q_{xx}$ ): The background amplitude remains constant, while the kink center shifts. A shelf of height $\epsilon/3$ forms.
- Linear Damping ( $F = -q$ ): The background amplitude decays exponentially, and the kink center shifts accordingly. The shelf height depends inversely on the square of the background amplitude.
Conservation Laws: A relationship is established between the evolution of the renormalized energy, the background amplitude, and the shelf height ( $E_T = q_{0T} - \alpha$ ), indicating that energy loss in dissipative systems manifests as shelf formation.

Significance and Claims
The authors claim that this work provides a rigorous foundation for analyzing the response of mKdV kinks to a broad class of perturbations, filling a gap in the literature regarding dark solitons and kinks on non-zero backgrounds. The paper emphasizes that the squared eigenfunction expansion method is particularly elegant for this problem as it diagonalizes the linearized flow and allows for closed-form expressions of scattering data corrections.

The study confirms that while the kink solution is structurally stable, perturbations generically induce a radiative shelf. The predictions derived from the perturbation theory are shown to be consistent with direct numerical simulations for both diffusion and damping cases, validating the theoretical framework. The paper concludes modestly, noting that the framework can be extended to dark solitons (which possess four discrete eigenvalues) and potentially to integrable fractional mKdV equations, but does not speculate on applications beyond these mathematical extensions.