Fock matrix dynamics

Abstract 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.

[1]  G. Scuseria,et al.  Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. II. Generalizations based on mass-weighting, idempotency, energy conservation and choice of initial conditions , 2001 .

[2]  Martin J. Field.  Constrained optimization of ab initio and semiempirical Hartree-Fock wave functions using direct minimization or simulated annealing , 1991 .

[3]  M. Parrinello,et al.  Surface solvation of halogen anions in water clusters: An ab initio molecular dynamics study of the Cl−(H2O)6 complex , 2001 .

[4]  J. Almlöf,et al.  Principles for a direct SCF approach to LICAO–MOab‐initio calculations , 1982 .

[5]  P. Pulay.  Improved SCF convergence acceleration , 1982 .

[6]  G. Scuseria,et al.  Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals , 2001 .

[7]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[8]  Matthias Krack,et al.  All-electron ab-initio molecular dynamics , 2000 .

[9]  Emily A. Carter,et al.  A COMPARISON OF CAR-PARRINELLO AND BORN-OPPENHEIMER GENERALIZED VALENCE BOND MOLECULAR DYNAMICS , 1995 .

[10]  Michele Parrinello,et al.  The Gaussian and augmented-plane-wave density functional method for ab initio molecular dynamics simulations , 1999 .

[11]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[12]  Shridhar R. Gadre,et al.  Structure and Stability of Water Clusters (H2O)n, n ) 8-20: An Ab Initio Investigation , 2001 .

[13]  Peter Pulay,et al.  Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules , 1969 .