Affiliation:
1. Department of Mathematics Iowa State University Ames Iowa USA
Abstract
AbstractNumerical solution of the time‐dependent Schrödinger equation with a position‐dependent effective mass is challenging to compute due to the presence of the non‐constant effective mass. To tackle the problem we present operator splitting‐based numerical methods. The wavefunction will be propagated either by the Krylov subspace method‐based exponential integration or by an asymptotic Green's function‐based time propagator. For the former, the wavefunction is given by a matrix exponential whose associated matrix–vector product can be approximated by the Krylov subspace method; and for the latter, the wavefunction is propagated by an integral with retarded Green's function that is approximated asymptotically. The methods have complexity per step with appropriate algebraic manipulations and fast Fourier transform, where is the number of spatial points. Numerical experiments are presented to demonstrate the accuracy, efficiency, and stability of the methods.
Funder
Division of Mathematical Sciences
Simons Foundation
Subject
Applied Mathematics,Computational Mathematics,Numerical Analysis,Analysis
Cited by
1 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献