Abstract:We introduce and analyze a symmetric low-regularity scheme for the nonlinear Schrödinger (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.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to perform reliable numerical approximation of the nonlinear Schrödinger equation (NLS) under low - regularity conditions. Specifically, the author introduces a new symmetric low - regularity integrator for approximating the nonlinear Schrödinger equation in a more general setting than the classical Fourier method. This new method allows for symmetric approximation under lower regularity assumptions and works on more general domains and with fractional convergence rates. Moreover, this method exhibits better structure - preserving properties and improved error constants under low - regularity conditions.
### Background and Problem Description of the Paper
**Nonlinear Schrödinger Equation (NLS)**:
\[ i\partial_t u(t, x) = -\Delta u(t, x) + |u(t, x)|^2 u(t, x), \quad (t, x) \in \mathbb{R} \times \Omega \]
where \(\Omega \subset \mathbb{R}^d\) or \(\Omega = \mathbb{T}^d\), and the initial condition is \(u|_{t = 0}=u_0\). When \(\partial\Omega\neq\emptyset\), assume that \(\Omega\) is a smooth bounded domain and is endowed with homogeneous boundary conditions.
**Research Objectives**:
1. **Low - regularity approximation**: Provide reliable numerical approximation when the initial data \(u_0\) is non - smooth (i.e., belongs to a low - order Sobolev space).
2. **Symmetry**: Design a symmetric numerical method to preserve time - reversibility and conservation properties.
3. **Generality**: The method is applicable not only to periodic boundary conditions but also to more general domains and boundary conditions.
### Main Contributions
1. **Symmetric low - regularity integrator**: A new symmetric low - regularity integrator is introduced, which exhibits good convergence and structure - preserving properties under low - regularity conditions.
2. **Convergence analysis**: Prove the fractional convergence rate of this method in the \(L^2\) norm, from first - order to second - order.
3. **Numerical experiments**: Verify the superior performance of the new method under low - regularity conditions through numerical experiments, especially in terms of structure - preserving properties and error constants.
### Mathematical Details
**Symmetric low - regularity integrator**:
\[ u^{n + 1}=\phi_\tau(u^n)=e^{i\tau\Delta}u^n - i\tau\left(\frac{1}{2}e^{i\tau\Delta}(u^n)^2\varphi_1(-i\tau\Delta)u^n+\frac{1}{2}(u^{n + 1})^2\varphi_1(i\tau\Delta)u^{n + 1}\right) \]
where \(\varphi_1(z)=\frac{e^z - 1}{z}\) and \(\tau\) is the time step size.
**Convergence Results**:
- On the periodic domain \(\mathbb{T}^d\), for \(u_0\in H^\alpha(\mathbb{T}^d)\) and \(\alpha\in\left[1+\frac{d}{4}, 3\right]\), the numerical solution \(u^n\) satisfies:
\[ \|u(n\tau)-u^n\|_{L^2}\leq C_T\left(\sup_{[0, T]}\|u(t)\|_{H^\alpha}\right)\tau^{1 + \gamma} \]
where \(\gamma\) satisfies specific conditions.
- On the smooth bounded domain \(\Omega\subset\mathbb{R}^d\), for \(u_0\in(H^{1 + d/4}\cap H_0^1)(\Omega)\), the numerical solution \(u^n\) satisfies first - order convergence; for \(u_0\in(H^\alpha\cap H_0^1)\)