Simulation study of ion pairing in concentrated aqueous salt solutions with a polarizable force field

The accuracy of empirical force fields is critical for meaningful molecular dynamics simulations of concentrated ionic solutions. Current models are typically developed on the basis of single ion properties such as the monohydrate energy in the gas phase, or the absolute hydration free energy at infinite dilution. However, the failure of these models to represent accurately the properties of concentrated solutions cannot be excluded. Here, these issues are illustrated for a polarizable potential based on classical Drude oscillators. To model accurately concentrated ionic solutions, the parameters of the potential functions are optimized to reproduce osmotic pressure data. The sodium-chloride potential of mean force in solution calculated from the empirically-adjusted model is consistent with the results from that calculated from ab initio CPMD simulations.