A splitting approach for the magnetic Schrödinger equation

详细信息    查看全文

• 作者：M. Caliari^a ; ^{marco.caliari@univr.it} ; A. Ostermann^b ; ^{alexander.ostermann@uibk.ac.at} ; C. Piazzola^b ; ^{chiara.piazzola@uibk.ac.at}
• 关键词：Magnetic Schr&ouml ; dinger equation ; Exponential splitting methods ; Convergence ; Fourier techniques ; Nonequispaced fast Fourier transform
• 刊名：Journal of Computational and Applied Mathematics
• 出版年：2017
• 出版时间：15 May 2017
• 年：2017
• 卷：316
• 期：Complete
• 页码：74-85
• 全文大小：502 K
• 卷排序：316

文摘

The Schrödinger equation in the presence of an external electromagnetic field is an important problem in computational quantum mechanics. It also provides a nice example of a differential equation whose flow can be split with benefit into three parts. After presenting a splitting approach for three operators with two of them being unbounded, we exemplarily prove first-order convergence of Lie splitting in this framework. The result is then applied to the magnetic Schrödinger equation, which is split into its potential, kinetic and advective parts. The latter requires special treatment in order not to lose the conservation properties of the scheme. We discuss several options. Numerical examples in one, two and three space dimensions show that the method of characteristics coupled with a nonequispaced fast Fourier transform (NFFT) provides a fast and reliable technique for achieving mass conservation at the discrete level.