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

Recognition: unknown

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

Wen-Xin Zhang , Yong Chen

Authors on Pith no claims yet

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 comments

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.

Desk Editor's Note private letter to a colleague

The paper gives a concrete NIST extension to the 3x3 coupled mKdV with working contour deformations, but the numerical validation stays thin on errors and comparisons.

read the letter

The paper takes the numerical inverse scattering transform and applies it to the coupled mKdV system, which produces a 3x3 matrix Riemann-Hilbert problem. That step is the main new piece: they compute scattering data via Chebyshev collocation on the direct problem, then deform the oscillatory RH problem using Deift-Zhou steepest descent. The phase has two symmetric stationary points, so they split the (x,t) plane into three regions and open lenses accordingly. This lets them evaluate the solution at fixed points without time marching, which is the practical payoff for long-time runs. The abstract indicates the numerics recover the main asymptotic features, and the setup looks like a direct, non-routine extension rather than a routine scalar copy-paste. If the full text shows the matrix jump conditions did not force extra g-function adjustments or break the decay estimates, that counts as solid technical work for people who already know the scalar NIST literature. The stress-test worry about extra phase contributions from the 3x3 structure does not seem to have blocked the three-region division, so the basic claim holds up on the description given. The soft spot is the evidence. The paper mentions experiments that illustrate performance, yet supplies no error tables, convergence plots, or side-by-side checks against known solutions or other solvers. Without those numbers it is hard to judge how well the method actually captures the asymptotics or how sensitive the contour choices are to the matrix entries. Reproducibility is also unclear since no code or data files are referenced. This work is aimed at the integrable-systems subgroup that builds numerical tools for nonlinear waves. A reader already comfortable with RH problems and steepest-descent analysis will pick up the specific 3x3 contour recipe and the region division quickly. It is not broad enough to interest people outside that niche, but inside it the contribution is focused and usable. The paper deserves a serious referee because the core construction is grounded in standard methods, the extension is presented as new, and the direct-evaluation advantage is real even if the error analysis needs tightening. I would send it out for review and ask the authors to add quantitative error metrics and at least a reproducibility statement.

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)

[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.
[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)

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

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
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)

Reference graph

Works this paper leans on

53 extracted references · 1 canonical work pages

[1]

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

1987
[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)

1974
[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)

1981
[4]

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

2017
[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)

2020
[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

2007
[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)

1999
[8]

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

2023
[9]

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

2024
[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)

2024
[11]

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

2019
[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)

1993
[13]

Higher order theory

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

1994
[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)

2001
[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)

2019
[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)

1967
[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)

1997
[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)

2026
[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)

2009
[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

2014
[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)

1917
[22]

molecule solutions

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

1997
[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)

2009
[24]

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

1997
[25]

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

1996
[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)

1998
[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)

2023
[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)

2018
[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]

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

2018
[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)

2019
[32]

Mathematics7(7), 573 (2019)

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

2019
[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)

1968
[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)

1984
[35]

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

2011
[36]

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

2012
[37]

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

1992
[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)

2013
[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)

2018
[40]

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

2013
[41]

SIAM, Philadelphia (2016)

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

2016
[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)

2012
[43]

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

1998
[44]

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

1973
[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)

2017
[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)

2018
[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)

2022
[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

2013
[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)

2015
[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)

2023
[51]

SIAM, Philadelphia (2010)

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

2010
[52]

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

1975
[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

2025