Thermodynamically Consistent Force Field for Coarse-Grained Modeling of Aqueous Electrolyte Solution

J Phys Chem B. 2019 Mar 14;123(10):2424-2431. doi: 10.1021/acs.jpcb.8b11190. Epub 2019 Mar 5.

Abstract

We propose a thermodynamically consistent methodology to parameterize interactions between charged particles inside the dissipative particle dynamics (DPD) formalism. We used osmotic pressure experimental data as a function of the salinity in order to optimize the interaction parameters. Results for NaCl aqueous solution show that both mean osmotic and activity coefficients of individual ions allowed the determination of Na+-water, Cl--water, and Na+-Cl- DPD repulsion parameters. A simple linear relationship between the hydration-free energies of ions and the ion-water repulsion parameters that allows the parameterization of the complete series of halide and alkaline ions is proposed. Two strategies have been used to obtain the anion-cation interaction parameters for halide and alkaline ions. In the first one, the parameters are obtained based on the numerical optimization of the anion-cation repulsion parameter with respect to the experimental osmotic pressure data (with mean average deviations <4%). Second, we propose a predictive approach based on the free-energy difference of hydration energies of anions and cations in the spirit of the law of matching water affinities [with a mean absolute relative deviation of about 13%, better than 6% if small ions (Li+ and F-) are removed].