{"id":"778a1240-d11f-4c75-a04b-0e8a5831e3c0","arxiv_id":"2509.08271","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The 2D cubic Klein-Gordon equation is approximated by the cubic NLS up to O(ε^2) error over times of order ε^{-2/(N+1)} with a (1+t)^N prefactor.","lead":"For the two-dimensional cubic Klein-Gordon equation, this paper proves that solutions converge to solutions of the cubic nonlinear Schrödinger equation as the speed of light tends to infinity, with error scaling like the square of the small parameter. The result justifies using the Schrödinger approximation in 2D numerical computations and gives explicit, if not sharp, long-time error bounds.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central error estimate rests on WKB structural lemmas (Prop. 3.4/3.5/Cor. 3.6) deferred to [7] and not verified for the 2D system; if they fail, the approximate solution U_a does not exist and Theorem 1.3 collapses.","rationale":"The paper's central claim is a long-time error estimate for the 2D cubic Klein-Gordon equation in the non-relativistic limit. The proof structure is: construct a WKB approximate solution U_a to arbitrary order, then prove nonlinear stability of U_a on a long time interval. The stability part (Section 4) is largely self-contained: the Duhamel formula, the Gronwall/bootstrap estimates, and the NLS decay estimates are written out in detail. The genuinely external input is the arbitrary-order WKB construction. Propositions 3.4, 3.5, and Corollary 3.6 are the inductive engine that produces the profile equations and the polynomial forms of all higher-order amplitudes; without them, U_a is only computed at low orders, and there is no way to bound the residual R_ε in Theorem 2.2. The manuscript explicitly defers these propositions to [7], but [7] is described as a three-dimensional result, and the current paper's state space has a different dimension (4 components instead of 5). This makes the reader's weakest assumption genuinely load-bearing: if the deferred identities do not carry over to 2D, the main theorems have no approximate solution to stabilize. I did not find a separate internal inconsistency that would be more serious. Minor issues, such as the Duhamel formula in (4.3) writing S(t)r_ε instead of εS(t)r_ε, are conservative and do not threaten the argument. Therefore the appropriate verdict remains CONDITIONAL, exactly as the reader recommended; the condition is that the deferred WKB structural lemmas are correct and applicable in 2D.","tokens_in":27450,"tokens_out":22094,"duration_ms":225907,"concrete_test":"Verify the deferred WKB construction directly for the 2D system. Concretely: (i) Starting from the displayed U0, U1, U2, U3 in (3.7), (3.13), (3.24), (3.36), use the induction formulas in Proposition 3.5 to compute U_{4,1}, U_{4,3}, U_{4,5} and check that they match the expressions obtained from (3.39)–(3.45); (ii) re-derive Proposition 3.4 by applying Π1 to the Φ_{n,1} equation for general n and confirm that (3.47) is obtained without any hidden dependence on spatial dimension beyond Δ_x; (iii) compare the proof of [7, Section 3] line-by-line and confirm that the change from R^3 to R^2 does not alter the harmonic sets H_n, the parity argument in Proposition 3.3, or the recursion (3.48). If these checks pass, the reader's concern is resolved; if any step fails, Theorems 2.2, 2.4, and 1.3 are not established as written.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The proof of Theorem 1.3 goes through the stability estimate Theorem 2.4, which presupposes the global WKB approximate solution U_a from Theorem 2.2. That U_a is built by induction, and the inductive backbone is Proposition 3.4 (form of U_{n,±1} and the profile equations (3.47)), Proposition 3.5 (form of U_{n,p} for |p|≥3 and the recursion (3.48)), and Corollary 3.6 (polynomial dependence). These are stated without proof and deferred to [7]. The 2D system has state vector U=(ε∇u, ε^2∂_t u, u) with ∇u in R^2, so the matrices A(∂x), A0, L_p, Π1 in Sections 2–3 are 4×4, whereas in the 3D paper [7] they are 5×5. The scalar profile equations may be identical, but that is an assertion rather than a demonstrated fact. All later regularity estimates in Section 3.8 and the stability analysis in Section 4 assume the exact algebraic form of U_a. If any of the structural identities in Propositions 3.4/3.5 or Corollary 3.6 fail in 2D, the approximate solution does not satisfy (2.7) to the claimed order, the residual bound in Theorem 2.2 is unavailable, and the ε^{K+1} stability estimate in Theorem 2.4 has no object to which it applies. This is load-bearing because it is the only place where arbitrary-order WKB existence is established; the rest of the paper controls perturbations of an approximate solution whose existence is taken as input.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the non-relativistic limit ε→0 for the two-dimensional cubic Klein-Gordon equation (1.1)–(1.2). It rewrites the equation as the symmetric hyperbolic system (2.2) and constructs WKB approximate solutions of arbitrary even order. The main results are stability estimates for these approximate solutions in H^{s-2K-4} over long time intervals, yielding error O(ε^{K+1}) with polynomial-in-time prefactors. Theorems 1.3 and 1.4 give the convergence to the cubic nonlinear Schrödinger equation with rates O(ε^2) and O(ε), with time spans of size O(ε^{-2/(N+1)}) and O(ε^{-α/2}), respectively. The proof extends the three-dimensional analysis in [7] to two dimensions; the new ingredients are the defocusing long-time stability argument and the translation of the WKB construction to the 4×4 first-order system. The exposition is generally clear, but several load-bearing structural propositions are deferred to [7] without a self-contained verification for the 2D system.","tokens_in":27889,"tokens_out":7042,"duration_ms":84825,"significance":"If the results are correct, they improve on the earlier L^2 bound of Wu and Lei by providing H^s error estimates with explicit algebraic-in-time prefactors, and they are consistent with numerical observations. The stability framework and the careful tracking of initial-data-dependent growth exponents are valuable. The paper also makes a sharp distinction between high- and low-regularity data and between focusing and defocusing cases. However, the core WKB existence result is not established inside the paper: Propositions 3.4, 3.5 and Corollary 3.6 are stated without proofs and are deferred to [7], a paper that treats the three-dimensional case. Since all subsequent error estimates are conditional on the exact algebraic form of the approximate solution, the current manuscript is not self-contained on this point.","major_comments":[{"comment":"The derivation of U_{4,3} and U_{4,5} is omitted with the note that the calculation is essentially the same as in [7]. These formulas feed into Corollary 3.6 and the regularity estimates, so they are part of the same gap. While it is reasonable to cite a prior paper for a routine computation, the dimension change makes it necessary to at least display the verification for the 2D case, or to state clearly that [7] already contains the 2D version. Without this, the induction is not independently checkable.","section":"§3.6, equations (3.44)–(3.45)"}],"minor_comments":[{"comment":"The deductions of Theorems 1.3–1.6 from Theorems 2.3–2.5 are said to be 'straightforward' and are omitted. Theorem 1.3 is shown explicitly, but a short paragraph for Theorems 1.4 and 1.6 would improve readability, particularly because the approximate solution u_a in Theorem 1.5 contains ε^K u_K and ε^{K+2}u_{K+2} and one must check that these extra terms are absorbed in the error.","section":"§4.4"},{"comment":"There are typos in the affiliation line: 'school of Mathmatics' and 'Sistrict' should be 'School of Mathematics' and 'District'.","section":"Title page and affiliation"},{"comment":"The last equality in (3.22) uses the Schrödinger equation for g0, but this is not stated at that point. Adding a parenthetical remark would help the reader.","section":"§3.3, equation (3.22)"}],"recommendation":"major_revision","confidential_remarks":"The heavy reliance on [7], a paper co-authored by the same first author, is also a novelty-disclosure issue: the genuinely new content appears to be the 2D adaptation and the stability estimates. I would ask for a self-contained proof of the WKB structural propositions or an explicit transfer theorem, rather than accepting the current level of deferral. The editor may also wish to consider whether, given the close relation to [7], the contribution is best framed as an extension paper with the deferred arguments included in an appendix."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new content here is real: 2D cubic Klein-Gordon with H^s error estimates of order C(1+t)^N ε^2 over times of order ε^{-2/(N+1)}, plus an arbitrary-order WKB approximation. That is not in the earlier 3D paper [7] (linear prefactor) nor in Wu-Lei's 2D L^2 result [22]. The defocusing-case long-time bounds and the explicit polynomial-in-time prefactors are the genuine contribution. The proof strategy—reformulate as a symmetric hyperbolic system, build a WKB approximate solution, then run a stability argument against the cubic NLS decay—is sensible and carried out in reasonable detail for the stability part.\n\nThe main soft spot is exactly what the stress test flags: Propositions 3.4, 3.5, and Corollary 3.6, which give the algebraic form of the profiles U_{n,p} and the recursion relations, are stated without proof and deferred to [7]. These lemmas are load-bearing: the approximate solution U_a is constructed by induction through them, and the residual estimates and stability theorems presuppose that construction. The 2D system has 4x4 matrices, whereas [7] works with 5x5; the scalar profile equations are likely the same, but the paper does not demonstrate the transfer. If any structural identity fails in 2D, Theorem 2.2 has no object to control and Theorem 1.3 collapses. That is a genuine gap in self-containedness, not a manufactured one.\n\nThat said, the gap is plausibly fillable. The computations in Sections 3.1–3.5 explicitly work out orders up to O(ε^2) in 2D, matching the pattern of [7], and the deferred lemmas look like the same induction in a different ambient dimension. The paper also has some minor missing details, like the 'straightforward' corollary proofs in Section 4.4, but those are genuinely routine. No fatal internal inconsistency. I disagree with any claim that the self-citation is itself disqualifying—[7] is a published, peer-reviewed paper, and deferring to it is legitimate, provided the transfer to 2D is made explicit.\n\nWho is this for? Researchers working on relativistic-to-nonrelativistic limits, especially numerical analysts who need rigorous error bounds matching observed (1+t)^N behavior. The paper deserves serious refereeing, but the referee should verify—or the authors should provide—the 2D versions of Propositions 3.4 and 3.5. I would recommend peer review with a request for revision: either prove the structural lemmas in an appendix or give a precise reduction to [7] showing the dimension change leaves the algebra unchanged.","headline":"New 2D result that extends the 3D non-relativistic-limit analysis, but it leans on deferred structural lemmas from the authors' prior paper and needs a revision to be fully convincing.","tokens_in":28412,"tokens_out":1060,"would_cite":true,"duration_ms":14181,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","35L05","35B25"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that the two-dimensional cubic Klein-Gordon equation is approximated by the cubic nonlinear Schrödinger equation with error O(ε²) up to times ε^{-2/(N+1)}.","keywords":["non-relativistic limit","cubic Klein-Gordon equation","cubic nonlinear Schrödinger equation","WKB expansion","geometric optics","error estimates","long-time stability","two-dimensional"],"falsifier":"Compute $u$ numerically for the defocusing 2D cubic Klein-Gordon equation with $\\varepsilon = 10^{-2}$ and smooth localized data, solving the cubic NLS for $g_0$ simultaneously, and measure $\\|u - (e^{i\\theta}g_0 + e^{-i\\theta}\\bar{g}_0)\\|_{H^{s-8}}$ up to $t = \\varepsilon^{-2/(N+1)}$; the theorem predicts this stays $O((1+t)^N \\varepsilon^2)$, so observing $O(\\varepsilon)$ growth, loss of $H^{s-8}$ regularity, or breakdown before $T_0 \\varepsilon^{-\\alpha}$ would refute it. As an analytic check, one can verify $\\Phi_{n,p} = 0$ in (2.5) for $n \\leq K$ using the explicit recursive $U_{n,p}$; any nonzero residual at order $\\leq K$ contradicts the construction.","tokens_in":27347,"feed_emoji":"⚛️","tokens_out":6018,"duration_ms":63485,"temperature":0.7,"texified_at":"2026-08-05T20:27:51.369525+00:00","pith_summary":"The paper analyzes the small-$\\varepsilon$ limit of the two-dimensional cubic Klein-Gordon equation, where $\\varepsilon$ is inversely proportional to the speed of light. It claims that in the defocusing case, the solution stays within $O(\\varepsilon^2)$ of the cubic nonlinear Schrödinger solution on a long time interval growing like $\\varepsilon^{-2/(N+1)}$, with an error prefactor that grows only algebraically in time. With less regular data, the rate slows to $O(\\varepsilon)$ on a shorter time interval. This establishes that the Schrödinger approximation is not merely formal but quantitatively accurate over the physically relevant non-relativistic timescale, in Sobolev norms rather than only $L^2$.","texify_model":"deepseek-v4-flash","texify_usage":{"total_tokens":5618,"prompt_tokens":947,"completion_tokens":4671,"prompt_tokens_details":{"cached_tokens":0},"prompt_cache_hit_tokens":0,"prompt_cache_miss_tokens":947,"completion_tokens_details":{"reasoning_tokens":3684}},"feed_headline":"2D Klein-Gordon matches cubic NLS to O(ε²) up to long times","feed_subtitle":"Error bound (1+t)^N ε² holds on the long-time interval ε^{-2/(N+1)}; lower regularity still gives O(ε).","key_machinery":"The engine is a WKB expansion with fast phase $\\theta = t/\\varepsilon^2$: the solution is written as a trigonometric polynomial $\\sum_n \\varepsilon^n \\sum_p e^{ip\\theta}U_{n,p}$, inserted into an equivalent symmetric hyperbolic system for $U = (\\varepsilon\\nabla u, \\varepsilon^2\\partial_t u, u)$. The lowest-order amplitudes satisfy the cubic NLS, and the higher profiles satisfy the linearized Schrödinger equations (3.47)–(3.48). The crucial mechanism is that the algebraic structure of these profiles cancels all resonant terms, leaving an error $\\varepsilon^{K+1}$ that is controlled by Duhamel stability together with the $t^{-1} L^\\infty$ decay of the defocusing cubic NLS.","core_discovery":"The central claim is that, for $\\lambda>0$ and initial data $(\\varphi,\\psi)$ in $H^s$ with $s>9$ plus weighted regularity, the solution $u$ of $\\varepsilon^2\\partial_{tt}u - \\Delta u + \\varepsilon^{-2}u + \\lambda u^3 = 0$ satisfies $\\|u - (e^{i\\theta}g_0 + e^{-i\\theta}\\bar{g}_0)\\|_{H^{s-8}} \\leq C_{\\varphi,\\psi}(1+t)^{N_{\\varphi,\\psi}} \\varepsilon^2$ for $t \\leq T_0 \\varepsilon^{-2/(N+1)}$, where $\\theta = t/\\varepsilon^2$ and $g_0$ solves the cubic NLS $2i\\partial_t g_0 - \\Delta g_0 + 3\\lambda|g_0|^2g_0 = 0$ with $g_0(0) = (\\varphi - i\\psi)/2$. For initial data with limited regularity, $s>5$, the error is $O((1+t) + (1+t)^{\\tilde N})\\varepsilon$ over a shorter time interval. The proof constructs WKB approximate solutions to arbitrary even order $K$, yielding error $\\varepsilon^{K+1}$ in $H^{s-2K-4}$ over the same long-time scale, and for the focusing case the approximation holds until the NLS solution's maximal ex","pith_inferences":["If the quoted structural WKB identities hold, the same stability framework should extend to related two-dimensional dispersive limits, such as Klein-Gordon–Zakharov systems, with the decay rates of the limit system entering the prefactor.","The algebraic growth (1+t)^N likely reflects the borderline t^{-1} decay of the two-dimensional cubic NLS; longer time scales would require additional structure such as scattering or weighted estimates.","The arbitrary-order result suggests a quantitative numerical recipe: schemes that include the ε²u2 correction should see the error drop by two powers of ε, a directly testable prediction.","Because the NLS decay prefactor depends on the profile of the initial datum, two initial data with identical norms can produce very different error constants, so numerical comparisons should report the full profile-dependence."],"forward_implications":["For defocusing smooth data, keeping only the leading profile already yields a rigorous O(ε²) approximation in H^{s-8} over times of order ε^{-2/(N+1)}, so reduced Schrödinger models inherit a concrete error bound.","Increasing regularity allows arbitrary-order approximations: including profiles up to order K lowers the error to ε^{K+1}, matching numerical observations of higher-order corrections.","The convergence holds in Sobolev spaces, not just L², so derivatives of the solution are also controlled by the Schrödinger approximation.","For focusing nonlinearity, the NLS description is validated up to the time when the Schrödinger solution itself may blow up, setting a fundamental limit on the non-relativistic approximation.","With less regular initial data the rate drops to O(ε), giving an explicit regularity-versus-rate tradeoff."],"supporting_citations":[{"why":"Supplies the algebraic structure of the WKB profiles (Propositions 3.4–3.6) and the three-dimensional convergence-rate method adapted here to two dimensions.","marker":"[7]"},{"why":"Provides local and global well-posedness of the cubic NLS, giving the existence of the leading profile g0.","marker":"[8]"},{"why":"Gives the L∞ decay ∥g0(t)∥∞ ≤ C(1+t)^{-1} for the defocusing cubic NLS, which is essential to the long-time stability estimate.","marker":"[12]"},{"why":"Earlier two-dimensional error estimate in L² norm; this paper extends the comparison to high-order Sobolev norms with longer time scales.","marker":"[22]"},{"why":"Introduces the geometric-optics stability framework and proves an O(ε) long-time error for quadratic nonlinearities, serving as the template for the stability analysis.","marker":"[13]"},{"why":"Numerical observation of O(ε²) and higher-order NLS approximation that motivates and corroborates the theorems.","marker":"[21]"}],"fun_headline_variants":["KG to NLS limit: ε² error, long-time validity","2D cubic Klein-Gordon matches NLS to O(ε²)","Non-relativistic limit: KG to NLS at ε² rate","For long times, KG error vs NLS is O(ε²)","Klein-Gordon to NLS: error bound (1+t)^N ε²"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The paper assumes, without proving here, that the WKB profiles have the exact algebraic form stated in Propositions 3.4, 3.5, and Corollary 3.6 and satisfy the profile equations (3.47)–(3.48); these structural identities are quoted from a prior paper. If they failed, the approximate solution would not exist and the error estimates would collapse.","fun_headline_variants_meta":{"raw":{"variants":["KG to NLS limit: ε² error, long-time validity","2D cubic Klein-Gordon matches NLS to O(ε²)","Non-relativistic limit: KG to NLS at ε² rate","For long times, KG error vs NLS is O(ε²)","Klein-Gordon to NLS: error bound (1+t)^N ε²"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000244,"raw_usage":{"total_tokens":1448,"prompt_tokens":903,"completion_tokens":545,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":647,"completion_tokens_details":{"reasoning_tokens":459}},"tokens_in":647,"tokens_out":545,"duration_ms":6587,"temperature":1.0,"reasoning_tokens":459,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T20:53:20.648117+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $u$ numerically for the defocusing 2D cubic Klein-Gordon equation with $\\varepsilon = 10^{-2}$ and smooth localized data, solving the cubic NLS for $g_0$ simultaneously, and measure $\\|u - (e^{i\\theta}g_0 + e^{-i\\theta}\\bar{g}_0)\\|_{H^{s-8}}$ up to $t = \\varepsilon^{-2/(N+1)}$; the theorem predicts this stays $O((1+t)^N \\varepsilon^2)$, so observing $O(\\varepsilon)$ growth, loss of $H^{s-8}$ regularity, or breakdown before $T_0 \\varepsilon^{-\\alpha}$ would refute it. As an analytic check, one can verify $\\Phi_{n,p} = 0$ in (2.5) for $n \\leq K$ using the explicit recursive $U_{n,p}$; any nonzero residual at order $\\leq K$ contradicts the construction.","supporting_citations":[],"review_version":1}