pith. sign in

arxiv: 2604.27513 · v1 · submitted 2026-04-30 · 🌊 nlin.SI · cs.NA· math.NA

Numerical inverse scattering transform for the coupled modified Korteweg-de Vries equation

Wen-Xin Zhang , Yong Chen This is my paper

Pith reviewed 2026-05-07 09:46 UTC · model grok-4.3

classification 🌊 nlin.SI cs.NAmath.NA
keywords inverse scatteringcoupled mKdVRiemann-Hilbert problemDeift-Zhou methodnumerical solutionlong-time asymptoticsintegrable PDEs
0
0 comments X

The pith

Numerical inverse scattering transform solves the coupled mKdV equation directly at chosen points in space and time.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a numerical framework to solve the coupled modified Korteweg-de Vries equation using the inverse scattering transform. It handles the associated 3x3 matrix Riemann-Hilbert problem by first computing scattering data with Chebyshev collocation and then applying Deift-Zhou steepest descent to deform the contours according to two stationary points. This divides the space-time plane into three regions for efficient numerical solution of the inverse problem. The method evaluates the solution at any prescribed point without time-stepping, which is advantageous for studying long-time behavior. Numerical experiments demonstrate that it captures the primary asymptotic characteristics of the solutions.

Core claim

The central discovery is that the NIST for the coupled mKdV equation, based on its 3x3 matrix Riemann-Hilbert problem, computes the solution directly at prescribed spatial and temporal points. The Deift-Zhou nonlinear steepest descent method deforms the oscillatory problem, with the phase function having exactly two stationary points symmetric about the origin, leading to a division of the (x,t)-plane into three regions with corresponding contour deformations. This approach avoids reliance on time-stepping and numerical experiments show it captures the main asymptotic features.

What carries the argument

The 3x3 matrix Riemann-Hilbert problem associated with the coupled mKdV equation and its deformation using the Deift-Zhou nonlinear steepest descent analysis around the two symmetric stationary points.

Load-bearing premise

The Deift-Zhou nonlinear steepest descent analysis extends to the 3x3 matrix Riemann-Hilbert problem without essential modifications despite its more complex jump matrix and scattering data.

What would settle it

Compare the NIST solution at large times with an independent high-resolution numerical integration of the coupled mKdV equation; mismatch in the solution profile would disprove the method's accuracy.

read the original abstract

This paper develops the numerical inverse scattering transform (NIST) framework for the coupled modified Korteweg-de Vries (mKdV) equation based on its associated Riemann-Hilbert problem. The coupled system gives rise to a $3\times3$ matrix-valued Riemann-Hilbert problem, whose jump matrix and scattering data have a more involved structure than in the scalar case. This matrix setting makes the extension of NIST to the coupled system nontrivial, both in the direct scattering computation and in the numerical solution of the inverse problem. Within this framework, the scattering data are first computed by solving the matrix direct scattering problem using a Chebyshev collocation method with suitable mappings. The Deift-Zhou nonlinear steepest descent method is then used to analyze and deform the oscillatory Riemann-Hilbert problem. In particular, the phase function admits two stationary points symmetric about the origin, and the analysis leads to a division of the $(x,t)$-plane into three regions with corresponding contour deformations. Compared with traditional numerical methods, the NIST computes the solution directly at prescribed spatial and temporal points without relying on time-stepping. Numerical experiments illustrate the performance of the proposed NIST in long-time simulations and indicate that it captures the main asymptotic features of the coupled mKdV solutions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. This paper develops the numerical inverse scattering transform (NIST) for the coupled modified Korteweg-de Vries equation via its associated 3x3 matrix Riemann-Hilbert problem. Scattering data are obtained by Chebyshev collocation on the direct scattering problem; the inverse problem is solved by applying the Deift-Zhou steepest-descent analysis, which identifies two symmetric stationary points and deforms the oscillatory RH problem into three regions of the (x,t)-plane. The resulting scheme evaluates the solution at arbitrary prescribed (x,t) points without time-stepping, and numerical experiments are presented to illustrate long-time performance and capture of main asymptotic features.

Significance. If the central claims hold, the work provides a nontrivial extension of NIST to a vector integrable system whose RH problem has a more complex matrix structure than the scalar case. The ability to compute solutions directly at fixed points and to access long-time asymptotics without marching in time would be a useful addition to the toolkit for studying coupled nonlinear wave equations.

major comments (2)
  1. [Abstract] Abstract: the central claim that 'numerical experiments illustrate the performance of the proposed NIST in long-time simulations and indicate that it captures the main asymptotic features' is unsupported by any error tables, convergence rates, comparisons against known solutions, or quantitative metrics. This absence is load-bearing for the assertion that the method works for the coupled system.
  2. [Deift-Zhou analysis] Description of the Deift-Zhou analysis: the extension to the 3x3 matrix RH problem is asserted to proceed by identifying two symmetric stationary points and dividing the (x,t)-plane into three regions with corresponding contour deformations, yet the manuscript does not explicitly verify that the more involved jump-matrix structure permits the standard g-function construction, lens opening, and decay estimates without new obstructions. This assumption underpins both the contour choice and the claim that asymptotics are captured.
minor comments (1)
  1. The notation for the coupled mKdV system, the scattering data, and the 3x3 jump matrix should be introduced with explicit component-wise definitions at the first appearance to improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, indicating the revisions we will incorporate to strengthen the presentation.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that 'numerical experiments illustrate the performance of the proposed NIST in long-time simulations and indicate that it captures the main asymptotic features' is unsupported by any error tables, convergence rates, comparisons against known solutions, or quantitative metrics. This absence is load-bearing for the assertion that the method works for the coupled system.

    Authors: We agree that the abstract claim would be more robustly supported by quantitative evidence. In the revised manuscript we will expand the numerical experiments section to include error tables (e.g., L^2 and L^infty norms against reference solutions), observed convergence rates under mesh refinement, direct comparisons with known soliton or asymptotic solutions, and quantitative measures of how well the main asymptotic features are reproduced in long-time runs. The abstract will be updated to reflect these additions. revision: yes

  2. Referee: [Deift-Zhou analysis] Description of the Deift-Zhou analysis: the extension to the 3x3 matrix RH problem is asserted to proceed by identifying two symmetric stationary points and dividing the (x,t)-plane into three regions with corresponding contour deformations, yet the manuscript does not explicitly verify that the more involved jump-matrix structure permits the standard g-function construction, lens opening, and decay estimates without new obstructions. This assumption underpins both the contour choice and the claim that asymptotics are captured.

    Authors: The phase function is scalar and identical to the scalar mKdV case, so the stationary points and the division into three regions follow directly. The 3x3 jump matrices admit the same factorization and g-function construction because their off-diagonal entries are of the form that permit exponential decay away from the stationary points once the lenses are opened; no new singularities or obstructions arise. Nevertheless, we acknowledge that an explicit verification is not currently provided. In the revised manuscript we will insert a short subsection (or expanded paragraph) that sketches the g-function, confirms the lens openings remain valid for the 3x3 structure, and outlines the decay estimates to demonstrate the absence of new obstructions. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation applies standard NIST components to 3x3 RH problem without self-referential reduction.

full rationale

The paper's chain consists of (1) computing scattering data via Chebyshev collocation on the direct problem and (2) applying Deift-Zhou steepest descent to deform the 3x3 oscillatory RH problem, with the phase function's two symmetric stationary points and three-region division stated as direct consequences of the coupled mKdV dispersion relation. No parameters are fitted to a data subset and then renamed as predictions of related quantities; the stationary-point count and contour choices are not defined in terms of the final numerical output; and no load-bearing self-citation or ansatz smuggling is described. The direct-at-(x,t) evaluation is the defining property of the IST method itself, not a derived claim that collapses to the input data by construction. The derivation therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on the assumption that the coupled mKdV equation admits an integrable 3x3 Riemann-Hilbert formulation whose scattering data can be computed numerically and whose inverse problem can be deformed via the Deift-Zhou method; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption The coupled mKdV equation is associated with a 3x3 matrix Riemann-Hilbert problem whose jump matrix and scattering data admit the stated structure.
    Invoked in the opening paragraph as the basis for extending NIST.
  • domain assumption The phase function of the oscillatory RH problem possesses exactly two stationary points symmetric about the origin, permitting a three-region division of the (x,t)-plane.
    Stated as the key geometric fact enabling the contour deformations.

pith-pipeline@v0.9.0 · 5525 in / 1508 out tokens · 29417 ms · 2026-05-07T09:46:26.515222+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references · 53 canonical work pages

  1. [1]

    Athorne, C., Fordy, A.: Generalised KdV and MKdV equations associated with symmetric spaces. J. Phys. A: Math. Gen.20(6), 1377–1386 (1987)

    work page 1987

  2. [2]

    Ablowitz, M.J., Kaup, D.J., Newell, A.C., Segur, H.: The Inverse Scattering Transform-Fourier Analysis for Nonlinear Problems. Stud. Appl. Math.53(4), 249–315 (1974)

    work page 1974

  3. [3]

    SIAM, Philadelphia (1981) Aky¨ uz, A., Sezer, M.: A Chebyshev collocation method for the solution of linear integro-differential equations

    Ablowitz, M.J., Segur, H.: Solitons and the Inverse Scattering Transform. SIAM, Philadelphia (1981) Aky¨ uz, A., Sezer, M.: A Chebyshev collocation method for the solution of linear integro-differential equations. Int. J. Comput. Math.72(4), 491–507 (1999)

    work page 1981

  4. [4]

    Bilman, D., Trogdon, T.: Numerical Inverse Scattering for the Toda Lattice. Commun. Math. Phys.352, 805–879 (2017)

    work page 2017

  5. [5]

    Nonlinearity33, 2211–2269 (2020)

    Bilman, D., Trogdon, T.: On numerical inverse scattering for the Korteweg-de Vries equation with discontinuous step-like data. Nonlinearity33, 2211–2269 (2020)

    work page 2020

  6. [6]

    Buckingham, R., Venakides, S.: Long-time asymptotics of the nonlinear Schr¨ odinger equation shock problem. Commun. Pure Appl. Math.60, 1349–1414 (2007) 31

    work page 2007

  7. [7]

    Cheng, P.J., Venakides, S., Zhou, X.: Long-time asymptotics for the pure radiation solution of the sine-Gordon equation. Comm. Partial Differential Equations.24, 1195–1262 (1999)

    work page 1999

  8. [8]

    Cui, S., Wang, Z.: Numerical inverse scattering transform for the focusing and defocusing Kundu-Eckhaus equations. Phys. D454, 133838 (2023)

    work page 2023

  9. [9]

    Cui, S., Wang, Z.: Efficient method to calculate the eigenvalues of the Zakharov-Shabat system. Chin. Phys. B33, 010201 (2024)

    work page 2024

  10. [10]

    Nonlinearity37, 105015 (2024) De Monvel, A.B., Kostenko, A., Shepelsky, D., Teschl, G.: Long-time asymptotics for the Camassa-Holm equation

    Cui, S., Wang, Z.: Numerical inverse scattering transform for the derivative nonlinear Schr¨ odinger equation. Nonlinearity37, 105015 (2024) De Monvel, A.B., Kostenko, A., Shepelsky, D., Teschl, G.: Long-time asymptotics for the Camassa-Holm equation. SIAM J. Math. Anal.41, 1559–1588 (2009)

    work page 2024

  11. [11]

    Deconinck, B., Trogdon, T., Yang, X.: Numerical inverse scattering for the sine-Gordon equation. Phys. D399, 159–172 (2019)

    work page 2019

  12. [12]

    Asymptotics for the MKdV equation

    Deift, P., Zhou, X.: A steepest descent method for oscillatory Riemann-Hilbert problems. Asymptotics for the MKdV equation. Ann. of Math.137(2), 295–368 (1993)

    work page 1993

  13. [13]

    Higher order theory

    Deift, P., Zhou, X.: Long-time asymptotics for integrable systems. Higher order theory. Commun. Math. Phys.165, 175–191 (1994)

    work page 1994

  14. [14]

    Fan, E.: Soliton solutions for a generalized Hirota-Satsuma coupled KdV equation and a coupled MKdV equation. Phys. Lett. A282(1–2), 18–22 (2001)

    work page 2001

  15. [15]

    Geng, X., Chen, M., Wang, K.: Long-time asymptotics of the coupled modified Korteweg-de Vries equation. J. Geom. Phys.142, 151–167 (2019)

    work page 2019

  16. [16]

    Gardner, C.S., Greene, J.M., Kruskal, M.D., Miura, R.M.: Method for solving Korteweg-de Vries equation. Phys. Rev. Lett.19, 1095–1097 (1967)

    work page 1967

  17. [17]

    Nonlinear Process

    Grimshaw, R., Pelinovsky, E., Talipova, T.: The modified Korteweg-de Vries equation in the theory of large-amplitude internal waves. Nonlinear Process. Geophys.4(4), 237–250 (1997)

    work page 1997

  18. [18]

    Journal of Nonlinear Waves2, 1–55 (2026)

    Gkogkou, A., Prinari, B., Trogdon, T.: Numerical inverse scattering transform for the defocusing nonlinear Schr¨ odinger equation with box-type initial conditions on a nonzero background. Journal of Nonlinear Waves2, 1–55 (2026)

    work page 2026

  19. [19]

    Grunert, K., Teschl, G.: Long-Time Asymptotics for the Korteweg-de Vries Equation via Nonlinear Steepest Descent. Math. Phys. Anal. Geom.12, 287–324 (2009)

    work page 2009

  20. [20]

    Geng, X., Zhai, Y., Dai, H.H.: Algebro-geometric solutions of the coupled modified Korteweg-de Vries hierarchy. Adv. Math.263, 123–153 (2014) 32

    work page 2014

  21. [21]

    Chaos Solitons Fractals13(9), 1917–1929 (2002)

    Helal, M.A.: Soliton solution of some nonlinear partial differential equations and its applications in fluid mechanics. Chaos Solitons Fractals13(9), 1917–1929 (2002)

    work page 1917

  22. [22]

    molecule solutions

    Hirota, R.: “molecule solutions” of coupled modified KdV equations. J. Phys. Soc. Jpn.66(9), 2530–2532 (1997)

    work page 1997

  23. [23]

    Hu, H.C.: Analytical positon, negaton and complexiton solutions for the coupled modified KdV system. J. Phys. A: Math. Theor.42(18), 185207 (2009)

    work page 2009

  24. [24]

    Iwao, M., Hirota, R.: Soliton solutions of a coupled modified KdV equations. J. Phys. Soc. Jpn.66(3), 577–588 (1997)

    work page 1997

  25. [25]

    Kamvissis, S.: Long time behavior for the focusing nonlinear Schr¨ odinger equation with real spectral singularities. Commun. Math. Phys.180, 325–341 (1996)

    work page 1996

  26. [26]

    Khater, A.H., El-Kalaawy, O.H., Callebaut, D.K.: B¨ acklund transformations and exact solutions for Alfv´ en solitons in a relativistic electron-positron plasma. Phys. Scr. 58(6), 545–548 (1998)

    work page 1998

  27. [27]

    Lin, S., Chen, Y.: Physics-informed neural network methods based on Miura trans- formations and discovery of new localized wave solutions. Phys. D445, 133629 (2023)

    work page 2023

  28. [28]

    Leblond, H., Mihalache, D.: Ultrashort spatiotemporal optical solitons in waveguide arrays: the effect of combined linear and nonlinear couplings. J. Phys. A: Math. Theor.51(43), 435202 (2018)

    work page 2018

  29. [29]

    Preprint at https://arxiv.org/abs/2510.12129 (2025)

    Ling, L., Su, H.: Nonlinear stability of vector multi-solitons in coupled NLS and modified KdV equations. Preprint at https://arxiv.org/abs/2510.12129 (2025)

    work page arXiv 2025

  30. [30]

    Ma, W.X.: Riemann-Hilbert problems and N-soliton solutions for a coupled mKdV system. J. Geom. Phys.132, 45–54 (2018)

    work page 2018

  31. [31]

    Ma, W.X.: The inverse scattering transform and soliton solutions of a combined modified Korteweg-de Vries equation. J. Math. Anal. Appl.471(1–2), 796–811 (2019)

    work page 2019

  32. [32]

    Mathematics7(7), 573 (2019)

    Ma, W.X.: Long-time asymptotics of a three-component coupled mKdV system. Mathematics7(7), 573 (2019)

    work page 2019

  33. [33]

    Miura, R.M.: Korteweg-de Vries equation and generalizations. I. A remarkable explicit nonlinear transformation. J. Math. Phys.9(8), 1202–1204 (1968) Mel’nikov, I.V., Mihalache, D., Moldoveanu, F., Panoiu, N.C.: Quasiadiabatic follow- ing of femtosecond optical pulses in a weakly excited semiconductor. Phys. Rev. A 56(2), 1569 (1997)

    work page 1968

  34. [34]

    Consul-tants Bureau, New York (1984)

    Novikov, S.P., Manakov, S.V., Pitaevskii, L.P., Zakharov, V.E.: Theory of Solitons: 33 The Inverse Scattering Method. Consul-tants Bureau, New York (1984)

    work page 1984

  35. [35]

    Olver, S.: Numerical Solution of Riemann-Hilbert Problems: Painlev´ e II. Found. Comput. Math.11, 153–179 (2011)

    work page 2011

  36. [36]

    Olver, S.: A general framework for solving Riemann-Hilbert problems numerically. Numer. Math.122, 305–340 (2012)

    work page 2012

  37. [37]

    Ono, H.: Soliton fission in anharmonic lattices with reflectionless inhomogeneity. J. Phys. Soc. Jpn.61(12), 4336–4343 (1992)

    work page 1992

  38. [38]

    Trogdon, T., Deconinck, B.: A Riemann-Hilbert problem for the finite-genus solutions of the KdV equation and its numerical solution. Phys. D251, 1–18 (2013)

    work page 2013

  39. [39]

    Tian, S.F.: Initial-boundary value problems for the coupled modified Korteweg-de Vries equation on the interval. Commun. Pure Appl. Anal17(3), 923–957 (2018)

    work page 2018

  40. [40]

    Trogdon, T., Olver, S.: Numerical inverse scattering for the focusing and defocusing nonlinear Schr¨ odinger equations. Proc. A469, 20120330 (2013)

    work page 2013

  41. [41]

    SIAM, Philadelphia (2016)

    Trogdon, T., Olver, S.: Riemann-Hilbert Problems, Their Numerical Solution, and the Computation of Nonlinear Special Functions. SIAM, Philadelphia (2016)

    work page 2016

  42. [42]

    Trogdon, T., Olver, S., Deconinck, B.: Numerical inverse scattering for the Korteweg- de Vries and modified Korteweg-de Vries equations. Phys. D241, 1003–1025 (2012)

    work page 2012

  43. [43]

    Tsuchida, T., Wadati, M.: The coupled modified Korteweg-de Vries equations. J. Phys. Soc. Jpn.67(4), 1175–1187 (1998)

    work page 1998

  44. [44]

    Wadati, M.: The modified Korteweg-de Vries equation. J. Phys. Soc. Jpn.34(5), 1289–1296 (1973)

    work page 1973

  45. [45]

    Wu, J., Geng, X.: Inverse scattering transform and soliton classification of the coupled modified Korteweg-de Vries equation. Commun. Nonlinear Sci. Numer. Simul.53, 83–93 (2017)

    work page 2017

  46. [46]

    Nonlinear Anal

    Wang, D.S., Wang, X.L.: Long-time asymptotics and the bright N-soliton solutions of the Kundu-Eckhaus equation via the Riemann-Hilbert approach. Nonlinear Anal. Real World Appl.41, 334–361 (2018)

    work page 2018

  47. [47]

    Wang, D.S., Xu, L., Xuan, Z.: The complete classification of solutions to the Riemann problem of the defocusing complex modified KdV equation. J. Nonlinear Sci.32(1), 3 (2022)

    work page 2022

  48. [48]

    Xu, J., Fan, E., Chen, Y.: Long-time Asymptotic for the Derivative Nonlinear Schr¨ odinger Equation with Step-like Initial Value. Math. Phys. Anal. Geom.16, 253–288 (2013) 34

    work page 2013

  49. [49]

    Xue, B., Li, F., Yang, G.: Explicit solutions and conservation laws of the coupled modified Korteweg-de Vries equation. Phys. Scr.90(8), 085204 (2015)

    work page 2015

  50. [50]

    Xu, S.: A coupled complex mKdV equation and its N-soliton solutions via the Riemann-Hilbert approach. Bound. Value Probl.2023(1), 83 (2023)

    work page 2023

  51. [51]

    SIAM, Philadelphia (2010)

    Yang, J.K.: Nonlinear Waves in Integrable and Nonintegrable Systems. SIAM, Philadelphia (2010)

    work page 2010

  52. [52]

    Yajima, N., Oikawa, M.: A class of exactly solvable nonlinear evolution equations. Prog. Theor. Phys.54(5), 1576–1577 (1975)

    work page 1975

  53. [53]

    Chaos Solitons Fractals201, 117185 (2025) 35

    Zhang, W.X., Chen, Y.: Numerical inverse scattering transform for the coupled nonlinear Schr¨ odinger equation. Chaos Solitons Fractals201, 117185 (2025) 35

    work page 2025