Nonlocal Modified Korteweg-de Vries Equation

• The nonlocal mKdV equation is an integrable nonlinear PDE that couples field values at (x,t) with those at (-x,-t) via PT-symmetric interactions.
• Its formulation uses Lax pair methods and modified inverse scattering techniques to derive explicit soliton, breather, and multi-soliton solutions.
• The nonlocality introduces novel spectral properties and practical implications for PT-symmetric optics and nonlinear wave theory.

The nonlocal modified Korteweg-de Vries (mKdV) equation refers to a class of integrable nonlinear partial differential equations in which the nonlinear term couples the field at $(x,t)$ to its “reflected” counterpart at $(-x, -t)$, typically manifesting parity-time (PT) symmetry. The canonical form is

$$q_t(x,t) + q_{xxx}(x,t) - 6\,q(x,t)\,q(-x,-t)\,q_x(x,t) = 0,$$

which generalizes the local mKdV by incorporating nonlocality via reverse space–time reflection. This nonlocal feature induces novel spectral and dynamical properties and has substantial impact on inverse scattering theory, soliton/breather classification, and multi-soliton interactions. Recent developments include rigorous IST formulations for decaying, periodic, and step-like backgrounds, as well as discrete and coupled extensions.

1. Equations, Nonlocal Reductions, and PT Symmetry

The nonlocal mKdV equation emerges from the cubic AKNS system under an Ablowitz–Musslimani type reduction. For real or complex fields $q(x,t)$ and reduction $r(x,t) = \pm q(-x, -t)$, the nonlocal form is

$q_t + q_{xxx} - 6\,q(x,t)\,q(-x,-t)\,q_x(x,t) = 0,$

where the nonlocal cubic interaction replaces the local $q^2 q_x$ by $q(x,t) q(-x,-t) q_x(x,t)$, enforcing invariance under $(x,t) \mapsto (-x,-t)$ (Ji et al., 2016, Khare et al., 2023). In its complex generalization, the PT-symmetric constraint may be combined with complex conjugation, $q^*(-x,-t)$, yielding a nonlocal cmKdV.

For discrete analogues, reverse-lattice involutions (e.g., $Q_{n,m}$ coupled to $Q_{-n,-m}$) encode the nonlocal property (Zhao et al., 2023, Hu et al., 2024). Coupled nonlocal mKdV systems extend this to multi-component fields with cross-terms involving both local and nonlocal interactions (Khare et al., 2023).

2. Lax Pair Formulations and Integrability

The nonlocal mKdV maintains complete integrability and admits a cubic Lax pair structure analogous to the local theory, but with nonlocal entries. The standard $2 \times 2$ Lax pair reads

$$\partial_x \varphi = (-ik\sigma_3 + Q(x,t))\,\varphi, \quad \partial_t \varphi = \left(-4ik^3\sigma_3 + 4k^2 Q - 2ik V_1 + V_2\right) \varphi,$$

where $Q(x,t)$ includes dependence on $q(x,t)$ and $q(-x,-t)$ (Ji et al., 2016, He et al., 2018). The nonlocal structure ensures the compatibility condition (zero curvature) reproduces the nonlocal mKdV equation.

Discrete and semi-discrete extensions carry integrability via lattice Lax pairs and associated Casoratian/Wronskian solution formulae (Zhao et al., 2023, Deng et al., 2024). Coupled and Alice–Bob variants realize nonlocality via symmetry-constrained off-diagonal blocks in the Lax pair, guaranteeing the existence of infinite conservation laws (Li et al., 2017).

3. Inverse Scattering Transform (IST), Riemann–Hilbert Problems, and Soliton Formulae

IST for nonlocal mKdV requires modifications of the direct/inverse scattering data to account for nonlocality. The spectral analysis involves a pair of Jost functions analytic in opposing half-planes, obeying additional symmetry relations (e.g., $a(k) = a^*(-k^*)$) and uncoupled discrete eigenvalues for $a(k)$ and $\bar a(k)$ (Ji et al., 2016, Zhang et al., 2021).

The reflectionless (pure soliton) case leads to finite-rank Gel'fand–Levitan–Marchenko systems with solutions expressible as quotients of determinants or sums over matrix elements; for example,

$$q(x,t) = -2\,\mathrm{tr}\left( [I_{\bar N}+EE^T]^{-1} \bar h\,\bar h^T \right),$$

with $E$ built from discrete eigenvalues and norming constants. The Riemann–Hilbert approach, including for NZBCs and step-like backgrounds, encodes the solution as an analytic matrix plus meromorphic corrections at discrete spectrum and background singularities (Zhang et al., 2018, Rybalko, 22 Jan 2026).

The general $N$-soliton solution takes determinant form, e.g.,

$$q(x,t) = -2i\, \frac{\det M^\natural(x,t)}{\det M(x,t)},$$

with $M^\natural$ and $M$ assembled from spectral data and norming constants, with signs modulating singularity and boundedness (Zhang et al., 2021).

4. Solution Taxonomy: Solitons, Breathers, Superposed Waves, and Backgrounds

Nonlocality generates a rich variety of solutions, including:

• Bright-like solitons: Exponentially localized pulses with amplitudes modulated by the difference in mirror eigenvalues; amplitude may grow/decay along rays (Ji et al., 2016).
• Singular solitons: Amplitude develops zeros or blow-ups at isolated times/positions, controlled by sign choices or norming constraints (Zhang et al., 2021).
• Kink/antikink and rational solutions: Step-like profiles connecting asymptotic plateaus, and power-law or rational-decaying solutions via degenerate eigenvalue limits (Khare et al., 2023, Khare et al., 2022).
• Breathers: Oscillatory-in-time structures associated to complex-conjugate spectral pairs, reducible to explicit hyperbolic/trigonometric expressions in special cases (Ji et al., 2016).
• Superposed periodic waves: Genuine cnoidal superpositions (e.g., kink–antikink or two pulses at shifted phases) are enabled by nonlocality and may be constructed via Jacobi–addition identities (Khare et al., 2022).
• Plane-wave modulated and coupled solitons: Nonlocal coupled systems admit complex plane-wave factors in solitary and periodic solutions, which are forbidden in the local systems (Khare et al., 2023).

On nonzero backgrounds (NZBCs), the IST produces soliton–kink, breathing bound states, and elastic/inelastic collision dynamics, with classification contingent on boundary phase and symmetry parameters (Zhang et al., 2018).

For step-like and oscillatory backgrounds, the Riemann–Hilbert formalism produces explicit two-soliton solutions, characterized by parameter regimes for amplitude and frequency ratios (e.g., $B < A/4$), and step oscillations persist asymptotically in the radiative regions (Rybalko, 22 Jan 2026).

5. Long-time Asymptotics, Dynamical Phenomena, and Collision Behavior

Asymptotic analysis for nonlocal mKdV employs the Deift–Zhou nonlinear steepest-descent method, with nuanced differences compared to the classical case. Leading order soliton/radiation decay rates are modulated by spectral symmetry factors (e.g., $t^{-\frac12+\nu}$), and the logarithmic phase in the solution reflects the extra nonlocal symmetry (He et al., 2018).

Collisions of solitons generally remain elastic for bounded solutions, with modified phase shifts and, in singular cases, periodic blow-up and collapse events. The ability to rotate characteristic lines alters propagation patterns and singularity locations, creating steerable collapse (Zhang et al., 2021). In discrete and semi-discrete cases, all soliton interactions are phase-shifting and preserve core amplitude and shape (Zhao et al., 2023, Deng et al., 2024).

On NZBCs and oscillatory backgrounds, long-time asymptotics in various similarity sectors display transitions: localized kinks/breathers emerge at critical velocities, followed by persistent background oscillations (Rybalko, 22 Jan 2026, Zhang et al., 2018).

6. Discrete, Semi-discrete, Coupled, and Alice–Bob Generalizations

Discrete, semi-discrete, and continuum nonlocal mKdV equations are unified via Casoratian/Wronskian constructions. For the fully discrete case, solution hierarchies are built from double Casoratians, recognizing parameter constraints induced by the reverse-space–time coupling (Zhao et al., 2023, Deng et al., 2024).

Coupled nonlocal mKdV equations utilize cross-couplings of field values at mirrored positions, yielding solution families not seen in local coupled mKdV—especially complex plane-wave modulated solitons and multi-component interplays (Khare et al., 2023).

Alice–Bob mKdV formalism generalizes nonlocality to shifted parity and delayed time reversals: $B(x,t) = -A(-x + x_0, -t + t_0),$ yielding integrable Darboux hierarchies with bright/dark solitons, complexitons, rogue waves, and multi-wave coherent structures, some breaking the base symmetry (Li et al., 2017).

7. Mathematical and Physical Significance, Applications, and Open Problems

The nonlocal mKdV equation expands the landscape of integrable nonlinear wave equations by allowing PT-symmetric interactions across spatial and temporal reflection. Mathematically, it yields new solution classes, parameter constraints, and IST techniques. Physically, implications arise for PT-symmetric field theories, complex nonlinear optics, waveguides with time-reversal invariance, and dispersive media with nonlocal effects (Ji et al., 2016, Khare et al., 2023).

Open problems include derivation and classification of full multi-soliton formulas for coupled systems, explicit Lax pairs for multi-component and discrete settings, conserved density hierarchies, and spectral stability under nonlocal perturbations. Moreover, discrete nonlocal versions offer rational, breatherlike, and rogue wave solutions absent singularities, suggesting further applications for lattice wave simulations and parity-time symmetric numerical benchmarking (Deng et al., 2024, Hu et al., 2024).