OSMOTIC AND ACTIVITY-COEFFICIENTS FROM EFFECTIVE POTENTIALS FOR HYDRATED IONS

Based on a method we previously suggested [Phys. Rev. E 52, 3730 (1995 )], effective interaction potentials between Na+ and Cl- ions have bee n derived from interionic radial distribution functions (RDF) in molec ular dynamics (MD) simulations of aqueous NaCl solution. The effective interaction potentials between the hydrated ions, which reproduce the original ion-ion RDF curves, can be used further to construct a corre sponding ionic solution in a much larger scale and to calculate any pr operties dependent on the structure of the electrolyte solution. In a subsequent Monte Carlo (MC) simulation, using the effective potentials , the osmotic and activity coefficients are calculated for the ions. C alculation of these properties directly from atomic MD or MC simulatio ns is beyond the capacity of the present computers due to the very lar ge number of molecules required in the simulations to obtain reliable results. A very good agreement with the experimental results is obtain ed. Effects of three-body interactions and concentration dependencies of the effective potentials are discussed.