Efficient Parallel Algorithms for Molecular Dynamics Simulations using Variable Charge Transfer Electrostatic Potentials

Document Type


Publication Date



We present an algorithm developed for the efficient computation on parallel computers of molecular dynamics simulations of materials with electrostatic potentials that allow for dynamic charge evolution. This algorithm combines a hierarchical cell multipole method for the fast evaluation of the electrostatic potential with a Broyden-Fletcher-Goldfarb-Shanno scheme for approximate solutions to arbitrary precision of the electronegativity equalization condition, which defines charge transfer in the material. We apply this algorithm to the simulation of a model -alumina slab at equilibrium. First, we demonstrate that simple fixed-charge models yield different bulk energies than do dynamic-charge models. Second, we compare our algorithm with two previously published techniques for the simulation of dynamic charges in electrostatic materials, discussing both computational speed and numerical accuracy. We show the relative computational speed-up of this scheme over a direct evaluation of the linear equations, and show that intrinsic flaws in other approximate dynamic-charge schemes can be avoided with this approach.

Files over 3MB may be slow to open. For best results, right-click and select "save as..."