Optimal error bounds on the exponential wave integrator for nonlinear Schrödinger equations with highly singular potential
Fuente:
arXiv
Enregistré dans:
| Auteurs principaux: | , , |
|---|---|
| Format: | Preprint |
| Publié: |
2026
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866909021264412672 |
|---|---|
| author | Bao, Weizhu Wang, Chushan Wu, Yifei |
| author_facet | Bao, Weizhu Wang, Chushan Wu, Yifei |
| contents | We establish error estimates of the first-order exponential wave integrator (EWI) for the nonlinear Schrödinger equation (NLSE) with a highly singular potential in $\mathbb{R}^d$ with $1\leq d \leq 3$. Our results deal with singular potentials in $L^p_\text{loc}(\mathbb{R}^d)$ with $p>\frac{d}{2}$ and $p\geq 1$, which is (almost) the weakest regularity of the potential required by the well-posedness of the NLSE. First, for $L^p_\text{loc}$-potentials with $p>2$, we establish an optimal first-order $L^2$-norm convergence for the EWI, with the convergence order slightly reduced to $1^-$ when $p=2$. To the best of our knowledge, the optimal first-order convergence for the three-dimensional $L^2$-potential is for the first time in the literature. The optimality of such an error bound is two-fold: (i) the first-order $L^2$-norm convergence is optimal for the EWI (and its higher-order versions) under the given $L^2$-regularity assumption on the potential, and (ii) to achieve the first-order $L^2$-norm convergence for the EWI, such an assumption is optimally weak. For more singular potentials in $L^p_\text{loc}(\mathbb{R}^d)$ with $\frac{d}{2} < p < 2$ and $p\geq 1$, we prove that the $L^2$-norm convergence is (almost) of $(1-α)$-order when $d=1,2$, and of $(1-\frac{3}{2}α)$-order when $d=3$, where $α:=d(1/p - 1/2)$ when $d =1,2,3$, $p>1$ and $α:=\frac{1}{2}^+$ when $d=1$, $p=1$. Notably, this result pushes the error estimate to the threshold regularity of the potential that matches the threshold regularity for the well-posedness of the NLSE, which is also for the first time. Two main ingredients are adopted in the proof: (i) the use of discrete space-time Lebesgue spaces together with discrete Strichartz estimates to establish the stability of the numerical scheme, and (ii) the use of normal form transformation and frequency decompositions to obtain optimal error bounds. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2605_03355 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Optimal error bounds on the exponential wave integrator for nonlinear Schrödinger equations with highly singular potential Bao, Weizhu Wang, Chushan Wu, Yifei Numerical Analysis 35Q55, 65M15, 65M70, 81-08 We establish error estimates of the first-order exponential wave integrator (EWI) for the nonlinear Schrödinger equation (NLSE) with a highly singular potential in $\mathbb{R}^d$ with $1\leq d \leq 3$. Our results deal with singular potentials in $L^p_\text{loc}(\mathbb{R}^d)$ with $p>\frac{d}{2}$ and $p\geq 1$, which is (almost) the weakest regularity of the potential required by the well-posedness of the NLSE. First, for $L^p_\text{loc}$-potentials with $p>2$, we establish an optimal first-order $L^2$-norm convergence for the EWI, with the convergence order slightly reduced to $1^-$ when $p=2$. To the best of our knowledge, the optimal first-order convergence for the three-dimensional $L^2$-potential is for the first time in the literature. The optimality of such an error bound is two-fold: (i) the first-order $L^2$-norm convergence is optimal for the EWI (and its higher-order versions) under the given $L^2$-regularity assumption on the potential, and (ii) to achieve the first-order $L^2$-norm convergence for the EWI, such an assumption is optimally weak. For more singular potentials in $L^p_\text{loc}(\mathbb{R}^d)$ with $\frac{d}{2} < p < 2$ and $p\geq 1$, we prove that the $L^2$-norm convergence is (almost) of $(1-α)$-order when $d=1,2$, and of $(1-\frac{3}{2}α)$-order when $d=3$, where $α:=d(1/p - 1/2)$ when $d =1,2,3$, $p>1$ and $α:=\frac{1}{2}^+$ when $d=1$, $p=1$. Notably, this result pushes the error estimate to the threshold regularity of the potential that matches the threshold regularity for the well-posedness of the NLSE, which is also for the first time. Two main ingredients are adopted in the proof: (i) the use of discrete space-time Lebesgue spaces together with discrete Strichartz estimates to establish the stability of the numerical scheme, and (ii) the use of normal form transformation and frequency decompositions to obtain optimal error bounds. |
| title | Optimal error bounds on the exponential wave integrator for nonlinear Schrödinger equations with highly singular potential |
| topic | Numerical Analysis 35Q55, 65M15, 65M70, 81-08 |
| url | https://arxiv.org/abs/2605.03355 |