Convergence of a moving window method for the Schrödinger equation with potential on $\mathbb{R}^d$
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2024
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866911993545359360 |
|---|---|
| author | Iserles, Arieh Li, Buyang Yao, Fangyan |
| author_facet | Iserles, Arieh Li, Buyang Yao, Fangyan |
| contents | We propose a novel framework, called moving window method, for solving the linear Schrödinger equation with an external potential in $\mathbb{R}^d$. This method employs a smooth cut-off function to truncate the equation from Cauchy boundary conditions in the whole space to a bounded window of scaled torus, which is itself moving with the solution. This allows for the application of established schemes on this scaled torus to design algorithms for the whole-space problem. Rigorous analysis of the error in approximating the whole-space solution by numerical solutions on a bounded window is established. Additionally, analytical tools for periodic cases are used to rigorously estimate the error of these whole-space algorithms. By integrating the proposed framework with a classical first-order exponential integrator on the scaled torus, we demonstrate that the proposed scheme achieves first-order convergence in time and $γ/2$-order convergence in space for initial data in $H^γ(\mathbb{R}^d) \cap L^2(\mathbb{R}^d;|x|^{2γ} dx)$ with $γ\geq 2$. In the case where $γ= 1$, the numerical scheme is shown to have half-order convergence under an additional CFL condition. In practice, we can dynamically adjust the window when waves reach its boundary, allowing for continued computation beyond the initial window. Extensive numerical examples are presented to support the theoretical analysis and demonstrate the effectiveness of the proposed method. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2408_09463 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Convergence of a moving window method for the Schrödinger equation with potential on $\mathbb{R}^d$ Iserles, Arieh Li, Buyang Yao, Fangyan Numerical Analysis 35Q41, 65M12, 65M70 We propose a novel framework, called moving window method, for solving the linear Schrödinger equation with an external potential in $\mathbb{R}^d$. This method employs a smooth cut-off function to truncate the equation from Cauchy boundary conditions in the whole space to a bounded window of scaled torus, which is itself moving with the solution. This allows for the application of established schemes on this scaled torus to design algorithms for the whole-space problem. Rigorous analysis of the error in approximating the whole-space solution by numerical solutions on a bounded window is established. Additionally, analytical tools for periodic cases are used to rigorously estimate the error of these whole-space algorithms. By integrating the proposed framework with a classical first-order exponential integrator on the scaled torus, we demonstrate that the proposed scheme achieves first-order convergence in time and $γ/2$-order convergence in space for initial data in $H^γ(\mathbb{R}^d) \cap L^2(\mathbb{R}^d;|x|^{2γ} dx)$ with $γ\geq 2$. In the case where $γ= 1$, the numerical scheme is shown to have half-order convergence under an additional CFL condition. In practice, we can dynamically adjust the window when waves reach its boundary, allowing for continued computation beyond the initial window. Extensive numerical examples are presented to support the theoretical analysis and demonstrate the effectiveness of the proposed method. |
| title | Convergence of a moving window method for the Schrödinger equation with potential on $\mathbb{R}^d$ |
| topic | Numerical Analysis 35Q41, 65M12, 65M70 |
| url | https://arxiv.org/abs/2408.09463 |