Ion Pairing in Molecular Simulations of Aqueous Alkali Halide Solutions

Abstract

Using classical molecular dynamics simulations, we study ion−ion interactions in water. We study the potentials of mean force (PMF) for the full set of alkali halide ion pairs, and in each case, we test different parameter sets for modeling both the water and the ions. Altogether, we compared 300 different PMFs. We also calculate association equilibrium constants (<i>K</i><sub>A</sub>) and compare them to two types of experiments. Of additional interest here was the proposition of Collins called the “law of matching water affinities”, where the relative affinity of ions in solution depends on the matching of cation and anion sizes. From observations on the relative depths of the free energies of the contact ion pair (CIP) and the solvent-shared ion pair (SIP), along with related solvent structure analyses, we find a good correlation with this proposition: small−small and large−large should associate in water, and small−large should be more dissociated.