Volume 13, Issue 1
Effect of the Reaction Field on Molecular Forces and Torques Revealed by an Image-Charge Solvation Model

Wei Song ,  Yuchun Lin ,  Andrij Baumketner ,  Shaozhong Deng ,  Wei Cai and Donald J. Jacobs

10.4208/cicp.290711.180711s

Commun. Comput. Phys., 13 (2013), pp. 129-149.

Preview Full PDF BiBTex 278 569
  • Abstract

We recently developed the Image-Charge Solvation Model (ICSM), which is an explicit/implicit hybrid model to accurately account for long-range electrostatic forces in molecular dynamics simulations [Linet al., J. Chem. Phys., 131, 154103,2009]. The ICSM has a productive spherical volume within the simulation cell for which key physical properties of bulk water are reproduced, such as density, radial distribution function, diffusionconstants and dielectric properties. Although the reaction´Čüeld (RF) is essential, it typically accounts for less than 2% of the total electrostatic force on a water molecule. This observation motivates investigating further the role of the RF within the ICSM. In this report we focus on distributions of forces and torques on water molecules as a function of distance from the origin and make extensive tests over a range of model parameters where Coulomb forces are decomposed into direct interactions from waters modeled explicitly and the RF. Molecular torques due to the RF typically account for 20% of the total torque, revealing why the RF plays an important role in the dielectric properties of simulated water. Moreover, it becomes clear that the buffer layer in the ICSM is essential to mitigate artifacts caused by the discontinuous change in dielectric constants at the explicit/implicit interface.


  • History

Published online: 2013-01

  • Keywords

  • AMS Subject Headings

  • Cited by