Uniform error bounds of an exponential wave integrator for the Klein–Gordon–Schrödinger system in the nonrelativistic and massless limit regime.

Saved in:
Bibliographic Details
Title: Uniform error bounds of an exponential wave integrator for the Klein–Gordon–Schrödinger system in the nonrelativistic and massless limit regime.
Authors: Li, Jiyong1 (AUTHOR) ljyong406@163.com, Yang, Minghui (AUTHOR)
Source: Mathematics & Computers in Simulation. Jul2025, Vol. 233, p237-258. 22p.
Subjects: Systems integrators, Oscillations, Wavelengths
Abstract: We propose an exponential wave integrator Fourier pseudo-spectral (EWI-FP) method and establish the uniform error bounds for the Klein–Gordon–Schrödinger system (KGSS) with ɛ ∈ (0 , 1 ]. In the nonrelativistic and massless limit regime (0 < ɛ ≪ 1), the solution of KGSS propagates waves with wavelength O (ɛ) in time and amplitude at O (ɛ α † ) where α † = min { α , β + 1 , 2 } with two parameters α and β. The parameters satisfy α ≥ 0 and β ≥ − 1. In this regime, due to the oscillation in time, it is very difficult to develop efficient schemes and make the corresponding error analysis for KGSS. In this paper, firstly, in order to overcome the difficulty of controlling the nonlinear terms, we transform the KGSS into a system with higher derivative. Then we construct an EWI-FP method and provide the error estimates with two bounds at O (h σ + 2 + min { τ / ɛ 1 − α ∗ , τ 2 / ɛ 2 − α † }) and O (h σ + 2 + τ 2 + ɛ α † ) , respectively, where α ∗ = min { 1 , α , 1 + β } , σ has to do with the smoothness of the solution in space, h is mesh size and τ is time step. From the two error bounds, we obtain the error estimates O (h σ + 2 + τ α † ) for α ≥ 0 and β ≥ − 1. Hence, we get uniform second-order error bounds at O (h σ + 2 + τ 2) in time when α ≥ 2 and β ≥ 1 , and uniformly accurate first-order error estimates for any α ≥ 1 and β ≥ 0. We also get uniformly accurate spatial spectral accuracy for any α ≥ 0 and β ≥ − 1. Our numerical results support our conclusions. [ABSTRACT FROM AUTHOR]
Copyright of Mathematics & Computers in Simulation is the property of Elsevier B.V. and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
Database: Engineering Source
Description
Abstract:We propose an exponential wave integrator Fourier pseudo-spectral (EWI-FP) method and establish the uniform error bounds for the Klein–Gordon–Schrödinger system (KGSS) with ɛ ∈ (0 , 1 ]. In the nonrelativistic and massless limit regime (0 < ɛ ≪ 1), the solution of KGSS propagates waves with wavelength O (ɛ) in time and amplitude at O (ɛ α † ) where α † = min { α , β + 1 , 2 } with two parameters α and β. The parameters satisfy α ≥ 0 and β ≥ − 1. In this regime, due to the oscillation in time, it is very difficult to develop efficient schemes and make the corresponding error analysis for KGSS. In this paper, firstly, in order to overcome the difficulty of controlling the nonlinear terms, we transform the KGSS into a system with higher derivative. Then we construct an EWI-FP method and provide the error estimates with two bounds at O (h σ + 2 + min { τ / ɛ 1 − α ∗ , τ 2 / ɛ 2 − α † }) and O (h σ + 2 + τ 2 + ɛ α † ) , respectively, where α ∗ = min { 1 , α , 1 + β } , σ has to do with the smoothness of the solution in space, h is mesh size and τ is time step. From the two error bounds, we obtain the error estimates O (h σ + 2 + τ α † ) for α ≥ 0 and β ≥ − 1. Hence, we get uniform second-order error bounds at O (h σ + 2 + τ 2) in time when α ≥ 2 and β ≥ 1 , and uniformly accurate first-order error estimates for any α ≥ 1 and β ≥ 0. We also get uniformly accurate spatial spectral accuracy for any α ≥ 0 and β ≥ − 1. Our numerical results support our conclusions. [ABSTRACT FROM AUTHOR]
ISSN:03784754
DOI:10.1016/j.matcom.2025.01.027