Optimal error bounds on the exponential wave integrator for the nonlinear Schrödinger equation with low regularity potential and nonlinearity

by   Weizhu Bao, et al.

We establish optimal error bounds for the exponential wave integrator (EWI) applied to the nonlinear Schrödinger equation (NLSE) with L^∞-potential and/or locally Lipschitz nonlinearity under the assumption of H^2-solution of the NLSE. For the semi-discretization in time by the first-order Gautschi-type EWI, we prove an optimal L^2-error bound at O(τ) with τ>0 being the time step size, together with a uniform H^2-bound of the numerical solution. For the full-discretization scheme obtained by using the Fourier spectral method in space, we prove an optimal L^2-error bound at O(τ + h^2) without any coupling condition between τ and h, where h>0 is the mesh size. In addition, for W^1, 4-potential and a little stronger regularity of the nonlinearity, under the assumption of H^3-solution, we obtain an optimal H^1-error bound. Furthermore, when the potential is of low regularity but the nonlinearity is sufficiently smooth, we propose an extended Fourier pseudospectral method which has the same error bound as the Fourier spectral method while its computational cost is similar to the standard Fourier pseudospectral method. Our new error bounds greatly improve the existing results for the NLSE with low regularity potential and/or nonlinearity. Extensive numerical results are reported to confirm our error estimates and to demonstrate that they are sharp.


page 1

page 2

page 3

page 4


Uniform error bound of an exponential wave integrator for the long-time dynamics of the nonlinear Schrödinger equation with wave operator

We establish the uniform error bound of an exponential wave integrator F...

Uniform L^∞-bounds for energy-conserving higher-order time integrators for the Gross-Pitaevskii equation with rotation

In this paper, we consider an energy-conserving continuous Galerkin disc...

Error Analysis and Correction for Weighted A*'s Suboptimality (Extended Version)

Weighted A* (wA*) is a widely used algorithm for rapidly, but suboptimal...

Error estimates of the time-splitting methods for the nonlinear Schrödinger equation with semi-smooth nonlinearity

We establish error bounds of the Lie-Trotter time-splitting sine pseudos...

Helmholtz FEM solutions are locally quasi-optimal modulo low frequencies

For h-FEM discretisations of the Helmholtz equation with wavenumber k, w...

On optimal zero-padding of kernel truncation method

The kernel truncation method (KTM) is a commonly-used algorithm to compu...

Sharp Analysis of Random Fourier Features in Classification

We study the theoretical properties of random Fourier features classific...

Please sign up or login with your details

Forgot password? Click here to reset