{"id":"f3790e06-a910-4d8f-804d-0bb7a47b479d","arxiv_id":"2412.11581","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Rescaled by 1/N, large-order rational solitons of the focusing complex mKdV equation converge to a universal limit function that solves c-mKdV in (X,T) and ODEs from the Painlevé-III hierarchy.","lead":"This paper works out what happens to very tall rational soliton waves of the complex modified Korteweg-de Vries equation as their order grows: in a rescaled frame they converge to a universal profile governed by the c-mKdV equation and by special equations from the Painlevé-III hierarchy. The new formulas describe the height, phase, and transition behavior of these extreme waves, which matter for predicting rogue-wave-like structures in optics and plasmas.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper asserts, but does not prove, that RHP 2 reproduces the explicit determinantal formula (1.19) for the multi-rational solitons; Theorem 1.2 as stated about the explicit family is therefore not yet established.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the asserted equivalence between the determinantal soliton formula (1.19) and RHP 2. I find no additional internal inconsistency in the asymptotic machinery; the small-norm RHP argument, the Laurent expansions, and the zero-curvature computations leading to the ODEs are algebraically consistent in the sections I spot-checked. The main unresolved point is a proof gap rather than a demonstrated contradiction: the paper does not show that the RHP 2 family coincides with the explicit formula (1.19) that Theorem 1.2 refers to. This warrants a conditional verdict, which is exactly what the reader issued. A concrete numerical or analytic check of this equivalence would settle the question and is inexpensive to perform. I therefore recommend keeping the verdict unchanged while requesting this verification.","tokens_in":50400,"tokens_out":13150,"duration_ms":116241,"concrete_test":"Solve RHP 2 numerically for k=1 and k=2 (n=1) using a standard RHP solver (e.g., Olver-Trogdon) on a grid of (X,T) in, say, [-5,5]^2, and compare the recovered q_1(x,t) and q_2(x,t) with formula (1.4) from Chen and Yan [21]. Agreement to numerical precision (e.g., relative error below 1e-8) would strongly support the identification; a systematic discrepancy would refute it. As an analytic complement, derive the Darboux/dressing transformation that converts the bare Jost-solution RHP into RHP 1 and verify that it reproduces the determinant formula (1.19) for k=1 and k=2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 1.2 states a large-order asymptotic for the kth-order rational solitons q_k given by the closed-form determinantal formula (1.19) imported from Chen and Yan [21]. The actual asymptotic analysis, however, is performed on the family defined by RHP 2 via Proposition 2.3, which gives q_k = 1 + 2i lim_{Lambda->∞} Lambda Mhat^(k)_12(Lambda;x,t). Proposition 2.3 derives this trace formula from RHP 1, but it does not show that RHP 1 (or RHP 2) produces the same solution as formula (1.19). The jump matrix of RHP 1 is written in terms of scattering data and a gauge transformation, yet the passage from the robust IST construction in [21] to this particular RHP is not demonstrated. If the RHP-generated family differs from (1.19) at finite order, then the O(1/n) convergence in Theorem 1.2 concerns a different object than the theorem's hypothesis. This gap is load-bearing because all downstream results—the limit c-mKdV equation, the Painleve-III ODEs, and the large-X asymptotics—describe the RHP-3 limit of the RHP-2 family. The gap is independent of the internal consistency of the asymptotic analysis, which otherwise appears coherent in the parts I checked, including the small-norm RHP argument and the zero-curvature computations in Section 3.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the large-order behavior of multi-rational solitons of the focusing complex mKdV equation with nonzero background. The authors construct a Riemann-Hilbert problem (RHP 1, then RHP 2 after removing the branch-cut jump) that is claimed to reproduce the kth-order rational solitons given by the closed-form determinant formula (1.19) imported from Chen and Yan. Through the rescaling X=nx, T=n^3 t, Λ=λ/n they derive a model RHP (RHP 3), prove existence and uniqueness, and show that the reconstruction q̂±(X,T) solves the c-mKdV equation in the rescaled variables. They then derive ODEs in X and T satisfied by q̂±, identify the T=0 X-ODE with the first member of the Painlevé-III hierarchy, and give large-X and transitional-asymptotic formulas, with a partial analysis for large T. The central claim, Theorem 1.2, is that (1/n) q_{2n}(X/n,T/n^3) and (1/n) q_{2n-1}(X/n,T/n^3) converge to q̂±(X,T) uniformly on compact sets.","tokens_in":50635,"tokens_out":9051,"duration_ms":87285,"significance":"If the identification between the explicit determinant solutions (1.19) and the RHP-generated family is established, the paper would provide a substantial extension to the complex mKdV equation of the large-order universality program developed for the focusing NLS equation by Bilman, Buckingham, Miller, and others. The asymptotic machinery is presented in detail: the small-norm RHP argument in Section 2, the zero-curvature computations in Section 3, and the steepest-descent parametrices in Section 4 all follow the established framework, and the asymptotic constants are derived from model problems rather than fitted to the target. The authors are also appropriately explicit that the large-T analysis is only partial. These strengths make the manuscript worth serious revision rather than rejection, but the missing equivalence proof for the object whose asymptotics are computed is load-bearing.","major_comments":[{"comment":"Theorem 1.2 is stated for the kth-order rational solitons q_k given by the explicit deterministic formula (1.19) from Chen and Yan [21], but the proof is carried out entirely for the family defined by RHP 2. Proposition 2.3 gives the reconstruction formula q_k = 1 + 2i lim_{Λ→∞} Λ M̂^(k)_{12}(Λ;x,t) from RHP 2, but the paper does not prove that this RHP family coincides with formula (1.19). In particular, RHP 1 is a pure discrete-spectrum RHP with jump data built from ((λ−i)/(λ+i))^{nσ_3}, yet the paper does not show that the multi-gauge-transformed IST solution in [21] has exactly this scattering data and no additional data. If the two families differ at finite order, then the O(1/n) convergence in Theorem 1.2 concerns a different sequence of functions than the one named in the theorem. This gap is load-bearing for Theorem 1.2 and, through it, for the interpretation of all downstream statements about q̂±. The authors should either supply a proof of equivalence or reformulate the theorem explicitly for the RHP-2 family and state the identification with (1.19) as a separate conjecture.","section":"§2, Theorem 1.2, Propositions 2.1–2.3"},{"comment":"The displayed derivation of the X-ODEs in Theorem 1.3 cannot be verified as printed. Equation (3.115) is not the X-derivative of Eq. (3.113): the expression contains a spurious '− +' before the term 6q̂_{XX}|q̂|², and the terms 6q̂²q̂* and 6|q̂|²q̂ appear without the derivative symbols that a differentiation of (3.113) would produce. Since Eqs. (3.117) and (3.120) are central to the identification of the Painlevé-III hierarchy, the authors should correct the displayed algebra and clarify each substitution used in the elimination leading to (3.120).","section":"§3.2.1, Eqs. (3.115)–(3.117)"}],"minor_comments":[{"comment":"There are numerous typographical errors that should be corrected, including 'gvien' after Eq. (1.5), 'with with' before Eq. (2.32), 'ﬁst' for 'first' in the discussion of Eq. (1.27), 'euqations' in Section 5, 'impies' in the proof of Proposition 3.5, and 'Propositio n' in the proof of Proposition 4.3.","section":"Throughout"},{"comment":"The large-T section is presented as 'part results', but it stops abruptly after the g-function construction and the statement that new model problems are needed. The section would be clearer if it explicitly summarized which statements are proved and which are left for future work.","section":"§4.2"},{"comment":"The phrase 'double simple real critical points' in items (ii) and (iii) is confusing; it should be 'two simple real critical points' or 'a pair of simple real critical points'.","section":"§4.1.1"},{"comment":"The uniform bound (1.20) for N± on |Λ|≠1 over compact sets K is used in Lemma 2.1, but the proof of Theorem 1.1 via the vanishing lemma establishes uniqueness and existence and does not explicitly justify this bound; a sentence explaining how the bound follows from continuity and the normalization at infinity would help.","section":"Theorem 1.1, Eq. (1.20)"}],"recommendation":"major_revision","confidential_remarks":"The central issue is not internal inconsistency of the asymptotic analysis but a missing identification: the paper proves asymptotics for the RHP-2 family while stating the theorem for the explicit determinant family (1.19). This is fixable within the scope of the manuscript by adding a proof or by restating the main theorems, so I do not recommend rejection. However, because Theorem 1.2 is the advertised main result and the gap propagates to the interpretation of the Painlevé-III and asymptotic statements, the revision must address it before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a genuine extension of the Bilman–Ling–Miller large-order soliton program to the focusing c-mKdV equation with nonzero background, and the asymptotic machinery is coherent. The one thing to watch is a load-bearing identification gap at the front of Section 2. Theorem 1.2 is stated for the explicit determinantal family (1.19), but the proof runs entirely on the RHP 2 family. Proposition 2.3 only converts RHP 1 into RHP 2; neither proposition shows that these RHPs produce the same solutions as formula (1.19) from Chen–Yan. The stress-test note is right about this, and it is not a nitpick: if the RHP family differed from (1.19) at finite order, the O(1/n) theorem would be about a different object. I think the identification is true and provable from the IST machinery the authors cite, so treat it as fixable, but it needs to be proven or the theorem restated.\n\nWhat is actually new: the phase function ΛX + 4Λ³T ± 2/Λ gives a five-region large-X structure rather than the NLS case's regions; the limit RHP 3 comes with an existence/uniqueness proof; q̂± solves the c-mKdV equation in rescaled variables; the X- and T-ODEs come out of zero-curvature with α²+β²=4; at T=0 the X-ODE is identified with the first member of the P-III hierarchy; and the transition at a = −1/96 gives P-II asymptotics. The constants—p = ln2/(2π), the modulus √(2p), the 64 in (3.124)—are derived, not fitted to the target asymptotics, and that is the strongest evidence the derivation is sound. The algebra I spot-checked (route to (3.118) and the elimination to (3.120)) is consistent.\n\nSofter spots, in order. (1) The identification gap above; it is the only issue that should block acceptance. (2) Zhou's vanishing lemma for RHP 3 is dispatched in a few lines; probably right, but compressed for a load-bearing step. (3) The abstract's claim about the Painlevé-III hierarchy is only demonstrated at T=0; the theorems themselves are careful about this, the abstract is not. (4) The paper itself flags what is left undone: the large-T section ends with \"new models needed,\" and numerics are deferred. A numerical check of Theorem 1.2, or of the large-X formulas against the RHP solution, would raise confidence but is not a blocker.\n\nWho it is for: researchers in Riemann–Hilbert asymptotics for integrable systems, particularly the large-order/infinite-order soliton program. It deserves a serious referee. Send it out; ask that the identification gap be closed (prove the equivalence or restate Theorem 1.2 for the RHP-generated family) before acceptance.","headline":"Genuine extension of the NLS large-order/rational-soliton program to the focusing c-mKdV equation, with a coherent asymptotic core; the real issue is a load-bearing identification gap between the determinantal family (1.19) and the RHP that is actually analyzed.","tokens_in":51334,"tokens_out":7757,"would_cite":true,"duration_ms":65555,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q51","35Q15","37K40","37K10"],"pacs":[],"model":"deepseek-v4-flash","headline":"Large-order rational solitons of the focusing complex mKdV equation converge, after rescaling by 1/n, to a single universal profile governed by the c-mKdV equation and a Painlevé-III hierarchy.","keywords":["Complex mKdV equation","Nonzero background","Lax pair","Inverse scattering transform","Riemann-Hilbert problem","Multi-rational solitons","Large-order asymptotics","Painlevé-III hierarchy"],"falsifier":"Evaluate the closed-form formula (1.19) for $n=10,20,40$ at a fixed compact set of $(X,T)$, compute $(1/n)q_k(X/n,T/n^3)$, and compare with a numerical solution of the model RHP 3: the difference should shrink like $O(1/n)$. A persistent nonzero difference would show that RHP 2 and formula (1.19) define different solution families, invalidating the application of the theorems to the explicit solitons.","tokens_in":50079,"feed_emoji":"🌊","tokens_out":10079,"duration_ms":84099,"temperature":0.7,"pith_summary":"This paper establishes a large-order limit for the multi-rational solitons of the focusing complex modified Korteweg–de Vries (c-mKdV) equation with nonzero background: after rescaling space and time by $X=nx$, $T=n^3t$ and dividing the soliton by $n$, the $k$th-order rational soliton converges, uniformly on compact sets, to a universal function $\\hat q_\\pm(X,T)$ that is independent of the fine details of the soliton formula. The limit function is obtained as the unique solution of a model Riemann–Hilbert problem and is itself a new global solution of the same c-mKdV equation in the rescaled variables. At fixed $T=0$, the spatial profile connects to the first member of the Painlevé-III hierarchy through an ODE for $f=\\partial_X\\ln\\hat q$. The paper also derives explicit large-$X$ asymptotics and a transitional asymptotics near the critical parameter $a=-1/96$, where a Painlevé-II function controls the profile. If correct, the result shows that arbitrarily high-order rational solitons of this integrable equation share a universal near-field shape, analogous to known large-order universality for NLS rogue waves.","feed_headline":"High-order rational solitons converge to one master profile","feed_subtitle":"After dividing by order, complex mKdV solitons approach one universal shape tied to Painlevé-III equations.","key_machinery":"The load-bearing object is a solvable $2\\times2$ matrix Riemann–Hilbert problem: find an analytic matrix $N^\\pm(\\Lambda;X,T)$ normalized to $I$ at infinity whose boundary values on the unit circle jump by $e^{-i(\\Lambda X+4\\Lambda^3T)\\sigma_3}Q e^{\\mp 2i\\Lambda^{-1}\\sigma_3}Q^{-1}e^{i(\\Lambda X+4\\Lambda^3T)\\sigma_3}$. The argument proceeds by writing the multi-rational soliton as a transformed RHP (RHP 2), rescaling $X=nx$, $T=n^3t$, $\\Lambda=\\lambda/n$, and showing the jump tends to the model jump; the vanishing lemma gives existence and uniqueness of $N^\\pm$, and $q$ is recovered as $2i\\lim_{\\Lambda\\to\\infty}\\Lambda N^\\pm_{12}(\\Lambda;X,T)$. A second transformation $H=D e^{-i(\\Lambda X+4\\Lambda^3T+2\\Lambda^{-1})\\sigma_3}$ converts the jump to the constant matrix $Q$, yielding a Lax pair in $X$ and $\\Lambda$ whose compatibility produces the ODEs; at $T=0$ this is the Lax pair of the first member of the Painlevé-III hierarchy. The asymptotic analysis uses nonlinear steepest descent with parabolic-cylinder and Painlevé-II model functions as local parametrices.","core_discovery":"The paper claims that the $k$th-order rational soliton $q_k$ of the focusing c-mKdV equation with unit background obeys, for $k=2n$ and $k=2n-1$ respectively, $\\frac1n q_{2n}(X/n,T/n^3)=\\hat q_+(X,T)+O(1/n)$ and $\\frac1n q_{2n-1}(X/n,T/n^3)=\\hat q_-(X,T)+O(1/n)$ uniformly on compact subsets of $\\mathbb R^2$, where $\\hat q_\\pm$ is reconstructed from a uniquely solvable model Riemann–Hilbert problem and satisfies the c-mKdV equation in $(X,T)$. The same limit functions satisfy two ordinary differential equations in $X$; at $T=0$, eliminating $\\hat q$ gives Eq. (1.27), the first member of the Painlevé-III hierarchy. Theorems 1.5 and 1.6 provide the large-$X$ asymptotic behavior: for $-1/96<a\\le 0$, $\\hat q_+(X,aX^2)\\sim \\sqrt{2p}\\cos(\\varphi(X,a))\\,X^{-3/4}/\\sqrt{6ab+b^{-3}}$, and as $a\\to -1/96$ with $a+1/96=O(X^{-1/3})$ the profile is $O(X^{-2/3})$ and is expressed through the Painlevé-II function $V_1$. The paper additionally proves $\\hat q_+(X,T)=-\\hat q_-(X,T)$, that both limits are real, and that they inherit the symmetry $\\hat q_\\pm(-X,-T)=\\hat q_\\pm(X,T)$.","pith_inferences":["A natural extension beyond the paper's claims is to push the small-norm expansion one order further: the jump in Lemma 2.1 is $I+O(1/n)$, so the next correction to $\\hat q_\\pm$ should be computable and would give an $O(1/n^2)$ refinement testable against the closed-form formula.","The self-similar master profile suggests that high-order rational solitons of c-mKdV, when generated in short-pulse optical systems with higher-order dispersion, should appear as a single rescaled waveform; this is an experimental consequence the paper does not assert.","The partial large-$T$ results indicate that the same steepest-descent machinery, with new local models at the four endpoints, should yield full Painlevé-type asymptotics for $T\\to\\infty$, giving a concrete next computation.","Because the model RHP has the same structural role as the one used for NLS large-order rogue waves, higher members of the Painlevé-III hierarchy may appear when the c-mKdV hierarchy is considered; the paper's method should transfer to those cases."],"forward_implications":["If the main theorem is correct, the $n$-scaled rational solitons $(1/n)q_k$ have a universal near-field limit $\\hat q_\\pm$ that no longer depends on the detailed determinantal formula used to generate them.","The limit profile solves the c-mKdV equation in the rescaled variables, so the large-order family is not an isolated approximation but an exact solution of the same dispersive PDE.","At $T=0$, the spatial profile is constrained by the first Painlevé-III hierarchy member, giving an explicit ODE classification of the near-field shape.","The explicit large-$X$ formulas in Theorems 1.5 and 1.6 and Corollary 4.1 turn the universal profile into concrete oscillatory decay laws with computable phase and amplitude, including a Painlevé-II transition when $a\\to -1/96$.","The paper states that the same RHP-based route extends to higher-order c-mKdV equations, the mKdV hierarchy, and the $(2+1)$-dimensional KP equation."],"supporting_citations":[{"why":"Supplies the closed-form determinantal $(2n-1,2n)$th-order rational soliton formula (1.19) that the RHP construction is claimed to reproduce.","marker":"[21]"},{"why":"Provides the robust inverse-scattering/RHP method used to build the solvable RHP for the c-mKdV equation with nonzero background.","marker":"[10]"},{"why":"Provides the vanishing lemma used to prove existence and uniqueness of the solution to the model RHP 3.","marker":"[96]"},{"why":"Provides the nonlinear steepest descent method used for the large-X and transitional asymptotic expansions.","marker":"[26]"},{"why":"Provides the conformal-mapping technique used near the double critical point $a=-1/96$ to set up the transitional parametrix.","marker":"[22]"},{"why":"Supplies the Painlevé-II model function used for the transitional asymptotics.","marker":"[66]"}],"fun_headline_variants":["Soliton cascade collapses to one Painlevé-III master profile","All high-order c-mKdV solitons share a single limit shape","Rational solitons converge to a universal Painlevé-III curve","Large-order solitons unify into one explicit profile"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"That the Riemann–Hilbert problem built in Proposition 2.3 reproduces exactly the same family of multi-rational solitons as the determinantal formula (1.19) imported from reference [21]; the large-order theorems are proved for the RHP-generated family, and any mismatch at finite order would change the object being approximated.","fun_headline_variants_meta":{"raw":{"variants":["Soliton cascade collapses to one Painlevé-III master profile","All high-order c-mKdV solitons share a single limit shape","Rational solitons converge to a universal Painlevé-III curve","Large-order solitons unify into one explicit profile"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000907,"raw_usage":{"total_tokens":4011,"prompt_tokens":1165,"completion_tokens":2846,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":781,"completion_tokens_details":{"reasoning_tokens":2780}},"tokens_in":781,"tokens_out":2846,"duration_ms":16999,"temperature":1.0,"reasoning_tokens":2780,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:49:56.996567+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the closed-form formula (1.19) for $n=10,20,40$ at a fixed compact set of $(X,T)$, compute $(1/n)q_k(X/n,T/n^3)$, and compare with a numerical solution of the model RHP 3: the difference should shrink like $O(1/n)$. A persistent nonzero difference would show that RHP 2 and formula (1.19) define different solution families, invalidating the application of the theorems to the explicit solitons.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the closed-form determinantal $(2n-1,2n)$th-order rational soliton formula (1.19) that the RHP construction is claimed to reproduce."},{"cited_title":"Zhou, The Riemann-Hilbert problem and inverse scatt ering, SIAM J","cited_arxiv_id":null,"evidence_quote":"Provides the vanishing lemma used to prove existence and uniqueness of the solution to the model RHP 3."},{"cited_title":"Deift, X","cited_arxiv_id":null,"evidence_quote":"Provides the nonlinear steepest descent method used for the large-X and transitional asymptotic expansions."},{"cited_title":"Chester, B","cited_arxiv_id":null,"evidence_quote":"Provides the conformal-mapping technique used near the double critical point $a=-1/96$ to set up the transitional parametrix."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Painlevé-II model function used for the transitional asymptotics."}],"review_version":1}