Abstract:We introduce a new non-resonant low-regularity integrator for the cubic nonlinear Schrödinger equation (NLSE) allowing for long-time error estimates which are optimal in the sense of the underlying PDE. The main idea thereby lies in treating the zeroth mode exactly within the discretization. For long-time error estimates, we rigorously establish the long-time error bounds of different low-regularity integrators for the nonlinear Schrödinger equation (NLSE) with small initial data characterized by a dimensionless parameter $\varepsilon \in (0, 1]$. We begin with the low-regularity integrator for the quadratic NLSE in which the integral is computed exactly and the improved uniform first-order convergence in $H^r$ is proven at $O(\varepsilon \tau)$ for solutions in $H^r$ with $r > 1/2$ up to the time $T_{\varepsilon } = T/\varepsilon $ with fixed $T > 0$. Then, the improved uniform long-time error bound is extended to a symmetric second-order low-regularity integrator in the long-time regime. For the cubic NLSE, we design new non-resonant first-order and symmetric second-order low-regularity integrators which treat the zeroth mode exactly and rigorously carry out the error analysis up to the time $T_{\varepsilon } = T/\varepsilon ^2$. With the help of the regularity compensation oscillation (RCO) technique, the improved uniform error bounds are established for the new non-resonant low-regularity schemes, which further reduce the long-time error by a factor of $\varepsilon^2$ compared with classical low-regularity integrators for the cubic NLSE. Numerical examples are presented to validate the error estimates and compare with the classical time-splitting methods in the long-time simulations.

Full- and low-rank exponential Euler integrators for the Lindblad equation

A Generic Method for Integrating Lindblad Master Equations

Quantum Simulation of Lindbladian Dynamics via Repeated Interactions

Taming the Bloch-Redfield equation: Recovering an accurate Lindblad equation for general open quantum systems

Low rank approximation for the numerical simulation of high dimensional Lindblad and Riccati equations

Stable and fast semi-implicit integration of the stochastic Landau-Lifshitz equation

A Quantum Kinetic Monte Carlo Method For Lindblad Equation

Long-time error bounds of low-regularity integrators for nonlinear Schrödinger equations

Wave Matrix Lindbladization II: General Lindbladians, Linear Combinations, and Polynomials

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

Optimal Convergence and Long-Time Conservation of Exponential Integration for Schrödinger Equations in a Normal or Highly Oscillatory Regime

Simulating Open Quantum Systems Using Hamiltonian Simulations

Explicit symmetric low-regularity integrators for the nonlinear Schrödinger equation

Solving the Lindblad equation with methods from computational fluid dynamics

Numerical integrators for confined Langevin dynamics

Optimal convergence and long-time conservation of exponential integration for Schr\"{o}dinger equations in a normal or highly oscillatory regime

Third order, uniform in low to high oscillatory coefficients, exponential integrators for Klein-Gordon equations

High-order linearly implicit structure-preserving exponential integrators for the nonlinear Schr\"odinger equation

Unconditionally stable exponential integrator schemes for the 2D Cahn-Hilliard equation

Efficient simulation of complex Ginzburg--Landau equations using high-order exponential-type methods