Scalable improvement of SPME multipolar electrostatics in anisotropic polarizable molecular mechanics using a general short‐range penetration correction up to quadrupoles. Issue 5 (15th February 2016)

• Record Type: Journal Article
• Authors: Narth, Christophe; Lagardère, Louis; Polack, Étienne; Gresh, Nohad; Wang, Qiantao; Bell, David R.; Rackers, Joshua A.; Ponder, Jay W.; Ren, Pengyu Y.; Piquemal, Jean‐Philip
• Is Part Of: Journal of computational chemistry. Volume 37:Issue 5(2016)
• Page Start: 494
• Page End: 506
• Publication Date: 2016-02-15
• Subjects: polarizable force fields -- electrostatics -- energy decomposition analysis; Chemistry -- Data processing -- Periodicals; 542.85
• DOI: 10.1002/jcc.24257
• Languages: English
• ISSNs: 0192-8651

Abstract

Abstract : We propose a general coupling of the Smooth Particle Mesh Ewald SPME approach for distributed multipoles to a short‐range charge penetration correction modifying the charge‐charge, charge‐dipole and charge‐quadrupole energies. Such an approach significantly improves electrostatics when compared to ab initio values and has been calibrated on Symmetry‐Adapted Perturbation Theory reference data. Various neutral molecular dimers have been tested and results on the complexes of mono‐ and divalent cations with a water ligand are also provided. Transferability of the correction is adressed in the context of the implementation of the AMOEBA and SIBFA polarizable force fields in the TINKER‐HP software. As the choices of the multipolar distribution are discussed, conclusions are drawn for the future penetration‐corrected polarizable force fields highlighting the mandatory need of non‐spurious procedures for the obtention of well balanced and physically meaningful distributed moments. Finally, scalability and parallelism of the short‐range corrected SPME approach are addressed, demonstrating that the damping function is computationally affordable and accurate for molecular dynamics simulations of complex bio‐ or bioinorganic systems in periodic boundary conditions. © 2016 Wiley Periodicals, Inc. Abstract : A general and scalable penetration correction has been implemented for distributed multipole electrostatics into Tinker to be used in conjunction to the smooth particle Mesh Ewald approach for periodic boundary conditions simulations with polarizable force fields such as AMOEBA and SIBFA.