Improved Uniform Error Bounds on Time-Splitting Methods for Long-Time Dynamics of the Nonlinear Klein--Gordon Equation with Weak Nonlinearity

We establish improved uniform error bounds on a second-order Strang time-splitting method which is equivalent to an exponential wave integrator for the long-time dynamics of the nonlinear Klein--Gordon equation (NKGE) with weak cubic nonlinearity, whose strength is characterized by ε² with 0 < ε ≤ 1 a dimensionless parameter. Actually, when 0 < ε 1, the NKGE with O(ε²) nonlinearity and $O(1)$ initial data is equivalent to that with $O(1)$ nonlinearity and small initial data, the amplitude of which is at O(ε). We begin with a semidiscretization of the NKGE by the second-order time-splitting method and derive a full-discretization by the Fourier spectral method in space. Employing the regularity compensation oscillation technique which controls the high frequency modes by the regularity of the exact solution and analyzes the low frequency modes by phase cancellation and energy method, we carry out the improved uniform error bounds at O(ε²τ²) and O(hᵐ+ε²τ²) for the second-order semidiscretization and full-discretization up to the long time T_ε = T/ε² with T fixed, respectively. Extensions to higher-order time-splitting methods and the case of an oscillatory complex NKGE are also discussed. Finally, numerical results are provided to confirm the improved error bounds and to demonstrate that they are sharp.