Acceleration schemes forab initiomolecular-dynamics simulations and electronic-structure calculations
- 15 October 1994
- journal article
- research article
- Published by American Physical Society (APS) in Physical Review B
- Vol. 50 (15), 10561-10573
- https://doi.org/10.1103/physrevb.50.10561
Abstract
We study the convergence and the stability of fictitious dynamical methods for electrons. first, we show that a particular damped second-order dynamics has a much faster rate of convergence to the ground state than first-order steepest-descent algorithms while retaining their numerical cost per time step. Our damped dynamics has efficiency comparable to that of conjugate gradient methods in typical electronic minimization problems. Then, we analyze the factors that limit the size of the integration time step in approaches based on plane-wave expansions. The maximum allowed time step is dictated by the highest frequency components of the fictitious electronic dynamics. These can result either from the larger wave vector components of the kinetic energy or from the small wave vector components of the Coulomb potential giving rise to the so called charge sloshing problem. We show how to eliminate large wave vector instabilities by adopting a preconditioning scheme in the context of Car-Parrinello ab initio molecular-dynamics simulations of the ionic motion. We also show how to solve the charge sloshing problem when this is present. We substantiate our theoretical analysis with numerical tests on a number of different silicon and carbon systems having both insulating and metallic character.Keywords
This publication has 17 references indexed in Scilit:
- Ab initiomolecular dynamics ford-electron systems: Liquid copper at 1500 KPhysical Review Letters, 1992
- Ab initiomolecular dynamics: Analytically continued energy functionals and insights into iterative solutionsPhysical Review Letters, 1992
- Adiabaticity in first-principles molecular dynamicsPhysical Review B, 1992
- Theory ofab initiomolecular-dynamics calculationsPhysical Review A, 1991
- Solution of Schrödinger’s equation for large systemsPhysical Review B, 1989
- Conjugate gradient minimization of the energy functional: A new method for electronic structure calculationPhysical Review B, 1989
- Molecular Dynamics andab initioTotal Energy CalculationsPhysical Review Letters, 1986
- Unified Approach for Molecular Dynamics and Density-Functional TheoryPhysical Review Letters, 1985
- A new method for diagonalising large matricesJournal of Physics A: General Physics, 1985
- Algorithm for Diagonalization of Large MatricesThe Journal of Chemical Physics, 1965