Acceleration schemes for ab initio molecular-dynamics simulations and electronic-structure calculations

F. Tassone, F. Mauri, Roberto Car

Research output: Contribution to journalArticle

150 Scopus citations

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.

Original languageEnglish (US)
Pages (from-to)10561-10573
Number of pages13
JournalPhysical Review B
Volume50
Issue number15
DOIs
StatePublished - Jan 1 1994

All Science Journal Classification (ASJC) codes

  • Condensed Matter Physics

Fingerprint Dive into the research topics of 'Acceleration schemes for ab initio molecular-dynamics simulations and electronic-structure calculations'. Together they form a unique fingerprint.

  • Cite this