Application of high-order SF-SFDTD scheme for solving a time-dependent Schrödinger equation
Xie Guo-Da,Pan Pan,Ren Xin-Gang,Feng Nai-Xing,Fang Ming,Li Ying-Song,Huang Zhi-Xiang,,
DOI: https://doi.org/10.7498/aps.73.20230771
IF: 0.906
2024-01-01
Acta Physica Sinica
Abstract:The finite-difference time-domain with second-order numerical accuracy in time and space (FDTD(2, 2)) method has been extensively employed in the field of quantum mechanics to solve the Schrödinger equation. Nevertheless, the presence of the Courant-Friedrichs-Lewy (CFL) condition imposes limitations on the grid size in the computational space, thereby constraining the admissible range of time steps. Accordingly, the efficiency of the FDTD(2, 2) method is significantly compromised. In addition, the second-order numerical accuracy of the FDTD(2, 2) method in both time and space domains often results in significant error accumulation during calculations, thereby undermining the fidelity of the simulation results. To surmount the constraints imposed by the CFL stability conditions and enhance the accuracy of computations, a novel approach termed the SF-SFDTD(3, 4) method has been proposed. This method combines spatial filtering (SF) with the high-order symplectic finite-difference time-domain (SFDTD) method (denoted by the numerical accuracies of 3 and 4 in time and space, respectively). Its primary objective is to solve the time-dependent Schrödinger equation while ensuring time stability and scalability. The SF-SFDTD(3, 4) method obviates the need for further derivation of the iterative formula employed in the conventional SFDTD(3, 4) method. Therefore, the method under consideration exhibits a remarkable degree of compatibility with its traditional counterpart. It merely necessitates the inclusion of a spatial filtering operation during each numerical iteration to eliminate spatial high-frequency components arising from the utilization of time step sizes that fail to satisfy the CFL stability condition, thereby ensuring the stability of the numerical scheme. Moreover, when the time step value satisfies the CFL stability condition, the amplitude of the high-frequency component approaches zero, thereby exerting minimal influence on the accuracy of the computational results. The adoption of time steps that do not meet the CFL stability conditions leads to an amplification in the amplitude of the high-frequency component. However, this finding solely affects the stability of the computational results, and the elimination of these unstable high-frequency components scarcely impacts the accuracy of the computational results. SF-SFDTD(3, 4) preserves the simplicity and efficacy inherent in the traditional SFDTD(3, 4) methods, while simultaneously enhancing computational efficiency. Additionally, the numerical stability and dispersion error of the SFSFDTD(3, 4) method were subjected to theoretical analysis. Finally, the validity and efficacy of the proposed method were corroborated through numerical illustrations.
physics, multidisciplinary