An explicit and symmetric exponential wave integrator for the nonlinear Schrödinger equation with low regularity potential and nonlinearity

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Bao, Weizhu, Wang, Chushan
Format: Preprint
Publié: 2023
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866910919636811776
author Bao, Weizhu
Wang, Chushan
author_facet Bao, Weizhu
Wang, Chushan
contents We propose and analyze a novel symmetric Gautschi-type exponential wave integrator (sEWI) for the nonlinear Schrödinger equation (NLSE) with low regularity potential and typical power-type nonlinearity of the form $ |ψ|^{2σ}ψ$ with $ ψ$ being the wave function and $ σ> 0 $ being the exponent of the nonlinearity. The sEWI is explicit and stable under a time step size restriction independent of the mesh size. We rigorously establish error estimates of the sEWI under various regularity assumptions on potential and nonlinearity. For ``good" potential and nonlinearity ($H^2$-potential and $σ\geq 1$), we establish an optimal second-order error bound in the $L^2$-norm. For low regularity potential and nonlinearity ($L^\infty$-potential and $σ> 0$), we obtain a first-order $L^2$-norm error bound accompanied with a uniform $H^2$-norm bound of the numerical solution. Moreover, adopting a new technique of \textit{regularity compensation oscillation} (RCO) to analyze error cancellation, for some non-resonant time steps, the optimal second-order $L^2$-norm error bound is proved under a weaker assumption on the nonlinearity: $σ\geq 1/2$. For all the cases, we also present corresponding fractional order error bounds in the $H^1$-norm, which is the natural norm in terms of energy. Extensive numerical results are reported to confirm our error estimates and to demonstrate the superiority of the sEWI, including much weaker regularity requirements on potential and nonlinearity, and excellent long-time behavior with near-conservation of mass and energy.
format Preprint
id arxiv_https___arxiv_org_abs_2310_20181
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle An explicit and symmetric exponential wave integrator for the nonlinear Schrödinger equation with low regularity potential and nonlinearity
Bao, Weizhu
Wang, Chushan
Numerical Analysis
35Q55, 65M15, 65M70, 81Q05
We propose and analyze a novel symmetric Gautschi-type exponential wave integrator (sEWI) for the nonlinear Schrödinger equation (NLSE) with low regularity potential and typical power-type nonlinearity of the form $ |ψ|^{2σ}ψ$ with $ ψ$ being the wave function and $ σ> 0 $ being the exponent of the nonlinearity. The sEWI is explicit and stable under a time step size restriction independent of the mesh size. We rigorously establish error estimates of the sEWI under various regularity assumptions on potential and nonlinearity. For ``good" potential and nonlinearity ($H^2$-potential and $σ\geq 1$), we establish an optimal second-order error bound in the $L^2$-norm. For low regularity potential and nonlinearity ($L^\infty$-potential and $σ> 0$), we obtain a first-order $L^2$-norm error bound accompanied with a uniform $H^2$-norm bound of the numerical solution. Moreover, adopting a new technique of \textit{regularity compensation oscillation} (RCO) to analyze error cancellation, for some non-resonant time steps, the optimal second-order $L^2$-norm error bound is proved under a weaker assumption on the nonlinearity: $σ\geq 1/2$. For all the cases, we also present corresponding fractional order error bounds in the $H^1$-norm, which is the natural norm in terms of energy. Extensive numerical results are reported to confirm our error estimates and to demonstrate the superiority of the sEWI, including much weaker regularity requirements on potential and nonlinearity, and excellent long-time behavior with near-conservation of mass and energy.
title An explicit and symmetric exponential wave integrator for the nonlinear Schrödinger equation with low regularity potential and nonlinearity
topic Numerical Analysis
35Q55, 65M15, 65M70, 81Q05
url https://arxiv.org/abs/2310.20181