Removing systematic errors in interionic potentials of mean force computed in molecular simulations using reaction-field-based electrostatics.
The performance of reaction-field methods to treat electrostatic interactions is tested in simulations of ions solvated in water. The potential of mean force between sodium chloride pair of ions and between side chains of lysine and aspartate are computed using umbrella sampling and molecular dynamics simulations. It is found that in comparison with lattice sum calculations, the charge-group-based approaches to reaction-field treatments produce a large error in the association energy of the ions that exhibits strong systematic dependence on the size of the simulation box. The atom-based implementation of the reaction field is seen to (i) improve the overall quality of the potential of mean force and (ii) remove the dependence on the size of the simulation box. It is suggested that the atom-based truncation be used in reaction-field simulations of mixed media.
Baumketner
Andrij
A
Department of Physics and Optical Science, University of North Carolina Charlotte, 9201 University City Blvd., Charlotte, North Carolina 28269, USA. abaumket@uncc.edu
D001224:Aspartic Acid; D003198:Computer Simulation; D007477:Ions / Q000737:chemistry*; D008239:Lysine; D008958:Models, Molecular; D012965:Sodium Chloride; D055672:Static Electricity; D014867:Water / Q000737:chemistry*
D007477:Ions; D014867:Water; D001224:Aspartic Acid; D012965:Sodium Chloride; D008239:Lysine