Journal
JOURNAL OF CHEMICAL PHYSICS
Volume 135, Issue 22, Pages -Publisher
AMER INST PHYSICS
DOI: 10.1063/1.3660689
Keywords
-
Funding
- U.S. Department of Energy (DOE) through the LANL LDRD
- Vetenskapsradet (VR), SeRC
- Kungliga Vetenskapsakademien (KVA)
Ask authors/readers for more resources
Geometric integration schemes for extended Lagrangian self-consistent Born-Oppenheimer molecular dynamics, including a weak dissipation to remove numerical noise, are developed and analyzed. The extended Lagrangian framework enables the geometric integration of both the nuclear and electronic degrees of freedom. This provides highly efficient simulations that are stable and energy conserving even under incomplete and approximate self-consistent field (SCF) convergence. We investigate three different geometric integration schemes: (1) regular time reversible Verlet, (2) second order optimal symplectic, and (3) third order optimal symplectic. We look at energy conservation, accuracy, and stability as a function of dissipation, integration time step, and SCF convergence. We find that the inclusion of dissipation in the symplectic integration methods gives an efficient damping of numerical noise or perturbations that otherwise may accumulate from finite arithmetics in a perfect reversible dynamics. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3660689]
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