On the calculation of solubilities via direct coexistence simulations: Investigation of NaCl aqueous solutions and Lennard-Jones binary mixtures

The dynamic and energetic properties of the alkali and halide ions were calculated using molecular dynamics (MD) and free energy simulations with various different water and ion force fields including our recently developed water-model-specific ion parameters. The properties calculated were activity coefficients, diffusion coefficients, residence times of atomic pairs, association constants, and solubility. Through calculation of these properties, we can assess the validity and range of applicability of the simple pair potential models and better understand their limitations. Due to extreme computational demands, the activity coefficients were only calculated for a subset of the models. The results qualitatively agree with experiment. Calculated diffusion coefficients and residence times between cation−anion, water−cation, and water−anion showed differences depending on the choice of water and ion force field used. The calculated solubilities of the alkali−halide salts were generally lower than the true solubility of the salts. However, for both the TIP4PEW and SPC/E water-model-specific ion parameters, solubility was reasonably well-reproduced. Finally, the correlations among the various properties led to the following conclusions: (1) The reliability of the ion force fields is significantly affected by the specific choice of water model. (2) Ion−ion interactions are very important to accurately simulate the properties, especially solubility. (3) The SPC/E and TIP4PEW water-model-specific ion force fields are preferred for simulation in high salt environments compared to the other ion force fields.