Symplectic Pseudospectral Time-Domain Scheme for Solving Time-Dependent Schrodinger Equation
This work provides a more accurate and stable numerical method for solving the Schrödinger equation, benefiting computational quantum mechanics simulations.
The authors developed a symplectic pseudospectral time-domain (SPSTD) scheme for solving the Schrödinger equation, combining FFT for spatial derivatives and high-order symplectic integrators for time evolution. Numerical tests on 1D quantum well and 3D harmonic oscillator showed advantages over traditional PSTD and FDTD methods, with infinite spatial accuracy and energy conservation for long-term simulations.
A symplectic pseudospectral time-domain (SPSTD) scheme is developed to solve Schrodinger equation. Instead of spatial finite differences in conventional finite-difference time-domain (FDTD) method, the fast Fourier transform is used to calculate the spatial derivatives. In time domain, the scheme adopts high-order symplectic integrators to simulate time evolution of Schrodinger equation. A detailed numerical study on the eigenvalue problems of 1D quantum well and 3D harmonic oscillator is carried out. The simulation results strongly confirm the advantages of the SPSTD scheme over the traditional PSTD method and FDTD approach. Furthermore, by comparing to the traditional PSTD method and the non-symplectic Runge-Kutta (RK) method, the explicit SPSTD scheme which is an infinite order of accuracy in space domain and energy-conserving in time domain, is well suited for a long-term simulation.