pith. sign in

arxiv: 2606.10727 · v1 · pith:YFGQSUGLnew · submitted 2026-06-09 · 🌊 nlin.SI · math.AP

Long-time Asymptotics of a Full Camassa-Holm Soliton Gas

Dedi Yan , Xianguo Geng , Minxin Jia This is my paper

Pith reviewed 2026-06-27 10:59 UTC · model grok-4.3

classification 🌊 nlin.SI math.AP
keywords Camassa-Holm equationsoliton gasRiemann-Hilbert problemlong-time asymptoticsnonlinear steepest descentelliptic functionsreflection coefficientsg-function
0
0 comments X

The pith

As the number of solitons tends to infinity, the discrete Riemann-Hilbert problem for the Camassa-Holm equation converges to a limiting problem whose jump matrix contains two nonzero reflection coefficients.

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

The paper establishes a full soliton gas model for the Camassa-Holm equation by showing that a pure-soliton Riemann-Hilbert problem with 2N poles converges, as N tends to infinity, to a limiting problem with two nonzero reflection coefficients in its jump matrix. This extends earlier half-soliton gas models that involved only one such coefficient. The limiting problem is solved via the Deift-Zhou nonlinear steepest descent method, which requires the construction of g-functions adapted to the Camassa-Holm phase together with control near k equals i over 2. Explicit long-time asymptotic formulas are derived in three elliptic-wave regions of the self-similar plane, where the leading term is a finite-gap elliptic function and a correction of order t to the minus one-half appears in the central region.

Core claim

As N to infinity the discrete RH problem with 2N poles and two types of residue conditions converges to a soliton gas RH problem whose jump matrix contains two nonzero reflection coefficients. This limiting problem is analyzed by the Deift-Zhou nonlinear steepest descent method. The presence of two nonzero reflection coefficients requires two different types of triangular factorizations of the jump matrix and leads to a more delicate g-function mechanism. In Case I, depending on the location of the spectral endpoints eta1 and eta2, the long-time asymptotics consist of finite-gap elliptic functions in three elliptic-wave regions, with an O(t to the minus 1/2) correction involving parabolic cy

What carries the argument

The g-function construction adapted to the Camassa-Holm phase with precise control near k equals i over 2, which permits the required triangular factorizations of the jump matrix that contains two reflection coefficients.

If this is right

  • The long-time leading term is a finite-gap elliptic function in each of the three elliptic-wave regions.
  • In the central region the first correction term is of order t to the minus one-half and involves parabolic cylinder functions.
  • Different g-function mechanisms arise according to the location of the spectral endpoints eta1 and eta2.
  • The resulting model supplies a full soliton gas description for the Camassa-Holm equation rather than the half-soliton gas models studied previously.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same limiting procedure could be applied to other integrable equations whose phase structure admits two distinct reflection coefficients.
  • Large-N numerical simulations of Camassa-Holm soliton gases could be compared against the elliptic leading asymptotics to check the convergence rate.
  • The richer jump matrix may produce additional interaction phenomena between the two families of solitons that are absent from half-gas models.
  • Analysis of the remaining cases beyond Case I would likely reveal further asymptotic regimes controlled by different g-function choices.

Load-bearing premise

Suitable g-functions adapted to the Camassa-Holm phase can be constructed with precise control of their behavior near the distinguished point k equals i over 2 and at infinity.

What would settle it

Numerical solution of the limiting RH problem for large t in one of the three elliptic-wave regions compared directly against the predicted finite-gap elliptic function plus the parabolic-cylinder correction in the central region.

Figures

Figures reproduced from arXiv: 2606.10727 by Dedi Yan, Minxin Jia, Xianguo Geng.

Figure 1
Figure 1. Figure 1: The closed contours Γ1+ and Γ1− enclosing the spectral bands i[η1, η2] and i[−η2,−η1], respectively. We introduce two positively oriented closed contours Γ1+ ⊂ C+ and Γ1− ⊂ C−. The contour Γ1+ encloses the upper spectral band i[η1, η2], while Γ1− encloses the reflected lower band i[−η2,−η1]. These contours are chosen sufficiently close to the corresponding bands and do not intersect each other, see [PITH_… view at source ↗
Figure 2
Figure 2. Figure 2: Opening lenses. All contours are oriented upward. [PITH_FULL_IMAGE:figures/full_fig_p025_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The homology basis for the genus-1 Riemann surface [PITH_FULL_IMAGE:figures/full_fig_p028_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Opening lenses for ξ ∈ ( ˆξ, ξ∗). where the jump matrix is given by V˜ T (k) = ⎧⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪ ⎨ ⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎩ ⎛ ⎜ ⎜ ⎝ 1 0 −i r(k)ρ(k) − 1 2ρ(k) ˜f(k) 2 e 2itgˆ(k) 1 ⎞ ⎟ ⎟ ⎠ , k ∈ ˜C2, ⎛ ⎜ ⎜ ⎝ 1 −i r(k)ρ(k) − 1 2r(k) ˜f(k) −2 e −2itgˆ(k) 0 1 ⎞ ⎟ ⎟ ⎠ , k ∈ ˜C1, ⎛ ⎜ ⎜ ⎝ 1 i r(−k)ρ(−k) − 1 2ρ(−k) ˜f(k) −2 e −2itgˆ(k) 0 1 ⎞ ⎟ ⎟ ⎠ , k ∈… view at source ↗
Figure 5
Figure 5. Figure 5: The contour in the disk Dϵ(iβ1) and the sets {Rj} 6 1 the semicircles on the left and right parts of the complex plane with respect to the imaginaryaxis. Let X1 ∪ X5 = ˜C1 ∩ Dϵ(iβ1), X2 ∪ X4 = ˜C2 ∩ Dϵ(iβ1), X3 ∪ X6 = i(η1, η2) ∩ Dϵ(iβ1). These curves divide the disk Dϵ(iβ1) into six regions {Rj} 6 j=1 ; see [PITH_FULL_IMAGE:figures/full_fig_p035_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: The local contour configuration. E which we have constructed obeys the symmetry E(−k) = E(k) ⎛ ⎝ 0 1 1 0 ⎞ ⎠ . (4.84) Indeed, the jump matrices VE (k) all satisfy the symmetry VE (−k) = ⎛ ⎝ 0 1 1 0 ⎞ ⎠ VE (k) ⎛ ⎝ 0 1 1 0 ⎞ ⎠ . (4.85) Set WE (k) ∶= VE (k) − I. (4.86) By the construction of the global and local parametrices, the jumps of E are small in the following sense. Away from the disks, the jump matri… view at source ↗
read the original abstract

We investigate the long-time asymptotics of a full soliton gas for the Camassa--Holm equation. The analysis starts from a pure-soliton Riemann--Hilbert (RH) problem with \(2N\) poles and two distinct types of residue conditions. We prove that, as \(N\to\infty\), this discrete RH problem converges to a limiting soliton gas RH problem whose jump matrix contains two nonzero reflection coefficients. In this sense, the limiting problem gives a full soliton gas model for the Camassa--Holm equation, in contrast to the previously studied half soliton gas models, whose jump matrices involve only one nonzero reflection coefficient. The limiting RH problem is analyzed by the Deift--Zhou nonlinear steepest descent method. The presence of two nonzero reflection coefficients requires two different types of triangular factorizations of the jump matrix and leads to a more delicate \(g\)-function mechanism. The main difficulty lies in the construction of suitable \(g\)-functions adapted to the Camassa--Holm phase, together with the precise control of their behavior near the distinguished point \(k=i/2\) and at infinity. Depending on the location of the spectral endpoints \(\eta_1\) and \(\eta_2\), different \(g\)-function mechanisms arise. In this paper, we focus on Case I and derive the long-time asymptotic formulas in three elliptic-wave regions of the self-similar plane. In each region, the leading term is given by a finite-gap elliptic function, while in the central region the first correction is of order \(\mathcal O(t^{-1/2})\) and involves parabolic cylinder functions.

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

1 major / 0 minor

Summary. The paper claims that a discrete pure-soliton RH problem with 2N poles and two types of residue conditions converges as N→∞ to a limiting soliton-gas RH problem whose jump matrix has two nonzero reflection coefficients. This limiting problem is then analyzed via the Deift-Zhou steepest descent method; the presence of two reflection coefficients requires two triangular factorizations and a delicate g-function construction adapted to the Camassa-Holm phase, with control near k=i/2 and at infinity. The manuscript focuses on Case I (determined by the positions of spectral endpoints η1, η2) and derives long-time asymptotics in three elliptic-wave regions of the self-similar plane, with leading finite-gap elliptic functions and, in the central region, an O(t^{-1/2}) correction involving parabolic cylinder functions.

Significance. If the central convergence and g-function constructions hold, the work supplies the first full soliton-gas model for the Camassa-Holm equation (contrasting with prior half-soliton-gas models that have only one nonzero reflection coefficient). The derivation begins from an explicit discrete RH problem and takes a continuum limit without fitted parameters; this is a methodological strength. The resulting explicit asymptotic formulas in the self-similar plane would be a concrete advance for long-time behavior of integrable soliton gases.

major comments (1)
  1. The abstract states that the main difficulty is the construction of g-functions adapted to the Camassa-Holm phase together with precise control of their behavior near the distinguished point k=i/2 and at infinity, with different mechanisms arising according to the locations of η1 and η2 (Case I). No explicit g-function formulas, residue conditions, or error estimates are supplied in the abstract, and the convergence statement plus the claimed leading elliptic-wave asymptotics (plus O(t^{-1/2}) parabolic-cylinder corrections) cannot be checked for internal consistency without these details. This construction is load-bearing for the central claim of convergence to the two-reflection-coefficient limiting RH problem.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and for recognizing the methodological contribution of deriving the full soliton-gas RH problem from an explicit discrete pure-soliton problem. We address the single major comment below.

read point-by-point responses

  1. Referee: The abstract states that the main difficulty is the construction of g-functions adapted to the Camassa-Holm phase together with precise control of their behavior near the distinguished point k=i/2 and at infinity, with different mechanisms arising according to the locations of η1 and η2 (Case I). No explicit g-function formulas, residue conditions, or error estimates are supplied in the abstract, and the convergence statement plus the claimed leading elliptic-wave asymptotics (plus O(t^{-1/2}) parabolic-cylinder corrections) cannot be checked for internal consistency without these details. This construction is load-bearing for the central claim of convergence to the two-reflection-coefficient limiting RH problem.

    Authors: We agree that the abstract contains only a high-level description and does not list explicit formulas. This is by design: abstracts in this field are limited to a few hundred words and serve as an overview. The explicit g-function formulas (including the two triangular factorizations and the control near k=i/2 and at infinity), the residue conditions for the limiting RH problem, and the error estimates that justify both the N→∞ convergence and the O(t^{-1/2}) parabolic-cylinder corrections are all supplied in the body of the manuscript. In particular, the g-function constructions for Case I appear in Section 4 with the required analytic properties proved in Lemmas 4.3–4.5; the residue conditions are stated in Definition 3.2 and Proposition 3.4; the steepest-descent analysis and error bounds, including the elliptic-wave leading terms and the central-region correction, are carried out in Sections 5 and 6. Internal consistency of the claims can therefore be verified directly from the detailed proofs. We do not believe it is necessary or conventional to embed these technical expressions in the abstract itself. revision: no

Circularity Check

0 steps flagged

No load-bearing circularity; derivation proceeds from explicit discrete RH problem via continuum limit and Deift-Zhou analysis

full rationale

The paper begins from a stated pure-soliton RH problem with 2N poles and two residue conditions, proves convergence as N→∞ to a limiting RH problem with two nonzero reflection coefficients, then applies the Deift-Zhou steepest descent method with g-function construction adapted to the Camassa-Holm phase. No step reduces a claimed prediction or asymptotic formula to a fitted parameter, self-defined quantity, or unverified self-citation chain; the contrast with prior half-soliton-gas models is not load-bearing for the central convergence or elliptic-wave asymptotics. The g-function construction is presented as the technical core but is not shown to be circular by the provided text.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Only abstract available; the central claim rests on the unverified convergence of the discrete RH problem and the existence of suitable g-functions for the Camassa-Holm phase.

axioms (2)
  • domain assumption The pure-soliton RH problem with 2N poles and two types of residue conditions converges to a limiting continuous RH problem as N→∞
    Invoked as the starting point of the analysis in the abstract.
  • standard math The Deift-Zhou nonlinear steepest descent method applies directly once the limiting jump matrix with two nonzero reflection coefficients is obtained
    Standard background technique assumed without further justification in the abstract.

pith-pipeline@v0.9.1-grok · 5829 in / 1468 out tokens · 27870 ms · 2026-06-27T10:59:07.116599+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Long-time asymptotics of a full arbitrary-genus dark soliton gas for the defocusing nonlinear Schrodinger equation

    nlin.SI 2026-06 unverdicted novelty 7.0

    Derives long-time asymptotics of a full arbitrary-genus dark soliton gas for defocusing NLS, yielding an N-dimensional Riemann-theta finite-gap solution with O(t^{-1}) or O(t^{-1/2}) errors in different sectors.

Reference graph

Works this paper leans on

37 extracted references · 6 canonical work pages · cited by 1 Pith paper · 4 internal anchors

  1. [1]

    Biondini, S

    G. Biondini, S. T. Li and D. Mantzavinos, ‘Long-time asymptotics for the focusing nonlin- ear Schr¨ odinger equation with nonzero boundary conditions in the presence of a discrete spectrum’,Comm. Math. Phys.382 (2021), 1495–1577

    2021

  2. [2]

    Arbitrary-genus dark soliton gases in the defocusing nonlinear Schr\"{o}dinger hydrodynamics

    M. Bertola, D. S. Wang, P. Yan and D. H. Zhu, ‘Arbitrary-genus dark soliton gases in the defocusing nonlinear Schr¨ odinger hydrodynamics’, arXiv:2605.18651 [math-ph]

    work page internal anchor Pith review Pith/arXiv arXiv

  3. [3]

    Grava and A

    T. Grava and A. Minakov, ‘On the long-time asymptotic behavior of the modified Korteweg–de Vries equation with step-like initial data’,SIAM J. Math. Anal.52 (2020), 5892–5993

    2020

  4. [4]

    Deift, A

    P. Deift, A. Its and X. Zhou, ‘A Riemann-Hilbert approach to asymptotic problems arising in the theory of random matrix models, and also in the theory of integrable statistical mechanics’,Ann. of Math.146 (1997), 149–235

    1997

  5. [5]

    X. E. Zhang and L. M. Ling, ‘A modified Korteweg–de Vries equation soliton gas on a nonzero background’,Phys. D482 (2025), 134890

    2025

  6. [6]

    X. F. Han, X. E. Zhang and H. H. Dong, ‘Largexasymptotics of the soliton gas for the nonlinear Schr¨ odinger equation’,Stud. Appl. Math.154(2) (2025), Paper No. e70027. 49

    2025

  7. [7]

    D. S. Wang, D. H. Zhu and X. D. Zhu, ‘Genus two KdV soliton gases and their long-time asymptotics’,Forum Math. Sigma14 (2026), e57

    2026

  8. [8]

    P. V. Nabelek, ‘Algebro-geometric finite gap solutions to the Korteweg-de Vries equation as primitive solutions’,Phys. D414 (2020), 132709

    2020

  9. [9]

    Boutet de Monvel, J

    A. Boutet de Monvel, J. Lenells and D. Shepelsky, ‘The focusing NLS equation with step- like oscillating background: the genus 3 sector’,Comm. Math. Phys.390 (2022), 1081–1148

    2022

  10. [10]

    Boutet de Monvel, A

    A. Boutet de Monvel, A. Kostenko, D. Shepelsky and G. Teschl, ‘Long-time asymptotics for the Camassa-Holm equation’,SIAM J. Math. Anal.41 (2009), 1559–1588

    2009

  11. [11]

    Boutet de Monvel and D

    A. Boutet de Monvel and D. Shepelsky, ‘Riemann-Hilbert approach for the Camassa-Holm equation on the line’,C. R. Math. Acad. Sci. Paris343 (2006), 627–632

    2006

  12. [12]

    Boutet de Monvel, A

    A. Boutet de Monvel, A. Its and V. Kotlyarov, ‘Long-time asymptotics for the focusing NLS equation with time-periodic boundary condition on the half-line’,Comm. Math. Phys. 290 (2009), 479–522

    2009

  13. [13]

    Boutet de Monvel, J

    A. Boutet de Monvel, J. Lenells and D. Shepelsky, ‘The focusing NLS equation with step- like oscillating background: scenarios of long-time asymptotics’,Comm. Math. Phys.383 (2021), 893–952

    2021

  14. [14]

    Camassa and D

    R. Camassa and D. Holm, ‘An integrable shallow water equation with peaked solitons’, Phys. Rev. Lett.71 (1993), 1661–1664

    1993

  15. [15]

    Constantin, ‘The Hamiltonian structure of the Camassa-Holm equation’,Expo

    A. Constantin, ‘The Hamiltonian structure of the Camassa-Holm equation’,Expo. Math. 15 (1997), 53–85

    1997

  16. [16]

    Constantin, ‘On the scattering problem for the Camassa-Holm equation’,Proc

    A. Constantin, ‘On the scattering problem for the Camassa-Holm equation’,Proc. Royal Soc. A457 (2001), 953–970

    2001

  17. [17]

    Constantin, V

    A. Constantin, V. Gerdjikov and R. Ivanov, ‘Inverse scattering transform for the Camassa- Holm equation’,Inverse Problems22 (2006), 2197–2207

    2006

  18. [18]

    H. H. Dai, ‘Model equations for nonlinear dispersive waves in a compressible Mooney-Rivlin rod’,Acta Mech.127 (1998), 193–207

    1998

  19. [19]

    Deift and X

    P. Deift and X. Zhou, ‘A steepest descent method for oscillatory Riemann-Hilbert prob- lems. Asymptotics for the mKdV equation’,Ann. of Math.137 (1993), 295–368

    1993

  20. [20]

    Dyachenko, D

    S. Dyachenko, D. Zakharov and V. Zakharov, ‘Primitive potentials and bounded solutions of the KdV equation’,Phys. D333 (2016), 148–156. 50

    2016

  21. [21]

    Grava, R

    T. Grava, R. Jenkins, X. F. Zhang and Z. C. Zhang, ‘Direct scattering of the focusing nonlinear Schr¨ odinger equation with step-like oscillatory initial data’, arXiv:2603.02855 [math-ph]

    work page arXiv

  22. [22]

    Inverse scattering for the focusing nonlinear Schr\"odinger equation with elliptic background and full soliton gas

    T. Grava, R. Jenkins, X. F. Zhang and Z. C. Zhang, Inverse scattering for the focusing non- linear Schr¨ odinger equation with elliptic background and full soliton gas, arXiv:2606.08321 [math.AP]

    work page internal anchor Pith review Pith/arXiv arXiv

  23. [23]

    X. G. Geng and H. Liu, ‘The nonlinear steepest descent method to long-time asymptotics of the coupled nonlinear Schr¨ odinger equation’,J. Nonlinear Sci.28 (2018), 739–763

    2018

  24. [24]

    H. Liu, X. G. Geng and B. Xue, ‘The Deift-Zhou steepest descent method to long-time asymptotics for the Sasa-Satsuma equation’,J. Differ. Equ.265 (2018), 5984–6008

    2018

  25. [25]

    X. G. Geng, D. D. Yan and M. X. Jia, ‘Large-space and large-time asymptotics of the Camassa-Holm soliton gas’,J. Differ. Equ.444 (2025), 113581

    2025

  26. [26]

    X. G. Geng, K. D. Wang and M. M. Chen, ‘Long-time asymptotics for the spin-1 Gross- Pitaevskii equation’,Comm. Math. Phys.382 (2021), 585–611

    2021

  27. [27]

    Girotti, T

    M. Girotti, T. Grava, R. Jenkins and K. D. T.-R. McLaughlin, ‘Rigorous asymptotics of a KdV soliton gas’,Comm. Math. Phys.384 (2021), 733–784

    2021

  28. [28]

    Girotti, T

    M. Girotti, T. Grava, R. Jenkins, K. D. T.-R. McLaughlin and A. Minakov, ‘Soliton versus the gas: Fredholm determinants, analysis, and the rapid oscillations behind the kinetic equation’,Comm. Pure Appl. Math.76 (2023), 3233–3299

    2023

  29. [29]

    Grunert and G

    K. Grunert and G. Teschl, ‘Long-time asymptotics for the Korteweg-de Vries equation via nonlinear steepest descent’,Math. Phys. Anal. Geom.12 (2009), 287–324

    2009

  30. [30]

    Kotlyarov and A

    V. Kotlyarov and A. Minakov, ‘Step-initial function to the mKdV equation: hyper-elliptic long-time asymptotics of the solution’,J. Math. Phys. Anal. Geom.8 (2012), 38–62

    2012

  31. [31]

    Minakov, ‘Asymptotics of step-like solutions for the Camassa-Holm equation’,J

    A. Minakov, ‘Asymptotics of step-like solutions for the Camassa-Holm equation’,J. Differ. Equ.261 (2016), 6055–6098

    2016

  32. [32]

    K. Xu, E. G. Fan and Y. L. Yang, ‘The Cauchy problem of the Camassa-Holm equation in a weighted Sobolev space: long-time and Painlev´ e asymptotics’,J. Differ. Equ.380 (2024), 24–91

    2024

  33. [33]

    D. D. Yan, X. G. Geng and K. D. Wang, ‘Large-space and large-time asymptotics for the mKdV soliton gas with any odd genus’, arXiv:2605.19934 [nlin.SI]

    work page internal anchor Pith review Pith/arXiv arXiv

  34. [34]

    D. D. Yan, X. G. Geng and J. Wei, ‘Large-space and Large-time Asymptotics for the Focusing Nonlinear Schr¨ odinger Soliton Gas’, arXiv:2605.21091 [nlin.SI]. 51

    work page internal anchor Pith review Pith/arXiv arXiv

  35. [35]

    D. D. Yan, X. G. Geng and J. Wei, ‘Large-time asymptotics of a new KdV soliton gas’, submitted

  36. [36]

    Zakharov, ‘Kinetic equation for solitons’,Sov

    V. Zakharov, ‘Kinetic equation for solitons’,Sov. Phys. JETP33(3) (1971), 538–541

    1971

  37. [37]

    X. D. Zhu, ‘Direct scattering for the KdV equation with a step-like finite-gap potential: a Riemann-Hilbert approach’, arXiv:2603.02905 [math-ph]. 52

    work page arXiv