{"id":"d47c4fb7-937d-4912-9c37-78928534aa78","arxiv_id":"2504.19663","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Any Schwartz-class solution of the bad Boussinesq equation on the half-line is recovered from its initial-boundary data via a 3x3 Riemann-Hilbert problem, under no-soliton and generic-behavior assumptions.","lead":"This paper proves that, assuming a smooth solution exists, the bad Boussinesq equation on a half-line can be reconstructed from its initial-boundary values by solving a matrix Riemann-Hilbert problem. This extends the Fokas method to a new integrable equation with a 3x3 Lax pair, opening a path to long-time analysis.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The recovery theorem rests on an unproved Schwartz-solution assumption that is non-generic for the linearly ill-posed bad Boussinesq equation, making the inverse result conditional on a possibly empty hypothesis.","rationale":"The reader identified the same weakest assumption: the paper assumes, rather than proves, the existence of a Schwartz-class solution, and the bad Boussinesq equation is linearly ill-posed. This is indeed the most load-bearing concern because the main theorems are conditional on it. In good faith, the paper is transparent about this conditionality, and the internal RH construction appears consistent: the direct problem builds spectral functions, the inverse problem constructs M, and Appendix A provides a detailed uniqueness argument. There are also many deferred proofs (e.g., Propositions 3.1, 3.4, 3.5, 4.1), but these are presented as standard or analogous to prior work and are not the primary weakness. The primary weakness is that the theorem's hypothesis may be vacuous or highly restrictive for the ill-posed equation, so the central claim is a conditional representation rather than an unconditional recovery result. A concrete linearized test can clarify whether nontrivial Schwartz solutions exist for t > 0. Since the reader's CONDITIONAL verdict already reflects this, no change is needed.","tokens_in":44313,"tokens_out":7468,"duration_ms":82889,"concrete_test":"Test the non-vacuity of Definition 2.1 at the linear level. Take u0(x) = e^{-x²}, u1(x) = 0 on R+ with compatible smooth boundary values, and solve the linearized bad Boussinesq equation utt = uxx + uxxxx on the half-line via Fourier/sine transform. The linear solution has Fourier modes with factor cosh(t sqrt(k⁴ − k²)), which grows like exp(t k²) for large k. Check whether the solution at any t > 0 is still in the Schwartz class S(R+). If it is not, then generic Schwartz initial data fail Definition 2.1, and the paper's recovery theorem applies only to a special, uncharacterized class of data; this would confirm that the existence hypothesis is a serious restriction and should be stated and proved as a separate theorem for a nontrivial class.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, Theorem 2.7, is conditional on the existence of a Schwartz-class solution, as stated in Definition 2.1 and the sentence 'we suppose that there exists a Schwartz class solution' in Section 2. This assumption is load-bearing because the bad Boussinesq equation is linearly ill-posed: the linearized dispersion relation ω² = k² − k⁴ gives high-frequency modes a growth factor like exp(t|k|²). Consequently, generic Schwartz initial data do not remain Schwartz for any t > 0, so the set of data satisfying Definition 2.1 together with Assumptions 2.2 and 2.4 may be very small or even empty for T > 0. The paper acknowledges this in the abstract ('Assuming that the solution exists'), so this is a limitation rather than a contradiction. However, it means the inverse scattering representation does not by itself establish that arbitrary initial-boundary data can be mapped to scattering data for a genuine solution; the entire construction applies only to solutions that are presupposed to exist. The RH problem itself is internally well-posed (existence by construction, uniqueness by Appendix A), but the theorem's practical content is weakened by the lack of any existence or non-vacuity result for the assumed solution class.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a direct and inverse scattering formalism for the Boussinesq equation on the half-line. Assuming the existence of a Schwartz-class solution with appropriate spectral assumptions, the authors define nine reflection coefficients from initial and boundary values via Volterra integral equations, and construct a 3×3 Riemann-Hilbert problem whose solution recovers u(x,t) through asymptotic formulas. The main results are Theorem 2.5 on the properties of the reflection coefficients and Theorem 2.7 on the recovery of u from the RH solution. The paper also proves an equivalence between the Boussinesq equation and a first-order system of Zakharov, and includes a uniqueness proof for the RH problem in Appendix A.","tokens_in":44628,"tokens_out":2363,"duration_ms":26019,"significance":"If the conditional results are accepted, the paper provides the first half-line inverse scattering representation for the bad Boussinesq equation, extending the Fokas method to a 3×3 Lax pair with a complicated 54-arc contour. The explicit jump matrices, the symmetry reductions, and the uniqueness argument for the RH problem are concrete technical achievements. The equivalence lemma between the scalar equation and the first-order system is cleanly proved. However, the central theorem is conditional on an unproved existence assumption that may be non-generic for this linearly ill-posed equation, which reduces the practical scope of the result. The paper is honest about this limitation in the abstract, but the manuscript would be substantially stronger with a non-vacuity check or a detailed discussion of the admissible data class.","major_comments":[{"comment":"The recovery theorem, Theorem 2.7, is conditional on the existence of a Schwartz-class solution of the half-line IBVP, and no existence or non-vacuity result is provided. This is load-bearing because the linearized bad Boussinesq equation has the dispersion relation ω² = k² − k⁴, so high-frequency modes grow like exp(t|k|²); generic Schwartz initial data therefore do not remain Schwartz for t > 0. The paper should either prove local existence in the Schwartz class for a nontrivial set of data satisfying Assumptions 2.2 and 2.4, or provide explicit nontrivial solutions of the assumed class. Without this, the inverse result applies only to solutions that are presupposed to exist, and the reader cannot tell whether the hypothesis is empty.","section":"Section 2, Definition 2.1 and the sentence \"we suppose that there exists a Schwartz class solution\""},{"comment":"These propositions contain the basic analytic properties of the eigenfunctions and of the sectionally defined function M, including domains of definition, boundedness, symmetries, and asymptotic behavior near k = ±1 and near the sixth roots of unity. Their proofs are omitted and referred to prior work or to 'similar' statements. Since the half-line setting requires a new analysis of three eigenfunctions and the contour has many intersection points, the omitted details are not routine for the reader. Please include complete proofs or provide precise theorem statements from the cited papers together with a verification that all hypotheses of those theorems are satisfied in the present setting.","section":"Propositions 3.1, 3.4, and 4.1"},{"comment":"The symmetry relations (2.13a)–(2.13e) are asserted to follow from a 'long but direct computation' that is not presented. These relations are essential: they are used in the proof of the jump condition for M (Lemma 4.5) and in the construction of the jump matrices. The computation should be included, at least in an appendix, because an error in these identities would invalidate the RH formulation.","section":"Theorem 2.5(iv), equations (2.13)"},{"comment":"The main theorems hold only under Assumption 2.2 (absence of solitons) and Assumption 2.4 (generic behavior at k = ±1). The manuscript does not discuss whether these assumptions are preserved by the time evolution or how restrictive they are. Because the inverse problem is stated for arbitrary initial-boundary data satisfying the existence assumption, the authors should clarify the admissible data class and, if possible, give examples of nontrivial data satisfying all hypotheses simultaneously.","section":"Assumptions 2.2 and 2.4"}],"minor_comments":[{"comment":"There are typographical errors in the abstract: 'Boussin esq' should be 'Boussinesq' and 't hat' should be 'that'; the phrase 'via the solution of a 3 × 3 Riemann-Hilbert problem' would read better without the extra spaces.","section":"Abstract"},{"comment":"The definition of the contour Γ and its orientation would benefit from a more explicit written description of the orientation of each arc and half-line, beyond Figure 1, since the jump condition M₊ = M₋ v depends on that orientation.","section":"Section 2.2, RH problem 2.6"},{"comment":"The functions h, g, h̃, g̃ are defined without motivation; a sentence explaining their role as auxiliary combinations of reflection coefficients would improve readability.","section":"Equation (2.15)"},{"comment":"The word 'interger' should be 'integer' in the sentence 'For each interger l ≥ −1'.","section":"Proposition 3.8(f)"},{"comment":"Reference [7] is cited as an arXiv preprint; if a published version exists, it should be cited to give the reader a stable reference.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper relies heavily on the author's own prior work for the proofs of several central propositions. This is not a circularity problem, since the reflection coefficients are defined independently of the desired output, but it does mean that the present manuscript's contribution is partly an assembly and adaptation of existing technical machinery. The deepest gap is the unproved existence of Schwartz-class solutions with the required spectral assumptions; this should be addressed head-on either by a new existence result or by a clear statement that the paper is a conditional representation theorem and that the data class may be small. The paper is likely to be of interest to the integrable-systems community if these issues are resolved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the first inverse-scattering representation for the bad Boussinesq equation on the half-line, and the result is honestly stated as conditional. The paper assumes a Schwartz-class solution exists and proves it can be recovered from the initial-boundary data via a 3x3 Riemann-Hilbert problem. The conditionality is explicit in the abstract and in Section 2, so it is a limitation the authors are not hiding.\n\nThe genuinely new piece is the adaptation of Lenells' 3x3 Fokas method to this equation: a 54-arc contour, 36 sectors, nine reflection coefficients, and the associated jump matrices. The machinery is organized carefully, and the uniqueness theorem for the RH problem in Appendix A is a real proof, not a sketch. Lemma 2.9, which shows the equivalence with Zakharov's system and cuts the boundary data down to four functions, is proved cleanly. Within the Charlier-Lenells program this is a natural and substantial step, and the references to their earlier papers are legitimate scaffolding rather than a substitute for the argument.\n\nThe soft spots, in proportion. The central assumption is load-bearing and, for this equation, non-trivial: bad Boussinesq is linearly ill-posed (the linearized dispersion gives growth like exp(t|k|²)), so the class of Schwartz solutions on [0,T] is plausibly small and is certainly not established to be non-empty. The stress-test concern lands here, although the paper's explicit caveat blunts the criticism. The conditional format is standard in this IBVP literature, but the paper provides no non-vacuity result or example showing the hypotheses are ever satisfied. A referee should push on this.\n\nTwo presentation gaps. Propositions 3.1, 3.4, and 4.1 are stated without proofs and deferred to earlier work; acceptable within a series, but it makes the paper hard to assess standalone. More pointedly, Theorem 2.5(iv) — the global relations in (2.13) that make the jump matrix consistent across the 54-arc contour — rests on 'a long but direct computation.' Enough depends on those identities that a referee should ask for the computation, or at least a summary of how it runs. Minor typos throughout are annoying but irrelevant.\n\nWho this is for: integrable-systems people working on Fokas-method IBVPs with 3x3 Lax pairs, and anyone planning long-time asymptotics for Boussinesq. The reader's conditional verdict is about right. I would send it to a serious referee, with instructions to focus on the existence-class question and the omitted global-relation computation.","headline":"First half-line inverse scattering for the bad Boussinesq equation: an honestly conditional Fokas-method construction that deserves a serious referee, with the unproved existence hypothesis and deferred computations as the main soft spots.","tokens_in":45058,"tokens_out":6215,"would_cite":true,"duration_ms":62280,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35G31","35Q15","37K15","76B15"],"pacs":[],"model":"deepseek-v4-flash","headline":"A 3x3 Riemann-Hilbert problem recovers half-line Boussinesq solutions from initial-boundary data.","keywords":["Riemann-Hilbert problem","direct and inverse scattering","initial-boundary value problem","Boussinesq equation","half-line","3x3 Lax pair","reflection coefficients"],"falsifier":"Take an explicit known half-line solution, for instance the zero solution with zero initial and boundary data, compute the nine reflection coefficients from (2.9)-(2.10), solve RH problem 2.6, and compare both recovery formulas in (2.25) with the known solution at many $(x,t)$; a single mismatch would show the claimed recovery is false.","tokens_in":44140,"feed_emoji":"🌊","tokens_out":10749,"duration_ms":100671,"temperature":0.7,"pith_summary":"The paper studies the Boussinesq equation $u_{tt}=u_{xx}+(u^2)_{xx}+u_{xxxx}$ on the half-line $x\\ge 0$, with initial data prescribed at $t=0$ and four boundary values prescribed at $x=0$. It proves that, assuming a smooth rapidly decaying solution exists, that solution is uniquely recovered from these initial-boundary values through a $3\\times3$ Riemann-Hilbert problem. The scattering data consist of nine reflection coefficients, and the reconstruction is given by two explicit limits of the Riemann-Hilbert solution. This transfers the inverse-scattering description of Boussinesq solutions from the whole line to the half-line. The statement is conditional on the assumed existence of the solution, which the paper does not establish.","feed_headline":"3x3 Riemann-Hilbert problem recovers half-line Boussinesq waves","feed_subtitle":"Nine reflection coefficients built from initial and boundary data determine the wave, provided a smooth solution exists.","key_machinery":"The load-bearing object is a $3\\times3$ Riemann-Hilbert problem for a matrix function $M(x,t,k)$ whose jump contour $\\Gamma$ consists of 18 arcs on the unit circle, 18 straight segments, and 18 half-lines. The problem is built from a Lax pair with diagonal matrices $L=\\operatorname{diag}(l_1,l_2,l_3)$ and $Z=\\operatorname{diag}(z_1,z_2,z_3)$, whose compatibility condition is the equation; Volterra integral equations along three contours define eigenfunctions $\\mu_1,\\mu_2,\\mu_3$ and adjoint eigenfunctions, and their values at the corners $(0,0)$, $(0,T)$, and $(+\\infty,t)$ produce spectral matrices $s,S,s^A,S^A$. Ratios of entries of these matrices define the nine reflection coefficients, and the jump matrix $v(x,t,k)=e^{x\\hat L+t\\hat Z}\\tilde v(k)$ is assembled from them. The recovery formulas come from the first two coefficients in the expansion of $M$ at $k=\\infty$, while the symmetries under $k\\mapsto\\omega k$ and $k\\mapsto k^{-1}$ and the prescribed pole behavior at $k=\\pm1$ make the problem uniquely solvable.","core_discovery":"The central claim is Theorem 2.7: if $u$ is a Schwartz-class solution of (1.1) on $[0,\\infty)\\times[0,T]$ with initial data $u_0,u_1\\in S(\\mathbb{R}_+)$ and boundary data $\\tilde u_0,\\tilde u_1,\\tilde u_2,\\tilde u_3\\in C^\\infty([0,T])$ satisfying Assumptions 2.2 (no solitons) and 2.4 (generic behavior at $k=\\pm1$), then the nine reflection coefficients defined by (2.9)-(2.10) give rise to a unique solution $M(x,t,k)$ of RH problem 2.6, and $u$ is recovered by $u(x,t)=-i\\sqrt{3}\\,\\partial_x\\lim_{k\\to\\infty} k[(M)_{33}-1]=\\frac{1-\\omega}{2}\\lim_{k\\to\\infty} k^2(M)_{32}$. The companion direct-scattering theorem establishes smoothness, pole and zero structure, boundary asymptotics, and algebraic identities for the reflection coefficients, and the paper shows the same results can be reformulated for the equivalent first-order system $v_t=u_x+(u^2)_x+u_{xxx}$, $u_t=v_x$.","pith_inferences":["If existence theory for this half-line problem is developed, the Riemann-Hilbert representation would immediately supply a route to long-time asymptotics, mirroring how the whole-line representation was used.","The identities (2.13) resemble admissibility conditions for boundary data; they could be tested numerically as constraints that any dataset supporting a solution must satisfy.","The reconstruction formulas give a concrete numerical algorithm: discretize the Riemann-Hilbert problem and compare the output with a direct finite-domain PDE solver, which would probe both the inverse-scattering construction and the assumed existence."],"forward_implications":["Any solitonless Schwartz-class half-line solution satisfying the two generic assumptions is completely determined by its initial-boundary data through the nine reflection coefficients.","The unique solution of RH problem 2.6 encodes the solution at every point $(x,t)$, so evaluating the two limits in (2.25) gives $u(x,t)$ directly.","The equivalence with the first-order system (2.26) means the same Riemann-Hilbert construction applies to $u$ and its integrated companion $v$.","The direct-scattering theorem records the pole and zero structure and the symmetries of the reflection coefficients, which is what makes the inverse problem well-posed.","The uniqueness proof for the Riemann-Hilbert problem means the map from initial-boundary data to scattering data is injective on the class of solutions considered."],"supporting_citations":[{"why":"supplies the line-case direct and inverse scattering analysis whose spectral functions this paper extends to the half-line.","marker":"[7]"},{"why":"provides the Lax pair for the equivalent first-order system, the compatibility condition underlying the whole construction.","marker":"[25]"},{"why":"gives the method for solving half-line initial-boundary value problems with 3x3 Lax pairs that organizes the Riemann-Hilbert problem.","marker":"[20]"},{"why":"introduces the unified transform approach for half-line integrable PDEs from which the boundary-value treatment grows.","marker":"[13]"},{"why":"supplies the Riemann-Hilbert and Volterra estimates for the Boussinesq equation that the technical propositions adapt.","marker":"[6]"},{"why":"gives the earlier half-line treatment of the related Boussinesq variant against which this 3x3 extension is positioned.","marker":"[16]"}],"fun_headline_variants":["Half-line Boussinesq solved via 3x3 Riemann-Hilbert","Boussinesq on half-line: recovery via 9 reflection coefficients","3x3 RH problem reconstructs half-line Boussinesq from data","New RH method for half-line Boussinesq with 18 arcs","Recovering half-line Boussinesq from boundary data via RH"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof assumes from the start that a Schwartz-class solution of the half-line initial-boundary value problem exists; if no such solution exists for the chosen data, the scattering data and the Riemann-Hilbert reconstruction have no object to describe.","fun_headline_variants_meta":{"raw":{"variants":["Half-line Boussinesq solved via 3x3 Riemann-Hilbert","Boussinesq on half-line: recovery via 9 reflection coefficients","3x3 RH problem reconstructs half-line Boussinesq from data","New RH method for half-line Boussinesq with 18 arcs","Recovering half-line Boussinesq from boundary data via RH"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000294,"raw_usage":{"total_tokens":1677,"prompt_tokens":878,"completion_tokens":799,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":494,"completion_tokens_details":{"reasoning_tokens":700}},"tokens_in":494,"tokens_out":799,"duration_ms":6353,"temperature":1.0,"reasoning_tokens":700,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:46:00.922920+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take an explicit known half-line solution, for instance the zero solution with zero initial and boundary data, compute the nine reflection coefficients from (2.9)-(2.10), solve RH problem 2.6, and compare both recovery formulas in (2.25) with the known solution at many $(x,t)$; a single mismatch would show the claimed recovery is false.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the Lax pair for the equivalent first-order system, the compatibility condition underlying the whole construction."},{"cited_title":"Lenells, Initial-boundary value problems for integrable evolut ion equations with 3 × 3 Lax pairs, Physica D 241 (2012), 857–875","cited_arxiv_id":null,"evidence_quote":"gives the method for solving half-line initial-boundary value problems with 3x3 Lax pairs that organizes the Riemann-Hilbert problem."},{"cited_title":"Fokas, A unified transform method for solving linear and cer tain nonlinear PDEs, Proc","cited_arxiv_id":null,"evidence_quote":"introduces the unified transform approach for half-line integrable PDEs from which the boundary-value treatment grows."},{"cited_title":"Charlier and J","cited_arxiv_id":null,"evidence_quote":"supplies the Riemann-Hilbert and Volterra estimates for the Boussinesq equation that the technical propositions adapt."},{"cited_title":"Himonas and D","cited_arxiv_id":null,"evidence_quote":"gives the earlier half-line treatment of the related Boussinesq variant against which this 3x3 extension is positioned."}],"review_version":1}