On Fourier time-splitting methods for nonlinear Schrödinger equations in the semi-classical limit. II. Analytic regularity (Q527828)

The paper is concerned with the analysis of the Lie-Trotter time-splitting scheme for the nonlinear Schrödinger equation \[ i \varepsilon \partial_t u^{\varepsilon}+\frac{\varepsilon^2}{2} \Delta u^{\varepsilon} = \lambda \left| u^{\varepsilon} \right|^{2 \sigma} u^{\varepsilon}, \quad \left(x,t\right) \in \left[0,T\right] \times \mathbb{R}^d, \] \(\lambda \in \mathbb{R}\), \(\sigma \in \mathbb{N}\), in the semi-classical limit \(\varepsilon \rightarrow 0\). It is assumed that the initial function has a WKB structure \[ u^{\varepsilon} \left( 0, x\right) = a_0 \left(x \right) e^{i \phi_0 \left(x\right) / \varepsilon}. \] The previous results concerning convergence and error estimates for the Lie-Trotter splitting are derived for the Schrödinger equation with nonlinearities of other types or under the assumption that the Sobolev norm of the exact solution satisfies \(\left\| u^{\varepsilon} \right\|_{H^s} \leq C\) with \(C\) independent on \(\varepsilon\). However, as \(\varepsilon \rightarrow 0\), the solution \(u^{\varepsilon}\) is oscillatory and \(\left\| u^{\varepsilon} \right\|_{H^s} \approx \varepsilon^{-s}\). In the paper, it is shown that the exact solution and the numerical solution preserve the WKB structure and error estimates are derived under the assumption that the Fourier transforms \(\hat{\phi}_0\) and \(\hat{a}_0\) of \(\phi_0\) and \(a_0\), respectively, satisfy \[ \int\limits_{\mathbb{R}^d} e^{\left(1+ \left|\xi \right|^2 \right)^{\left(1+\delta\right)/2}} \left( \left| \hat{\phi}_0 \left( \xi \right)\right|^2 + \left| \hat{a}_0 \left( \xi \right) \right|^2 \right) d \xi < \infty \] for some \(\delta > 0\).