Prediction of the concentration dependence of the surface tension and density of salt solutions: atomistic simulations using Drude oscillator polarizable and nonpolarizable models

Molecular simulations using Drude oscillator polarizable and nonpolarizable models for water and ions are carried out to predict the dependence of the surface tension on salt concentration. The polarizable water and ion models are based only on the classical Drude oscillators. The temperature dependence of the surface tension of water is examined for different water models. The dependence of salt densities on salt concentration is investigated through the nonpolarizable and Drude oscillator polarizable models. Finally, the reproduction of the surface tension of salt solution over a large range of concentrations is analyzed through a number of combinations between ions and water force fields. The structure of the interface is then discussed as a function of polarization effects. We establish here the inability of the Drude oscillator polarizable force fields to reproduce the salt concentration dependence of surface tension of NaCl aqueous solutions.