REVIEW 2 major objections 5 minor 48 references
Asymptotic stability of solutions to the good Boussinesq equation in dispersive wave region
T0 review · 2 major / 5 minor · reviewed 2026-07-13 · grok-4.5
Pith's one-line read The good Boussinesq equation has asymptotically stable dispersive waves when reflection coefficients lie in a weighted Sobolev space, with leading term of order t^{-1/2} and error O(t^{-3/4}).
desk verdict Solid Dbar extension of good-Boussinesq asymptotics to H^{3,4} data, with the usual open RH-existence gap for the new class. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The ¯∂-steepest-descent method applied to the 3×3 Riemann-Hilbert problem associated with the Lax pair: successive contour deformations and continuous extensions of the reflection coefficients reduce the problem to a model problem near the three critical points plus two small-norm ¯∂-problems whose solutions contribute only O(t^{-3/4}) after reconstruction.
What would settle it
Exhibit a pair of reflection coefficients in H^{3,4} satisfying the paper’s assumptions for which either no RH solution exists or the reconstructed u(x,t) fails to obey the stated O(t^{-3/4}) bound for large t on some compact interval of positive x/t.
Extended reading notes
Core claim
For reflection coefficients r1 in H^{3,4}(0,∞) and r2 in H^{3,4}(-∞,0) satisfying the no-soliton and double-pole assumptions, the solution u(x,t) of the good Boussinesq equation satisfies |u(x,t) + (3^{5/4} k0 √ν / √(2t)) sin(phase involving ν ln(t k0^{2}), arg r1(k0), arg Γ(iν) and a logarithmic integral)| ≤ O(t^{-3/4}) uniformly for ζ = x/t in compact subsets of (0,∞), where k0 = ζ/2 and ν = - (1/(2π)) ln(1-|r1(k0)|^{2}).
Load-bearing premise
A solution to the underlying 3×3 Riemann-Hilbert problem is assumed to exist for the given weighted-Sobolev reflection coefficients; only uniqueness is proved.
Editorial extensions
If this is right
- The same leading-term formula and error rate hold, after the substitution r1 ↔ r2, for the left half-plane x < 0.
- Asymptotic stability persists for initial data whose Fourier-like scattering data decay only like weighted L^{2} rather than Schwartz.
- The error O(t^{-3/4}) is optimal for the ¯∂-method under the stated Sobolev regularity, so further improvement would require stronger assumptions on r1,r2.
- The reconstruction formulae for both u and w remain valid on open subsets of the half-plane once the RH solution exists.
Reading between the lines
- The same contour-opening and local-parametrix construction should extend without essential change to other 3×3 integrable systems whose phase functions possess three stationary points of the same type.
- If the double-pole condition at k=0 is relaxed, the expansion near the origin would acquire an additional residue contribution that could alter the constant term but not the t^{-1/2} decay.
- The weighted-Sobolev threshold H^{3,4} appears dictated by the need for four derivatives after reconstruction; a lower-regularity version might still yield a weaker but still decaying error.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper applies the ¯∂-steepest descent method to the 3×3 Riemann–Hilbert problem for the good Boussinesq equation and derives the long-time asymptotic expansion of the solution u(x,t) in the dispersive region ζ=x/t∈I⊂(0,∞). Under the assumptions that the reflection coefficients r1∈H^{3,4}(0,∞), r2∈H^{3,4}(-∞,0) satisfy Assumptions 2.1–2.2 (no solitons, double poles of s11 and sA11 at k=0), Theorem 2.1 states that |u+(3^{5/4}k0√ u/√(2t))sin(phase)| = O(t^{-3/4}), with the phase involving u, arg r1(k0), arg Γ(i u) and a principal-value integral of d ln(1-|r1|^{2}). The argument proceeds by three contour deformations (Sections 3.1–3.3), construction of local parametrices at the three critical points ω^j k0 via a model RH problem (Section 4), small-norm estimates for the residual RH problem, and L^∞ bounds of order t^{-1/4} (and solid-Cauchy contributions of order t^{-3/4}) for the two pure ¯∂ problems (Section 5). The reconstruction formula (9) then yields the claimed expansion. The main claimed advance relative to Charlier–Lenells–Wang [13] is the passage from Schwartz to weighted-Sobolev reflection coefficients and the resulting asymptotic stability statement.
Significance. If the existence gap is closed (or the result is carefully restated as conditional on RH existence), the work supplies a usable template for ¯∂ analysis of 3×3 Lax pairs and improves the admissible data class for the good Boussinesq equation. The error tracking through the solid Cauchy operators is careful and the O(t^{-3/4}) rate is optimal for the method. The explicit phase formula and the left-half-plane counterpart (Remark 2.4) are concrete and potentially useful for comparison with numerics or other asymptotic regimes. The manuscript therefore has clear technical value for the integrable-systems community working on higher-order spectral problems.
major comments (2)
- Remark 2.1 and Proposition 2.1 establish uniqueness of a solution to RH problem 2.1 only if a solution exists; existence for r1,r2∈H^{3,4} is left open (the vanishing lemma is unavailable under the good-Boussinesq symmetries, and the global well-posedness result of [48] is for the PDE, not for the RH map). Theorem 2.1 and the reconstruction (9) therefore rest on an unproved existence hypothesis. All subsequent deformations, local parametrices and error estimates inherit this conditionality. The claimed extension beyond Schwartz data is incomplete until existence is proved or the theorem is restated as conditional on the existence of a solution to RH problem 2.1.
- The justification for the precise weight H^{3,4} (Remark 2.1, Remark 2.3) is only sketched: the authors invoke integration by parts on Volterra equations and the need for fourth-order derivatives of the reflection coefficients when the asymptotic formula is substituted back into the PDE. A short, self-contained argument that H^{3,4} is both necessary and sufficient for the uniform convergence of the expansions (6)–(7) and for the L^{2} bounds used in the ¯∂ estimates would strengthen the paper; without it the choice of space remains somewhat ad hoc.
minor comments (5)
- Section 6 title: “Long-time saymptotics” → “asymptotics”.
- Abstract and Theorem 2.1: “the long-time asymptotic expansion of the solution are obtained” → subject–verb agreement.
- Fig. 1 and subsequent contour figures: the labels of the six sectors D_j and of the critical rays are dense; a clearer legend or a single composite figure with all intermediate contours would help the reader follow the three transformations.
- Notation for the continuous extensions G_j and R_j (Propositions 3.1, 3.3) is introduced without a uniform convention for the angular cut-offs; a short table or a sentence listing the support of each cut-off function would improve readability.
- References [12] and [13] are cited heavily for the original RH construction and the model problem; a one-paragraph summary of the jump matrices and of the model solution m^X would make the paper more self-contained.
Circularity Check
No circularity: long-time expansion is computed from explicit contour deformations, model parametrices and Dbar estimates on a pre-existing RH problem, not defined by or fitted to the claimed formula.
full rationale
The derivation chain begins from the 3 imes3 RH problem 2.1 (constructed in the independent reference [12]) whose jump matrices are built from the given reflection coefficients r1,r2. Three successive transformations (G(1), riangle, R(3)) open lenses and factor the problem into a pure RH problem plus two pure Dbar problems. Local parametrices at the three stationary points are taken from the model problem already solved in [13]; the resulting small-norm RH problem and the two solid Cauchy operators are estimated directly in L^p and L^ o L^ o norms, producing the O(t^{-3/4}) remainder. Reconstruction formula (9) then yields Theorem 2.1. None of these steps defines the target asymptotic expression in terms of itself, fits a free parameter to data that is later “predicted,” or rests on a uniqueness theorem proved only by the same authors for the same claim. The open existence question for the RH map under H^{3,4} data (Remark 2.1) is a completeness gap, not a circular reduction. Self-citations supply only the starting RH problem and the model parametrix; the long-time analysis performed here is independent of those inputs.
Assumptions & free parameters
assumptions (4)
- domain assumption A solution to RH problem 2.1 exists for r1∈H^{3,4}(0,∞), r2∈H^{3,4}(-∞,0) satisfying the stated limits at k=0.
- domain assumption s11(k) and sA11(k) are nonzero on the closed sectors D1 and D4 except at the origin, and both have double poles at k=0 (Assumptions 2.1–2.2).
- domain assumption The direct scattering map sends the initial data into reflection coefficients belonging to H^{3,4} and the reconstruction formula (9) recovers the PDE solution.
- standard math Standard L^p boundedness of Cauchy projectors on Carleson contours and the Beals-Coifman theory for small-norm RH problems.
Cite this review
Pith. "Pith review of Asymptotic stability of solutions to the good Boussinesq equation in dispersive wave region." pith.science (2026). https://pith.science/paper/D3IORFCC
@misc{pith2026260709050,
author = {Pith},
title = {Pith review of: Asymptotic stability of solutions to the good Boussinesq equation in dispersive wave region},
year = {2026},
howpublished = {\url{https://pith.science/paper/D3IORFCC}},
note = {Machine review of arXiv:2607.09050}
}
abstract
This work studies the asymptotic stability of solutions to the good Boussinesq equation in dispersive wave region when the reflection coefficients associated with the initial data belong to weighted Sobolev space. The Dbar-steepest descent method is applied to the Riemann-Hilbert problem and the long-time asymptotic expansion of the solution are obtained up to an optimal error of order $\mathcal{O}(t^{-3/4})$. Compared with previous results, we extend the initial data from the rapidly decaying Schwartz space to a weighted Sobolev space, and prove the asymptotic stability of the solution in dispersive wave region.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[13]
S.The “good” Boussinesq equation: long-time asymp- totics
Charlier, C., Lenells, J., and Wang, D. S.The “good” Boussinesq equation: long-time asymp- totics. Anal. PDE16(2023), no. 5, 1351–1388
2023
-
[48]
Bona, J. L. and Sachs, R. L. Global existence of smooth solutions and stability of solitary waves for a generalized Boussinesq equation. Commun. Math. Phys.118(1988), 15–29. 52
1988
-
[1]
S.Integrable nonlinear evolution equations on the half-line
Fokas, A. S.Integrable nonlinear evolution equations on the half-line. Commun. Math. Phys. 230(2002), no. 1, 1–39
2002
-
[2]
H.New integrable nonlinear evolution equations
Wadati, M., Konno, K., and Ichikawa, Y. H.New integrable nonlinear evolution equations. J. Phys. Soc. Japan47(1979), no. 5, 1698–1700
1979
-
[3]
E., Spector, M
Fal’kovich, G. E., Spector, M. D., and Turitsyn, S. K.Destruction of stationary solutions and collapse in the nonlinear string equation. Phys. Lett. A99(1983), no. 6, 271–274
1983
-
[4]
S., Mitchell, A
Manoranjan, V. S., Mitchell, A. R., and Morris, J. L.Numerical solutions of the good Boussi- nesq equation. SIAM J. Sci. Stat. Comput.5(1984), no. 4, 946–957
1984
-
[5]
M.Pseudospectral method for the “good” Boussinesq equation
De Frutos, J., Ortega, T., and Sanz-Serna, J. M.Pseudospectral method for the “good” Boussinesq equation. Math. Comp.57(1991), no. 195, 109–122
1991
-
[6]
S., Ortega, T., and Sanz-Serna, J
Manoranjan, V. S., Ortega, T., and Sanz-Serna, J. M.Soliton and antisoliton interactions in the “good” Boussinesq equation. J. Math. Phys.29(1988), no. 9, 1964–1968
1988
Show all 48 references
-
[7]
Ablowitz, M. J. and Segur, H.Solitons and the Inverse Scattering Transform. SIAM, Philadel- phia, 1981
1981
-
[8]
J., Kaup, D
Ablowitz, M. J., Kaup, D. J., Newell, A. C., and Segur, H.The inverse scattering transform- Fourier analysis for nonlinear problems. Stud. Appl. Math.53(1974), no. 4, 249–315
1974
-
[9]
C.The inverse scattering transform
Newell, A. C.The inverse scattering transform. In:Solitons, 177–242, Springer, Berlin, Heidelberg, 1980
1980
-
[10]
McLaughlin, K. T. R. and Miller, P. D.The ¯∂steepest descent method and the asymptotic behavior of polynomials orthogonal on the unit circle with fixed and exponentially varying nonanalytic weights. Int. Math. Res. Pap.2006(2006), 48673. 49
2006
-
[11]
McLaughlin, K. T. R. and Miller, P. D.The ¯∂steepest descent method for orthogonal poly- nomials on the real line with varying weights. Int. Math. Res. Not. IMRN2008(2008), rnn075
2008
-
[12]
and Lenells, J.The “good” Boussinesq equation: a Riemann-Hilbert approach
Charlier, C. and Lenells, J.The “good” Boussinesq equation: a Riemann-Hilbert approach. Indiana Univ. Math. J.71(2022), no. 4, 1505–1562
2022
-
[14]
and Lenells, J.Miura transformation for the “good” Boussinesq equation
Charlier, C. and Lenells, J.Miura transformation for the “good” Boussinesq equation. Stud. Appl. Math.152(2024), no. 1, 73–110
2024
-
[15]
Darrigol, O.Joseph Boussinesq’s legacy in fluid mechanics. C. R. Mec.345(2017), no. 7, 427–445
2017
-
[16]
B., Kirby, J
Kennedy, A. B., Kirby, J. T., and Chen, Q.Boussinesq-type equations with improved nonlin- ear performance. Wave Motion33(2001), no. 3, 225–243
2001
-
[17]
M.Multiple-soliton solutions for the Boussinesq equation
Wazwaz, A. M.Multiple-soliton solutions for the Boussinesq equation. Appl. Math. Comput. 192(2007), no. 2, 479–486
2007
-
[18]
Deift, P., Tomei, C., and Trubowitz, E.Inverse scattering and the Boussinesq equation. Comm. Pure Appl. Math.35(1982), no. 5, 567–628
1982
-
[19]
B., Alharbi, A
Almatrafi, M. B., Alharbi, A. R., and Tun¸ c, C.Constructions of the soliton solutions to the good Boussinesq equation. Adv. Difference Equ.2020(2020), no. 1, Paper No. 629
2020
-
[20]
Clarkson, P. A. and Dowie, E.Rational solutions of the Boussinesq equation and applications to rogue waves. Trans. Math. Appl.1(2017), no. 1, tnx003
2017
-
[21]
Bona, J. L. and Sachs, R. L.Global existence of smooth solutions and stability of solitary waves for a generalized Boussinesq equation. Comm. Math. Phys.118(1988), no. 1, 15–29
1988
-
[22]
and Zhou, X.A steepest descent method for oscillatory Riemann-Hilbert problems: Asymptotics for the MKdV equation
Deift, P. and Zhou, X.A steepest descent method for oscillatory Riemann-Hilbert problems: Asymptotics for the MKdV equation. Ann. of Math. (2)137(1993), no. 2, 295–368
1993
-
[23]
Q., Tian, S
Li, Z. Q., Tian, S. F., and Yang, J. J.On the soliton resolution and the asymptotic stability of N-soliton solution for the Wadati-Konno-Ichikawa equation with finite density initial data in space-time solitonic regions. Adv. Math.409(2022), 108639
2022
-
[24]
and Yan, Z
Zhang, M H. and Yan, Z. Y.The Sasa-Satsuma equation with high-order discrete spectra in space-time solitonic regions: soliton resolution via the mixed ¯∂-Riemann-Hilbert problem. Commun. Theor. Phys.76(2024), no. 6, 065002
2024
-
[25]
G., Wang, J., Wang, K
Geng, X. G., Wang, J., Wang, K. D., and Li, R. M.Soliton resolution for the complex short-pulse positive flow with weighted Sobolev initial data in the space-time soliton regions. J. Differ. Equ.386(2024), 214–268. 50
2024
-
[26]
Q., Tian, S
Li, Z. Q., Tian, S. F., and Yang, J. J.Soliton resolution for the Wadati-Konno-Ichikawa equation with weighted Sobolev initial data. Ann. Henri Poincar´ e23(2022), no. 7, 2611– 2655
2022
-
[27]
R.The Riemann-Hilbert problem and integrable systems
Its, A. R.The Riemann-Hilbert problem and integrable systems. Notices Amer. Math. Soc. 50(2003), no. 11, 1389–1400
2003
-
[28]
Zhou, X.The Riemann-Hilbert problem and inverse scattering. SIAM J. Math. Anal.20 (1989), no. 4, 966–986
1989
-
[29]
Kashiwara, M.The Riemann-Hilbert problem for holonomic systems. Publ. Res. Inst. Math. Sci.20(1984), no. 2, 319–365
1984
-
[30]
Borghese, M., Jenkins, R., and McLaughlin, K. D. T.-R.Long time asymptotic behavior of the focusing nonlinear Schr¨ odinger equation. Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire 35(2018), no. 4, 887–920
2018
-
[31]
J., and Guo, B
Liu, N., Chen, M. J., and Guo, B. L.Long-time asymptotic behavior of the fifth-order modified KdV equation in low regularity spaces. Stud. Appl. Math.147(2021), no. 1, 230–299
2021
-
[32]
and McLaughlin, K
Dieng, M. and McLaughlin, K. D. T.-R.Long-time asymptotics for the NLS equation via dbar methods. Preprint, arXiv:0805.2807
-
[33]
Wang, D. S. and Zhu, X. D.Long-time asymptotics of the good Boussinesq equation with qxx-term and its modified version. J. Math. Phys.63(2022), no. 12, 123501
2022
-
[34]
S., and Zhu, X
Huang, L., Wang, D. S., and Zhu, X. D.Long-time asymptotics of the Tzitz´ eica equation on the line. Preprint, arXiv:2404.04999
-
[35]
Yang, Y. L. and Fan, E. G.On the long-time asymptotics of the modified Camassa-Holm equation in space-time solitonic regions. Adv. Math.402(2022), 108340
2022
-
[36]
Wang, Z. Y. and Fan, E. G.The defocusing nonlinear Schr¨ odinger equation with a nonzero background: Painlev´ e asymptotics in two transition regions. Comm. Math. Phys.402(2023), no. 3, 2879–2930
2023
-
[37]
Nonlinearity30(2017), no
Giavedoni, P.Long-time asymptotic analysis of the Korteweg-de Vries equation via the dbar steepest descent method: the soliton region. Nonlinearity30(2017), no. 3, 1165–1196
2017
-
[38]
Yang, Y. L. and Fan, E. G.Soliton resolution and large time behavior of solutions to the Cauchy problem for the Novikov equation with a nonzero background. Adv. Math.426(2023), 109088
2023
-
[39]
Q., Tian, S
Li, Z. Q., Tian, S. F., and Yang, J. J.On the asymptotic stability ofN-soliton solution for the short pulse equation with weighted Sobolev initial data. J. Differ. Equ.377(2023), 121–187
2023
-
[40]
¯∂-problem for the focusing nonlinear Schr¨ odinger equation and soliton shielding
Bertola, M., Grava, T., and Orsatti, G. ¯∂-problem for the focusing nonlinear Schr¨ odinger equation and soliton shielding. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci.481(2025), no. 2310, 20240764. 51
2025
-
[41]
Q., Tian, S
Li, Z. Q., Tian, S. F., and Yang, J. J.On the long-time asymptotic for a coupled generalized nonlinear Schr¨ odinger equations with weighted Sobolev initial data. Phys. D456(2023), 133915
2023
-
[42]
and Liu, J
Chen, G. and Liu, J. Q.Long-time asymptotics of the modified KdV equation in weighted Sobolev spaces. Forum Math. Sigma10(2022), e66
2022
-
[43]
and Wang, L.The coupled Hirota equations with a3×3Lax pair: Painlev´ e-type asymptotics in transition zone
Zhao, X. and Wang, L.The coupled Hirota equations with a3×3Lax pair: Painlev´ e-type asymptotics in transition zone. Stud. Appl. Math.153(2024), no. 3, e12745
2024
-
[44]
and Fan, E
Xu, J. and Fan, E. G.Initial-boundary value problem for integrable nonlinear evolution equa- tion with3×3Lax pairs on the interval. Stud. Appl. Math.136(2016), no. 3, 321–354
2016
-
[45]
G.The initial boundary value problems for a nonlinear integrable equation with3×3Lax pair on the finite interval
Xiao, Y., Xu, J., and Fan, E. G.The initial boundary value problems for a nonlinear integrable equation with3×3Lax pair on the finite interval. Acta Math. Sci. Ser. B (Engl. Ed.)41 (2021), no. 5, 1733–1748
2021
-
[46]
and Coifman, R
Beals, R. and Coifman, R. R.Scattering and inverse scattering for first order systems. Comm. Pure Appl. Math.37(1984), no. 1, 39–90
1984
-
[47]
Lenells, J.Matrix RH problems with jumps across Carleson contours. Monatsh. Math.186 (2018), no. 1, 111–152
2018
Reviewed July 13, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.