A symmetric low-regularity integrator for the nonlinear Schrödinger equation

We introduce and analyze a symmetric low-regularity scheme for the nonlinear Schr\"odinger (NLS) equation beyond classical Fourier-based techniques. We show fractional convergence of the scheme in $L^2$-norm, from first up to second order, both on the torus $\mathbb{T}^d$ and on a smooth bounded domain $\Omega \subset \mathbb{R}^d$, $d\le 3$, equipped with homogeneous Dirichlet boundary condition. The new scheme allows for a symmetric approximation to the NLS equation in a more general setting than classical splitting, exponential integrators, and low-regularity schemes (i.e. under lower regularity assumptions, on more general domains, and with fractional rates). We motivate and illustrate our findings through numerical experiments, where we witness better structure preserving properties and an improved error-constant in low-regularity regimes.