Computer simulations of NaCl association in polarizable water

Abstract

Classical molecular dynamics computer simulations have been used to investigate the thermodynamics and kinetics of sodium chloride association in polarizable water. The simulations make use of the three-site polarizable water model of Dang [J. Chem. Phys. 97, 2659 (1992)], which accurately reproduces many bulk water properties. The model’s static dielectric constant and relaxation behavior have been calculated and found to be in reasonable agreement with experimental results. The ion–water interaction potentials have been constructed through fitting to both experimental gas-phase binding enthalpies for small ion–water clusters and to the measured structures and solvation enthalpies of ionic solutions. Structural properties and the potential of mean force for sodium chloride in water have been calculated. In addition, Grote–Hynes theory has been used to predict dynamical features of contact ion-pair dissociation. All of the calculated ionic solution properties have been compared with results from simulations using the extended simple point charge (SPC/E), nonpolarizable water model [J. Phys. Chem. 91, 6296 (1987)]. The dependence on polarizability is found to be small, yet measurable, with the largest effects seen in the solvation structure around the highly polarizable chlorine anion. This work validates the use of some nonpolarizable water models in simulations of many condensed-phase systems of chemical and biochemical interest.