High-Accuracy Numerical Methods and Convergence Analysis for Schrödinger Equation with Incommensurate Potentials

Numerical solving the Schrödinger equation with incommensurate potentials presents a great challenge since its solutions could be space-filling quasiperiodic structures without translational symmetry nor decay. In this paper, we propose two high-accuracy numerical methods to solve the time-dependent quasiperiodic Schrödinger equation. Concretely, we discretize the spatial variables by the quasiperiodic spectral method and the projection method, and the time variable by the second-order operator splitting method. The corresponding convergence analysis is also presented and shows that the proposed methods both have spectral convergence rates in space and second order accuracy in time, respectively. Meanwhile, we analyse the computational complexity of these numerical algorithms. One- and two-dimensional numerical results verify these convergence conclusions, and demonstrate that the projection method is more efficient.