In this paper, we study the Schr\"{o}dinger equation in the semiclassical regime and with multiscale potential function. We develop the so-called constraint energy minimization generalized multiscale finite element method (CEM-GMsFEM), in the framework of Crank-Nicolson (CN) discretization in time. The localized multiscale basis functions are constructed by addressing the spectral problem and a constrained energy minimization problem related to the Hamiltonian norm. A first-order convergence in the energy norm and second-order convergence in the $L^2$ norm for our numerical scheme are shown, with a relation between oversampling number in the CEM-GMsFEM method, spatial mesh size and the semiclassical parameter provided. Furthermore, we demonstrate the convergence of the proposed Crank-Nicolson CEM-GMsFEM scheme with $H/\sqrt{\Lambda}$ sufficiently small (where $H$ represents the coarse size and $\Lambda$ is the minimal eigenvalue associated with the eigenvector not included in the auxiliary space). Our error bound remains uniform with respect to $\varepsilon$ (where $0 < \varepsilon\ll 1$ is the Planck constant). Several numerical examples including 1D and 2D in space, with high-contrast potential are conducted to demonstrate the efficiency and accuracy of our proposed scheme.
翻译:暂无翻译