Acceleration Schemes for Ab-Initio Molecular Dynamics and Electronic Structure Calculations

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 analyse 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 large wavevector components of the kinetic energy or from the small wavevector components of the Coulomb potential giving rise to the so called {\it charge sloshing} problem. We show how to eliminate large wavevector instabilities by adopting a preconditioning scheme that is implemented here for the first-time 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.

RevTex, 9 figures available upon request, to appear in Phys. Rev. B

References in corpus (1)

• Electronic structure calculations and molecular dynamics simulations with linear system-size scaling

Cited by in corpus (30)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Density-functional theory calculations for poly-atomic systems: Electronic structure, static and elastic properties and ab initio molecular dynamics
• Ab initio theory and modeling of water
• Real-Space Mesh Techniques in Density Functional Theory
• The Individual and Collective Effects of Exact Exchange and Dispersion Interactions on the Ab Initio Structure of Liquid Water
• Why does hydronium diffuse faster than hydroxide in liquid water?
• Order-N implementation of exact exchange in extended systems
• First principles study of adsorbed Cu_n (n=1-4) microclusters on MgO(100): structural and electronic properties
• Pressure and stress tensor of complex anharmonic crystals within the stochastic self-consistent harmonic approximation
• How well do Car-Parrinello simulations reproduce the Born-Oppenheimer surface ? Theory and Examples
• Local Structure Analysis in $Ab$ $Initio$ Liquid Water
• Density functional theory description of hole-trapping in SiO${}_2$: a successful self-interaction-corrected approach
• First-principle molecular dynamics with ultrasoft pseudopotentials: parallel implementation and application to extended bio-inorganic system
• Structural, Electronic, and Dynamical Properties of Liquid Water by ab initio Molecular Dynamics based on SCAN Functional within the Canonical Ensemble
• Comparison of Global and Local Adaptive Coordinates for Density Functional Calculations
• Isotope Effects in Liquid Water via Deep Potential Molecular Dynamics
• First-principle study of excitonic self-trapping in diamond
• Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance
• Electronic structure and Peierls instability in graphene nanoribbons sculpted in graphane
• Effectiveness of smearing and tetrahedron methods: best practices in DFT codes
• Spontaneous formation and stability of small GaP fullerenes
• Stabilized Quasi-Newton Optimization of Noisy Potential Energy Surfaces
• Accelerated ab-initio Molecular Dynamics: probing the weak dispersive forces in dense liquid hydrogen
• Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles
• A Real-Space Full Multigrid study of the fragmentation of Li11+ clusters
• Second generation wave-function thermostat for ab-initio molecular dynamics
• Accelerated First-Order Methods: Differential Equations and Lyapunov Functions
• Efficient Langevin dynamics for "noisy" forces
• Inflationary dynamics for matrix eigenvalue problems
• A concise introduction to molecular dynamics simulation: theory and programming