Efficient Parallel Algorithms for Molecular Dynamics Simulations using Variable Charge Transfer Electrostatic Potentials
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.
AU: D. J. Keffer AU: J. W. Mintmire TI: Efficient parallel algorithms for molecular dynamics simulations using variable charge transfer electrostatic potentials SO: International Journal of Quantum Chemistry VL: 80 NO: 4-5 PG: 733-742 YR: 2000 CP: Copyright © 2000 John Wiley & Sons, Inc. ON: 1097-461X PN: 0020-7608 AD: Code 6189, U.S. Naval Research Laboratory, Washington, DC 20375-5342 DOI: 10.1002/1097-461X(2000)80:4/5<733::AID-QUA23>3.0.CO;2-Q US: http://dx.doi.org/10.1002/1097-461X(2000)80:4/5<733::AID-QUA23>3.0.CO;2-Q