{"id":"8280e38a-98b9-401c-b919-d34242e16a12","arxiv_id":"2505.06768","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Line solitons of the 2D Toda lattice are shown to be linearly stable in an exponentially weighted space, with perturbations converging to a damped-wave profile.","lead":"This paper proves that the moving straight-soliton solutions of the two-dimensional Toda lattice are linearly stable: small perturbations decay exponentially in a suitable weighted space, once a few symmetry-related secular modes are removed. The result matters because the 2D Toda lattice is a canonical integrable discrete wave equation, and this is the first transverse stability proof for its line solitons.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The α2-endpoint estimate (102) for the mid-frequency band is not justified for arbitrary data: (97) requires secular conditions that P0(η1,η2)U f does not satisfy, so the interpolation step (99)+(102)→(103) lacks a valid endpoint.","rationale":"The reader's weakest-assumption flag correctly identifies the interpolation step between (99)/(102) and (103) as the most load-bearing point in the proof of Theorem 1.1. My stress-test sharpens this: the interpolation machinery itself (Riesz–Thorin on the scale ℓ^2_αH^1×ℓ^2_αL^2) is standard and would work if both endpoint estimates were valid for the same operator on the respective endpoint spaces. The real problem is that (102), the α2-endpoint, is asserted for P0(η1,η2)U(t,s)f while the only stated justification, (97), requires the secular orthogonality conditions (78) and (84), which the Fourier restriction P0 does not impose. This is not a mere typo: without an additional decomposition and an estimate for the secular projection in the [η1,η2] band, the mid-frequency band is uncontrolled in the target weight α. I considered other possible concerns — the apparent missing e^{αct} factor in Lemma 2.4, the α∈(0,κ) restriction in Lemma 5.3, and the typo in Lemma 5.9's kernel vector — but these are explicit errors that are easy to repair and do not threaten the overall strategy. The (102) gap is more substantial because it sits at the junction of the high-frequency, low-frequency, and interpolation arguments; if it cannot be filled, Theorem 1.1 does not follow. Because the needed estimate is plausibly obtainable by projecting and using the explicit temporal decay of g±, the appropriate verdict remains CONDITIONAL, as the reader already concluded.","tokens_in":28829,"tokens_out":28577,"duration_ms":281928,"concrete_test":"Check the derivation of (102) by taking an initial datum supported only on frequencies η∈(η1,η2) with a nonzero component of g+(η0) for some η0∈(η1,η2), and verify whether the resulting Q'_η satisfies (78) and (84). Constructively, replace P0(η1,η2) in (102) by P0(η1,η2)(I−P1(t,η2)) and add the missing bound for {P1(t,η2)−P1(t,η1)}U(t,s)f using formulas (88)–(91) and Lemma 2.6; check whether this yields a decay rate δ_R(η1)>0 in the e^{-α2ct} scale. Only if both pieces satisfy the claimed bound does (102) hold for arbitrary f, and only then can interpolation produce (103).","verdict_should_be":"UNCHANGED","load_bearing_attack":"In the proof of Theorem 1.1, the band [η1,η2] is controlled by applying complex interpolation to (99) and (102) to obtain (103). The α1-endpoint (99) is obtained from the high-frequency estimate (95), which holds for all data. The α2-endpoint (102), however, is claimed to follow from (97), which is proved only under the secular orthogonality conditions (78) and (84). The operator P0(η1,η2)U(t,s)f does not enforce those conditions: for η∈(η1,η2), its Fourier components include the secular modes g±(η), which belong to ℓ^2_{α2} because η2<η*(α2). Thus (102) is not a direct consequence of (97), and the sentence \"Similarly, it follows from (97) and Lemma 2.6 that (102)\" is not justified as written. To obtain (102) one would need to decompose P0(η1,η2) as P0(η1,η2)(I−P1(t,η2)) on the band plus the projection remainder P1(t,η2)−P1(t,η1), estimating the first term via (97) and the second via the explicit decay of g1,g2 in the e^{-α2ct} scale; the latter estimate is absent. Since (103) is essential for the mid-frequency band, this is a genuine gap in the proof of Theorem 1.1. The gap is concrete but likely repairable, so the correct disposition remains conditional.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the transverse linear stability of 1-line solitons for the two-dimensional Toda lattice. The linearized equation around the soliton is rewritten in terms of the variable Q', and the author constructs Jost functions and dual Jost functions that generate the secular modes of the linearized problem. Using Darboux transformations, the linearized problem around the soliton is related to the linearized problem around the zero solution, for which exponential decay in exponentially weighted spaces is proved by a Fourier argument. Theorem 1.1 claims that, under a secular-term orthogonality condition, every perturbation decays exponentially in a weighted space with weight e^{αn}, α∈(0,2κ). Theorem 1.2 claims that the dominant part of the solution is described by a damped one-dimensional wave equation in the transverse variable, with an explicit convolution profile and an O(t^{-1/4}) error.","tokens_in":29115,"tokens_out":7055,"duration_ms":67140,"significance":"If the proof is completed, the paper gives the first linear stability result for line solitons of the 2D Toda lattice for solitons of any size, with explicit constants and decay rates and no fitted parameters. The Darboux-transformation approach is well suited to the integrable structure and the paper contains detailed spectral and Jost-function estimates. The central idea of removing the secular modes before applying the Darboux correspondence is natural and, if the gaps described below are fixed, would be an important step toward nonlinear stability of 1-line solitons. However, the manuscript as written has three load-bearing gaps: the mid-frequency band estimate (102) is not justified for arbitrary data, the low-frequency Darboux lemmas are stated only for α<κ while the main theorem claims α<2κ, and the proof of Theorem 1.2 stops before deriving the claimed asymptotic profile.","major_comments":[{"comment":"The estimate for the mid-frequency band is not derived for arbitrary data. The sentence 'Similarly, it follows from (97) and Lemma 2.6 that (102)' is not justified, because (97) is proved only for solutions Q'_η of (34) that satisfy the secular orthogonality conditions (78) and (84). The Fourier-localized component P0(η1,η2)U(t,s)f does not satisfy those conditions. Since (103) is obtained by complex interpolation between (99) and (102), the control of the band [η1,η2] in the target weight is not established as written. A repair would require decomposing P0(η1,η2) into P0(η1,η2)(I−P1(t,η2)) on the band plus P1(t,η2)−P1(t,η1), estimating the first term via (97) and the second via the explicit decay of g1 and g2 in the e^{−α2ct} scale; the latter estimate is absent.","section":"Section 6, Eq. (102)"},{"comment":"Lemma 5.3 is stated only for α∈(0,κ), but Lemma 5.6, and through it the low-frequency estimates (96) and (97), are used for weights up to 2κ. In the proof of Theorem 1.1 one chooses α1<α<α2<2κ, which may give α2≥κ; the Darboux-inversion argument in Lemma 5.6 relies on Lemma 5.3 and is therefore not proved in that range. Either Lemma 5.3 must be extended to α∈(0,2κ), or Theorem 1.1 must be restricted to α<κ. As written, the case α∈[κ,2κ) is not covered by the proof.","section":"Section 5, Lemmas 5.3 and 5.6"},{"comment":"The proof of Theorem 1.2 ends immediately after the assertion that f_j∈L^1(R_y). It does not define the function f, does not derive the convolution profile H_t*W_t*f, and does not prove the claimed O(t^{-1/4}) bound for the full difference in the ℓ^2_αH^1×ℓ^2_αL^2 norm. In addition, the displayed line 'f̂(η) + +O(η)' contains a typo, and the operator P(t,η0) appearing in the first line of the proof is not defined (it should presumably be P1(t,η0)). The asymptotic statement is therefore not proven as written.","section":"Section 6, proof of Theorem 1.2"}],"minor_comments":[{"comment":"'equaitons' should be 'equations'.","section":"Section 1, Remark 1.2"},{"comment":"'Daroboux transformations' should be 'Darboux transformations'.","section":"Section 3 heading"},{"comment":"'respectivey' should be 'respectively'.","section":"Lemma 5.2"},{"comment":"The complex interpolation step should explicitly identify the interpolation couple, for example (ℓ^2_{α1}H^1×ℓ^2_{α1}L^2, ℓ^2_{α2}H^1×ℓ^2_{α2}L^2), and state why the operator P0(η1,η2)U(t,s) interpolates between the two endpoint estimates; this is currently only asserted.","section":"Around Eq. (103)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a serious contribution and the main theorem is likely correct, but the proof as written has three technical gaps that need to be addressed before publication: the interpolation endpoint for the mid-frequency band, the range of α in the low-frequency Darboux lemmas, and the incomplete proof of Theorem 1.2. These are repairable within the manuscript's framework, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this is a genuine new result—transverse linear stability of 1-line solitons for 2D Toda for all amplitudes, with a damped-wave asymptotic description—and the Darboux machinery is the right tool. I believe the main theorem is true, but the written proof has a real gap in the mid-frequency band, and Theorem 1.2 is sketched rather than proved.\n\nWhat's new and good: the result itself is not in any cited paper; the constants are explicit; the reduction from the linearized 1-soliton equation to the free equation via Darboux transformations is carefully set up; Lemmas 5.3–5.9 and the secular-mode orthogonality are substantial and mostly checkable. The paper also states clearly that the linearized equation is non-autonomous and that slowly decaying modes have to be projected out. That is honest and correct.\n\nSoft spots, in order of seriousness.\n\n1. The stress-test note lands. Equation (102) is claimed to follow from (97), but (97) is only proved under the secular orthogonality conditions (78),(84). The operator P0(η1,η2)U(t,s)f does not enforce those conditions. For each η in the band, the projection contains components of g±(η), which are nonzero in ℓ²_{α2} and are precisely the secular modes. So the estimate (102) is not justified as written. You need a decomposition of P0(η1,η2) into P0(η1,η2)(I−P1(t,η2)) plus a remainder, plus an explicit bound on the remainder in the e^{-α2 ct} scale. That bound is absent. This is a genuine gap in the proof of Theorem 1.1, but it looks repairable with the paper's own tools.\n\n2. Theorem 1.2's proof stops after defining f_j and asserting f_j∈L^1. The convolution profile and the O(t^{-1/4}) bound for the full solution are asserted rather than derived. This is an omission, not obviously an error; a referee should ask for the missing page.\n\n3. Minor: Remark 1.2 mixes the nκ−t sinhκ and nκ+t sinhκ conventions. Typo-level.\n\nVerdict: conditional. The central argument is coherent, the gap is concrete but likely fixable, and the result deserves referee time. Send it to a serious referee, with instructions to check the interpolation endpoint and demand the missing derivation in Theorem 1.2. I would not desk-reject.","headline":"New result with a repairable mid-frequency gap: it deserves refereeing, not rejection.","tokens_in":29693,"tokens_out":2739,"would_cite":true,"duration_ms":27521,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B35","37K40","35Q51"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves exponential linear stability of 1-line solitons of the 2D Toda lattice in a weighted space, for any soliton size, after removing the slowly decaying secular modes.","keywords":["2D Toda","line solitary waves","transverse linear stability","Darboux transformation","Bäcklund transformation","exponentially weighted space","secular modes","KP-II equation"],"falsifier":"Compute the projected evolution $P_0(\\eta_1,\\eta_2)U(t,s)$ on the mid-frequency band and test whether the bound (103) holds with a decay rate independent of $\\alpha$ between $\\alpha_1$ and $\\alpha_2$; a counterexample, such as initial data in this band whose $\\ell^2_\\alpha$-weighted norm decays only polynomially, would disprove Theorem 1.1. The specific place to look is the complex interpolation step between (99) and (102), where the paper does not verify that the projectors, the solution operator, and the decay estimates interpolate to the target weight.","tokens_in":28567,"feed_emoji":"🌊","tokens_out":9679,"duration_ms":86232,"temperature":0.7,"pith_summary":"The paper aims to prove that the one-dimensional line soliton of the two-dimensional Toda lattice—the integrable semi-discrete wave equation whose continuous limit is KP-II—is linearly stable against transverse perturbations, and to identify the asymptotic mechanism of that stability. The main theorem states that any solution of the linearized equation whose initial data are orthogonal to the slowly decaying secular modes decays exponentially in an exponentially weighted energy space, uniformly for solitons of every amplitude (every $\\kappa>0$). The proof works by a Darboux transformation that maps the linearized problem around the soliton onto the linearized problem around the zero solution, where decay is elementary. The paper further shows that the surviving part of a generic perturbation is a modulation of the time derivative of the soliton, with the modulation governed by a damped wave equation in the transverse variable, mirroring the known picture for KP-II line solitons. This linear result is positioned as the key step toward a full nonlinear stability theorem for 2D Toda line solitons.","feed_headline":"2D Toda line solitons are linearly stable at any size","feed_subtitle":"Darboux transformations show perturbations decay exponentially, leaving a damped-wave modulation.","key_machinery":"The load-bearing object is the Darboux transformation, expressed through Jost functions and dual Jost functions for the Lax pair of the 2D Toda equation. For the zero background the Jost functions are $\\Phi^0_n(\\beta)=\\beta^n e^{\\beta x-s/\\beta}$; for the one-soliton background they are modified by the $\\tau$-function of the soliton. Products $\\Phi(\\beta_1)\\Phi^*(\\beta_2)$ solve the linearized equation, and the transformation (41) connects solutions $Q'$ of the linearized equation around the soliton to solutions $Q$ of the linearized equation around zero. The secular modes $g_\\pm(\\eta)$, $g_{\\pm,*}(\\eta)$ are built from these Jost functions; the orthogonality condition in Theorem 1.1 removes exactly the slowly decaying modes that would otherwise spoil uniform exponential decay. The exponential weight $e^{2\\alpha n}$ in the lattice direction is what makes the soliton outrun the perturbations, the idea drawn from weighted-space stability theory for solitary waves.","core_discovery":"On its own terms, the paper's central claim is Theorem 1.1: fix $\\kappa>0$, $\\alpha\\in(0,2\\kappa)$, $c=\\sinh\\kappa/\\kappa$ and $\\eta_0\\in(0,\\eta_*(\\alpha))$, where $\\eta_*(\\alpha)=\\tanh(\\kappa+\\alpha)\\sqrt{\\sinh\\alpha\\,\\sinh(2\\kappa+\\alpha)}$. If $R'$ is a solution of the linearized equation (9) in $C(\\mathbb{R};\\ell^2_\\alpha H^1)\\cap C^1(\\mathbb{R};\\ell^2_\\alpha L^2)$ and if at time $t_0$ its weighted pairing with each secular mode $\\tilde{g}^{\\pm,*}(\\eta)$ vanishes for all $|\\eta|\\le\\eta_0$, then for every $t\\ge t_0$, $e^{-\\alpha ct}(\\|R'(t)\\|_{\\ell^2_\\alpha H^1}+\\|\\partial_t R'(t)\\|_{\\ell^2_\\alpha L^2})$ is bounded by $K e^{-b(t-t_0)}e^{-\\alpha ct_0}$ times the same norm at $t_0$. In words: once the slowly decaying secular component is removed, every perturbation decays exponentially in the weighted norm, at a rate independent of the initial time. Theorem 1.2 sharpens this by identifying the asymptotic profile: the dominant part is a convolution of the initial amplitude with a heat kernel $H_t$ and a finite-speed wave kernel $W_t$, multiplied by $(\\partial_t R^\\kappa,\\partial_t^2 R^\\kappa)$, with error $O(t^{-1/4})$.","pith_inferences":["If the interpolation step flagged below can be justified, the same Darboux machinery should give linear stability for multi-line solitons by iterating the transformation, with the secular condition becoming a finite set of orthogonality constraints per line.","Because the secular modes are indexed by a continuous parameter $\\eta$, the linearized equation has an infinite-dimensional critical manifold; the theorem suggests the 1-line soliton is 1-codimensionally stable in a critical space, and a nonlinear analogue would track a phase shift in both the lattice and transverse directions.","A direct numerical check of Theorem 1.2 is feasible: initialize (9) with a compactly supported transverse profile orthogonal to $g_{\\pm,*}$, and compare the solution at large $t$ with the predicted $(H_t*W_t*f)(y)\\partial_t R^\\kappa$; agreement to $O(t^{-1/4})$ would confirm the damped-wave description.","The $O(t^{-1/4})$ rate comes from the $L^2$ norm of $e^{-t\\lambda_2\\eta^2}$, so one could test whether the rate improves to $t^{-1/2}$ in $L^1$ or in sup-norm-type weights, an extension the paper does not address."],"forward_implications":["Any perturbation of a 1-line soliton that is orthogonal to the secular modes decays exponentially in the weighted energy norm; the line soliton is therefore linearly stable in $\\ell^2_\\alpha H^1\\times\\ell^2_\\alpha L^2$.","For data with additional $\\ell^2_\\alpha L^1$ regularity, the asymptotic shape is explicit: the amplitude is $(H_t*W_t*f)(y)$, so the transverse spreading is diffusive at rate $t^{-1/4}$ in suitable norms, just as for KP-II.","The stability statement holds uniformly for solitons of any size $\\kappa>0$, not only small-amplitude ones.","In the opposite weight $\\ell^2_{-\\alpha}$, the same theorem states that the line soliton is linearly stable, which is the weight relevant for perturbations behind the soliton.","The linear decay estimate is the natural input for a B\\\"acklund-based nonlinear stability argument, the route the paper identifies for proving nonlinear stability of 2D Toda line solitons."],"supporting_citations":[{"why":"Supplies the $\\tau$-function bilinear formalism and the B\\\"acklund transformations that the Darboux maps are built from.","marker":"[15]"},{"why":"Establishes the Darboux correspondence for linearized Toda solitons; Lemma 5.2 of the present paper follows its Lemma 5.","marker":"[31]"},{"why":"Introduces the exponentially weighted space strategy that makes the soliton outrun perturbations, which the paper adapts to 2D Toda.","marker":"[36]"},{"why":"Provides the analogous weighted-space linear stability result for KP-II line solitons that the present theorem extends to the discrete 2D Toda setting.","marker":"[27]"},{"why":"Uses Darboux transformations to prove asymptotic stability of Toda m-solitons, the methodological antecedent for the linearized argument.","marker":"[4]"},{"why":"Gives the canonical/B\\\"acklund transformation for the one-dimensional Toda lattice that underlies the modified-Toda equations (37)-(40).","marker":"[42]"}],"fun_headline_variants":["2D Toda line solitons stable at any size","Exponential decay confirms 2D Toda soliton stability","Darboux transform locks line solitons in 2D Toda","Perturbations fade: 2D Toda line solitons stable","Any-size line solitons proven stable in 2D Toda"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The decay estimate for the middle range of transverse frequencies is obtained by interpolating two bounds that hold in different exponential weights, and the proof does not verify that the projected solution operator and its decay rate interpolate to the target weight.","fun_headline_variants_meta":{"raw":{"variants":["2D Toda line solitons stable at any size","Exponential decay confirms 2D Toda soliton stability","Darboux transform locks line solitons in 2D Toda","Perturbations fade: 2D Toda line solitons stable","Any-size line solitons proven stable in 2D Toda"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000636,"raw_usage":{"total_tokens":2967,"prompt_tokens":1013,"completion_tokens":1954,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":629,"completion_tokens_details":{"reasoning_tokens":1859}},"tokens_in":629,"tokens_out":1954,"duration_ms":15339,"temperature":1.0,"reasoning_tokens":1859,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:34:42.883408+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the projected evolution $P_0(\\eta_1,\\eta_2)U(t,s)$ on the mid-frequency band and test whether the bound (103) holds with a decay rate independent of $\\alpha$ between $\\alpha_1$ and $\\alpha_2$; a counterexample, such as initial data in this band whose $\\ell^2_\\alpha$-weighted norm decays only polynomially, would disprove Theorem 1.1. The specific place to look is the complex interpolation step between (99) and (102), where the paper does not verify that the projectors, the solution operator, and the decay estimates interpolate to the target weight.","supporting_citations":[{"cited_title":"Hirota , The direct method in soliton theory","cited_arxiv_id":null,"evidence_quote":"Supplies the $\\tau$-function bilinear formalism and the B\\\"acklund transformations that the Darboux maps are built from."},{"cited_title":"Mizumachi and R","cited_arxiv_id":null,"evidence_quote":"Establishes the Darboux correspondence for linearized Toda solitons; Lemma 5.2 of the present paper follows its Lemma 5."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the exponentially weighted space strategy that makes the soliton outrun perturbations, which the paper adapts to 2D Toda."},{"cited_title":"Mizumachi , Stability of line solitons for the KP-II equation in ^2 , Mem","cited_arxiv_id":null,"evidence_quote":"Provides the analogous weighted-space linear stability result for KP-II line solitons that the present theorem extends to the discrete 2D Toda setting."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Uses Darboux transformations to prove asymptotic stability of Toda m-solitons, the methodological antecedent for the linearized argument."},{"cited_title":"Toda and M","cited_arxiv_id":null,"evidence_quote":"Gives the canonical/B\\\"acklund transformation for the one-dimensional Toda lattice that underlies the modified-Toda equations (37)-(40)."}],"review_version":1}