Premium
Efficient parallel algorithms for molecular dynamics simulations using variable charge transfer electrostatic potentials
Author(s) -
Keffer D. J.,
Mintmire J. W.
Publication year - 2000
Publication title -
international journal of quantum chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.484
H-Index - 105
eISSN - 1097-461X
pISSN - 0020-7608
DOI - 10.1002/1097-461x(2000)80:4/5<733::aid-qua23>3.0.co;2-q
Subject(s) - charge (physics) , electrostatics , multipole expansion , molecular dynamics , computation , algorithm , statistical physics , computer science , physics , computational chemistry , chemistry , quantum mechanics
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. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 80: 733–742, 2000