Solubility of NaCl in water by molecular simulation revisited.

In this paper, the solubility of NaCl in water is evaluated by using computer simulations for three different force fields. The condition of chemical equilibrium (i.e., equal chemical potential of the salt in the solid and in the solution) is obtained at room temperature and pressure to determine the solubility of the salt. We used the same methodology that was described in our previous work [E. Sanz and C. Vega, J. Chem. Phys. 126, 014507 (2007)] although several modifications were introduced to improve the accuracy of the calculations. It is found that the predictions of the solubility are quite sensitive to the details of the force field used. Certain force fields underestimate the experimental solubility of NaCl in water by a factor of four, whereas the predictions of other force fields are within 20% of the experimental value. Direct coexistence molecular dynamic simulations were also performed to determine the solubility of the salt. Reasonable agreement was found between the solubility obtained from free energy calculations and that obtained from direct coexistence simulations. This work shows that the evaluation of the solubility of salts in water can now be performed in computer simulations. The solubility depends on the ion-ion, ion-water, and water-water interactions. For this reason, the prediction of the solubility can be quite useful in future work to develop force fields for ions in water.

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

[2]  J. Skinner,et al.  Water structure, dynamics, and vibrational spectroscopy in sodium bromide solutions. , 2009, The Journal of chemical physics.

[3]  Erik L. G. Wernersson,et al.  Effect of Water Polarizability on the Properties of Solutions of Polyvalent Ions: Simulations of Aqueous Sodium Sulfate with Different Force Fields. , 2010, Journal of chemical theory and computation.

[4]  P. Shepson,et al.  Halide and sodium ion parameters for modeling aqueous solutions in TIP5P-Ew water , 2010 .

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

[6]  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 .

[7]  P. Shepson,et al.  Ions at the ice/vapor interface , 2007 .

[8]  R. Netz,et al.  Ionic force field optimization based on single-ion and ion-pair solvation properties: going beyond standard mixing rules. , 2012, The Journal of chemical physics.

[9]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

[10]  T. Cheatham,et al.  Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations , 2008, The journal of physical chemistry. B.

[11]  S. Hanlon,et al.  Structural transitions of deoxyribonucleic acid in aqueous electrolyte solutions. II. The role of hydration. , 1975, Biochemistry.

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

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

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

[15]  C. Vega,et al.  Computing the free energy of molecular solids by the Einstein molecule approach: ices XIII and XIV, hard-dumbbells and a patchy model of proteins. , 2008, The Journal of chemical physics.

[16]  William R. Smith,et al.  Molecular simulations of aqueous electrolyte solubility: 1. The expanded-ensemble osmotic molecular dynamics method for the solution phase. , 2005, The journal of physical chemistry. B.

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

[18]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[19]  Julian W. Vincze,et al.  The nonmonotonic concentration dependence of the mean activity coefficient of electrolytes is a result of a balance between solvation and ion-ion correlations. , 2010, The Journal of chemical physics.

[20]  G. S. Manning.  The molecular theory of polyelectrolyte solutions with applications to the electrostatic properties of polynucleotides , 1978, Quarterly Reviews of Biophysics.

[21]  Fernando Bresme,et al.  The short range anion-H interaction is the driving force for crystal formation of ions in water. , 2009, The Journal of chemical physics.

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

[23]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

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

[25]  A. Panagiotopoulos,et al.  Simulations of phase transitions and free energies for ionic systems , 2008 .

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

[27]  J. Anwar,et al.  Uncovering molecular processes in crystal nucleation and growth by using molecular simulation. , 2011, Angewandte Chemie.

[28]  William R. Smith,et al.  Molecular simulation of aqueous electrolyte solubility. 2. Osmotic ensemble Monte Carlo methodology for free energy and solubility calculations and application to NaCl. , 2011, The journal of physical chemistry. B.

[29]  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 .

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

[31]  D. Frenkel,et al.  Rate of homogeneous crystal nucleation in molten NaCl. , 2005, The Journal of chemical physics.

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

[33]  P. Madden,et al.  Calculation of activities of ions in molten salts with potential application to the pyroprocessing of nuclear waste. , 2008, The journal of physical chemistry. B.

[34]  A. Ladd,et al.  Interfacial and co-existence properties of the Lennard-Jones system at the triple point , 1978 .

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

[36]  P. Jungwirth,et al.  Molecular dynamics simulations of freezing of water and salt solutions , 2007 .

[37]  P. Gallo,et al.  A route to explain water anomalies from results on an aqueous solution of salt. , 2010, The Journal of chemical physics.

[38]  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.

[39]  T. Steitz,et al.  The contribution of metal ions to the structural stability of the large ribosomal subunit. , 2004, RNA.

[40]  C. C. Hardin,et al.  Monovalent cation induced structural transitions in telomeric DNAs: G-DNA folding intermediates. , 1991, Biochemistry.

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

[42]  M. Coletta,et al.  Anion concentration modulates the conformation and stability of the molten globule of cytochrome c , 2003, JBIC Journal of Biological Inorganic Chemistry.

[43]  M. W. Chase.  NIST-JANAF thermochemical tables , 1998 .

[44]  J. Hansen,et al.  Ions in water: from ion clustering to crystal nucleation. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[45]  Ivo Nezbeda,et al.  The Lennard-Jones fluid: an accurate analytic and theoretically-based equation of state , 1994 .

[46]  D. Cole,et al.  Ion-specific effects under confinement: the role of interfacial water. , 2010, ACS nano.

[47]  Thomas E. Cheatham,et al.  Molecular Dynamics Simulations of the Dynamic and Energetic Properties of Alkali and Halide Ions Using Water-Model-Specific Ion Parameters , 2009, The journal of physical chemistry. B.

[48]  P. Ulbig,et al.  Monitoring Ice Nucleation in Pure and Salty Water via High-Speed Imaging and Computer Simulations , 2008 .

[49]  P. Gallo,et al.  Effect of concentration on the thermodynamics of sodium chloride aqueous solutions in the supercooled regime. , 2009, The Journal of chemical physics.

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

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

[52]  C. Vega,et al.  Calculation of the melting point of alkali halides by means of computer simulations. , 2012, The Journal of chemical physics.

[53]  Maria M. Reif,et al.  Computation of methodology-independent single-ion solvation properties from molecular simulations. IV. Optimized Lennard-Jones interaction parameter sets for the alkali and halide ions in water. , 2011, The Journal of chemical physics.

[54]  P. Gallo,et al.  Structural properties of high and low density water in a supercooled aqueous solution of salt. , 2011, Journal of Physical Chemistry B.

[55]  Andrew S. Paluch,et al.  A method for computing the solubility limit of solids: application to sodium chloride in water and alcohols. , 2010, The Journal of chemical physics.

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

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

[58]  C. Vega,et al.  The melting point of ice Ih for common water models calculated from direct coexistence of the solid-liquid interface. , 2006, The Journal of chemical physics.

[59]  C. Vega,et al.  Solubility of KF and NaCl in water by molecular simulation. , 2007, The Journal of chemical physics.

[60]  M. Lísal,et al.  Molecular simulation of aqueous electrolyte solubility. 3. Alkali-halide salts and their mixtures in water and in hydrochloric acid. , 2012, The journal of physical chemistry. B.

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

[62]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[63]  A. D. J. Haymet,et al.  The ice/water interface: A molecular dynamics simulation study , 1988 .

[64]  P. Shepson,et al.  Molecular dynamics simulation of ice growth from supercooled pure water and from salt solution , 2006, Annals of Glaciology.

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

[66]  C. Vega,et al.  Absence of superheating for ice Ih with a free surface: a new method of determining the melting point of different water models , 2006, 1304.5078.

[67]  C. Vega,et al.  Revisiting the Frenkel-Ladd method to compute the free energy of solids: the Einstein molecule approach. , 2007, The Journal of chemical physics.

[68]  Ignacio Tinoco,et al.  Measurement of the effect of monovalent cations on RNA hairpin stability. , 2007, Journal of the American Chemical Society.

[69]  P. Gallo,et al.  Molecular dynamics studies on the thermodynamics of supercooled sodium chloride aqueous solution at different concentrations , 2010, Journal of physics. Condensed matter : an Institute of Physics journal.

[70]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[71]  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 .

[72]  A. Ladd,et al.  Triple-point coexistence properties of the lennard-jones system , 1977 .