Journal
JOURNAL OF CHEMICAL PHYSICS
Volume 135, Issue 4, Pages -Publisher
AMER INST PHYSICS
DOI: 10.1063/1.3605303
Keywords
carbon compounds; convergence of numerical methods; density functional theory; HF calculations; molecular dynamics method; organic compounds; SCF calculations; tight-binding calculations; water
Funding
- National Institutes of Health
- Office of Science of the U.S. Department of Energy (DOE) [DE-AC02-05CH11231]
- U.S. Department of Energy
Ask authors/readers for more resources
An important element determining the time requirements of Born-Oppenheimer molecular dynamics (BOMD) is the convergence rate of the self-consistent solution of Roothaan equations (SCF). We show here that improved convergence and dynamics stability can be achieved by use of a Lagrangian formalism of BOMD with dissipation (DXL-BOMD). In the DXL-BOMD algorithm, an auxiliary electronic variable (e. g., the electron density or Fock matrix) is propagated and a dissipative force is added in the propagation to maintain the stability of the dynamics. Implementation of the approach in the self-consistent charge density functional tight-binding method makes possible simulations that are several hundred picoseconds in lengths, in contrast to earlier DFT-based BOMD calculations, which have been limited to tens of picoseconds or less. The increase in the simulation time results in a more meaningful evaluation of the DXL-BOMD method. A comparison is made of the number of iterations (and time) required for convergence of the SCF with DXL-BOMD and a standard method (starting with a zero charge guess for all atoms at each step), which gives accurate propagation with reasonable SCF convergence criteria. From tests using NVE simulations of C2F4 and 20 neutral amino acid molecules in the gas phase, it is found that DXL-BOMD can improve SCF convergence by up to a factor of two over the standard method. Corresponding results are obtained in simulations of 32 water molecules in a periodic box. Linear response theory is used to analyze the relationship between the energy drift and the correlation of geometry propagation errors. (C) 2011 American Institute of Physics.[doi:10.1063/1.3605303]
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available