High order finite difference algorithms for solving the Schrödinger equation in molecular dynamics

The view of considering global Pseudospectral methods (Sinc and Fourier) as the infinite order limit of local finite difference methods, and vice versa, finite difference as a certain sum acceleration of the pseudospectral methods is exploited to investigate high order finite difference algorithms f...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:The Journal of chemical physics 1999-12, Vol.111 (24), p.10827-10835
Hauptverfasser: Guantes, Raul, Farantos, Stavros C.
Format: Artikel
Sprache:eng
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:The view of considering global Pseudospectral methods (Sinc and Fourier) as the infinite order limit of local finite difference methods, and vice versa, finite difference as a certain sum acceleration of the pseudospectral methods is exploited to investigate high order finite difference algorithms for solving the Schrödinger equation in molecular dynamics. A Morse type potential for iodine molecule is used to compare the eigenenergies obtained by a Sinc Pseudospectral method and a high order finite difference approximation of the action of the kinetic energy operator on the wave function. Two-dimensional and three-dimensional model potentials are employed to compare spectra obtained by fast Fourier transform techniques and variable order finite difference. It is shown that it is not needed to employ very high order approximations of finite differences to reach the numerical accuracy of pseudospectral techniques. This, in addition to the fact that for complex configuration geometries and high dimensionality, local methods require less memory and are faster than pseudospectral methods, put finite difference among the effective algorithms for solving the Schrödinger equation in realistic molecular systems.
ISSN:0021-9606
1089-7690
DOI:10.1063/1.480446