Fock matrix dynamics

An efficient method is suggested for direct ab initio Born–Oppenheimer molecular dynamics. It is based on extrapolating the Fock matrices from previous time steps to provide a start for the SCF procedure. Several polynomial extrapolation schemes are compared; low-degree schemes work best. SCF conver...

Full description

Saved in:
Bibliographic Details
Published in:Chemical physics letters Vol. 386; no. 4; pp. 272 - 278
Main Authors: Pulay, Peter, Fogarasi, Géza
Format: Journal Article
Language:English
Published: Elsevier B.V 11-03-2004
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:An efficient method is suggested for direct ab initio Born–Oppenheimer molecular dynamics. It is based on extrapolating the Fock matrices from previous time steps to provide a start for the SCF procedure. Several polynomial extrapolation schemes are compared; low-degree schemes work best. SCF convergence can be achieved in 2–3, sometimes less than two cycles. We have identified a phenomenon, hysteresis, a violation of the energy conservation and reversibility, due to systematic errors arising from incomplete SCF convergence. It can be eliminated by rescaling the velocities. Examples include water clusters and hydrated halogen ions, running simulations for several picoseconds.
ISSN:0009-2614
1873-4448
DOI:10.1016/j.cplett.2004.01.069