Solubility of KF and NaCl in water by molecular simulation.

The solubility of two ionic salts, namely, KF and NaCl, in water has been calculated by Monte Carlo molecular simulation. Water has been modeled with the extended simple point charge model (SPC/E), ions with the Tosi-Fumi model and the interaction between water and ions with the Smith-Dang model. The chemical potential of the solute in the solution has been computed as the derivative of the total free energy with respect to the number of solute particles. The chemical potential of the solute in the solid phase has been calculated by thermodynamic integration to an Einstein crystal. The solubility of the salt has been calculated as the concentration at which the chemical potential of the salt in the solution becomes identical to that of the pure solid. The methodology used in this work has been tested by reproducing the results for the solubility of KF determined previously by Ferrario et al. [J. Chem. Phys. 117, 4947 (2002)]. For KF, it was found that the solubility of the model is only in qualitative agreement with experiment. The variation of the solubility with temperature for KF has also been studied. For NaCl, the potential model used predicts a solubility in good agreement with the experimental value. The same is true for the hydration chemical potential at infinite dilution. Given the practical importance of solutions of NaCl in water the model used in this work, whereas simple, can be of interest for future studies.

[1]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[2]  E. Wang,et al.  Dissolution dynamics of NaCl nanocrystal in liquid water. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  M. P. Tosi,et al.  Ionic sizes and born repulsive parameters in the NaCl-type alkali halides—I: The Huggins-Mayer and Pauling forms , 1964 .

[4]  L. Dang.  Fluoride—fluoride association in water from molecular dynamics simulations , 1992 .

[5]  J. Mayer,et al.  Interatomic Distances in Crystals of the Alkali Halides , 1933 .

[6]  Kenneth S. Pitzer,et al.  Thermodynamics of electrolytes. I. Theoretical basis and general equations , 1973 .

[7]  L. Smits,et al.  Transport Numbers of Concentrated Sodium Chloride Solutions at 25 , 1966 .

[8]  A D J Haymet,et al.  Free energy of solvation of simple ions: molecular-dynamics study of solvation of Cl- and Na+ in the ice/water interface. , 2005, The Journal of chemical physics.

[9]  E. Sanz,et al.  The range of meta stability of ice-water melting for two simple models of water , 2005, 0902.3966.

[10]  K. Higashitani,et al.  Molecular Dynamics Simulations of Water at NaCl(001) and NaCl(011) Surfaces , 1998 .

[11]  H. Uchida,et al.  MD simulation of crystal growth of NaCl from its supersaturated aqueous solution , 2005 .

[12]  J. Rasaiah,et al.  Computer simulation studies of aqueous sodium chloride solutions at 298 K and 683 K , 2000 .

[13]  D. Zahn.  Atomistic mechanism of NaCl nucleation from an aqueous solution. , 2004, Physical review letters.

[14]  S. K. Prasad,et al.  Recent developments in discotic liquid crystals , 1999 .

[15]  D. Williams.  Metals, ligands, and cancer. , 1972, Chemical reviews.

[16]  C. Vega,et al.  The fluid–solid equilibrium for a charged hard sphere model revisited , 2003 .

[17]  C. Cavazzoni,et al.  Structure of NaCl and KCl concentrated aqueous solutions by ab initio molecular dynamics , 2004 .

[18]  David E. Smith,et al.  Computer simulations of NaCl association in polarizable water , 1994 .

[19]  C. Vega,et al.  A general purpose model for the condensed phases of water: TIP4P/2005. , 2005, The Journal of chemical physics.

[20]  J. Brodholt.  Molecular dynamics simulations of aqueous NaCl solutions at high pressures and temperatures , 1998 .

[21]  J. Mayer,et al.  Dispersion and Polarizability and the van der Waals Potential in the Alkali Halides , 1933 .

[22]  T. Straatsma,et al.  Free energy of ionic hydration: Analysis of a thermodynamic integration technique to evaluate free energy differences by molecular dynamics simulations , 1988 .

[23]  Richard M. Noyes,et al.  Thermodynamics of Ion Hydration as a Measure of Effective Dielectric Properties of Water , 1962 .

[24]  R. C. Weast.  CRC Handbook of Chemistry and Physics , 1973 .

[25]  Aatto Laaksonen,et al.  Concentration Effects in Aqueous NaCl Solutions. A Molecular Dynamics Simulation , 1996 .

[26]  D. Frenkel,et al.  Calculation of the melting point of NaCl by molecular simulation , 2003 .

[27]  C. Sagui,et al.  Calculation of ionic charging free energies in simulation systems with atomic charges, dipoles, and quadrupoles , 2003 .

[28]  I. R. Mcdonald,et al.  CORRIGENDUM: Rigid-ion models of the interionic potential in the alkali halides , 1974 .

[29]  C. Vega,et al.  Solid–fluid equilibria for quadrupolar hard dumbbells via Monte Carlo simulation , 1995 .

[30]  Lloyd L. Lee,et al.  Molecular Thermodynamics of Nonideal Fluids , 1988 .

[31]  J. Martí,et al.  A molecular dynamics simulation study of hydrogen bonding in aqueous ionic solutions , 2005 .

[32]  Daan Frenkel,et al.  New Monte Carlo method to compute the free energy of arbitrary solids. Application to the fcc and hcp phases of hard spheres , 1984 .

[33]  G. Ciccotti,et al.  Solubility of KF in water by molecular dynamics using the Kirkwood integration method , 2002 .

[34]  P. Jungwirth,et al.  Brine rejection from freezing salt solutions: a molecular dynamics study. , 2005, Physical review letters.

[35]  J. Rasaiah,et al.  Using simulation to study solvation in water , 2001 .

[36]  L. Degrève,et al.  LARGE IONIC CLUSTERS IN CONCENTRATED AQUEOUS NACL SOLUTION , 1999 .

[37]  J. Rasaiah,et al.  Computer simulation studies of aqueous solutions at ambient and supercritical conditions using effective pair potential and polarizable potential models for water , 2001 .

[38]  M. P. Tosi,et al.  Ionic sizes and born repulsive parameters in the NaCl-type alkali halides—II: The generalized Huggins-Mayer form☆ , 1964 .

[39]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[40]  C. Vega,et al.  The melting temperature of the most common models of water. , 2005, The Journal of chemical physics.

[41]  Alexander P. Lyubartsev,et al.  Solvation free energies of methane and alkali halide ion pairs: An expanded ensemble molecular dynamics simulation study , 1998 .

[42]  C. Vega,et al.  Phase diagram of water from computer simulation. , 2004, Physical review letters.

[43]  Aatto Laaksonen,et al.  Determination of solvation free energies by adaptive expanded ensemble molecular dynamics. , 2004, The Journal of chemical physics.

[44]  L. Dang,et al.  Molecular dynamics simulations of aqueous ionic clusters using polarizable water , 1993 .

[45]  C. Vega,et al.  Solid–fluid equilibrium for a molecular model with short ranged directional forces , 1998 .

[46]  C. Vega,et al.  A potential model for the study of ices and amorphous water: TIP4P/Ice. , 2005, The Journal of chemical physics.

[47]  Berend Smit,et al.  Accelerating Monte Carlo Sampling , 2002 .

[48]  Gerhard Hummer,et al.  Free Energy of Ionic Hydration , 1996 .

[49]  J. S. Rowlinson,et al.  Molecular Thermodynamics of Fluid-Phase Equilibria , 1969 .

[50]  C. Vega,et al.  Tracing the phase diagram of the four-site water potential (TIP4P). , 2004, The Journal of chemical physics.

[51]  John A. Zollweg,et al.  The Lennard-Jones equation of state revisited , 1993 .

[52]  L. Degrève,et al.  Structure of concentrated aqueous NaCl solution: A Monte Carlo study , 1999 .

[53]  E. Mastny,et al.  Direct calculation of solid-liquid equilibria from density-of-states Monte Carlo simulations. , 2005, The Journal of chemical physics.

[54]  R. C. Weast.  Handbook of chemistry and physics , 1973 .

[55]  L. Dang,et al.  Mechanism and Thermodynamics of Ion Selectivity in Aqueous Solutions of 18-Crown-6 Ether: A Molecular Dynamics Study , 1995 .