Efficient numerical method for the Schr\"odinger equation with high-contrast potentials
Abstract
In this paper, we study the Schr\"odinger 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 L2 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. The convergence requires H/=O(54), t=O(54) if ≤ δ; while if δ<, the convergence requires H/=O(14δ), t=O(δ23/4) (where H represents the maximum diameter of coarse elements, is the minimal eigenvalue associated with the eigenvector not included in the auxiliary space, t is the time step, 0 < 1 is the Planck constant and δ describes the multiscale structure of the potential).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.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.