ICSM: An order N method for calculating electrostatic interactions added to TINKER

We present an order N method for calculating electrostatic interactions that has been integrated into the molecular dynamics portion of the TINKER Molecular Modeling package. This method, introduced in a previous paper [J. Chem. Phys. 131 (2009) 154103] and termed the Image-Charge Solvation Model (ICSM), is a hybrid electrostatic approach that combines the strengths of both explicit and implicit representations of the solvent. A multiple-image method is used to calculate reaction fields due to the implicit part while the Fast Multipole Method (FMM) is used to calculate the Coulomb interactions for all charges, including the explicit part. The integrated package is validated through test simulations of liquid water. The results are compared with those obtained by the Particle Mesh Ewald (PME) method that is built in the TINKER package. Timing performance of TINKER with the integrated ICSM is benchmarked on bulk water as a function of the size of the system. In particular, timing analysis results show that the ICSM outperforms the PME for sufficiently large systems with the break-even point at around 30,000 particles in the simulated system.

[1]  Donald J. Jacobs,et al.  An image-based reaction field method for electrostatic interactions in molecular dynamics simulations of aqueous solutions. , 2009, The Journal of chemical physics.

[2]  Jay W. Ponder,et al.  Analysis and Application of Potential Energy Smoothing and Search Methods for Global Optimization , 1998 .

[3]  Bo Zhang,et al.  FMM-Yukawa: An adaptive fast multipole method for screened Coulomb interactions , 2009, Comput. Phys. Commun..

[4]  Jay W. Ponder,et al.  Algorithms for calculating excluded volume and its derivatives as a function of molecular conformation and their use in energy minimization , 1991 .

[5]  D. Zorin,et al.  A kernel-independent adaptive fast multipole algorithm in two and three dimensions , 2004 .

[6]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[7]  Wilfred F. van Gunsteren,et al.  Biomolecular Modeling: Goals, Problems, Perspectives , 2006 .

[8]  Wei Cai,et al.  Extending the fast multipole method for charges inside a dielectric sphere in an ionic solvent: High-order image approximations for reaction fields , 2007, J. Comput. Phys..

[9]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

[10]  J. Ponder,et al.  The NMR solution structure of intestinal fatty acid-binding protein complexed with palmitate: application of a novel distance geometry algorithm. , 1996, Journal of molecular biology.

[11]  Bo Zhang,et al.  FMM-Yukawa: An adaptive fast multipole method for screened Coulomb interactions , 2009, Comput. Phys. Commun..

[12]  WALTER GAUTSCHI.  Algorithm 726: ORTHPOL–a package of routines for generating orthogonal polynomials and Gauss-type quadrature rules , 1994, TOMS.

[13]  J. Ponder,et al.  An efficient newton‐like method for molecular mechanics energy minimization of large molecules , 1987 .

[14]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[15]  Harold L. Friedman,et al.  Image approximation to the reaction field , 1975 .

[16]  D. Jacobs,et al.  Ionic solvation studied by image-charge reaction field method. , 2011, The Journal of chemical physics.

[17]  L. Greengard.  The Rapid Evaluation of Potential Fields in Particle Systems , 1988 .

[18]  Wei Cai,et al.  Extending the fast multipole method to charges inside or outside a dielectric sphere , 2007, J. Comput. Phys..

[19]  L. Greengard,et al.  Regular Article: A Fast Adaptive Multipole Algorithm in Three Dimensions , 1999 .

[20]  L. Greengard,et al.  A new version of the Fast Multipole Method for the Laplace equation in three dimensions , 1997, Acta Numerica.

[21]  W. Cai,et al.  Discrete Image Approximations of Ionic Solvent Induced Reaction Field to Charges , 2007 .

[22]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[23]  L. Greengard,et al.  A new version of the fast multipole method for screened Coulomb interactions in three dimensions , 2002 .