Molecular Simulation Study on Hofmeister Cations and the Aqueous Solubility of Benzene

We study the ion-specific salting-out process of benzene in aqueous alkali chloride solutions using Kirkwood-Buff (KB) theory of solutions and molecular dynamics simulations with different empirical force field models for the ions and benzene. Despite inaccuracies in the force fields, the simulations indicate that the decrease of the Setchenow salting-out coefficient for the series NaCl > KCl > RbCl > CsCl is determined by direct benzene-cation correlations, with the larger cations showing weak interactions with benzene. Although ion-specific aqueous solubilities of benzene may be affected by indirect ion-ion, ion-water, and water-water correlations, too, these correlations are found to be unimportant, with little to no effect on the Setchenow salting-out coefficients of the various salts. We further considered LiCl, which is experimentally known to be a weaker salting-out agent than NaCl and KCl and, therefore, ranks at an unusual position within the Hofrneister cation series. The simulations indicate that hydrated Li+ ions can take part of the benzene hydration shell while the other cations are repelled by it. This causes weaker Li+ exclusion around the solute and a resulting, weaker salting-out propensity of LiCl compared to that of the other salts. Removing benzene-water and benzene-salt electrostatic interactions in the simulations does not affect this mechanism, which may therefore also explain the smaller effect of LiCl, as compared to that of NaCl or KCl, on aqueous solvation and hydrophobic interaction of nonpolar molecules.