A multiscale finite element method for the Schrödinger equation with multiscale potentials
Fuente:
arXiv
Enregistré dans:
| Auteurs principaux: | , , |
|---|---|
| Format: | Preprint |
| Publié: |
2019
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866916532879097856 |
|---|---|
| author | Chen, Jingrun Ma, Dingjiong Zhang, Zhiwen |
| author_facet | Chen, Jingrun Ma, Dingjiong Zhang, Zhiwen |
| contents | In recent years, an increasing attention has been paid to quantum heterostructures with tailored functionalities, such as heterojunctions and quantum matematerials, in which quantum dynamics of electrons can be described by the Schrödinger equation with multiscale potentials. The model, however, cannot be solved by asymptoics-based approaches where an additive form of different scales in the potential term is required to construct the prescribed approximate solutions. In this paper, we propose a multiscale finite element method to solve this problem in the semiclassical regime. The localized multiscale basis are constructed using sparse compression of the Hamiltonian operator, and thus are "blind" to the specific form of the potential term. After an one-shot eigendecomposition, we solve the resulting system of ordinary differential equations explicitly for the time evolution. In our approach, the spatial mesh size is $ H=\mathcal{O}(ε) $ where $ε$ is the semiclassical parameter and the time stepsize $ k$ is independent of $ε$. Numerical examples in one dimension with a periodic potential, a multiplicative two-scale potential, and a layered potential, and in two dimension with an additive two-scale potential and a checkboard potential are tested to demonstrate the robustness and efficiency of the proposed method. Moreover, first-order and second-order rates of convergence are observed in $H^1$ and $L^2$ norms, respectively. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_1901_00343 |
| institution | arXiv |
| publishDate | 2019 |
| record_format | arxiv |
| spellingShingle | A multiscale finite element method for the Schrödinger equation with multiscale potentials Chen, Jingrun Ma, Dingjiong Zhang, Zhiwen Numerical Analysis 65M60, 74Q10, 35J10 In recent years, an increasing attention has been paid to quantum heterostructures with tailored functionalities, such as heterojunctions and quantum matematerials, in which quantum dynamics of electrons can be described by the Schrödinger equation with multiscale potentials. The model, however, cannot be solved by asymptoics-based approaches where an additive form of different scales in the potential term is required to construct the prescribed approximate solutions. In this paper, we propose a multiscale finite element method to solve this problem in the semiclassical regime. The localized multiscale basis are constructed using sparse compression of the Hamiltonian operator, and thus are "blind" to the specific form of the potential term. After an one-shot eigendecomposition, we solve the resulting system of ordinary differential equations explicitly for the time evolution. In our approach, the spatial mesh size is $ H=\mathcal{O}(ε) $ where $ε$ is the semiclassical parameter and the time stepsize $ k$ is independent of $ε$. Numerical examples in one dimension with a periodic potential, a multiplicative two-scale potential, and a layered potential, and in two dimension with an additive two-scale potential and a checkboard potential are tested to demonstrate the robustness and efficiency of the proposed method. Moreover, first-order and second-order rates of convergence are observed in $H^1$ and $L^2$ norms, respectively. |
| title | A multiscale finite element method for the Schrödinger equation with multiscale potentials |
| topic | Numerical Analysis 65M60, 74Q10, 35J10 |
| url | https://arxiv.org/abs/1901.00343 |