"Optimizing the accuracy and efficiency of fast hierarchical multipole expansions for MD simulations" Konstantin Lorenzen, Magnus Schwörer, Philipp Tröster, Simon Mates, and Paul Tavan
J. Chem. Theory Comput. 8, 36283636 (2012).
Abstract: Based on p’th order Cartesian Taylor expansions of Coulomb interactions and on hierarchical decompositions of macromolecular simulation systems into a hierarchy of nested,
structureadapted, and adaptively formed clusters of increasing size, fast multipole methods are constructed for rapid and accurate calculations of electrostatic interactions. These socalled
SAMM_{p} algorithms are formulated through totally symmetric and traceless tensors describing the multipole moments and the coefficients of local Taylor expansions. Simple recursions for the efficient evaluation and shifting of multipole moments are given. The required tensors are explicitly given up to order p = 4. The SAMM_{p} algorithms are shown to guarantee the reaction principle. For systems with periodic boundaries a reaction field (RF) correction is applied, which introduces at distances beyond the “minimum image convention” boundary a dielectric continuum surrounding each cluster at the toplevel of coarse graining. The correctness of the present SAMM_{p} implementation is demonstrated by analyzing the scaling of the residuals and by checking the numerical accuracy of the reaction principle for a pair of distant molecular ions in vacuum. Molecular dynamics (MD) simulations of pure water and of aqueous solutions
containing artificial ions, which are enclosed by periodic boundaries, demonstrate the stability
and lownoise behavior of SAMM_{p}/RF.
BMO authors (in alphabetic order): Konstantin Lorenzen Magnus Schwörer Paul Tavan Philipp Troester
Assoziierte Projekte: Polarizable force fields for molecular mechanics from first principles Longrange electrostatics in molecular dynamics simulations
WWWVersion
