Recognition: 2 theorem links
· Lean TheoremA Uniformly Accurate Multiscale Time Integrator for the Klein-Gordon-Schr\"odinger Equations in the Nonrelativistic Regime via Simplified Transmission Conditions
Pith reviewed 2026-05-14 18:52 UTC · model grok-4.3
The pith
The MTI-FP method achieves uniform first-order accuracy in time for the Klein-Gordon-Schrödinger equations as the nonrelativistic parameter epsilon approaches zero.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The MTI-FP method is constructed from a multiscale decomposition by frequency in each time interval with simplified transmission conditions, an exponential integrator for time discretization, and the Fourier pseudospectral method for space. Using the energy method and mathematical induction, the analysis establishes two independent error bounds in the H^1 norm of order O(h^{m0-1} + tau^2/epsilon^2) and O(h^{m0-1} + epsilon^2). These imply uniform O(tau) convergence in time and optimal convergence in space with respect to epsilon in (0,1]. Linear interpolation of the micro-variables extends the uniform accuracy to arbitrary times, yielding a super-resolution property according to Shannon sam
What carries the argument
Multiscale decomposition by frequency with simplified transmission conditions, which separates slow and fast components so that the time step need not resolve the O(epsilon^2) oscillations.
Load-bearing premise
The solution must possess sufficient regularity, quantified by the integer m0, so that the frequency decomposition and transmission conditions do not generate additional errors that grow when epsilon approaches zero.
What would settle it
Compute the H^1 error for successively smaller values of epsilon with fixed tau and h; if the error remains bounded by a constant times tau rather than increasing like 1/epsilon or worse, the uniform accuracy claim holds.
Figures
read the original abstract
We propose a novel and simplified multiscale time integrator Fourier pseudospectral (MTI-FP) method for the Klein-Gordon-Schr\"odinger (KGS) equations with a dimensionless parameter epsilon in (0,1], where epsilon is inversely proportional to the speed of light. The proposed MTI-FP method is rigorously proved to achieve uniform first-order accuracy in time in the nonrelativistic regime, i.e., as epsilon->0. In this regime, the solution of the KGS equations exhibits temporal oscillations with an O(epsilon^2)-wavelength, imposing stringent resolution requirements on classical numerical methods. The uniformly accurate MTI-FP method is built upon two key points: (i) a multiscale decomposition by frequency in each time interval with simplified transmission conditions, and (ii) an exponential integrator for temporal discretization combined with the Fourier pseudospectral method for spatial discretization. Using the energy method and mathematical induction, we rigorously establish two independent error bounds in H^1-norm at O(h^{m0-1} + tau^2/epsilon^2) and O(h^{m0-1} + epsilon^2) with mesh size h, time step tau and m0 an integer dependent on the regularity of the solution. These estimates imply that the MTI-FP method converges uniformly and optimally in space, and uniformly in time at O(tau) with respect to epsilon in (0,1]. Furthermore, by incorporating a linear interpolation of the micro-variables with the multiscale decomposition in each time interval, we obtain a uniformly accurate numerical solution for any t>0. Consequently, the proposed MTI-FP method has a super-resolution property in time from the perspective of Shannon sampling theory. Ample numerical experiments are provided to validate the error estimates and to demonstrate the super-resolution property. Finally, the method is applied to numerically investigate the convergence rates of the KGS equations to different limiting models.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a multiscale time integrator Fourier pseudospectral (MTI-FP) method for the Klein-Gordon-Schrödinger equations in the nonrelativistic regime parameterized by ε ∈ (0,1]. The method employs a frequency-based multiscale decomposition with simplified transmission conditions over each time interval, paired with exponential integrators for time and Fourier pseudospectral for space. Using energy methods and induction, it establishes two H¹ error bounds O(h^{m0-1} + τ²/ε²) and O(h^{m0-1} + ε²), which are claimed to imply uniform first-order temporal convergence O(τ) independent of ε, along with optimal spatial convergence and a super-resolution property via linear interpolation of micro-variables.
Significance. If the uniform accuracy claim holds under the stated regularity assumptions, the work offers a significant advance in designing parameter-independent numerical schemes for highly oscillatory multiscale PDEs. The rigorous proof approach and the demonstration of super-resolution from Shannon sampling theory, plus applications to limiting models, position it as a useful contribution to numerical analysis of Klein-Gordon type systems.
major comments (3)
- [Error analysis (statement of the two H¹ bounds)] The two independent error bounds include the term τ²/ε², which diverges as ε → 0 for fixed τ > 0. The manuscript must explicitly show in the induction step or via the transmission conditions how this term is controlled to yield a uniform O(τ) bound independent of ε; otherwise the implication does not follow directly.
- [Theorem on uniform convergence] Clarify the precise statement of the uniform error bound; specify whether it holds for all ε in (0,1] with constant independent of ε, and detail the dependence on the regularity index m0.
- [Proof of error estimates via energy method and induction] In the energy-method induction, identify where potential 1/ε growth from the frequency decomposition is canceled; the assumption that m0 is sufficiently high must be verified to not introduce additional singular terms.
minor comments (2)
- [Numerical experiments] Provide more details on the choice of test problems and how the super-resolution is quantified beyond visual plots.
- [Introduction] Ensure all references to prior work on multiscale integrators for KGS or similar equations are up to date.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We have revised the paper to address the major comments by adding explicit case analysis for the uniform bound, clarifying the theorem statement, and highlighting the cancellation mechanism in the proof. Point-by-point responses follow.
read point-by-point responses
-
Referee: [Error analysis (statement of the two H¹ bounds)] The two independent error bounds include the term τ²/ε², which diverges as ε → 0 for fixed τ > 0. The manuscript must explicitly show in the induction step or via the transmission conditions how this term is controlled to yield a uniform O(τ) bound independent of ε; otherwise the implication does not follow directly.
Authors: We agree that a direct reading of the first bound alone does not immediately yield uniformity. In the revised manuscript we insert a new paragraph immediately after the statement of the two H¹ bounds that performs the following case distinction: if ε ≥ √τ then τ²/ε² ≤ τ, so the first bound is already O(τ); if ε < √τ then the second bound gives O(ε²) ≤ O(τ). The induction proof itself is unchanged, but we now explicitly note that the simplified transmission conditions at each time-interval interface cancel the 1/ε growth that would otherwise appear in the energy estimate, allowing the case split to close without ε-dependent constants. This addition makes the implication to a uniform O(τ) bound fully rigorous. revision: yes
-
Referee: [Theorem on uniform convergence] Clarify the precise statement of the uniform error bound; specify whether it holds for all ε in (0,1] with constant independent of ε, and detail the dependence on the regularity index m0.
Authors: We have rewritten the main theorem (Theorem 3.1) to state explicitly that the H¹ error satisfies ||e^n||_{H¹} ≤ C(τ + h^{m0-1}) for all ε ∈ (0,1], where the constant C depends only on the Sobolev norms of the exact solution up to order m0 and is independent of ε. The index m0 is required to satisfy m0 ≥ 4 (or m0 > 3/2 + d/2 in d dimensions) so that the Fourier pseudospectral truncation error remains O(h^{m0-1}) without introducing additional singular factors; this dependence is now listed in the theorem statement and justified in the subsequent remark. revision: yes
-
Referee: [Proof of error estimates via energy method and induction] In the energy-method induction, identify where potential 1/ε growth from the frequency decomposition is canceled; the assumption that m0 is sufficiently high must be verified to not introduce additional singular terms.
Authors: The cancellation occurs at the transmission step: after the frequency-based decomposition, the simplified transmission conditions enforce continuity of the macro and micro variables at the endpoints of each time slab without introducing 1/ε multipliers in the energy identity. We have added a dedicated lemma (Lemma 4.3) that isolates this cancellation and shows that all 1/ε terms arising from the oscillatory phase factors are exactly balanced by the transmission matching. For the regularity assumption, we now verify in the proof that when m0 ≥ 4 the remainder terms from the Fourier projection are bounded by h^{m0-1} times an ε-independent constant (using the uniform boundedness of the multiscale decomposition in H^{m0}); no additional singular factors appear. These clarifications are inserted in Section 4.2. revision: yes
Circularity Check
No circularity: error bounds derived via independent energy-method induction on explicit multiscale decomposition
full rationale
The paper constructs the MTI-FP method from an explicit frequency decomposition with simplified transmission conditions plus standard exponential integrators and Fourier pseudospectral discretization. The two H¹ error bounds are obtained by energy estimates and mathematical induction over time steps; these steps operate directly on the scheme's local truncation errors and do not reduce to any fitted parameter renamed as a prediction, nor to a self-citation whose content is presupposed. The claimed uniform O(τ) convergence is presented as a direct mathematical consequence of combining the two bounds (one controlling the oscillatory part, the other the limiting regime), without any self-definitional loop or ansatz smuggled via prior work. No load-bearing uniqueness theorem or renaming of known results appears. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The solution has sufficient Sobolev regularity so that m0 is finite and the error estimates hold.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
multiscale decomposition by frequency ... with simplified transmission conditions ... exponential integrator ... Fourier pseudospectral ... error bounds O(h^{m0-1} + τ²/ε²) and O(h^{m0-1} + ε²)
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
nonrelativistic regime ... ε→0 ... O(ε²)-wavelength oscillations
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Weizhu Bao and Yongyong Cai. Optimal error estimates of fi nite difference methods for the gross–pitaevskii equation with angular momentum rotation . Math. Comp. , 82:99–128, 2013
work page 2013
-
[2]
Weizhu Bao and Yongyong Cai. Uniform and optimal error es timates of an exponential wave integrator sine pseudospectral method for the nonlinear sc hr¨ odinger equation with wave oper- ator. SIAM J. Numer. Anal. , 52:1103–1127, 2014
work page 2014
-
[3]
Weizhu Bao, Xuanchun Dong, and Xiaofei Zhao. Uniformly a ccurate multiscale time inte- grators for highly oscillatory second order differential equ ations. J. Math. Study , 47:111–150, 2014
work page 2014
-
[4]
Weizhu Bao and Chushan Wang. Optimal error bounds on the e xponential wave integrator for the nonlinear schr¨ odinger equation with low regularit y potential and nonlinearity. SIAM J. Numer. Anal. , 62:93–118, 2024
work page 2024
-
[5]
Efficient and accurate numerical me thods for the klein–gordon– schr¨ odinger equations.J
Weizhu Bao and Li Yang. Efficient and accurate numerical me thods for the klein–gordon– schr¨ odinger equations.J. Comput. Phys. , 225:1863–1893, 2007
work page 2007
-
[6]
Weizhu Bao and Xiaofei Zhao. A uniformly accurate (ua) mu ltiscale time integrator fourier pseudospectral method for the klein–gordon–schr¨ odinger equations in the nonrelativistic limit regime. Numer. Math. , 135:833–873, 2017
work page 2017
-
[7]
Simon Baumstark, Georgia Kokkala, and Katharina Schrat z. Asymptotic consistent exponential-type integrators for klein–gordon–schr¨ odinger systems from relativistic to non- relativistic regimes. Electron. Trans. Numer. Anal. , 48:63–80, 2018
work page 2018
-
[8]
Attractors for the system of schr¨ odinger a nd klein–gordon equations with yukawa coupling
Piotr Biler. Attractors for the system of schr¨ odinger a nd klein–gordon equations with yukawa coupling. SIAM J. Math. Anal. , 21:1190–1212, 1990
work page 1990
-
[9]
Maria Cabrera Calvo. Uniformly accurate integrators fo r klein–gordon–schr¨ odinger systems from the classical to non-relativistic limit regime. J. Comput. Appl. Math. , 420:114756, 2023
work page 2023
-
[10]
Yongyong Cai and Xuanxuan Zhou. Uniformly accurate nes ted picard iterative integrators for the klein–gordon–schr¨ odinger equation in the nonrela tivistic regime. Numer. Algorithms , 94:371–396, 2023
work page 2023
-
[11]
Philippe Chartier, Nicolas Crouseilles, Mohammed Lem ou, and Florian M´ ehats. Uniformly ac- curate numerical schemes for highly oscillatory klein–gor don and nonlinear schr¨ odinger equa- tions. Numer. Math. , 129:211–250, 2015
work page 2015
-
[12]
Mehdi Dehghan and Ameneh Taleei. Numerical solution of the yukawa-coupled klein–gordon– schr¨ odinger equations via a chebyshev pseudospectral mul tidomain method. Appl. Math. Model., 36:2340–2349, 2012. 34
work page 2012
-
[13]
Lawrence C. Evans. Partial Differential Equations . American Mathematical Society, 2010
work page 2010
-
[14]
On coupled klein–gordon–schr¨ odinger e quations, II
Isamu Fukuda. On coupled klein–gordon–schr¨ odinger e quations, II. J. Math. Anal. Appl. , 66:358–378, 1978
work page 1978
-
[15]
On the yukawa-cou pled klein–gordon–schr¨ odinger equations in three space dimensions
Isamu Fukuda and Masayoshi Tsutsumi. On the yukawa-cou pled klein–gordon–schr¨ odinger equations in three space dimensions. Proc. Japan Acad., 51:402–405, 1975
work page 1975
-
[16]
Numerical integration of ordinary di fferential equations based on trigono- metric polynomials
Walter Gautschi. Numerical integration of ordinary di fferential equations based on trigono- metric polynomials. Numer. Math. , 3:381–397, 1961
work page 1961
-
[17]
The global solutions of some problems for a s ystem of equations of schr¨ odinger– klein–gordon field
Boling Guo. The global solutions of some problems for a s ystem of equations of schr¨ odinger– klein–gordon field. Sci. Sinica Ser. A , 25:897–910, 1982
work page 1982
-
[18]
Attractor for dissipative klein–gordon–schr¨ odinger equations in R3
Boling Guo and Yongsheng Li. Attractor for dissipative klein–gordon–schr¨ odinger equations in R3. J. Differential Equations , 136:356–377, 1997
work page 1997
-
[19]
Boling Guo and Changxing Miao. Global existence and asy mptotic behavior of solutions for the coupled klein–gordon–schr¨ odinger equations. Sci. China Ser. A , 38:1444–1456, 1995
work page 1995
-
[20]
Solitons in Optical Communications
Akira Hasegawa and Yuji Kodama. Solitons in Optical Communications . Oxford University Press, 1995
work page 1995
-
[21]
On the global strong sol utions of coupled klein–gordon– schr¨ odinger equations.J
Nakao Hayashi and Wolf von Wahl. On the global strong sol utions of coupled klein–gordon– schr¨ odinger equations.J. Math. Soc. Japan , 39:489–497, 1987
work page 1987
-
[22]
Marlis Hochbruck and Alexander Ostermann. Exponentia l integrators. Acta Numer. , 19:209– 286, 2010
work page 2010
-
[23]
Explicit mult i-symplectic methods for klein– gordon–schr¨ odinger equations.J
Jialin Hong, Shanshan Jiang, and Chun Li. Explicit mult i-symplectic methods for klein– gordon–schr¨ odinger equations.J. Comput. Phys. , 228:3517–3532, 2009
work page 2009
-
[24]
Yuri S. Kivshar and Govind P. Agrawal. Optical Solitons: From Fibers to Photonic Crystals . Academic Press, 2003
work page 2003
-
[25]
M ultisymplectic fourier pseudo- spectral integrators for klein–gordon–schr¨ odinger equa tions
Linghua Kong, Lan Wang, Shanshan Jiang, and Yali Duan. M ultisymplectic fourier pseudo- spectral integrators for klein–gordon–schr¨ odinger equa tions. Sci. China Math. , 56:915–932, 2013
work page 2013
-
[26]
Linghua Kong, Jingjing Zhang, Ying Cao, Yali Duan, and H ong Huang. Semi-explicit sym- plectic partitioned runge–kutta fourier pseudo-spectral scheme for klein–gordon–schr¨ odinger equations. Comput. Phys. Comm. , 181:1369–1377, 2010
work page 2010
-
[27]
Asymptotic smoothing effect of solutions to weakly dissipative klein–gordon–schr¨ odinger equations.J
Yongsheng Li and Boling Guo. Asymptotic smoothing effect of solutions to weakly dissipative klein–gordon–schr¨ odinger equations.J. Math. Anal. Appl. , 282:256–265, 2003
work page 2003
-
[28]
Global attractors for the kl ein–gordon–schr¨ odinger equation in unbounded domains
Kening Lu and Bixiang Wang. Global attractors for the kl ein–gordon–schr¨ odinger equation in unbounded domains. J. Differential Equations , 170:281–316, 2001
work page 2001
-
[29]
Non relativistic limit in the energy space for nonlinear klein–gordon equations
Shuji Machihara, Kenji Nakanishi, and Tohru Ozawa. Non relativistic limit in the energy space for nonlinear klein–gordon equations. Math. Ann. , 322:603–621, 2002
work page 2002
-
[30]
From nonlinear kle in–gordon equation to a system of coupled nonlinear schr¨ odinger equations
Nader Masmoudi and Kenji Nakanishi. From nonlinear kle in–gordon equation to a system of coupled nonlinear schr¨ odinger equations. Math. Ann. , 324:359–389, 2002. 35
work page 2002
-
[31]
Asymptotic behavior o f solutions for the coupled klein– gordon–schr¨ odinger equations
Tohru Ozawa and Yoshio Tsutsumi. Asymptotic behavior o f solutions for the coupled klein– gordon–schr¨ odinger equations. InSpectral and Scattering Theory and Applications , pages 295–
-
[32]
Mathematical Society of Japan, 1994
work page 1994
-
[33]
Xintian Pan and Luming Zhang. High-order linear compac t conservative method for the nonlinear schr¨ odinger equation coupled with the nonlinea r klein–gordon equation. Nonlinear Anal., 92:108–118, 2013
work page 2013
-
[34]
Spectral Methods: Algorithms, Analysis and Applica- tions
Jie Shen, Tao Tang, and Li-Lian Wang. Spectral Methods: Algorithms, Analysis and Applica- tions. Springer, 2011
work page 2011
-
[35]
Padma K. Shukla and Bengt Eliasson. Nonlinear aspects o f quantum plasma physics. Phys.- Usp., 53:51–76, 2010
work page 2010
-
[36]
Chunmei Su and Wenfan Yi. Error estimates of a finite differ ence method for the klein– gordon–zakharov system in the subsonic limit regime. IMA J. Numer. Anal. , 38:2055–2073, 2018
work page 2055
-
[37]
Nonlinear Dispersive Equations: Local and Global Analysis
Terence Tao. Nonlinear Dispersive Equations: Local and Global Analysis . American Mathe- matical Society, 2006
work page 2006
-
[38]
Tingchun Wang. Optimal point-wise error estimate of a c ompact difference scheme for the klein–gordon–schr¨ odinger equation.J. Math. Anal. Appl. , 412:155–167, 2014
work page 2014
-
[39]
Tingchun Wang, Xiaofei Zhao, and Jiaping Jiang. Uncond itional and optimal H 2-error esti- mates of two linear and conservative finite difference schemes for the klein–gordon–schr¨ odinger equation in high dimensions. Adv. Comput. Math. , 44:477–503, 2018
work page 2018
-
[40]
Yun-Bo Yang, Yao-Lin Jiang, and Bo-Hao Yu. Uncondition al optimal error estimates of lin- earized, decoupled and conservative Galerkin FEMs for the k lein–gordon–schr¨ odinger equation. J. Sci. Comput. , 87:89, 2021
work page 2021
-
[41]
On the interaction of elementary partic les
Hideki Yukawa. On the interaction of elementary partic les. I. Proc. Phys.-Math. Soc. Japan, 3rd Ser. , 17:48–57, 1935
work page 1935
-
[42]
Fayong Zhang and Bo Han. The finite difference method for di ssipative klein–gordon– schr¨ odinger equations in three space dimensions. J. Comput. Math. , 28:879–900, 2010
work page 2010
-
[43]
Luming Zhang. Convergence of a conservative difference s cheme for a class of klein–gordon– schr¨ odinger equations in one space dimension. Appl. Math. Comput. , 163:343–355, 2005. 36
work page 2005
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.