Fast Computation Of The Fractional Schrodinger Equation On An Unbounded Domain

A fast and accurate numerical scheme is presented and analyzed for the computa-tion of the time fractional Schrodinger equation on an unbounded domain.The main idea consists of two parts.Firstly,we use artificial boundary methods to equivalent-ly reformulate the unbounded problem into an initial-boundary value(IBV)problem.Secondly,we present two numerical schemes for the IBV problem:a direct scheme and a fast scheme.The direct scheme stands for the direct discretization of the Caputo fractional derivative by using Ll-formula.The fast scheme means that the sum-of-exponentials approximation is used to speed up the evaluation of the Caputo fractional derivative.The resulting fast algorithm significantly reduces the storage requirement and the overall computational cost,and keeps almost the same accuracy as the direct scheme.Furthermore,the corresponding stability analysis and error estimates of two schemes are carried out,and numerical examples are provided to verify the performance of our approach and demonstrate that the fast scheme is more efficient than the direct scheme.