ICSM: An order [formula omitted] method for calculating electrostatic interactions added to TINKER

Abstract

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 [Y. Lin, A. Baumketner, S. Deng, Z. Xu, D. Jacobs, W. Cai, An image-based reaction field method for electrostatic interactions in molecular dynamics simulations of aqueous solutions, 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.