Ionic force field optimization based on single-ion and ion-pair solvation properties

Molecular dynamics simulations of ionic solutions depend sensitively on the force fields employed for the ions. To resolve the fine differences between ions of the same valence and roughly similar size and in particular to correctly describe ion-specific effects, it is clear that accurate force fields are necessary. In the past, optimization strategies for ionic force fields either considered single-ion properties (such as the solvation free energy at infinite dilution or the ion-water structure) or ion-pair properties (in the form of ion-ion distribution functions). In this paper we investigate strategies to optimize ionic force fields based on single-ion and ion-pair properties simultaneously. To that end, we simulate five different salt solutions, namely CsCl, KCl, NaI, KF, and CsI, at finite ion concentration. The force fields of these ions are systematically varied under the constraint that the single-ion solvation free energy matches the experimental value, which reduces the two-dimensional $\{\sigma,\epsilon\}$ parameter space of the Lennard Jones interaction to a one dimensional line for each ion. From the finite-concentration simulations, the pair-potential is extracted and the osmotic coefficient is calculated, which is compared to experimental data. We find a strong dependence of the osmotic coefficient on the force field, which is remarkable as the single-ion solvation free energy and the ion-water structure remain invariant under the parameter variation. Optimization of the force field is achieved for the cations Cs$^+$ and K$^+$, while for the anions I$^-$ and F$^-$ the experimental osmotic coefficient cannot be reached. This suggests that in the long run, additional parameters might have to be introduced into the modeling, for example by modified mixing rules.