含时Schrödinger方程的高阶辛FDTD算法研究

提出了一种新的算法——高阶辛时域有限差分法(SFDTD(3, 4): symplectic finite-difference time-domain)求解含时薛定谔方程.在时间上采用三阶辛积分格式离散, 空间上采用四阶精度的同位差分格式离散, 建立了求解含时薛定谔方程的高阶离散辛框架;探讨了高阶辛算法的稳定性及数值色散性.通过理论上的分析及数值算例表明:当空间采用高阶同位差分格式时, 辛积分可提高算法的稳定度;SFDTD(3, 4)法和FDTD(2, 4)法较传统的FDTD(2, 2)法数值色散性明显改善.对二维量子阱和谐振子的仿真结果表明: SFDTD(3, 4)法较传统的FDTD(2, 2)法及高阶FDTD(2, 4)法有着更好的计算精度和收敛性, 且SFDTD(3, 4)法能够保持量子系统的能量守恒, 适用于长时间仿真.

High-oder symplectic FDTD scheme for solving time-dependent Schrödinger equation

Using three-order symplectic integrators and fourth-order collocated spatial differences, a high-order symplectic finite-difference time-domain (SFDTD(3, 4)) scheme is proposed to solve the time-dependent Schrödinger equation. First, high-order symplectic framework for discretizing the Schrödinger equation is described. The numerical stability and dispersion analyses are provided for the FDTD(2, 2), FDTD(2, 4) and SFDTD(3, 4) schemes. The results are demonstrated in terms of theoretical analyses and numerical simulations. The spatial high-order collocated difference reduces the stability that can be improved by the high-order symplectic integrators. The SFDTD(3, 4) scheme and FDTD(2, 4) approach show better numerical dispersion than the traditional FDTD(2, 2) method. The simulation results of a two-dimensional quantum well and harmonic oscillator strongly confirm the advantages of the SFDTD(3, 4) scheme over the traditional FDTD(2, 2) method and other high-order approaches. The explicit SFDTD(3, 4) scheme, which is high-order-accurate and energy-conserving, is well suited for long-term simulation.