matlab逆变换法产生随机数_分裂算符法数值求解含时薛定谔方程

author{Qinghua Ke}
date{2019, 09, 07}
上一篇文章我们已经简单介绍了虚时间演化法求解本征态，从而可以求得初态。

qinghuake：虚时间演化法 (imaginary-time propagation method) 计算本征态​zhuanlan.zhihu.com

接下来，在原子或分子上加上外电场，用分裂算符的方法 [1] 数值求解含时薛定谔方程。原子或分子外的电子在强激光场下二维的含时薛定谔方程（单电子近似），

其中哈密顿量（长度规范）

,

是电子的坐标，软核势

,

是软核参数。对于氩原子，

，计算得到相应的电离能

如果

很小，

时刻与

波函数的关系为

知道

时刻的波函数

（初始波函数由虚时演化得到），将波函数傅立叶变换到动量空间，乘以动量算符，然后逆变换到坐标空间，乘以势能项，接着变换到动量空间，乘以动量算符，最后变换回来，即可得到

时刻的波函数

. 当然，实际计算过程中，要将坐标空间分为内外两个部分。把每个时刻到外部空间的波函数变换到动量空间，演化到激光结束，将这些波函数相干叠加，取模的平方，就可以得到动量谱。

[1] M. D. Feit, J. A. Fleck JR, A Steiger , ``Solution of the Schrödinger equation by a spectral method'', J. Comput. Phys. 47, 412 (1982)